Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 22 November 2019
Sec. Ethnopharmacology
This article is part of the Research Topic Network Pharmacology and Traditional Medicine View all 31 articles

Comparative Network Pharmacology Analysis of Classical TCM Prescriptions for Chronic Liver Disease

  • 1Department of Gastroenterology, Shanghai Tenth People's Hospital, School of Life Sciences and Technology, Tongji University, Shanghai, China
  • 2Institute of Interdisciplinary Integrative Medicine Research, Shanghai University of Traditional Chinese Medicine, Shanghai, China
  • 3Department of Pharmacology, Shanghai University of Traditional Chinese Medicine, Shanghai, China
  • 4Key Laboratory of Liver and Kidney Diseases of Ministry of Education of China, Institute of Liver Diseases, Shuguang Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China

Chronic liver disease (CLD) has become a major global health problem while herb prescriptions are clinically observed with significant efficacy. Three classical Traditional Chinese Medicine (TCM) formulae, Yinchenhao Decoction (YCHT), Huangqi Decoction (HQT), and Yiguanjian (YGJ) have been widely applied in China to treat CLD, but no systematic study has yet been published to investigate their common and different mechanism of action (MOA). Partial limitation may own to deficiency of effective bioinformatics methods. Here, a computational framework of comparative network pharmacology is firstly proposed and then applied to herbal recipes for CLD disease. The analysis showed that, the three formulae modulate CLD mainly through functional modules of immune response, inflammation, energy metabolism, oxidative stress, and others. On top of that, each formula can target additional unique modules. Typically, YGJ ingredients can uniquely target the ATP synthesis and neurotransmitter release cycle. Interestingly, different formulae may regulate the same functional module in different modes. For instance, YCHT and YGJ can activate oxidative stress-related genes of SOD family while HQT are found to inhibit SOD1 gene. Overall, our framework of comparative network pharmacology proposed in our work may not only explain the MOA of different formulae treating CLD, but also provide hints to further investigate the biological basis of CLD subtypes.

Introduction

With increasing alcohol consumption and viral infection, chronic liver disease (CLD) has become a worldwide health concern. Chronic liver disease in the clinical context is a disease process of the liver that involves a process of progressive destruction and regeneration of the liver parenchyma leading to fibrosis and cirrhosis. “Chronic liver disease” refers to disease of the liver that lasts over a period of six months. It consists of a wide range of liver pathologies covering inflammation (chronic hepatitis), liver cirrhosis, and hepatocellular carcinoma. According to causes, chronic liver diseases can be classified as viral disease, toxic and drug-related disease, metabolic disease, autoimmune-response disease, and so on. While the nucleotide analogues (NUCs) and pegylated-interferon (Peg-IFN) therapies are effective, patients may not tolerate their adverse effects in lifelong treatment. In the meantime, most of CLD patients seek medical help from herbal medicine for its significant efficacy and low side effects. Recently, several meta-studies of randomized, controlled, clinical trials (RCTs) suggested that herbal remedies have a similar and even better effects compared to western drug therapy in dealing with CLD (Peng et al., 2016; Chen et al., 2019).

According to Guidelines for the diagnosis and treatment of liver fibrosis in integrative medicine practice of China, three classical formulae, Yinchenhao Decoction (YCHT), Huangqi Decoction (HQT), and Yiguanjian (YGJ) have been recommended to treat different subtypes of CLD according to “Zheng” differentiation, which classifies the pathological state of a patient at system level leading to individualized treatments. Syndrome Differentiation and Treatment (SDAT) is the basis of clinical application of traditional Chinese Medicine prescription. A chronic liver disease can be treated with the different efficacy TCM prescription based on its different TCM syndrome. Formula YCHT is suggested to treat “damp-heat” type of patients featured by jaundice, inflammation, and abnormal fatty acids metabolism (Yan et al., 2017), and the therapeutic characteristic of Yinchenhao decoction is to clear heat and drain dampness. Clinically, it is applied mainly to treat patients who have acute jaundice hepatitis, hepatitis B, cholelithiasis, cholecystitis, leptospirosis, typhoid fever, pneumonia, hemolytic jaundice caused by favism, and so on (Yan et al., 2017), while formula HQT, with the proportion of Huangqi [Astragalus mongholicus Bge.] and Zhigancao [Glycyrrhiza uralensis Fisch.] at 6:1, is recommended for those “Qi-deficiency” patients with weak body, black and yellowish complexion, palpitations and vexations in the chest, dry lips and mouths, pale complexion and lack of appetite. The therapeutic characteristic of HQT is tonifying qi. HQT is applied clinically to treat CLD patients who have liver fibrosis or cirrhosis, heart failure, constipation, type 2 diabetes, and diabetic peripheral neuropathy with “Qi-blood deficiency” (Wang et al., 2015). Additionally, HQT alleviated DMN-induced liver fibrosis (Song et al., 2016). As the therapeutic characteristic of YGJ is nourishing liver and kidney, it is applied clinically to treat those who have chronic liver diseases, chronic gastritis, gastric and duodenal ulcer, intercostal neuralgia, and neurosis with “Yin-deficiency” (Wang et al., 2015).

