Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2019 / Article

Research Article | Open Access

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

Guangyao Li, Yuan Zhou, Daniel Man-Yuen Sze, Chao Liu, Qianru Zhang, Zihao Wang, Hua Yu, Ging Chan, Zhongdao Wu, Shibing Su, Yuanjia Hu, "Active Ingredients and Action Mechanisms of Yi Guan Jian Decoction in Chronic Hepatitis B Patients with Liver Fibrosis", Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 2408126, 13 pages, 2019. https://doi.org/10.1155/2019/2408126

Active Ingredients and Action Mechanisms of Yi Guan Jian Decoction in Chronic Hepatitis B Patients with Liver Fibrosis

Academic Editor: Victor Kuete
Received21 May 2019
Accepted29 Jul 2019
Published03 Sep 2019


Background and Aim. The progression of liver fibrosis in chronic hepatitis B (CHB) patients is currently insufficiently controlled worldwide. The Yi Guan Jian decoction (YGJD) has been widely used in the treatment of liver fibrosis in CHB cases. Although animal studies have reported the antifibrotic effects of the decoction, the active ingredients of the YGJD remain unknown. This study aimed at identifying the potential active ingredients and exploring the mechanisms of action (MOA) of the decoction when treating CHB patients with fibrosis. Methods. Using data mining techniques and a structural clustering analysis, the potential active ingredients were determined. A network analysis of the differentially expressed genes was conducted to identify the potential targets. Selected compounds were docked to the potential targets for the compound-target interaction simulation. In vitro validation, including a cell proliferation assay and Western blot analysis, was conducted to evaluate the prediction results. Results. In the microarray data, 224 differentially expressed genes related to liver fibrosis were considered to be potential targets. Thirty active ingredients of the YGJD and 15 main targets and relevant pathways were identified. Among them, two active ingredients, methylophiopogonone A and 8-geranyloxypsoralen, were validated as exhibiting antifibrotic effects on hepatic stellate cells. Conclusions. We identified the potential active ingredients of the YGJD and proposed the possible explanation for the MOA in the treatment of CHB patients with liver fibrosis. Moreover, this study provides a methodological reference for the systematic investigation of the bioactive compounds and related MOA of a traditional Chinese medicine formula in a clinical context.

1. Introduction

The hepatitis B virus (HBV) is a major health problem worldwide. Indeed, the World Health Organization has estimated that some 257 million people are infected with the HBV globally, with the virus having been associated with nearly one million deaths in 2015 alone, mainly due to complications such as cirrhosis and liver cancer [1]. In the absence of curative treatment, the most common treatments for the HBV include IFN-α and nucleotide analogues, which often lead to severe side effects and drug resistance issues, while having little effect in terms of treating liver fibrosis, cirrhosis, or liver cancer [2].

Due to the issues associated with mainstream treatments, the use of herbal products as a form of complementary medicine has long been widely practiced in Asian countries. In China, over 90% of chronic hepatitis B (CHB) patients received the traditional Chinese medicine (TCM) therapy [3]. Studies have shown various types of TCM to target the HBV [4], liver fibrosis [5], or hepatocellular carcinoma [6]. Further, a meta-analysis of clinical trials found that TCM had an equivalent or better effect than interferon or lamivudine in the treatment of CHB [7]. Moreover, randomized controlled clinical studies have reported the effects of different TCM formulae with regard to improving the clinical parameters of patients with liver fibrosis [8].

As a TCM formula listed in the “Catalogue of Classical TCM Prescriptions” issued by the State Administration of TCM of the People’s Republic of China, the Yi Guan Jian decoction (YGJD) is known to have hepatoprotective effects [9]. The YGJD consists of six herbal medicines, namely, Rehmanniae Radix (Rehmannia glutinosa), Glehniae Radix (Glehnia littoralis), Ophiopogonis Radix (Ophiopogon japonicus), Angelica Sinensis Radix (Angelica sinensis), Lycii Fructus (Lycium barbarum), and Toosendan Fructus (Melia toosendan). Three possible mechanisms of action (MOA) have been suggested for the antifibrotic effects of the YGJD: first, the inhibition of hepatic stellate cell (HSC) activation and hepatocytes apoptosis [10]; second, antiangiogenesis via the HIF-1α/VEGF signaling pathway [11]; and finally, anti-inflammation [12]. However, the bioactive ingredients of the YGJD remain unknown, while the MOA at the molecular level require further investigation.

In clinical practice involving TCM, the YGJD has been widely used as a basic formula in the treatment of liver-kidney Yin deficiency syndrome (LKYDS) in patients with CHB [13]. LKYDS is mainly seen in CHB patients with severe liver fibrosis [14]. We hypothesized that the beneficial clinical effect of the YGJD in terms of reducing liver fibrosis stems from the oral bioavailable chemical compounds found in the decoction, which are related to the changes in the gene transcription of the blood leukocytes seen in CHB patients with LKYDS. This study hence had two key aims: first, to identify the potential active ingredients of the YGJD with regard to reducing liver fibrosis, and second, to delineate the MOA of the YGJD in relation to the effective treatment of CHB patients with liver fibrosis.

2. Materials and Methods

To achieve the abovementioned research aims, this study gathered relevant bioinformatic data concerning both the YGJD and CHB patients with LKYDS. Further, the study comprehensively employed various methods of analysis, including network pharmacology [15], microarray analysis, molecular docking, and in vitro experiments, to analyze the collected data. The research roadmap is shown in Figure 1, while detailed descriptions of the utilized materials and methods are given below.

