Skip to main content

Integrated proteomics and metabolomics analysis of lumbar in a rat model of osteoporosis treated with Gushukang capsules



Gushukang (GSK) capsules are a Chinese patented medicine that is widely used in clinics for the treatment of osteoporosis (OP). Animal experiments have revealed that the bone mineral density of osteoporotic rats increase after treatment with GSK capsules. However, the specific mechanism and target of GSK in the treatment of osteoporosis are unclear. Further studies are needed.


Metabolomics (GC/MS) and proteomics (TMT-LC-MC/MC) with bioinformatics (KEGG pathway enrichment), correlation analysis (Pearson correlation matrix), and joint pathway analysis (MetaboAnalyst) were employed to determine the underlying mechanisms of GSK. The differential expression proteins were verified by WB experiment.


The regulation of proteins, i.e., Cant1, Gstz1, Aldh3b1, Bid, and Slc1a3, in the common metabolic pathway of differential proteins and metabolites between GSK/OP and OP/SHAM was corrected in the GSK group. The regulation of 12 metabolites (tyramine, thymidine, deoxycytidine, cytosine, L-Aspartate, etc.) were differential in the common enrichment metabolic pathway between GSK /OP and OP/SHAM. Differential proteins and metabolites jointly regulate 11 metabolic pathways, such as purine metabolism, pyrimidine metabolism, histidine metabolism, beta-alanine metabolism, and so on.


GSK may protect bone metabolism in osteoporotic rats by affecting nucleotide metabolism, amino acid metabolism, and the immune system.

Peer Review reports


Osteoporosis (OP) is a systemic metabolic disease characterized by low bone mass, microstructural damage of bone tissue, decreased bone strength, and increased bone fragility [1]. Its clinical manifestations include primarily pain and spinal curvature [2]. For thousands of years, traditional Chinese medicine (TCM) has been used in both China and other Asian countries [3]. At present, TCM has obvious curative effects and benefits in the treatment of osteoporosis. These can better improve the side effects and safety of standard treatment. TCM formulas not only reduce bone loss by decreasing bone resorption but also increase bone formation in the multicomponent and multitarget pattern [4]. It can also regulate the overall function of the human body and relieve back pain and lumbago [5].

The GSK capsule is a Chinese patented medication listed in the Chinese Pharmacopoeia that is used for treating osteoporosis. Its main effective components include Epimedii folium, Rehmannia radix praeparata, Astragalus membranaceus (Fisch.)Bge. and Drynaria rhizoma. Chinese scholars have established an HPLC-CAD method to determine the content of the main components of the GSK capsules, which provides a reference for the quality control of the GSK capsules [6]. There is sufficient evidence that the main components of the GSK capsules, such as Epimedium and Drynaria, have antiosteoporosis function. Epimedium has been proved to improve bone mineral density or bone metabolism [7]. Rhizoma Drynariae can promote osteogenesis differentiation by activating the Wnt/β-catenin signaling pathways, which play a role in antiosteoporosis activity [8]. Animal experiments have shown that Rehmannia radix praeparata could significantly increases the expression levels of Akt1 and MAPK1 mRNA in the bone tissue to increase bone density [9]. Astragalus saponin IV (AS- IV) extracted from Astragalus membranaceus (Fisch.)Bge., promotes the osteogenic differentiation of the bone marrow mesenchymal stem cells via the MIR-21/NGF/BMP2/Runx2 pathway [10]. GSK capsules provide convenient administration, stable curative effects, and few side effects. In the guideline for the clinical diagnosis and treatment of primary osteoporosis in China, GSK is selected as the recommended Chinese patent medicine [11], and it is widely used in the clinical settings [12]. Clinical research shows that GSK is effective in the treatment of primary osteoporosis, by adjusting the whole-body function. GSK can inhibit the formation of osteoclasts and stimulate the formation of osteoblasts in ovariectomized mice [13]. It also significantly increases the production of vitamin D and calcium [14]. The BMP/Smad signaling pathway is important in bone remodeling. Interestingly, a study found that GSK significantly enhanced the BMP-2/Smad signaling pathway and improved bone microstructure by upregulating relevant osteogenic factors in osteoporotic rats [15]. Furthermore, H-type blood vessels have shown the ability to induce angiogenesis and bone formation [16]. Recent studies demonstrated that GSK can enhance hypoxia-inducible factor-1α to induce H-type vascular formation as well as bone formation [17].

In recent years, multiomics integrated analysis technology has been widely used to study various disease states [18]. Integrating multiomic data analyses can make up for data loss, noise, and other factors in single omic data analysis; integrate various interactions isolated at the gene level or protein level, various metabolic, and regulatory pathways; and jointly clarify the overall state of the biological system [19]. The combined application of proteomic and metabolomic technology has undeniable advantages [20]. Data from metabolomic analysis can provide information on molecular function for the results of proteomic research and provide clues to explore the regulatory relationship between them. The resulting protein–metabolite interaction network intuitively shows the involved, maladjusted pathways in the disease state or the overall changes after drug treatment. Proteomics and metabolomics provide us with vast information of differential proteins and metabolites. In turn, this provides a new research strategy for the target research of TCM in the treatment of diseases. However, the specific mechanism and target of GSK in OVX rats remain unclear. Therefore, based on the integrated analysis of proteomics and metabolomics, this study used GSK capsules to further explore the mechanism of GSK in the treatment of osteoporotic rats. The research flowchart is shown in Fig. 1.

Fig. 1
figure 1

A graphic abstract of the article


Drugs and reagents