Studies have investigated the liver-protective activity of above formula. For instance, being composed of Yinchen [Artemisia capillaris Thunb.], Zhizi [Gardenia jasminoides J. Ellis], and Dahuang [Rheum palmatum L.], YCHT demonstrated variable ability of liver protection in treating cholestasis, liver fibrosis, hepatitis, biliary cirrhosis, and cholesteric liver diseases (Yan et al., 2017). The liver protection of YCHT was also confirmed by histopathology and biochemical experiments (Sun et al., 2019), and YCHD treatment could reverse the damage by DMN-induced in liver function (Cai et al., 2018). According to literature validation, pharmacological research revealed that YCHT showed function of anti-oxidative stress, inhibiting apoptosis. Also, it was validated to regulate cellular inflammation and lipid metabolism (Wang et al., 2015). Kupffer cells (KC) activation of inflammatory response and abnormal fatty acid metabolism were found to be the main pathological basis in the early inflammation for YCHT Syndrome (Xu-Dong et al., 2014).

Through literature searching, HQT decoction was reported to downregulate the expressions of PDGF genes, collagen genes (COL1A1, COL1A2, COL5A2), and THBS1, inhibiting TGF-beta and PDGF signaling pathways, followed by verification via qRT-PCR (Zhang et al., 2015). Also, it was noticed to promote Kupffer cell activation (Du et al., 2012; Liu et al., 2012; Zhang et al., 2015; Song et al., 2016), inhibit Notch signaling pathway (Wang et al., 2015), and alleviate oxidative stress and lipid peroxidation injury. In an animal experiment, HQT was found to significantly reduce alpha-naphthylisothiocyanate-induced cholestasis in mice (Wu et al., 2017). At herb level, herb Huangqi in HQT were reported with anti-oxidative and immune regulating effects (Shahzad et al., 2016; Li et al., 2019).

YGJ decoction has 6 herbs: Beishashen [Glehnia littoralis F. Schmidt ex Miq.], Maidong [Ophiopogon japonicus (Thunb.) Ker Gawl.], Danggui [Angelica sinensis (Oliv.) Diels], Shengdihuang [Rehmannia glutinosa (Gaertn.) DC.], Gouqizi [Lycium barbarum L.], and Chuanlianzi [Melia azedarach L.]. Research indicates the hepatoprotective, anti-fibrogenic, anti-angiogenesis effects of YGJ in hepatic injury mice models (Tian et al., 2016). In rat model of liver fibrosis, YGJ was found to alleviate the hepatic collagen hyperplasia and deteriorated hepatic function (Tao et al., 2009). Also, YGJ can reduce the liver oxidative stress and lipid peroxidation injury while inhibiting angiogenesis and induce cell differentiation (Zhou et al., 2015). In a rat model of liver cirrhosis, YGJ was reported with repairing function of liver cirrhosis, and the key mechanism was suggested as relating to the regulation of macrophage activation state (Xu et al., 2018).

Above studies showed the biological activity of formulae, but it would be interesting to know the detailed MOA difference or similarities between the three formulae under the same background of CLD disease. As TCM formulae are believed to produce efficacy in a holistic way, the difficulty in detecting the MOA similarities and differences lie in the inherent complexity of multi-ingredients and multi-targets for TCM formulae.

With the development of systems biology and accumulated herbal targets information, the emerging method of network pharmacology makes it possible for MOA analysis of herbal ingredients. Coupled with transcriptomics validation, a latest article provided an excellent model adopting network pharmacology to predict the active components and regulation mechanisms of YCHT against liver fibrosis (Cai et al., 2019). Until now, network pharmacology has been applied to study individual ingredients, herbs, or formulae. One reason causing the difficulty is the lack of effective bioinformatics methods. Here, we developed a computational framework aiming to compare the MOA between formulae for same disease, taking the YCHT, HQT, and YGJ as an example.

The results indicated that although different formulae have different molecular targeting profile, they are commonly involved in a group of functional modules, indicating the CLD disease background. In addition to that, each formula was detected with unique function modules in MOA analysis. More interestingly, even on the same functional module, different regulating modes were found for different formulae agreeing with observed “Zheng” classification of CLD subtyping.

Materials and Methods

Dataset

Formulae, Herbs, Ingredients, and Targets

For each of the three formulae, herb components, ingredients and targets of each ingredients were collected from HIT (Ye et al., 2011), NPASS (Zeng et al., 2018b), TCMDB (Chen, 2011) and TCM-ID (Huang et al., 2018) databases respectively.

Method

Disease Network, Disease Modules and Function Annotation

Targets of herbal active ingredients were obtained from HIT (Ye et al., 2011) and NPASS (Zeng et al., 2018b) database for target network construction. Potential protein-protein interaction (PPI) network of disease is generated via Reactome database (Croft et al., 2014) based on combined unique targets from 3 formulae via step size 1 (one bridging node). Then, the disease network was modularized through ReactomeFIViz (Wu and Haw, 2017). Genecards database (Safran et al., 2010) was used for function annotation of each module because of its high level of function annotation. Disease network and modularization were displayed via Cytoscape (Shannon et al., 2003). KEGG (Kanehisa, 2002) and GO Biological Process (Ashburner et al., 2000; Mi et al., 2017; The Gene Ontology, 2017) databases were used for biological function annotation for targets of each formula.

