Identified the Synergistic Mechanism of <i>Drynariae Rhizoma</i> for Treating Fracture Based on Network Pharmacology

Lin, Haixiong; Wang, Xiaotong; Wang, Ligang; Dong, Hang; Huang, Peizhen; Cai, Qunbin; Mo, Yingjie; Huang, Feng; Jiang, Ziwei

doi:https://doi.org/10.1155/2019/7342635

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 7342635 | https://doi.org/10.1155/2019/7342635

Identified the Synergistic Mechanism of Drynariae Rhizoma for Treating Fracture Based on Network Pharmacology

Haixiong Lin,¹Xiaotong Wang,²Ligang Wang,³Hang Dong,⁴Peizhen Huang,⁴Qunbin Cai,⁴Yingjie Mo,⁵Feng Huang,^1,4and Ziwei Jiang⁴

Academic Editor: Simona Martinotti

Received22 Jun 2019

Revised14 Sept 2019

Accepted20 Sept 2019

Published20 Oct 2019

Abstract

Background. Drynariae Rhizoma (DR) has been widely used in the prevention and treatment of various fractures. However, the specific mechanisms of DR’s active ingredients have not been elucidated. The purpose of this study was to explore the synergistic mechanisms of DR for treating fracture. Methods. A network pharmacology approach integrating ingredient screening, target exploration, active ingredients-gene target network construction, protein-protein interaction network construction, molecular docking, gene-protein classification, gene ontology (GO) functional analysis, KEGG pathway enrichment analysis, and signaling pathway integration was used. Results. This approach identified 17 active ingredients of DR, interacting with 144 common gene targets and 143 protein targets of DR and fracture. NCOA1, GSK3B, TTPA, and MAPK1 were identified as important gene targets. Five most important protein targets were also identified, including MAPK1, SRC, HRAS, RXRA, and NCOA1. Molecular docking found that DR has a good binding potential with common protein targets. GO functional analysis indicated that common genes involve multiple processes, parts and functions in biological process, cellular component, and molecular function, including positive regulation of transcription from RNA polymerase II promoter, signal transduction, cytosol, extracellular exosome, cytoplasm, and protein binding. The KEGG pathway enrichment analysis indicated that common gene targets play a role in repairing fractures in multiple signaling pathways, including MAPK, PI3K/AKT, Ras, and VEGF signaling pathways. MAPK and PI3K/AKT signaling pathways were involved in osteoblast formation, Ras signaling pathway was involved in enhancing mesenchymal stromal cell migration, and VEGF signaling pathway was involved in angiogenesis. Conclusion. The study revealed the correlation between DR and fracture and the potential synergistic mechanism of different targets of DR in the treatment of fractures, which provides a reference for the development of new drugs.

1. Introduction

Fracture is a common and frequent disease that occurs in patients with various injuries or osteoporosis [1]. In China, the population-weighted incidence of traumatic fractures of the legs, arms, or trunk in 2014 was 3.21 per 1,000 people (95% CI 2.83–3.59) [2]. Osteoporotic fractures are estimated to account for half of all fractures by 2050, and the estimated cost of osteoporotic hip fractures worldwide may reach $131 billion [3]. Therefore, the study of drugs for the prevention and treatment of fractures plays an important role in promoting patient health and reducing family economic pressure.

Recently, DR, one of the plants from Davalliaceae and Davallia Sm., has been widely used in the prevention and treatment of various fractures due to excellent treatment, low side effects, extensive use, and safety [4]. Animal experiments have confirmed that DR could alter the bone histomorphology and increase the number of trabeculae by 10% [5], and its osteogenesis is related to Runx2 and BMP-2 signaling pathways [6]. In addition, it is believed that the various ingredients contained in an herb could regulate multiple targets in different signaling pathways and produce synergistic therapeutic effects [7]. However, such research has not been carried out in the treatment of fractures with DR.

Network pharmacology based on systems biology and polypharmacology has achieved a paradigm shift from “one drug, one goal” to “multi-ingredient therapy, biological network,” which has attracted the attention of Chinese medicine researchers and has been recognized as an effective tool for elucidating multiple components, targets, synergistic effects, and mechanisms of Chinese medicine [8–10]. It is reported that network pharmacology predicts the clinical efficacy, pathways, and side effects of drugs by constructing drug-drug networks, disease-drug networks, and disease-disease networks, providing valuable information for improving the clinical efficacy, reducing toxicity, and elucidating multimechanisms of drugs [11]. For example, Wang Nani found that Er-Xian Decotion has 13 main components closely related to 65 osteoporosis-related targets by using network pharmacology, thereby constructing Er-Xian Decotion component-osteoporosis target network and potential antiosteoporosis mechanism [12]. Yueying et al. identified 108 compounds, 86 potential targets, and 47 signal transduction pathways that Danshiliuhao Granule regulates liver fibrosis by the network pharmacology method, which reflects the multicomponent, multitarget, and multichannel characteristics of Chinese herbal medicine in antiliver fibrosis [13]. Therefore, in order to reveal the relationship between fracture and the active ingredients involved in the DR, we conducted network pharmacology to achieve this goal from protein and gene level. We collected the information of targets from active ingredients in DR and targets of fracture from several databases, respectively, and used network pharmacology to explore the potential synergistic mechanisms of DR for treating fracture.

