Abstract
Background. Eucommiae Cortex is a Chinese herbal medicine with bone protective effects and treats osteoporosis. This study aimed to explore the pharmacological mechanisms of this complex mixture. Methods. The active compounds and disease targets involved in the study were obtained from publicly available websites and databases. Core target genes were identified by protein-protein interaction (PPI) network and topology analysis and mapped to critical components by a “component-target” regulatory network. Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway analysis and Gene Ontology (GO)analysis were performed to analyze the biological processes of target genes. Moreover, we carried out molecular docking, cell experiments, and quantitative real-time PCR to propel the research forward. Results. Eucommiae Cortex contained 28 active ingredients and 85 potential targets for antiosteoporosis. PPI networks and topology analysis screened 17 core target genes. The therapeutic mechanism involves a series of biological reactions, including host-virus response, chemical stress, oxidative stress, and cell cycle control. The KEGG enrichment illustrated that the MAPK signaling pathway might play a significant role. The final experiment detected ten genes (EGF, AKT1, JUN, MAPK8, MAPK1, CASP3, FOS, VEGFA, EGFR, and MYC) and three compounds (quercetin, kaempferol, and beta- carotene) in the MAPK pathway. Firstly, the CCK-8 and ALP activity test results showed that three compounds could enhance the proliferation and differentiation of osteoblast-like MC3T3-E1 cells. Secondly, molecular docking confirmed the favorable binding potential. Subsequently, we observed that adding mol/L quercetin, mol/L kaempferol, and mol/L beta-carotene activated the ERK/JNK cascades and the heterodimer complex AP-1(Fos/Jun) in the MAPK pathway. Conclusion. MAPK pathway might provide an essential mechanism for the antiosteoporosis effect of Eucommiae Cortex. Eucommiae Cortex and its active ingredients have the potential to treat osteoporosis.
1. Introduction
Osteoporosis is a common chronic metabolic bone disease with considerable social and health implications. The prevalence of osteoporosis increased with age, bringing safety hazards to the elderly, including signs of augmented fracture [1]. According to the latest statistics of the International Osteoporosis Foundation, one-third of women over 50 years old and one-fifth of men will encounter osteoporosis fractures in their life [2]. Osteoporotic fractures account for 0.83% of the global burden of noncommunicable diseases, and this share is likely to rise [3]. Unfortunately, osteoporosis requires lifelong medication to prevent fractures. A recent meta-analysis found that bisphosphonate benefited 1% of potential female patients after a minimum of 12.4 months of continuous use [4]. Although the advent of newer drugs such as teriparatide and denosumab ushered at the dawn of antiosteoporosis treatment, evidence of long-term benefit is absent [5]. The long-term efficacy and safety of first-line therapy for osteoporosis remain an issue of concern. There is a strong demand to develop drugs with more moderate side effects and boost bone metabolism over the long term.
Eucommiae Cortex, a traditional Chinese medicine (TCM), is widely used to treat osteoporosis, lumbago, and hypertension in TCM clinics. In the Compendium of Materia Medica, it is recorded that Eucommia can nourish the liver and kidney and treat waist and knee pain. Eucommiae Cortex 60% ethanol extract can remarkably inhibit the osteopenia caused by tail suspension, protect the microstructure of femoral trabecular bone, and improve the biomechanical properties of rat femur [6]. In particular, Eucommiae cortex could induce osteogenic differentiation of bone marrow mesenchymal stem cells through Wnt/β-catenin, MAPK, and Rho A/ROCK signaling pathways [7]. Eucommiae Cortex has a variety of active components and complex mechanisms of action. It is indispensable for screening and validation before Eucommiae Cortex is considered a latent medication. Network pharmacology can establish systematic networks of drug components, targets, and diseases at different levels for screening. The research perspective of systems biology and the network pharmacology methods are more aligned with the characteristics of natural herbal medicine.
Based on network pharmacology and cell experiments, this study explored the important mechanism and targets of Eucommiae Cortex against osteoporosis. After obtaining the targets of Eucommiae Cortex and osteoporosis in authoritative public databases, we constructed a “component-target” regulatory network. The primary pathways are analyzed with topology analysis and biological function analysis. Subsequently, in silico and in vitro experiments, we verified that relative components and core target genes of the primary pathway play significant roles in the proliferation and differentiation of osteoblasts.
2. Materials and Methods
2.1. Acquisition of Active Drug Ingredients and Their Targets
TCMSP platform is a platform for drug screening and evaluation of herbal medicine systems. It provides information on essential characteristics of ADME, such as half-life, oral bioavailability (OB), drug-like properties (DL), and blood-brain barrier. In previous studies, OB and DL were two crucial criteria for drug screening [8]. Through the TCMSP platform, “du Zhong” and “Eucommiae Cortex”were used as the keywords for retrieval. The limiting conditions were set as and , and obtained the information of all the chemical components of Eucommiae Cortex. The active ingredient targets were obtained from the Drugbank database and published literature. In addition, we searched the protein database ID from the Uniprot database.
2.2. Acquisition of Potential Therapeutic Targets of Eucommiae Cortex for Osteoporosis
The keyword “Osteoporosis” was sought in the GeneCards database and online human Mendelian inheritance (OMIM) database. Further, obtained common target genes of drugs and diseases serve as potential therapeutic targets for antiosteoporosis.
2.3. Construction of a “Component-Target” Regulatory Network
The active ingredients and potential antiosteoporosis targets were imported into Cytoscape 3.7.2 software to build a “component-target” regulatory network. This network can guide each compound and target to find the corresponding match.
2.4. Protein-Protein Interaction (PPI) Network and Topology Analysis
Potential therapeutic targets were retrieved from the STRING database and obtained the interaction relationships. The protein interaction network was drawn by Cytoscape software. The TSV files exported in the STRING database are imported into the Cytonca plugin. A total of six degrees are available for comprehensive screening of target genes, including betweenness (BC), clotonca, degree (DC), eigenvector (EC), local average connectivity- (LAC-) based method, and network (NC) [9]. The six-degree values represent the influence of nodes at different levels and choose the median value of each degree as the filtering condition of topology analysis.
2.5. Analysis of Biological Processes and Pathways of Core Targets
R software was used to study the core targets of Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway analysis and Gene Ontology (GO) analysis, with a value <0.05 as the screening condition. Bar chart and bubble chart show the top 20 results.
2.6. Molecular Docking
The compounds’ structures were obtained from PubChem’s bioassay database. The compounds were executed to energy minimization using the software Chem3D (version 17.0) and saved in the PDB format. We set rotatable keys of the ligands and saved them as a PDBQT file for semiflexible docking using the software AutoDockTools-1.5.6. We downloaded the structures in the PDB format of the targets from the Uniport database. The receptor proteins were processed by Pymol software to delete the ligands bound initially to the receptors. We selected the docking region of the proteins, including the central site (X, Y, and Z) and the docking size on the X, Y, and Z axis. To display the docking results, open the receptor PDB file and the ligand PDBQT file by Pymol software simultaneously. AutoDock Vina was used to conduct molecular docking.
2.7. Cell Experiment
2.7.1. Materials and Samples
Cell Counting Kit-8 (CCK-8) was purchased from Beijing Solar-bio Science & Technology Co., Ltd. (Beijing, China). ALP kit was bought from Jiancheng Bioengineering Institute (Nanjing, China). MC3T3-E1 cells were purchased from the Cell Culture Collection Committee of the Chinese Academy of Sciences (Shanghai, China). DMEM and FBS were purchased from Life-iLab (Shanghai, China). TRIzol reagent was ordered from Invitrogen (Carlsbad, CA, USA). SYBR Green Pro Tap was brought from AG (Shanghai, China). Primer was synthesized by Sangong Bioengineering (Shanghai, China). Quercetin, kaempferol, and beta-carotene were ordered from Sigma-Aldrich, Merck KGaA (Darmstadt, Germany). The molecular formula is shown in Figure 1.
(a)
(b)
(c)
2.7.2. Cell Culture
MC3T3-E1 cells, established from C57BL/6 mouse calvaria, were outsourced from the Cell Culture Collection Committee of the Chinese Academy of Sciences (Shanghai, China). These cells were cultured in a DMEM medium (Gibco, Grand Island, NY, USA), supplemented with 10% fetal bovine serum (FBS, Gibco, Grand Island, NY, USA) at 37°C in a chamber with 5% CO2.
2.7.3. CCK-8 and ALP Activity Assay
As for CCK-8, MC3T3-E1 cells were cultured for 24 h in 96-well plates at a density of cells/well, and then treated with quercetin at concentrations of 10−9, 10−8, 10−7, 10−6, 10-5, and 0 mol/L for 24, 48, and72 h in DMEM. Then, 100 μL fresh medium containing 10 μL CCK-8 solution was added to each well, and the cells were incubated at 37°C for one h. The study was using a microplate reader to measure the absorbance at 450 nm.
As for the ALP activity assay, MC3T3-E1 cells were cultured for 24 h in 24-well plates at a density of cells/well and then treated with quercetin at concentrations of 10−9, 10−8, 10−7, 10−6, 10-5, and 0 mol/L for 24, 72, and 120 h in DMEM. Each hole cell was gently washed thrice with phosphate-buffered saline and then lysed with 500 μL 0.1% Triton X-100. The ALP activity assay kit was used to determine ALP activity in the lysate. The micromoles of p-nitrophenol liberated per-nanogram protein expressed the ALP activity. Kaempferol and beta-carotene adopt the same strategy.
2.7.4. Total RNA Isolation and Quantitative Real-Time PCR Analysis
Total RNA was isolated from osteoblast cell cultures using TRIzol reagent. Then, RNA samples were reverse-transcribed into cDNA. cDNA is used as the template for qRT-PCR analysis of target genes using SYBR Green Pro Tap (AG, China). Specific primers are listed in Table 1. The following parameters for PCR amplification were used: 94°C 30 sec, 30 cycles of 98°C for 10 sec, 55°C for 30 sec, and 72°C for 60 sec, and a final extension step of 72°C 120 sec. Relative mRNA expression levels were calculated by the 2-ΔΔCt method. Each sample was assayed at least five times.
2.8. Statistical Analysis
Statistical analysis was processed with Prism 8 software and R4.1.0software. Data were expressed as the and analyzed using Student t-test or one-way ANOVA. Differences between groups were considered statistically significant if values of .
3. Results
3.1. Screening of Active Ingredients and Target Genes
In the TCMSP database, Eucommiae Cortex contains 28 active ingredients in total, with and as the conditions (Table 2). Then, the DrugBank database and Uniprot database predicted the targets. Eventually, a total of 181 drug targets were obtained.
3.2. Eucommiae Cortex Potential Therapeutic Targets
A total of 932 targets of osteoporosis were obtained in the Genecards database and OMIM database. After comparison with drug targets, 85 intersection targets were derived (Table 3, Figure 2(a)).
(a)
(b)
(c)
(d)
3.3. A “Component-Target” Regulatory Network
The “component-target” regulatory network represents the correspondence between active compounds and antiosteoporosis targets (Figure 2(b)). Among them, quercetin and kaempferol have the most related marks.
3.4. PPI Network and Topological Analysis
PPI network showed the interaction between potential therapeutic targets, with 84 nodes (one of the targets is hidden because of no protein interaction) and 1516 edges (Figure 2(c)). A subsequent topology analysis filtered 84 nodes. The median value of BC, CC, DC, LAC, and NC are taken for filtering (BC: 21.594879345; CC: 0.6264240145; DC: 34; EC: 0.0984259805; LAC: 26.26559715; NC: 28.311158725). 36 nodes and 563 edges were obtained. Then, take the median filter again (BC: 3.2453257765; CC: 0.897435897; DC: 31; EC: 0.1647364195; LAC: 27.51612903; NC: 29.57936678); 17 nodes and 136 edges were eventually identified (Figure 2(d)). According to the “component-target” regulatory network, quercetin, kaempferol, and beta-carotene were most associated with 17 core targets, which were selected for subsequent verification.
3.5. GO and KEGG Analysis
GO analysis illustrates gene function on three levels: biological process (BP), cellular component (CC), and molecular function (MF) (Figure 3(a)). BP mainly involves cellular response to chemical stress, oxidative stress, and metal ion. CC is mainly related to the vesicle lumen, platelet alpha granule lumen, and secretory granule lumen. MF is involved primarily in cytokine receptor, growth factor receptor, activating transcription factor binding, MAP kinase activity, and R−SMAD binding. According to KEGG enrichment results, the mechanisms of Eucommiae Cortex of antiosteoporosis are mainly concentrated in Kaposi sarcoma-associated herpesvirus infection, MAPK signaling pathway, and hepatitis B signaling pathway (Figure 3(b)). Combined with GO, KEGG results, and disease-related biological knowledge, the MAPK signaling pathway was most valuable in developing osteoporosis. Ten genes from the MAPK signaling pathway were selected (Figure 3(c)).
(a)
(b)
(c)
3.6. Molecular Docking
Molecular docking mimics the interaction of small ligands with large receptor protein molecules. The affinity can be predicted by calculating the binding energy. Smaller binding energy indicates a more stable structure. Quercetin, kaempferol, and beta-carotene have good binding potential with MAPK pathway target genes (Table 4, Figure 4).
(a)
(b)
(c)
(d)
(e)
(f)
3.7. The Proliferation of MC3T3-E1 Cells
As shown in Figures 5(a)–5(c), quercetin, kaempferol, and beta-carotene could vastly boost the proliferation in MC3T3-E1 cells. Moreover, it is proportional to time. When the concentration of quercetin, kaempferol, and beta-carotene was mol/L, mol/L, and mol/L, respectively, the maximum proliferation capacity was obtained.
(a)
(b)
(c)
(d)
(e)
(f)
3.8. The ALP Activity of MC3T3-E1 Cells
ALP is a glycoprotein secreted by osteoblasts, indirectly reflecting differentiation and maturity. As shown in Figures 5(d)–5(f), at 24, 48, and 72 h, the activity of ALP increased with time. As for quercetin, the concentration of reached the ALP activity peak, and differentiation was the highest. At 10-5, it was flat or slightly decreased. The optimal concentration of was selected. As for kaempferol and beta-carotene, the concentration gradient of is the optimum selection.
3.9. qPCR Validation of the Core Genes
The CCK-8 and ALP activity test results showed quercetin at mol/L, kaempferol at mol/L, and beta-carotene at mol/L be selected to verify. As shown in Figure 6, the addition of the three compounds amplified the transcription factor complex (Jun/FOS), thus playing a role in promoting proliferation and differentiation. Quercetin and beta-carotene mainly passed through the MAPK1/ERK2 pathway, while kaempferol passed through the MAPK8/JNK1 pathway.
(a)
(b)
(c)
4. Discussion
After entering old age, the body will slow down the production of new bone tissue. The etiology and pathogenesis of osteoporosis are frustratingly complex. Some of the leading causes may be the low level of sex hormones, the lack of differentiation of bone marrow stromal cells into osteoblasts, and the breaking of the dynamic balance between osteoclast bone resorption and osteoblast bone remodeling [10]. A recent review highlighted type I collagen processing, Wnt signaling, TGF-β signal transduction, RANKL rank system, and mechanical injury of bone cells [11]. Our study found that the MAPK pathway has substantial value in treating osteoporosis, particularly related to the pharmacodynamic effect of Eucommiae Cortex.
Network pharmacology combines molecular biology, genetics, and computer science to collect, analyze, mine, and deploy biological information. Behind diseases and natural drugs backgrounds, this technology can better predict and classify various therapeutical targets and interactions [8]. In previous studies, Eucommiae Cortex presented a significant bone protective effect on rats induced by lead acetate, promoting bone formation and inhibiting bone resorption [12]. However, the specific mechanisms and targets of action remain unclear. We explore 85 potential therapeutical targets through the “component-target” regulatory network, including apoptosis-related protein, growth factor, inflammatory factors, and protein kinase. Further, PPI analysis and topology analysis screened 17 core targets. According to the enrichment of biological functions and the supplement of existing literature, core targets are widely involved in the host-virus response, chemical stress, oxidative stress, and cell cycle control. KEGG analysis also demonstrated that the antiosteoporosis effect of Eucommiae Cortex was connected to virus infection and cell signal transduction, among which the MAPK pathway was the most valuable. There are three active compounds that correlated with the MAPK pathway—quercetin, kaempferol, and beta-carotene.
Mitogen-activated protein kinase (MAPK) cascades are essential signaling pathways that regulate various cellular processes, including proliferation, differentiation, and apoptosis [13]. According to our analysis results, Eucommiae Cortex is mainly related to the two branches, extracellular signal-related kinase (ERK) and Jun amino-terminal kinase (JNK). ERK signal is an indispensable member of cell development. It responds to various stimuli such as aging, trauma, and infection and drives biochemical activities for organ regeneration [14]. Growth factors can activate the ERK2/MAPK1 to induce cell proliferation and differentiation [15]. Our qPCR experiments reflected that quercetin simultaneously activated pro-epidermal growth factor (EGF), vascular endothelial growth factor A (VEGFA), and ERK2 genes. While beta-carotene only activated VEGFA and ERK2, JNK activity is crucial in the late differentiation of osteoblasts [16]. According to transcriptomic analysis, JNK1/MAPK8 was desired to up-regulate several osteoblast-derived angiogenic factors and was the critical mediator of osteoblast function [17]. Only kaempferol-stimulated osteoblasts showed JNK1 up-regulation.
Activator protein 1 (AP-1) transcription factor is a dimer complex composed of Jun and Fos proteins [18]. It is the downstream executive element of the MAPK pathway. AP-1 acts as a valid regulator in osteoblast differentiation and proliferation; under physiological conditions, the activity of AP-1 is determined by parathyroid hormone and transforming growth factor-β and 1,25-dihydroxy vitamin D [19]. Rainer Zenz [20] reviewed the capacity of AP-1 in bone immunological therapy. The therapeutic interventions on AP-1 signal transduction might provide potent weapons for treating low bone mass diseases. After three compounds stimulated osteoblasts, c-jun and c-fos were amplified in our study. Binding properties in vitro also demonstrated that the three compounds have superior binding force with Jun and FOS. This research provides indirect evidence for the study of targeting AP-1 transcription factors.
Quercetin and kaempferol are natural flavonoids found in a variety of phytomedicine. They participate in bone remodeling through the methods of inhibiting adipogenesis, inflammation, oxidative stress, and osteoblast apoptosis [21]. Beta-carotene is likewise salutary for maintaining bone health. Recent clinical studies have shown that higher beta-carotene intake was associated with a lower risk of osteoporosis and fractures [22]. Beta-carotene treatment also increased the expression of osteopontin and ALP, markers of bone differentiation [23]. Although the values of the three compounds as osteoprotective agents have received some attention, their effects on osteoblast proliferation and differentiation have not been carefully evaluated. Our research has proved that quercetin, kaempferol, and carotene can promote the proliferation and differentiation of osteoblasts. After being stimulated by the three compounds, the proliferation rate and ALP activity of MC3T3-E1 cells increased gradually at 24, 48, and 72 hours. It is worth mentioning that all three compounds reduced the expression of apoptosis factor CASP3 in MC3T3-E1 cells. Therefore, we speculate that the three compounds can inhibit osteoblast apoptosis.
In summary, we found 17 core targets, and the MAPK signaling pathway contributed to Eucommiae Cortex in treating osteoporosis. Subsequently, cell experiments and qPCR detections confirmed that active components (quercetin, kaempferol, and beta-carotene) promote the proliferation and differentiation of MC3T3-E1 cells. Molecular docking simulates the binding state at the computer level. Based on the above evidence, we speculate that Eucommiae Cortex can be used as a potential drug to treat osteoporosis. Our research experiment still has many shortcomings, but it can provide a basis for developing therapeutic drugs for low bone mass diseases and bone repair.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there is no conflict for the publication of this paper.