Statistic Enrichment of Functional Modules

For the given number of formula targets, we tested whether it is significantly enriched in a specific module A, comparing to a randomly picked module with equal size of A, from nodes of disease network background. To make the variation comparable in statistics, we repeated the picking for 5000 times. And each time, half the number of given targets were randomly picked. T-test statistics was made between the hitting of tested targets on specific module and the random module equally sized. According to above, different number of targets for each formula were enriched respectively into the functional modules of disease.

Effect of Kaempferol on the Absorption of Astragalosides

To investigate the effect of kaempferol on the absorption of astragalosides, experiments of verted rat sac was performed as previously described (Yan et al., 2015). Briefly, the everted intestinal sacs were prepared by rapidly removing the small intestine from starved rats euthanized under carbon dioxide (CO2) anesthesia. Then, four intestinal segments (the duodenum, upper jejunum, bottom jejunum, and ileum) were excised, flushed several times with saline solution at room temperature, and then placed immediately into oxygenated buffer solution at 37°C. Intestine segments from four rats were divided into four groups treated with astragalosides without or with kaempferol at three different concentrations (1, 3, 10M) using a 4 × 4 Latin square design (Gressley et al., 2016).

www.frontiersin.org

The concentration of Astragaloside solutions was set at 8 mg/mL based on the gastrointestinal concentration after oral administration of Huangqi [Astragalus mongholicus Bge.] decoction to rats, and the concentration of the additional kaempferol was 1, 3, and 10 µM. Then, aliquots of the sac fluid sample (200 µL) were removed, and the same volume of buffer solution was added at 0, 0.5, 1, and 2 h. The area of each sac was calculated closely. Furthermore, each sac was weighed before and after fluid collection to calculate accurately the volume inside the sac. The samples were analyzed using UPLC-LTQ-Orbitrap in both positive and negative ionization mode as previously described. The proposed method could be used to determine the 10 compounds in astragalosides simultaneously with high precision, sensitivity, and accuracy (Zeng et al., 2018a). A multivariate analysis of variance (ANOVA) was performed using the SPSS ver. 13.0 software while p-values < 0.05 were considered statistically significant.

Results

Computational Framework of Comparative Network Pharmacology

The idea of comparative network pharmacology is to construct a background disease network covering all the targets of different formulae, and then analyze the common and different modules or regulations. The framework of comparative network pharmacology is proposed as below (Figure 1): 1. Disease network construction and modularization of combined targets; 2. Comparative analysis of formulae on disease module network; 3. Detailed regulation comparison between formulae.

FIGURE 1
www.frontiersin.org

Figure 1 Computational Framework of comparative network pharmacology.

Accordingly, the computational framework works like this: 1. Retrieve targets of each formula for same disease; 2. Construct the disease network based on combined targets of all formulae; 3. Partition and annotate the disease network into function modules; 4. Calculate the enriched modules for common targets among all formulae and unique targets for each formula; 5. Analyze the MOA similarities and differences.

Disease Network and Function Modules

The known ingredients with experimental target evidence from literatures were retrieved from databases and summarized into Supplementary Table 1 for each herb and formula. Though different formulae are composed of different number of herbs, the overall known targets for each formula is roughly similar ranging from 120 to 140 genes. A total of 293 non-redundant targets were combined for all 3 formulae.

Targets of each formula was firstly enriched into KEGG and GO database for function annotation. The results show that all the three formulae are primarily targeting pathway of Hepatitis B, suggesting their common therapeutic effects in liver disease, but the enriched pathway lists are heavily overlapping (Supplementary Table 2). Similarly, mapping to GO databases also shows similar pathways profiles and GO terms. This suggests that the commonly used bioinformatics analysis of KEGG and GO annotation is good at inferring the MOA similarities, rather than the MOA differences for different formulae.

After running the pipeline of comparative network pharmacology, we proposed the biggest connected subnetwork containing 244 targets (83.28% of total) as the CLD disease background network.

This CLD network was initially clustered into 16 modules by ReactomeFIViz (Wu and Haw, 2017). After removing the two modules fully composed of background bridging nodes, we kept the remaining 14 modules for further analysis.

Common and Unique Functional Modules for TCM Formulae

Since active ingredient list is often an arguable issue for herbs, consensus has not been reached particularly for those commonly found among herbs, such as kaempferol, vitamins, or quercetins. Despite the molecular activity in vitro, the low bioavailability was often reported (Barve et al., 2009). Under the complexity of multi-ingredients context of herbs, their potential ADME effect to known active ingredients deserves further investigation. Taking kaempferol as an example, we tested the effects of kaempferol on the absorption of astragalosides in Huangqi herb via verted rat sac method under 4 different dosages. In the preliminary test, it was found that there were significant differences in the absorption of astragalosides in different rats and different gut segments. Then, Latin square was adopted to exclude the influence of different rats and intestinal segments. As the detailed results shown in Supplementary Table 3, under the concentration range of 1-10µM, kaempferol was observed with no significant effect on the absorption of 10 components in astragalosides in verted rat sac model.

