Exploring the Effect of Jiawei Buguzhi Pills on TGF-β-Smad Pathway in Postmenopausal Osteoporosis Based on Integrated Pharmacological Strategy

Abstract

Objective. To explore the effect of Jiawei Buguzhi Pills (JWBGZP) on the TGF-β-Smad pathway in postmenopausal osteoporosis (PMO) based on integrated pharmacological strategy. Method. The ETCM database was used to collect JWBGZP. GeneCards and OMIM databases were utilized to obtain PMO-related genes. Cytoscape was used for network construction and analysis, and DAVID was used for GO and KEGG enrichment analysis of key targets. Animal experiments and cell experiments were conducted to further explore the mechanism. The bone mass density was detected by dual-energy X-ray bone densitometer. The TGF-β1 and Smad4 mRNA in bone tissue were detected by RT-qPCR. The TGF-β1 and Smad4 protein in bone tissue were detected by the western blot. The TGF-β1 and Smad4 protein in osteoblasts were determined by immunohistochemistry. Result. A total of 721 JWBGZP potential targets and 385 PMO-related genes were obtained. The enrichment analysis showed that JWBGZP may regulate the TGF-beta signaling pathway, oxidation-reduction process, aging, response to hypoxia, response to ethanol, negative regulation of cell proliferation, PI3K-Akt, HIF-1, and other signaling pathways. The animal experiments showed that compared with the model group, the femoral bone mineral density and lumbar bone mineral density of the JWBGZP group increased (P < 0.05); the expression levels of TGF-β1 and Smad mRNA and proteins in the JWBGZP group were significantly higher (P < 0.05). The cell experiment results showed a large number of osteoblast stained blue-purple and orange-red calcified nodules. The expression levels of TGF-β1 and Smad proteins in the JWBGZP group were significantly higher than those in the blank control group and the sham operation group, and the protein expression levels in the model group were the lowest (P < 0.05). Conclusion. JWBGZP may be involved in PI3K-Akt, HIF-1, estrogen, prolactin, and other signaling pathways and regulate MAPK1, AKT1, PIK3CA, JAK2, and other gene targets, regulate bone metabolism, and thereby treat PMO.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Postmenopausal osteoporosis (PMO) is a metabolic bone disease that usually occurs 3–4 years after menopause. The sharp drop in estrogen levels after menopause leads to bone resorption greater than bone formation, resulting in bone loss. The patient usually has no obvious symptoms, but due to the loss of bone mass, it is prone to fragility fractures [1, 2]. Clinically, PMO often presents lower back pain or leg pain, which is chronic and persistent and can be aggravated by coughing and bending over [3]. In severe cases, there may be a series of symptoms such as hunchback or shortening of height due to vertebral compression fractures. Epidemiological studies have shown that China has over 210 million people over 60 years old, and there are approximately 140 million people aged 65 and over [4]. In China, among the elderly, 60–69 years old, osteoporosis in women accounts for about 50% to 70%, and about 30% in men [5]. By 2050, medical expenses for fractures in Europe will reach approximately 7.5 billion euros [6]. A series of complications such as bone deformation and fracture caused by PMO not only affect the physical and mental health and quality of life of postmenopausal women, but also cause a heavy family and social burden [6]. In postmenopausal women, due to the decrease in estrogen secretion, oxidative stress response occurs in the body. The large amount of active oxygen not only inhibits bone formation, but also promotes bone resorption, which eventually causes bone loss [7]. At present, the treatment and prevention drugs for PMO mainly include calcium, vitamin D, antibone resorption drugs, and bone active drugs [6, 8–10]. Although they play an important role in the treatment of osteoporosis, there are also problems such as drug side effects and insufficient patient compliance. On the other hand, Chinese medicine (TCM) is based on the whole and seeks the root of the disease. TCM is gradually playing more roles in the treatment of osteoporosis [11, 12].

Buguzhi pills from “Compendium of Materia Medica” is a classic prescription for strengthening the kidney and strengthening the bones. Current studies have shown that Psoralea corylifolia Linn has an estrogen-like effect, maintains bone density, inhibits bone resorption, and promotes bone formation in patients with PMO, and may improve the body's bone metabolism balance [13–16]. Jiawei Buguzhi pills (JWBGZP) are modified from Buguzhi Pills, in order to further explore the potential mechanism of JWBGZP in the treatment of PMO. Network pharmacology is based on the theory of system biology and the network analysis of biological systems [17, 18]. Network pharmacology emphasizes the multichannel regulation of signal pathways because it is consistent with the characteristics of the overall concept of Chinese medicine (TCM). It provides an effective research approach for predicting the core targets and pathways of drug treatment of diseases and is especially suitable for the research of complex systems such as TCM or chronic diseases [19,20]. At present, through the application of system biology strategies in the analysis of TCM compound prescriptions, not only has the TCM drug compatibility law been verified at the system level, but it has also played an important role in finding potential active compounds [21, 22]. This study uses network pharmacology methods to explore the potential effective components and targets of JWBGZP in the treatment of PMO and analyze the key targets of gene ontology functions and signal pathways to provide references for future anti-PMO treatments and in-depth research.

2. Materials and Methods

2.1. JWBGZP Potential Targets and PMO Gene Collection

