Systems Pharmacology Uncovers the Multiple Mechanisms of Xijiao Dihuang Decoction for the Treatment of Viral Hemorrhagic Fever

Liu, Jianling; Pei, Tianli; Mu, Jiexin; Zheng, Chunli; Chen, Xuetong; Huang, Chao; Fu, Yingxue; Liang, Zongsuo; Wang, Yonghua

doi:https://doi.org/10.1155/2016/9025036

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2016 | Article ID 9025036 | https://doi.org/10.1155/2016/9025036

Systems Pharmacology Uncovers the Multiple Mechanisms of Xijiao Dihuang Decoction for the Treatment of Viral Hemorrhagic Fever

Jianling Liu,^1,2Tianli Pei,^1,2Jiexin Mu,^1,2Chunli Zheng,²Xuetong Chen,²Chao Huang,²Yingxue Fu,²Zongsuo Liang,³and Yonghua Wang²

Academic Editor: Hyunsu Bae

Received21 Dec 2015

Revised17 Mar 2016

Accepted23 Mar 2016

Published27 Apr 2016

Abstract

Background. Viral hemorrhagic fevers (VHF) are a group of systemic diseases characterized by fever and bleeding, which have posed a formidable potential threat to public health with high morbidity and mortality. Traditional Chinese Medicine (TCM) formulas have been acknowledged with striking effects in treatment of hemorrhagic fever syndromes in China’s history. Nevertheless, their accurate mechanisms of action are still confusing. Objective. To systematically dissect the mechanisms of action of Chinese medicinal formula Xijiao Dihuang (XJDH) decoction as an effective treatment for VHF. Methods. In this study, a systems pharmacology method integrating absorption, distribution, metabolism, and excretion (ADME) screening, drug targeting, network, and pathway analysis was developed. Results. 23 active compounds of XJDH were obtained and 118 VHF-related targets were identified to have interactions with them. Moreover, systematic analysis of drug-target network and the integrated VHF pathway indicate that XJDH probably acts through multiple mechanisms to benefit VHF patients, which can be classified as boosting immune system, restraining inflammatory responses, repairing the vascular system, and blocking virus spread. Conclusions. The integrated systems pharmacology method provides precise probe to illuminate the molecular mechanisms of XJDH for VHF, which will also facilitate the application of traditional medicine in modern medicine.

1. Introduction

Viral hemorrhagic fevers (VHF) are a group of systemic diseases caused by certain viruses, such as Ebola, Lassa, Dengue, and Crimean-Congo hemorrhagic fever viruses. Patients with VHF show the common cardinal symptoms, including fever, hemorrhages, and shock [1]. Data obtained over the past years indicate that these diseases are characterized by intense inflammatory responses with generalized signs of increased vascular permeability, severely impaired immune functions, diffuse vascular dysregulation, and coagulation abnormalities [2, 3]. VHF are generally prevalent in developing countries, which have posed a serious public health threat with high mortality, morbidity, and infectivity in recent years [4]. Currently, many large pharmaceutical companies are pursuing an effective antiviral therapy for VHF. Although the broad-spectra antiviral drug ribavirin is approved for treatment of several types of VHF, there remains a need for a safe and more effective medication to replace the antiviral drug [5].

Traditional Chinese Medicine (TCM) formulas consisting of complex mixtures of multiple plants play an outstanding role in the treatment of various acute infectious diseases because of the pharmacological and pharmacokinetic synergistic effects of the abundant bioactive ingredients [6]. A series of TCM prescriptions for hemorrhagic fever syndromes have been described in history [7, 8]. For example, XJDH is a famous TCM formula for treating hemorrhagic fever syndromes [9]. XJDH originally comes from “Prescriptions Worth A Thousand Gold” which is written by the “Medicinal King” Sun Simiao in the Tang Dynasty (around 700 AD) [10, 11]. The components of the formula include Rhino horn (substituted by Buffalo Horn now, Shui Niujiao in Chinese), Rehmannia dried rhizome (Sheng Dihuang in Chinese), Paeonia lactiflora Pall. (Shao Yao in Chinese), and Paeonia suffruticosa Andr. (Mu Danpi in Chinese). Actually, XJDH has been normally used for cooling the blood for hemostasis, stopping bleeding accompanied with fever, removing toxic substances, and treating the cases of high fever and sweating, spontaneous bleeding, hemoptysis, and nosebleeds [12, 13]. Although the therapeutic efficiency of XJDH in the treatment of VHF is attractive, several fundamental questions are still unclear. What are the potential active ingredients of XJDH? What are the underlying molecular mechanisms of action of the formula in the treatment of VHF? What are the precise targets of these medicines? Since the multiple components-multiple targets interaction model of TCM formulas, traditional experimental research methods show up the shortcomings of long-term investment.

Fortunately, as an emerging discipline, systems pharmacology provides a new way to solve the complex pharmacological problems [14]. Systems pharmacology integrates pharmacokinetic data (ADME/T characteristics of a drug) screening together with targets prediction, networks, and pathways analyses to explore the drug actions from molecular and cellular levels to tissue and organism levels. It also provides an analysis platform for decoding molecular mechanisms of TCM formulas. In our previous work, a series of systems pharmacology methods have been exploited to uncover the underlying mechanisms of action of TCM formulas for cardiovascular diseases, depression, and cancer [15–17].

The purpose of the present study is to investigate the underlying molecular mechanisms of XJDH in treating VHF based on systems pharmacology method. Firstly, four pharmacokinetic models, including oral bioavailability (OB), drug-likeness (DL), Caco-2 permeability, and drug half-life (HL), were employed to filter out the potential active ingredients with favorable ADME profiles from XJDH. Then, based on an integrated target prediction method which combined the biological and mathematical models, the corresponding targets of these active ingredients were identified. Finally, the network pharmacology and VHF-related signaling pathways analysis was carried out to systematically disclose the underlying interactions between drugs, target proteins, and pathways. The detailed flowchart of the systems pharmacology method is shown in Figure 1.