As such, the active ingredient list was perturbed from full to shortlist by removing kaempferol, vitamin C, and vitamin E one by one. The module enrichment results were summarized in Supplementary Table 1 (Sheet name: results of enrichment). Despite the perturbation, all formulae share 3 consensus modules and one conditional module, and each formula own two unique modules. Ten modules covered 236 targets, and their functional annotation are shown in Table 1 by GeneCards (Safran et al., 2010).

TABLE 1
www.frontiersin.org

Table 1 List of functional modules for CLD disease network derived from three TCM formulae.

Through the computational framework, the MOA similarity and differences are illustrated in Figure 2 as shown below in the background of disease modules. The four common modules cover 167 targets, mainly involving oxidative stress, drug metabolism, DNA damage, NFAT and Cardiac Hypertrophy, Jak-STAT signaling and inflammation, etc. Further details of modules are shown in Figure 3. These common processes of oxidative stress, inflammation, and corresponding signaling pathways are highly involved in the pathological basis of CLD. For instance, oxidative stress, caused by imbalance between production and consumption by antioxidants, is a prominent feature in the pathophysiology of CLD (Vuppalanchi et al., 2011). Once oxidative stress is initiated, continuous cycle of cellular damage and release of proinflammatory cytokines will move on leading to hepatic inflammation, fibrosis, and cirrhosis (Comporti et al., 2009).

FIGURE 2
www.frontiersin.org

Figure 2 Common and unique functional modules between the three formulae. Each cycle represents a functional module, and the size of each cycle is correlated to the number of genes in each module. Modules 8, 9, and 10 are illustrated as common module, while module 7 is meant as conditional common module. Herbs are labelled as well indicating that their ingredients can target unique module.

FIGURE 3
www.frontiersin.org

Figure 3 The overall network modularization. Green nodes indicate unique targets of YCHT; blue nodes indicate unique targets of HQT; pink nodes indicate unique targets of YGJ; and yellow nodes indicate overlapping targets from 2 or 3 formulae.

In view of literature report, consistent mechanistic evidences have been found between our computing analysis and independent validation. Modules 8 and 7 are involved with anti-oxidative stress. All three formulae can target these two common modules, indicating their common anti-oxidative function. Furthermore, module 10 is mainly involved with inflammation, and all three formulae were predicted to interact with a wide variety of inflammation factors. For instance, YCHT can target IL1B, IL4, and IL6; HQT can target IL6 and IL10; and YGJ can target more factors of IL1B, IL2, IL4, IL5, IL6, IL10, and IL13. At molecular level, our predictive results gained support from previous validation in that TGF-beta was targeted by HQT, and collagen genes (COL1A1, COL4A2, and COL7A1) were targeted by YGJ.

Meanwhile, each formula can significantly target two unique modules as illustrated in Figure 2. YCHT can significantly target modules 1 and 2. Module 1 is specifically involved with immune response of lectin induced complement pathway, where Dahuang herb can uniquely interact with it. Module 2 is related to leukocyte trans-endothelial migration and cell membrane interaction or adhesion, where all herbs in YCHT can target this module. HQT compounds can interact with CDK-mediated phosphorylation and removal of Cdc6, together with ATP binding in LKB1 signaling events (modules 3 and 4). Most interestingly, YGJ ingredients can directly and uniquely target the ATP synthesis (Kamanna et al., 2013) and neurotransmitter release cycle (modules 5 and 6), implying the potential modulation of ATP energy and neuron alertness.

Different Regulation Modes in the Same Functional Module

It is noted that, despite the unique modules for different formulae, the majority of targets are still involving the same functional modules. We investigated the regulation modes of herbal ingredients in each formula on the same module and found that a substantial number of targets are differently regulated by different formulae. It’s realized that the information of regulation modes on targets is highly precious, where only a few databases such as HIT collected the sparse evidences. Details of regulation modes on three consensus modules can be found in Table 2.

TABLE 2
www.frontiersin.org

Table 2 Different regulation modes of different formulae on the same functional module. Information of regulation modes is collected from HIT database (Ye et al., 2011).

For the module of DNA Damage and Oxidative Stress (Module 8), we can see that YCHT and YGJ ingredients can activate the antioxidant-related genes SOD families, while HQT compounds inhibit SOD1 genes. This actually agrees well with previous literature results for HQT (Wang et al., 2008). On Module 9 of NFAT and Cardiac Hypertrophy, YCHT ingredient can inhibit ESR1 target and HQT ingredient has inverse action, while YGJ shows mixed function. The 10th Module contains a list of inflammation factors. Generally, all formulae inhibit a wide range of pro-inflammation factors, such as IL1, IL2, IL4, IL5, etc.

