# A Network Pharmacology-Based Study on Vital Pharmacological Pathways and Targets of Eucommiae Cortex Acting on Osteoporosis

Authors: Libo Zhou; Tao Wu
Journal: BioMed Research International
 (2022)

---

## 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 1∗10−6 mol/L quercetin, 1∗10−5 mol/L kaempferol, and 1∗10−5 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.


---

## Body

## 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 theCompendium 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 OB>30% and DL>0.18, 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 aP 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 Figure1.Figure 1
Chemical formula of (a) quercetin, (b) kaempferol, and (c) beta-carotene.
(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 of5×103 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 of2×104 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 Table1. 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.Table 1
Quantitative real-time PCR primer sequences.
GeneFORREVEGF5′-TTGGTATGCAAGGATGTGTC-3′5′-CCACTTTGCGAAGTAACTTGGTA-3′EGFR5′-GCCATCTGGGCCAAAGATACC-3′5′-GTCTTCGCATGAATAGGCCAAT-3′VEGFA5′-CCCGGGCCTCGGTTCCAG-3′5′-GTCGTGGGTGCAGCCTGGG-3′MAPK1/ERK25′-GCTGAATCACATCCTGGGTA-3′5′-CTGTTCCACGGGCACCTTATT-3′MAPK8/JNK15′-GTGGAATCAAGCACCTTCACT-3′5′-TCCTCGCCAGTCCAAAATCAA-3′CASP35′-GAAATTGTGGAATTGATGCGTGA-3′5′-CTACAACGATCCCCTCTGAAAAA-3′AKT15′-GGACTACTTGCACTCCGAGAA-3′5′-CATAGTGGCACCGTCCTTGATC-3′c-JUN5′-ACTCG-GACCTCCTCACCTCG-3′5′-TGTTTAAGCTGTGCCACCTGTT-3′FOS5′-GCGAGCAACTGAGAAGAC-3′5′-TTGAAACCCGAGAACATC-3′c-MYC5′-AACTTACAATCTGCGAGCCA-3′5′-AGCAGCTCGAATTTCTTCCAGATAT-3′β-Actin5′-GTACTCTGTGTGGATCGGTGG-3′5′-AACGCAGCTCAGTAACAGTCC-3′


### 2.8. Statistical Analysis

Statistical analysis was processed with Prism 8 software and R4.1.0software. Data were expressed as themean±SD and analyzed using Student t-test or one-way ANOVA. Differences between groups were considered statistically significant if values of P<0.05.


## 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 OB>30% and DL>0.18, 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 aP 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 Figure1.Figure 1
Chemical formula of (a) quercetin, (b) kaempferol, and (c) beta-carotene.
(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 of5×103 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 of2×104 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 Table1. 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.Table 1
Quantitative real-time PCR primer sequences.
GeneFORREVEGF5′-TTGGTATGCAAGGATGTGTC-3′5′-CCACTTTGCGAAGTAACTTGGTA-3′EGFR5′-GCCATCTGGGCCAAAGATACC-3′5′-GTCTTCGCATGAATAGGCCAAT-3′VEGFA5′-CCCGGGCCTCGGTTCCAG-3′5′-GTCGTGGGTGCAGCCTGGG-3′MAPK1/ERK25′-GCTGAATCACATCCTGGGTA-3′5′-CTGTTCCACGGGCACCTTATT-3′MAPK8/JNK15′-GTGGAATCAAGCACCTTCACT-3′5′-TCCTCGCCAGTCCAAAATCAA-3′CASP35′-GAAATTGTGGAATTGATGCGTGA-3′5′-CTACAACGATCCCCTCTGAAAAA-3′AKT15′-GGACTACTTGCACTCCGAGAA-3′5′-CATAGTGGCACCGTCCTTGATC-3′c-JUN5′-ACTCG-GACCTCCTCACCTCG-3′5′-TGTTTAAGCTGTGCCACCTGTT-3′FOS5′-GCGAGCAACTGAGAAGAC-3′5′-TTGAAACCCGAGAACATC-3′c-MYC5′-AACTTACAATCTGCGAGCCA-3′5′-AGCAGCTCGAATTTCTTCCAGATAT-3′β-Actin5′-GTACTCTGTGTGGATCGGTGG-3′5′-AACGCAGCTCAGTAACAGTCC-3′


## 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 Figure1.Figure 1
Chemical formula of (a) quercetin, (b) kaempferol, and (c) beta-carotene.
(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 of5×103 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 of2×104 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 Table1. 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.Table 1
Quantitative real-time PCR primer sequences.
GeneFORREVEGF5′-TTGGTATGCAAGGATGTGTC-3′5′-CCACTTTGCGAAGTAACTTGGTA-3′EGFR5′-GCCATCTGGGCCAAAGATACC-3′5′-GTCTTCGCATGAATAGGCCAAT-3′VEGFA5′-CCCGGGCCTCGGTTCCAG-3′5′-GTCGTGGGTGCAGCCTGGG-3′MAPK1/ERK25′-GCTGAATCACATCCTGGGTA-3′5′-CTGTTCCACGGGCACCTTATT-3′MAPK8/JNK15′-GTGGAATCAAGCACCTTCACT-3′5′-TCCTCGCCAGTCCAAAATCAA-3′CASP35′-GAAATTGTGGAATTGATGCGTGA-3′5′-CTACAACGATCCCCTCTGAAAAA-3′AKT15′-GGACTACTTGCACTCCGAGAA-3′5′-CATAGTGGCACCGTCCTTGATC-3′c-JUN5′-ACTCG-GACCTCCTCACCTCG-3′5′-TGTTTAAGCTGTGCCACCTGTT-3′FOS5′-GCGAGCAACTGAGAAGAC-3′5′-TTGAAACCCGAGAACATC-3′c-MYC5′-AACTTACAATCTGCGAGCCA-3′5′-AGCAGCTCGAATTTCTTCCAGATAT-3′β-Actin5′-GTACTCTGTGTGGATCGGTGG-3′5′-AACGCAGCTCAGTAACAGTCC-3′


## 2.8. Statistical Analysis

Statistical analysis was processed with Prism 8 software and R4.1.0software. Data were expressed as themean±SD and analyzed using Student t-test or one-way ANOVA. Differences between groups were considered statistically significant if values of P<0.05.


## 3. Results

### 3.1. Screening of Active Ingredients and Target Genes

In the TCMSP database, Eucommiae Cortex contains 28 active ingredients in total, withOB≥30% and DL≥0.18 as the conditions (Table 2). Then, the DrugBank database and Uniprot database predicted the targets. Eventually, a total of 181 drug targets were obtained.Table 2
The available compounds of Eucommiae Cortex.
MOl IDMolecule nameOB (%)DLMOL00205840957-99-157.20.62MOL000211Mairin55.380.78MOL000358Beta-sitosterol36.910.75MOL000422Kaempferol41.880.24MOL004367Olivil62.230.41MOL000443Erythraline49.180.55MOL005922Acanthoside B43.350.77MOL006709AIDS21463492.430.55MOL0070593-Beta-hydroxymethyllenetanshiquinone32.160.41MOL000073Ent-epicatechin48.960.24MOL007563Yangambin57.530.81MOL009007Eucommin A30.510.85MOL009009(+)-Medioresinol87.190.62MOL009015(-)-Tabernemontanine58.670.61MOL009027Cyclopamine55.420.82MOL009029Dehydrodiconiferyl alcohol 4,gamma’-di-O-beta-D-glucopyanoside_qt51.440.4MOL009030Dehydrodieugenol30.10.24MOL009031Cinchonan-9-al, 6'-methoxy-, (9R)-68.220.4MOL009038GBGB45.580.83MOL009042Helenalin77.010.19MOL009047(+)-Eudesmin33.290.62MOL0090534-[(2S,3R)-5-[(E)-3-Hydroxyprop-1-enyl]-7-methoxy-3-methylol-2,3-dihydrobenzofuran-2-yl]-2-methoxy-phenol50.760.39MOL009055Hirsutin_qt49.810.37MOL009057Liriodendrin_qt53.140.8MOL000098Quercetin46.430.28MOL002773Beta-carotene37.180.58MOL008240(E)-3-[4-[(1R,2R)-2-Hydroxy-2-(4-hydroxy-3-methoxy-phenyl)-1-methylol-ethoxy]-3-methoxy-phenyl]acrolein56.320.36MOL011604Syringetin36.820.37


### 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 (Table3, Figure 2(a)).Table 3
Eucommiae Cortex bioactive compound-related osteoporosis targets.
GENE NAMEProtein nameENTRYAHRAryl hydrocarbon receptorP35869AKT1RAC-alpha serine/threonine-protein kinaseP31749ALBSerum albuminP02768ALOX5Arachidonate 5-lipoxygenaseP09917ARAndrogen receptorP10275BAXApoptosis regulator BAXQ07812BCL2Apoptosis regulator Bcl-2P10415BCL2L1Bcl-2-like protein 1Q07817BIRC5Baculoviral IAP repeat-containing protein 5O15392CASP3Caspase-3P42574CASP8Caspase-8Q14790CASP9Caspase-9P55211CAV1Caveolin-1Q03135CCL2C-C motif chemokine 2P13500CCNB1G2/mitotic-specific cyclin-B1P14635CCND1G1/S-specific cyclin-D1P24385CD40LGCD40 ligandP29965CDKN1ACyclin-dependent kinase inhibitor 1P38936CHRM3Muscarinic acetylcholine receptor M3P20309IKKAInhibitor of nuclear factor kappa-B kinase subunit alphaO15111CLDN4Claudin-4O14493CRPC-reactive proteinP02741CTNNB1Catenin beta-1P35222CXCL2C-X-C motif chemokine 2P19875CXCL8Interleukin-8P10145CYP1A1Cytochrome P450 1A1P04798CYP1A2Cytochrome P450 1A2P05177CYP3A4Cytochrome P450 3A4P08684DUOX2Dual oxidase 2Q9NRD8EGFPro-epidermal growth factorP01133EGFREpidermal growth factor receptorP00533ERBB2Receptor tyrosine-protein kinase erbB-2P04626ERBB3Receptor tyrosine-protein kinase erbB-3P21860F3Tissue factorP13726FOSProto-oncogene c-FosP01100GJA1Gap junction alpha-1 proteinP17302GSK3BGlycogen synthase kinase-3 betaP49841GSTM1Glutathione S-transferase Mu 1P09488GSTP1Glutathione S-transferase PP09211HIF1AHypoxia-inducible factor 1-alphaQ16665HSPB1Heat shock protein beta-1P04792ICAM1Intercellular adhesion molecule 1P05362IFNGInterferon gammaP01579IGF2Insulin-like growth factor IIP01344IKBKBInhibitor of nuclear factor kappa-B kinase subunit betaO14920IL10Interleukin-10P22301IL1AInterleukin-1 alphaP01583IL1BInterleukin-1 betaP01584IL2Interleukin-2P60568IL6Interleukin-6P05231JUNTranscription factor AP-1P05412MAPK1Mitogen-activated protein kinase 1P28482MAPK14Mitogen-activated protein kinase 14Q16539MAPK8Mitogen-activated protein kinase 8P45983MMP1Interstitial collagenaseP03956MMP272 kDa type IV collagenaseP08253MMP3Stromelysin-1P08254MMP9Matrix metalloproteinase-9P14780MPOMyeloperoxidaseP05164MYCMyc proto-oncogene proteinP01106NFE2L2Nuclear factor erythroid 2-related factor 2Q16236NFKBIANF-kappa-B inhibitor alphaP25963NOS2Nitric oxide synthase, inducibleP35228NOS3Nitric oxide synthase, endothelialP29474NQO1NAD(P)H dehydrogenase [quinone] 1P15559ODC1Ornithine decarboxylaseP11926PARP1Poly [ADP-ribose] polymerase 1P09874PGRProgesterone receptorP06401PLATTissue-type plasminogen activatorP00750PLAUUrokinase-type plasminogen activatorP00749PPARGPeroxisome proliferator-activated receptor gammaP37231PTGS1Prostaglandin G/H synthase 1P05979PTGS2Prostaglandin G/H synthase 2P35354RASSF1Ras association domain-containing protein 1Q9NS23RB1Retinoblastoma-associated proteinP06400SELEE-selectinP16581SERPINE1Plasminogen activator inhibitor 1P05121SLC6A4Sodium-dependent serotonin transporterP31645SLPIAntileukoproteinaseP03973SOD1Superoxide dismutase [Cu-Zn]P00441SPP1OsteopontinP10451THBDThrombomodulinP07204TP53Cellular tumor antigen p53P04637VCAM1Vascular cell adhesion protein 1P19320VEGFAVascular endothelial growth factor AP15692Figure 2
(a) Venn diagram. It displays the interset of compounds targets related to osteoporosis targets. (b) “Component-target” regulatory network. The regulatory network shows the matched relationship between the active components and the osteoporosis genes. The blue symbol in the inner ring represents herbal Eucommiae cortex. The rose symbols in the outer circle represent 28 active ingredients. The red symbol represents osteoporosis, and the green signs in the outer ring represent 85 intersection genes. (c).PPI network. (d). Topological analysis.
(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 (Figure2(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 (Figure2(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) (Figure3(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)).Figure 3
(a) The BP, CC, and MF column color symbolizes the enrichment significance based on the correctedP value. The horizontal axis represents the number of genes enriched on each item. (b) The vertical axis of the KEGG bubble diagram represents the enrichment degree according to the corrected P value. Moreover, the horizontal axis shows the gene proportion enriched in each entry. (c) MAPK signaling pathway. The relevant targets of the MAPK signaling pathway are marked red.
(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 (Table4, Figure 4).Table 4
Binding energy of the selected compounds to the core genes.
QuercetinKaempferolBeta-caroteneEGF-8.6-5.4-6.4VEGFA-10.2-7.4-10.2EGFR-9.8-6.5-7.3MAPK1/ERK2-7.9-7.4-8.8C-MYC-7.2-8.9-11C-FOS-10.4-9.2-8.4CASP3-9.6-10.1-9.8AKT1-8.1-8.3-8.7MAPK8/JNK1-6.3-10.2-7.7C-JUN-8.9-8.5-8.3Figure 4
Molecular docking: (a) quercetin-ERK2, (b) quercetin-EGF, (c) quercetin-C-FOS, (d) kaempferol-AKT, (e) kaempferol-CASP3, and (f) kaempferol-MAPK8.
(a)(b)(c)(d)(e)(f)


### 3.7. The Proliferation of MC3T3-E1 Cells

As shown in Figures5(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 1∗10−6 mol/L, 1∗10−5 mol/L, and 1∗10−5 mol/L, respectively, the maximum proliferation capacity was obtained.Figure 5
(a)–(c). CCK-8 detect cell viability: (a) quercetin, (b) kaempferol, and (c) beta-carotene. (d)–(f) ALP activity detect cell differentiation: (d) quercetin, (e) kaempferol, and (f) beta-carotene.n=5; Student t-test vs. 0 mol/L; ∗P<0.05, ∗∗P<0.01, ∗∗∗P<0.001.
(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 Figures5(d)–5(f), at 24, 48, and 72 h, the activity of ALP increased with time. As for quercetin, the concentration of 1∗10−6 reached the ALP activity peak, and differentiation was the highest. At 10-5, it was flat or slightly decreased. The optimal concentration of 1∗10−6 was selected. As for kaempferol and beta-carotene, the concentration gradient of 1∗10−5 is the optimum selection.


### 3.9. qPCR Validation of the Core Genes

The CCK-8 and ALP activity test results showed quercetin at1∗10−6 mol/L, kaempferol at 1∗10−5 mol/L, and beta-carotene at 1∗10−5 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.Figure 6
The effects of quercetin, kaempferol, and beta-carotene treatment on the mRNA expression of osteoblast: (a) quercetin, (b) kaempferol, and (c) beta-carotene.n=5 vs. MC3T3-E1; ns: no significant; ∗P<0.05, ∗∗P<0.01, ∗∗∗P<0.001.
(a)(b)(c)


## 3.1. Screening of Active Ingredients and Target Genes

In the TCMSP database, Eucommiae Cortex contains 28 active ingredients in total, withOB≥30% and DL≥0.18 as the conditions (Table 2). Then, the DrugBank database and Uniprot database predicted the targets. Eventually, a total of 181 drug targets were obtained.Table 2
The available compounds of Eucommiae Cortex.
MOl IDMolecule nameOB (%)DLMOL00205840957-99-157.20.62MOL000211Mairin55.380.78MOL000358Beta-sitosterol36.910.75MOL000422Kaempferol41.880.24MOL004367Olivil62.230.41MOL000443Erythraline49.180.55MOL005922Acanthoside B43.350.77MOL006709AIDS21463492.430.55MOL0070593-Beta-hydroxymethyllenetanshiquinone32.160.41MOL000073Ent-epicatechin48.960.24MOL007563Yangambin57.530.81MOL009007Eucommin A30.510.85MOL009009(+)-Medioresinol87.190.62MOL009015(-)-Tabernemontanine58.670.61MOL009027Cyclopamine55.420.82MOL009029Dehydrodiconiferyl alcohol 4,gamma’-di-O-beta-D-glucopyanoside_qt51.440.4MOL009030Dehydrodieugenol30.10.24MOL009031Cinchonan-9-al, 6'-methoxy-, (9R)-68.220.4MOL009038GBGB45.580.83MOL009042Helenalin77.010.19MOL009047(+)-Eudesmin33.290.62MOL0090534-[(2S,3R)-5-[(E)-3-Hydroxyprop-1-enyl]-7-methoxy-3-methylol-2,3-dihydrobenzofuran-2-yl]-2-methoxy-phenol50.760.39MOL009055Hirsutin_qt49.810.37MOL009057Liriodendrin_qt53.140.8MOL000098Quercetin46.430.28MOL002773Beta-carotene37.180.58MOL008240(E)-3-[4-[(1R,2R)-2-Hydroxy-2-(4-hydroxy-3-methoxy-phenyl)-1-methylol-ethoxy]-3-methoxy-phenyl]acrolein56.320.36MOL011604Syringetin36.820.37


## 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 (Table3, Figure 2(a)).Table 3
Eucommiae Cortex bioactive compound-related osteoporosis targets.
GENE NAMEProtein nameENTRYAHRAryl hydrocarbon receptorP35869AKT1RAC-alpha serine/threonine-protein kinaseP31749ALBSerum albuminP02768ALOX5Arachidonate 5-lipoxygenaseP09917ARAndrogen receptorP10275BAXApoptosis regulator BAXQ07812BCL2Apoptosis regulator Bcl-2P10415BCL2L1Bcl-2-like protein 1Q07817BIRC5Baculoviral IAP repeat-containing protein 5O15392CASP3Caspase-3P42574CASP8Caspase-8Q14790CASP9Caspase-9P55211CAV1Caveolin-1Q03135CCL2C-C motif chemokine 2P13500CCNB1G2/mitotic-specific cyclin-B1P14635CCND1G1/S-specific cyclin-D1P24385CD40LGCD40 ligandP29965CDKN1ACyclin-dependent kinase inhibitor 1P38936CHRM3Muscarinic acetylcholine receptor M3P20309IKKAInhibitor of nuclear factor kappa-B kinase subunit alphaO15111CLDN4Claudin-4O14493CRPC-reactive proteinP02741CTNNB1Catenin beta-1P35222CXCL2C-X-C motif chemokine 2P19875CXCL8Interleukin-8P10145CYP1A1Cytochrome P450 1A1P04798CYP1A2Cytochrome P450 1A2P05177CYP3A4Cytochrome P450 3A4P08684DUOX2Dual oxidase 2Q9NRD8EGFPro-epidermal growth factorP01133EGFREpidermal growth factor receptorP00533ERBB2Receptor tyrosine-protein kinase erbB-2P04626ERBB3Receptor tyrosine-protein kinase erbB-3P21860F3Tissue factorP13726FOSProto-oncogene c-FosP01100GJA1Gap junction alpha-1 proteinP17302GSK3BGlycogen synthase kinase-3 betaP49841GSTM1Glutathione S-transferase Mu 1P09488GSTP1Glutathione S-transferase PP09211HIF1AHypoxia-inducible factor 1-alphaQ16665HSPB1Heat shock protein beta-1P04792ICAM1Intercellular adhesion molecule 1P05362IFNGInterferon gammaP01579IGF2Insulin-like growth factor IIP01344IKBKBInhibitor of nuclear factor kappa-B kinase subunit betaO14920IL10Interleukin-10P22301IL1AInterleukin-1 alphaP01583IL1BInterleukin-1 betaP01584IL2Interleukin-2P60568IL6Interleukin-6P05231JUNTranscription factor AP-1P05412MAPK1Mitogen-activated protein kinase 1P28482MAPK14Mitogen-activated protein kinase 14Q16539MAPK8Mitogen-activated protein kinase 8P45983MMP1Interstitial collagenaseP03956MMP272 kDa type IV collagenaseP08253MMP3Stromelysin-1P08254MMP9Matrix metalloproteinase-9P14780MPOMyeloperoxidaseP05164MYCMyc proto-oncogene proteinP01106NFE2L2Nuclear factor erythroid 2-related factor 2Q16236NFKBIANF-kappa-B inhibitor alphaP25963NOS2Nitric oxide synthase, inducibleP35228NOS3Nitric oxide synthase, endothelialP29474NQO1NAD(P)H dehydrogenase [quinone] 1P15559ODC1Ornithine decarboxylaseP11926PARP1Poly [ADP-ribose] polymerase 1P09874PGRProgesterone receptorP06401PLATTissue-type plasminogen activatorP00750PLAUUrokinase-type plasminogen activatorP00749PPARGPeroxisome proliferator-activated receptor gammaP37231PTGS1Prostaglandin G/H synthase 1P05979PTGS2Prostaglandin G/H synthase 2P35354RASSF1Ras association domain-containing protein 1Q9NS23RB1Retinoblastoma-associated proteinP06400SELEE-selectinP16581SERPINE1Plasminogen activator inhibitor 1P05121SLC6A4Sodium-dependent serotonin transporterP31645SLPIAntileukoproteinaseP03973SOD1Superoxide dismutase [Cu-Zn]P00441SPP1OsteopontinP10451THBDThrombomodulinP07204TP53Cellular tumor antigen p53P04637VCAM1Vascular cell adhesion protein 1P19320VEGFAVascular endothelial growth factor AP15692Figure 2
(a) Venn diagram. It displays the interset of compounds targets related to osteoporosis targets. (b) “Component-target” regulatory network. The regulatory network shows the matched relationship between the active components and the osteoporosis genes. The blue symbol in the inner ring represents herbal Eucommiae cortex. The rose symbols in the outer circle represent 28 active ingredients. The red symbol represents osteoporosis, and the green signs in the outer ring represent 85 intersection genes. (c).PPI network. (d). Topological analysis.
(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 (Figure2(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 (Figure2(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) (Figure3(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)).Figure 3
(a) The BP, CC, and MF column color symbolizes the enrichment significance based on the correctedP value. The horizontal axis represents the number of genes enriched on each item. (b) The vertical axis of the KEGG bubble diagram represents the enrichment degree according to the corrected P value. Moreover, the horizontal axis shows the gene proportion enriched in each entry. (c) MAPK signaling pathway. The relevant targets of the MAPK signaling pathway are marked red.
(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 (Table4, Figure 4).Table 4
Binding energy of the selected compounds to the core genes.
QuercetinKaempferolBeta-caroteneEGF-8.6-5.4-6.4VEGFA-10.2-7.4-10.2EGFR-9.8-6.5-7.3MAPK1/ERK2-7.9-7.4-8.8C-MYC-7.2-8.9-11C-FOS-10.4-9.2-8.4CASP3-9.6-10.1-9.8AKT1-8.1-8.3-8.7MAPK8/JNK1-6.3-10.2-7.7C-JUN-8.9-8.5-8.3Figure 4
Molecular docking: (a) quercetin-ERK2, (b) quercetin-EGF, (c) quercetin-C-FOS, (d) kaempferol-AKT, (e) kaempferol-CASP3, and (f) kaempferol-MAPK8.
(a)(b)(c)(d)(e)(f)


## 3.7. The Proliferation of MC3T3-E1 Cells

As shown in Figures5(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 1∗10−6 mol/L, 1∗10−5 mol/L, and 1∗10−5 mol/L, respectively, the maximum proliferation capacity was obtained.Figure 5
(a)–(c). CCK-8 detect cell viability: (a) quercetin, (b) kaempferol, and (c) beta-carotene. (d)–(f) ALP activity detect cell differentiation: (d) quercetin, (e) kaempferol, and (f) beta-carotene.n=5; Student t-test vs. 0 mol/L; ∗P<0.05, ∗∗P<0.01, ∗∗∗P<0.001.
(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 Figures5(d)–5(f), at 24, 48, and 72 h, the activity of ALP increased with time. As for quercetin, the concentration of 1∗10−6 reached the ALP activity peak, and differentiation was the highest. At 10-5, it was flat or slightly decreased. The optimal concentration of 1∗10−6 was selected. As for kaempferol and beta-carotene, the concentration gradient of 1∗10−5 is the optimum selection.


## 3.9. qPCR Validation of the Core Genes

The CCK-8 and ALP activity test results showed quercetin at1∗10−6 mol/L, kaempferol at 1∗10−5 mol/L, and beta-carotene at 1∗10−5 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.Figure 6
The effects of quercetin, kaempferol, and beta-carotene treatment on the mRNA expression of osteoblast: (a) quercetin, (b) kaempferol, and (c) beta-carotene.n=5 vs. MC3T3-E1; ns: no significant; ∗P<0.05, ∗∗P<0.01, ∗∗∗P<0.001.
(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.



---

*Source: 8510842-2022-03-31.xml*