2. Materials and Methods

2.1. Active Compounds Database

All chemicals of these four medicines in XJDH were manually collected from a wide-scale text mining and our in-house developed database: the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP, http://lsp.nwsuaf.edu.cn/tcmsp.php) [18]. In order to obtain the potential active compounds from these medicines, we applied a method incorporating OB, DL, Caco-2 permeability, and drug HL evaluation in this work.

2.1.1. OB Prediction

OB is defined as “the ratio of how many active components absorbed into the circulatory system to play a role at the site of action.” OB is one of the vital pharmacokinetic profiles in active compounds screening processes. In this work, the OB screening was calculated by a robust in-house system, OBioavail1.1 [19], and components with OB ≥ 30% were selected as the candidate molecules for further study. The following two basic sections describe the design principles of the threshold: () information from the studied medicines is obtained as much as possible using the least number of compounds and () the established model can be elucidated within reason by the reported pharmacological data [6].

2.1.2. DL Prediction

DL generally means “molecule which holds functional groups and/or has physical properties consistent with the majority of known drugs” [20]. In this study, we performed a self-constructed model pre-DL (predicts drug-likeness) based on the molecular descriptors and Tanimoto coefficient [21]. The DL index of the compounds was calculated by Tanimoto coefficient defined aswhere is the molecular descriptors of herbal compounds, represents the average molecular properties of all compounds in DrugBank database (http://www.drugbank.ca/) [22]. The DL ≥ 0.18 (average value for DrugBank) was defined as the criterion to select those drug-like compounds which are chemically suitable for drugs.

2.1.3. Caco-2 Permeability Prediction

The majority of orally administered drugs absorption occurs in the small intestine where the surface absorptivity greatly improves with the presence of villi and microvilli [23]. Previously, researchers have developed a quantity of in silico drug absorption models using in vitro Caco-2 permeability in drug discovery and development processes [24]. In this study, based on 100 drug molecules with satisfactory statistical results, a robust in silico Caco-2 permeability prediction model pre-Caco-2 (predicts Caco-2 permeability) was employed to predict the compound’s intestinal absorption [25]. Finally, on the account of the fact that compounds with Caco-2 value less than 0 are not permeable, in the study, the threshold of Caco-2 permeability was set to 0.

2.1.4. Drug HL Prediction

An in silico pre-HL (predicts half-life) has been developed in our previous work to calculate the drug HL by using the C-partial least square (C-PLS) algorithm which is supported by 169 drugs with known half-life values [26–28]. HL evaluates the time needed for compounds in the body to fall by half, and components with long HL were selected as the candidate molecules.

In order to obtain the potential active ingredients, the screening principle was defined as follows: OB ≥ 30%; Caco-2 ≥ 0; DL ≥ 0.18; or long HL.

2.2. Drug Targeting

Apart from screening out the active compounds, the therapeutic targets exploration is also a vital stage. Firstly, the potential targets exploration was fulfilled based on the systematic drug targeting tool (SysDT) as described in our previous work. Based on two mathematical tools, Random Forest (RF) and Support Vector Machine (SVM), the method can comprehensively ascertain the compound-target interaction profiles [29]. These two models exert great property of predicting the drug-target mutual effects with a concordance of 82.83%, a sensitivity of 81.33%, and a specificity of 93.62%. In this work, the compound-target interactions with SVM score ≥ 0.8 and the RF score ≥ 0.7 were selected for further research. Secondly, a recently developed computational model named weighted ensemble similarity (WES) was also introduced to detect drug direct targets [30]. For internal validation, this model performed remarkably well in predicting the binding (average sensitivity 72%, SEN) and the nonbinding (average specificity 82%, SPE) patterns, with the average areas under the receiver operating curves (ROC, AUC) of 85.2% and an average concordance of 77.5%. Thirdly, the obtained protein targets were mapped to the database UniProt (http://www.uniprot.org/) for normalization [31]. Finally, in order to identify and analyze the specific biological properties of the potential targets, the Gene Ontology (GO) biological processes were introduced to dissect target genes in a hierarchically structured way based on biological terms [32]. The GlueGO, a Cytoscape plug-in, was utilized to interpret the biology processes of large lists of genes.

2.3. Network Construction and Analysis

In order to explore the multiple mechanisms of action of XJDH for VHF, currently we analyzed the relationship between candidate compounds and potential targets by constructing the drug-target network (D-T network), in which all active compounds are connected to their targets. The network was generated by Cytoscape 2.8.1 [33]. In the network, compounds and targets are represented by nodes, while the interactions between them are represented by edges. In addition, a vital topological parameter, namely, degree was analyzed by the plugin NetworkAnalyzer of Cytoscape [34]. The degree of a node is defined as the number of edges connected to the node.

2.4. Pathway Construction and Analysis

At the pathway level, in order to probe into the action mechanisms of the formula for VHF, an incorporated “VHF pathway” was established based on the current knowledge of VHF pathology. Firstly, the obtained human target proteins were collected to be input into the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/) database to acquire the information of pathways. Then, based on the obtained information of basic pathways, we assembled an incorporated “VHF pathway” by picking out closely linked pathways related with VHF pathology.

3. Results

3.1. Active Compounds Screening

We employed four ADME parameters to screen out the potential active components of XJDH. As a result, from the 136 compounds of XJDH (as shown in supporting information, Table S1, in Supplementary Material available online at http://dx.doi.org/10.1155/2016/9025036), a total of 20 active compounds pass through the criteria of OB ≥ 30%, Caco-2 ≥ 0, DL ≥ 0.18, or long HL. Besides, in order to obtain a more accurate result, some certain rejected compounds, which have relatively poor pharmacokinetic properties, but are the most abundant and active ingredients of certain herbs, were also selected as the active components for further research. For example, although catalpol (MOL108) has poor OB, Caco-2, and HL properties, it has been reported to be rich in the roots of Rehmannia dried rhizome. And rehmaglutin D (MOL116) and paeoniflorin (MOL046) with poor OB (14.43%) take a large proportion in Paeonia lactiflora Pall. Thus, the three compounds were also retained for further analysis. Finally, a total of 23 active ingredients were obtained in this study (as shown in Table 1).

For all these 23 ingredients, many of them have been reported to demonstrate significant biologic activity including anti-inflammatory, antivirus, antipyretic, and immune-regulatory activities and protection effect of vascular endothelial cell. For instance, methyl gallate (MOL006, OB = 30.91%, Caco-2 = 0.26, and long HL), obtained from Paeonia lactiflora Pall., shows antivirus activity by interacting with virus proteins and altering the adsorption and penetration of the virion [35]. Salicylic acid (MOL045, OB = 32.13%, Caco-2 = 0.63, and long HL) and paeoniflorin (MOL046) with poor OB from Paeonia lactiflora Pall. exhibit antipyretic, anti-inflammatory, and immune-regulatory activities [36, 37]. In addition, kaempferol (MOL060, OB = 69.61%, Caco-2 = 0.15, DL = 0.24, and long HL), paeonol (MOL072, OB = 30.98%, Caco-2 = 0.91, long HL), and eugenol (MOL070, OB = 44.47%, Caco-2 = 1.36, and long HL), the main active compounds of the radix of Paeonia suffruticosa Andr., have been reported to have potential therapeutic effect for inflammation and vascular injury disorders [38–40]. Besides, it is worth noting that β-sitosterol (MOL018, OB = 36.9%, DL = 0.75) is a common ingredient of Rehmannia dried rhizome, Paeonia lactiflora Pall., and Paeonia suffruticosa Andr., indicating that these active compounds may show synergetic pharmacological effects on VHF.

3.2. Drug Targeting and Functional Analysis

Traditional information retrieval approaches of therapeutic targets of drugs are expatiatory and complicated [41]. To overcome this barrier, we introduced our previous developed target prediction model [29, 30] to dissect interactions between drugs and proteins. As a result, 23 candidate compounds are linked with 118 candidate targets (as shown in Table 2). The results show that many components simultaneously can act on more than one target and many targets can connect to all of the four medicines, demonstrating the promiscuous actions and analogous pharmacological effects of the bioactive molecules. For instance, kaempferol (MOL060) not only serves as the restrainer of Prostaglandin G/H synthase 2 [42] but also acts as the inhibitor of tumor necrosis factor [43]. And β-sitosterol (MOL018), which is shared by Rehmannia dried rhizome, Paeonia lactiflora Pall., and Paeonia suffruticosa Andr., acts as the activator of estrogen receptor [44] and transcription factor AP-1 [45]. Meanwhile, the results show that different drugs in XJDH can immediately impact on the common targets such as DNA ligase 1 (LIG1), indicating the synergism or cumulative effects of the drug molecules.

In general, vascular system, particularly the endothelium, plays a key role in VHF development [3]. A strong inflammatory response characterized by high circulating concentrations of cytokines and chemokines occurs early during the VHF infectious process [46]. And the patients’ immune functions might also be severely impaired; innate defenses are further hindered by the loss of natural killer cells [47]. The relevant biological processes of above targets were revealed by GlueGO (as shown in supporting information, Table S2). Figure 2 provides primary biological processes of these targets by cluster analysis. It is interesting to note that these targets are involved in a variety of biological processes including regulation of macrophage derived from cell differentiation, regulation of vasoconstriction and vasodilation, nitric oxide biosynthetic process, and vascular process in circulatory system. These biological processes largely fall into three groups: controlling inflammation response, modulating the immune system, and accommodation of vascular system. For example, peroxisome proliferator-activated receptor gamma (PPARγ), tumor necrosis factor receptor superfamily member 1A (TNFRSF1A), beta-2 adrenergic receptor (ADRB2), and so forth are involved in the regulation of acute inflammatory response. Nitric oxide synthase, inducible (NOS2), dipeptidyl peptidase 4 (DPP4), and so forth are associated with regulation of immune effector process, while nitric oxide synthase, endothelial (NOS3), NOS2, Krüppel-like factor 5 (KLF5), and so forth are directly connected to blood vessel remodeling and blood vessel morphogenesis. These suggest that XJDH might exert the therapeutic effect on VHF mainly through anti-inflammation, enhancing immunity and vascular repair therapy.

3.3. Drug-Target Network Construction and Analysis

As shown in Figure 3, D-T network is constructed including 141 nodes (23 active compounds and 118 potential targets) and 382 edges. The degrees of the candidate compounds are shown in Table 1; this provides us with an intuitionistic concept to distinguish those highly connected vital compounds or targets from the others in the network. The results of network analysis show that 18 out of 23 candidate compounds are linked with more than ten targets, among which kaempferol (MOL060) displays the highest number of target interactions (degree = 41), followed by eugenol (MOL070, degree = 34) and paeonol (MOL072, degree = 27). This confirms the multitarget properties of herbal compounds. We speculate that the top three ingredients might be the crucial elements in the treatment of VHF. For instance, kaempferol (MOL060) is predicted to interact with 41 targets like calcitonin gene-related peptide type 1 receptor (MAPK14), PPARγ, and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit gamma isoform (PIK3CG). MAPK14 takes part in the vascular endothelial growth factor (VEGF) synthesis through the mediation of angiotensin II [48]. VEGF can induce angiogenesis and improve the increased vascular permeability, so as to prevent bleeding in patients with VHF [49]. Besides, kaempferol is also found to significantly upregulate the transcriptional activity of PPARγ, which acts as an inhibitor of inflammatory gene expression and vandalizes proinflammatory transcription factor signaling pathways in vascular cells [50]. Additionally, previous finding suggests that PIK3CG interacting with kaempferol can participate in inflammation processes and influence the innate immune system [51]. Thus, these key active ingredients of XJDH work mainly by modulating inflammatory factor, innate immune system, and VEGF.

Meanwhile, the results also show that one target can be hit by multiple compounds from different medicines, indicating synergism or summation effects of the formula. According to the D-T network analysis, 64 out of the 118 targets have at least two links with the components of different herbs. XJDH exerts its therapeutic effect for VHF by binding and regulating particular protein targets. For instance, prostaglandin G/H synthase 2 (PTGS2) is simultaneously targeted by 10 active compounds including eugenol (MOL070, from Paeonia suffruticosa Andr.), paeonol (MOL072, from Paeonia suffruticosa Andr.), and salicylic acid (MOL045, from Paeonia lactiflora Pall.). PTGS2 is possibly an effective marker of platelet dysfunction; a reduced PTGS2 expression in the VHF primate model cells could directly result in platelet dysfunction [52]. Fortunately, in agreement with our study, previous findings suggest that eugenol and paeonol can control the expression of PTGS2 through the suppression of NF-κB in macrophage [53–55], so as to recover the function of thrombocyte. Study shows that severe disseminated intravascular coagulation is the mechanisms of bleeding in all VHF [3]. Prothrombin (F2), the precursor substances of clotting enzyme, plays a crucial role in optimizing the procoagulant activity through controlling the anticoagulant function of meizothrombin [56]. Thus, these 10 ingredients such as methyl gallate (MOL006, from Paeonia lactiflora Pall.), salicylic acid (MOL045, from Paeonia lactiflora Pall.), and apocynin (MOL053, Paeonia suffruticosa Andr.) interacting with F2 may be the key factors in the treatment of bleeding in patients with VHF. By analyzing the above D-T network, we can conclude that XJDH produces the healing efficacy for VHF probably by three different ways, intervening in the process of inflammation, boosting immune reaction, and repairing vascular system.

3.4. Pathway Analysis

To explore the integral regulation of XJDH for the treatment of VHF, we assembled an integrated “VHF pathway” (Figure 4) on the basis of the current knowledge of VHF pathogenesis. By means of inputting the obtained human target proteins into KEGG pathway database, result shows that 110 of the 118 targets can be mapped to the KEGG pathways, including NF-κB signaling pathway, AMPK pathway, and PI3K-AKT signaling pathway. Now, three detailed therapeutic modules are provided (inflammation module, angiogenesis module, and virus spreading module).

3.4.1. Inflammation Module

As shown in Figure 4, 7 key proteins targeted by XJDH are mapped onto a key inflammation process, namely, NF-κB signaling pathway, indicating the anti-inflammatory action may play a vital role in the treatment of VHF. Patients with VHF have a strong inflammatory response with high inflammatory cytokines and chemokines levels such as IL-1β, TNF, and IL-6 in the early phase of VHF [46]. The expression of inflammatory cytokines (TNF-α, IL-1β, IL-6, and IL-8) is mediated by NF-κB [57], and NF-κB is one of the most important regulators of proinflammatory gene expression. The result demonstrates that paeoniflorin (MOL046) from Paeonia lactiflora Pall. and Paeonia suffruticosa Andr. can regulate transcription factor NF-κB activity. Meanwhile, other researchers have verified that paeoniflorin can restrain the activation of the NF-κB pathway via inhibiting IκB kinase [58]. In addition, vascular adhesion molecule 1 (VCAM-1), a cell adhesion molecule, also plays an important role in the pathogenesis of inflammatory and immune processes [59]. Our work indicates that kaempferol (MOL060) and β-sitosterol (MOL018) can regulate the expression of inflammatory cytokines by targeting vascular adhesion molecule 1 (VCAM-1). Consequently, the foregoing analysis shows that XJDH has the effect of ameliorating the symptoms of inflammation disorders of patients with VHF.

3.4.2. Virus Spreading Module

In Figure 4, the phosphoinositide-3 kinase (PI3K) pathway is a significant cell signaling pathway that regulates diverse cellular activities including cell proliferation, differentiation, apoptosis, and vesicular trafficking. Notably, our research shows that, in line with previous studies, kaempferol (MOL060) is predicted to modulate PI3Ks activity, and AKT is also a target for kaempferol (MOL060) and paeonol (MOL072) [60–62]. Besides, PI3Ks, a family of lipid kinases, can prevent hemorrhagic fever virus entry into host cells by regulating cellular activities of vesicular trafficking [63]. AKT is a major downstream effector of the PI3K pathway, and this target protein can control the expression of many molecules directly or indirectly. Moreover, evidence also suggests that activity of PI3K/Akt pathway is required for hemorrhagic fever virus intruding into the host cells [63]. Therefore, depressors of PI3K and AKT dramatically reduced the risk of hemorrhagic fever virus infection at an early step during the replication cycle. These above analyses show that XJDH could make effective control of hemorrhagic fever virus entry into cells, thus blocking virus spread by interfering with the PI3K-AKT signaling pathway.

3.4.3. Angiogenesis Module

VHF is a severe multisystem syndrome characterized by diffuse vascular damage. The vascular system, particularly the vascular endothelium, seems to be directly and indirectly targeted by hemorrhagic fever viruses [3]. In the VHF pathway shown in Figure 4, PI3K pathway and AMPK pathway are involved in regulating the angiogenesis progress. We find out that apocynin (MOL053), kaempferol (MOL060), methyl salicylate (MOL066), and eugenol (MOL070) can affect the activity of endothelial NO synthase (eNOS) and then bring about NO production changes in endothelial cell. Moreover, at present, a large number of researches indicate that the synthesis of bioactive endothelium-derived NO is required for the progress of angiogenesis [64–66]. Therefore, the evidence presented enables us to reasonably conclude that the XJDH takes part in regulation of angiogenesis progress through PI3K signaling pathway and AMPK pathway.

4. Discussion

Actually, XJDH has been normally used for cooling the blood for hemostasis, stopping bleeding accompanied with fever, removing toxic substances, and treating the cases of high fever and sweating, spontaneous bleeding, hemoptysis, and nosebleeds [13, 67]. Although XJDH has been used historically for treating hemorrhagic fever syndromes, the specific bioactive molecules responsible for VHF and their precise mechanisms of action are still unclear. Thus, in this work, a systems pharmacology method combining the screening active components, drug targeting, network, and pathway analysis was carried out, so as to uncover the active ingredients, targets, and pathways of XJDH and systematically decipher its therapeutic mechanism of actions.

Our results show that 23 active ingredients were obtained from XJDH, and 118 potential targets were predicted. These manifest that the characteristics of XJDH are multicomponent botanical therapeutics and multitargets synergetic therapeutic effects. The GO analysis of targets and integrated D-T network analysis demonstrate the synergistic effects of XJDH for the treatment of VHF mainly through boosting of immune system, inhibiting inflammatory response, and repairing vascular system. Meanwhile, the integrated “VHF pathway” analysis in our work shows that XJDH might simultaneously regulate multitargets/pathways coupled with a range of therapeutic modules, for example, anti-inflammation, antivirus, and angiogenesis.

Now most researchers believe that VHF can be attributed to the simultaneous occurrence of multiple pathogenic mechanisms. They are mainly as follows: hemorrhagic fever virus infection stimulates macrophages to release cytokines, chemokines, and other mediators, causing fever, malaise, alterations in vascular function, and a shift in the coagulation system toward a procoagulant state, and immune functions might also be severely impaired [2]. Besides, hemorrhagic fever virus can target the vascular system directly and indirectly and cause endothelial activation and dysfunction [3]. In this study, we show here for the first time using GO enrichment analysis, network analysis, and integrated pathway analysis that XJDH significantly enriches target genes involved in reducing the inflammation response, enhancing immunity, combating the spreading virus, and preventing vascular dysfunction. And more experiments are needed to verify the validity of the results in further research works.

5. Conclusions

The result of this study provides bioactive ingredients, vital targets, and pathways of XJDH. We have come to the conclusion that the action mechanisms of XJDH for VHF mainly include restoring the immune system and enhancing immune response, ameliorating the symptoms of inflammation disorders, improving their vascular endothelial dysfunction, and combating the spreading virus. The systems pharmacology method established in our work provides preliminary clues that the multilayer networks of drug-target paradigm may be valuable for the modernization of TCM formulas at molecular level and then push forward their acceptance into mainstream medicine.

Competing Interests

The authors declare that there are no competing interests with respect to this study.

Authors’ Contributions

Jianling Liu and Tianli Pei contributed equally to this work.

Acknowledgments

This work was supported by Grants from Northwest A&F University, National Natural Science Foundation of China (31170796 and 81373892), and Science and Technology Department of Shaanxi Province, Natural Science Fund Projects (2014JM3063). It was also supported in part by China Academy of Chinese Medical Sciences (ZZ0608).

Supplementary Materials

Table S1: the 136 compounds information of XJDH.

Table S2: the relevant biological processes of the candidate targets.

Supplementary Tables

References

P. B. Jahrling, A. M. Marty, and T. W. Geisbert, “Viral hemorrhagic fevers,” in Medical Aspects of Biological Warfare, pp. 271–310, Office of the Surgeon General, United States Army, and Borden Institute, Walter Reed Army Medical Center, Washington, DC, USA, 2007.
View at: Google Scholar
M. Bray, “Pathogenesis of viral hemorrhagic fever,” Current Opinion in Immunology, vol. 17, no. 4, pp. 399–403, 2005.
View at: Publisher Site | Google Scholar
H.-J. Schnittler and H. Feldmann, “Viral hemorrhagic fever-a vascular disease?” Thrombosis and Haemostasis, vol. 89, no. 6, pp. 967–972, 2003.
View at: Google Scholar
L. Borio, T. Inglesby, C. J. Peters et al., “Hemorrhagic fever viruses as biological weapons: medical and public health management,” JAMA, vol. 287, no. 18, pp. 2391–2405, 2002.
View at: Publisher Site | Google Scholar
M. Bray, “Highly pathogenic RNA viral infections: challenges for antiviral research,” Antiviral Research, vol. 78, no. 1, pp. 1–8, 2008.
View at: Publisher Site | Google Scholar
X. Wang, X. Xu, Y. Li et al., “Systems pharmacology uncovers Janus functions of botanical drugs: activation of host defense system and inhibition of influenza virus replication,” Integrative Biology, vol. 5, no. 2, pp. 351–371, 2013.
View at: Publisher Site | Google Scholar
Y.-P. Zhu and H. J. Woerdenbag, “Traditional Chinese herbal medicine,” Pharmacy World and Science, vol. 17, no. 4, pp. 103–112, 1995.
View at: Publisher Site | Google Scholar
Z.-A. Yina, Y.-Y. He, L. Long et al., “Prevention and treatment idea of Chinese medicine on ebola hemorrhagic fever,” Journal of Acupuncture and Herbs, vol. 2, no. 2014, pp. 55–62, 2014.
View at: Google Scholar
S. Dharmananda, “Treatment of leukemia using integrated Chinese and western medicine,” International Journal of Oriental Medicine, vol. 22, pp. 169–183, 1997.
View at: Google Scholar
J. Cai and Y. Zhen, “Medicine in ancient China,” in Medicine Across Cultures, pp. 49–73, Springer, Berlin, Germany, 2003.
View at: Google Scholar
S. Dharmananda, Sun Simiao, vol. 1, Deutsche Übersetzung in ZTCM, 2004.
S. Jing and Z. Xiaojiang, “Professor Zhang Shiqing's experience on treating pediatric allergic purpura with Xijiao Dihuang decoction,” Journal of Gansu College of Traditional Chinese Medicine, vol. 3, article 3, 2005.
View at: Google Scholar
B.-R. Zheng, J.-P. Shen, H.-F. Zhuang, S.-Y. Lin, Y.-P. Shen, and Y.-H. Zhou, “Treatment of severe aplastic anemia by immunosuppressor anti-lymphocyte globulin/anti-thymus globulin as the chief medicine in combination with Chinese drugs,” Chinese Journal of Integrative Medicine, vol. 15, no. 2, pp. 145–148, 2009.
View at: Publisher Site | Google Scholar
S. I. Berger and R. Iyengar, “Network analyses in systems pharmacology,” Bioinformatics, vol. 25, no. 19, pp. 2466–2472, 2009.
View at: Publisher Site | Google Scholar
C. Huang, C. Zheng, Y. Li, Y. Wang, A. Lu, and L. Yang, “Systems pharmacology in drug discovery and therapeutic insight for herbal medicines,” Briefings in Bioinformatics, vol. 15, no. 5, Article ID bbt035, pp. 710–733, 2013.
View at: Publisher Site | Google Scholar
X. Li, X. Xu, J. Wang et al., “A system-level investigation into the mechanisms of Chinese traditional medicine: compound danshen formula for cardiovascular disease treatment,” PLoS ONE, vol. 7, no. 9, Article ID e43918, 2012.
View at: Publisher Site | Google Scholar
W. Tao, X. Xu, X. Wang et al., “Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease,” Journal of Ethnopharmacology, vol. 145, no. 1, pp. 1–10, 2013.
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, article 13, 2014.
View at: Publisher Site | Google Scholar
X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Sciences, vol. 13, no. 6, pp. 6964–6982, 2012.
View at: Publisher Site | Google Scholar
W. P. Walters and M. A. Murcko, “Prediction of ‘drug-likeness’,” Advanced Drug Delivery Reviews, vol. 54, no. 3, pp. 255–271, 2002.
View at: Publisher Site | Google Scholar
H. Liu, J. Wang, W. Zhou, Y. Wang, and L. Yang, “Systems approaches and polypharmacology for drug discovery from herbal medicines: an example using licorice,” Journal of Ethnopharmacology, vol. 146, no. 3, pp. 773–793, 2013.
View at: Publisher Site | Google Scholar
D. S. Wishart, C. Knox, A. C. Guo et al., “DrugBank: a comprehensive resource for in silico drug discovery and exploration,” Nucleic Acids Research, vol. 34, supplement 1, pp. D668–D672, 2006.
View at: Publisher Site | Google Scholar
K. S. Pang, “Modeling of intestinal drug absorption: roles of transporters and metabolic enzymes (for the gillette review series),” Drug Metabolism and Disposition, vol. 31, no. 12, pp. 1507–1519, 2003.
View at: Publisher Site | Google Scholar
A. Kam, K. M. Li, V. Razmovski-Naumovski et al., “The protective effects of natural products on blood-brain barrier breakdown,” Current Medicinal Chemistry, vol. 19, no. 12, pp. 1830–1845, 2012.
View at: Publisher Site | Google Scholar
L. Li, Y. Li, Y. Wang, S. Zhang, and L. Yang, “Prediction of human intestinal absorption based on molecular indices,” Journal of Molecular Science, vol. 23, pp. 286–291, 2007.
View at: Google Scholar
A.-L. Boulesteix, “PLS dimension reduction for classification with microarray data,” Statistical Applications in Genetics and Molecular Biology, vol. 3, no. 1, 2004.
View at: Publisher Site | Google Scholar | MathSciNet
D. Chung and S. Keles, “Sparse partial least squares classification for high dimensional data,” Statistical Applications in Genetics and Molecular Biology, vol. 9, no. 1, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
C. Knox, V. Law, T. Jewison et al., “DrugBank 3.0: a comprehensive resource for ‘omics’ research on drugs,” Nucleic Acids Research, vol. 39, supplement 1, pp. D1035–D1041, 2011.
View at: Publisher Site | Google Scholar
H. Yu, J. Chen, X. Xu et al., “A systematic prediction of multiple drug-target interactions from chemical, genomic, and pharmacological data,” PLoS ONE, vol. 7, no. 5, Article ID e37608, 2012.
View at: Publisher Site | Google Scholar
C. Zheng, Z. Guo, C. Huang et al., “Large-scale direct targeting for drug repositioning and discovery,” Scientific Reports, vol. 5, 2015.
View at: Publisher Site | Google Scholar
C. H. Wu, R. Apweiler, A. Bairoch et al., “The Universal Protein Resource (UniProt): an expanding universe of protein information,” Nucleic Acids Research, vol. 34, supplement 1, pp. D187–D191, 2006.
View at: Publisher Site | Google Scholar
G. Bindea, B. Mlecnik, H. Hackl et al., “ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks,” Bioinformatics, vol. 25, no. 8, pp. 1091–1093, 2009.
View at: Publisher Site | Google Scholar
M. E. Smoot, K. Ono, J. Ruscheinski, P.-L. Wang, and T. Ideker, “Cytoscape 2.8: new features for data integration and network visualization,” Bioinformatics, vol. 27, no. 3, Article ID btq675, pp. 431–432, 2011.
View at: Publisher Site | Google Scholar
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, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
C. J. M. Kane, J. H. Menna, C.-C. Sung, and Y.-C. Yeh, “Methyl gallate, methyl-3,4,5-trihydroxybenzoate, is a potent and highly specific inhibitor of herpes simplex virus in vitro. II. Antiviral activity of methyl gallate and its derivatives,” Bioscience Reports, vol. 8, no. 1, pp. 95–102, 1988.
View at: Publisher Site | Google Scholar
T. Grosser, E. Smyth, and G. FitzGerald, “Antiinflammatory, antipyretic and analgesic agents,” in Goodman and Gilman's the Pharmacological Basis of Therapeutics, L. Bruton, B. A. Chabner, and B. C. Knollmann, Eds., p. 985, McGraw-Hill Education, 2011.
View at: Google Scholar
X.-X. Xu, X.-M. Qi, W. Zhang et al., “Effects of total glucosides of paeony on immune regulatory toll-like receptors TLR2 and 4 in the kidney from diabetic rats,” Phytomedicine, vol. 21, no. 6, pp. 815–823, 2014.
View at: Publisher Site | Google Scholar
L. L. Pan and M. Dai, “Paeonol from Paeonia suffruticosa prevents TNF-α-induced monocytic cell adhesion to rat aortic endothelial cells by suppression of VCAM-1 expression,” Phytomedicine, vol. 16, no. 11, pp. 1027–1032, 2009.
View at: Publisher Site | Google Scholar
P. Prakash and N. Gupta, “Therapeutic uses of Ocimum sanctum Linn (Tulsi) with a note on eugenol and its pharmacological actions: a short review,” Indian Journal of Physiology and Pharmacology, vol. 49, no. 2, pp. 125–131, 2005.
View at: Google Scholar
K. Markowitz, M. Moynihan, M. Liu, and S. Kim, “Biologic properties of eugenol and zinc oxide-eugenol: a clinically oriented review,” Oral Surgery, Oral Medicine, Oral Pathology, vol. 73, no. 6, pp. 729–737, 1992.
View at: Publisher Site | Google Scholar
B. Lomenick, R. W. Olsen, and J. Huang, “Identification of direct protein targets of small molecules,” ACS Chemical Biology, vol. 6, no. 1, pp. 34–46, 2010.
View at: Publisher Site | Google Scholar
K. M. Lee, K. W. Lee, S. K. Jung et al., “Kaempferol inhibits UVB-induced COX-2 expression by suppressing Src kinase activity,” Biochemical Pharmacology, vol. 80, no. 12, pp. 2042–2049, 2010.
View at: Publisher Site | Google Scholar
J. L. Pang, D. A. Ricupero, S. Huang et al., “Differential activity of kaempferol and quercetin in attenuating tumor necrosis factor receptor family signaling in bone cells,” Biochemical Pharmacology, vol. 71, no. 6, pp. 818–826, 2006.
View at: Publisher Site | Google Scholar
C. Shi, F. Wu, X. Zhu, and J. Xu, “Incorporation of β-sitosterol into the membrane increases resistance to oxidative stress and lipid peroxidation via estrogen receptor-mediated PI3K/GSK3β signaling,” Biochimica et Biophysica Acta (BBA)—General Subjects, vol. 1830, no. 3, pp. 2538–2544, 2013.
View at: Publisher Site | Google Scholar
S. Loizou, I. Lekakis, G. P. Chrousos, and P. Moutsatsou, “β-Sitosterol exhibits anti-inflammatory activity in human aortic endothelial cells,” Molecular Nutrition and Food Research, vol. 54, no. 4, pp. 551–558, 2010.
View at: Publisher Site | Google Scholar
E. M. Leroy, S. Baize, V. E. Volchkov et al., “Human asymptomatic Ebola infection and strong inflammatory response,” The Lancet, vol. 355, no. 9222, pp. 2210–2215, 2000.
View at: Publisher Site | Google Scholar
D. S. Reed, L. E. Hensley, J. B. Geisbert, P. B. Jahrling, and T. W. Geisbert, “Depletion of peripheral blood T lymphocytes and NK cells during the course of Ebola hemorrhagic fever in cynomolgus macaques,” Viral Immunology, vol. 17, no. 3, pp. 390–400, 2004.
View at: Publisher Site | Google Scholar
Y. S. Kang, Y. G. Park, B. K. Kim et al., “Angiotensin II stimulates the synthesis of vascular endothelial growth factor through the p38 mitogen activated protein kinase pathway in cultured mouse podocytes,” Journal of Molecular Endocrinology, vol. 36, no. 2, pp. 377–388, 2006.
View at: Publisher Site | Google Scholar
A. Srikiatkhachorn, C. Ajariyakhajorn, T. P. Endy et al., “Virus-induced decline in soluble vascular endothelial growth receptor 2 is associated with plasma leakage in dengue hemorrhagic fever,” Journal of Virology, vol. 81, no. 4, pp. 1592–1600, 2007.
View at: Publisher Site | Google Scholar
L. A. Moraes, L. Piqueras, and D. Bishop-Bailey, “Peroxisome proliferator-activated receptors and inflammation,” Pharmacology & Therapeutics, vol. 110, no. 3, pp. 371–385, 2006.
View at: Publisher Site | Google Scholar
K. Wojciechowska-Durczynska, K. Krawczyk-Rusiecka, A. Cyniak-Magierska, A. Zygmunt, S. Sporny, and A. Lewinski, “The role of phosphoinositide 3-kinase subunits in chronic thyroiditis,” Thyroid Research, vol. 5, no. 1, article 22, 4 pages, 2012.
View at: Publisher Site | Google Scholar
S. Fisher-Hoch, J. B. McCormick, D. Sasso, and R. B. Craven, “Hematologic dysfunction in Lassa fever,” Journal of Medical Virology, vol. 26, no. 2, pp. 127–135, 1988.
View at: Publisher Site | Google Scholar
M.-S. Kim and S.-H. Kim, “Inhibitory effect of astragalin on expression of lipopolysaccharideinduced inflammatory mediators through NF-κB in macrophages,” Archives of Pharmacal Research, vol. 34, no. 12, pp. 2101–2107, 2011.
View at: Publisher Site | Google Scholar
A. N. Daniel, S. M. Sartoretto, G. Schmidt, S. M. Caparroz-Assef, C. A. Bersani-Amado, and R. K. N. Cuman, “Anti-inflammatory and antinociceptive activities of eugenol essential oil in experimental animal models,” Revista Brasileira de Farmacognosia, vol. 19, no. 1, pp. 212–217, 2009.
View at: Publisher Site | Google Scholar
T.-C. Chou, “Anti-inflammatory and analgesic effects of paeonol in carrageenan-evoked thermal hyperalgesia,” British Journal of Pharmacology, vol. 139, no. 6, pp. 1146–1152, 2003.
View at: Publisher Site | Google Scholar
J. P. Wood, J. R. Silveira, N. M. Maille, L. M. Haynes, and P. B. Tracy, “Prothrombin activation on the activated platelet surface optimizes expression of procoagulant activity,” Blood, vol. 117, no. 5, pp. 1710–1718, 2011.
View at: Publisher Site | Google Scholar
P. P. Tak and G. S. Firestein, “NF-κB: a key role in inflammatory diseases,” The Journal of Clinical Investigation, vol. 107, no. 1, pp. 7–11, 2001.
View at: Publisher Site | Google Scholar
W.-L. Jiang, X.-G. Chen, H.-B. Zhu, Y.-B. Gao, J.-W. Tian, and F.-H. Fu, “Paeoniflorininhibits systemic inflammation and improves survival in experimental sepsis,” Basic & Clinical Pharmacology and Toxicology, vol. 105, no. 1, pp. 64–71, 2009.
View at: Publisher Site | Google Scholar
M. I. Cybulsky, J. W. U. Fries, A. J. Williams et al., “Gene structure, chromosomal location, and basis for alternative mRNA splicing of the human VCAM1 gene,” Proceedings of the National Academy of Sciences of the United States of America, vol. 88, no. 17, pp. 7859–7863, 1991.
View at: Publisher Site | Google Scholar
G. Marfe, M. Tafani, M. Indelicato et al., “Kaempferol induces apoptosis in two different cell lines via Akt inactivation, bax and SIRT3 activation, and mitochondrial dysfunction,” Journal of Cellular Biochemistry, vol. 106, no. 4, pp. 643–650, 2009.
View at: Publisher Site | Google Scholar
W.-T. Kuo, Y.-C. Tsai, H.-C. Wu et al., “Radiosensitization of non-small cell lung cancer by kaempferol,” Oncology Reports, vol. 34, no. 5, pp. 2351–2356, 2015.
View at: Publisher Site | Google Scholar
H. S. Lee, H. J. Cho, G. T. Kwon, and J. H. Park, “Kaempferol downregulates insulin-like growth factor-I receptor and ErbB3 signaling in HT-29 human colon cancer cells,” Journal of Cancer Prevention, vol. 19, no. 3, pp. 161–169, 2014.
View at: Publisher Site | Google Scholar
M. F. Saeed, A. A. Kolokoltsov, A. N. Freiberg, M. R. Holbrook, and R. A. Davey, “Phosphoinositide-3 kinase-akt pathway controls cellular entry of ebola virus,” PLoS Pathogens, vol. 4, no. 8, Article ID e1000141, 2008.
View at: Publisher Site | Google Scholar
J. P. Cooke, “NO and angiogenesis,” Atherosclerosis Supplements, vol. 4, no. 4, pp. 53–60, 2003.
View at: Publisher Site | Google Scholar
D. Fulton, J.-P. Gratton, T. J. McCabe et al., “Regulation of endothelium-derived nitric oxide production by the protein kinase Akt,” Nature, vol. 399, no. 6736, pp. 597–601, 1999.
View at: Publisher Site | Google Scholar
S. Babaei, K. Teichert-Kuliszewska, Q. Zhang, N. Jones, D. J. Dumont, and D. J. Stewart, “Angiogenic actions of angiopoietin-1 require endothelium-derived nitric oxide,” The American Journal of Pathology, vol. 162, no. 6, pp. 1927–1936, 2003.
View at: Publisher Site | Google Scholar
S. Jing and Z. Xiaojiang, “Professor Zhang Shiqing's experience on treating pediatric allergic purpura with Xijiao Dihuang decoction,” Journal of Gansu College of Traditional Chinese Medicine, vol. 3, article 003, 2005.
View at: Google Scholar

Copyright

Copyright © 2016 Jianling Liu 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

Views

4449

Downloads

1460

Citations