Previous studies have shown that YCHT can regulate Kupffer cells; inhibit the release of proinflammatory cytokines, such as TNF-α, IL-1β, and IL-6; and promote the release of anti-inflammatory and anti-fibrogenic factors, IFN-γ and IL-10 (Zhang et al., 2013). Palmitic acid contained in herb Zhizi can inhibit the cytokines involved in various inflammatory processes as an antagonist of CASP1 (Xu et al., 2016). Further, the active ingredients in YCHT, such as salicylic acid, rottlerin, and garlic acid can also reduce the expression of NOS2, TNF, and RELA, and then inhibit the inflammation (Hee Sun et al., 2008). Among the cytokines, IL10 is well known for its balanced function of inflammation-inhibition. On this particular gene of IL10, HQT shows negative effects, while YGT ingredients show positive effects. It is noted that only individual ingredients were reported with activating/inhibiting effect on molecular targets; the final in vivo effects may depend on the relative content and bioavailability of multiple ingredients in formula.

Discussion

Highly personalized prescriptions of different formulae are frequently applied to treat one same disease in different states. Comparing the MOA similarities and differences between formulae for same disease will provide us a good opportunity to study the pathogenesis of disease subtypes. However, as shown in this study, the commonly simple adoption of KEGG or GO annotation often gives more similarities than differences because the pathway maps or GO terms have been strictly predefined while they are actually densely cross-talked. We thus set up a computational framework of comparative network pharmacology for this purpose and illustrated an example to reveal MOA details of 3 herbal recipes for CLD disease.

In applying this framework to other formulae, there are some suggestions to get better understanding of the method. First, instead of the predicted targets, we strongly support to collect those clearly validated ingredients and targets with regulation modes. One reason is that the prediction error may expand and mess up the modularization part during the subsequent analysis. Another reason is, regulation modes are normally deficient in the results of computational prediction. That is why we have to abandon several elegant TCM databases, such as SymMap (Wu et al., 2018), TCMSP (Ru et al., 2014), or YaTCM (Li et al., 2018). Secondly, other databases or household PPI list can also be adopted to generate PPI network if the information in Reactome is not enough. Lastly, the statistics of module enrichment is elaborated here by comparing with a randomly picked module with equal size of targeting groups in disease network background.

The comparative analysis finds the MOA difference correlates well with the Zheng symptoms between the formulae. As we know, YCHT is used to treat CLD patients with jaundice or inflammation. The unique function module of YCHT is just related to lectin-induced complement pathway and leukocyte trans-endothelial migration, which might mobilize extra anti-inflammation factors to deal with “damp-heat” symptoms. Another interesting finding here is that YGJ can directly target ATP synthesis and neurotransmitter release, which may promote ATP supply and increase the alertness of neuro systems. In our previous study, YGJ was found with ability to modulate the abnormal energy metabolism (Yan et al., 2017). This unique MOA of YGJ agrees well with its clinical application to treat “Yin deficiency” patients with extreme fatigue and weakness.

Yet, most of the time, disease subtypes of TCM Zheng are classified by highly summarized concepts of tongue color, pulse feeling, and personal overall phenotypes. While the Zheng research is still at beginning stage, our comparison framework from formula's perspective may provide alternative views to study their material basis. In addition to above, more consistent evidence may be proposed in future with both the progress in herbal target information and the TCM Zheng study.

Lastly, interaction between compounds and protein targets plays an important role in modulating the biological process and activity. We collected only those compounds with target information and activity evidence for network analysis. Yet, what ingredients are actually pharmacologically active deserves further investigation for these formulae. Correspondingly, perturbation tests are necessary to reach consensus conclusion. In fact, the list of active compounds may be affected by multiple factors such as the quantities in herbs, herb dosage, conjugated forms, compound interactions, biotransformation, dietary usage, bioavailability, other ADME issues (Dabeek and Marra, 2019).

To summarize, the comparative network pharmacology can detect a list of functional modules inferring the potential mechanisms of different formulae and provide functional insights for the subsequent exploration for not only formulae, but also disease subtyping. It's aware that all the information in the framework is collected from public database. Thus, the analysis results may be influenced by the abundancy and accuracy of related information. In the future, with more herbal experimental evidences accumulated, MOA comparison between different formulae is expected to be refreshed quantitatively and highlighted.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: www.ncbi.nlm.nih.gov/pmc/articles/PMC3013727/.

Ethics Statement

The experimental procedures involving the use of animals were approved by the Committee on the Use of Live Animals for Teaching and Research of the Shanghai University of Traditional Chinese Medicine and all experiments were performed in accordance with the guidelines of this committee.

Author Contributions

ZKC conducted the project and wrote the manuscript. YL and YW finished the experimental verification. XW, KT, and DW collected the data and interpreted the results. WZ modified the manuscript. YM, PL, and ZWC designed and supervised the project. All authors read, critically reviewed, and approved the final manuscript.

Funding

This work was supported in part by National Key R&D Program of China (grant number 2017YFC1700200 and 2017YFC0908400), the National Natural Science Foundation of China (31671379, 81573873, 81774196).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.01353/full#supplementary-material

Supplementary Table 1 | The detailed results of targets and ingredients.

Supplementary Table 2 | The detailed results of Enrichment in KEGG and GO.

Supplementary Table 3 | The detailed results of the everted intestinal sacs experiment.