2.1. Compound Collection and Oral Bioavailability Prediction

The relevant chemical information concerning the YGJD was retrieved from four TCM databases in December 2017, including the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP) [16], the Herbal Ingredients Targets Database (HIT) [17], the Psoriasis Database of Traditional Chinese Medicine (PDTCM) [18], and the Potential Target Database of Traditional Chinese Medicine (TCM-PTD). The three-dimensional structures were downloaded from the corresponding database or PubChem.

We used the two parameters proposed by Veber for the preliminary screening of the collected compounds. Veber’s filter was used for the oral bioavailability prediction [19], which meant that the number of rotatable bonds of a given compound had to be ≤10, while the topological polar surface area had to be ≤140 Å2 [20]. Further, DataWarrior was employed to calculate the molecular properties of the compounds.

2.2. Structural Clustering Analysis

The ChemmineR package run in Rstudio was used to compute the fingerprints and the Tanimoto coefficients among the compounds. The obtained cross-similarity matrix was subsequently applied to build the structural similarity network of the compounds with Tanimoto coefficients ≥0.6 [21]. The MCODE plugin in Cytoscape was used to partition the clusters for the network.

2.3. Microarray Data Analysis

The microarray data concerning the leukocyte samples obtained from three CHB patients with LKYDS and three healthy donors were provided by the Research Center for TCM Complexity System, Shanghai University of TCM [22]. The differentially expressed genes were identified under the conditions of a fold-change >2 and a value <0.05. The collected genes were then submitted to the STRING database to establish the protein-protein interaction (PPI) network, while the edges with confidence scores greater than 0.7 were retained for further analysis.

The liver fibrosis-related genes were collected from DisGeNET. Using the search term of “Fibrosis, Liver” (UMLS CUI: C0239946), a total of 354 genes were selected from among the results. The CHB with LKYDS genes that exhibited high correlations (a confidence score greater than 0.9 in the STRING database) with the liver fibrosis-related genes were considered to be potential targets for the YGJD. PPI network analysis via pathway enrichment and topological parameter calculation was performed to verify the importance of the potential targets [23].

2.4. Molecular Docking

A cross-docking analysis involving sets of compounds and targets was conducted using AutoDock Vina. The crystal structures of the proteins were downloaded from the Research Collaboratory for Structural Bioinformatics Protein Data Bank, and only those X-ray structures with a resolution of less than 3 Å were selected for the analysis (Table S1). The binding sites were obtained from the original ligand binding locations of the crystal structures, or the reported pockets, or else were predicted by the CASTp 3.0 web server [24]. The gird boxes were adjusted to cover the entire pocket, the majority of which were 22 × 22 × 22 Å3. To eliminate any bias caused by the properties of the protein pockets, all the docking scores for each protein were normalized by the median number, as shown in equation (1) [25]. The receiver operating characteristic (ROC) curves between the normalized scores of the active ligands and the ingredients in the YGJD were used to determine the cutoff values [26]:where is the original docking score and is the normalized score, i and j are the indexes of the corresponding ligand and protein, respectively, and and are the median and standard deviation, respectively, of the docking scores of protein j.

2.5. C-T and C-T-P Network Analysis

Cytoscape was used to visualize the compound-target (C-T) network and the compound-target-pathway (C-T-P) network for the YGJD in the treatment of fibrosis caused by CHB. Two main parameters, namely, the degree centrality (DC) and the betweenness centrality (BC), were used to determine the key nodes in the network [27]. The pathway enrichment of the targets was performed using the STRING database with a false discover rate (FDR) <0.05. In the C-T-P network, the number of shortest paths from all the active ingredients to a biological pathway was exploited to rank the pathways, and the utilized calculation method was as shown in the following equation [28]:where is the number of shortest paths from the active ingredients to pathway m, is the in-DC of the target n in the C-T-P network, N represents the number of targets in the network, and is a dummy variable. Further, when there is a path between the target m and the pathway n; otherwise, .

2.6. Cell Proliferation Assay

The LX-2 immortalized human hepatic stellate cell line was purchased from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China) [29]. The cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, USA), which was supplemented with 10% fetal bovine serum (FBS, Gibco, USA), 100 U/mL penicillin, and 100 U/mL streptomycin and then incubated at 37°C in humidified air containing 5% CO2.

After being seeded overnight in 96-well plates with DMEM containing 2% FBS, the LX-2 cells were treated with serially diluted compounds for 24 h, 48 h, and 72 h. The cell number was measured by means of a Cell Counting Kit-8 (Dojindo, Kumamoto, Japan) assay. The absorbance was monitored at 450 nm using a microplate reader (PE Envision). Further, the viability percentage was calculated using the following equation:where T and C refer to the absorbance of the treatment group and the control group, respectively.

2.7. Western Blot Analysis

The cellular content of the collagen I and the α-smooth muscle actin (α-SMA) in the LX-2 cells was determined by means of Western blot analysis. The LX-2 cells were treated under different conditions with the active ingredients in 5 ng/mL of TGF-β1 for 24 h. The lysates of the denatured cells were separated on NuPAGE 4–12% Bis-Tris gel (Thermo Fisher Scientific), and the proteins were transferred to Hybond ECL membranes (GE Healthcare). The membranes were incubated overnight at 4°C with the primary antibodies against collagen type I (Abcam), α-SMA (Abcam), and β-actin (Sigma) as the endogenous controls, followed by the corresponding secondary peroxidase-coupled antibodies (Sigma). After the chemiluminescence had been enhanced (ECL, Pierce), the digital detection was evaluated using Image Lab Software (ChemiDoc, Bio-Rad).