Gushukang capsule was supplied by Liaoning Kangchen Pharmaceutical Co., Ltd. (Z20060270) (Dandong, China). BCA quantitative kit was purchased from Beyotime Biotech Inc. (Shanghai, China). SDT lysate, TMT labeling kit were from Thermo Fisher Scientific Inc. (Waltham, USA). Acetonitrile was from Merck & Co Inc. (Kenilworth, USA), Ammonium acetate was from Sigma (Aldrich, USA). Anti-Slc1a3 rabbit polyclonal antibody [EPR12686] was from Abcam Plc ((Cambridge, USA), anti-Aldh3b1 rabbit polyclonal antibody [TA323583S] was from OriGene Technologies, Inc. (Rockwell, USA), anti-Gstz1 rabbit polyclonal antibody [No.Ag6676], anti-Cant1 rabbit polyclonal antibody [No.12164–1-AP], anti-Bid rabbit polyclonal antibody [No.10988–1-AP] were from Proteintech Group, Inc. (Rosement, USA). Anti-β-actin rabbit polyclonal antibody [bs-0061R], Goat anti-rabbit IgG H&L HRP [bs-0295G-HRP] were from Beijing Biosynthesis Biotech Co., Ltd. (Beijing, China). MonPro™ ECL ultrasensitive substrate pro [PW30701] was from Monad Biotech Co., Ltd. (Wuhan, China).

Animal model establishment

All experimental procedures complied with the revised Animals (Scientific Procedures) Act 1986 and the ARRIVE (Animal Research: Reporting of In Vivo Experiments) guidelines, and were approved by the Animal Ethics Committee of Fujian Academy of Traditional Chinese Medicine (Fujian, China) (FJATCM-IAEC2018034). Anesthesia methods that we employed in the experiment complied with American Veterinary Medical Association (AVMA) Guidelines for the Euthanasia of Animals (2020). Thirty-six, 3-month-old female SPF Sprague Dawley rats were obtained from Shanghai Slake Laboratory Animal Co., Ltd. (Shanghai, China) [SCXK(Shanghai)2007–0005], and raised to 6 months old in the Experimental Center of Comparative Medicine, Fujian Academy of Chinese Medical Sciences (Fujian, China) [SYXK(Fujian)2009–000]. All rats were maintained in a 12-h light/dark cycle, with controlled temperature (22 °C–24 °C) and humidity (50–60%). After 1 week of adaptive feeding, 6-month-old rats (n = 36) were randomly divided into three groups: the GSK group (n = 12), the model group (n = 12), and the sham operation group (n = 12). An osteoporosis model was established by removing the bilateral ovaries of rats. Approximately 1 g of adipose tissue adjacent to the ovary was removed in the sham group [21]. Intramuscular injection of penicillin at 800,000 units/day was used to prevent infection for 3 days after the ovariectomy.

GSK treatment

The GSK capsules comprise seven medicinal compositions as shown in Table 1. GSK capsules were opened and mixed with normal saline, and administered by gavage once a day (0.32 g/(kg × d)). Both the sham group and model group were given an equal volume of normal saline once a day for 12 weeks. After 12 weeks of treatment, rats in each group were anesthetized with 2% pentobarbital sodium. The first and second lumbar vertebrae were separated and quickly put into liquid nitrogen. The left tibia was used for micro-CT.

Table 1 Ingredients of Gushukang capsules

Micro-CT bone analysis

Micro-computed tomography (micro-CT) scanning and morphometric analysis were performed on a micro-CT imaging system (ZKKS-MCT-Sharp, Zhongkekaisheng Co., China), which operated at a voltage of 70 kV and an electric current of 100 μA. The voxel size after reconstruction was 10 × 10 × 10 μm. Based on the micro-CT results, three-dimensional images were reconstructed by micro-CT reconstruction. The bone morphometric parameters of the left tibia were measured by analysis software, including bone mineral density (BMD), bone volume/total volume (BV/TV), trabecular separation (Tb.Sp) and trabecular number (Tb.N). The data are displayed as means ± standard errors of the mean (SEM). All statistical analyses were performed on Prism 8.0 (GraphPad Prism Software, USA).

Proteomic sample processing

Lumbar samples were added with appropriate SDT lysate. After boiling in a water bath, the supernatant was centrifuged to extract the protein. The lumbar vertebrae proteins of every four rats were combined into one sample for proteomic analysis. Protein was quantified using the BCA method. An appropriate amount of protein from each sample was trypsinized using the filter-aided proteome preparation (FASP) method [22], and the peptides were quantified (OD280). Approximately 100 UG peptide segments were taken from each sample and labeled according to the instructions of the TMT labeling kit.

Liquid chromatography tandem mass spectrometry (LC-MS/MS)

Each fractionated sample was separated by the HPLC liquid-phase system EASY-nLC (Thermo Fisher Scientific, USA) with nanoliter flow rate. After chromatographic separation, the samples were analyzed by a QExactive mass spectrometer (Thermo Fisher Scientific, USA). Positive ion detection method, parent ion scanned range 300–1800 m/Z. Mass charge ratio collection method of polypeptide and polypeptide fragments: collect 20 fragment maps (MS2 scan, HCD) after each full scan.

Proteomic data analysis

The original data was analyzed by LC-MS/MS software, Mascot 2.2, and Proteome Discoverer 1.4. The protein database is UniProt_ Rattus_ norvegicus_36107_ 20,190,524. The quantification method was based on the median quantitative value of the unique peptide. The maximum number of missed cuts is two. The variable modifications were oxidation (m) and TMT6/10 plex (y). The fixed modifications include carbamide methyl (c), TMT6/10 plex (N-term), and TMT6/10 plex (k). The screening criteria of credible peptides were ≤ 0.01.

Screening of differential proteins

Differentially expressed proteins were screened using a multiplier variation that was higher than 1.2 times (up to 1.2 times or less than 0.83) and with a p-value of < 0.05. Simultaneously, hierarchical clustering was used to analyze differentially expressed proteins in the comparison group, after which obtained data were displayed on a heatmap.

Metabonomic sample processing

Approximately 80-mg samples were added to 200-ul water. Samples were then homogenated with an MP homogenizer and vortexed for 60 s. Approximately 800-ul methanol acetonitrile solution (1:1, v/v) was added, and the mixture was vortexed for another 60 s. The suspension was ultrasonicated at low temperature for 30 min twice and placed at − 20 °C for 1 h to precipitate the protein. The suspension was then centrifuged at 4 °C, 14000 RCF, for 20 min. Finally, the supernatant was freeze-dried and stored at − 80 °C.

Gas chromatography mass spectrometry

The samples were separated in an HILIC column by 1290 Infinity LC (Agilent, USA) ultra-performance liquid chromatography (UHPLC). The column temperature was 25 °C, and the flow rate was 0.3 ml/min. The mobile phase consisted of either water + 25 mM ammonium acetate + 25 mm ammonia water (A) or acetonitrile, gradient elution (B). QC samples were inserted into the sample queue to evaluate both the stability of the system and the reliability of experimental data. The positive and negative ion modes of electrospray ionization (ESI) were detected. Samples were separated by UHPLC and then analyzed by a TripleTOF 5600 mass spectrometer (AB SCIEX). Second-order mass spectrometry was obtained by information-dependent acquisition (IDA) with the high sensitivity mode. The following settings were used: declustering potential of ±60 V (positive and negative modes) and collision energy of 35 ± 15 eV. IDA was set as follows: Exclude isotopes within 4 Da Candidate ions to monitor per cycle: 6.

Metabolomic data processing

The original data was converted into an mzXML format by ProteoWizard. Peak alignment, retention time correction, and peak area extraction were performed by the XCMS program. Metabolite structure identification uses accurate mass matching (< 25 ppm) and secondary spectrum matching methods while searching the laboratory’s self-built database. After pareto-scaling pretreatment, the data was analyzed using unsupervised principal component analysis (PCA).

Screening of differential metabolites

We used the PLS-DA data model to further screen differential metabolites among our three groups to obtain variable importance in projection (VIP) under both positive and negative ion mode conditions. Afterward, metabolites with a multidimensional statistical analysis of VIP > 1 and univariate statistical analysis p-value < 0.05 were selected to identify potential metabolites with significant differences. Subsequently, MetaboAnalyst 5.0 was used to analyze metabolic pathways involved in differential metabolites, and then to evaluate their importance. When the impact value of the understudied metabolic pathways was > 0.2, we considered these pathways as potential target metabolic pathways in this experiment.

Integration analysis of differential proteins and metabolites and the common pathway enrichment

For pathway enrichment analysis, dysregulated metabolites and proteins were selected and mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [23,24,25] (

Correlation analysis between metabolomics and proteomics

The correlation coefficient between significantly differential proteins and metabolites was calculated using Pearson’s correlation, after which results were displayed in the form of a correlation coefficient matrix heatmap (the software was R3.4.2). To intuitively reflect the difference in expression patterns of the significant differences obtained between differential protein and metabolites, we then conducted Pearson’s correlation hierarchical clustering (the software was a R3.4.2 heatmap). Subsequently, all nodes were loaded onto Cytoscape for network construction based on the correlation data. Afterward, we screened metabolites and proteins with significant differences at critical nodes in the network.

Detection of differential proteins by western blotting

Total protein in the lumbar was extracted with RIPA buffer. A BCA protein assay kit was used for concentration measurement. 4–12% PAGE gel was used for electrophoresis. The protein was transferred to PVDF membrane using eBlot L1 Fast Wet Transfer System (Genescript biotech Corp., China). After blocking with 5% skim milk dissolved in 1× Tris-buffered saline and 0.1% Tween 20, the PVDF membranes were cut prior to hybridisation with antibodies. The cropped PVDF membranes were incubated with the diluted primary antibody [Anti-Slc1a3 (11000 dilution), anti-Aldh3b1 (11000), anti-Cant1 (11000), anti-Gstz1(11000), anti-Bid (1: 1000) and anti-β-actin (14000) antibodies] at 4 °C overnight. The membrane was then washed with 1 × TBST three times and incubated with secondary antibody for 1 h. After washing three times with 1 × TBST, MonPro™ ECL ultrasensitive substrate pro was added dropwise and exposed in the dark room for 1 min. Results were visualized using the Fluor Chem M Digital Darkroom (ProteinSimple Ltd., USA).. Stained bands were analyzed with Image J to calculate the gray value. Gray values of aim protein/β-actin were taken as the relative value of differential protein expression in the lumbar. All statistical analyses were performed on Prism 8.0.


Effect of the GSK capsules on bone trabecular in osteoporotic rats

The proximal end of the left tibia was evaluated with a micro-CT plain scan (Fig. 2A, B, C). From the analysis results of bone trabecular morphological parameters, the BMD (P < 0.001, Fig. 2D), BV/TV (P < 0.01, Fig. 2E), BS/TV (P < 0.05, Fig. 2F) of the OP group decreased significantly, while Tb.Pf (P < 0.05, Fig. 2G) and Tb.Sp (P < 0.01, Fig. 2I) increased significantly compared with the SHAM group. In the OP group, Tb.N tended to decrease, but there was no statistical difference (P > 0.05, Fig. 2H).

Fig. 2
figure 2

Bone plain scan and bone mineral density (BMD) of the trabecular bone in the proximal left tibia. Representative plain scan images of the proximal tibia were shown (sagittal) in A SHAM, B OP, and C GSK. The BMD indices and parameters of micro CT are shown as the mean ± SEM (n = 3 per group). ***p < 0.001 vs. the sham group; **p < 0.01 vs. the sham group; *p < 0.05 vs. the sham group; # p < 0.05 vs. the OP group

In the GSK group, the BMD, BV/TV and BS/TV increased significantly compared with the OP group (P < 0.05, resptively, Fig. 2D, E, F), while Tb.Pf decreased significantly (P < 0.05, Fig. 2G). Tb.Sp showed a decreasing trend in GSK group, but there was no significant difference compared with OP group (P = 0.09, Fig. 2I).

Results of differential protein analysis

We identified 3614 proteins in the three groups. Subsequently, we then compared the differentially expressed proteins between these groups using the differential protein screening criteria mentioned in the above experimental method (Fig. 3), followed by a statistical analysis of the expression of differential proteins among GSK, OP, and SHAM groups. As observed, differentially expressed proteins in OVX rats were corrected in the GSK group (top 5) (Table 2). Table 3 lists the five proteins with the most significant upregulation and the five with the greatest downregulation in the OP group.

Fig. 3
figure 3

Different expression trends of differential proteins in each group (A). Cluster analysis of differentially expressed proteins in each group (B) (C). Hierarchical clustering results were represented as a tree heat map. Each row represents a protein, each column represents a set of samples, and Log2expression values of significant differentially expressed proteins in different samples are presented in different colors on the heat map. Furthermore, red represents significantly upregulated proteins, blue represents significantly downregulated proteins, and gray represents no quantitative protein information

Table 2 Statistical results of upregulated and downregulated proteins between GSK and OVX (top 5)
Table 3 Statistics of upregulated and downregulated proteins between OVX and SHAM (top5)

Results of differential metabolite analysis

Figure 4A and B shows the number of abnormal metabolites in positive and negative ion modes. PCA results indicated that the QC and other samples were well aggregated within the group, indicating that the overall quality of the experimental data met the analysis requirements. Significant differences were observed in metabolite levels among the GSK, SHAM, and OP groups (Fig. 4C, D). The significantly different metabolites (Fig. 5A) were mainly involved in the following pathways (Fig. 5B): phenylalanine, tyrosine, tryptophan, glycine, serine and threonine metabolism, phenylalanine metabolism, synthesis and degradation of ketone bodies, beta-alanine metabolism, cysteine and methionine metabolism, and pyruvate metabolism.

Fig. 4
figure 4

Differing expression trends of differential metabolites in each group (A) (B). PCA map showing metabolomics results in positive and negative ion modes (C) (D). Each point in the figure represents a sample. Different colors represent different groups. Furthermore, closer samples represent closer metabolite expression patterns, which means more minor differences between the groups. We used the PCA method to evaluate the aggregation degree and the overall distribution trend of samples within the group

Fig. 5
figure 5

Cluster analysis showing differentially expressed metabolites in the three groups (A). Map showing significant metabolic pathways involved in differential metabolites (B)

The number of common metabolic pathways of differential proteins and differential metabolites

Differential proteins between GSK and OP groups were involved in 70 metabolic pathways (Fig. 6A). However, although these differential metabolites were involved in 73 metabolic pathways, two parts jointly participated in 45 pathways. The results also showed that differential proteins between OP and SHAM groups were involved in 93 metabolic pathways (Fig. 6B). Nevertheless, while differential metabolites were involved in 65 metabolic pathways, two parts were jointly involved in 26 pathways. Figure 6C, D shows the top 10 pathways.

Fig. 6
figure 6

Venn diagrams showing the number of metabolic pathways involved in differential proteins and differential metabolites in GSK/OP and OP/SHAM (A) (B). The number of differential proteins and differential metabolites contained in these metabolic pathways (TOP 10) (C) (D)

Correlation analysis of significantly different proteins and metabolites

We constructed a correlation network of proteins and metabolites (|r| ≥ 0.5 and P < 0.05) (Table 4), used matrix heatmap and hierarchical heatmaps to represent their correlation coefficients (Fig. 7A-D), and found that the metabolites involved in histidine metabolism and β-alanine metabolism were strongly related to the differential proteins. In the correlation network of GSK/OP, seven proteins and eight metabolites within nodes included two correlations and six negative correlations (Fig. 7E). In the correlation network of OP/SHAM, two proteins and two metabolites within nodes included two negative correlations (Fig. 7F).

Table 4 Correlation statistics of significant differential proteins and metabolites between GSK/OVX and OVX/SHAM
Fig. 7
figure 7

Matrix heatmap showing the correlation coefficient matrix of significantly different proteins expressing significantly different metabolites (A) (C). r was identified between −1 and + 1. The correlation coefficient (r) for proteins and metabolites is indicated using colors. While r > 0 indicates a positive correlation and is represented in red, r < 0 indicates a negative correlation and is represented in blue. Darker colors indicate stronger correlations. Hierarchical clustering heatmap showing Pearson’s correlation analysis of differential proteins and metabolites (B) (D). Each row indicates a significantly different metabolite, whereas each column indicates a significantly different protein. The correlation analysis network between GSK/OP and OP/SHAM for significantly different proteins and metabolites (E) (F)

Integrated analysis of metabolic pathways involved in differential proteins and differential metabolites

According to the statistics of the common metabolic pathways after GSK treatment and castration, differential proteins and metabolites jointly regulate 11 metabolic pathways. These include purine metabolism, pyrimidine metabolism, histidine metabolism, beta-alanine metabolism, inflammatory mediator regulation of TRP channels, platelet activation, tyrosine metabolism, cancer pathways, phenylalanine metabolism, glutamatergic synapses, and gap junctions, among others (Table 5). Among them, purine metabolism is correlated with histidine metabolism, beta-alanine metabolism is corelated with pyrimidine metabolism, and tyrosine metabolism is correlated with phenylalanine metabolism (Fig. 8).

Table 5 Differential proteins and metabolites enriched by the common metabolic pathways in GSK/OVX and OVX/SHAM. * represents a protein or metabolite with differential expression trends between GSK/OP and OP/SHAM
Fig. 8
figure 8

Figure showing common metabolic pathways of differential proteins and metabolites between GSK/OP and OP/SHAM. The boxes contain differentially expressed proteins. Red represents the upregulated expression of differential proteins or metabolites. Green shows the downregulated expression of proteins or metabolites. However, yellow depicts metabolic pathways involved in both differential metabolites and proteins

Protein levels of Slc1a3, Aldh3b1, Cant1, Gstz1, Bid decreased after treatment with GSK

Based on the combined analysis of proteomics and metabonomics, we carried out western blot to verify the expression changes of five differential proteins. There were differences in levels of Slc1a3, Aldh3b1, Cant1, Gstz1, Bid among the groups (Fig. 9A,B). The expression levels of Slc1a3, Aldh3b1, Cant1, Bid in the OP group were increased significantly (P < 0.05 or P < 0.01) compared with those in the SHAM, and decreased significantly after 12 weeks of GSK treatment (GSK vs OP) (P < 0.05 or P < 0.01), which were consistent with the results of integrated analysis. In addition, GSK almost downregulated Gstz1, but the results showed no statistical significance. (P = 0.061).

Fig. 9
figure 9

Expression of differential proteins in lumbar samples from the three groups. A Western blotting analyses of protein levels in lumbar in the SHAM, OP and GSK group. The image has been cropped and the original images can be found in supplementary material (Fig. 9C, Fig. 9D). The samples derived from the same experiment and blots were processed in parallel. B Protein levels in the lumbar of different rat groups. **p < 0.01 vs. the sham group; *p < 0.05 vs. the sham group; ## p < 0.01 vs. the OP group; # p < 0.05 vs. the OP group


With the increasing incidence of osteoporosis each year, the disease has become a worldwide public health and social problem [26]. Exploring the mechanism of GSK and other TCM can provide greater advantages and theoretical supplementation for TCM treatment and prevention of osteoporosis. Clearly, proteins and small molecule metabolites are closely associated with the occurrence of diseases. Thus, the combination of proteomics and metabonomics deserves our matched attention.

In this study, micro scanning CT data of bone trabecula showed that GSK capsules increased bone density and improved bone microarchitecture in OP rats. The specific molecular mechanism of GSK in the prevention and treatment of OP has not yet to be explained. Based on the combined analysis of proteomics and metabonomics, the results showed that the expression trends of differential proteins, i.e., Cant1, Gstz1, Aldh3b1, Bid, and Slc1a3, in the common metabolic pathway of differential proteins and metabolites between GSK/OP and OP/SHAM were different. Importantly, they were corrected in the GSK group, and we verified them through experiments. Cant1 is involved in purine and pyrimidine metabolism. Aldh3b1 is involved in the metabolism of tyrosine, histidine, beta-alanine, and phenylalanine. After GSK treatment, Slc1a3 was strongly correlated with L-histidine (P = 0.045), and Aldh3b1 was significantly correlated with 1-stearoyl-sn-glycerol 3-phosphocholine (P = 0.045). Combined with proteomic analysis, Gstz1 and Aldh3b1 were located in the functional nodes of the protein–protein interaction network. Additionally, the expression trend of 12 metabolites (tyramine, thymidine, deoxycytidine, cytosine, L-aspartate, ergothioneine, L-histidine, L-anserine, pantothenate, arachidonic acid (peroxide-free), benzoic acid, and DL-phenylalanine) differed in the common enrichment metabolic pathway between GSK/OP and OP/SHAM. The above metabolic pathways are primarily involved in nucleotide metabolism, amino acid metabolism, the immune system, cell processes, and other major signaling pathways, indicating that GSK may treat osteoporosis in rats through these.

Uric acid (UA) is the final product of purine metabolism in the human body. It is mainly decomposed from nucleic acids by enzymes or from other purine compounds metabolized by cells. UA is a reducing substance in the human body, participating in redox reactions and scavenging oxygen free radicals, thereby inhibiting oxidative stress [27, 28]. Several studies have reported that serum UA is involved in the pathogenesis of OP by affecting oxidative stress and inflammatory cascades [29]. Higher serum UA levels appear to be protective from bone loss in peri- and postmenopausal women [30]. Cant1, a calcium-activated nucleotide, is essential for glycosaminoglycan synthesis in cartilage [31]. In this study, Cant1, the key differential protein, was involved in purine metabolism, which downregulated after GSK treatment.

Amino acid metabolism is crucial for life activities. As the basic unit of macromolecular proteins, amino acids are essential for the normal metabolism of the human body [32]. Previous studies showed that amino acid metabolism was important to osteoporosis [33, 34]. In the correlation analysis, the significant differential proteins and metabolites indicate a strong correlation with histidine metabolism and β-alanine metabolism. During childhood, histidine is an essential amino acid in the human body. In adulthood, histidine can be synthesized by the human body, thereby making it nonessential. Histamine is formed under the action of histidine decarboxylase and is a factor that regulates inflammation and allergic reactions [35]. An increase in the amount of bone cortex and bone minerals was found in mice deficient in histamine synthesis, revealing that a lack of histamine may increase the bone density and bone formation of histamine-deficient mice [36].

The results showed that after treatment with GSK, Aldh3b1 in the histidine metabolism pathway was downregulated. L-aspartate was also downregulated. Ergothioneine, L-histidine, and L-anserine were upregulated. β-alanine plays a key role in bone metabolism, mainly by improving the production of insulin and insulin-like growth factor-1, as well as the synthesis of collagen and muscle proteins. Studies have shown that β-alanine may be a potential biomarker of osteoporosis, as β-alanine levels increase significantly in postmenopausal osteoporosis patients with low bone density [37]. Tyrosine is an important amino acid in the body and is the main raw material for the synthesis of thyroxine [38]. Previous studies have suggested that GSK may participate in bone metabolism by regulating tyrosine metabolism and differential proteins. It is reported that tyrosine metabolism is significantly related to osteoporosis [39]. In addition, Aldh3b1 is involved in the metabolism of tyrosine, histidine, beta-alanine, and phenylalanine. We speculate that Aldh3b1 may be closely related to the metabolism of amino acids. Aldh3b1 is a member of the ALDH family. Its function is related to the scavenging of reactive oxygen species [40]. In postmenopausal women, the level of estrogen and antioxidants decreased with age. Meanwhile, the ROS that accumulated in the body could not be cleared in time, leading to oxidative stress and the damage of osteoblasts and osteocytes [41, 42]. ALDH can reduce oxidative stress through a variety of pathways involved in aldehyde metabolism [43].

In this study, the expression pattern of arachidonic acid (AA) in the platelet activation pathway and inflammatory mediator regulation of TRP channels changed after GSK treatment. Platelets are known to play an important role in inflammation and thrombosis [44]. Higher platelet counts are significantly associated with osteopenia and osteoporosis [45]. AA and its products play important roles in regulating the inflammatory response [46], vascular elasticity, platelet activation, and bone remodeling. AA belongs to the polyunsaturated fatty acid (PUFA) family, one of the indispensable free fatty acids in the human body [47]. According to the position of the first double bond, PUFAs can be divided into n-3 PUFAs and n-6 PUFAs. n-6 PUFAs mainly include linoleic acid and AA. Recent research highlights the potential role of PUFA in the bone protection [48,49,50]. The combined use of n-3 PUFA and E2 exerted synergistic bone-protective efficacy through the upregulation of RUNX2, an essential transcription factor for bone formation, as well as the suppression of bone-resorbing cytokine IL-1β [51].


This study reports the significant changes of proteins and metabolites in osteoporotic rats treated with GSK capsules. Integrated proteomic and metabonomic analyses showed that the GSK capsules may regulate differential proteins and differential metabolites to participate in nucleotide metabolism, amino acid metabolism, immune processes, and general cellular processes. They affect bone metabolism via the above mentioned pathways and play a role in bone protection. This study provides a reference for the molecular mechanism of the GSK capsules in the treatment of osteoporosis. Follow-up studies are also needed to further explore the differential proteins and metabolites in these metabolic pathways in order to find new drug targets and mechanisms of action.

Availability of data and materials

The data analyzed during the study are available from the corresponding author on reasonable request.







Traditional Chinese medicine




Variable Importance for the Projection


Information dependent acquisition


Partial least squares discriminant analysis


Orthogonal partial least squares


Bone mineral density


Principal component analysis


Uric acid


Arachidonic acid


Polyunsaturated fatty acids


  1. Armas LA, Recker RR. Pathophysiology of osteoporosis: new mechanistic insights. Endocrinol Metab Clin N Am. 2012;41(3):475–86.

    Article  CAS  Google Scholar 

  2. Men Z, Huang C, Xu M, Ma J, Wan L, Huang J, et al. Zhuanggu Zhitong Capsule alleviates postmenopausal osteoporosis in ovariectomized rats by regulating autophagy through AMPK/mTOR signaling pathway. Ann Transl Med. 2022;10(16):900.

    Article  CAS  Google Scholar 

  3. Jin R, Lin ZJ, Xue CM, Zhang B. An improved association-mining research for exploring Chinese herbal property theory: based on data of the Shennong's Classic of Materia Medica. J Integrative Med. 2013;11(5):352–65.

    Article  Google Scholar 

  4. Mukwaya E, Xu F, Wong MS, Zhang Y. Chinese herbal medicine for bone health. Pharm Biol. 2014;52(9):1223–8.

    Article  Google Scholar 

  5. Zhang ND, Han T, Huang BK, Rahman K, Jiang YP, Xu HT, et al. Traditional Chinese medicine formulas for the treatment of osteoporosis: Implication for antiosteoporotic drug discovery. J Ethnopharmacol. 2016;189:61–80.

    Article  Google Scholar 

  6. Naiqiang Z, Jinxin L, Guiyun M, Jingyi H. Simultaneous Determination of Six Ingredients in Gushukang Capsule by HPLC-DAD. Chin J Exp Tradit Med Formulae. 2018;19:42–6.

    Article  Google Scholar 

  7. Shi S, Wang F, Huang Y, Chen B, Pei C, Huang D, et al. Epimedium for Osteoporosis Based on Western and Eastern Medicine: An Updated Systematic Review and Meta-Analysis. Front Pharmacol. 2022;13:782096.

    Article  CAS  Google Scholar 

  8. Li S, Zhou H, Hu C, Yang J, Ye J, Zhou Y, et al. Corrigendum: Total Flavonoids of Rhizoma Drynariae Promotes Differentiation of Osteoblasts and Growth of Bone Graft in Induced Membrane Partly by Activating Wnt/β-Catenin Signaling Pathway. Front Pharmacol. 2021;12:726831.

    Article  Google Scholar 

  9. Ou L, Kang W, Liang Z, Gao F, Dong T, Wei P, et al. Investigation of anti-osteoporosis mechanisms of Rehmanniae Radix Preparata based on network pharmacology and experimental verification. J Orthop Surg Res. 2021;16(1):599.

    Article  Google Scholar 

  10. Wu XH, Dou B, Sun NY, Gao J, Liu XL. Astragalus saponin IV promotes osteogenic differentiation of bone marrow mesenchymal stem cells via miR-21/NGF/BMP2/Runx2 pathway. Acta Histochem. 2020;122(4):151549.

    Article  CAS  Google Scholar 

  11. Xie YM, Yuwen Y, Dong FH, Sun SC, Wang HM, Liu QS, et al. Clinical practice guideline of traditional medicine for primary osteoporosis. Chin J Integrative Med. 2011;17(1):52–63.

    Article  Google Scholar 

  12. Wang G, Huo L, Chen G, He H. Effects of Gushukang for postmenopausal osteoporosis: A protocol for systematic review and meta-analysis. Medicine. 2020;99(27):e20908.

    Article  Google Scholar 

  13. Wang Q, Zhao Y, Sha N, Zhang Y, Li C, Zhang H, et al. The systemic bone protective effects of Gushukang granules in ovariectomized mice by inhibiting osteoclastogenesis and stimulating osteoblastogenesis. J Pharmacol Sci. 2018;136(3):155–64.

    Article  CAS  Google Scholar 

  14. Li XL, Wang L, Bi XL, Chen BB, Zhang Y. Gushukang exerts osteopreserve effects by regulating vitamin D and calcium metabolism in ovariectomized mice. J Bone Miner Metab. 2019;37(2):224–34.

    Article  CAS  Google Scholar 

  15. Chai S, Wan L, Wang JL, Huang JC, Huang HX. Gushukang inhibits osteocyte apoptosis and enhances BMP-2/Smads signaling pathway in ovariectomized rats. Phytomedicine. 2019;64:153063.

    Article  CAS  Google Scholar 

  16. Peng Y, Wu S, Li Y, Crane JL. Type H blood vessels in bone modeling and remodeling. Theranostics. 2020;10(1):426–36.

    Article  CAS  Google Scholar 

  17. Li W, Zhou X, Jiang T, He H, Wen T. Positive Effect of Gushukang on Type-H Vessel and Bone Formation. Front Cell Dev Biol. 2020;8:265.

    Article  Google Scholar 

  18. Palmieri O, Mazza T, Castellana S, Panza A, Latiano T, Corritore G, et al. Inflammatory Bowel Disease Meets Systems Biology: A Multi-Omics Challenge and Frontier. Omics. 2016;20(12):692–8.

    Article  CAS  Google Scholar 

  19. Sun YV, Hu YJ. Integrative Analysis of Multi-omics Data for Discovery and Functional Studies of Complex Human Diseases. Adv Genet. 2016;93:147–90.

    Article  CAS  Google Scholar 

  20. Shao WH, Chen JJ, Fan SH, Lei Y, Xu HB, Zhou J, et al. Combined Metabolomics and Proteomics Analysis of Major Depression in an Animal Model: Perturbed Energy Metabolism in the Chronic Mild Stressed Rat Cerebellum. Omics. 2015;19(7):383–92.

    Article  CAS  Google Scholar 

  21. Thompson DD, Simmons HA, Pirie CM, Ke HZ. FDA Guidelines and animal models for osteoporosis. Bone. 1995;17(4 Suppl):125s–33s.

    Article  CAS  Google Scholar 

  22. Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359–62.

    Article  CAS  Google Scholar 

  23. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  Google Scholar 

  24. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  Google Scholar 

  25. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  Google Scholar 

  26. Sietsema DL. Fighting the Epidemic: Bone Health and Osteoporosis. Nurs Clin North Am. 2020;55(2):193–202.

    Article  Google Scholar 

  27. Pasalic D, Marinkovic N, Feher-Turkovic L. Uric acid as one of the important factors in multifactorial disorders--facts and controversies. Biochem Med (Zagreb). 2012;22(1):63–75.

    Article  CAS  Google Scholar 

  28. Cheng S, Yang Y, Zhou Y, Xiang W, Yao H, Ma L. Influence of different concentrations of uric acid on oxidative stress in steatosis hepatocytes. Exp Ther Med. 2018;15(4):3659–65.

    Article  CAS  Google Scholar 

  29. Kim SC, Paik JM, Liu J, Curhan GC, Solomon DH. Gout and the Risk of Non-vertebral Fracture. J Bone Mineral Res. 2017;32(2):230–6.

    Article  CAS  Google Scholar 

  30. Makovey J, Macara M, Chen JS, Hayward CS, March L, Seibel MJ, et al. Serum uric acid plays a protective role for bone loss in peri- and postmenopausal women: a longitudinal study. Bone. 2013;52(1):400–6.

    Article  CAS  Google Scholar 

  31. Kodama K, Takahashi H, Oiji N, Nakano K, Okamura T, Niimi K, et al. CANT1 deficiency in a mouse model of Desbuquois dysplasia impairs glycosaminoglycan synthesis and chondrocyte differentiation in growth plate cartilage. FEBS Open Bio. 2020;10(6):1096–103.

    Article  CAS  Google Scholar 

  32. Wu G. Amino acids: metabolism, functions, and nutrition. Amino Acids. 2009;37(1):1–17.

    Article  CAS  Google Scholar 

  33. Torricelli P, Fini M, Giavaresi G, Giardino R. Human osteopenic bone-derived osteoblasts: essential amino acids treatment effects. Artif Cells Blood Substit Immobil Biotechnol. 2003;31(1):35–46.

    Article  CAS  Google Scholar 

  34. Jennings A, MacGregor A, Spector T, Cassidy A. Amino Acid Intakes Are Associated With Bone Mineral Density and Prevalence of Low Bone Mass in Women: Evidence From Discordant Monozygotic Twins. J Bone Mineral Res. 2016;31(2):326–35.

    Article  CAS  Google Scholar 

  35. Gibbs NK. l-Histidine Supplementation in Adults and Young Children with Atopic Dermatitis (Eczema). J Nutr. 2020;150(Suppl 1):2576s–9s.

    Article  Google Scholar 

  36. Fitzpatrick LA, Buzas E, Gagne TJ, Nagy A, Horvath C, Ferencz V, et al. Targeted deletion of histidine decarboxylase gene in mice increases bone formation and protects against ovariectomy-induced bone loss. Proc Natl Acad Sci U S A. 2003;100(10):6027–32.

    Article  CAS  Google Scholar 

  37. Yu L, Qi H, An G, Bao J, Ma B, Zhu J, et al. Association between metabolic profiles in urine and bone mineral density of pre- and postmenopausal Chinese women. Menopause (New York, NY). 2019;26(1):94–102.

    Article  Google Scholar 

  38. den Hartog MT, Sijmons CC, Bakker O, Ris-Stalpers C, de Vijlder JJ. Importance of the content and localization of tyrosine residues for thyroxine formation within the N-terminal part of human thyroglobulin. Eur J Endocrinol. 1995;132(5):611–7.

    Article  Google Scholar 

  39. Ling CW, Miao Z, Xiao ML, Zhou H, Jiang Z, Fu Y, et al. The Association of Gut Microbiota With Osteoporosis Is Mediated by Amino Acid Metabolism: Multiomics in a Large Cohort. J Clin Endocrinol Metab. 2021;106(10):e3852–64.

    Article  Google Scholar 

  40. Marchitti SA, Brocker C, Orlicky DJ, Vasiliou V. Molecular characterization, expression analysis, and role of ALDH3B1 in the cellular protection against oxidative stress. Free Radic Biol Med. 2010;49(9):1432–43.

    Article  CAS  Google Scholar 

  41. Li M, Hao L, Li L, Liu L, Chen G, Jiang W, et al. Cinnamtannin B-1 Prevents Ovariectomy-Induced Osteoporosis via Attenuating Osteoclastogenesis and ROS Generation. Front Pharmacol. 2020;11:1023.

    Article  CAS  Google Scholar 

  42. Hyeon S, Lee H, Yang Y, Jeong W. Nrf2 deficiency induces oxidative stress and promotes RANKL-induced osteoclast differentiation. Free Radic Biol Med. 2013;65:789–99.

    Article  CAS  Google Scholar 

  43. Li H, Toth E, Cherrington NJ. Alcohol Metabolism in the Progression of Human Nonalcoholic Steatohepatitis. Toxicol Sci. 2018;164(2):428–38.

    Article  CAS  Google Scholar 

  44. Koupenova M, Clancy L, Corkrey HA, Freedman JE. Circulating Platelets as Mediators of Immunity, Inflammation, and Thrombosis. Circ Res. 2018;122(2):337–51.

    Article  CAS  Google Scholar 

  45. Kim J, Kim HS, Lee HS, Kwon YJ. The relationship between platelet count and bone mineral density: results from two independent population-based studies. Arch Osteoporos. 2020;15(1):43.

    Article  Google Scholar 

  46. Calder PC. Polyunsaturated fatty acids and inflammatory processes: New twists in an old tale. Biochimie. 2009;91(6):791–5.

    Article  CAS  Google Scholar 

  47. Wang L, Zhang C, Liang H, Zhou N, Huang T, Zhao Z, et al. Polyunsaturated Fatty Acids Level and Bone Mineral Density: A Two-Sample Mendelian Randomization Study. Front Endocrinol (Lausanne). 2022;13:858851.

    Article  Google Scholar 

  48. Thanabalan A, Kiarie EG. Influence of Feeding Omega-3 Polyunsaturated Fatty Acids to Broiler Breeders on Indices of Immunocompetence, Gastrointestinal, and Skeletal Development in Broiler Chickens. Front Vet Sci. 2021;8:653152.

    Article  Google Scholar 

  49. Cugno C, Kizhakayil D, Calzone R, Rahman SM, Halade GV, Rahman MM. Omega-3 fatty acid-rich fish oil supplementation prevents rosiglitazone-induced osteopenia in aging C57BL/6 mice and in vitro studies. Sci Rep. 2021;11(1):10364.

    Article  CAS  Google Scholar 

  50. Longo AB, Ward WE. PUFAs, Bone Mineral Density, and Fragility Fracture: Findings from Human Studies. Adv Nutr. 2016;7(2):299–312.

    Article  CAS  Google Scholar 

  51. Jin Y, Lee M, Park Y. Synergistic attenuation of ovariectomy-induced bone loss by combined use of fish oil and 17β-oestradiol. Br J Nutr. 2017;117(4):479–89.

    Article  CAS  Google Scholar 

Download references


We gratefully appreciate all laboratory personnel for their experimental guidance and Shanghai Applied Protein Technology Co. (Shanghai, China) for their assistance with bioinformatics analysis techniques.


The research was supported by the National Natural Sciences Foundations of China (NO. 81774350) and Basic research projects of provincial scientific research institutions for public welfare in Fujian province (NO. 2019R1003–1).

Author information

Authors and Affiliations



RL wrote the first draft of the manuscript and did the formal analysis. BX and LX investigated the experiments or contributed reagents and tools. SL designed the research project and was responsible for funding acquisition. JG administered and supervised the project. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shengqiang Li.

Ethics declarations

Ethics approval and consent to participate

The presented research work involves the use of SD rats and all experimental protocols were complied with the the revised Animals (Scientific Procedures) Act 1986 and the ARRIVE (Animal Research: Reporting of In Vivo Experiments) guidelines, and were approved by the Animal Ethics Committee of Fujian Academy of Traditional Chinese Medicine (Fujian, China) [SYXK(Fujian)2009–000] (FJATCM-IAEC2018034). All animals purchase from Shanghai Slake Laboratory Animal Co., Ltd. (Shanghai, China) [SCXK (Shanghai) 2007–0005], and were raised in a SPF environment in the Experimental Center of Comparative Medicine, Fujian Academy of Chinese Medical Sciences. All methods were performed in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information


Additional file 1.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lin, R., Xie, B., Xie, L. et al. Integrated proteomics and metabolomics analysis of lumbar in a rat model of osteoporosis treated with Gushukang capsules. BMC Complement Med Ther 22, 333 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Proteomics
  • Metabonomics
  • Gushukang capsules
  • Osteoporosis
  • Integration
  • Bioinformatics