- Open Access
Integrated proteomics and metabolomics analysis of lumbar in a rat model of osteoporosis treated with Gushukang capsules
BMC Complementary Medicine and Therapies volume 22, Article number: 333 (2022)
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.
Osteoporosis (OP) is a systemic metabolic disease characterized by low bone mass, microstructural damage of bone tissue, decreased bone strength, and increased bone fragility . Its clinical manifestations include primarily pain and spinal curvature . For thousands of years, traditional Chinese medicine (TCM) has been used in both China and other Asian countries . 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 . It can also regulate the overall function of the human body and relieve back pain and lumbago .
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 . 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 . Rhizoma Drynariae can promote osteogenesis differentiation by activating the Wnt/β-catenin signaling pathways, which play a role in antiosteoporosis activity . 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 . 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 . 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 , and it is widely used in the clinical settings . 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 . It also significantly increases the production of vitamin D and calcium . 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 . Furthermore, H-type blood vessels have shown the ability to induce angiogenesis and bone formation . Recent studies demonstrated that GSK can enhance hypoxia-inducible factor-1α to induce H-type vascular formation as well as bone formation .
In recent years, multiomics integrated analysis technology has been widely used to study various disease states . 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 . The combined application of proteomic and metabolomic technology has undeniable advantages . 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.
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 . Intramuscular injection of penicillin at 800,000 units/day was used to prevent infection for 3 days after the ovariectomy.
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.
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 , 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] (http://home.jp/kegg/).
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 (1∶1000 dilution), anti-Aldh3b1 (1∶1000), anti-Cant1 (1∶1000), anti-Gstz1(1∶1000), anti-Bid (1: 1000) and anti-β-actin (1∶4000) 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).
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.
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.
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.
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).
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).
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).
With the increasing incidence of osteoporosis each year, the disease has become a worldwide public health and social problem . 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 . Higher serum UA levels appear to be protective from bone loss in peri- and postmenopausal women . Cant1, a calcium-activated nucleotide, is essential for glycosaminoglycan synthesis in cartilage . 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 . 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 . 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 .
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 . Tyrosine is an important amino acid in the body and is the main raw material for the synthesis of thyroxine . 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 . 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 . 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 .
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 . Higher platelet counts are significantly associated with osteopenia and osteoporosis . AA and its products play important roles in regulating the inflammatory response , 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 . 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β .
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
Polyunsaturated fatty acids
Armas LA, Recker RR. Pathophysiology of osteoporosis: new mechanistic insights. Endocrinol Metab Clin N Am. 2012;41(3):475–86. https://doi.org/10.1016/j.ecl.2012.04.006.
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. https://doi.org/10.21037/atm-22-3724.
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. https://doi.org/10.3736/jintegrmed2013051.
Mukwaya E, Xu F, Wong MS, Zhang Y. Chinese herbal medicine for bone health. Pharm Biol. 2014;52(9):1223–8. https://doi.org/10.3109/13880209.2014.884606.
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. https://doi.org/10.1016/j.jep.2016.05.025.
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. https://doi.org/10.13422/j.cnki.syfjx.20181420.
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. https://doi.org/10.3389/fphar.2022.782096.
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. https://doi.org/10.3389/fphar.2021.726831.
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. https://doi.org/10.1186/s13018-021-02751-5.
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. https://doi.org/10.1016/j.acthis.2020.151549.
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. https://doi.org/10.1007/s11655-011-0613-6.
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. https://doi.org/10.1097/md.0000000000020908.
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. https://doi.org/10.1016/j.jphs.2018.01.007.
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. https://doi.org/10.1007/s00774-018-0924-1.
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. https://doi.org/10.1016/j.phymed.2019.153063.
Peng Y, Wu S, Li Y, Crane JL. Type H blood vessels in bone modeling and remodeling. Theranostics. 2020;10(1):426–36. https://doi.org/10.7150/thno.34126.
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. https://doi.org/10.3389/fcell.2020.00265.
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. https://doi.org/10.1089/omi.2016.0147.
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. https://doi.org/10.1016/bs.adgen.2015.11.004.
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. https://doi.org/10.1089/omi.2014.0164.
Thompson DD, Simmons HA, Pirie CM, Ke HZ. FDA Guidelines and animal models for osteoporosis. Bone. 1995;17(4 Suppl):125s–33s. https://doi.org/10.1016/8756-3282(95)00285-l.
Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359–62. https://doi.org/10.1038/nmeth.1322.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51. https://doi.org/10.1002/pro.3715.
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. https://doi.org/10.1093/nar/gkaa970.
Sietsema DL. Fighting the Epidemic: Bone Health and Osteoporosis. Nurs Clin North Am. 2020;55(2):193–202. https://doi.org/10.1016/j.cnur.2020.02.002.
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. https://doi.org/10.11613/bm.2012.007.
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. https://doi.org/10.3892/etm.2018.5855.
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. https://doi.org/10.1002/jbmr.2978.
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. https://doi.org/10.1016/j.bone.2012.10.025.
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. https://doi.org/10.1002/2211-5463.12859.
Wu G. Amino acids: metabolism, functions, and nutrition. Amino Acids. 2009;37(1):1–17. https://doi.org/10.1007/s00726-009-0269-0.
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. https://doi.org/10.1081/bio-120018002.
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. https://doi.org/10.1002/jbmr.2703.
Gibbs NK. l-Histidine Supplementation in Adults and Young Children with Atopic Dermatitis (Eczema). J Nutr. 2020;150(Suppl 1):2576s–9s. https://doi.org/10.1093/jn/nxaa200.
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. https://doi.org/10.1073/pnas.0934373100.
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. https://doi.org/10.1097/gme.0000000000001158.
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. https://doi.org/10.1530/eje.0.1320611.
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. https://doi.org/10.1210/clinem/dgab492.
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. https://doi.org/10.1016/j.freeradbiomed.2010.08.004.
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. https://doi.org/10.3389/fphar.2020.01023.
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. https://doi.org/10.1016/j.freeradbiomed.2013.08.005.
Li H, Toth E, Cherrington NJ. Alcohol Metabolism in the Progression of Human Nonalcoholic Steatohepatitis. Toxicol Sci. 2018;164(2):428–38. https://doi.org/10.1093/toxsci/kfy106.
Koupenova M, Clancy L, Corkrey HA, Freedman JE. Circulating Platelets as Mediators of Immunity, Inflammation, and Thrombosis. Circ Res. 2018;122(2):337–51. https://doi.org/10.1161/CIRCRESAHA.117.310795.
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. https://doi.org/10.1007/s11657-020-0700-4.
Calder PC. Polyunsaturated fatty acids and inflammatory processes: New twists in an old tale. Biochimie. 2009;91(6):791–5. https://doi.org/10.1016/j.biochi.2009.01.008.
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. https://doi.org/10.3389/fendo.2022.858851.
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. https://doi.org/10.3389/fvets.2021.653152.
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. https://doi.org/10.1038/s41598-021-89827-8.
Longo AB, Ward WE. PUFAs, Bone Mineral Density, and Fragility Fracture: Findings from Human Studies. Adv Nutr. 2016;7(2):299–312. https://doi.org/10.3945/an.115.009472.
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. https://doi.org/10.1017/s0007114517000344.
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).
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
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.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
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). https://doi.org/10.1186/s12906-022-03807-7
- Gushukang capsules