3. Results

3.1. Identification of the Representative Compounds of the YGJD

A total of 836 compounds contained within the YGJD were acquired from the databases, with 526 being predicted to exhibit better oral bioavailability. A structural similarity network was constructed to identify the representative compounds among the 526 ingredients (Figure 2). The network was divided into 53 clusters. Further investigation of the network showed that the compounds contained within clusters 1–3 were mainly cholestanes from Lycii Fructus, while those within clusters 6–12 were mostly limonoids from Toosendan Fructus, and those within clusters 13–15 were mostly isoflavones from Ophiopogonis Radix.

According to the theory of structure-activity relationship, structurally similar compounds tend to perform similar biological activities [30]. The representative compounds were identified after ranking all the compounds based on their DC and BC scores. The compounds with the top scores in each cluster were selected, giving a total of 57 representative compounds. The representative compounds, in addition to the compounds shown as isolated nodes within the network, were chosen for further docking analysis. In total, 278 compounds (Table S2) were docked to the selected targets.

3.2. The Pool of Potential Targets for the YGJD

The CHB with LKYDS genes were mapped by means of analyzing the microarray data concerning the CHB with LKYDS patients who had been diagnosed by senior TCM doctors at the Shanghai Longhua Hospital. A PPI network for the differently expressed genes was constructed (Figure 3(a)). In order to obtain the genes highly associated with liver fibrosis and to reduce the potential bias of the microarray data, liver fibrosis-related genes were collected from the DisGeNET database. Some 224 CHB with LKYDS genes that exhibited a high correlation with liver fibrosis were selected as potential targets.

To verify the importance of the selected targets, the different gene sets were subjected to a pathway enrichment analysis (Figure 3(b)) and a topological parameter analysis (Figure 3(c)). After manually deleting those pathways that were apparently unrelated diseases, the potential targets were found to be involved in 64 pathways, while all 740 genes were involved in ten pathways, and those genes not selected as potential targets were only correlated with one pathway. Further, it was found that the selected targets exhibited significantly higher DC and BC scores than the unselected ones. Taken together, the network analysis results indicated the successful exclusion of interfering genes as well as the high correlation between LKYDS and liver fibrosis in CHB patients.

3.3. Compound-Target Interaction Simulation Using In Silico Molecular Docking

The selected compounds from Figure 2 were docked with the potential targets one to one with AutoDock Vina. The active ligands collected from the protein crystal structures, the DrugBank, or the Therapeutic Target Database were used to determine the cutoff values for the docking scores. The distribution of the normalized docking scores of the active ligands (Table S3), as well as the selected compounds from the YGJD, is illustrated in Figure 4(a).

The ROC curve presented in Figure 4(b) depicts the sensitivity/specificity trade-off for the threshold values ranging from −4 to 1.5. A cutoff value of −1.75 was selected, with a sensitivity of 0.722 and a specificity of 0.952, and 4.79% of the compound-target interactions were determined to fall within that range. The C-T network was subsequently constructed and visualized based on the predicted compound-target interactions.

3.4. Identification of the Active Ingredients in the YGJD That Act on CHB-Related Liver Fibrosis

The C-T network consisted of selected compounds and their respective targets, as visualized in Figure 5. The compounds that exhibited the top 30 out-DC scores were considered to be the active ingredients, and they are listed in Table 1.

IDOut DCCompound nameOriginMolecular formulaCluster

YGJ51065DiosgeninOphiopogonis RadixC27H42O33
YGJ00960BelladonnineLycii FructusC34H42N2O4Isolated
YGJ386596-Deacetyloxy-7-deacetylchisochetonToosendan FructusC26H36O4Isolated
YGJ7755524-Methylenecycloartanol ferulateLycii FructusC41H60O4Isolated
YGJ06154Methylophiopogonanone AOphiopogonis RadixC19H16O615
YGJ23050Mesendanin MToosendan FructusC30H44O46
YGJ83346MarmesininGlehniae RadixC20H24O9Isolated
YGJ26044Meliasenin XToosendan FructusC30H48O57
YGJ303433,3′Z-6.7′,7.6′-diligustilideAngelica Sinensis RadixC24H28O4Isolated
YGJ20542Mesendanin JToosendan FructusC28H40O848
YGJ37742Trichilinin EToosendan FructusC35H42O812
YGJ81241Lantadene ALycii FructusC35H54O5Isolated
YGJ01940Ophiopogonone AOphiopogonis RadixC18H14O614
YGJ48539LupeolGlehniae Radix;
Lycii Fructus
YGJ05139N-p-coumaroyltyramineOphiopogonis RadixC18H15NO4Isolated
YGJ25338Meliasenin WToosendan FructusC30H50O49
YGJ16038Meliasenin PToosendan FructusC31H50O48
YGJ32035Nimbolinin CToosendan FructusC38H46O911
YGJ60232Levistolide AAngelica Sinensis RadixC24H28O430
YGJ23531Mesendanin UToosendan FructusC32H50O6Isolated
YGJ012308-GeranyloxypsoralenGlehniae RadixC21H22O435
YGJ62428Campesteryl ferulateLycii FructusC38 H56 O453
YGJ5662224-MethylenecycloartanolLycii FructusC31H52O5
YGJ6902112-O-nicotinoylisolineoloneLycii FructusC27H35NO6Isolated
YGJ72320IsotetandrineAngelica Sinensis RadixC38H42N2O6Isolated
YGJ48618SitoglusideRehmanniae Radix;
Glehniae Radix;
Angelica Sinensis Radix
YGJ41418Toosendanic acid AToosendan FructusC30H48O46
YGJ005172′-HydroxymatteucinolOphiopogonis RadixC18H18O6Isolated
YGJ00616Taurochenodeoxycholic acidLycii FructusC26H45NO6SIsolated
YGJ358162-Hydroxy-3-(3-methylbut-2-enyl) furo[3,2 g] chromen-7-oneGlehniae RadixC16H14O4Isolated