2. Materials and Methods

2.1. Screening of Active Ingredients of Drynariae Rhizoma

Traditional Chinese Medicine Systems Pharmacology (TCMSP, http://lsp.nwu.edu.cn/, Version 2.3) Database and Analysis Platform includes chemicals, targets, and drug-target-disease networks, as well as pharmacokinetic properties involving oral bioavailability, druglikeness, blood-brain-barrier, and so on [14]. There were 71 compounds of DR which were obtained from the TCMSP. The potential active ingredients of DR for treating fracture were screened according to their oral bioavailability (OB) ≥30% and druglikeness (DL) ≥0.18 recommended by TCMSP.

2.2. Obtaining the Chemical Structure of Active Ingredients

The structure of the potential active ingredients of DR was downloaded from TCMSP and stored in mol2 format. If there was no chemical structure, the PubChem compound was put into the PubChem (https://pubchem.ncbi.nlm.nih.gov/) to download a chemical structure and save it in sdf format, or the PubChem compound was put into the Zinc database (https://zinc.docking.org/) to download a chemical structure and save it in mol2 format. The related SMILES of potential active ingredients was received from TCMSP or PubChem or Zinc database. Then, the SMILES was put into the Swiss Target Prediction database (http://www.swisstargetprediction.ch/) to obtain the related drug target and save it.

2.3. Gene Targets of Drynariae Rhizoma

The DRAR-CPI server (http://cpi.bio-x.cn/drar, update in 2017-7-26) has a collection of drug molecules and targetable human proteins [15]. When submitting a drug molecule, the server docks the drug uploaded by users with the three-dimensional structure of all protein targets in the database, scores, and ranks them with the affinity scoring function based on the protein-ligand interaction, thereby predicting the potential protein targets of human-targetable drugs [15, 16]. This affinity score is called Z-score in the DRAR-CPI server [17]. Protein-ligand interaction with Z-score <−0.5 was recommended by DRAR-CPI as a potential protein target for human-targetable drugs [16]. We uploaded the potential active ingredients of DR in mol2 or sdf format and used Z-score <−0.5 to select potential protein targets for DR. A total of 1760 proteins with Z-score <−0.5 and 355 protein targets were obtained after deletion of the duplicate data. The PDB ID of the protein targets were inputted into UniProt KB (http://www.uniprot.org/uniprot/) of the UniProt database, and the “popular organisms” was selected as human to obtain the gene targets associated with the potential active ingredient of DR.

2.4. Gene Target Prediction for Drynariae Rhizoma to Treat Fractures

The following electronic databases were searched to identify the genes related to fractures: Genetic Association Database (https://geneticassociationdb.nih.gov/), Therapeutic Targets Database (http://bidd.nus.edu.sg/BIDD-Databases/TTD/TTD.asp), PharmGkb database (https://www.pharmgkb.org/), GeneCards database (http://www.genecards.org/), and OMIM database (http://www.ncbi.nlm.nih.gov/omim). Then, the duplicate data and false-positive genes were deleted. Finally, the Venny tool (http://bioinfogp.cnb.csic.es/tools/venny/index.html, Version 2.1) was used to identify the common gene targets of DR and fracture, which may be the potential targets for DR to treat fractures.

2.5. Constructing the Ingredient-Target Network of Drynariae Rhizoma

The common gene targets of DR and fracture were introduced into the Cytoscape software (Version 3.4.0) to construct an ingredient-target network of Drynariae Rhizoma and analyze the topology properties of the network, including degree, betweenness centrality, and closeness centrality [18]. The degree describes the number of connections to a node in the network, indicating interaction with other nodes in the network. Betweenness centrality measures the proportion of a node between shortest paths among other nodes, suggesting the importance of nodes in maintaining network tightness. Closeness centrality indicates the degree of nodes close to the “center” of the network. A node with high degree, betweenness centrality, and closeness centrality values means that it plays a very important role in the network [18].

2.6. Constructing Protein-Protein Interaction (PPI) of Drynariae Rhizoma

The String database (https://string-db.org/, Version 10.5) is a database containing known and predicted PPIs, which collect and integrate a large number of protein interactions involving 9,643,763 proteins and 1,380,838,440 interactions, including experimental data and interactive prediction data derived from bioinformatic methods [19]. Common gene targets of DR and fracture were imported into the STRING database, and the species were set to humans for PPIs. Then, the highest confidence was set to 0.9 in the minimum required interaction score and the results were updated. The TSV format of the updated results were downloaded. Then, node1, node2, and combined scores were extracted and imported into the Cytoscape software to create a PPI network, and the network was analyzed as follows: Step 1: analyze the topology properties of the network: cytoscape ⟶ tools ⟶ network analyzer ⟶ network analysis ⟶ analyze network, save the CSV format of the network result and extract the degree value. Step 2: create a network map according to the degree: cytoscape ⟶ tool ⟶ network analyzer ⟶ network analysis ⟶ generate style from statistics ⟶ map node size to degree ⟶ map node color to degree and save the PPI network map.

2.7. Molecular Docking

SystemsDock (http://systemsdock.unit.oist.jp, Version 2.0) is a web server for network pharmacology-based prediction and analysis that could be used to illustrate the role of ligands on a complex molecular network [20]. It evaluates the protein-ligand binding potential of molecular docking by combining docking with the intelligence (dock-IN) score. The dock-IN score is the negative logarithm of the experimental dissociation/inhibition constant (pKd/pKi), which ranges from 0 to 10, indicating weak to strong binding [20]. It is believed that the docK-IN score above 4.25 indicates a slight binding potential between the protein and ligand; a value greater than 5.0 indicates a moderate binding potential, and a value greater than 7.0 indicates a strong binding potential [16]. We extracted the top 5 proteins with the highest degree value in the PPI network. The proteins that were recognized by systemsDock docked with the potential active ingredients of DR to receive the dock-IN score. The results were saved, and their dock-IN score was analyzed to assess the binding potentials between the potential active ingredients of DR and protein targets.

2.8. GO Functional Analysis and KEGG Pathway Enrichment Analysis

GO (http://www.geneontology.org) is widely used for annotation of gene function, providing detailed annotations of gene function in terms of biological process (BP), cellular component (CC), and molecular function (MF), respectively [21]. Database for annotation, visualization, and integrated discovery (David, https://david.ncifcrf.gov/, Version 6.8) is a functional genomic annotation database that provides bioinformatics annotation for genes or proteins based on the gene annotation function of the GO database and the signaling pathway information of the KEGG database [22]. We performed GO functional analysis and KEGG pathway enrichment analysis in the David database. The procedure was as follows: Step 1: paste the common gene targets of DR and fracture list. Step 2: select “OFFICIAL_GENE_SYMBOL” in “Select Identifier.” Step 3: select “Gene List” in “List Type.” Step 4: select “Homo sapiens” in species. Step 5: submit list. Step 6: download the results of BP, CC, and MF in the gene ontology. Step 7: download the results of KEGG pathway in the pathways. Step 8: targets with were screened and sorted by count (number of targets), and the top-ranked biological processes or KEGG pathways were extracted. Step 9: BP, CC, and MF were designed using GraphPad Prism 5.0 software. The KEGG pathways were designed by the advanced bubble chart of the omicshare tool (http://omicshare.com/tools/Home/Soft/getsoft/type/index).

2.9. Collect Protein Class Corresponding to Common Gene Targets

DisGeNET (http://www.disgenet.org/web/DisGeNET/menu, Version 5.0) is a discovery platform that contains one of the largest publicly available genes and variants associated with human disease. It could be used to analyze the properties of disease genes and investigate the molecular basis of specific diseases and their comorbidities, as well as adverse drug reactions [23]. We used the search function of the DisGeNET platform to retrieve the protein class corresponding to common gene targets.

2.10. Pathway Integration

We used the KEGG Mapper tool in the KEGG database (http://www.kegg.jp/) to retrieve some pathways of DR for fractures and then integrate into a final pathway map. The procedure was as follows: Step 1: used the UniProt KB search function of the UniProt database to retrieve the UniprotID of the common gene targets. Step 2: import the UniProt ID of the common gene targets. Step 3: set the parameters: search against: hsa, primary ID: NCBI-UniProt ID, and examples: Homo sapiens pathway. Step 4: download the PI3K-AKT, MAPK, Ras, and VEGF signaling pathways. Step 5: integrate the signal path.

3. Results

3.1. Active Ingredients of Drynariae Rhizoma

A total of 71 ingredients of DR were retrieved from TCMSP, and 18 active ingredients were screened according to the biological functions of DR. However, marioside_qt (Molecule ID: MOL009087) was removed because it could not be recognized by the PubChem or Zinc database. The remaining 17 active ingredients are shown in Table 1, including (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one, aureusidin, eriodictyol (flavanone), stigmasterol, beta-sitosterol, kaempferol, naringenin, (+)-catechin, eriodictyol, digallate, luteolin, 22-stigmasten-3-one, cyclolaudenol acetate, cycloartenone, cyclolaudenol, davallioside A_qt, and xanthogalenol.

3.2. Gene Target Prediction

A total of 303 gene targets associated with the potential active ingredients of DR were retrieved in the UniProt database. A total of 3,173 fracture-related genes were received, and 3,054 genes remained after deletion of the duplicate and false-positive genes. Common gene target screening for fracture and DR are shown in Figure 1. A total of 144 common gene targets of DR and fracture were received, indicating the potential targets for DR to treat fractures, as shown in Table 2.

3.3. Ingredient-Target Network of Drynariae Rhizoma

The active ingredients and gene targets of DR was inputted into Cytoscape software to construct the ingredient-target network, as shown in Figure 2. In the network, the pink oval nodes represent the main active ingredients of DR, and the light green rectangle nodes represent the potential gene targets for DR to treat fractures. The line represents the correlation between the active ingredients of DR and the gene targets. There are 161 nodes and 774 lines in the network. An active ingredient could be linked to different gene targets, and a gene target could be linked to different active ingredients, suggesting the multicomponent and multitarget characteristics of DR. The topology properties of active ingredients of DR are shown in Supplementary Table 1. Both cyclolaudenol and cycloartenone were linked to a maximum number of gene targets, with 70 (48.61%) different gene targets. Cyclolaudenol acetate was linked to 54 (37.5%) different gene targets. Xanthogalenol was linked to a minimum number of gene targets for a total of 29 (20.14%). In addition, there were four gene targets with the top four degree values, betweenness centrality, and closeness centrality at the same time (Table 3), which were NCOA1, GSK3B, TTPA, and MAPK1.

3.4. Protein-Protein Interaction of Drynariae Rhizoma

The PPI network of DR is shown in Figure 3. In the network, the node represents the protein, and the size and color of the node represent the value of the degree. The larger the node and the brighter the color (yellow to blue), the greater the value of the degree. The line indicates the association between proteins. Results showed that there were 143 nodes and 315 lines.

Degree in the network indicates the number of proteins that a protein has interacting with. In other words, top-degree protein targets screened in PPI plays a pivotal role in the treatment of fractures with DR. Five important protein targets with top degree of DR were identified in the PPI network and are shown in Table 4. They were MAPK1, SRC, HRAS, RXRA, and NCOA1.

3.5. Molecular Docking

Three important protein targets with top degree of DR were identified by SystemsDock, including SRC, RXRA, and NCOA1. Dock-IN score of these three proteins docked with 17 active ingredients of DR are shown in Table 5. Molecular docking results showed that there were 17 (33.33%) with a dock-IN score greater than 7.0, 24 (47.06%) with a dock-IN score between 7.0 and 5.0, 8 (15.69%) with a dock-IN score between 5.0 and 4.25, and 2 (3.92%) with a dock-IN score less than 4.25.

3.6. Gene Ontology (GO) Functional Analysis and KEGG Pathway Enrichment Analysis

Enriched gene ontology terms for BP, CC, and MF of potential therapeutic fracture targets from the main active ingredients of DR are shown in Figure 4. In the BP (Figure 4(a)), positive regulation of transcription from RNA polymerase II promoter involved 33 (22.92%) potential therapeutic fracture targets, signal transduction involved 30 (20.84%) potential therapeutic fracture targets, negative regulation of transcription from RNA polymerase II promoter involved 20 (13.89%) potential therapeutic fracture targets, positive regulation of transcription and DNA-template involved 19 (13.19%) potential therapeutic fracture targets, and transcription initiation from RNA polymerase II promoter involved 18 (12.5%) potential therapeutic fracture targets. In the CC (Figure 4(b)), cytosol involved 63 (43.75%) potential therapeutic fracture targets, extracellular exosome involved 60 (41.67%) potential therapeutic fracture targets, cytoplasm involved 58 (40.28%) potential therapeutic fracture targets, nucleus involved 56 (38.89%) potential therapeutic fracture targets, and plasma membrane involved 53 (36.81%) potential therapeutic fracture targets. In the MF (Figure 4(c)), protein binding involved 110 (76.39%) potential therapeutic fracture targets, zinc ion binding involved 31 (21.53%) potential therapeutic fracture targets, identical protein binding involved 29 (20.14%) potential therapeutic fracture targets, ATP binding involved 27 (18.75%) potential therapeutic fracture targets, and enzyme binding involved 23 (15.97%) potential therapeutic fracture targets.

(a)

(b)

(c)

Enriched KEGG pathways of potential targets for treating fracture from the main active ingredients of DR are shown in Figure 5. The MAPK signaling pathway was identified as an important signaling pathway involving 17 (11.81%) potential therapeutic fracture targets with . The PI3K-Akt signaling pathway involved 17 (11.81%) potential therapeutic fracture targets, the Rap1 signaling pathway involved 14 (9.72%) potential therapeutic fracture targets, the Ras signaling pathway involved 14 (9.72%) potential therapeutic fracture targets, and the signaling pathways regulating pluripotency of stem cells involved 12 (9.03%) potential therapeutic fracture targets.

3.7. Protein Class Corresponding to Common Gene Targets

The protein class corresponding to potential targets for treating fracture from the main active ingredients of DR is presented in Table 6. The results showed that DR treatment of the fracture process involved a variety of substances, such as signaling molecule, transcription factor, receptor, enzyme modulator, chaperone, cell adhesion molecule, protein (transporter, transfer protein, carrier protein, calcium-binding protein, defense protein, and immune protein), enzyme modulator, and enzymes (oxidoreductase, kinase, phosphatase, hydrolase, ligase, protease, isomerase, lyase, enzyme regulator, and transferase).

3.8. Signaling Pathway Integration

Four pathways associated with the potential targets of DR main active ingredients for treating fracture are presented in Figure 6. The arrow (⟶) indicates the promoting effect, the T-arrows (⊣) indicate the inhibition, and the arrows of different colors represent different signaling pathways. The targets of the signaling pathway were marked as light blue, and the potential targets of DR main active ingredients for treating fracture were marked as dark blue. There were 21 (14.58%) potential targets of main active ingredients of DR for treating fracture in the PI3K-AKT, MAPK, Ras, and VEGF signaling pathways, indicating that the fracture targets play a role in these signaling pathways. In addition, some targets play a role in a variety of signaling pathways, such as Ras, RafB, AKT/PKB, PI3K, ERK, and JNK.

4. Discussion

In order to reveal the relationship between fracture and the active ingredients involved in the DR, we predicted the mechanism of DR treatment fractures by constructing a biological network of interactions between active ingredients and common gene targets and common protein targets from a molecular level. A total of 17 active ingredients of DR were received in our study, including (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one, aureusidin, eriodictyol (flavanone), stigmasterol, beta-sitosterol, kaempferol, naringenin, (+)-catechin, eriodictyol, digallate, luteolin, 22-stigmasten-3-one, cyclolaudenol acetate, cycloartenone, cyclolaudenol, davallioside A_qt, and xanthogalenol. Most of them were polyphenolic compounds, which are also called flavonoids. Flavonoids are considered to be the main active ingredients of DR and have been reported to reduce bone loss in ovariectomized rats [24]. In addition, Kang Suk-Nam finds that the total phenolics and flavonoids of DR are better extracted with 70% ethanol instead of water, and this ethanol extraction method also makes these extracts have higher antioxidant activity and in vitro antiosteoporosis effect [25]. In the ingredient-target network, all active ingredients were also identified to bind well to the fracture gene targets, binding to at least 29 (20.14%) different gene targets. Therefore, the 17 active ingredients of DR may have the effect of reducing bone loss and promoting fracture healing.

In our study, 144 common gene targets of DR and fracture were received, and 774 interactions between the active ingredients of DR and common gene targets were found. Some gene targets have been confirmed by clinical trials or animal experiments. For example, Guimarães et al. found that polymorphisms in the FGFR1 and BMP4 genes were associated with fracture nonunion in patients [26]. And our team's previous study also found that the total flavonoids of DP could promote osteogenesis and mineralization in rats with tibial defects by increasing the gene expression of BMP2, BMP4, BMPR1A, and Smadl [27]. In the ingredient-target network, NCOA1, GSK3B, TTPA, and MAPK1 were identified as important gene targets based on degree values, betweenness centrality, and closeness centrality. Qin et al. found that NCOA1 promotes angiogenesis by upregulating HIF1α- and AP-1-mediated VEGFa transcription [28]. Galli et al. demonstrated by cell experiments that inhibition of GSK3B could increase cytoplasmic availability of b-catenin, thereby enhancing Wnt classical signaling and osteoblastic differentiation [29]. Fujita et al. found that mice deficient in TTPA developed a high bone mass phenotype in vertebrae and long bones due to lower bone resorption [30]. Matsushita et al. confirmed that MAPK1 (also called ERK2) plays an important role in osteoblast differentiation and osteoclastogenesis [31]. These gene targets are involved in vascularization, osteoblast differentiation, and osteoclastogenesis in fracture repair. Besides, we found that one active ingredient can interact with different gene targets, and one gene target can interact with different active ingredients, which is consistent with the modern drug theory of “multi-ingredient, multitarget” [9].

To identify the interactions of proteins corresponding to common genes, we conducted a PPI network. A total of 143 common protein targets for DR and fracture were received, with 315 PPIs. In addition, MAPK1, SRC, HRAS, RXRA, and NCOA1 were identified as the five most important target proteins. Previous studies have found that MAPK1 and SRC could promote proliferation and differentiation of myeloid cells and inhibit apoptosis [32, 33]. Clinical cases have found that elevated levels of fibroblast growth factor 23 in patients with dysplasia are associated with HRAS mutations [34]. RXRA is an essential cofactor in the action of 1,25-dihydroxyvitamin D, and umbilical cord RXRA methylation was inversely related to offspring bone mineral content [35]. Coronnello et al. found that NCOA1 modulate the estrogen effects in bone, and miR-488-5p overexpression reduces NCOA1 protein levels, thereby reducing bone mineral density [36]. These protein targets are associated with bone growth and angiogenesis in fracture repair. At the same time, we docked SRC, RXRA, and NCOA1 with 17 potential active ingredients of DR and found that 41 (80.39%) had moderate binding potential, suggesting that DR could bind well to fracture-related protein targets.

In order to identify the function of the common gene, we performed GO functional analysis on these genes. The results showed that the common gene involves multiple processes, parts and functions in BP, CC, and MF, which was consistent with existing studies about DR and fracture repair. For example, in the BP, 33 (22.92%) gene targets were involved in positive regulation of transcription from the RNA polymerase II promoter, and 30 (20.84%) gene targets were involved in signal transduction. Previous studies have shown that the promoter activates the polymerase to bind precisely to the template DNA and has the specificity of transcription initiation [37]. The RNA polymerase II promoter responsible for mRNA transcription is the largest and most important class of promoters [37]. This provides conditions for DR to initiate osteogenic targets. Besides, some signal transduction genes have been found in experiments. Song Nan found that VEGFR-2 may play a signal transduction role for naringin, one ingredient of DR, to stimulate angiogenesis and promote fracture healing [38]. In the CC, 63 (43.75%) gene targets were involved in cytosol, 60 (41.67%) gene targets were involved in extracellular exosome, and 58 (40.28%) gene targets were involved in cytoplasm. This indicates that the recovery of the fracture requires the support of various components in the cell, which is consistent with previous studies [39]. In the MF, 110 (76.39%) gene targets were involved in protein binding, suggesting that mutual recognition between proteins has good gene regulation conditions. This is consistent with the protein class corresponding to the potential target. These results were further validated in the protein class corresponding to the common gene. In the protein class, all of these common genes have been found to regulate a variety of fracture-related molecules, such as transcription factors, receptors, enzyme regulators, molecular chaperones, cell adhesion molecules, enzyme, and so on.

In order to identify the synergistic mechanism of DR for fracture, we performed KEGG pathway enrichment analysis and summarized some important signaling pathways, which provides direction for future research. In the KEGG pathway enrichment analysis, 17 (11.81%) gene targets were involved in MAPK signaling pathway, 17 (11.81%) gene targets were involved in PI3K-Akt signaling pathway, 14 (9.72%) gene targets were involved in Ras signaling pathway, and 6 (4.17%) gene targets were involved in VEGF signaling pathway, which suggest that common gene targets play a role in repairing fractures in multiple signaling pathways. MAPK and PI3K/AKT signaling pathways have been demonstrated to promote osteoblastic bone formation [40]. Zhang et al. confirmed that total flavonoids from DR promote the osteogentic differentiation of ciliary neurotrophic factor-modified myoblasts by activating p38 MAPK signaling pathway [41]. Moreover, total flavonoids of DR could promote osteogenic differentiation of rat dental pulp stem cells via the PI3K/Akt pathway [42]. Lin et al. found that the effect of naringin on the healing of fracture may be related to the promotion of the synthesis and secretion of cellular chemokines (CXCL5, CXCL6) and enhancement of mesenchymal stromal cell migration through Ras signaling pathway [43]. In addition, naringin stimulates angiogenesis by regulating the VEGF/VEGFR-2 signaling pathway in rats, thereby promoting fracture healing [38]. However, the mechanism of some active ingredients of DR in the treatment of fractures has not yet been verified. Therefore, we integrated MAPK, PI3K/AKT, Ras, and VEGF signaling pathways to provide a reference for researchers to verify the mechanism of other DR active ingredients in the treatment of fractures.

5. Conclusion

We collected the gene and protein targets of fractures and active ingredients of DR and then used network pharmacology to reveal the correlation between drugs and diseases and the potential synergistic mechanism of different targets of DR in the treatment of fractures, which provides a reference for the development of new drugs.

Abbreviations

DR:	Drynariae Rhizoma
GO:	Gene ontology
TCMSP:	Traditional Chinese Medicine Systems Pharmacology
OB:	Oral bioavailability
DL:	Druglikeness
CPI:	Chemical-protein interactome
PPI:	Protein-protein interaction
dock-IN:	Combining docking with intelligence
BP:	Biological process
CC:	Cellular component
MF:	Molecular function
DAVID:	Database for Annotation, Visualization, and Integrated Discovery.

Data Availability

All data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Haixiong Lin and Xiaotong Wang conceived and designed the study. Ziwei Jiang and Feng Huang revised the protocol. Haixiong Lin and Xiaotong Wang extracted the data. Ligang Wang, Hang Dong, Peizhen Huang, Qunbin Cai, and Yingjie Mo checked the data. Haixiong Lin, Xiaotong Wang, and Ligang Wang performed statistical analysis and wrote the manuscript. Haixiong Lin, Xiaotong Wang, Ligang Wang, Hang Dong, Peizhen Huang, Qunbin Cai, Yingjie Mo, Feng Huang, and Ziwei Jiang interpreted the results. Ziwei Jiang and Feng Huang reviewed and proposed advice. All authors contributed to constructive comments on the paper. Haixiong Lin, Xiaotong Wang, and Ligang Wang contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 81774337) and Inheritance Studio Construction Project of Prestigious TCM Doctors (Feng Huang) of Guangdong Province (no. YZYBH[2019]5, Index no. 006939748/2017-00583).

Supplementary Materials

Supplementary Table 1: the topology properties of active ingredients of DR. (Supplementary Materials)

References

G. Zhang, B. Hou, L. Shan, J. Wang, W. Yang, and C. Liang, “A survey on the incidence of hip fractures in middle-aged and old population from Changning District, Shanghai,” Chinese Journal of Tissue Engineering Research, vol. 19, no. 37, pp. 6055–6059, 2015.
View at: Google Scholar
W. Chen, H. Lv, S. Liu et al., “National incidence of traumatic fractures in China: a retrospective survey of 512 187 individuals,” The Lancet Global Health, vol. 5, no. 8, pp. E807–E817, 2017.
View at: Publisher Site | Google Scholar
H.-M. Huang, X.-L. Li, S.-Q. Tu, X.-F. Chen, C.-C. Lu, and L.-H. Jiang, “Effects of roughly focused extracorporeal shock waves therapy on the expressions of bone morphogenetic protein-2 and osteoprotegerin in osteoporotic fracture in rats,” Chinese Medical Journal, vol. 129, no. 21, pp. 2567–2575, 2016.
View at: Publisher Site | Google Scholar
Y. Zhang, J. Jiang, H. Shen, Y. Chai, X. Wei, and Y. Xie, “Total flavonoids from Rhizoma Drynariae (Gusuibu) for treating osteoporotic fractures: implication in clinical practice,” Drug Design, Development and Therapy, vol. 11, pp. 1881–1890, 2017.
View at: Publisher Site | Google Scholar
R. W. K. Wong, B. Rabie, M. Bendeus, and U. Hagg, “The effects of rhizoma curculiginis and rhizoma drynariae extracts on bones,” Chinese Medicine, vol. 2, no. 1, p. 13, 2007.
View at: Publisher Site | Google Scholar
W. Yao, H. Zhang, X. Jiang et al., “Effect of total flavonoids of rhizoma drynariae on tibial dyschondroplasia by regulating BMP-2 and Runx2 expression in chickens,” Frontiers in Pharmacology, vol. 9, p. 1251, 2018.
View at: Publisher Site | Google Scholar
B. Li, J. Rui, X. Ding, and X. Yang, “Exploring the multicomponent synergy mechanism of Banxia Xiexin decoction on irritable bowel syndrome by a systems pharmacology strategy,” Journal of Ethnopharmacology, vol. 233, pp. 158–168, 2019.
View at: Publisher Site | Google Scholar
M. N. Trame, K. Biliouris, L. J. Lesko, and J. T. Mettetal, “Systems pharmacology to predict drug safety in drug development,” European Journal of Pharmaceutical Sciences, vol. 94, no. S1, pp. 93–95, 2016.
View at: Publisher Site | Google Scholar
S. Sheng, J. Wang, L. Wang et al., “Network pharmacology analyses of the antithrombotic pharmacological mechanism of Fufang Xueshuantong capsule with experimental support using disseminated intravascular coagulation rats,” Journal of Ethnopharmacology, vol. 154, no. 3, pp. 735–744, 2014.
View at: Publisher Site | Google Scholar
X.-x. Xu, J.-p. Bi, L. Ping, P. Li, and F. Li, “A network pharmacology approach to determine the synergetic mechanisms of herb couple for treating rheumatic arthritis,” Drug Design, Development and Therapy, vol. 12, pp. 967–979, 2018.
View at: Publisher Site | Google Scholar
A. Liu and G. Du, “Network pharmacology: new guidelines for drug discovery,” Acta Pharmaceutica Sinica, vol. 45, no. 12, pp. 1472–1477, 2010.
View at: Google Scholar
N. Wang, P. Xu, X. Wang et al., “Integrated pathological cell fishing and network pharmacology approach to investigate main active components of Er-Xian decotion for treating osteoporosis,” Journal of Ethnopharmacology, vol. 241, Article ID 111977, 2019.
View at: Publisher Site | Google Scholar
Y. Tao, K. Tian, J. Chen et al., “Network pharmacology-based prediction of the active compounds, potential targets, and signaling pathways involved in Danshiliuhao Granule for treatment of liver fibrosis,” Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 2630357, 14 pages, 2019.
View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014.
View at: Publisher Site | Google Scholar
H. Luo, J. Chen, L. Shi et al., “DRAR-CPI: a server for identifying drug repositioning potential and adverse drug reactions via the chemical-protein interactome,” Nucleic Acids Research, vol. 39, no. 2, pp. W492–W498, 2011.
View at: Publisher Site | Google Scholar
D. Wu, Y. Gao, H. Xiang et al., “Exploration into mechanism of antidepressant of Bupleuri radix based on network pharmacology,” Acta Pharmaceutica Sinica, vol. 53, no. 2, pp. 210–219, 2018.
View at: Google Scholar
T. J. A. Ewing, S. Makino, A. G. Skillman, and I. D. Kuntz, “DOCK 4.0: search strategies for automated molecular docking of flexible molecule databases,” Journal of Computer-Aided Molecular Design, vol. 15, no. 5, pp. 411–428, 2001.
View at: Google Scholar
X. Wang, X. Li, and G. Chen, Complex Network Theory and its Application, Tsinghua University Press, Beijing, China, 2006.
D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic Acids Research, vol. 45, no. D1, pp. D362–D368, 2017.
View at: Publisher Site | Google Scholar
K.-Y. Hsin, Y. Matsuoka, Y. Asai et al., “systemsDock: a web server for network pharmacology-based prediction and analysis,” Nucleic Acids Research, vol. 44, no. W1, pp. W507–W513, 2016.
View at: Publisher Site | Google Scholar
M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
View at: Publisher Site | Google Scholar
G. Dennis, B. T. Sherman, D. A. Hosack et al., “DAVID: database for annotation, visualization, and integrated discovery,” Genome Biology, vol. 4, no. 9, p. R60, 2003.
View at: Publisher Site | Google Scholar
J. Piñero, À. Bravo, N. Queralt-Rosinach et al., “DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants,” Nucleic Acids Research, vol. 45, no. D1, pp. D833–D839, 2017.
View at: Publisher Site | Google Scholar
S. Song, Z. Gao, X. Lei et al., “Total flavonoids of drynariae rhizoma prevent bone loss induced by hindlimb unloading in rats,” Molecules, vol. 22, no. 7, p. 1033, 2017.
View at: Publisher Site | Google Scholar
S.-N. Kang, J. Lee, J.-H. Park et al., “In vitro anti-osteoporosis properties of diverse Korean drynariae rhizoma phenolic extracts,” Nutrients, vol. 6, no. 4, pp. 1737–1751, 2014.
View at: Publisher Site | Google Scholar
J. M. Guimarães, I. C. d. V. Guimarães, M. E. L. Duarte et al., “Polymorphisms in BMP4 and FGFR1 genes are associated with fracture non-union,” Journal of Orthopaedic Research, vol. 31, no. 12, pp. 1971–1979, 2013.
View at: Publisher Site | Google Scholar
J. Zeng, “A preliminary study to explore the effect of total flavonoids from Rhizoma drynariae on distraction osteogenesis based on the BMP-smad signal path,” Guangzhou University of Chinese Medicine, Guangzhou, China, 2016, Doctoral thesis.
View at: Google Scholar
L. Qin, Y. Xu, Y. Xu et al., “NCOA1 promotes angiogenesis in breast tumors by simultaneously enhancing both HIF1 alpha- and AP-1-mediated VEGFa transcription,” Oncotarget, vol. 6, no. 27, pp. 23890–23904, 2015.
View at: Publisher Site | Google Scholar
C. Galli, M. Piemontese, S. Lumetti, E. Manfredi, G. M. Macaluso, and G. Passeri, “GSK3b-inhibitor lithium chloride enhances activation of Wnt canonical signaling and osteoblast differentiation on hydrophilic titanium surfaces,” Clinical Oral Implants Research, vol. 24, no. 8, pp. 921–927, 2013.
View at: Publisher Site | Google Scholar
K. Fujita, M. Iwasaki, H. Ochi et al., “Vitamin E decreases bone mass by stimulating osteoclast fusion,” Nature Medicine, vol. 18, no. 4, pp. 589–594, 2012.
View at: Publisher Site | Google Scholar
T. Matsushita, Y. Y. Chan, A. Kawanami, G. Balmes, G. E. Landreth, and S. Murakami, “Extracellular signal-regulated kinase 1 (ERK1) and ERK2 play essential roles in osteoblast differentiation and in supporting osteoclastogenesis,” Molecular and Cellular Biology, vol. 29, no. 21, pp. 5843–5857, 2009.
View at: Publisher Site | Google Scholar
S. A. Bashanfer, M. Saleem, O. Heidenreich, E. Moses, and N. Yusoff, “Disruption of MAPK1 expression in the ERK signalling pathway and the RUNX1-RUNX1T1 fusion gene attenuate the differentiation and proliferation and induces the growth arrest in t (8;21) leukaemia cells,” Oncology Reports, vol. 41, no. 3, pp. 2027–2040, 2018.
View at: Publisher Site | Google Scholar
Q. Liu, Y. Zhou, and Z. Li, “PDGF-BB promotes the differentiation and proliferation of MC3T3-E1 cells through the Src/JAK2 signaling pathway,” Molecular Medicine Reports, vol. 18, no. 4, pp. 3719–3726, 2018.
View at: Publisher Site | Google Scholar
E. Avitan-Hersh, S. Tatur, M. Indelman et al., “Postzygotic HRAS mutation causing both keratinocytic epidermal nevus and thymoma and associated with bone dysplasia and hypophosphatemia due to elevated FGF23,” The Journal of Clinical Endocrinology & Metabolism, vol. 99, no. 1, pp. E132–Ed136, 2014.
View at: Publisher Site | Google Scholar
N. C. Harvey, A. Sheppard, K. M. Godfrey et al., “Childhood bone mineral content is associated with methylation status of the RXRA promoter at birth,” Journal of Bone and Mineral Research, vol. 29, no. 3, pp. 600–607, 2014.
View at: Publisher Site | Google Scholar
C. Coronnello, R. Hartmaier, A. Arora et al., “Novel modeling of combinatorial miRNA targeting identifies SNP with potential role in bone density,” PLoS Computational Biology, vol. 8, no. 12, Article ID e1002830, 2012.
View at: Publisher Site | Google Scholar
H. Zhi, “Research on human POL II promoter recognizing,” Tongji University, Shanghai, China, 2008, Master thesis.
View at: Google Scholar
N. Song, Z. Zhao, X. Ma et al., “Naringin promotes fracture healing through stimulation of angiogenesis by regulating the VEGF/VEGFR-2 signaling pathway in osteoporotic rats,” Chemico-Biological Interactions, vol. 261, pp. 11–17, 2017.
View at: Publisher Site | Google Scholar
Y. Moedinger, A. E. Rapp, A. Vikman et al., “Reduced terminal complement complex formation in mice manifests in low bone mass and impaired fracture healing,” The American Journal of Pathology, vol. 189, no. 1, pp. 147–161, 2019.
View at: Publisher Site | Google Scholar
Y. Xiao, J. Zeng, L. Jiao, and X. Xu, “Review for treatment effect and signaling pathway regulation of kidney-tonifying traditional Chinese medicine on osteoporosis,” China Journal of Chinese Materia Medica, vol. 43, no. 1, pp. 21–30, 2018.
View at: Google Scholar
L. Zhang, B. Xu, J. Li, and W. Wang, “Influence of mitogen-activated protein kinase inhibitor in the process of assemble flavone of rhizome drynaria promoting the osteogentic differentiation of myoblasts,” Chinese Journal of Tissue Engineering Research, vol. 20, no. 51, pp. 7622–7627, 2016.
View at: Google Scholar
X. Huang, S. Yuan, and C. Yang, “Effects of drynaria total flavonoid on osteogenic differentiation of rat dental pulp stem cells via PI3K/Akt pathway,” Chinese Journal of Tissue Engineering Research, vol. 17, no. 1, pp. 92–97, 2013.
View at: Google Scholar
F. Lin, Y. Zhu, and G. Hu, “Naringin promotes cellular chemokine synthesis and potentiates mesenchymal stromal cell migration via the Ras signaling pathway,” Experimental and Therapeutic Medicine, vol. 16, no. 4, pp. 3504–3510, 2018.
View at: Publisher Site | Google Scholar

Copyright

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

PDF Download Citation

Download other formats

Order printed copies