References

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology Consortium. Nat. Genet. 25 (1), 25–29. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Barve, A., Chen, C., Hebbar, V., Desiderio, J., Saw, C. L., Kong, A. N. (2009). Metabolism, oral bioavailability and pharmacokinetics of chemopreventive kaempferol in rats. Biopharm. Drug Dispos. 30 (7), 356–365. doi: 10.1002/bdd.677

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, F. F., Wu, R., Song, Y. N., Xiong, A. Z., Chen, X. L., Yang, M. D., et al. (2018). Yinchenhao decoction alleviates liver fibrosis by regulating bile acid metabolism and TGF-beta/Smad/ERK signalling pathway. Sci. Rep. 8 (1), 15367. doi: 10.1038/s41598-018-33669-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, F. F., Bian, Y. Q., Wu, R., Sun, Y., Chen, X. L., Yang, M. D., et al. (2019). Yinchenhao decoction suppresses rat liver fibrosis involved in an apoptosis regulation mechanism based on network pharmacology and transcriptomic analysis. BioMed. Pharmacother. 114, 108863. doi: 10.1016/j.biopha.2019.108863

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Huang, C., Lu, J., Huang, R., Zhang, J., Chen, Z., et al. (2019). Efficacy and safety of Yinchen Sini decoction in treating biliary atresia patients after Kasai portoenterostomy: a systematic review and meta-analysis protocol. Med. (Baltimore) 98 (2), e13935. doi: 10.1097/MD.0000000000013935

CrossRef Full Text | Google Scholar

Chen, C. Y. (2011). TCM Database@Taiwan: the world's largest traditional Chinese medicine database for drug screening in silico. PloS One 6 (1), e15939. doi: 10.1371/journal.pone.0015939

PubMed Abstract | CrossRef Full Text | Google Scholar

Comporti, M., Arezzini, B., Signorini, C., Vecchio, D., Gardi, C. (2009). Oxidative stress, isoprostanes and hepatic fibrosis. Histol. Histopathol. 24 (7), 893–900. doi: 10.14670/HH-24.893

PubMed Abstract | CrossRef Full Text | Google Scholar

Croft, D., Mundo, A. F., Haw, R., Milacic, M., Weiser, J., Wu, G., et al. (2014). The Reactome pathway knowledgebase. Nucleic Acids Res. 42 (Database issue), D472–D477. doi: 10.1093/nar/gkt1102

PubMed Abstract | CrossRef Full Text | Google Scholar

Dabeek, W. M., Marra, M. V. (2019). Dietary Quercetin and Kaempferol: bioavailability and potential cardiovascular-related bioactivity in humans. Nutrients 11 (10), 2288. doi: 10.3390/nu11102288

CrossRef Full Text | Google Scholar

Du, J. X., Sun, M. Y., Du, G. L., Li, F. H., Liu, C., Mu, Y. P., et al. (2012). Ingredients of Huangqi decoction slow biliary fibrosis progression by inhibiting the activation of the transforming growth factor-beta signaling pathway. BMC Complement. Altern. Med. 12, 33. doi: 10.1186/1472-6882-12-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Gressley, T. F., Davison, K. A., Macies, J., Leonardi, C., McCarthy, M. M., Nemec, L. M., et al. (2016). Effect of abomasal carbohydrates and live yeast on measures of postruminal fermentation. J. Anim. Sci. 94 (1), 284–296. doi: 10.2527/jas.2015-9324

PubMed Abstract | CrossRef Full Text | Google Scholar

Hee Sun, B., Minho, W., Kyeong Ah, P., Young-Rae, K., Byung Lyul, C., Hyunji, L., et al. (2008). Prevention of TNF-induced necrotic cell death by rottlerin through a Nox1 NADPH oxidase. Exp. Mol. Med. 40 (2), 186. doi: 10.3858/emm.2008.40.2.186

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, L., Xie, D., Yu, Y., Liu, H., Shi, Y., Shi, T., et al. (2018). TCMID 2.0: a comprehensive resource for TCM. Nucleic Acids Res. 46 (D1), D1117–D1120. doi: 10.1093/nar/gkx1028

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamanna, V. S., Ganji, S. H., Kashyap, M. L. (2013). Recent advances in niacin and lipid metabolism. Curr. Opin. In Lipidol. 24 (3), 239–245. doi: 10.1097/MOL.0b013e3283613a68

CrossRef Full Text | Google Scholar