3.5. MOA for the YGJD When Treating CHB-Related Liver Fibrosis

The C-T-P network was constructed (Figure 6) to further explore the MOA of the YGJD in the treatment of liver fibrosis caused by CHB. For each active ingredient, the targets that exhibited the top ten highest binding affinities to the ingredients were retained and shown in the graph. The links between the targets and the pathways were drawn based on the results obtained in the pathway enrichment analysis described above.

The BC was used to measure the importance of the targets within the C-T-P network. Those targets with high BC scores were the intermediaries that controlled the flow of interactions within the network. Based on the BC score, the top 15 targets, which are listed in Table 2, were screened as the main targets of YGJD treatment.

SymbolEntrez IDProtein nameBC (10−3)

RELA5970Transcription factor p65, NF-κB subunit5.28
EGFR1956Epidermal growth factor receptor2.66
RAC35881Ras-related C3 botulinum toxin substrate 31.44
FOS2353Proto-oncogene c-Fos, AP-1 transcription factor subunit1.25
IFNA13439Interferon-alpha 1/131.05
JUN3725Transcription factor AP-10.89
MAPK126300Mitogen-activated protein kinase 120.61
INPP5B3633Type II inositol 1,4,5-trisphosphate 5-phosphatase0.56
TGFB37043Transforming growth factor-beta 30.47
RXRG6258Retinoic acid receptor RXR-gamma0.45
CYP2C91559Cytochrome P450 2C90.43
OCRL4952Inositol polyphosphate 5-phosphatase OCRL-10.36
HGS9146Hepatocyte growth factor-regulated tyrosine kinase substrate0.32

Next, based on the C-T-P network, we calculated the shortest paths from the active ingredients to each pathway. By ranking the scores for this indicator, the key biological pathways influenced by the YGJD were identified. The top ten pathways included MAPK signaling pathway, cytokine-cytokine receptor interaction, Toll-like receptor signaling pathway, hepatitis C, endocytosis, PI3K-Akt signaling pathway, hepatitis B, metabolic pathways, chemokine signaling pathway, and osteoclast differentiation. Details concerning these key pathways are provided in Table S4.

3.6. Antifibrotic Effects of Methylophiopogonone A and 8-Geranyloxypsoralen on LX-2 Cells

In the 30 active ingredients in Table 1, the 17 representative compounds were selected based on their topological features in Figure 2. The 13 compounds in 17 were considered for the further experimental validation, except for 4 compounds reported to show antifibrotic or hepatoprotective effects in existing literature. In consideration of the availability of these sample compounds, 4 representative compounds in Table 1 were selected for further pharmacological experiments. As shown in Figure 7(a), two ingredients in four, namely, methylophiopogonone A from cluster 15 and 8-geranyloxypsoralen from cluster 35, were found to inhibit the proliferation of the LX-2 cells in a time- and dose-dependent manner. However, there was inadequate evidence to confirm the role of sitogluside and lupeol in inhibiting the proliferation of the LX-2 cells, which indicates the need for further validation of the other potential MOA of the YGJD predicted using our method.

Additionally, two of the main markers of HSC activation, namely, collagen I and α-SMA, were determined by means of Western blot analysis after the cells were treated with the compounds for 24 h. As shown in Figure 7(b), the methylophiopogonone A and 8-geranyloxypsoralen inhibited the TGF-β1-induced collagen I and α-SMA expression.

4. Discussion

The YGJD is an herbal formula that is widely used in clinical practice involving TCM for the treatment of CHB patients with LKYDS, a syndrome known to be related to liver fibrosis. Modern pharmacological studies have confirmed the antifibrotic and anti-inflammatory effects of YGJD [31]. In this study, we explored the effects of YGJD on liver fibrosis of CHB by using the LKYDS as a clue to collect microarray data in the specific patient pool of CHB. The microarray data analysis further provided a potential target set in the treatment of liver fibrosis with YGJD from a clinical perspective. The interactions between the compounds of YGJD and the potential targets were revealed by means of molecular docking. Moreover, the method of network pharmacology was used to predict the bioactive ingredients and biological basis of YGJD acting on CHB patients with liver fibrosis.

For the 30 active ingredients predicted by our method, some have been reported to exhibit effects against liver fibrosis. Diosgenin from Ophiopogonis Radix and levistolide A from Angelica Sinensis Radix were observed to inhibit the HSC activation and proliferation [32, 33]. Lupeol, lantadene A, and two trans-ferulic acids from Lycii Fructus, namely, 24-methylenecycloartanol ferulate and campesteryl ferulate, were reported to exhibit hepatoprotective effects in mice with toxically induced liver injuries [3436]. N-p-coumaroyltyramine from Ophiopogonis Radix, 8-geranyloxypsoralen from Glehniae Radix, and sitogluside from Rehmanniae Radix were all reported to be involved in anti-inflammatory activities [3739]. In addition, several limonoids from Toosendan Fructus, such as mesendanin M, trichilinin E, and nimbolinin C, exhibited high out-DC scores in the C-T network. It has previously been reported that limonoids extracted from Meliaceae possess hepatoprotective functions and could hence alleviate liver fibrosis by reducing the TGF-β1 and collagen levels [40, 41]. Thus, the limonoids listed in Table 1 might be worthy of further investigation.