JWBGZP consists of Fructus Psoraleae (Bu Gu Zhi) 15 g, Epimedii Folium (Yin Yang Huo) 15 g, Lycii Fructus (Gou Qi) 10 g, Cyathulae Radix (Chuan Niu Xi) 10 g, Fructus Ligustri Lucidi (Nv Zhen Zi) 10 g, Rhizoma Cibotii (Gou Ji) 10 g, and Rhizoma Drynariae (Gu Sui Bu) 10 g. The herbs of JWBGZP were used as a keyword to search in ETCM (http://www.tcmip.cn/ETCM/index.php/Home/Index/) [23] (Table S1). “Postmenopausal osteoporosis” was used as a search term to search for PMO-related genes in the GeneCards [24] and OMIM databases [25], respectively. The search results of the two data are summarized, and duplicate genes are removed, which is the target of the PMO disease (Table S2).

2.2. Network Construction and Analysis

JWBGZP potential targets and PMO genes were input into STRING database (https://string-db.Org/), organization is set to “Homo sapiens” to obtain protein-protein interaction (PPI) data. The PPI data is imported into Cytoscape to construct the network and perform topology analysis. JWBGZP potential targets and PMO genes were imported into the DAVID 6.8 database (https://david.ncifcrf.gov/) [26], and the gene ontology (GO) function enrichment analysis and KEGG pathway enrichment analysis were performed, respectively. Among them, GO analysis includes three contents: biological process (BP), molecular function (MF), and cellular component (CC).

2.3. Experimental Materials

2.3.1. Experimental Animals

(1) Animal for animal experiments: 75 SD female rats, body weight (200 ± 10) g, were provided by Guangdong Medical Laboratory Animal Center (license number: SCXK (Guangdong) 2018-0002).

(2) Animal for cell experiments: 60 SD female rats, body weight (200 ± 10) g, were provided by Guangdong Medical Laboratory Animal Center (license number: SCXK (Guangdong) 2018-0002). 10 SD rats born within 72 hours were provided by the Experimental Animal Center of Hunan University of Chinese Medicine (Certificate Number: SCXK (Xiang) 2019-0004). The experiment was approved by the Animal Ethics Committee of Hunan University of Chinese Medicine (Ethics approval number: LLBH-202103050004).

2.3.2. Instruments and Reagents

JWBGZP extract (prepared by the Department of Natural Medicine, Central South University) is a mixture of oil phase and water phase. Psoralen (batch number: 110739-201817, purity: 99.7%) and isopsoralen (batch number: 110738-201814, purity: 99.9%) were purchased from China Institute for Food and Drug Control.

DEME low-sugar solution, penicillin-streptomycin mixed solution, and type II collagenase (Shanghai Lifei Biotechnology Company) were used. The serum E2 kit was purchased from Wuhan Boster Biology Co. Ltd. Fetal bovine serum was purchased from Zhejiang Tianhang Biological Technology Co. Ltd; alkaline phosphatase staining solution and Alizarin red staining solution were purchased from Progema Biotechnology Company. TGF-β1 ELISA kit was purchased from Blue Gene Biotech Co. (lot number: EK3812/2); TGF-β1-anti (Santa Cruz Biotechnology, lot number: sc-130348). HRP-labeled goat anti-rabbit IG was purchased from Abcam, UK). RIPA protein lysate and protease inhibitor were purchased from Beyotime Biotechnology (product number P0013B (P1030)). Goat anti-rabbit IgG was purchased from Zhongshan Jinqiao (batch number: ZB-2301). β-Actin antibody (anti-beta actin antibody [mAbcam 8226]-loading control (ab8226)) was used. PCR primers were provided by Shanghai Bioengineering Co. Ltd. Animal tissue total RNA extraction kit (TianGEN, lot number: Q5301) and 2 × SYBR Premix ExTaqTM (TianGEN, lot number: Q6220) were used. 1260 type high-performance liquid chromatography (HPLC), including G1329B autosampler, G1315D diode array detector, G1311C quaternary gradient elution pump, G1316A column thermostat, and 1260 chromatography workstation (Agilent, USA) was used; low-temperature centrifuge, 37°C thermostat, DNP-9162 electric heating thermostat incubator, and LEICA DM LB2 binocular microscope were used; dual energy X-ray bone densitometer (Unigamma X-ray Plus), Motic B microscopic camera system, and MIAS medical image analysis system were also used.

2.4. JWBGZP Quality Control by HPLC

2.4.1. Solution Preparation

Reference solution: an appropriate amount of psoralen and isopsoralen reference substance was accurately weighed and placed in the same 50 mL measuring flask, to which methanol-0.1% acetic acid solution (1 : 1, V/V) was added to dissolve and dilute to volume. Then, 1 ml of it was taken and placed in a 10 ml measuring flask, and methanol-0.1% acetic acid solution (1 : 1, V/V) was added to make the volume constant. JWBGZP solution was prepared from the JWBGZP extract.

2.4.2. HPLC Condition

HPLC conditions were as follows: column: Kinetex-C1 (250 mm × 4.6 mm, 5 um); mobile phase: methanol (A)-0.1% acetic acid solution (B), gradient elution 1; flow rate: 1.0 mL/min; detection wavelength: 248 nm; column temperature: 30°C; and injection volume: 20 uL. The results are shown in Figure 1.

Table 1

The primer sequences.

Gene	Sequence forward	Sequence reverse
TGF-β1	5′-CCAAGGAGACGGAATACAGG3′	5′-GTGTTGGTTGTAGAGGGCAAG-3′
Smad4	5′-TCAGCCAGCTACTTACCACCA-3′	5′-ACACGTCCCCTTCACCTTTAC-3′
β-Actin	5′-TGTATGCCTCTGGTCGTACCAC-3′	5′-ACAGAGTACTTGCGCTCAGGAG3′

[figures omitted; refer to PDF]

2.5. Experimental Methods for Animal Experiments

2.5.1. Animal Grouping, Modeling, and Intervention

The female rats were randomly divided into sham operation group, model group, JWBGZP low dose, medium dose, and high dose groups with 15 rats each. Modeling was carried out according to the literature method [3]. The rats were anesthetized with 2% sodium pentobarbital (0.2 mL/100 g), and the ovaries were harvested by making approximately 1 cm incisions on both sides of the lumbar spine. In the sham operation group, only the same amount of fat was cut, and the abdominal wall was sutured in layers and sterilized with ethanol. 800,000 U of penicillin was given daily within 3 days after treatment. After the operation, the rats received vaginal smears every day to observe the changes in the estrus cycle. Rats with estrous cycle disorders are considered to be successful in ovariectomy.

The bone density was measured by using the dual-energy X-ray absorption method, and the rats with bone density reduction were selected as the successful model of postmenopausal osteoporosis. 12 weeks after modeling, rats were given intragastric administration. The dosage is calculated based on the body surface area of humans and rats. The high dose of JWBGZP is equivalent to 13.2 g/kg of crude drug, the medium dose is equivalent to 6.6 g/kg of crude drug, and the low dose is equivalent to 3.3 g/kg of the crude drug. The administration volume was 1 mL/100g, and the model group and the sham operation group were given equal volumes of purified water. The administration lasted for 13 weeks.

2.5.2. Serum E2 Level Detection

After 13 weeks of administration, rats in each group were anesthetized by intraperitoneal injection of 2% sodium pentobarbital (0.2 mL/100 g), blood was collected from the abdominal aorta, and serum was separated. Serum E2 was measured by using radioimmunoassay.

2.5.3. Femur Tissue mRNA Expression Detected by RT-qPCR

The femoral tissue was ground in liquid nitrogen, and 1 ml TRlzol was added to extract total RNA from the tissue. NanoDrop 2000 spectrophotometer was used to detect the purity and concentration of total RNA; PrimeScriptMRT kit and gDNA Eraser were used to reverse transcribe RNA into cDNA. Amplification conditions were as follows: 95°C predenaturation for 5 min: 95°C denaturation for 15 s, 60°C extension for 30 s, repeat 40 cycles. Finally, the expression levels of TGF-β1, Smad4, and β-actin mRNA were calculated according to the 2-△△ method. The RT-PCR primer sequences are shown in Table 1.

2.5.4. Femur Tissue Protein Expression Detected by Western Blot

50 mg of femoral bone tissue was extracted by RIPA protein lysate according to the manufacturer's instructions, and the protein was quantified using the BCA method. After the separation of the proteins by SDS-PAGE electrophoresis, the proteins were transferred to a polyvinylidene (PVDF) membrane and blocked with 5% skimmed milk powder for 60 min. Afterwards, the PVDF membrane was incubated with primary antibodies (TGF-β1 Smad4 and β-actin, diluted 1 : 1000) at 4°C overnight, and the secondary antibody (HRP-labeled goat anti-IG, diluted 1 : 10000) was incubated for 90 min. Finally, the ECL reagent was used for chemiluminescence development, and Image J was used to analyze the gray value.

2.6. Experimental Methods for Cell Experiments

2.6.1. Animal Grouping and Modeling

The female rats were randomly divided into blank control group, sham operation group, model group, and JWBGZP group with 15 rats each. Modeling and model assessment were carried out according to Section 2.5.1.

2.6.2. Extraction and Culture of Osteoblasts

Ten SD rats born within 72 hours were anesthetized and sacrificed by cervical dislocation. They were immersed in absolute ethanol for 20 s, and the skulls were cut off aseptically in PBS buffer and rinsed. Then, it was cut into 1 mm × 1 mm bone pieces and placed in 4 mL 0.25% trypsin for digestion at 37°C for 30 min. The bone fragments were transferred into a centrifuge tube containing 0.1% type I collagenase, digested in a constant temperature water bath at 37°C for 1 h, and centrifuged at 1,000 r/min for 10 min twice. Then, they were added into the culture medium in aliquots, placed in a volume fraction of 5% CO₂, and incubated at 37°C. The medium was changed for the first time after 3 days, and the medium was changed every 3 days thereafter. Then, it was digested with 0.25% trypsin, the cells were passaged in a 1 : 1 manner, and the 3rd generation cells were used for subsequent experiments.

2.6.3. Drug-Containing Serum Preparation

Intervention began in the 3rd week after modeling. The JWBGZP group was given 1.5 mL of JWBGZP every day, and the other groups were given 1.5 mL of normal saline for 12 consecutive weeks. The dose of the drug was corrected once a week according to body weight. 24 hours after the last administration, 2% sodium pentobarbital was injected intraperitoneally, and 8 mL of blood was taken from the main abdominal artery. The blood sample was placed at 4°C for 4 h, centrifuged at 3,000 r/min for 30 min, the supernatant was taken, and inactivated at 56°C for 30 min. The sera of the same group were mixed, filtered through a 0.22 μL sterile filter and aliquoted, and stored at −20°C.

2.6.4. Drug-Containing Serum Intervention

The third-generation osteoblasts were digested with 0.25% trypsin, and the cell suspension was inoculated into a 24-well culture plate at a final concentration of 2 × 10⁴/mL, divided into 4 × 4 wells. After 24 hours of cell attachment, 1 mL of their respective serum was added, and the medium was changed for 3 days during the culture. The culture supernatant was collected when the medium was changed for the first time and stored at −20°C for later use for immunohistochemical detection.

2.7. Statistical Analysis

The SPSS 20.0 software was used for statistical analysis, and the measurement data were expressed as “mean ± SD”. The comparison was performed by the t-test, and P < 0.05 was considered statistically significant.

3. Results and Discussion

3.1. JWBGZP Potential Targets and PMO-Related Genes

Bu Gu Zhi is warm and pungent; it returns to the kidney meridian. Yin Yang Huo is warm, sweet, mildly bitter, and pungent; it returns to the liver meridian, triple burner meridian, and kidney meridian. Chuan Niu Xi is even, sweet, bitter, and sour; it returns to the kidney meridian and liver meridian. Gou Qi Zi is even and sweet; it returns to kidney meridian and liver meridian. Nv Zhen Zi is cool, sweet, and bitter; it returns to kidney meridian and liver meridian. Gu Sui Bu is warm and bitter; it returns to kidney meridian and liver meridian (Figure 2(a)). A total of 150 JWBGZP components, 721 JWBGZP potential targets, and 385 PMO-related genes were obtained. The JWBGZP herbs, components, and potential targets are shown in Figures 2(b) and 2(c). The JWBGZP potential target, PMO-related gene, and JWBGZP-PMO PPI data were input into Cytoscape 3.8.2 to construct the JWBGZP-PMO PPI network. This network consists of 610 JWBGZP target nodes, 285 PMO gene nodes, 55 JWBGZP-PMO target nodes, and 20658 edges (Figure 3). The top 30 targets were divided into 3 categories: (1) JWBGZP target set: GAPDH (350 edges), DECR1 (184 edges), CASP3 (179 edges), HRAS (178 edges), and CREB1 (167 edges); (2) PMO gene set: ALB (305 edges), IL6 (296 edges), MYC (226 edges), EGF (217 edges), STAT3 (203 edges), IGF1 (200 edges), JUN (199 edges), FOS (192 edges), CXCL8 (190 edges), CTNNB1 (187 edges), IL1B (181 edges), IL10 (177 edges), MMP9 (175 edges), PTEN (174 edges), CAT (173 edges), and LEP (166 edges); (3) JWBGZP-PMO target set: INS (355 edges), AKT1 (327 edges), TP53 (291 edges), VEGFA (256 edges), MAPK3 (254 edges), TNF (249 edges), SRC (213 edges), ESR1 (187 edges), and PTGS2 (170 edges).

[figures omitted; refer to PDF]

[figure omitted; refer to PDF]

Through PPI target network analysis, the important targets of JWBGZP for the treatment of PMO are obtained: MAPK1, AKT1, PIK3CA, PIK3R1, MAPK8, JAK2, etc. Among them, MAPK1 and MAKE8 are members of the MAPK signaling pathway, which can preclude osteoblast apoptosis through the mitochondrial pathway [27]. MAPK1, also known as ERK2, can mediate osteoblast differentiation through the ERK1/2 pathway [28, 29]. AKT1 is a member of the gene family of AKT (protein kinase B). Studies have found that after the AKT1 gene knockout treatment in rats, the rats have restricted bone development. It shows that AKTI may play a role in promoting bone growth in bone development [30], and at the same time, AKT1, PIK3CA, and PIK3R1 jointly act as members of the PI3K-Akt signaling pathway, and play a regulatory role in the process of bone metabolism through this pathway [31, 32]. JAK2 is a member of the JAK/STAT signal pathway, which can cause bone destruction by participating in the inflammatory response [33]. The TGF-β1/Smads signaling pathway is an important signaling pathway between bone resorption and bone formation in the process of bone metabolism, and it plays an important role in the balance regulation of bone metabolism [34]. The smad protein family is the only transduction molecule that conducts signals in TGF-β1 cells. Smads can combine with other transcription molecules to transfer the response signal from the cell membrane to the nucleus. Among them, Smad4 is a common molecule required for TGF-β1 signal transduction [35].

The results of GO enrichment analysis can infer that these key targets mainly play a role in the plasma membrane, membrane raft, and cell nucleus. JWBGZP may exert molecular functions such as kinase activity, protein binding, and protein kinase binding through positive regulation of nitric oxide biosynthesis, signal transduction, protein phosphorylation, and other biological processes. The PI3K-Akt signaling pathway is an important signaling pathway in the process of bone metabolism. It has the effects of mediating osteoblast proliferation [36, 37], osteoblast differentiation [38], and interfering with osteoclast apoptosis [39]. It is worth mentioning that the key targets AKT1, PIK3CA, and PIK3R1 selected in this study are all members of this pathway, and the research results are mutually verified. The HIF-1 signaling pathway has previously been proven to promote bone formation [40]. In animal experiments, it was found that the HIF-1 signaling pathway was inhibited, and the degree of osteoporosis in experimental mice was improved. While the osteoclast index of bone metabolism decreased, the index of bone metabolism osteogenesis also decreased [40]. The HIF-1 signaling pathway may be mainly involved in the regulation of osteoclast activity and exert its anti-PMO effect. This research also involves three hormonal signaling pathways. Among them, the estrogen signaling pathway promotes the proliferation and differentiation of osteoblasts by up-regulating the expression of OPG/RANKL through estrogen receptors [41], and at the same time, up-regulating CFTR to induce the apoptosis of osteoclasts, which has a two-way regulation effect on bone formation and osteoclast. The prolactin signaling pathway acts on the reproductive axis to inhibit the synthesis of steroid hormones and mediate bone metabolism [42]. Thyroxine directly or indirectly acts on chondrocytes, osteoclasts, and osteoblasts to achieve bone transformation [43].

In summary, the treatment of PMO by JWBGZP is to act on multiple gene targets through multiple active ingredients in JWBGZP. The functions of each gene are complex and diverse and regulate bone metabolism through different signaling pathways. This reflects the characteristics of multiple components, multiple targets, and multiple pathways in the treatment of PMO by JWBGZP and provides a reference for in-depth research on the mechanism of JWBGZP in the treatment of PMO in the future.

3.3. Serum E2 and Bone Density of Ovariectomized Female Rats

After castration, the serum E2 of the rats decreased significantly, which was statistically significant compared with the nonexperimental group (P < 0.01). After treatment with JWBGZP, E2 content increased significantly (P < 0.05) (Figure 8). Compared with the sham operation group, the femoral bone mineral density and lumbar bone mineral density of the model group were significantly lower (P < 0.05), indicating that the model was successful. Compared with the model group, the femoral bone mineral density and lumbar bone mineral density of the JWBGZP group increased (P < 0.05) (Figure 9).

[figure omitted; refer to PDF]

The TGF-β-Smad signaling transduction pathway consists of the TGF-β signal molecule and the Smad protein. The TGF-β signal molecule conducts signal transduction through the transmembrane receptor complex [44]. There are mainly two types of TGF-β receptors (TβR), type I and type II, which are necessary for its conduction process [45], and its functions involve extracellular matrix deposition, growth inhibition, and apoptosis. The TGF-β gene exerts antiosteoporosis by regulating the function of bone deposition and osteoclasts, and the TGF-β protein expressed by it promotes the proliferation and differentiation of osteoblasts and promotes the synthesis of extracellular matrix to affect the formation of bone [46]. In addition, TGF-β also has an active bone repair function, and studies have found that it is highly correlated with the concentration of estrogen [47]. Smad signaling protein is a transcription factor that plays an important role in the intracellular signal transduction of members of the TGF-β superfamily. All TGF-β signaling molecules must be mediated by Smad4 protein to enter the nucleus, and Smad6 and 7 are inhibitors of this pathway [48, 49]. This study found that under the action of JWBGZP, the expression level of the TGF-β-Smad signal transduction pathway in osteoblasts was significantly higher than that of other groups. This proves that TGF-β-Smad is an important pathway for JWBGZP to treat postmenopausal osteoporosis, and JWBGZP can increase the expression level of this pathway.

This research explored the molecular mechanism of JWBGZP in the treatment of osteoporosis through network pharmacology and found that JWBGZP may improve the pathological process of osteoporosis through TGF-β-related signaling pathways through in vivo and in vitro experiments. However, this study still has shortcomings. For example, in this study, the exploration of the downstream signal pathway transduction of the TGF-β signal pathway is still not deep enough. In the future, we hope to design inhibitors and siRNAs targeting the TGF-β signaling pathway to further study in vivo and in vitro osteoporosis models.

4. Conclusion

JWBGZP may be involved in PI3K-Akt, HIF-1, estrogen, prolactin, and other signaling pathways, regulate MAPK1, AKT1, PIK3CA, JAK2, and other gene targets, regulate bone metabolism, and thereby treat PMO.

Disclosure

Xiao Yuan, Yonghe Wu and Kailin Yang should be considered joint first author.

Authors’ Contributions

Kailin Yang, Xiao Yuan, Yonghe Wu, Huiping Liu, and Guomin Zhang participated in the concept and design. Kailin Yang, Xiao Yuan, and Yonghe Wu were responsible for data analysis and interpretation in the network pharmacological section; Huiping Liu and Guomin Zhang were responsible for data analysis and interpretation of experiments; Kailin Yang, and Yonghe Wu drafted the paper; Guomin Zhang and Huiping Liu supervised the study; all authors participated in the analysis and interpretation of data and approved the final paper.

References

[1] W. Mysiw, R. Jackson, "Insights into the epidemiology of postmenopausal osteoporosis: the women's health initiative," Seminars in Reproductive Medicine, vol. 32 no. 06, pp. 454-462, DOI: 10.1055/s-0034-1384629, 2014.

[2] L. A. G. Armas, R. R. Recker, "Pathophysiology of osteoporosis," Endocrinology and Metabolism Clinics of North America, vol. 41 no. 3, pp. 475-486, DOI: 10.1016/j.ecl.2012.04.006, 2012.

[3] J. M. Lane, L. Russell, S. N. Khan, "Osteoporosis," Clinical Orthopaedics and Related Research, vol. 372 no. 372, pp. 139-150, DOI: 10.1097/00003086-200003000-00016, 2000.

[4] Q. Zeng, N. Li, Q. Wang, J. Feng, D. Sun, Q. Zhang, J. Huang, Q. Wen, R. Hu, L. Wang, Y. Ma, X. Fu, S. Dong, X. Cheng, "The prevalence of osteoporosis in China, a nationwide, multicenter DXA survey," Journal of Bone and Mineral Research, vol. 34 no. 10, pp. 1789-1797, DOI: 10.1002/jbmr.3757, 2019 Oct.

[5] Osteoporosis and Bone Mineral Salt Disease Branch of Chinese Medical Association, "Osteoporosis in China: an epidemiological survey," Chinese Journal of Osteoporosis and Bone Mineral Salt Diseases, vol. 12 no. 04, pp. 317-318, 2019. (in chinese)

[6] J. A. Kanis, C. Cooper, C. Cooper, R. Rizzoli, J.-Y. Reginster, "European guidance for the diagnosis and management of osteoporosis in postmenopausal women," Osteoporosis International, vol. 30 no. 1,DOI: 10.1007/s00198-018-4704-5, 2019.

[7] R. Rizzoli, H. Bischoff-Ferrari, B. Dawson-Hughes, C. Weaver, "Nutrition and bone health in women after the menopause," Women's Health, vol. 10 no. 6, pp. 599-608, DOI: 10.2217/whe.14.40, 2014.

[8] J. A. Kanis, C. Cooper, C. Cooper, R. Rizzoli, J.-Y. Reginster, "Executive summary of the European guidance for the diagnosis and management of osteoporosis in postmenopausal women," Calcified Tissue International, vol. 104 no. 3, pp. 235-238, DOI: 10.1007/s00223-018-00512-x, 2019.

[9] D. Gatti, A. Fassio, "Pharmacological management of osteoporosis in postmenopausal women: the current state of the art," Journal of Population Therapeutics and Clinical Pharmacology, vol. 26 no. 4, pp. e1-e17, DOI: 10.15586/jptcp.v26.i4.646, 2019.

[10] L. Cui, T. He, Y. Jiang, M. Li, O. Wang, R. Jiajue, Y. Chi, Q. Xu, X. Xing, W. Xia, "Predicting the intervention threshold for initiating osteoporosis treatment among postmenopausal women in China: a cost-effectiveness analysis based on real-world data," Osteoporosis International, vol. 31 no. 2, pp. 307-316, DOI: 10.1007/s00198-019-05173-6, 2020.

[11] T. Wang, Q. Liu, W. Tjhioe, J. Zhao, A. Lu, G. Zhang, R. X. Tan, M. Zhou, J. Xu, H. T. Feng, "Therapeutic potential and outlook of alternative medicine for osteoporosis," Current Drug Targets, vol. 18 no. 9, pp. 1051-1068, DOI: 10.2174/1389450118666170321105425, 2017.

[12] N.-D. Zhang, T. Han, B.-K. Huang, K. Rahman, Y.-P. Jiang, H.-T. Xu, L.-P. Qin, H.-L. Xin, Q.-Y. Zhang, Y.-m. Li, "Traditional Chinese medicine formulas for the treatment of osteoporosis: implication for antiosteoporotic drug discovery," Journal of Ethnopharmacology, vol. 189, pp. 61-80, DOI: 10.1016/j.jep.2016.05.025, 2016.

[13] K. Li, J.-jin Wang, N. Zhou, Li-jie Du, Z.-ling Zhang, "Discussion on the rationality of Buguzhi pill formulation and salt preparation based on efficacy evaluation," Chinese Traditional and Herbal Medicine, vol. 50 no. 22, pp. 5508-5514, 2019. (in chinese)

[14] W. D. Li, C. P. Yan, Y. Wu, Z. B. Weng, F. Z. Yin, G. M. Yang, B. C. Cai, Z. P. Chen, "Osteoblasts proliferation and differentiation stimulating activities of the main components of Fructus Psoraleae corylifoliae," Phytomedicine, vol. 21 no. 4, pp. 400-405, DOI: 10.1016/j.phymed.2013.09.015, 2014.

[15] X.-y. Cai, Z.-j. Zhang, J.-l. Xiong, M. Yang, Z.-t. Wang, "Experimental and molecular docking studies of estrogen-like and anti-osteoporosis activity of compounds in Fructus Psoraleae," Journal of Ethnopharmacology, vol. 276,DOI: 10.1016/j.jep.2021.114044, 2021.

[16] A. X. D. Yu, M. L. Xu, P. Yao, K. K. L. Kwan, Y. X. Liu, R. Duan, T. T. X. Dong, R. K. M. Ko, K. W. K. Tsim, "Corylin, a flavonoid derived from Psoralea Fructus, induces osteoblastic differentiation via estrogen and Wnt/ β ‐catenin signaling pathways," The FASEB Journal, vol. 34 no. 3, pp. 4311-4328, DOI: 10.1096/fj.201902319RRR, 2020.

[17] L. Zeng, K. Yang, "Exploring the pharmacological mechanism of Yanghe decoction on HER2-positive breast cancer by a network pharmacology approach," Journal of Ethnopharmacology, vol. 199, pp. 68-85, DOI: 10.1016/j.jep.2017.01.045, 2017.

[18] D. Zhao, Y. Yi, Q. He, S. Wang, K. Yang, J. Ge, "Exploring the regulatory mechanism of Nao Tai Fang on vascular Dementia's biological network based on cheminformatics and transcriptomics strategy," Journal of Ethnopharmacology, vol. 274,DOI: 10.1016/j.jep.2021.114065, 2021.

[19] K. Yang, L. Zeng, A. Ge, Y. Shi, X. Zhu, W. Liu, J. Ge, "The effect of hedysarum multijugum maxim.-chuanxiong rhizoma compound on ischemic stroke: a research based on network and experimental pharmacology," Oxidative Medicine and Cellular Longevity, vol. 2020,DOI: 10.1155/2020/6072380, 2020.

[20] K. Yang, L. Zeng, A. Ge, Y. Yi, S. Wang, J. Ge, "Exploring the oxidative stress mechanism of buyang huanwu decoction in intervention of vascular dementia based on systems biology strategy," Oxidative Medicine and Cellular Longevity, vol. 2021,DOI: 10.1155/2021/8879060, 2021.

[21] L. Zeng, K. Yang, L. Liu, T. Zhang, H. Liu, Z. Tan, L. Lei, "Systematic biological and proteomics strategies to explore the regulation mechanism of Shoutai Wan on recurrent spontaneous Abortion's biological network," Journal of Ethnopharmacology, vol. 263,DOI: 10.1016/j.jep.2020.113156, 2020.

[22] K. Yang, L. Zeng, A. Ge, T. Bao, T. Xu, X. Xie, L. Liu, "Exploring the regulation mechanism of xihuang pill, olibanum and β -boswellic acid on the biomolecular network of triple-negative breast cancer based on transcriptomics and chemical informatics methodology," Frontiers in Pharmacology, vol. 11,DOI: 10.3389/fphar.2020.00825, 2020.

[23] H.-Y. Xu, Y.-Q. Zhang, Z.-M. Liu, T. Chen, C.-Y. Lv, S.-H. Tang, X.-B. Zhang, W. Zhang, Z.-Y. Li, R.-R. Zhou, H.-J. Yang, X.-J. Wang, L.-Q. Huang, "ETCM: an encyclopaedia of traditional Chinese medicine," Nucleic Acids Research, vol. 47 no. D1, pp. D976-D982, DOI: 10.1093/nar/gky987, 2018.

[24] G. Stelzer, R. Rosen, I. Plaschkes, S. Zimmerman, M. Twik, S. Fishilevich, T. Iny Stein, R. Nudel, I. Lieder, Y. Mazor, S. Kaplan, D. Dahary, D. Warshawsky, Y. Guan - Golan, A. Kohn, N. Rappaport, M. Safran, D. Lancet, "The GeneCards suite: from gene data mining to disease genome sequence analysis," Current Protocols in Bioinformatics, vol. 54, pp. 1301-1303, DOI: 10.1002/cpbi.5, 2016.

[25] A. Hamosh, A. F. Scott, J. S. Amberger, C. A. Bocchini, V. A. McKusick, "“Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders," Nucleic Acids Research, vol. 33 no. Database issue, pp. D514-D517, DOI: 10.1093/nar/gki033, 2005.

[26] D. W. Huang, B. T. Sherman, R. A. Lempicki, "Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources," Nature Protocols, vol. 4 no. 1, pp. 44-57, DOI: 10.1038/nprot.2008.211, 2009.

[27] Weng Wenyu Hao Ye Hua Zhang Yusheng, "The role of mitochondrial pathway in apoptosis of osteoblasts," Chinese Journal of Osteoporosis, vol. 24 no. 12, pp. 1665-1670, 2018.

[28] Z. Y. Cheng, T. Ye, Q. Y. Ling, T. Wu, G. Y. Wu, G. J. Zong, "Parathyroid hormone promotes osteoblastic differentiation of endothelial cells via the extracellular signal‑regulated protein kinase 1/2 and nuclear factor‑ κ B signaling pathways," Experimental and Therapeutic Medicine, vol. 15 no. 2, pp. 1754-1760, DOI: 10.3892/etm.2017.5545, 2018.

[29] T. Noguchi, K. Ebina, M. Hirao, S. Otsuru, A. J. Guess, R. Kawase, T. Ohama, S. Yamashita, Y. Etani, G. Okamura, H. Yoshikawa, "Apolipoprotein E plays crucial roles in maintaining bone mass by promoting osteoblast differentiation via ERK1/2 pathway and by suppressing osteoclast differentiation via c-Fos, NFATc1, and NF- κ B pathway," Biochemical and Biophysical Research Communications, vol. 503 no. 2, pp. 644-650, DOI: 10.1016/j.bbrc.2018.06.055, 2018.

[30] X.-d. Peng, P. Z. Xu, M. L. Chen, A. Hahn-Windgassen, J. Skeen, J. Jacobs, D. Sundararajan, W. S. Chen, S. E. Crawford, K. G. Coleman, N. Hay, "Dwarfism, impaired skin development, skeletal muscle atrophy, delayed bone development, and impeded adipogenesis in mice lacking Akt1 and Akt2," Genes & Development, vol. 17 no. >11, pp. 1352-1365, DOI: 10.1101/gad.1089403, 2003.

[31] J. Xi, Q. Li, X. Luo, Y. Wang, J. Li, L. Guo, G. Wu, "Celastrol inhibits glucocorticoid‑induced osteoporosis in rat via the PI3K/AKT and Wnt signaling pathways," Molecular Medicine Reports, vol. 18 no. 5, pp. 4753-4759, DOI: 10.3892/mmr.2018.9436, 2018.

[32] Y. Zhao, Y. Xu, H. Zheng, N. Lin, "QingYan formula extracts protect against postmenopausal osteoporosis in ovariectomized rat model via active ER-dependent MEK/ERK and PI3K/Akt signal pathways," Journal of Ethnopharmacology, vol. 268,DOI: 10.1016/j.jep.2020.113644, 2021.

[33] C.-h. Li, L.-l. Xu, L.-l. Jian, R.-h. Yu, J.-x. Zhao, L. Sun, G.-h. Du, X.-y. Liu, "Static inhibits RANKL-mediated osteoclastogenesis by suppressing activation of STAT3 and NF- κ B pathways," International Immunopharmacology, vol. 58, pp. 136-144, DOI: 10.1016/j.intimp.2018.03.021, 2018.

[34] X. Li, X. Sui, Q. Yang, Y. Li, N. Li, X. Shi, D. Han, Y. Li, X. Huang, P. Yu, X. Qu, "Oviductus Ranae protein hydrolyzate prevents menopausal osteoporosis by regulating TGF β /BMP2 signaling," Archives of Gynecology and Obstetrics, vol. 299 no. 3, pp. 873-882, DOI: 10.1007/s00404-018-5033-9, 2019 Mar.

[35] L. Wu, Y. Wu, B. Gathings, M. Wan, X. Li, W. Grizzle, Z. Liu, C. Lu, Z. Mao, X. Cao, "Smad4 as a transcription corepressor for estrogen receptor α," Journal of Biological Chemistry, vol. 278 no. 17, pp. 15192-15200, DOI: 10.1074/jbc.M212332200, 2003.

[36] D. Liang, M. Yang, B. Guo, J. Cao, L. Yang, X. Guo, Y. Li, Z. Gao, "Zinc inhibits H2O2-induced MC3T3-E1 cells apoptosis via MAPK and PI3K/AKT pathways," Biological Trace Element Research, vol. 148 no. 3, pp. 420-429, DOI: 10.1007/s12011-012-9387-8, 2012 Sep.

[37] J.-C. Xi, H.-Y. Zang, L.-X. Guo, H.-B. Xue, X.-D. Liu, Y.-B. Bai, Y.-Z. Ma, "The PI3K/AKT cell signaling pathway is involved in regulation of osteoporosis," Journal of Receptors and Signal Transduction, vol. 35 no. 6, pp. 640-645, DOI: 10.3109/10799893.2015.1041647, 2015.

[38] Y.-X. Gu, J. Du, M.-S. Si, J.-J. Mo, S.-C. Qiao, H.-C. Lai, "The roles of PI3K/Akt signaling pathway in regulating MC3T3-E1 preosteoblast proliferation and differentiation on SLA and SLActive titanium surfaces," Journal of Biomedical Materials Research Part A, vol. 101A no. 3, pp. 748-754, DOI: 10.1002/jbm.a.34377, 2013 Mar.

[39] D. Yu, D. Zhao, X. Yao, "Effects of PI3-K/Akt signaling pathway on osteoclast apoptosis," Chinese Journal of Osteoporosis, vol. 24 no. 7, pp. 874-878, 2018.

[40] J. Xu, Y. Sun, T. Wu, Y. Liu, L. Shi, J. Zhang, Q. Kang, Y. Chai, G. Li, "Enhancement of bone regeneration with the accordion technique via HIF‐1 α /VEGF activation in a rat distraction osteogenesis model," Journal of Tissue Engineering and Regenerative Medicine, vol. 12 no. 2, pp. e1268-e1276, DOI: 10.1002/term.2534, 2018 Feb.

[41] L. J. Sun, C. Li, X. H. Wen, L. Guo, Z. F. Guo, L. Q. Liao, Y. Guo, "Icariin stimulates hFOB 1.19 osteoblast proliferation and differentiation via OPG/RANKL mediated by the estrogen receptor," Current Pharmaceutical Biotechnology, vol. 22 no. 1, pp. 168-175, DOI: 10.2174/1389201021666200123102550.PMID:31971108, 2021.

[42] V. Bernard, J. Young, N. Binart, "Prolactin - a pleiotropic factor in health and disease," Nature Reviews Endocrinology, vol. 15 no. 6, pp. 356-365, DOI: 10.1038/s41574-019-0194-6, 2019.

[43] F. Mirza, E. Canalis, "Management of endocrine disease: secondary osteoporosis: pathophysiology and management," European Journal of Endocrinology, vol. 173 no. 3, pp. R131-R151, DOI: 10.1530/EJE-15-0118, 2015.

[44] J. Yu, L. Xu, K. Li, N. Xie, Y. Xi, Y. Wang, X. Zheng, X. Chen, M. Wang, X. Ye, "Zinc-modified calcium silicate coatings promote osteogenic differentiation through TGF- β /smad pathway and osseointegration in osteopenic rabbits," Scientific Reports, vol. 7 no. 1,DOI: 10.1038/s41598-017-03661-5, 2017.

[45] B. T. Hu, W. Z. Chen, "MOTS-c improves osteoporosis by promoting osteogenic differentiation of bone marrow mesenchymal stem cells via TGF- β /Smad pathway," European Review for Medical and Pharmacological Sciences, vol. 22 no. 21, pp. 7156-7163, DOI: 10.26355/eurrev_201811_16247, 2018.

[46] X. X. Li, N. Li, C. J. Ban, M. Zhu, B. Xiao, H. P. Dai, "Idiopathic pulmonary fibrosis in relation to gene polymorphisms of transforming growth factor- β 1 and plasminogen activator inhibitor 1," Chinese Medical Journal, vol. 124 no. 13, pp. 1923-1927, 2011.

[47] M. A. Estai, F. Suhaimi, S. Das, A. N. Shuid, Z. Mohamed, I.-N. Soelaiman, "Expression of TGF- β 1 in the blood during fracture repair in an estrogen-deficient rat model," Clinics, vol. 66 no. 12, pp. 2113-2119, DOI: 10.1590/s1807-59322011001200018, 2011.

[48] A. Moustakas, S. Souchelnytskyi, C. H. Heldin, "Smad regulation in TGF-beta signal transduction," Journal of Cell Science, vol. 114 no. Pt 24, pp. 4359-4369, DOI: 10.1242/jcs.114.24.4359, 2001.

[49] A. J. McCarthy, R. Chetty, "Smad4/DPC4," Journal of Clinical Pathology, vol. 71 no. 8, pp. 661-664, DOI: 10.1136/jclinpath-2018-205095, 2018.