Application of Natural Products in Experimental, and Cheminformatics-based Therapy for Neurological DisordersView this Special Issue
Research Article | Open Access
Zhencheng Xiong, Feng Yang, Wenhao Li, Xiangsheng Tang, Haoni Ma, Ping Yi, "Deciphering Pharmacological Mechanism of Buyang Huanwu Decoction for Spinal Cord Injury by Network Pharmacology Approach", Evidence-Based Complementary and Alternative Medicine, vol. 2021, Article ID 9921534, 20 pages, 2021. https://doi.org/10.1155/2021/9921534
Deciphering Pharmacological Mechanism of Buyang Huanwu Decoction for Spinal Cord Injury by Network Pharmacology Approach
Objective. The purpose of this study was to investigate the mechanism of action of the Chinese herbal formula Buyang Huanwu Decoction (BYHWD), which is commonly used to treat nerve injuries, in the treatment of spinal cord injury (SCI) using a network pharmacology method. Methods. BYHWD-related targets were obtained by mining the TCMSP and BATMAN-TCM databases, and SCI-related targets were obtained by mining the DisGeNET, TTD, CTD, GeneCards, and MalaCards databases. The overlapping targets of the abovementioned targets may be potential therapeutic targets for BYHWD anti-SCI. Subsequently, we performed protein-protein interaction (PPI) analysis, screened the hub genes using Cytoscape software, performed Gene Ontology (GO) annotation and KEGG pathway enrichment analysis, and finally achieved molecular docking between the hub proteins and key active compounds. Results. The 189 potential therapeutic targets for BYHWD anti-SCI were overlapping targets of 744 BYHWD-related targets and 923 SCI-related targets. The top 10 genes obtained subsequently included AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN. Fifteen signaling pathways were also screened out after enrichment analysis and literature search. The results of molecular docking of key active compounds and hub target proteins showed a good binding affinity for both. Conclusion. This study shows that BYHWD anti-SCI is characterized by a multicomponent, multitarget, and multipathway synergy and provides new insights to explore the specific mechanisms of BYHWD against SCI.
Spinal cord injury (SCI) is a group of disorders in which the integrity and continuity of the spinal cord are disrupted by trauma, tumor, hypoxia, inflammation, and other factors . According to the World Health Organization (WHO) report, about 250,000 to 500,000 people worldwide suffer from SCI every year, and the population is mainly under 30 years old . The Beijing SCI Epidemiology Survey Group reported in 2002 that the incidence of SCI in Beijing was about 60/106, an increase of nearly 10 times compared to the late 1980s, and that the per capita hospitalization cost was 27819.3 CNY . The high disability rate and high consumption of SCI have imposed a heavy burden on the families of patients and society, while how to protect and restore injured neurons after SCI has been considered a great challenge in clinical and experimental research . There is still a lack of effective treatment for SCI, but it is known that traditional Chinese medicine (TCM) has a long history and rich practical experience, which may be explored for the treatment of SCI .
Modern Chinese medicine scholars classify SCI as “Ti Duo” and “Wei Zheng” based on the TCM classics “Ling Shu Jing-Han Re Bing” and “Su Wen-Wei Lun” . Based on the analysis of the etiology and pathogenesis of SCI by TCM, as well as the pathological changes of the damaged Governor Vessel and the symptoms corresponding to SCI, the main treatment of TCM is to unblock Governor Vessel and strengthen “Yang Qi,” supplemented by tonifying the kidney and filling the essence, strengthening the tendons and bones, activating blood circulation and removing blood stasis, and warming the meridians . Buyang Huanwu Decoction (BYHWD), recorded in “Yi Lin Gai Cuo” by Wang Qingren in the Qing Dynasty, has been used in China for hundreds of years to treat stroke-related disabilities and improve the recovery of neurological function . The whole formula of BYHWD consists of seven herbs: Chishao (Radix Paeoniae Rubra), Chuanxiong (Rhizoma Chuanxiong), Danggui (Radix Angelicae Sinensis), Dilong (Pheretima), Huangqi (Radix Astragali), Honghua (Flos Carthami), and Taoren (Semen Persicae) . Studies have shown that BYHWD may promote the repair of SCI by upregulating the expression of Notch1 gene , and BYHWD combined with bone marrow mesenchymal stem cells (BMSCs) transplantation can promote the recovery after SCI by saving axotomized red nucleus neurons in rats . Although some progress has been made in the study of BYHWD for SCI, the complexity of herbal prescriptions with multiple components, targets, and pathways makes the study of the specific mechanism of action of BYHWD for SCI difficult to clarify.
In the last decade, with the development of modern multiomics theories such as genomics, proteomics, and metabolomics, the introduction of a systems biology perspective, and the application of bioinformatics, the concept of network pharmacology has emerged. According to the network basis of drug action, network pharmacology methods, including protein-protein interaction (PPI) network, Gene Ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, are applied to infer the specific pathways of multimolecular drug action in the human body, which is consistent with the multicomponent, multitarget, and multipathway synergy of Chinese herbal medicine and helps to clarify the mechanism of action of BYHWD (Figure 1).
2. Materials and Methods
2.1. Acquisition of BYHWD-Related Targets and SCI-Related Targets
The compounds related to BYHWD were obtained by searching “Chishao or CHI SHAO or Radix Paeoniae Rubra,” “Chuanxiong or CHUAN XIONG or Rhizoma Chuanxiong,” “Danggui or DANG GUI or Radix Angelicae Sinensis,” “Dilong or DI LONG or Pheretima,” “Huangqi or HUANG QI or Radix Astragali,” “Honghua or HONG HUA or Flos Carthami,” and “Taoren or TAO REN or Semen Persicae” in Traditional Chinese Medicine Systems Pharmacology (TCMSP, version 2.3, https://tcmspw.com/tcmsp.php) database , and then the active compounds that met the criteria (oral bioavailability (OB) ≥ 30 and drug-likeness (DL) ≥ 0.18) were further screened by the absorption, distribution, metabolism, and excretion (ADME) principle [8, 9], and the corresponding target protein information of the active compounds was also obtained. Subsequently, by using the UniProt database (https://www.uniprot.org/) (species should be set to “Homo sapiens”), the target protein name is converted into the corresponding gene name and UniProt ID . The compounds and targets related to BYHWD were also obtained directly by searching the above keywords in the BATMAN-TCM platform (http://bionet.ncpsb.org/batman-tcm/) according to the set criteria (score cutoff >25, adjusted P value < 0.05) . Finally, the results of the two databases are summarized, integrated, and deduplicated to obtain BYHWD-related targets.
DisGeNet database (https://www.disgenet.org/, v7.0) (score ≥ 0.01) , GeneCards (https://www.genecards.org/) (relevance score ≥ 20) , MalaCards (https://www.malacards.org/) , Comparative Toxicogenomics Database (CTD) (http://ctdbase.org/, last update by June 2020) (inference score ≥ 20) , and Therapeutic Target Database (TTD) (http://db.idrblab.net/ttd/, last update by June 1, 2020)  were used to screen out potential targets related to SCI by searching for the keyword “Spinal cord injury or Spinal cord diseases.” The results of the five databases are summarized, integrated, and deduplicated to obtain SCI-related targets.
2.2. Acquisition of Potential Therapeutic Targets for BYHWD Anti-SCI
Potential therapeutic targets of BYHWD in the treatment of SCI should be derived from the overlapping targets of BYHWD-related targets and SCI-related targets obtained from the Venn online tool (http://www.bioinformatics.com.cn/). Subsequently, we used Cytoscape software (version 3.7.2) to construct a network diagram of “BYHWD-active compounds-target genes-SCI” .
2.3. Construction of PPI Network and Acquisition of Hub Genes
PPI is composed of proteins that interact with each other to participate in various aspects of life processes such as biological signaling, regulation of gene expression, energy and material metabolism, and cell cycle regulation . The potential therapeutic targets of BYHWD against SCI were analyzed by a PPI network using the STRING database (http://string-db.org/; version 11) with the species set to “Homo sapiens” and a confidence level greater than 0.4 . Subsequently, the PPI network was visualized by Cytoscape software (version 3.7.2).
To obtain the hub genes for BYHWD anti-SCI, we use the CytoHubba plug-in in Cytoscape software, which provides 12 topological analysis methods, including Degree, Maximal Clique Centrality (MCC), BottleNeck, Betweenness, ClusteringCoefficient, Closeness, Density of Maximum Neighborhood Component (DMNC), EcCentricity, Edge Percolated Component (EPC), Maximum Neighborhood Component (MNC), Radiality, and Stress [20, 21].
2.4. GO and KEGG Pathway Enrichment Analysis of Potential Therapeutic Genes for BYHWD Anti-SCI
R is commonly used for statistical analysis and plotting software, and its open-source and free features have been favored by many programming enthusiasts, resulting in the birth of many convenient and efficient packages, such as the ClusterProfiler package for GO annotation and KEGG pathway enrichment analysis [22–24]. In this study, we used ClusterProfiler package (adjusted P value < 0.05) to explore the biological processes (BP), cellular components (CC), molecular functions (MF), and signaling pathways associated with BYHWD anti-SCI.
2.5. Molecular Docking of Hub Genes and Key Active Compounds
Based on the previously mined BYHWD-related targets and related active compounds, the active compounds and herbs corresponding to the top 10 hub genes need to be sorted out. The interrelationship between “herbs in BYHWD-active compounds-top 10 hub genes” was visualized by a Sankey diagram (http://sankeymatic.com/). In the following, we will select key active compounds as small molecule ligands and corresponding large molecule receptors (proteins encoded by the top 10 hub genes) for molecular docking analysis. The raw files of key active compounds (MOL2 format) can be downloaded from PubChem database (https://pubchem.ncbi.nlm.nih.gov/), and the raw files of the top 10 core genes (PDB format) can be downloaded from the RCSB protein data (http://www.rcsb.org/), and after corresponding processing by AutoDock Tool , they are finally converted to PDBQT format for molecular docking in Pymol software (https://pymol.org/2/; version 2.4.1) . The molecular docking score reflects the binding affinity, with smaller values representing higher affinity for the binding of small molecule ligands and large molecule receptors .
3.1. Acquisition Results of BYHWD-Related Targets and SCI-Related Targets
A total of 775 compounds (BATMAN-TCM: 219; after deduplication: 189) in BYHWD were preliminarily retrieved from the TCMSP database and BATMAN-TCM platform (score cutoff >25; adjusted P value < 0.05), of which 119 (BATMAN-TCM: 8) were from Chishao, 189 (87) were from Chuanxiong, 125 (77) were from Danggui, 87 (23) were from Huangqi, 189 (22) were from Honghua, and 66 (2) were from Taoren. Relevant compounds and targets’ information for Dilong was not retrieved in either database. After setting the ADME criteria (OB ≥ 30% and DL ≥ 0.18), a total of 103 active compounds (after removing duplicates and nontarget compounds: 60) were screened from the TCMSP database, of which 29 (After removing nontarget compounds: 14) were from Chishao, 7 (6) were from Chuanxiong, 2 (2) were from Danggui, 20 (17) were from Huangqi, 22 (17) were from Honghua, and 23 (19) were from Taoren. The BATMAN-TCM platform provides the gene names directly, while the TCMSP database first provides the target protein names of the active compounds and therefore needs to be converted to gene names via the UniProt database. Finally, after integrating the results of the two databases, we obtained 244 BYHWD-related active compounds and 744 BYHWD-related targets.
After deduplication, a total of 923 SCI-related targets were obtained by mining the DisGeNET (number: 84), TTD (number: 5), CTD (number: 553), GeneCards (number: 432), and MalaCards (number: 19) databases.
3.2. Acquisition Results of Potential Therapeutic Targets for BYHWD Anti-SCI
The overlapping targets of BYHWD-related targets and SCI-related targets were considered as potential therapeutic targets for BYHWD anti-SCI. As shown in Figure 2, we screened the overlapping targets of both by constructing a Venn diagram, and a total of 189 targets were obtained. Through multiple databases mining and construction of the Venn diagram, we obtained BYHWD-related compounds and potential therapeutic targets for BYHWD anti-SCI and finally used Cytoscape software to construct a “BYHWD-active compounds-target genes-SCI” network (Figure 3). In Figure 3, the lines between two nodes represent the existence of mutual relationships, and the larger the node is, the more relationships exist.
3.3. Acquisition Results of Hub Genes for BYHWD Anti-SCI and PPI Network Construction
In the STRING database, we first set the species to “Homo sapiens” and then entered 189 potential therapeutic targets for BYHWD anti-SCI to obtain the PPI network (Figure 4(a)), which involved 189 nodes and 3717 edges. The obtained TSV files were then imported into Cytoscape software (version 3.7.2) for further analysis and visualization (Figure 4(b)). To obtain the hub genes for BYHWD anti-SCI, based on the above the PPI network, we used the CytoHubba plug-in of Cytoscape software, which currently contains 12 topological analysis methods. We targeted the top 10 hub genes for BYHWD anti-SCI, and the results of the 12 algorithms each contained 10 genes. We integrated and analyzed the results to obtain a total of 45 different genes, which were then sorted by the number of algorithms to which they were attributed. Table 1 shows the basic information on hub genes for BYHWD anti-SCI. The top 10 hub genes (AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN) sorted are completely consistent with some algorithms (Degree, Closeness, and MNC) (Figure 4(c)).
(1): Degree, (2): Maximal Clique Centrality (MCC), (3): Betweenness, (4): BottleNeck, (5): Closeness, (6): ClusteringCoefficient, (7): Density of Maximum Neighborhood Component (DMNC), (8): Edge Percolated Component (EPC), (9): EcCentricity, (10): Maximum Neighborhood Component (MNC), (11): Radiality, and (12): Stress
3.4. Results of GO and KEGG Pathway Enrichment Analysis
GO and KEGG pathway enrichment analysis are important tools to explore the biological processes and signal pathways involved in clinical drug therapy for diseases, as well as key genes. Through the ClusterProfiler package in R, we performed GO enrichment analysis on the 189 potential therapeutic targets for BYHWD anti-SCI and obtained 2842 GO items (after adjustment, P < 0.05), including 2567 BP items, 99 CC items, and 176 MF items. We combined the top 10 enrichment results for each of GO-BP, GO-CC, and GO-MF by a bar graph (Figure 5). The basic information of the top 10 GO enrichment items is shown in Table 2.
In the same way, through the ClusterProfiler package in R, we performed KEGG pathway enrichment analysis on the 189 potential therapeutic targets for BYHWD anti-SCI and obtained 177 KEGG pathways (after adjustment, P < 0.05). The top 20 KEGG enrichment pathways are listed below, including AGE-RAGE signaling pathway in diabetic complications (hsa04933), Fluid shear stress and atherosclerosis (hsa05418), Prostate cancer (hsa05215), IL-17 signaling pathway (hsa04657), Hepatitis B (hsa05161), TNF signaling pathway (hsa04668), Kaposi sarcoma-associated herpesvirus infection (hsa05167), Human cytomegalovirus infection (hsa05163), Hepatitis C (hsa05160), Human T-cell leukemia virus 1 infection (hsa05166), Small cell lung cancer (hsa05222), Pancreatic cancer (hsa05212), Measles (hsa05162), Epstein-Barr virus infection (hsa05169), Nonalcoholic fatty liver disease (hsa04932), Influenza A (hsa05164), Endocrine resistance (hsa01522), Bladder cancer (hsa05219), Platinum drug resistance (hsa01524), and Apoptosis (hsa04210). The visualization of the abovementioned KEGG enrichment pathway is realized in Figure 6. More KEGG pathways were enriched, and another fine screening of the pathways was required. Therefore, the SCI-related literature was searched in the PubMed database, and then the retrieved potentially relevant pathways were compared with the enriched 177 pathways, and finally, 15 relatively relevant pathways for SCI were obtained. Subsequently, we used Cytoscape software to visualize the network relationship between potential therapeutic targets for BYHWD anti-SCI and pathways (Figure 7). Table 3 shows the basic information of 15 pathways that may be relevant for SCI.
3.5. Results of Molecular Docking of Hub Genes and Key Active Compounds
Sankey diagram was used to establish one-to-one correspondence between the top 10 hub genes (AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN) for BYHWD anti-SCI and BYHWD-related compounds, providing the basis for the next molecular docking analysis. Among them, the gene that targets the most active compounds is TNF, and the key active compounds that target the most hub genes are quercetin and luteolin (Figure 8). As seen in the Sankey diagram, 6 of the 33 active compounds (quercetin, luteolin, beta-carotene, kaempferol, baicalein, and beta-sitosterol) targeted more hub genes, so we decided to perform the molecular docking between these 6 active compounds and top 10 hub genes. The results of the molecular docking score are presented as a heat map (Figure 9), with a range of −6.4 to −10.9 kcal·mol−1, representing the large molecule protein receptor (hub genes) that binds well to the small molecule ligand (key active compounds). As shown in Figure 10, the molecular docking process of each key active compound and its corresponding protein encoded by the hub gene with the best docking affinity was visualized using Pymol. Table 4 shows the basic information of active compounds targeting hub genes in BYHWD.
OB, oral bioavailability; DL, drug-likeness.
According to the report of the National SCI Statistics Center, there are about 17,000 new cases of SCI in the United States each year, and the annual incidence of SCI is about 54 cases per million people . Over the past century, advances in understanding the mechanisms of SCI have changed the clinical management strategies for SCI, including surgical procedures, supportive measures, and rehabilitative training that have improved neurological outcomes and reduced morbidity in SCI patients . However, there is still a lack of effective treatments for SCI and thus has been stimulating the desire of researchers to explore and try to find effective treatments, including drugs, surgical methods, and rehabilitation programs. As we all know, TCM has a long history of being used to treat various diseases in China and other Asian countries . Among them, BYHWD, Jisuikang, single herb, and Governor Vessel electroacupuncture have all been shown in relevant studies to be used in the treatment of SCI [4, 5, 28–30].
BYHWD, which comes from Qing Dynasty medical classics, consists of seven herbs: Chishao, Chuanxiong, Danggui, Dilong, Huangqi, Honghua, and Taoren . Based on the theory of TCM, Huangqi is the main ingredient in this formula, which has the function of tonifying Qi and nourishing blood, eliminating blood stasis without harming the proper function; Danggui is the secondary herb, which has the effect of invigorating blood and harmonizing blood; Chishao, Chuanxiong, Dilong, Honghua, and Taoren are the adjuvant herbs in this formula; all the five herbs have the function of activating blood circulation and removing blood stasis, promoting the flow of Qi and blood circulation; the combination of the five herbs not only activates blood circulation and removes blood stasis but also helps the Qi and blood that the main herbs benefit to reach the whole body [3, 31]. Studies have shown that BYHWD in combination with neural stem cells (NSCs) attenuates the process of demyelination or improves the recovery of myelin and exerts a synergistic effect on the neurological function recovery . It has been shown that BYHWD restores hindlimb motor function in rats with SCI, and the neuroprotective effect is associated with the regulation of apoptosis-related protein expression . BYHWD rescues axotomized neurons and promotes functional recovery after spinal cord injury in rats. A study published in 2008 showed that BYHWD protected the injured neurons and promoted functional recovery after SCI in rats . However, since the specific mechanism of BYHWD for the treatment of SCI is still unclear due to its complex composition and wide range of effects, we decided to use a network pharmacology approach to investigate in this study and to lay the foundation for the next in-depth study.
In the first step, BYHWD-related compounds and related targets were obtained. We obtained relevant compounds and targets by mining two herbal-related databases, TCMSP and BATMAN-TCM. The mining methods of these two databases are different: TCMSP database (ADME principle: OB ≥ 30 and DL ≥ 0.18) first obtains BYHWD-related compounds, then obtains related target proteins, and finally needs to be converted into standard gene names by UniProt database; BATMAN-TCM platform (score cutoff >25, adjusted P value <0.05) directly obtains BYHWD-related compounds and target genes. The screening process of TCMSP database is divided into five steps: (1) the relevant compounds of each component of BYHWD are directly obtained (Chishao: 119, Chuanxiong: 189, Danggui: 125, Huangqi: 87, Honghua: 189, Taoren: 66, and Dilong: 0); (2) the relevant compounds of each component are screened by ADME criteria (Chishao: 29, Chuanxiong: 7, Danggui: 2, Huangqi: 20, Honghua: 22, Taoren: 23, and Dilong: 0); (3) the above compounds are screened again by removing the null targets (Chishao: 14, Chuanxiong: 6, Danggui: 2, Huangqi: 17, Honghua: 17, Taoren: 19, and Dilong: 0); (4) the above compounds are finally screened by removing duplicates (total number of active compounds: 60); and (5) the target proteins corresponding to the obtained active compounds are converted into gene names (total number of targets: 227). And the screening process of BATMAN-TCM platform is divided into two steps: (1) the relevant compounds and targets of each component of BYHWD are directly obtained (Chishao: 8, Chuanxiong: 87, Danggui: 77, Huangqi: 23, Honghua: 22, Taoren: 2, and Dilong: 0), and (2) the relevant compounds of each component are screened after deduplication (total number of active compounds: 189; the total number of targets: 586). The results of the two databases were integrated and deduplicated to finally obtain 244 BYHWD-related active compounds and 744 BYHWD-related targets. In the second step, a total of 923 SCI-related targets were acquired by mining 5 disease databases, including DisGeNET (number: 84), TTD (number: 5), CTD (number: 553), GeneCards (number: 432), and MalaCards (number: 19). Next, the Venn diagram of BYHWD-related targets and SCI-related targets was plotted by the website online tool to obtain 189 potential therapeutic targets of BYHWD against SCI.
In total, Cytoscape software was used four times in this study: (1) the “BYHWD-active compounds-target genes-SCI” network was constructed, involving 410 nodes and 1571 interactions; (2) the PPI network obtained from the STRING database was reconstructed here, involving 189 nodes and 3717 interactions; (3) the top 10 hub genes (AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN) of BYHWD against SCI were obtained based on 12 topological analysis methods (Degree, MCC, Stress, ClusteringCoefficient, EcCentricity, BottleNeck, Closeness, Radiality, Betweenness, EPC, DMNC, and MNC) contained in the CytoHubba plug-in; (4) the “potential therapeutic targets-pathways” network was constructed, involving 347 nodes and 2650 interactions. Studies have shown that the protein encoded by TP53 and other factors regulate the regeneration, germination, and functional recovery of axons after central nervous system injury ; the rapid upregulation of CASP3 mRNA seen after SCI in rats may be associated with cell death in the spinal cord ; inflammatory mediators such as TNF-α and IL6 mediate the recruitment of inflammatory cells to the site of injury and by targeting these cytokines may be a potential strategy to reduce secondary injury in SCI [36, 37]; by targeting MAPK1, the overexpression of miRNA-433-5p protects motor dysfunction and inflammation after SCI .
The Sankey diagram was used once, mainly to show the correspondence between the top 10 hub genes and the corresponding active compounds contained in BYHWD (quercetin, luteolin, beta-carotene, kaempferol, baicalein, beta-sitosterol, carvacrol, O-cresol, anisic acid, 1-methyl-2-dodecyl-4-(1h)-quinolone, 4-ethylresorcinol, angelicin, astragaloside I, astragaloside II, caffeic acid dimethyl ether, cibarian, cordycepin, dihydropinosylvin, dodecenoic acid, ellagic acid, ethyl-P-methoxycinnamate, formononetin, linoleic acid, M-cresol, M-ethylphenol, O-ethylphenol, oleic acid, paeoniflorin, P-cresol, P-ethylphenol, sebiferic acid, sucrose, and thymol). Among them, quercetin targeted six hub genes, luteolin targeted six hub genes, beta-carotene targeted five hub genes, kaempferol targeted four hub genes, baicalein targeted three hub genes, beta-sitosterol targeted two hub genes, anisic acid targeted two hub genes, and the others targeted one hub gene. The experimental results showed that quercetin exerts a protective effect on the spinal cord by inhibiting the activation of the p38MAPK/iNOS signaling pathway in SCI rats and subsequently regulating secondary oxidative stress ; the combination of luteolin and palmitoylethanolamide reduced autophagy in SCI ; beta-carotene effectively reduced the course of secondary damage events after SCI by blocking NF-κB pathway activation ; baicalein may alleviate the harm caused by SCI by activating PI3K and inducing autophagy to reduce neuronal apoptosis .
Based on the results of the Sankey diagram, we performed molecular docking of the top 10 genes and the top 6 compounds in terms of the number of targeted genes (quercetin, luteolin, beta-carotene, kaempferol, baicalein, and beta-sitosterol). The results of molecular docking of 60 cohorts of hub target proteins and key active compounds showed a good binding affinity for both. We used R software twice, once to present the scores of the above molecular docking as a heat map, and another time to perform GO and KEGG pathway enrichment analysis using the ClusterProfiler package.
The top 10 results of the GO enrichment analysis of 189 potential therapeutic targets for BYHWD anti-SCI, based on the adjusted P value (from small to large), are as follows: response to metal ion (GO-BP:0010038), response to nutrient levels (GO-BP:0031667), response to lipopolysaccharide (GO-BP:0032496), response to molecule of bacterial origin (GO-BP:0002237), cellular response to chemical stress (GO-BP: 0062197), response to antibiotic (GO-BP:0046677), response to oxidative stress (GO-BP: 0006979), response to reactive oxygen species (GO-BP:0000302), response to nutrient (GO-BP:0007584), cellular response to oxidative stress (GO-BP: 0034599), membrane raft (GO-CC:0045121), membrane microdomain (GO-CC:0098857), membrane region (GO-CC:0098589), vesicle lumen (GO-CC: 0031983), protein kinase complex (GO-CC: 1902911), RNA polymerase II transcription regulator complex (GO-CC: 0090575), serine/threonine-protein kinase complex (GO-CC:1902554), secretory granule lumen (GO-CC: 0034774), cytoplasmic vesicle lumen (GO-CC: 0060205), transcription regulator complex (GO-CC: 0005667), RNA polymerase II-specific DNA-binding transcription factor binding (GO-MF:0061629), DNA-binding transcription factor binding (GO-MF:0140297), tetrapyrrole binding (GO-MF:0046906), heme binding (GO-MF:0020037), electron transfer activity (GO-MF: 0009055), nuclear receptor activity (GO-MF:0004879), ligand-activated transcription factor activity (GO-MF:0098531), ubiquitin-like protein ligase binding (GO-MF: 0044389), steroid hormone receptor activity (GO-MF:0003707), and scaffold protein binding (GO-MF:0097110). The first 20 results of the (adjusted P value based) KEGG pathway enrichment analysis of 189 potential therapeutic targets of BYHWD against SCI have been presented, and we screened the enriched pathways again to find possible relevant pathways for SCI. By searching the relevant literature, we screened a total of 15 pathways with strong relevance, as follows: AGE-RAGE signaling pathway in diabetic complications (hsa04933) , TNF signaling pathway (hsa04668), Apoptosis (hsa04210) , Pathways of neurodegeneration-multiple diseases (hsa05022), Toll-like receptor signaling pathway (hsa04620) , PI3K-Akt signaling pathway (hsa04151) , HIF-1 signaling pathway (hsa04066) , MAPK signaling pathway (hsa04010) , NF-kappa B signaling pathway (hsa04064) , EGFR tyrosine kinase inhibitor resistance (hsa01521) , JAK-STAT signaling pathway (hsa04630) , Wnt signaling pathway (hsa04310) , AMPK signaling pathway (hsa04152) , and mTOR signaling pathway (hsa04150) . The above pathways and hub genes are not in a disconnected relationship; rather, many hub genes are important components of some of the pathways. Taking PI3K-Akt signaling pathway as an example, it contains AKT1, MAPK1, IL6, TP53, and VEGFA. Based on the results of the above network pharmacological analysis, we hypothesize that the key active compounds of BYHWD may exert anti-SCI effects through the above pathways and hub genes.
This study still has some limitations. First, the corresponding experimental validation is lacking. Second, the compounds, targets, and pathways contained in these databases may not be exhaustive. Finally, compounds screened based on ADME principles may be missing other important compounds.
4.2. Future Perspectives
This study uncovered the multicompound, multitarget, and multipathway characteristics of BYHWD through a network pharmacology approach. Among them, the key compounds include quercetin, luteolin, beta-carotene, kaempferol, baicalein, and beta-sitosterol; the hub genes include AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN; and the important pathways are the 15 related pathways mentioned above. These results clearly give a general direction for SCI-related researchers, and there is no doubt that the anti-SCI mechanism of BYHWD is much more than the above-mentioned research directions, and we will find more accurate targets and pathways while exploring them one by one.
Using a network pharmacology approach, we explored the potential mechanism of BYHWD as a classical herbal formula for treating nerve injury in the treatment of SCI. The top 10 core genes among the potential therapeutic targets for BYHWD anti-SCI include AKT1, IL6, MAPK1, TNF, TP53, VEGFA, CASP3, ALB, MAPK8, and JUN. Key active compounds of BYHWD against SCI include quercetin, luteolin, beta-carotene, kaempferol, baicalein, and beta-sitosterol. In addition, combining the KEGG pathway enrichment results and reviewing SCI-related literature, we obtained 15 SCI potentially relevant pathways, mainly including PI3K-Akt signaling pathway, MAPK signaling pathway, NF-kappa B signaling pathway, TNF signaling pathway, Apoptosis, Toll-like receptor signaling pathway, HIF-1 signaling pathway, JAK-STAT signaling pathway, Wnt signaling pathway, AMPK signaling pathway, and mTOR signaling pathway. The above analysis results demonstrated the multitarget, multicomponent, and multipathway characteristics of BYHWD in treating SCI, which laid the foundation for us to carry out specific experimental validation in the next step, and provided new ideas for researchers dedicated to the treatment of SCI with herbal extracts.
The data used to support the findings of this study are included within the article.
Zhencheng Xiong and Feng Yang are joint first authors.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
This work was supported by the National Natural Science Foundation of China (no. 81873141).
- G. Courtine and M. V. Sofroniew, “Spinal cord repair: advances in biology and technology,” Nature Medicine, vol. 25, no. 6, p. 898, 2019.
- Y. Zheng, S. Qi, F. Wu et al., “Chinese herbal medicine in treatment of spinal cord injury: a systematic review and meta-analysis of randomized controlled trials,” The American Journal of Chinese Medicine, vol. 48, no. 7, p. 1593, 2020.
- “Current status of traditional Chinese medicine research on spinal cord injury,” World Journal of Integrated Traditional and Western Medicine, vol. 12, no. 3, p. 440, 2017.
- P. Yang, A. Chen, Y. Qin et al., “Buyang huanwu decoction combined with BMSCs transplantation promotes recovery after spinal cord injury by rescuing axotomized red nucleus neurons,” Journal of Ethnopharmacology, vol. 228, p. 123, 2019.
- M. Zhang, Y. Chai, T. Liu, N. Xu, and C. Yang, “Synergistic effects of Buyang Huanwu decoction and embryonic neural stem cell transplantation on the recovery of neurological function in a rat model of spinal cord injury,” Experimental and Therapeutic Medicine, vol. 9, no. 4, p. 1141, 2015.
- Z. P. Guo, M. N. Huang, A. Q. Liu, Y. J. Yuan, J. B. Zhao, and X. F. Mei, “Buyang Huanwu decoction up-regulates Notch1 gene expression in injured spinal cord,” Neural Regeneration Research, vol. 10, no. 8, p. 1321, 2015.
- 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. 13, 2014.
- K. Tsaioun, B. Blaauboer, and T. Hartung, “Evidence-based absorption, distribution, metabolism, excretion (ADME) and its interplay with alternative toxicity methods,” ALTEX, vol. 33, no. 4, p. 343, 2016.
- S. Sheng, Z. Yang, F. Xu, and Y. Huang, “Network pharmacology-based exploration of synergistic mechanism of guanxin II formula (II) for coronary heart disease,” Chinese Journal of Integrative Medicine, vol. 27, no. 2, pp. 106–114, 2021.
- T. UniProt Consortium, “UniProt: the universal protein knowledgebase,” Nucleic Acids Research, vol. 46, no. 5, p. 2699, 2018.
- Z. Liu, F. Guo, Y. Wang et al., “BATMAN-TCM: a bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine,” Scientific Reports, vol. 6, p. 21146, 2016.
- J. Piñero, J. M. Ramírez-Anguita, J. Saüch-Pitarch et al., “The DisGeNET knowledge platform for disease genomics: 2019 update,” Nucleic Acids Research, vol. 48, no. D1, p. D845, 2020.
- G. Stelzer, N. Rosen, I. Plaschkes et al., “The GeneCards suite: from gene data mining to disease genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, pp. 1–30, 2016.
- N. Rappaport, M. Twik, I. Plaschkes et al., “MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search,” Nucleic Acids Research, vol. 45, no. D1, p. D877, 2017.
- A. Davis, C. Grondin, R. Johnson et al., “Comparative Toxicogenomics database (CTD): update 2021,” Nucleic Acids Research, vol. 49, no. D1, 2021.
- Y. Wang, S. Zhang, F. Li et al., “Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics,” Nucleic Acids Research, vol. 48, no. D1, p. D1031, 2020.
- P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, p. 2498, 2003.
- K. Yang, L. Zeng, A. Ge et al., “Hedysarum multijugum MaximThe effect of .- compound on ischemic stroke: a research based on network and experimental pharmacology,” Oxidative Medicine and Cellular Longevity, vol. 2020, Article ID 6072380, 2020.
- 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, p. D362, 2017.
- C. Xu, R. Li, and J. Wu, “Effects of Yuanhu- Zhitong tablets on alcohol-induced conditioned place preference in mice,” Biomedicine & Pharmacotherapy = Biomedecine & Pharmacotherapie, vol. 133, p. 110962, 2020.
- C. Chin, S. Chen, H. Wu, C. Ho, M. Ko, and C. Lin, “cytoHubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. S11, 2014.
- T. G. O. Consortium, “Gene Ontology consortium: going forward,” Nucleic Acids Research, vol. 43, p. D1049, 2015.
- M. Kanehisa, M. Araki, S. Goto et al., “KEGG for linking genomes to life and the environment,” Nucleic Acids Research, vol. 36, p. D480, 2008.
- G. Yu, L.-G. Wang, Y. Han, and Q.-Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, p. 284, 2012.
- O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, p. 455, 2010.
- N. Zhu and J. Hou, “Exploring the mechanism of action Xianlingubao Prescription in the treatment of osteoporosis by network pharmacology,” Computational Biology and Chemistry, vol. 85, p. 107240, 2020.
- J. Chang, L. Liu, Y. Wang, Y. Sui, H. Li, and L. Feng, “Investigating the multitarget mechanism of traditional Chinese medicine prescription for cancer-related pain by using network pharmacology and molecular docking approach. Evidence-based complementary and alternative medicine,” eCAM, vol. 7617261, p. 2020, 2020.
- Y. Lu, J. Yang, X. Wang et al., “Research progress in use of traditional Chinese medicine for treatment of spinal cord injury,” Biomedicine & Pharmacotherapy, vol. 127, p. 110136, 2020.
- C. Wu, Y. Zhou, P. Tu et al., “Jisuikang promotes the repair of spinal cord injury in rats by regulating NgR/RhoA/ROCK signal pathway. Evidence-based complementary and alternative medicine,” eCAM, vol. 9542359, p. 2020, 2020.
- H. Xu, Y. Yang, Q. Deng et al., “Governor Vessel electro-acupuncture promotes the intrinsic growth ability of spinal neurons through activating calcitonin gene-related peptide/α-calcium/calmodulin-dependent protein kinase/neurotrophin-3 pathway after spinal cord injury,” Journal of Neurotrauma, vol. 38, no. 6, pp. 734–745, 2020.
- Y.-S. Lee, S.-C. Woo, S.-Y. Kim, and J.-Y. Park, “Understanding the multi-herbal composition of Buyang Huanwu Decoction: a review for better clinical use,” Journal of Ethnopharmacology, vol. 255, Article ID 112765, 2020.
- X. Dong, X. He, and W. Gao, “Neuroprotective effects of the Buyang Huanwu decoction on functional recovery in rats following spinal cord injury,” The Journal of Spinal Cord Medicine, vol. 39, no. 1, p. 85, 2016.
- A. Chen, H. Wang, J. Zhang et al., “BYHWD rescues axotomized neurons and promotes functional recovery after spinal cord injury in rats,” Journal of Ethnopharmacology, vol. 117, no. 3, p. 451, 2008.
- Y. Joshi, M. G. Sória, G. Quadrato et al., “The MDM4/MDM2-p53-IGF1 axis controls axonal regeneration, sprouting and functional recovery after CNS injury,” Brain, vol. 138, no. 7, p. 1843, 2015.
- B. A. Citron, P. M. Arnold, C. Sebastian et al., “Rapid upregulation of caspase-3 in rat spinal cord after injury: mRNA, protein, and cellular localization correlates with apoptotic cell death,” Experimental Neurology, vol. 166, no. 2, p. 213, 2000.
- E. Esposito and S. Cuzzocrea, “Anti-TNF therapy in the injured spinal cord,” Trends in Pharmacological Sciences, vol. 32, no. 2, p. 107, 2011.
- J. Kim, H. P. Joshi, S. H. Sheen et al., “Resolvin D3 promotes inflammatory resolution, neuroprotection, and functional recovery after spinal cord injury,” Molecular Neurobiology, vol. 58, no. 1, p. 424, 2021.
- C. L. Zhou, F. Li, X. W. Wu et al., “Overexpression of miRNA-433-5p protects acute spinal cord injury through activating MAPK1,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 6, p. 2829, 2020.
- Y. Song, J. Liu, F. Zhang, J. Zhang, T. Shi, and Z. Zeng, “Antioxidant effect of quercetin against acute spinal cord injury in rats and its correlation with the p38MAPK/iNOS signaling pathway,” Life Sciences, vol. 92, no. 24-26, p. 1215, 2013.
- R. Siracusa, I. Paterniti, G. Bruschetta et al., “The association of palmitoylethanolamide with Luteolin decreases autophagy in spinal cord injury,” Molecular Neurobiology, vol. 53, no. 6, p. 3783, 2016.
- L. Zhou, L. Ouyang, S. Lin et al., “Protective role of β-carotene against oxidative stress and neuroinflammation in a rat model of spinal cord injury,” International Immunopharmacology, vol. 61, p. 92, 2018.
- Y. Li, S. Lin, C. Xu, P. Zhang, and X. Mei, “Triggering of autophagy by Baicalein in response to apoptosis after spinal cord injury: possible involvement of the PI3K activation,” Biological and Pharmaceutical Bulletin, vol. 41, no. 4, p. 478, 2018.
- H. Wang, Z. Zhao, C. Liu et al., “Receptor for advanced glycation end-products (RAGE) blockade do damage to neuronal survival via disrupting Wnt/β-catenin signaling in spinal cord injury,” Neurochemical Research, vol. 43, no. 7, p. 1405, 2018.
- K. Gao, J. Niu, and X. Dang, “Neuroprotection of melatonin on spinal cord injury by activating autophagy and inhibiting apoptosis via SIRT1/AMPK signaling pathway,” Biotechnology Letters, vol. 42, no. 10, p. 2059, 2020.
- S. Xu, J. Wang, J. Jiang et al., “TLR4 promotes microglial pyroptosis via lncRNA-F630028O10Rik by activating PI3K/AKT pathway after spinal cord injury,” Cell Death & Disease, vol. 11, no. 8, p. 693, 2020.
- Z. Ding, C. Dai, L. Zhong et al., “Neuregulin-1 converts reactive astrocytes toward oligodendrocyte lineage cells via upregulating the PI3K-AKT-mTOR pathway to repair spinal cord injury,” Biomedicine & Pharmacotherapy, vol. 134, p. 111168, 2021.
- H. Chen, J. Li, S. Liang et al., “Effect of hypoxia-inducible factor-1/vascular endothelial growth factor signaling pathway on spinal cord injury in rats,” Experimental and Therapeutic Medicine, vol. 13, no. 3, p. 861, 2017.
- J. Ye, R. Xue, Z. Y. Ji et al., “Effect of NT-3 on repair of spinal cord injury through the MAPK signaling pathway,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 5, p. 2165, 2020.
- C. Wang, T. Xu, B. Lachance et al., “Critical roles of sphingosine kinase 1 in the regulation of neuroinflammation and neuronal injury after spinal cord injury,” Journal of Neuroinflammation, vol. 18, no. 1, p. 50, 2021.
- W. Xue, Y. Zhao, Z. Xiao et al., “Epidermal growth factor receptor-extracellular-regulated kinase blockade upregulates TRIM32 signaling cascade and promotes neurogenesis after spinal cord injury,” Stem Cells, vol. 38, no. 1, p. 118, 2020.
- K. Alexander, H. Tseng, W. Fleming et al., “Inhibition of JAK1/2 tyrosine kinases reduces neurogenic heterotopic ossification after spinal cord injury,” Frontiers in Immunology, vol. 10, p. 377, 2019.
Copyright © 2021 Zhencheng Xiong 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.