Interestingly, 17 of the 30 active ingredients were found to belong to 57 representative compounds in Figure 2, whereas only 13 of the 221 isolated nodes performed well in the docking analysis, which indicated that the highly connected nodes in the compound similarity network had greater opportunity to be biologically active. The compounds with similar structures might represent the driving force within the YGJD, since they could work synergistically to achieve the identified therapeutic effects.

To reveal the MOA of the YGJD in relation to the effective treatment of CHB patients with liver fibrosis, the main targets and the relevant pathways were obtained using the C-T-P network analysis. The findings reported in the literature regarding these principal targets support the suggestion that the YGJD attenuated or even reversed the progression of liver fibrosis among CHB patients using three main approaches. The first approach was the inhibition of the fibrosis-related HSC activation and proliferation. RELA (also known as nuclear factor-κB p65), which has previously been reported to be downregulated in YGJD treatment [42], was found to be a significant factor against the apoptosis of the HSC [43]. Additionally, it has been reported that the inhibition of EGFR could decrease the activation of the HSC [44]. Further, OCRL and INPP5B could regulate HSC proliferation by modulating the AKT signaling pathway [45]. HGS might affect HSC activation via the modulation of both Smad2 and PDGFR [46]. The second approach concerned the fact that some of the main targets were involved in preventing excessive extracellular matrix (ECM) deposition. Fibronectin was part of the ECM [5]. JUN and FOS were reported to influence several genes related to the ECM in the HSC [47]. Moreover, MAPK12 might be involved in regulating α1(I) collagen, a key part of the ECM in fibrosis, via the p38 MAPK signaling pathway [48]. The third approach involved regulating the inflammatory response to the HBV by regulating both IL-1B and IFNA1. It has previously been reported that the YGJD could decrease the levels of IL-1B so as to exert anti-inflammatory effects [12]. Further, interferon-alpha encoded by IFNA1 can inhibit the replication of the HBV, thereby leading to the remission of liver diseases [49]. All of the three summarized approaches are consistent with the antifibrotic effects of YGJD reported in the existing literature. This further confirms the accuracy of the target prediction in this study.

The key pathways were also confirmed to be involved in the above-mentioned three approaches accordingly. The MAPK signaling pathway and the PI3K-Akt signaling pathway both played an important role in the entire pathological process of liver fibrosis. Moreover, it was found that the active substances targeting the two pathways could hinder the HSC activation and exert an antifibrotic effect [50, 51]. Cytokine-cytokine receptor interaction, the chemokine signaling pathway, and the Toll-like receptor signaling pathway were all highly related to the immunological process associated with the hepatitis B infection [52]. As for the metabolic pathways, both the retinol metabolism and fatty acid metabolism pathways shown in the C-T-P network were found to be associated with HSC activities as well as with antifibrotic effects [53, 54].

However, a research limitation should be noted. In this study, only two active ingredients were verified to exhibit antifibrotic effects on the HSC by significantly inhibiting the cell proliferation. Further animal experiments and clinical trials on the two ingredients and more potential active ingredients may be performed to explore the scientific nature of YGJD.

5. Conclusions

This study predicted 224 potential targets for treating CHB-caused liver fibrosis using microarray data analysis and PPI network analysis and also identified 278 candidates of active ingredients of YGJD through network cluster analysis. The interaction between the compounds and potential targets was revealed by molecular docking analysis and led to a total of 30 active ingredients. Among these active ingredients, methylophiopogonone A and 8-geranyloxypsoralen were found to inhibit HSC proliferation in LX-2 cells in vitro. And almost half of the active ingredients have been reported to have antifibrotic or anti-inflammatory effects.

As for the MOA, fifteen genes related to the liver fibrosis caused by the CHB were predicted to be the main targets of YGJD. Literature findings regarding these main targets support YGJD attenuated or even reversed the progression of liver fibrosis among CHB patients by inhibiting fibrosis-related HSC activation and proliferation via RELA, EGFR, OCRL, INPP5B, and HGS, preventing excessive ECM deposition via FN1, JUN, FOS, and MAPK12, and regulating inflammatory response to HBV by regulating both IL-1B and IFNA1. Furthermore, the network analysis revealed the pathways involved in the YGJD treatment, including MAPK signaling, PI3K-Akt signaling, cytokine-cytokine receptor interaction, Toll-like receptor signaling, endocytosis, and metabolism pathway.

In summary, we identified the active ingredients in the YGJD and proposed the possible MOA for the YGJD in the treatment of liver fibrosis among CHB patients. Our study not only provided a direction for the development of new drugs to treat liver fibrosis in CHB patients, but also offered an alternative approach for the systematic investigation of the MOA of a TCM formula by using the YGJD as an example.

Data Availability

