Tao Li1†
Tingting Tan
Cheng Xiang
Zhengxiao Ouyang- 1Department of Orthopedics, The Second Xiangya Hospital, Central South University, Changsha, Hunan, China
- 2Department of Immunology, School of Basic Medical Science, Central South University, Changsha, Hunan, China
Background: Osteoporosis is a prevalent bone metabolism disease characterized by a reduction in bone density, leading to several complications that significantly affect patients’ quality of life. The Achyranthes bidentata–Dipsacus asper (AB–DA) herb pair is commonly used in Traditional Chinese Medicine (TCM) to treat osteoporosis. This study aimed to investigate the therapeutic compounds and potential mechanisms of AB–DA using network pharmacology, molecular docking, molecular dynamics simulation, and experimental verification.
Methods: Identified compounds of AB–DA were collected from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), Traditional Chinese Medicine Information Database (TCM-ID), TCM@Taiwan Database, BATMAN-TCM, and relevant literature. The main bioactive ingredients were screened based on the criteria of “OB (oral bioavailability) ≥ 30, DL (drug-likeness) ≥ 0.18.” Potential targets were predicted using the PharmMapper and SwissTargetPrediction websites, while disease (osteoporosis)-related targets were obtained from the GeneCards, DisGeNET, and OMIM databases. The PPI network and KEGG/GO enrichment analysis were utilized for core targets and pathway screening in the STRING and Metascape databases, respectively. A drug–compound–target–pathway–disease network was constructed using Cytoscape software to display core regulatory mechanisms. Molecular docking and dynamics simulation techniques explored the binding reliability and stability between core compounds and targets. In vitro and in vivo validation experiments were utilized to explore the anti-osteoporosis efficiency and mechanism of sitogluside.
Results: A total of 31 compounds with 83 potential targets for AB–DA against osteoporosis were obtained. The PPI analysis revealed several hub targets, including AKT1, CASP3, EGFR, IGF1, MAPK1, MAPK8, and MAPK14. GO/KEGG analysis indicated that the MAPK cascade (ERK/JNK/p38) is the main pathway involved in treating osteoporosis. The D–C–T–P–T network demonstrated therapeutic compounds that mainly consisted of iridoids, steroids, and flavonoids, such as sitogluside, loganic acid, and β-ecdysterone. Molecular docking and dynamics simulation analyses confirmed strong binding affinity and stability between core compounds and targets. Additionally, the validation experiments showed preliminary evidence of antiosteoporosis effects.
Conclusion: This study identified iridoids, steroids, and flavonoids as the main therapeutic compounds of AB–DA in treating osteoporosis. The underlying mechanisms may involve targeting core MAPK cascade (ERK/JNK/p38) targets, such as MAPK1, MAPK8, and MAPK14. In vivo experiments preliminarily validated the anti-osteoporosis effect of sitogluside. Further in-depth experimental studies are required to validate the therapeutic value of AB–DA for treating osteoporosis in clinical practice.
1 Introduction
Osteoporosis is a prevalent bone metabolism disease that affects over 200 million people worldwide (McDonald et al., 2021; Grewe et al., 2022). It is characterized by a reduction in bone density, which greatly increases the risk of fractures (NIH Consensus Development Panel on Osteoporosis Prevention et al., 2001). Osteoporotic fractures and associated complications can have a significant and lasting impact on patients’ quality of life, sometimes even threatening their lives, and place a considerable cost on society and individuals (Compston et al., 2019; Liang et al., 2022). With the aging of the global population, osteoporosis has become a pressing health problem (Jiang et al., 2020).
Anti-osteoporosis drugs such as estrogen, raloxifene, bisphosphonates, calcitonin, and parathyroid hormone (PTH) are commonly used in clinical practice (Chen et al., 2022). These medications, however, have side effects and severe responses that restrict their long-term usage. For example, bisphosphonates may cause jaw osteonecrosis and renal impairment (Li et al., 2021). There is, therefore, an urgent need to identify potential anti-osteoporosis drugs that are both more effective and safer.
Traditional Chinese Medicine (TCM) has a rich history and has been widely used in Asia to treat various diseases, including osteoporosis (Huang et al., 2022). TCM is cost-effective and has fewer side effects than chemosynthetic drugs, making it more suitable for long-term use (Mukwaya et al., 2014). The application of TCM in modern society provides a new pathway for complementary and alternative medicine (CAM) treatment (Chu et al., 2022). Over the years, many TCM treatments and prescriptions have been used to treat various orthopedic diseases, especially osteoporosis and fractures, with great success (Suvarna et al., 2018; Peng et al., 2022). TCM’s prescription to treat osteoporosis can also play a comprehensive role in regulating body function and relieving pain (Feng et al., 2022). A Chinese herb pair, generally composed of two kinds of herbal medicine, is the essence of TCM prescriptions. Compared with all herbs in prescriptions, studying and elucidating the complex pharmacological mechanism of herb pairs is simpler and more beneficial (Liu et al., 2020).
Achyranthes bidentata (AB), also known as Niu Xi, is included in the Chinese, Japanese, and Korean pharmacopoeia (He et al., 2017). Additionally, its dried roots are regularly used in TCM for osteoporosis (Yan et al., 2019). A number of biological activities, including anti-osteoporosis (He et al., 2010; Zhang et al., 2012; Jiang et al., 2014; Suh et al., 2014; Zhang M. et al., 2018; Zhang S. et al., 2018), anti-tumor (Jin et al., 2007), and anti-oxidant (Huang et al., 2015), have been demonstrated by contemporary pharmacological research on AB extracts. Dipsacus asper (DA), also known in Chinese as Xu Duan —meaning “to rebuild fractures and unite bones”—is discussed in Shennong’s Classic of Material Medicine, which is the earliest source (Tao et al., 2020). DA can be used to treat muscle pain and bone repair, golden sores, and collapses (Wu et al., 2022). According to modern pharmacological investigations, numerous disorders, including osteoporosis and osteoarthritis, have been successfully treated with DA (Liu et al., 2019; Yu et al., 2019; Zhang et al., 2019). Jiegudan capsules, which contain AB and DA, are a common traditional Chinese medicine prescription to treat osteoporosis. Although many compounds have been isolated from AB and DA, the potential pharmacological mechanisms of AB–DA herb pairs and their interactions with osteoporosis-related targets and pathways remain unclear and need further exploration.
In recent years, the use of network pharmacology has become increasingly popular for exploring the interaction network of TCM therapy (Shuai et al., 2020). Molecular docking, a virtual screening technology that simulates the behavior of small-molecule ligands at the binding sites of receptor proteins, has also gained popularity for developing novel drugs (Pagadala et al., 2017). This research aims to elucidate the potential mechanism in TCM of the AB–DA herb pair for treating osteoporosis, bioinformatics prediction by network pharmacology, molecular docking, and molecular dynamics simulation, and verify these via alkaline phosphatase (ALP) activity, osteoblast mineralization assays, Western blot, q-PCR, and an ovariectomy (OVX) osteoporosis mouse model. A flow chart outlining the study’s approach is presented in Figure 1.
FIGURE 1. Study flow chart to investigate the potential underlying mechanisms for AB–DA treatment of osteoporosis.
2 Materials and methods
2.1 Screening of bioactive compounds of AB–DA
The TCM@Taiwan (http://tcm.cmu.edu.tw/zh-tw/) Database (Chen, 2011), Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP http://lsp.nwu.edu.cn/tcmsp.php) (Ru et al., 2014), BATMAN-TCM (L et al., 2016) (a bioinformatics analytical tool for the molecular mechanisms of TCM: http://bionet.ncpsb.org.cn/), Traditional Chinese Medicine Information Database (Kang et al., 2013) (TCM-ID http://bidd.group), and relevant literature were utilized to acquire all identified AB–DA compounds. ADME (absorption, distribution, metabolism, and excretion) properties were applied to screen bioactive ingredients, and the screening criteria were set as “oral bioavailability (OB)≥ 30, drug-likeness (DL)≥ 0.18” (Gao et al., 2022) to screen the compounds from the TCMSP database. Similarly, compounds from different sources were screened in the SwissADME (Daina et al., 2017) database using their pharmacokinetic properties (http://www.swissadme.ch).
2.2 Relevant targets of AB–DA compound and osteoporosis
To predict potential targets based on their spatial configuration, the compounds generated in the previous step were imported into the SwissTargetPrediction (Daina et al., 2019) (http://www.swisstargetprediction.ch/) and PharmMapper databases (Liu et al., 2010) (http://www.lilab-ecust.cn/pharmmapper/). The UniProt ID of the target was converted into a standardized gene name using the UniProt database (Holzhüter and Geertsma, 2022) (https://www.uniprot.org). The keyword “osteoporosis” was searched in the GeneCards (Barshir et al., 2021) (https://www.genecards.org/), DisGeNET (Piñero et al., 2020) (https://www.disgenet.org), and Online Mendelian Inheritance in Man (Li et al., 2012) (OMIM, https://omim.org) databases to obtain relevant targets. Then, the overlapping targets identified by Venn diagram were considered as targets of AB–DA for the treatment of osteoporosis after merging and removing duplicates.
2.3 PPI network for core target selection
To identify potential hub genes, a protein–protein interaction (PPI) network was established in the STRING (Szklarczyk et al., 2021) (http://string-db.org, Version 11.5) database, with a focus on the anti-osteoporosis efficacy of AB–DA in Homo sapiens and an interaction score threshold of 0.4. Topological analysis was performed, and the core targets of AB–DA for treating osteoporosis were accurately selected by using the Cytoscape plug-ins MCODE (molecular complex detection) and CytoHubba (Ye et al., 2022).
2.4 GO and KEGG pathway enrichment analyses for the core pathways
Analysis of Gene Ontology (GO) functions, including cellular component (CC), molecular function (MF), and biological process (BP), and Kyoto Encyclopedia of Genome and Genome (KEGG) pathway enrichment analysis, was utilized to clarify the key anti-osteoporosis mechanism of the AB–DA herb pair. When entering the targets into the Metascape database (http://www.metascape.org/), the cut-off p-value, minimum overlap value, and concentration value were set to 0.01, 3, and 1.5, respectively (Zhou et al., 2019). False-positive rate (FPR) analysis was eliminated using the Benjamini–Hochberg method with a q-value of 0.05 or lower (Zou et al., 2016). The enriched findings were displayed as bar and bubble plots on the bioinformatics website using the R package (http://www.bioinformatics.com.cn/). Comprehensive information on the most significantly enriched pathway was then extracted and colored (Kanehisa and Sato, 2020). Finally, a herb–compound–target–pathway–disease network was created using Cytoscape software (v.3.9.1, https://cytoscape.org/) to present the complicated network of the AB–DA herb pair in the treatment of osteoporosis.
2.5 Molecular docking to validate binding affinity
The SwissDock platform (Grosdidier et al., 2011) (http://www.swissdock.ch/) is an online molecular docking (MD) tool to determine the binding affinity from each binding site between small molecule ligands and receptor proteins. The X-ray diffraction of the protein crystal structure of key targets were downloaded from the Protein Data Bank (PDB) database (www.rcsb.org) (Nakamura et al., 2022). The binding sites were ranked based on their binding affinity scores, with the site having the smallest score considered the best binding site. Discovery Studio 2019 software (https://www.3ds.com) was used to visualize the binding details (Sultana et al., 2022).
2.6 Molecular dynamics simulation to validate binding stability
To investigate the stability of the complexes between small-molecule ligands and proteins, molecular dynamics simulations (MDS) were performed using the Standard Dynamics Cascade subunit of the Discovery Studio 2019 software. The ligand–protein complex with the lowest binding affinity score according to molecular docking analysis was selected (Hu et al., 2022). In this simulation system, water molecules are used to fill the solvent chamber, and Cl and Na+ ions are used to maintain an electrically neutral state. The simulation time was set as 300 ps, and the heating, balancing, and manufacturing phases were carried out after the system was balanced by an NPT ensemble, which fixed the pressure, temperature, and particle number. The analysis of the locus was performed using root mean square fluctuation (RMSF), root mean square deviation (RMSD), and hydrogen bond properties to produce the results.
2.7 Reagent and cell culture
Sitogluside, identified as one of the most promising bioactive compounds in the AB–DA herb pair, was further investigated for its anti-osteoporotic effects and associated mechanisms. Sitogluside was procured from the Dalian Meilunbio company and solubilized in dimethyl sulfoxide (DMSO). Human fetal osteoblast (hFOB) cells were obtained from the American Type Culture Collection (ATCC) and cultured in six-well plates using Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS, Gibco, United States), 1% penicillin–streptomycin (P/S), and 0.3 mg/mL Geneticin (G418). The JNK-specific inhibitor (SP600125) and p38-specific inhibitor (SB203580) were acquired from Absin (Shanghai). The cells were maintained in a humidified sterile atmosphere at 34 °C. They were subsequently treated with sitogluside when they reached a confluence of 60%–80% per well.
2.8 Alkaline phosphatase activity, Alizarin red staining mineralization, and osteoclast differentiation assays
For alkaline phosphatase (ALP) activity analysis, the BCIP/NBT Alkaline Phosphatase color development kit (Beyotime Institute of Biotechnology, Shanghai, China) was used according to the manufacturer’s procedure. Briefly, osteoblast precursor cells were seeded at 3 × 104 cells/well in 24-well plates and grown for 14 days in osteogenic media (DMEM +10% FBS +1% P/S +100 M ascorbic acid +2 mM 2-glycerophosphate +10 nM dexamethasone). The stained culture plates were photographed using a microscope (Leica image analysis system, Q500MC) and quantified using ImageJ software (National Institutes of Health, Bethesda, MD, United States). In addition, a 2% ARS reagent (Beyotime Institute of Biotechnology) was used to detect matrix mineralization, with the same protocol as ALP assay except the dye. In addition, we also investigated the effect of sitogluside on the osteoclast to fully illustrate the anti-OP effects. The RAW264.7 (osteoclast precursor) cells were treated with a different concentration of sitogluside and supplemented with 50 ng/mL receptor activator of NF-κB ligand (RANKL) as an osteoclast formation stimulator. Cells were then fixed in paraformaldehyde and stained with tartrate-resistant acid phosphatase (TRAP) after the intervention.
2.9 Ovariectomized mouse model
This study received ethical approval from the Animal Care Committee of the Second Xiangya Hospital of Central South University. A total of 30 10-week-old female C57BL/6 mice were procured from SLAC Laboratory Animal Co. Ltd. (SLACCAS, Shanghai, China). They were acclimatized in specific pathogen-free (SPF) cages for 1 week, during which measures were taken to minimize animal suffering through the use of anesthesia and sterile techniques during the surgical procedures. Bilateral ovariectomy (OVX) or sham surgery (retroperitoneal incision without ovariectomy) was performed based on group assignment, as described in the following paragraph.
Following the surgical procedures, the mice were randomly assigned to one of three groups: sham group (non-OVX mice, n = 10), vehicle group (OVX mice, n = 10), and sitogluside group (OVX mice intraperitoneally injected with sitogluside at a dose of 10 mg/kg/day, n = 10).
After 12 weeks, all mice were euthanized by cervical dislocation, and the right femurs were dissected and fixed in 4% paraformaldehyde (PFA) for 48 h. High-resolution micro-computed tomography (μCT40, Scanco, Zurich, Switzerland) was employed for bone analysis at the following parameters: scanning voltage = 80 kV, electric current = 80 μA, and resolution = 10 μm. The relevant trabecular bone volume fractions (BV/TV), trabecular number (Tb. N), trabecular thickness (Tb. Th), and trabecular separation (Tb. Sp) were subsequently calculated to assess the protective efficacy of sitogluside in OVX mice. In addition, to investigate the biosafety of sitogluside in the OVX model, the main organs were also obtained and detected using hematoxylin-eosin (H&E) staining. The heart, liver, spleen, lung, and kidney were hence fixed with formalin and embedded in paraffin cut to a 4 μm section. They were dewaxed in xylene, rehydrated with concentration gradient ethanol, and then stained with H&E for histological examinations and morphometric analysis. The serum of the mice was also collected to examine the biomarkers of alanine aminotransferase (ALT), creatine kinase (CK), and blood urea nitrogen (BUN).
2.10 Quantitative real-time PCR analysis
Human osteoblast (hFOB) cells were treated with 40 μM sitogluside or combined with the specific inhibitor of JNK and p38, depending on the groups, with ascorbic acid added in the osteogenic media. The cells were then harvested using the RNeasy Mini kit (QIAGEN, CA, United States) to extract total RNA following the manufacturer’s protocol. Subsequently, cDNA synthesis was performed using the reverse transcriptase kit (Takara Biotechnology, Japan). Real-time PCR analysis was performed using the SYBR Premix Ex Taq kit (Takara Biotechnology, Japan). The PCR parameters were set as follows: 40 cycles (denaturation at 95 °C for 10 s and amplification at 60 °C for 30 s). The resultant data were recorded as cycle threshold (Ct) values, and the 2−ΔΔCT method was employed for further analysis of RNA expression. In addition, to determine the modulated effect of sitogluside with JNK and osteogenic genes, the knockdown and activation of JNK expression were applied. ShRNA (shGnai3) was thus used to downregulate the expression level of JNK, and ASM was used as an activator to increase the p-JNK level (Meng et al., 2021). The alterations of osteogenic biomarkers were detected via q-PCR.
2.11 Western blot analysis
The hFOB cells were harvested by trypsin and then lysed in RIPA for 30 min on ice. Cell lysates were centrifuged at 12,000 g for 15 min at 4 C; the supernatant was collected, and protein content was quantified via the BSA protein assay kit following the manufacturer’s instruction. Proteins were separated by electrophoresis on 10%–12% SDS-PAGE at 100 V for 1.5 h and transferred onto a 0.45 μm polyvinylidene difluoride (PVDF) membrane at 250 mA for 1 h. The PVDF membrane was blocked with 5% non-fat milk in TBST buffer for 1 h at room temperature and incubated with primary antibody at 4 °C overnight. They were then incubated with secondary antibody for 1 h at room temperature and detected using the Chemiluminescence Kit.
2.12 Statistical analysis
Statistical analyses were conducted using GraphPad Prism 8.0.2 software (San Diego, United States). The data are presented as means ± standard deviation (SD). Data comparisons were performed using one-way analysis of variance (ANOVA), and statistical significance was determined by a p-value of less than 0.05.
3 Results
3.1 Relevant targets of AB–DA compound and osteoporosis
According to TCM databases, 185 components of AB and 82 compounds of DA were obtained in this work (Figures 3A, B). After filtration by screening criteria, 19 potential bioactive compounds of AB were obtained: arjunolic acid, baicalein, baicalin, berberine, chondrillasterol, coptisine, delta-7-stigmastenol, epiberberine, inophyllum E, kaempferol, oleanol, palmatine, quercetin, sitogluside, spinasterol, spinoside A, stigmasterol, wogonin, and β-ecdysterone. The 12 compounds from DA were 2,6-dihydroxycinnamic acid, caffeate, cauloside A, gentisin, isochlorogenic acid A, japonine, loganetin, loganic acid, loganin, sweroside, and sylvestroside III. Sitosterol is a common compound of both AB and DA (details shown in Supplementary Table S1; structures shown in Figure 2).
FIGURE 2. Chemical structure of 31 screened bioactive compounds. (A) Chemical structure of 19 bioactive compounds derived from AB. (B) Chemical structure of 12 bioactive compounds derived from DA; sitosterol is common to both.
After merging and duplicating results, 413 potential targets associated with 31 bioactive compounds were obtained, and 827 relevant targets of osteoporosis were acquired. Some 83 overlapped genes on Venn between AB–DA ingredients and osteoporosis were regarded as potential therapeutic targets (Figure 3C; details shown in Supplementary Table S2).
FIGURE 3. Data collection and hub gene screening for AB–DA against osteoporosis. (A) Venn diagram of the identified compounds of AB; data derived from TCMSP, TCM@Taiwan, BATMAN-TCM, and TCM-ID databases. (B) Venn diagram of the identified compounds of DA; data derived from the TCMSP, TCM-ID, and BATMAN-TCM databases and relevant literature. (C) Overlapping targets of AB–DA and osteoporosis, representing potential therapeutic targets. (D) PPI network analysis applied in the STRING database. (E) Plug-ins MCODE and CytoHubba of Cytoscape software to screen hub genes. (F) Degree value of top 10 hub genes.
3.2 PPI network of AB–DA against osteoporosis
A total of 83 targets of AB–DA against osteoporosis were imported to the STRING database; after deleting the disconnect targets, a PPI network with 78 nodes and 702 edges was constructed. Cytoscape software was utilized for further visualization, and plug-ins of MCODE and CytoHubba based on the topological parameters were applied to screen the hub genes, including AKT1, IGF1, CASP3, MMP9, EGFR, PPARG, ESR1, MAPK1, MAPK8, and MAPK14 (Figures 3D–F).
3.3 GO and KEGG pathway enrichment analyses of AB–DA against osteoporosis
GO and KEGG pathway enrichment analyses were performed in the Metascape database, and 83 targets with 875 GO items were enriched. It contained 772 biological processes (BP), 29 cellular components (CC), and 74 molecular functions (MF) items (p < 0.01, adjusted q < 0.05) (details shown in Figure 4A).
FIGURE 4. GO/KEGG enrichment analysis and D–C–T–P–D network of AB–DA in the treatment of osteoporosis. (A) Top 10 enriched GO items; left to right are biological process (BP), cellular compound (CC), and molecular function (MF), respectively. (B) Bubble diagram of top 20 enriched KEGG pathways. (C) GO chord chart presenting the corresponding relationship between core targets and pathways. (D) D–C–T network map; dark green diamond represents drugs (AB and DA), light green hexagon represents 31 bioactive compounds of AB–DA, and circles colored orange to red represent targets with low to high degrees. (E) MAPK signaling pathway mapped and colored by KEGG Mapper database; therapeutic targets of AB–DA are shown in red, targets of AB–DA but without therapeutic effect on osteoporosis are shown in blue, and other targets of osteoporosis are shown in yellow. (F) D–C–T–P–D network exhibits the regulatory mechanisms for AB–DA in the treatment of osteoporosis; dark green diamond represents drugs (AB and DA), light green hexagon represents compounds of AB–DA, circles colored orange to red represent targets with low to high degrees, blue V icon represents enriched core pathways, and yellow rectangle indicates the disease (osteoporosis). (G) the degree values of top 12 bioactive compounds.
A total of 135 KEGG pathways were significantly enriched, which mainly included the MAPK signaling pathway, osteoclast differentiation, PI3K-Akt signaling pathway, and endocrine resistance (Figures 4B, C). We further mapped and colored the regulation details of the MAPK cascade in the KEGG mapper database; in that map, red objects represent targets of AB–DA against osteoporosis, blue objects show the targets of AB–DA without the therapeutic effects of osteoporosis, and the other untargeted targets of osteoporosis are colored yellow (Figure 4E).
3.4 Drug–compound–target–pathway–disease network analysis
A drug–compound–target network was constructed using Cytoscape to illustrate core compounds and targets (Figure 4D); the core bioactive compounds included sitogluside, arjunolic acid, chondrillasterol, stigmasterol, spinasterol, spinoside A, cauloside A, sylvestroside III, β-ecdysterone, sitosterol, oleanol, and baicalin (ranked by degree value) (Figure 4G). A drug–compound–target–pathway–disease (D–C–T–P–D) network was then constructed to exhibit the complex molecular mechanisms of AB–DA anti-osteoporosis with multi-compound, multi-target, and multi-pathway characteristics. The dark green diamond represents the drugs (AB and DA), the light green hexagon represents the compounds of AB–DA, the circles colored orange to red represent targets with low to high degrees, the blue V icon represents the core enriched pathways, and the yellow rectangle indicates the disease (osteoporosis) (Figure 4F).
3.5 Molecular docking
The binding affinity between the core compounds and core targets are shown by heatmap (Figure 5A). According to relevant theories of molecular docking, the results of binding affinity < −5.0 kcal/mol suggest that there is a good spontaneous binding activity between molecule ligands and protein receptors, and results < −7.0 kcal/mol are stronger. Our research results showed that all the compounds had good binding activity with core targets, with binding affinities ranging from −5.64 kcal/mol to −9.43 kcal/mol. Sitogluside has the highest binding activity with IGF1. As shown in Figures 5B–D, the molecular interaction forces between IGF1 and sitogluside include π–donor hydrogen bond, π–alkyl bond, conventional hydrogen bond, carbon–hydrogen bond, and alkyl bond. The distances between the sitogluside atoms and amino acid residues of IGF1 range from 1.72 Å (number 1133, histidine residue) to 5.43 Å (number 1154, phenylalanine residue).
FIGURE 5. Molecular docking between bioactive compounds and core targets. (A) Binding affinity heatmap of compound ligand–protein receptor complexes, showing stronger binding affinity. (B) Binding details of the sitogluside–IGF1 complex (3D). (C) Binding details of the sitogluside–IGF1 complex (spatial structure). (D) Binding details of intermolecular force types of the sitogluside–IGF1 complex (2D).
3.6 Molecular dynamics simulation
The previous analysis of molecular docking showed the strong binding affinity between AB–DA compounds and core targets, and molecular dynamics simulation was utilized to identify the stability of the ligand–protein complex after docking by conformation alternation with potential energy under Newton’s law of motion. After the 300 ps simulation, the energy and temperature alternate tendency of the ensemble, hydrogen bond, RMSD, and RMSF changes of ligand–receptor interaction were calculated for stability analysis. RMSD was used to analyze the conformational alternation of receptors made by the ligand, and the results showed that curves and fluctuations only occurred at the beginning of the 80 ps simulation, and then tended to be stable (Figure 6A). The RMSF curve was utilized to monitor the conformational alternation of amino acid residues, and results show that the whole process is stable with only some small random fluctuations—which also reflects the whole ensemble’s stability to some extent (Figure 6B). The same tendency is also observed in the hydrogen bond heatmap (Figure 6C). The energies and temperature alternation of the whole ensemble were also stable and controlled. Therefore, the overall results exhibited good stability between the small-molecule ligand–protein receptor complex.
FIGURE 6. The molecular dynamic simulations to calculate the binding stability of sitogluside–IGF1 complex. (A), The total energy (a) and potential energy (b) curves of the whole ensemble in 300 ps simulation, showed the stability. (B), The temperature alteration of whole ensemble is controllable in 300 ps simulation. (C), The RMSD curve present the conformational alternation of receptor made by ligand. (D), The RMSF curve showed the conformational alternation of amino acid residues. (E), The hydrogen bond heat map of sitogluside–IGF1 complex.
3.7 Sitogluside promotes mineralization and ALP activity in osteoblast without any effect on osteoclast
To further investigate the impact of sitogluside on alkaline phosphatase (ALP) activity and mineralization in osteoblasts, we conducted ALP activity and ARS experiments. As illustrated in Supplementary Figure S3D, the sitogluside group exhibited increased ALP activity and mineralization compared to the control group. However, when shGnai3 intervention was introduced, the stimulatory effect of sitogluside on ALP activity and mineralization in osteoblasts was partially diminished. Subsequently, with the addition of the p-JNK activator ASM, the promotion of ALP activity and mineralization was restored. Additionally, our results indicated that sitogluside might not have a significant effect on osteoclast formation (Supplementary Figure S2).
3.8 Sitogluside administration protects against OVX-induced bone loss
The OVX model was employed to further investigate the anti-osteoporotic effects of sitogluside in vivo. As depicted in Figure 7A, the μCT scan results demonstrated a significant loss of bone in the vehicle group (OVX mice) compared to the sham group, confirming the successful establishment of the OVX model. Additionally, the sitogluside-treated group exhibited a mitigating effect on bone loss compared to the vehicle group. Specifically, the sitogluside-treated OVX mice demonstrated increased bone mineral density (BMD), trabecular thickness (Tb. Th), bone volume fraction (BV/TV), and trabecular number (Tb. N), while exhibiting a decreased bone-surface-to-bone-volume ratio (BS/BV) and trabecular separation (Tb. Sp). These findings collectively indicate the anti-osteoporotic efficacy of sitogluside in the OVX mouse model (Figure 7B). Furthermore, histological examination of major organs in the drug-treated group, including the heart, liver, spleen, lungs, and kidneys, was conducted using (H&E) staining. Additionally, mouse serum was analyzed to assess cardiac, hepatic, and renal functions. As depicted in Supplementary Figure S1, the experimental results indicate that sitogluside did not exert significant toxic effects on mouse organs.
FIGURE 7. In vivo OVX model and q-PCR experiments to validate the potential anti-OP effect of sitogluside. (A) μCT scanning of mice tibia showing the bone loss alleviated by sitogluside’s efficiency. (B) Statistical results of osteoporosis phenotype parameters of OVX mice, including BMD, Tb. Th, BV/TV, and Tb. N. (C) q-PCR results of osteogenic biomarker, which is shown to promote the osteogenic effect of sitogluside. (D) q-PCR results of OSX, which, treated with sitogluside or combined with the specific inhibitor of JNK (SP600125) and p38 (SB203580), indicates that sitogluside could target the JNK pathway to promote osteogenic genes in the treatment of OP.
3.9 Effects of sitogluside on osteoblast-related genes
To gain further insights into the underlying mechanism of sitogluside’s anti-osteoporotic effects, we examined the expression of osteogenic-related genes using q-PCR and Western blot. As shown in Figure 7, treatment with sitogluside for 48 h resulted in a significant increase in the expression of osteogenic markers, such as Osterix (OSX) and osteocalcin (OCN). These findings suggest that sitogluside has the potential to promote osteogenic activity. Additionally, we investigated the core target and pathway associated with sitogluside in the treatment of osteoporosis via blocking the JNK and p38 cascade (Figure 7D). Our results revealed that targeting JNK could reduce the osteogenic efficacy of sitogluside while blocking the p38 pathway without significant expression changes of OSX. To further investigate whether sitogluside exerts its osteogenic effects through the JNK pathway, we employed shRNA to silence the JNK pathway. As shown in Supplementary Figure S3C, compared to the sitogluside group, silencing the JNK pathway with shRNA (shGnai3) resulted in a partial reduction in the expression levels of osteogenic markers, including Runx2, OSX, OCN, and ALP. However, when treated with the p-JNK activator anisomycin (ASM), these osteogenic markers increase. Moreover, as shown in Supplementary Figures S3A, B, the phosphorylation level of p-JNK increases with sitogluside intervention. These further suggest that sitogluside may promote osteoblast differentiation through the JNK signaling pathway. We also investigated the effect of sitogluside on the Smad 1/5/8 protein in the BMP signaling pathway, which also showed increasing expression (Supplementary Figures S3A, B).
4 Discussion
Osteoporosis is regarded as a silent disease without clinical symptoms before complications are apparent (Anthamatten and Parish, 2019). Current conventional therapy mainly focuses on symptom prevention with long-term supplements of calcium and on intervention to regulate bone metabolism. The curative effects depend on the individual response and have mild or severe side effects (Ensrud and Crandall, 2017). AB and DA have been two of the most important herbs for bone diseases in TCM therapy for more than 2,000 years (He et al., 2017; Tao et al., 2020), and play an important role in many classic anti-osteoporosis drugs. Recently, researchers have isolated more than 100 ingredients and verified numerous bioactivities (He et al., 2017; Tao et al., 2020). Our research obtained 31 bioactive compounds of AB–DA based on the screening criteria and applied the main therapeutic compounds to anti-osteoporosis. Compounds of AB–DA such as baicalin, kaempferol, oleanol, quercetin, and sitosterol are widespread in herbal medicine, and many studies have shown their anti-osteoporosis activity. It has been reported that kaempferol could ameliorate the inhibitory effects of osteogenesis by activating JNK and p38 pathways in the glucocorticoid-induced and ovariectomy-induced osteoporosis model (Wong et al., 2019). Vakili et al. (2021) also constructed an OVX osteoporosis model treated with quercetin; results showed that it might modulate cell autophagy and apoptosis to alleviate osteoporosis; the potential mechanisms mainly involve Wnt, NF-κB, and MAPK cascades. Moreover, asperosaponin VI is the quality indicator of DA and could induce the differentiation of osteoblastic cells by increased expression of BMP2, promote osteogenesis and angiogenesis via regulating the OPG/RANKL signaling pathway, and inhibit the differentiation of osteoclast (Chen et al., 2022). β-Ecdysterone is an iconic ingredient of AB and could upregulate the activity of alkaline phosphatase in mesenchymal stem cells by modulating the expression of estrogen receptors (Gao et al., 2008) and could also inhibit apoptosis and autophagy induced by the dexamethasone of osteoblast cells in vivo and in vitro (Tang et al., 2018). This body of evidence suggests that the numerous bioactive compounds retrieved from AB–DA have certain anti-OP effects and partly confirm our findings; however, more accurate and in-depth study is still necessary to explore the AB–DA herb pair.
After screening the core therapeutic ingredients, a PPI network was constructed to determine core targets; our results showed that AKT1, CASP3, EGFR, ESR1, IGF1, MAPK1, and MAPK14 are important for AB–DA treatment of osteoporosis. It is widely reported that RAC-alpha serine/threonine-protein kinase (AKT1) regulates a series of biological processes, including cell proliferation, growth, metabolism, and angiogenesis (Heron-Milhavet et al., 2011). A targeted knockdown AKT1 mouse model showed that AKT1 deficiency would induce osteoclast-osteogenesis disorder and diminish the RANKL (NF-κB ligand) and MCSF (macrophage colony-stimulating factor) receptors on multinucleated osteoclasts. This is evidence for AKT1 as an intermediator to regulate osteoblast and osteoclast differentiation (Mukherjee and Rotwein, 2012). Wang et al. (2022) also showed that modulating the expression of AKT1 could relieve osteoporosis. Epidermal growth factor receptor (EGFR) could activate the downstream of MEK-ERK, PI3K-AKT, and NF-κB signaling, transferring extracellular cues into cellular response (Runkle et al., 2016). It was reported that promoting the phosphorylation of EGFR and ERK1/2 could alleviate the apoptosis induced by H2O2 of MC3T3-E1 (osteoblast cell) (Yang et al., 2019), exhibiting the potential targeted therapeutic value of EGFR on osteoporosis. Insulin-like growth factor I (IGF1) is the most abundant growth promotor in the bone matrix and also regulates glycogen synthesis in osteoblasts; it plays an important role in bone homeostasis maintenance and osteoblast differentiation by mediating the mTOR (mammalian target of rapamycin) signaling pathway (Xian et al., 2012). The results of molecular docking and molecular dynamics simulation also showed a stronger binding affinity between AB–DA compounds and IGF1. All the aforementioned evidence shows that the hub genes regulated by AB–DA compounds are important and meaningful in the treatment of osteoporosis.
The GO/KEGG pathway enrichment analysis of 83 potential therapeutic targets showed that AB–DA anti-osteoporosis mainly involved the MAPK signaling pathway, osteoclast differentiation, and the PI3K-Akt signaling pathway. As shown in Figure 4E, the colored target details of the MAPK pathway—ERK (MAPK1), JNK(MAPK8), and p38 (MAPK14)—are three classical cascades of the MAPK pathway, all of which were involved in the potential mechanism of AB–DA treating osteoporosis. ERK cascade mediated cell growth and differentiation via cytoskeletal rearrangement, and upregulating the phosphorylation of ERK (p-ERK) might promote osteoblast differentiation (Jing et al., 2019). The JNK and p38 cascades would be activated by extracellular stimulation, such as pro-inflammatory and physical stress. Lee et al. (2019) showed that downregulated osteoclast-related gene expression was associated with JNK cascade inhibition and that suppressing the p38 cascade would also relieve osteoporosis (Wang et al., 2018). Liu et al. (2022) indicated that vitexin could act against osteoporosis by promoting osteogenesis and angiogenesis in an ovariectomized rat model; the potential underlying mechanism might upregulate the PI3K-AKT cascade. Regulated PI3K-AKT signaling could also mediate the biological function of osteoclast (Jiang et al., 2022). Thus, our research mined the herb databases and screened therapeutic compounds of AB–DA with appropriate pharmacokinetic properties. In summary, the 31 therapeutic compounds have different targets and regulate different signaling pathways with a synergistic effect against osteoporosis, showing the complex molecular mechanisms with “multi-compound,” “multi-target,” and “multi-pathway.”
Furthermore, in vitro and in vivo experimental validation were both employed to explore the anti-osteoporotic effects of sitogluside (also known as daucosterol) and its underlying molecular mechanism. Some research has indicated the regulation of both the osteoblast and osteoclast formation of AB–DA (He et al., 2017; Tao et al., 2020). In this present research, our findings show the potential effects of sitogluside with osteoblast differentiation and mineralization, but the effects on osteoclast were not significant; this was also verified by the collaborative therapeutic effects of AB–DA compounds. Moreover, q-PCR and Western blot analyses showed that sitogluside might upregulate the JNK cascade to promote osteogenics, such as Runx2, Osx, and OCN, and knockdown or block with the inhibitor could partly reduce its efficacy. Previous research has also reported that daucosterol could increase the p-JNK expression to exert an anti-prostate cancer effect (Gao et al., 2019). Huang et al. reported that the JNK kinase pathway with downstream OSX belonged to the non-canonical Smad-independent BMP signaling pathway to promote osteogenics. Hence, our study illustrated that a potential mechanism for sitogluside in the treatment of osteoporosis was to promote the JNK pathway and non-canonical BMP signaling to regulate downstream osteogenic genes. In addition, the OVX mouse model also showed anti-osteoporosis efficacy in vivo without observable toxicity. Our results therefore highlight the potential therapeutic value of sitogluside in the treatment of osteoporosis. However, further systematic and in-depth research is required to explore other bioactive compounds of the AB–DA combination.
There were some limitations to our research. First, it is based on the bioinformatic analysis of network pharmacology, molecular docking, and molecular dynamics simulation. The results illustrated the core therapeutic ingredients, core targets, and core signaling pathways of AB–DA in the treatment of osteoporosis, although we have preliminarily verified the anti-osteoporotic effect of sitogluside. Our findings thus need more wet experiments in vitro and in vivo for corroboration. The second limitation is that all the bioactive compounds were filtered by the ADME properties; the content of specific ingredients were not considered, and the toxicological information was ignored. Therefore, the next task for us is to comprehensively verify our network pharmacological findings with experiments and evaluate the biosafety of these ingredients for optimally utilizing AB–DA treatment of osteoporosis during clinical practice (Zou et al., 2021).
5 Conclusion
This study represents the first comprehensive investigation into the bioactive compounds of the AB–DA herb pair. Our findings illustrate that the primary therapeutic compounds responsible for treating osteoporosis are iridoids, steroids, and flavonoids. Additionally, we propose that the underlying mechanisms of action may involve targeting key core targets, including MAPK1, MAPK8, and MAPK14, to modulate the MAPK cascade (ERK/JNK/p38). Furthermore, in vitro and in vivo experiments have provided preliminary validation of the anti-osteoporotic effect of the most potent bioactive compound, sitogluside. Nevertheless, further in-depth experiments are needed to fully harness the therapeutic potential of AB–DA for treating osteoporosis in clinical practice.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The animal study was approved by the Animal Care Committee of the Second Xiangya Hospital of Central South University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
ZO designed this study; CX wrote the manuscript; TL collected the data and performed statistical analysis; WL and TT revised the manuscript; XG created all tables, figures, and prepared the Supplementary Material. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the National Natural Science Foundation of China (No. 82371600), the Natural Science Foundation of Hunan Province (No. 2021JJ40810 and 2023JJ40833), the Natural Science Foundation of Hunan Province of Outstanding Young Scholars (Grant No. 2021JJ20086), the Scientific Research Projects of the Health Commission of Hunan Province (project number B202304077119), the Natural Science Foundation of Hunan Province (grant number 2021JJ30930), and the Fundamental Research Funds for the Central Universities of Central South University (grant number 2023ZZTS0558).
Acknowledgments
The authors sincerely appreciate the TCMSP, TCM@Taiwan, BATMAN-TCM, TCM-ID, PharmMapper, Swiss Institute of Bioinformatics, UniProt, GeneCards, DisGeNET, OMIM, STRING, and Metascape databases and Cytoscape and Discovery Studio 2019 software for data collection, statistical analysis, and figure plotting.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1242194/full#supplementary-material
SUPPLEMENTARY TABLE S1 | Detailed information of AB–DA herb pair bioactive compounds.
SUPPLEMENTARY TABLE S2 | Protein name and Entrez ID of all targets.
References
Anthamatten, A., and Parish, A. (2019). Clinical update on osteoporosis. J. Midwifery Womens Health 64, 265–275. doi:10.1111/jmwh.12954
Barshir, R., Fishilevich, S., Iny-Stein, T., Zelig, O., Mazor, Y., Guan-Golan, Y., et al. (2021). GeneCaRNA: A comprehensive gene-centric database of human non-coding RNAs in the GeneCards suite. J. Mol. Biol. 433, 166913. doi:10.1016/j.jmb.2021.166913
Chen, C. Y. (2011). TCM Database@Taiwan: the world's largest traditional Chinese medicine database for drug screening in silico. PLoS One 6, e15939. doi:10.1371/journal.pone.0015939
Chen, F., Liang, Q., Mao, L., Yin, Y., Zhang, L., Li, C., et al. (2022). Synergy effects of Asperosaponin VI and bioactive factor BMP-2 on osteogenesis and anti-osteoclastogenesis. Bioact. Mater 10, 335–344. doi:10.1016/j.bioactmat.2021.09.001
Chu, H., Moon, S., Park, J., Bak, S., Ko, Y., and Youn, B. Y. (2022). The use of artificial intelligence in complementary and alternative medicine: A systematic scoping review. Front. Pharmacol. 13, 826044. doi:10.3389/fphar.2022.826044
Compston, J. E., McClung, M. R., and Leslie, W. D. (2019). Osteoporosis. Lancet 393, 364–376. doi:10.1016/S0140-6736(18)32112-3
Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 7, 42717. doi:10.1038/srep42717
Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic acids Res. 47, W357–w364. doi:10.1093/nar/gkz382
Ensrud, K. E., and Crandall, C. J. (2017). Osteoporosis. Ann. Intern Med. 167, Itc17–itc32. doi:10.7326/aitc201708010
Feng, S., Wang, T., Fan, L., An, X., Ding, X., Wang, M., et al. (2022). Exploring the potential therapeutic effect of Eucommia ulmoides-Dipsaci Radix herbal pair on osteoporosis based on network pharmacology and molecular docking technology. RSC Adv. 12, 2181–2195. doi:10.1039/d1ra05799e
Gao, L., Cai, G., and Shi, X. (2008). Beta-ecdysterone induces osteogenic differentiation in mouse mesenchymal stem cells and relieves osteoporosis. Biol. Pharm. Bull. 31, 2245–2249. doi:10.1248/bpb.31.2245
Gao, P., Huang, X., Liao, T., Li, G., Yu, X., You, Y., et al. (2019). Daucosterol induces autophagic-dependent apoptosis in prostate cancer via JNK activation. Biosci. Trends 13, 160–167. doi:10.5582/bst.2018.01293
Gao, Y., Ji, W., Lu, M., Wang, Z., Jia, X., Wang, D., et al. (2022). Systemic pharmacological verification of Guizhi Fuling decoction in treating endometriosis-associated pain. J. Ethnopharmacol. 297, 115540. doi:10.1016/j.jep.2022.115540
Grewe, J. M., Knapstein, P. R., Donat, A., Jiang, S., Smit, D. J., Xie, W., et al. (2022). The role of sphingosine-1-phosphate in bone remodeling and osteoporosis. Bone Res. 10, 34. doi:10.1038/s41413-022-00205-0
Grosdidier, A., Zoete, V., and Michielin, O. (2011). SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic acids Res. 39, W270–W277. doi:10.1093/nar/gkr366
He, C. C., Hui, R. R., Tezuka, Y., Kadota, S., and Li, J. X. (2010). Osteoprotective effect of extract from Achyranthes bidentata in ovariectomized rats. J. Ethnopharmacol. 127, 229–234. doi:10.1016/j.jep.2009.11.016
He, X., Wang, X., Fang, J., Chang, Y., Ning, N., Guo, H., et al. (2017). The genus Achyranthes: A review on traditional uses, phytochemistry, and pharmacological activities. J. Ethnopharmacol. 203, 260–278. doi:10.1016/j.jep.2017.03.035
Heron-Milhavet, L., Khouya, N., Fernandez, A., and Lamb, N. J. (2011). Akt1 and Akt2: Differentiating the aktion. Histol. Histopathol. 26, 651–662. doi:10.14670/HH-26.651
Holzhüter, K., and Geertsma, E. R. (2022). Uniport, not proton-symport, in a non-mammalian SLC23 transporter. J. Mol. Biol. 434, 167393. doi:10.1016/j.jmb.2021.167393
Hu, Y., Wu, Y., Jiang, C., Wang, Z., Shen, C., Zhu, Z., et al. (2022). Investigative on the molecular mechanism of licorice flavonoids anti-melanoma by network pharmacology, 3D/2D-QSAR, molecular docking, and molecular dynamics simulation. Front. Chem. 10, 843970. doi:10.3389/fchem.2022.843970
Huang, A. Y., Xiong, Z., Liu, K., Chang, Y., Shu, L., Gao, G., et al. (2022). Identification of kaempferol as an OSX upregulator by network pharmacology-based analysis of qianggu Capsule for osteoporosis. Front. Pharmacol. 13, 1011561. doi:10.3389/fphar.2022.1011561
Huang, Y., Li, N., Wan, J. B., Zhang, D., and Yan, C. (2015). Structural characterization and antioxidant activity of a novel heteropolysaccharide from the submerged fermentation mycelia of Ganoderma capense. Carbohydr. Polym. 134, 752–760. doi:10.1016/j.carbpol.2015.08.067
Jiang, M., Liu, R., Liu, L., Kot, A., Liu, X., Xiao, W., et al. (2020). Identification of osteogenic progenitor cell-targeted peptides that augment bone formation. Nat. Commun. 11, 4278. doi:10.1038/s41467-020-17417-9
Jiang, T., Gu, H., and Wei, J. (2022). Echinacoside inhibits osteoclast function by down-regulating PI3K/akt/C-fos to alleviate osteolysis caused by periprosthetic joint infection. Front. Pharmacol. 13, 930053. doi:10.3389/fphar.2022.930053
Jiang, Y., Zhang, Y., Chen, W., Liu, C., Li, X., Sun, D., et al. (2014). Achyranthes bidentata extract exerts osteoprotective effects on steroid-induced osteonecrosis of the femoral head in rats by regulating RANKL/RANK/OPG signaling. J. Transl. Med. 12, 334. doi:10.1186/s12967-014-0334-7
Jin, L. Q., Zheng, Z. J., Peng, Y., Li, W. X., Chen, X. M., and Lu, J. X. (2007). Opposite effects on tumor growth depending on dose of Achyranthes bidentata polysaccharides in C57BL/6 mice. Int. Immunopharmacol. 7, 568–577. doi:10.1016/j.intimp.2006.12.009
Jing, Z., Wang, C., Yang, Q., Wei, X., Jin, Y., Meng, Q., et al. (2019). Luteolin attenuates glucocorticoid-induced osteoporosis by regulating ERK/Lrp-5/GSK-3β signaling pathway in vivo and in vitro. J. Cell Physiol. 234, 4472–4490. doi:10.1002/jcp.27252
Kanehisa, M., and Sato, Y. (2020). KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. a Publ. Protein Soc. 29, 28–35. doi:10.1002/pro.3711
Kang, H., Tang, K., Liu, Q., Sun, Y., Huang, Q., Zhu, R., et al. (2013). HIM-herbal ingredients in-vivo metabolism database. J. cheminformatics 5, 28. doi:10.1186/1758-2946-5-28
Liu, Z., Guo, F., Wang, Y., Li, C., Zhang, X., Li, H., et al. (2016). BATMAN-TCM: A bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine. Sci. Rep. 6, 21146. doi:10.1038/srep21146
Lee, K. M., Lee, C. Y., Zhang, G., Lyu, A., and Yue, K. K. M. (2019). Methylglyoxal activates osteoclasts through JNK pathway leading to osteoporosis. Chem. Biol. Interact. 308, 147–154. doi:10.1016/j.cbi.2019.05.026
Li, T., Jiang, G., Hu, X., Yang, D., Tan, T., Gao, Z., et al. (2021). Punicalin attenuates breast cancer-associated osteolysis by inhibiting the NF-κB signaling pathway of osteoclasts. Front. Pharmacol. 12, 789552. doi:10.3389/fphar.2021.789552
Li, Z., Ying, B., Liu, X., Zhang, X., and Yu, H. (2012). An examination of the OMIM database for associating mutation to a consensus reference sequence. Protein & Cell 3, 198–203. doi:10.1007/s13238-012-2037-2
Liang, B., Burley, G., Lin, S., and Shi, Y. C. (2022). Osteoporosis pathogenesis and treatment: Existing and emerging avenues. Cell Mol. Biol. Lett. 27, 72. doi:10.1186/s11658-022-00371-3
Liu, K., Liu, Y., Xu, Y., Nandakumar, K. S., Tan, H., He, C., et al. (2019). Asperosaponin VI protects against bone destructions in collagen induced arthritis by inhibiting osteoclastogenesis. Phytomedicine Int. J. phytotherapy Phytopharm. 63, 153006. doi:10.1016/j.phymed.2019.153006
Liu, X., Ouyang, S., Yu, B., Liu, Y., Huang, K., Gong, J., et al. (2010). PharmMapper server: A web server for potential drug target identification using pharmacophore mapping approach. Nucleic acids Res. 38, W609–W614. doi:10.1093/nar/gkq300
Liu, Y., Xue, Q., Li, A., Li, K., and Qin, X. (2020). Mechanisms exploration of herbal pair of HuangQi-DanShen on cerebral ischemia based on metabonomics and network pharmacology. J. Ethnopharmacol. 253, 112688. doi:10.1016/j.jep.2020.112688
Liu, Y., Zhu, S., Liu, J., Chen, Y., Zhong, S., Lian, D., et al. (2022). Vitexin regulates angiogenesis and osteogenesis in ovariectomy-induced osteoporosis of rats via the VDR/PI3K/AKT/eNOS signaling pathway. J. Agric. Food Chem. 71, 546–556. doi:10.1021/acs.jafc.2c07005
McDonald, M. M., Khoo, W. H., Ng, P. Y., Xiao, Y., Zamerli, J., Thatcher, P., et al. (2021). Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption. Cell 184, 1940–1347.e13. doi:10.1016/j.cell.2021.03.010
Meng, L., Yuan, L., Ni, J., Fang, M., Guo, S., Cai, H., et al. (2021). Mir24-2-5p suppresses the osteogenic differentiation with Gnai3 inhibition presenting a direct target via inactivating JNK-p38 MAPK signaling axis. Int. J. Biol. Sci. 17, 4238–4253. doi:10.7150/ijbs.60536
Mukherjee, A., and Rotwein, P. (2012). Selective signaling by Akt1 controls osteoblast differentiation and osteoblast-mediated osteoclast development. Mol. Cell Biol. 32, 490–500. doi:10.1128/MCB.06361-11
Mukwaya, E., Xu, F., Wong, M. S., and Zhang, Y. (2014). Chinese herbal medicine for bone health. Pharm. Biol. 52, 1223–1228. doi:10.3109/13880209.2014.884606
Nakamura, R., Ogawa, S., Takahashi, Y., and Fujishiro, T. (2022). Cycloserine enantiomers inhibit PLP-dependent cysteine desulfurase SufS via distinct mechanisms. FEBS J. 289, 5947–5970. doi:10.1111/febs.16455
NIH Consensus Development Panel on Osteoporosis Prevention, Diagnosis, and Therapy (2001). Osteoporosis prevention, diagnosis, and therapy. Jama 285, 785–795. doi:10.1001/jama.285.6.785
Pagadala, N. S., Syed, K., and Tuszynski, J. (2017). Software for molecular docking: A review. Biophys. Rev. 9, 91–102. doi:10.1007/s12551-016-0247-1
Peng, Z., Xu, R., and You, Q. (2022). Role of traditional Chinese medicine in bone regeneration and osteoporosis. Front. Bioeng. Biotechnol. 10, 911326. doi:10.3389/fbioe.2022.911326
Piñero, J., Ramírez-Anguita, J. M., Saüch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic acids Res. 48, D845–d855. doi:10.1093/nar/gkz1021
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: A database of systems pharmacology for drug discovery from herbal medicines. J. cheminformatics 6, 13. doi:10.1186/1758-2946-6-13
Runkle, K. B., Kharbanda, A., Stypulkowski, E., Cao, X. J., Wang, W., Garcia, B. A., et al. (2016). Inhibition of DHHC20-mediated EGFR palmitoylation creates a dependence on EGFR signaling. Mol. Cell 62, 385–396. doi:10.1016/j.molcel.2016.04.003
Shuai, Y., Jiang, Z., Yuan, Q., Tu, S., and Zeng, F. (2020). Deciphering the underlying mechanism of eucommiae cortex against osteoporotic fracture by network pharmacology. eCAM 2020, 7049812. doi:10.1155/2020/7049812
Suh, K. S., Lee, Y. S., and Choi, E. M. (2014). The protective effects of Achyranthes bidentata root extract on the antimycin A induced damage of osteoblastic MC3T3-E1 cells. Cytotechnology 66, 925–935. doi:10.1007/s10616-013-9645-4
Sultana, N., Chung, H. J., Emon, N. U., Alam, S., Taki, M. T. I., Rudra, S., et al. (2022). Biological functions of Dillenia pentagyna roxb. Against pain, inflammation, fever, diarrhea, and thrombosis: Evidenced from in vitro, in vivo, and molecular docking study. Front. Nutr. 9, 911274. doi:10.3389/fnut.2022.911274
Suvarna, V., Sarkar, M., Chaubey, P., Khan, T., Sherje, A., Patel, K., et al. (2018). Bone health and natural products- an insight. Front. Pharmacol. 9, 981. doi:10.3389/fphar.2018.00981
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic acids Res. 49, D605–d612. doi:10.1093/nar/gkaa1074
Tang, Y. H., Yue, Z. S., Xin, D. W., Zeng, L. R., Xiong, Z. F., Hu, Z. Q., et al. (2018). β-Ecdysterone promotes autophagy and inhibits apoptosis in osteoporotic rats. Mol. Med. Rep. 17, 1591–1598. doi:10.3892/mmr.2017.8053
Tao, Y., Chen, L., and Yan, J. (2020). Traditional uses, processing methods, phytochemistry, pharmacology and quality control of dipsacus asper wall. Ex C.B. Clarke: A review. J. Ethnopharmacol. 258, 112912. doi:10.1016/j.jep.2020.112912
Vakili, S., Zal, F., Mostafavi-Pour, Z., Savardashtaki, A., and Koohpeyma, F. (2021). Quercetin and vitamin E alleviate ovariectomy-induced osteoporosis by modulating autophagy and apoptosis in rat bone cells. J. Cell Physiol. 236, 3495–3509. doi:10.1002/jcp.30087
Wang, W. D., Kang, W. B., Zhou, X. Q., Yao, G. F., and Wang, X. J. (2018). Mitochondrial protein OPA mediates osteoporosis induced by radiation through the P38 signaling pathway. Eur. Rev. Med. Pharmacol. Sci. 22, 8091–8097. doi:10.26355/eurrev_201812_16499
Wang, Y., Liu, L., Qu, Z., Wang, D., Huang, W., Kong, L., et al. (2022). Tanshinone ameliorates glucocorticoid-induced bone loss via activation of AKT1 signaling pathway. Front. Cell Dev. Biol. 10, 878433. doi:10.3389/fcell.2022.878433
Wong, S. K., Chin, K. Y., and Ima-Nirwana, S. (2019). The osteoprotective effects of kaempferol: The evidence from in vivo and in vitro studies. Drug Des. Devel Ther. 13, 3497–3514. doi:10.2147/DDDT.S227738
Wu, H., Lv, Y., Wei, F., Li, C., Ge, W., and Du, W. (2022). Comparative analysis of anti-osteoporosis efficacy in Radix Dipsaci before and after processing with salt based on spectrum-effect relationship. J. Pharm. Biomed. analysis 221, 115078. doi:10.1016/j.jpba.2022.115078
Xian, L., Wu, X., Pang, L., Lou, M., Rosen, C. J., Qiu, T., et al. (2012). Matrix IGF-1 maintains bone mass by activation of mTOR in mesenchymal stem cells. Nat. Med. 18, 1095–1101. doi:10.1038/nm.2793
Yan, C., Zhang, S., Wang, C., and Zhang, Q. (2019). A fructooligosaccharide from Achyranthes bidentata inhibits osteoporosis by stimulating bone formation. Carbohydr. Polym. 210, 110–118. doi:10.1016/j.carbpol.2019.01.026
Yang, F., Lin, Z. W., Huang, T. Y., Chen, T. T., Cui, J., Li, M. Y., et al. (2019). Ligustilide, a major bioactive component of Angelica sinensis, promotes bone formation via the GPR30/EGFR pathway. Sci. Rep. 9, 6991. doi:10.1038/s41598-019-43518-7
Ye, X. W., Wang, H. L., Cheng, S. Q., Xia, L. J., Xu, X. F., and Li, X. R. (2022). Network pharmacology-based strategy to investigate the pharmacologic mechanisms of coptidis rhizoma for the treatment of alzheimer's disease. Front. aging Neurosci. 14, 890046. doi:10.3389/fnagi.2022.890046
Yu, J. H., Yu, Z. P., Wang, Y. Y., Bao, J., Zhu, K. K., Yuan, T., et al. (2019). Triterpenoids and triterpenoid saponins from Dipsacus asper and their cytotoxic and antibacterial activities. Phytochemistry 162, 241–249. doi:10.1016/j.phytochem.2019.03.028
Zhang, M., Wang, Y., Zhang, Q., Wang, C., Zhang, D., Wan, J. B., et al. (2018a). UPLC/Q-TOF-MS-based metabolomics study of the anti-osteoporosis effects of Achyranthes bidentata polysaccharides in ovariectomized rats. Int. J. Biol. Macromol. 112, 433–441. doi:10.1016/j.ijbiomac.2018.01.204
Zhang, R., Hu, S. J., Li, C., Zhang, F., Gan, H. Q., and Mei, Q. B. (2012). Achyranthes bidentata root extract prevent OVX-induced osteoporosis in rats. J. Ethnopharmacol. 139, 12–18. doi:10.1016/j.jep.2011.05.034
Zhang, S., Zhang, Q., Zhang, D., Wang, C., and Yan, C. (2018b). Anti-osteoporosis activity of a novel Achyranthes bidentata polysaccharide via stimulating bone formation. Carbohydr. Polym. 184, 288–298. doi:10.1016/j.carbpol.2017.12.070
Zhang, W., Xue, K., Gao, Y., Huai, Y., Wang, W., Miao, Z., et al. (2019). Systems pharmacology dissection of action mechanisms of Dipsaci Radix for osteoporosis. Life Sci. 235, 116820. doi:10.1016/j.lfs.2019.116820
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10, 1523. doi:10.1038/s41467-019-09234-6
Zou, M. L., Chen, Z. H., Teng, Y. Y., Liu, S. Y., Jia, Y., Zhang, K. W., et al. (2021). The Smad dependent TGF-β and BMP signaling pathway in bone remodeling and therapies. Front. Mol. Biosci. 8, 593310. doi:10.3389/fmolb.2021.593310
Keywords: osteoporosis, Achyranthes bidentata, Dipsacus asper, network pharmacology, MAPK
Citation: Li T, Li W, Guo X, Tan T, Xiang C and Ouyang Z (2023) Unraveling the potential mechanisms of the anti-osteoporotic effects of the Achyranthes bidentata–Dipsacus asper herb pair: a network pharmacology and experimental study. Front. Pharmacol. 14:1242194. doi: 10.3389/fphar.2023.1242194
Received: 18 June 2023; Accepted: 07 September 2023;
Published: 02 October 2023.
Edited by:
Xinhua Qu, Shanghai Jiao Tong University, ChinaReviewed by:
Jun Wan, The University of Queensland, AustraliaVenkatesh Katari, University of Toledo, United States
Copyright © 2023 Li, Li, Guo, Tan, Xiang and Ouyang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhengxiao Ouyang, ouyangzhengxiao@csu.edu.cn; Wenzhao Li, liwenzhao@csu.edu.cn
†These authors have contributed equally to this work