Kanehisa, M. (2002). The KEGG database. Novartis Found Symp. 247, 91–101. doi: 10.1002/0470857897.ch8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Ma, C., Zhao, X., Hu, Z., Du, T., Xu, X., et al. (2018). YaTCM: yet another traditional Chinese medicine database for drug discovery. Comput. Struct. Biotechnol. J. 16, 600–610. doi: 10.1016/j.csbj.2018.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W. K., Wang, G. F., Wang, T. M., Li, Y. Y., Li, Y. F., Lu, X. Y., et al. (2019). Protective effect of herbal medicine Huangqi decoction against chronic cholestatic liver injury by inhibiting bile acid-stimulated inflammation in DDC-induced mice. Phytomedicine 62, 152948. doi: 10.1016/j.phymed.2019.152948

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Wang, G., Chen, G., Mu, Y., Zhang, L., Hu, X., et al. (2012). Huangqi decoction inhibits apoptosis and fibrosis, but promotes Kupffer cell activation in dimethylnitrosamine-induced rat liver fibrosis. BMC Complement. Altern. Med. 12, 51. doi: 10.1186/1472-6882-12-51

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, H., Huang, X., Muruganujan, A., Tang, H., Mills, C., Kang, D., et al. (2017). PANTHER version 11: expanded annotation data from Gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 45 (D1), D183–D189. doi: 10.1093/nar/gkw1138

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, H., He, Y., Zheng, G., Zhang, W., Yao, Z., Xie, W. (2016). Meta-analysis of traditional herbal medicine in the treatment of nonalcoholic fatty liver disease. Cell Mol. Biol. (Noisy-le-grand) 62 (4), 88–95. doi: 10.14715/cmb/2016.62.4.16

PubMed Abstract | CrossRef Full Text | Google Scholar

Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 13. doi: 10.1186/1758-2946-6-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Safran, M., Dalah, I., Alexander, J., Rosen, N., Iny Stein, T., Shmoish, M., et al. (2010). GeneCards Version 3: the human gene integrator. Database (Oxford) 2010, baq020. doi: 10.1093/database/baq020

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahzad, M., Shabbir, A., Wojcikowski, K., Wohlmuth, H., Gobe, G. C. (2016). The antioxidant effects of Radix Astragali (Astragalus membranaceus and related species) in protecting tissues from injury and disease. Curr. Drug Targets 17 (12), 1331–1340. doi: 10.2174/1389450116666150907104742

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Y. N., Zhang, G. B., Lu, Y. Y., Chen, Q. L., Yang, L., Wang, Z. T., et al. (2016). Huangqi decoction alleviates dimethylnitrosamine-induced liver fibrosis: an analysis of bile acids metabolic mechanism. J. Ethnopharmacol. 189, 148–156. doi: 10.1016/j.jep.2016.05.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, H., Zhang, A. H., Yang, L., Li, M. X., Fang, H., Xie, J., et al. (2019). High-throughput chinmedomics strategy for discovering the quality-markers and potential targets for Yinchenhao decoction. Phytomedicine 54, 328–338. doi: 10.1016/j.phymed.2018.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, Q., Sun, M. Y., Feng, Q. (2009). Syndrome identification of CCl4 induced liver fibrosis model rats based on syndrome detecting from recipe used. Chin. J. Integr. Tradit. Western Med. 29 (3), 246–250. doi: CNKI:SUN:ZZXJ.0.2009-03-020

Google Scholar

The Gene Ontology, C. (2017). Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. 45 (D1), D331–D338. doi: 10.1093/nar/gkw1108

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, M., Liu, W., You, H., Zhao, Q., Ouyang, L., Gao, B., et al. (2016). Protective effect of Yiguanjian decoction against DNA damage on concanavalin A-induced liver injury mice model. J. Tradit. Chin. Med. 36 (4), 471–478. doi: 10.1016/s0254-6272(16)30064-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuppalanchi, R., Juluri, R., Bell, L. N., Ghabril, M., Kamendulis, L., Klaunig, J. E., et al. (2011). Oxidative stress in chronic liver disease: relationship between peripheral and hepatic measurements. Am. J. Med. Sci. 342 (4), 314–317. doi: 10.1097/MAJ.0b013e31821d9905

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Liu, P., Wang, C. S. (2008). Effects of 5 classical recipes on anti-oxidative stress in rat liver with cirrhosis. Zhongguo Zhong Xi Yi Jie He Za Zhi 28 (5), 435. doi: CNKI:SUN:ZZXJ.0.2008-05-023

PubMed Abstract | Google Scholar

Wang, X., Xie, G., Wang, X., Zhou, M., Yu, H., Lin, Y., et al. (2015). Urinary metabolite profiling offers potential for differentiation of liver-kidney yin deficiency and dampness-heat internal smoldering syndromes in posthepatitis B cirrhosis patients. Evid Based Complement. Alternat. Med. 2015, 464969. doi: 10.1155/2015/464969

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G., Haw, R. (2017). Functional interaction network construction and analysis for disease discovery. Methods Mol. Biol. 1558, 235–253. doi: 10.1007/978-1-4939-6783-4_11

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J. S., Li, Y. F., Li, Y. Y., Dai, Y., Li, W. K., Zheng, M., et al. (2017). Huangqi decoction alleviates Alpha-Naphthylisothiocyanate induced intrahepatic cholestasis by reversing disordered bile acid and glutathione homeostasis in mice. Front. Pharmacol. 8, 938. doi: 10.3389/fphar.2017.00938

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Zhang, F., Yang, K., Fang, S., Bu, D., Li, H., et al. (2018). SymMap: an integrative database of traditional Chinese medicine enhanced by symptom mapping. Nucleic Acids Res. 47 (D1), D1110–D1117. doi: 10.1093/nar/gky1021