The data that support the findings of this study are available in the Supplementary Materials and openly available in the TCMSP (http://lsp.nwu.edu.cn/tcmsp.php), the HIT (http://lifecenter.sgst.cn/hit/), the PDTCM (http://cadd.gdhtcm.com:2180/PDTCM/), the TCMPTD (http://pharminfo.zju.edu.cn/ptd), the PubChem (https://pubchem.ncbi.nlm.nih.gov/), and the DisGeNET (http://www.disgenet.org/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Guangyao Li and Yuan Zhou contributed equally to this work.


We thank the Science and Technology Development Fund of Macao SAR and the University of Macau and Guangdong-Macao Traditional Chinese Medicine Technology Industrial Park for financial support for the projects 013-2015-A1, MYRG2016-00144-ICMS-QRCM and D-Pro-02742017 for this research.

Supplementary Materials

Table S1: potential targets for docking analysis. Table S2: selected compounds for docking analysis. Table S3: active ligands as reference for the docking analysis. Table S4: key pathways involved with the YGJD acting on CHB-related liver fibrosis. Table S5: raw data of the docking results. (Supplementary Materials)


  1. WHO, WHO Guidelines Approved by the Guidelines Review Committee, World Health Organization, Geneva, Switzerland, 2015.
  2. N. A. Terrault, N. H. Bzowej, K.-M. Chang, J. P. Hwang, M. M. Jonas, and M. H. Murad, “AASLD guidelines for treatment of chronic hepatitis B,” Hepatology, vol. 63, no. 1, pp. 261–283, 2016. View at: Publisher Site | Google Scholar
  3. X. Zheng, J. Wang, and D. Yang, “Antiviral therapy for chronic hepatitis B in China,” Medical Microbiology and Immunology, vol. 204, no. 1, pp. 115–120, 2015. View at: Publisher Site | Google Scholar
  4. X. Cui, Y. Wang, N. Kokudo, D. Fang, and W. Tang, “Traditional Chinese medicine and related active compounds against hepatitis B virus infection,” BioScience Trends, vol. 4, no. 2, pp. 39–47, 2010. View at: Google Scholar
  5. R. Bataller and D. A. Brenner, “Liver fibrosis,” Journal of Clinical Investigation, vol. 115, no. 2, pp. 209–218, 2005. View at: Publisher Site | Google Scholar
  6. B. Hu, S. S. Wang, and Q. Du, “Traditional Chinese medicine for prevention and treatment of hepatocarcinoma: from bench to bedside,” World Journal of Hepatology, vol. 7, no. 9, pp. 1209–1232, 2015. View at: Publisher Site | Google Scholar
  7. L. Zhang, G. Wang, W. Hou, P. Li, A. Dulin, and H. L. Bonkovsky, “Contemporary clinical research of traditional Chinese medicines for chronic hepatitis B in China: an analytical review,” Hepatology, vol. 51, no. 2, pp. 690–698, 2010. View at: Publisher Site | Google Scholar
  8. L. Zhang and D. Schuppan, “Traditional Chinese medicine (TCM) for fibrotic liver disease: hope and hype,” Journal of Hepatology, vol. 61, no. 1, pp. 166–168, 2014. View at: Publisher Site | Google Scholar
  9. R. Weiskirchen, A. Mahli, S. Weiskirchen, and C. Hellerbrand, “The hop constituent xanthohumol exhibits hepatoprotective effects and inhibits the activation of hepatic stellate cells at different levels,” Frontiers in Physiology, vol. 6, p. 303, 2015. View at: Publisher Site | Google Scholar
  10. Y. Mu, P. Liu, G. Du et al., “Action mechanism of Yi Guan Jian decoction on CCl4 induced cirrhosis in rats,” Journal of Ethnopharmacology, vol. 121, no. 1, pp. 35–42, 2009. View at: Publisher Site | Google Scholar
  11. Y. N. Zhou, Y. P. Mu, W. W. Fu et al., “Yiguanjian decoction and its ingredients inhibit angiogenesis in carbon tetrachloride-induced cirrhosis mice,” BMC Complementary and Alternative Medicine, vol. 15, no. 1, p. 342, 2015. View at: Publisher Site | Google Scholar
  12. S. Shui, S. Shen, R. Huang, B. Xiao, and J. Yang, “Metabonomic analysis of biochemical changes in the plasma and urine of carrageenan-induced rats after treatment with Yi-Guan-Jian decoction,” Journal of Chromatography B, vol. 1033-1034, pp. 80–90, 2016. View at: Publisher Site | Google Scholar
  13. W. Jia, W.-Y. Gao, Y.-Q. Yan et al., “The rediscovery of ancient Chinese herbal formulas,” Phytotherapy Research, vol. 18, no. 8, pp. 681–686, 2004. View at: Publisher Site | Google Scholar
  14. X.-X. Zeng, Z.-X. Bian, T.-X. Wu, S.-F. Fu, E. Ziea, and W. T. C. Woon, “Traditional Chinese medicine syndrome distribution in chronic hepatitis B populations: a systematic review,” The American Journal of Chinese Medicine, vol. 39, no. 6, pp. 1061–1074, 2011. View at: Publisher Site | Google Scholar
  15. S. Li, “Mapping ancient remedies: applying a network approach to traditional Chinese medicine,” Science, vol. 350, no. 6262, pp. S72–S74, 2015. View at: Google Scholar
  16. J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014. View at: Publisher Site | Google Scholar
  17. H. Ye, L. Ye, H. Kang et al., “HIT: linking herbal active ingredients to targets,” Nucleic Acids Research, vol. 39, pp. D1055–D1059, 2011. View at: Publisher Site | Google Scholar
  18. D. Wang, J. Gu, W. Zhu et al., “PDTCM: a systems pharmacology platform of traditional Chinese medicine for psoriasis,” Annals of Medicine, vol. 49, no. 8, pp. 652–660, 2017. View at: Publisher Site | Google Scholar
  19. C. A. Lipinski, “Lead- and drug-like compounds: the rule-of-five revolution,” Drug Discovery Today: Technologies, vol. 1, no. 4, pp. 337–341, 2004. View at: Publisher Site | Google Scholar
  20. D. F. Veber, S. R. Johnson, H.-Y. Cheng, B. R. Smith, K. W. Ward, and K. D. Kopple, “Molecular properties that influence the oral bioavailability of drug candidates,” Journal of Medicinal Chemistry, vol. 45, no. 12, pp. 2615–2623, 2002. View at: Publisher Site | Google Scholar
  21. J. S. Delaney, “Assessing the ability of chemical similarity measures to discriminate between active and inactive compounds,” Molecular Diversity, vol. 1, no. 4, pp. 217–222, 1996. View at: Publisher Site | Google Scholar
  22. Y. Y. Lu, Q. L. Chen, Y. Guan et al., “Study of ZHENG differentiation in hepatitis B-caused cirrhosis: a transcriptional profiling analysis,” BMC Complementary and Alternative Medicine, vol. 14, no. 1, p. 371, 2014. View at: Publisher Site | Google Scholar
  23. X. He and J. Zhang, “Why do hubs tend to be essential in protein networks?” PLoS Genetics, vol. 2, no. 6, p. e88, 2006. View at: Publisher Site | Google Scholar
  24. W. Tian, C. Chen, X. Lei, J. Zhao, and J. Liang, “CASTp 3.0: computed atlas of surface topography of proteins,” Nucleic Acids Research, vol. 46, no. W1, pp. W363–w367, 2018. View at: Publisher Site | Google Scholar
  25. Q. Luo, L. Zhao, J. Hu, H. Jin, Z. Liu, and L. Zhang, “The scoring bias in reverse docking and the score normalization strategy to improve success rate of target fishing,” PLoS One, vol. 12, no. 2, Article ID e0171433, 2017. View at: Publisher Site | Google Scholar
  26. M. W. Chang, W. Lindstrom, A. J. Olson, and R. K. Belew, “Analysis of HIV wild-type and mutant structures via in silico docking against diverse ligand libraries,” Journal of Chemical Information and Modeling, vol. 47, no. 3, pp. 1258–1262, 2007. View at: Publisher Site | Google Scholar
  27. M. E. J. Newman, “The structure of scientific collaboration networks,” Proceedings of the National Academy of Sciences, vol. 98, no. 2, pp. 404–409, 2001. View at: Publisher Site | Google Scholar
  28. Z. Xiang, H. Sun, X. Cai, D. Chen, and X. Zheng, “The study on the material basis and the mechanism for anti-renal interstitial fibrosis efficacy of rhubarb through integration of metabonomics and network pharmacology,” Molecular BioSystems, vol. 11, no. 4, pp. 1067–1078, 2015. View at: Publisher Site | Google Scholar
  29. L. Xu, A. Y. Hui, E. Albanis et al., “Human hepatic stellate cell lines, LX-1 and LX-2: new tools for analysis of hepatic fibrosis,” Gut, vol. 54, no. 1, pp. 142–151, 2005. View at: Publisher Site | Google Scholar
  30. P. G. Mezey, “Functional groups in quantum chemistry,” Advances in Quantum Chemistry, vol. 27, pp. 163–222, 1996. View at: Publisher Site | Google Scholar
  31. W. L. Yao, W. W. Fan, X. L. Chen et al., “Researching prospect of correspondence between prescription and syndromes on the pathologic basis of yin deficiency syndrome of liver and kidney and effective mechanism of decoction of Yiguan Jian,” Journal of Developing Drugs, vol. 06, no. 02, pp. 1–6, 2017. View at: Publisher Site | Google Scholar
  32. W. L. Xie, R. Jiang, X. L. Shen, Z. Y. Chen, and X. M. Deng, “Diosgenin attenuates hepatic stellate cell activation through transforming growth factor-beta/Smad signaling pathway,” International Journal of Clinical and Experimental Medicine, vol. 8, no. 11, pp. 20323–20329, 2015. View at: Google Scholar
  33. T.-F. Lee, Y.-L. Lin, and Y.-T. Huang, “Studies on antiproliferative effects of phthalides from Ligusticum chuanxiong in hepatic stellate cells,” Planta Medica, vol. 73, no. 6, pp. 527–534, 2007. View at: Publisher Site | Google Scholar
  34. S. P. Preetha, M. Kanniappan, E. Selvakumar, M. Nagaraj, and P. Varalakshmi, “Lupeol ameliorates aflatoxin B1-induced peroxidative hepatic damage in rats,” Comparative Biochemistry and Physiology Part C: Toxicology & Pharmacology, vol. 143, no. 3, pp. 333–339, 2006. View at: Publisher Site | Google Scholar
  35. C. Grace-Lynn, Y. Chen, L. Latha et al., “Evaluation of the hepatoprotective effects of Lantadene A, a pentacyclic triterpenoid of Lantana plants against acetaminophen-induced liver damage,” Molecules, vol. 17, no. 12, pp. 13937–13947, 2012. View at: Google Scholar
  36. C. Chotimarkorn and H. Ushio, “The effect of trans-ferulic acid and gamma-oryzanol on ethanol-induced liver injury in C57BL mouse,” Phytomedicine, vol. 15, no. 11, pp. 951–958, 2008. View at: Publisher Site | Google Scholar
  37. L.-C. Lin, C.-C. Shen, Y.-C. Shen, and T.-H. Tsai, “Anti-inflammatory neolignans from Piper kadsura,” Journal of Natural Products, vol. 69, no. 5, pp. 842–844, 2006. View at: Publisher Site | Google Scholar
  38. D.-Y. Shen, Y.-Y. Chan, T.-L. Hwang et al., “Constituents of the roots of Clausena lansium and their potential anti-inflammatory activity,” Journal of Natural Products, vol. 77, no. 5, pp. 1215–1223, 2014. View at: Publisher Site | Google Scholar
  39. Y. Ding, C. Liang, J. H. Kim et al., “Triterpene compounds isolated from Acer mandshuricum and their anti-inflammatory activity,” Bioorganic & Medicinal Chemistry Letters, vol. 20, no. 5, pp. 1528–1531, 2010. View at: Publisher Site | Google Scholar
  40. M. P. Germanò, V. D’Angelo, R. Sanogo, A. Morabito, S. Pergolizzi, and R. De Pasquale, “Hepatoprotective activity of Trichilia roka on carbon tetrachloride-induced liver damage in rats,” Journal of Pharmacy and Pharmacology, vol. 53, no. 11, pp. 1569–1574, 2001. View at: Publisher Site | Google Scholar
  41. S. Fan, H.-N. Chen, C.-J. Wang, W.-C. Tseng, H.-K. Hsu, and C.-F. Weng, “Toona sinensis Roem (Meliaceae) leaf extract alleviates liver fibrosis via reducing TGFβ1 and collagen,” Food and Chemical Toxicology, vol. 45, no. 11, pp. 2228–2236, 2007. View at: Publisher Site | Google Scholar
  42. Y. Xiang, B.-Y. Pang, Y. Zhang et al., “Effect of Yi Guan Jian decoction on differentiation of bone marrow mesenchymalstem cells into hepatocyte-like cells in dimethylnitrosamine-induced liver cirrhosis in mice,” Molecular Medicine Reports, vol. 15, no. 2, pp. 613–626, 2017. View at: Publisher Site | Google Scholar
  43. M. Pinzani and K. Rombouts, “Liver fibrosis: from the bench to clinical targets,” Digestive and Liver Disease, vol. 36, no. 4, pp. 231–242, 2004. View at: Publisher Site | Google Scholar
  44. B. C. Fuchs, Y. Hoshida, T. Fujii et al., “Epidermal growth factor receptor inhibition attenuates liver fibrosis and development of hepatocellular carcinoma,” Hepatology, vol. 59, no. 4, pp. 1577–1590, 2014. View at: Publisher Site | Google Scholar
  45. M. Bohdanowicz, D. M. Balkin, P. De Camilli, and S. Grinstein, “Recruitment of OCRL and Inpp5B to phagosomes by Rab5 and APPL1 depletes phosphoinositides and attenuates Akt signaling,” Molecular Biology of the Cell, vol. 23, no. 1, pp. 176–187, 2012. View at: Publisher Site | Google Scholar
  46. A. A. Farooqi, S. Waseem, A. M. Riaz et al., “PDGF: the nuts and bolts of signalling toolbox,” Tumor Biology, vol. 32, no. 6, pp. 1057–1070, 2011. View at: Publisher Site | Google Scholar
  47. J. E. Poulos, J. D. Weber, J. M. Bellezzo et al., “Fibronectin and cytokines increase JNK, ERK, AP-1 activity, and transin gene expression in rat hepatic stellate cells,” American Journal of Physiology-Gastrointestinal and Liver Physiology, vol. 273, no. 4, pp. G804–G811, 1997. View at: Publisher Site | Google Scholar
  48. C. J. Parsons, M. Takashima, and R. A. Rippe, “Molecular mechanisms of hepatic fibrogenesis,” Journal of Gastroenterology and Hepatology, vol. 22, no. s1, pp. S79–S84, 2007. View at: Publisher Site | Google Scholar
  49. A. S. F. Lok and B. J. McMahon, “Chronic hepatitis B,” Hepatology, vol. 45, no. 2, pp. 507–539, 2007. View at: Publisher Site | Google Scholar
  50. W. He, F. Shi, Z. W. Zhou et al., “A bioinformatic and mechanistic study elicits the antifibrotic effect of ursolic acid through the attenuation of oxidative stress with the involvement of ERK, PI3K/Akt, and p38 MAPK signaling pathways in human hepatic stellate cells and rat liver,” Drug Design, Development and Therapy, vol. 9, pp. 3989–4104, 2015. View at: Publisher Site | Google Scholar
  51. J. Wang, E. S. H. Chu, H. Y. Chen et al., “microRNA-29b prevents liver fibrosis by attenuating hepatic stellate cell activation and inducing apoptosis through targeting PI3K/AKT pathway,” Oncotarget, vol. 6, no. 9, pp. 7325–7338, 2015. View at: Publisher Site | Google Scholar
  52. M.-C. Jung and G. R. Pape, “Immunology of hepatitis B infection,” The Lancet Infectious Diseases, vol. 2, no. 1, pp. 43–50, 2002. View at: Publisher Site | Google Scholar
  53. R. Blomhoff and K. Wake, “Perisinusoidal stellate cells of the liver: important roles in retinol metabolism and fibrosis,” The FASEB Journal, vol. 5, no. 3, pp. 271–277, 1991. View at: Publisher Site | Google Scholar
  54. X. Gou, Q. Tao, Q. Feng et al., “Urine metabolic profile changes of CCl4-liver fibrosis in rats and intervention effects of Yi Guan Jian decoction using metabonomic approach,” BMC Complementary and Alternative Medicine, vol. 13, no. 1, p. 123, 2013. View at: Publisher Site | Google Scholar

Copyright © 2019 Guangyao Li 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.