CrossRef Full Text | Google Scholar

Xu, W., Guo, Y. B., Li, X., He, M. R., Liu, S. D. (2016). [Palmitic acid induces hepatocellular oxidative stress and activation of inflammasomes]. Nan Fang Yi Ke Da Xue Xue Bao 36 (5), 655–659. doi: 10.3969/j.issn.1673-4254.2016.05.11

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Fan, W. W., Xu, W., Jiang, S. L., Chen, G. F., Liu, C., et al. (2018). Yiguanjian decoction enhances fetal liver stem/progenitor cell-mediated repair of liver cirrhosis through regulation of macrophage activation state. World J. Gastroenterol. 24 (42), 4759–4772. doi: 10.3748/wjg.v24.i42.4759

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu-Dong, H. U., Ting-Jie, Y. E., Wang, X. L., Bo-Qin, L. I., Liu, P., Biology, D. O. (2014). The prevention and treatment effect of Yin-chen-hao decoction on Concanavalin A-induced acute mice hepatitis by inhibiting the activation of Kupffer cells. Lishizhen Med. Materia Med. Res. 25 (06), 1340-1343. doi: 10.3969/j.issn.1008-0805.2014.06.022

CrossRef Full Text | Google Scholar

Yan, D., Ma, B., Shi, R., Wang, T., Ma, Y. (2015). Involvement of herb-herb interactions in the influences of Radix Scutellaria and Coptis Chinensis on the bioavailability of the anthraquinones form Rhei Rhizoma in rats. Eur. J. Drug Metab. Pharmacokinet. 40 (1), 103–110. doi: 10.1007/s13318-014-0188-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, J., Xie, G., Liang, C., Hu, Y., Zhao, A., Huang, F., et al. (2017). Herbal medicine Yinchenhaotang protects against α-naphthylisothiocyanate-induced cholestasis in rats. Sci. Rep. 7 (1), 4211. doi: 10.1038/s41598-017-04536-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, H., Ye, L., Kang, H., Zhang, D., Tao, L., Tang, K., et al. (2011). HIT: linking herbal active ingredients to targets. Nucleic Acids Res. 39 (Database issue), D1055–D1059. doi: 10.1093/nar/gkq1165

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, J. K., Li, Y. Y., Wang, T. M., Zhong, J., Wu, J. S., Liu, P., et al. (2018a). Simultaneous quantification of multiple components in rat plasma by UPLC-MS/MS and pharmacokinetic study after oral administration of Huangqi decoction. BioMed. Chromatogr. 32 (5), e4178. doi: 10.1002/bmc.4178

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, X., Zhang, P., He, W., Qin, C., Chen, S., Tao, L., et al. (2018b). NPASS: natural product activity and species source database for natural product research, discovery and tool development. Nucleic Acids Res. 46 (D1), D1217–D1222. doi: 10.1093/nar/gkx1026

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. Q., Tang, X. D., Wang, F. Y., Yang, B., Liu, Y. L., Guo, P., et al. (2013). [Effect of ronggan mixture on immunoregulation and hepatocyte apoptosis-related factors in concanavalin A induced acute immunological liver injury mice]. Zhongguo Zhong Xi Yi Jie He Za Zhi 33 (11), 1500–1506. doi: 10.7661/CJIM.2013.11.1500

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G. B., Song, Y. N., Chen, Q. L., Dong, S., Lu, Y. Y., Su, M. Y., et al. (2015). Actions of Huangqi decoction against rat liver fibrosis: a gene expression profiling analysis. Chin. Med. 10, 39. doi: 10.1186/s13020-015-0066-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y. N., Mu, Y. P., Fu, W. W., Ning, B. B., Du, G. L., Chen, J. M., et al. (2015). Yiguanjian decoction and its ingredients inhibit angiogenesis in carbon tetrachloride-induced cirrhosis mice. BMC Complement. Altern. Med. 15, 342. doi: 10.1186/s12906-015-0862-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: network comparative analysis, TCM, chronic liver disease, Yinchenhao Decoction, Huangqi Decoction, Yiguanjian

Citation: Chen Z, Wang X, Li Y, Wang Y, Tang K, Wu D, Zhao W, Ma Y, Liu P and Cao Z (2019) Comparative Network Pharmacology Analysis of Classical TCM Prescriptions for Chronic Liver Disease. Front. Pharmacol. 10:1353. doi: 10.3389/fphar.2019.01353

Received: 19 February 2019; Accepted: 25 October 2019;
Published: 22 November 2019.

Edited by:

Shi-Bing Su, Shanghai University of Traditional Chinese Medicine, China

Reviewed by:

Jianxin Chen, Beijing University of Chinese Medicine, China
Yiyang Hu, Shanghai University of Traditional Chinese Medicine, China
Yonghua Wang, Northwest A&F University, China

Copyright © 2019 Chen, Wang, Li, Wang, Tang, Wu, Zhao, Ma, Liu and Cao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yueming Ma, mayueming@shutcm.edu.cn; Ping Liu, liuliver@vip.sina.com; Zhiwei Cao, zwcao@tongji.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.