Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 4863015 | https://doi.org/10.1155/2020/4863015

Shanghui Liu, Run Wang, Yan Lou, Jia Liu, "Uncovering the Mechanism of the Effects of Pien-Tze-Huang on Liver Cancer Using Network Pharmacology and Molecular Docking", Evidence-Based Complementary and Alternative Medicine, vol. 2020, Article ID 4863015, 12 pages, 2020. https://doi.org/10.1155/2020/4863015

Uncovering the Mechanism of the Effects of Pien-Tze-Huang on Liver Cancer Using Network Pharmacology and Molecular Docking

Academic Editor: Victor Kuete
Received08 Jun 2020
Revised20 Aug 2020
Accepted26 Aug 2020
Published07 Sep 2020


Pien-Tze-Huang (PTH) has a long history in the treatment of liver cancer. However, its molecular mechanism of action remains unclear. TCMSP and TCM were used to collect the active ingredients. Bioactive compounds targets were predicted by reverse pharmacophore models. The antiliver cancer targets of PTH were selected by gene comparison of liver cancer in the GEO database. Molecular docking was used to verify the binding activity of the targets and the active ingredients. The DAVID was used to analyze the gene function and signal pathway. A model was built with Cytoscape. The core genes were obtained by PPI network. We screened the 4 main medicinal ingredients of PTH to obtain 16 active ingredient, 190 potential targets, and 6 core genes. We found that active small molecules exert anticancer effects by multiple pathways. The core genes were involved in multiple biological processes. We also found that eight chemical components play a greater role in inhibiting liver cancer. PTH achieves the effect of inhibiting liver cancer through the synergistic effect of multiple components, multiple targets, and multiple pathways. This study provides a potential scientific basis for further elucidating the molecular mechanism of action of PTH against liver cancer.

1. Introduction

Liver cancer is a common clinical malignant tumor of liver. The incidence ranks the 5th in the world [1]. The fatality rate ranks the 2nd in the world [2]. Currently surgical resection is the main treatment. This is because most of the patients are already in the late stage once diagnosed, so that they can only be treated by drug instead of surgery. Clinically, first-line treatment drugs such as sorafenib are mainly used, but long-term use makes it easy for patients to have tolerance and a variety of adverse reactions [3, 4]. Therefore, exploring effective drugs to prevent and treat liver cancer is becoming an urgent issue. Many clinical studies have proved that Chinese medicine has good effects in improving the symptoms of patients with advanced liver cancer [5], reducing tumor recurrence [6], and controlling disease progression [7].

Pien-Tze-Huang (PTH), a classic Chinese medicine prescription, is made from Chinese ingredients such as Rhizoma Notoginseng (Sanqi in Chinese), Moschus (Shexiang in Chinese), Calculus Bovis (Niuhuang in Chinese), and Snake Gall (Shedan in Chinese) [8]. It is a proprietary Chinese medicine with the function of detoxification, anti-inflammation, and immune regulation [9, 10]. PTH has a history of more than 500 years and has been used as a folk medicine in China and Southeast Asia to treat various inflammation-related diseases such as hepatitis. In China, PTH is clinically used to protect the liver [11, 12], treat hepatitis [13], improve liver fibrosis [14], and treat liver cancer [1517]. Studies on the effects of PTH on lung cancer [18], intestinal cancer [19, 20], osteosarcoma [21], and liver cancer [22] have been reported, but the active ingredients of PTH on liver cancer and the relevant molecular mechanisms are still unclear.

Network pharmacology was developed on the basis of the rapid development of system biology, revealing the complex network relationship between drug-gene-target-disease, systematically observing the intervention and influence of drugs on the disease network, and revealing the molecular mechanism of drug therapy [23, 24]. Traditional Chinese medicine compounds are difficult to study because of their multiple herbs, complex components, and unclear mechanism of action. Network pharmacology provides a new way of thinking for the study of complex traditional Chinese medicine systems [25]. Molecular docking technology is to place small drug molecules (ligands) in the binding area of macromolecular targets (receptors) through computer simulation, predict the binding energy (binding affinity) and binding mode (conformation) of the two by calculating physical and chemical parameters, and then find the lowest energy conformation when the ligands and receptors combine in their active areas [26]. Low binding energy is the basis of stable binding between molecules. The aim of molecular docking is to find the best binding position when the binding energy is the lowest.

In this paper, with the help of network pharmacology, the bioactive compounds and possible molecular network mechanism of PTH against liver cancer are analyzed from a systematic and molecular level. The inhibition of PTH on liver cancer is verified by microarray data of liver cancer in GEO database and molecular docking technology. At the same time, it provides a method to identify the active molecules and target genes involved in complex diseases.

2. Materials and Methods

2.1. Screening of Active Ingredients in PTH

By searching TCM Database@Taiwan (http://tcm.cmu.edu.tw/), TCMSP (http://tcmspw.com/tcmsp.php), and related literature, we collected the main chemical components of Calculus Bovis, Moschus, Rhizoma Notoginseng, and Snake Gall in PTH, respectively. The oral bioavailability (OB ≥ 30%), drug-likeness (DS ≥ 0.18), and clear biological activity report [13, 27] were employed to identify the possible bioactive compounds in PTH. 2D and 3D structures and structural graphics of the above compounds were obtained using the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) and saved in mol2 format.

2.2. Prediction of Bioactive Compounds Targets for Liver Cancer

Through PharmMapper database (http://www.lilab-ecust.cn/pharmmapper/), we uploaded the mol2 format file of main active ingredient of PTH and used reverse pharmacophore method to get the information of interaction with compound and target protein. The target protein structures were downloaded and accurate information about proteins was further extracted from UniProt database (https://sparql.uniprot.org/), which is a comprehensive resource of protein sequence and rich functional information resources.

RNA sequencing data of liver cancer were obtained from the GSE76427 microarray in the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The mRNA sequencing data contained 115 liver cancer tissues and 52 normal liver tissues. The data used in this article was from GEO open databases and does not require informed consent. Differential analysis was performed using the limma package in R to obtain DEmRNA with thresholds of |log2 fold change (FC)| > 1 and value < 0.05.

Venn diagrams were used to further calculate the intersection of DEmRNA and the potential target genes, which were considered to be the target genes for predicting the active ingredients of PTH.

2.3. Active Ingredient-Target Protein Molecular Docking

To investigate the interaction relationship and mechanisms of action between candidate compounds and targets, we conducted a molecular docking study to explore binding activity of the two. We used the RCSB PDB (http://www.rcsb.org/pdb/home/home.do) database to retrieve and download the 3D structure file of the target protein as the receptor. The ligands were prepared based on ChemDraw Ultra 8.0 software for constructing two-dimensional structures and saved in mol2 format. Both of them were coimported into AutoDock Vina for molecular docking calculation. The docking results use binding energy as a reference to screen for the most active ligand molecules and target genes. At present, there is no unified standard for binding energy screening. Binding energy less than 0 means that ligand and receptor can bind spontaneously. The smaller the binding energy is, the more stable the binding between ligand and receptor is. After consulting the literature, it was considered that PTH had better binding activity with binding energy less than −5.0 kJ·mol−1 [28], which was the basis for screening candidate target of active ingredient.

2.4. Go and Pathways Analysis of Candidate Targets

Potential targets with binding energies less than 5.0 kJ·mol−1 were uploaded to the DAVID database (https://david.ncifcrf.gov/) for Gene Ontology (GO) and pathway enrichment analyses (KEGG). The target gene name list was defined as Homo sapiens, and the threshold of significant difference was set at . The results were plotted in R language.

2.5. Component-Target-Pathway Network Construction

To explore the mechanism of PTH with multiple ingredients, targets, and pathways, active ingredient-target-pathway network was built by Cytoscape 3.7.2, which has good binding power of active ingredient and target. Nodes of different colors indicate molecules, target proteins, and associated annotation pathways, and edges indicate molecule-target protein and target protein-pathway interrelationships.

2.6. Construction of PPI Network and Key Subnetwork

Based on target protein of PTH, the String (https://string-db.org/) database was used to construct a protein PPI network. Using Cytoscape 3.7.2 software, we selected core genes through maximal clique centrality (MCC) topological analysis method in cytoHubba and built the key subnetwork of PTH target protein.

2.7. Visualization of Docking between Core Genes and Active Ingredients

The target protein represented by the core genes and the active ingredient of interaction were analyzed to find the active ingredient with the minimum binding energy. Molecular visualization software PyMOL can be used to show docking effect conformation.

3. Results

3.1. Screening Results of Active Ingredients in PTH

By mining TCMSP and TCM, the main chemical components, such as Calculus Bovis, Moschus, Rhizoma Notoginseng, and Snake Gall, in PTH were identified. During the screening process, it was found that some of the components failed to meet the requirements of OB ≥ 30% and DL ≥ 0.18. However, those components in PTH were relatively high and had been proved to have clearly biological activity by literatures; they were also included in the study of active components. At last, 18 PTH active ingredients were collected, including 6 Rhizoma Notoginseng, 7 Calculus Bovis, 1 Moschus, and 4 Snake galls. 16 remained after weight removal. Detailed data is listed in Table 1.

Traditional Chinese medicine (TCM)No.Active ingredientChemical Molecular Formula

Rhizoma Notoginseng1Notoginsenoside R1C47H80O18
2Ginsenoside ReC48H82O18
3Ginsenoside Rg1C42H72O14
4Ginsenoside Rb1C54H92O23
5Astragaloside AC41H68O14
6Ginsenoside RdC48H82O18

Calculus Bovis1TaurineC2H7NO3S
2Taurocholic acidC26H45NO7S
3Glycocholic acid hydrateC26H43NO6
4Cholic acidC24H40O5
5Glycodeoxycholic acid hydrateC26H43NO5
6Chenodeoxycholic acidC24H40O4
7Deoxycholic acidC24H40O4


Snake galls1Taurocholic acidC26H45NO7S
2Glycocholic acid hydrateC26H43NO6
3Sodium taurodeoxylateC26H44NO6S
4Taurochenodeoxycholic acidC26H45NO6S

3.2. Prediction Results of Bioactive Compounds Targets in PTH for Liver Cancer

Through the PharmMapper, 4505 potential targets were obtained by analysing 16 active ingredients in PTH. The target names were input into UniProt database, and all the target names were corrected to the gene names of the target; 986 target genes were obtained. 307 target genes were obtained after further gene deduplication. Then, 307 target genes and 15385 differentially expressed mRNA (DEmRNA) of GEO were intersected, and 197 potential liver cancer related targets were screened by the intersection of Venn (Figure 1).

3.3. Molecular Docking Results

With 16 active ingredients as ligands and 197 target protein genes as receptors, 614 groups of corresponding relationships were formed. The binding energy of each group can be calculated by molecular docking, and the binding energy of docking results can be displayed by heatmap (Figure 2). Red indicates that the binding energy is negative. Green indicates positive. Gray indicates that protein gene is not the target of active ingredient. Being dominated by red means that most of the active ingredients in PTH had good interaction and binding activity with target. The 190 target genes with better binding activity were obtained by screening out targets with binding energy less than −5.

3.4. GO and KEGG Analysis of Target Genes

We used DAVID database to analyze the process of GO biological function and the enrichment of KEGG metabolic pathway of 190 target genes with good binding activity. The threshold was set to . There were 71, 45, and 25 enrichment items of biological processes (BP), cell component (CC), and molecular function (MF), respectively. Top 10 records were screened from small to large values. Items related to BP mainly focused on transcription initiation from RNA polymerase II promoter, steroid hormone mediated signaling pathway, protein autophosphorylation, positive regulation of apoptotic process, protein phosphorylation, response to lipopolysaccharide, etc. (Figure 3). The items related to CC mainly involved cytosol, cytoplasm, nucleoplasm, extracellular exosome, nucleus, etc. (Figure 4). MF related items mainly involved protein binding, same protein binding, zinc ion binding, protein homopolymer activity, ATP binding, protein kinase activity, etc. (Figure 5).

KEGG signaling pathway analysis is a process in which multiple proteins interact to regulate cell function and metabolism. 12 signal pathways with mainly include metabolic pathways, pathways in cancer, biosynthesis of antibiotics, proteoglycans in cancer, focal adhesion, purine metabolism, etc. (Figure 6).

3.5. Construction Results of Component-Target-Pathway Network

The component-target-pathway network was constructed by using Cytoscape 3.7.2. Green dot indicated the active component of PTH. Blue was the target protein of the active component, and red was KEGG pathway (Figure 7). It can be seen from the figure that 16 active components of PTH act on 190 targets and distribute in 12 different metabolic pathways. The network degree value and the number of target genes of each active component in PTH were calculated (Table 2). The degree of 15 compounds was more than 30. In particular, the numbers of target genes targeted by sodium taurodeoxylate, taurochenodeoxycholic acid, and glycocholic acid hydrate were 48, 46, and 43. Only the number of target genes targeted by Taurine was as few as 2, and the binding energies were −4.3 and −3.2. The target factor less than −5 was 0. PTH has the mechanism of multicomponent, multitarget, and multichannel.

No.Chemical compoundDegree

1Sodium taurodeoxylate48
2Taurochenodeoxycholic acid46
3Glycocholic acid hydrate43
4Glycodeoxycholic acid hydrate41
5Ginsenoside Re40
6Ginsenoside Rg140
7Notoginsenoside R140
8Chenodeoxycholic acid39
9Deoxycholic acid38
10Astragaloside A34
11Cholic acid37
12Ginsenoside Rb137
13Ginsenoside Rd37
15Taurocholic acid34

3.6. Construction Results of PPI Network and Key Subnetwork

The PPI network was constructed by analysing 190 target proteins in STRING database. The minimum required interaction score of PPI network is 0.4 by default. The disconnected nodes are hidden. After screening, 167 targets can interact with proteins, 492 edges represent the interaction relationship between proteins and the average degree of each node is 5.9. PPI network of PTH is shown in Figure 8. Using the MCC algorithm of cytoHubba plug-in in Cytoscape, a PPI key subnetwork composed of six genes with the highest score is obtained. The PPI key subnetwork with ESR1, AR, ALB, NOTCH1, ERBB2, and IGF1R has the most connectivity and contains more genes than other subnets, which can avoid missing the key factors related to disease (Figure 9).

3.7. Visual Analysis of Docking between Core Genes and Active Ingredients

Molecular docking of the core genes to the active ingredients with high mutual affinity results is shown in Table 3. The corresponding molecular docking effects are shown in Figure 10. It shows the best conformation when the binding energy between the small molecule and the target protein is the lowest and shows the interaction mode between the drug molecule and the core target.

Core geneActive ingredientCombined energy

ESR1Glycocholic acid hydrate−6.4
ARChenodeoxycholic acid−7.3
ALBChenodeoxycholic acid−9.5
NOTCH1Ginsenoside Rb1−7.4
ERBB2Taurochenodeoxycholic acid−7.4
IGF1RSodium taurodeoxylate−7.1

4. Discussion

There are 16 main chemical components in PTH, which act on 307 target genes in human body. 197 target genes related to liver cancer were obtained by crossing with hepatoma genes in GEO database. In order to verify the binding activity of chemical components and target genes, 16 chemical components were combined with 197 target genes for molecular docking, and 190 target genes with binding energy less than −5 were finally obtained. Then, GO enrichment analysis, KEGG pathway analysis, PPI network, and key subnet construction as well as the establishment of component-target-pathway network were carried out, and the synergistic mechanism of PTH multicomponent, multitarget, and multipathway was obtained.

PTH has many active target genes against liver cancer. Six core active target genes (ESR1, ERBB2, IGF1R, AR, NOTCH1, and ALB) were obtained by constructing PPI subnetwork. The antihepatocellular carcinoma mechanism of the six core genes was obtained by the comprehensive analysis of the component-target-pathway network, the relevant literature retrieval, and the screening of chemical components with minimum binding energy.

ESR1 is an estrogen receptor, which has a wide range of effects on the heterogeneous metabolism of human liver [29]. Liver is the target organ of sex hormone. The sex hormone level and the variants of sex hormone specific gene expression are supposed to be the main cause of liver disease [30]. It has been found that ESR1 is a direct downstream target of mir-393-3p, and mir-939-3p/ESR1 axis may be a potential new target for the treatment of liver cancer [31]. Verrier et al. found that ESR1 is one of the key host factors of the liver virus and a target for the inhibition of the liver virus [32]. Wang et al. also revealed that ESR1 is probably associated with HCC survival [33]. In this study, it was shown by target prediction and pathway enrichment that the bioactive compounds, glycolic acid hydrate, ginsenoside Rg1, and notoginsenoside R1, all act on the proteoglycans in cancer signal pathway and inhibit the progress of liver cancer by inhibiting the ESR1 transduction on the pathway.

Tyrosine kinase receptor 2 (ERBB2) is a member of the epidermal growth factor receptor (EGFR) family. Patients with high ERBB2 expression are prone to tumor metastasis and have a short survival period, which has become an ideal target of tumor immunobiotherapy and a hot molecule in the field of tumor therapy. It has been found that mir-296-5p inhibits EMT related metastasis of liver cancer through ERBB2 signaling pathway [34, 35]. Mir-375 upregulates cell proliferation and apoptosis by targeting ERBB2 and inhibits the growth of human hepatoma cells [36]. Yap directly regulates the transcription of ERBB2 and influences the proliferation of hepatocytes through ERBB2 signaling pathway [37]. Target prediction and pathway enrichment showed that the bioactive compounds, muscone and taurochenodeoxycholic acid, target ERBB2 through five signal pathways (pathways in cancer, prostate cancer, proteoglycans in cancer, focal adhesion, and bladder cancer) to inhibit the growth of liver cancer cells.

The insulin-like growth factor 1 receptor (IGF1R) is a transmembrane tyrosine protein located on chromosome 15q25-q26. In some malignancies, overexpression of IGF1R promotes tumor growth, invasion, and metastasis [38, 39]. Fan et al. believed that targeted IGF1R may be a promising therapeutic method for liver cancer [40]. Lin et al. believed that IGF1R and GRK2 are negatively correlated in liver cancer, and IGF1R may be a potential marker of poor prognosis in liver cancer [41]. Ye et al. found that the decrease of inhibitory effect on driving oncogene IGF1R and MAPK1 gene may trigger EMT and stem cell transformation of HCC [42]. According to target prediction and pathway enrichment, the bioactive compound stadium oxygen can act on the signal pathways of prostate cancer, pathways in cancer, proteoglycans in cancer, and focal adhesion and controls the progress of liver cancer by inhibiting the overexpression of IGF1R.

Androgen receptor (AR) may play a key role in influencing the progression of hepatocellular carcinoma at different stages [43, 44]. AR may be a marker for the response of patients with liver cancer to sorafenib [45]. AR signal transduction is involved in many aspects of the pathogenesis of liver cancer. The expression of AR-SV in liver cancer leads to the progress of liver cancer and the resistance to traditional AR antagonists. The successful therapy targeting AR-SVs is beneficial to liver cancer [46]. Target prediction and pathway enrichment revealed that 14 bioactive compounds (notoginsenoside R1, ginsenoside Re, ginsenoside Rg1, ginsenoside Rb1, ginsenoside Rd, chenodeoxycholic acid, cholic acid, glycodeoxycholic acid hydrate, deoxycholic acid, taurocholic acid, astragaloside A, ginsenoside Rd, sodium taurodeoxylate, and taurochenodeoxycholic acid) act on pathways in cancer signaling pathway and play a therapeutic role in liver cancer by regulating the expression of AR.

Notch signaling pathway is highly conserved and widely expressed in almost all tissues and organs, which plays a regulatory role in cell generation, development, differentiation, and proliferation. It has been reported that abnormal activation of Notch signal is related to the metastasis of hepatocellular carcinoma, and NOTCH1 can activate RNF187 promoter to promote the invasion and metastasis of hepatocellular carcinoma [47]. Inhibition of NOTCH1 signaling pathway may inhibit the development of HBx-induced hepatocellular carcinoma, and NOTCH1 may become a therapeutic target for hepatocellular carcinoma [48]. The bioactive compounds ginsenoside Rb1 and ginsenoside Rd inhibit the development of hepatocellular carcinoma by inhibiting the activation of NOTCH1 signal.

Albumin (ALB) is the most important protein in human plasma, which can maintain the nutrition and osmotic pressure and also reflect the synthesis function of the liver. It is found that the activity of ALB in patients with hepatitis, cirrhosis, and liver cancer is decreased significantly. ALB is considered to play a pivotal role in the inflammatory process and has been known as a prognostic indicator in several types of cancer [49]. Peng et al. found that ALB is a biomarker, which is related to the early diagnosis, metastasis, prognosis, or treatment of HCC [50]. The levels of ALB and GGT before operation are potential biomarkers to predict the prognosis of patients with liver cancer undergoing radical resection [51]. The evaluation of liver fibrosis degree is an important basis for the clinical diagnosis and treatment of HCC. Also, ALB is valuable while evaluating the liver fibrosis degree [52]. Chenodeoxycholic acid and cholic acid have good binding activity with ALB target genes, which can effectively regulate the activity of ALB, so as to improve the treatment status of patients with liver cancer.

Our KEGG pathway annotation analysis showed that some pathways were closely related to liver cancer [53]. The cyclin dependent kinase CDK2 is the target of tumor treatment [54]. The enrichment of 17 target genes in the pathways in cancer is primarily associated with liver cancer, such as Hsp90AB1e; the p53 signaling pathway may be an important molecular mechanism for the occurrence and development of liver cancer, which is consistent with the research results of Meng et al. [55] and Fan et al. [56]. The amplification of CCNE1 in this pathway is significantly related to the metastasis of liver cancer, which has been identified as an important target for tumor cells [57]. He et al. found that downregulation of CBL expression can inhibit the development of HCC [58]. In our study, the target genes of liver cancer are most abundant in metabolic pathways. Many studies have confirmed that this pathway is closely related to the occurrence and progress of HCC [59]. Sun et al. found that targeting specific metabolic pathways can prevent recurrence of liver cancer [60]. There are many target genes enriched in focal adhesion pathway, in which focal adhesion kinase and β-catenin collaboratively induce the production of HCC [61]. Also, many liver cancer target genes are enriched in the “proteoglycans in cancer” pathway. Proteoglycans are attractive biomarkers and therapeutic targets in hepatocellular carcinoma [62]. Purine metabolism pathway is also related to liver cancer. For example, the loss of dual specific tyrosine phosphorylation regulated kinase 3 can activate purine metabolism, which can promote the progress of liver cancer [63].

The degree of active ingredient indicates the number of correlations between the ingredient and target. According to the component-target-pathway network, sodium taurodeoxylate, taurochenodeoxycholic acid, glycocholic acid hydrate, chenodeoxycholic acid, ginsenoside Rd, muscone, cholic acid, and ginsenoside Rb1 are 8 bioactive compounds with high degree. They act on 6 core target genes and perform well in molecular docking. These provide a basis for further research on the bioactive compounds of antihepatoma.

5. Conclusions

PHT is mainly made from four compatible medicinal materials, which are Calculus Bovis, Moschus, Rhizoma Notoginseng, and Snake Gall. Its 16 chemical components act on different pathways by regulating different target proteins and collaboratively inhibit liver cancer, which shows the integrity and systematicness of Chinese traditional medicine prescriptions. ① Active small molecule ligands regulate liver cancer target gene receptors and inhibit the occurrence, development, and metastasis of liver cancer cells by interfering in tumor-related signaling pathways such as pathways in cancer, P53, proteoglycans in cancer, focal adhesion, metabolic and purine metabolism. ② Core target genes, ESR1, ERBB2, IGF1R, AR, NOTCH1,and ALB, are mainly involved in the transcription of RNA polymerase II promoters, steroid hormone mediated signaling pathway, protein autophosphorylation, positive regulation of apoptosis process, protein phosphorylation, transmembrane receptor protein tyrosine kinase signaling pathway, response to lipopolysaccharide, and other biological processes, which inhibits liver cancer. ③ Eight bioactive compounds (sodium taurodeoxylate, taurochenodeoxycholic acid, glycocholic acid hydrate, chenodeoxycholic acid, ginsenoside Rd, muscone, cholic acid, and ginsenoside Rb1) are effective to inhibit liver cancer. This study provides a basis for the molecular mechanism and biological process about how PHT inhibits liver cancer and lays a foundation for subsequent studies.

Data Availability

The datasets generated for this study can be found at NCBI using accession number GSE76427 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE76427).

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

S. H. L. conceptualized the study, S. H. L. and J. L. contributed to methodology, R. W. and J. L. validated the study, S. H. L. was responsible for formal analysis, and S. H. L. investigated the study. S. H. L. and J. L. were responsible for writing—original draft preparation, and S. H. L. contributed to writing—review and editing. Visualization was performed by J. L., and supervision was performed by Y. L. All authors have read and agreed the published version of the manuscript.


  1. X. Sun, L. Li, H.-g. Ma et al., “Bisindolylmaleimide alkaloid BMA-155Cl induces autophagy and apoptosis in human hepatocarcinoma HepG-2 cells through the NF-κB p65 pathway,” Acta Pharmacologica Sinica, vol. 38, no. 4, pp. 524–538, 2017. View at: Publisher Site | Google Scholar
  2. J. Zhou, H. Sun, Z. Wang et al., “Guidelines for diagnosis and treatment of primary liver cancer in China: 2017 edition,” Liver Cancer, vol. 7, pp. 235–260, 2018. View at: Google Scholar
  3. M. Santoni, A. Conti, U. De Giorgi et al., “Risk of gastrointestinal events with sorafenib, sunitinib and pazopanib in patients with solid tumors: a systematic review and meta-analysis of clinical trials,” International Journal of Cancer, vol. 135, no. 4, pp. 763–773, 2014. View at: Publisher Site | Google Scholar
  4. H. C. Randrup, D. Grimm, J. Bauer et al., “Effects and side effects of using sorafenib and sunitinib in the treatment of metastatic renal cell carcinoma,” International Journal of Molecular Sciences, vol. 18, p. E461, 2017. View at: Google Scholar
  5. K.-y. Tang, S.-l. Du, Q.-l. Wang, Y.-f. Zhang, and H.-y. Song, “Traditional Chinese medicine targeting cancer stem cells as an alternative treatment for hepatocellular carcinoma,” Journal of Integrative Medicine, vol. 18, no. 3, pp. 196–202, 2020. View at: Publisher Site | Google Scholar
  6. Y. Cheng, S. Ni, Y. Chen, Q. Ling, and J. Chen, “Erzhu Qinggan Jiedu recipe improves the clinical outcome of hepatocellular cancer after surgical resection: a case-control retrospective study,” Internal Medicine Journal, 2020. View at: Publisher Site | Google Scholar
  7. X. Wang, N. Wang, F. Cheung, L. Lao, C. Li, and Y. Feng, “Chinese medicines for prevention and treatment of human hepatocellular carcinoma: current progress on pharmacological actions and mechanisms,” Journal of Integrative Medicine, vol. 13, no. 3, pp. 142–164, 2015. View at: Publisher Site | Google Scholar
  8. Pharmacopoeia Commission, Pharmacopoeia of the People’s Republic of China, Chinese Medical Science and Technology Press, Beijing, China, 2015.
  9. X. Qiu, H. Luo, X. Liu et al., “Therapeutic potential of Pien Tze Huang on experimental autoimmune encephalomyelitis rat,” Journal of Immunology Research, vol. 2018, Article ID 2952471, p. 10, 2018. View at: Publisher Site | Google Scholar
  10. Y. Deng, H. Luo, J. Shu et al., “Pien Tze Huang alleviate the joint inflammation in collagen-induced arthritis mice,” Chinese Medicine, vol. 15, no. 1, p. 30, 2020. View at: Publisher Site | Google Scholar
  11. K. K. H. Lee, W. H. Kwong, F.-t. Chau, D. T. Yew, and W. Y. Chan, “Pien Tze Huang protects the liver against carbon tetrachloride-induced damage,” Pharmacology and Toxicology, vol. 91, no. 4, pp. 185–192, 2002. View at: Publisher Site | Google Scholar
  12. J. Y. Zhao, D. X. Chen, S. Lin et al., “Hepatoprotective activity of Pien Tze Huang Gan Bao in CCl4-induced chronic liver injury models,” International Journal of Clinical and Experimental Medicine, vol. 12, pp. 11324–11332, 2019. View at: Google Scholar
  13. L. Huang, Y. Zhang, X. Zhang et al., “Therapeutic potential of Pien-Tze-Huang: a review on its chemical composition, pharmacology, and clinical application,” Molecules, vol. 24, no. 18, p. 3274, 2019. View at: Publisher Site | Google Scholar
  14. H. Zheng, X. Wang, Y. Zhang, L. Chen, L. Hua, and W. Xu, “Pien-Tze-Huang ameliorates hepatic fibrosis via suppressing NF-κB pathway and promoting HSC apoptosis,” Journal of Ethnopharmacology, vol. 244, Article ID 111856, 2019. View at: Publisher Site | Google Scholar
  15. Y. Y. Xu and E. X. Yu, “Clinical analysis of the effect of Pien Tze Huang in treatment of 42 patients with moderate or advanced liver cancer,” Shanghai Journal of Traditional Chinese Medicine, vol. 12, pp. 4-5, 1994. View at: Google Scholar
  16. S. L. Zhao and J. Pan, “A clinical trial of combined use of Pien Tze Huang and chemotherapy in the treatment of primary liver cancer,” Medical World, pp. 49–51, 2006. View at: Google Scholar
  17. L. P. Meng and X. F. Gu, “Clinical observation of the changes of immune indexes of primary liver cancer treated with Pien Tze Huang and TACE,” Jiujiang Medical Journal, vol. 1, pp. 31–33, 2008. View at: Google Scholar
  18. Y. H. Hu, Q. L. Cui, and D. Y. Ma, “Inhibitory effect and mechanism of Pien Tze Huang on proliferation, migration, and invasion of lung cancer cells,” Chinese Traditional and Herbal Drugs, vol. 50, pp. 5527–5531, 2019. View at: Google Scholar
  19. J. Lin, J. Feng, Y. Jin, Z. Yan, Z. Lai, and J. Peng, “Pien Tze Huang suppresses VEGF-C-mediated lymphangiogenesis in colorectal cancer,” Oncology Reports, vol. 36, no. 6, pp. 3568–3576, 2016. View at: Publisher Site | Google Scholar
  20. Z. Chen, A. Shen, L. Liu et al., “Pien Tze Huang induces apoptosis and inhibits proliferation of 5-fluorouracil-resistant colorectal carcinoma cells via increasing miR-22 expression,” Experimental and Therapeutic Medicine, vol. 14, no. 4, pp. 3533–3540, 2017. View at: Publisher Site | Google Scholar
  21. Y. Fu, L. Zhang, Z. Hong et al., “Methanolic extract of Pien Tze Huang induces apoptosis signaling in human osteosarcoma MG63 cells via multiple pathways,” Molecules, vol. 21, no. 3, p. 283, 2016. View at: Publisher Site | Google Scholar
  22. F. Qi, S. Zhou, L. Li et al., “Pien Tze Huang inhibits the growth of hepatocellular carcinoma cells by upregulating miR-16 expression,” Oncology Letters, vol. 14, no. 6, pp. 8132–8137, 2017. View at: Google Scholar
  23. M. Kibble, N. Saarinen, J. Tang, K. Wennerberg, S. Mäkelä, and T. Aittokallio, “Network pharmacology applications to map the unexplored target space and therapeutic potential of natural products,” Natural Product Reports, vol. 32, no. 8, pp. 1249–1266, 2015. View at: Publisher Site | Google Scholar
  24. S. Zhao and R. Iyengar, “Systems pharmacology: network analysis to identify multiscale mechanisms of drug action,” Annual Review of Pharmacology and Toxicology, vol. 52, no. 1, pp. 505–521, 2012. View at: Publisher Site | Google Scholar
  25. S. Li and B. Zhang, “Traditional Chinese medicine network pharmacology: theory, methodology and application,” Chinese Journal of Natural Medicines, vol. 11, no. 2, pp. 110–120, 2013. View at: Publisher Site | Google Scholar
  26. D. B. Kitchen, H. Decornez, J. R. Furr, and J. Bajorath, “Docking and scoring in virtual screening for drug discovery: methods and applications,” Nature Reviews Drug Discovery, vol. 3, no. 11, pp. 935–949, 2004. View at: Publisher Site | Google Scholar
  27. Q. M. Huang, Y. N. Tai, Y. L. Zhu, S. J. Shen, H. T. Guo, and M. Q. Huang, “Simultaneous determination of thirteen constituents in Pientzehuang by UPLCQQQ-MS,” Chinese Traditional Patent Medicine, vol. 40, pp. 101–105, 2018. View at: Google Scholar
  28. K.-Y. Hsin, Y. Matsuoka, Y. Asai et al., “SystemsDock: a web server for network pharmacology-based prediction and analysis,” Nucleic Acids Research, vol. 44, no. W1, pp. W507–W513, 2016. View at: Publisher Site | Google Scholar
  29. D. Wang, R. Lu, G. Rempala, and W. Sadee, “Ligand-free estrogen receptor α (ESR1) as master regulator for the expression of CYP3A4 and other cytochrome P450 enzymes in the human liver,” Molecular Pharmacology, vol. 96, no. 4, pp. 430–440, 2019. View at: Publisher Site | Google Scholar
  30. P. Kur, A. Kolasa-Wołosiuk, K. Misiakiewicz-Has, and B. Wiszniewska, “Sex hormone-dependent physiology and diseases of liver,” International Journal of Environmental Research and Public Health, vol. 17, no. 8, p. 2620, 2020. View at: Publisher Site | Google Scholar
  31. F. Chen, X. Y. Ni, L. X. Chen et al., “miR-939-3p promotes epithelial-mesenchymal transition and may be used as a prognostic marker in hepatocellular carcinoma,” Oncology Letters, vol. 19, pp. 2727–2732, 2020. View at: Google Scholar
  32. E. R. Verrier, A. Weiss, C. Bach et al., “Combined small molecule and loss-of-function screen uncovers estrogen receptor alpha and CAD as host factors for HDV infection and antiviral targets,” Gut, vol. 69, no. 1, pp. 158–167, 2020. View at: Publisher Site | Google Scholar
  33. X. Wang, X. Liao, K. Huang et al., “Clustered microRNAs HSA-miR-221-3p/HSA-miR-222-3p and their targeted genes might be prognostic predictors for hepatocellular carcinoma,” Journal of Cancer, vol. 10, no. 11, pp. 2520–2533, 2019. View at: Publisher Site | Google Scholar
  34. D.-M. Shi, L.-X. Li, X.-Y. Bian et al., “miR-296-5p suppresses EMT of hepatocellular carcinoma via attenuating NRG1/ERBB2/ERBB3 signaling,” Journal of Experimental & Clinical Cancer Research, vol. 37, no. 1, p. 294, 2018. View at: Publisher Site | Google Scholar
  35. D.-M. Shi, X.-L. Shi, K.-L. Xing, H.-X. Zhou, L.-L. Lu, and W.-Z. Wu, “miR-296-5p suppresses stem cell potency of hepatocellular carcinoma cells via regulating Brg1/Sall4 axis,” Cellular Signalling, vol. 72, Article ID 109650, 2020. View at: Publisher Site | Google Scholar
  36. L. Li, L. Jia, and Y. Ding, “Upregulation of miR-375 inhibits human liver cancer cell growth by modulating cell proliferation and apoptosis via targeting ErbB2,” Oncology Letters, vol. 16, no. 3, pp. 3319–3326, 2018. View at: Google Scholar
  37. E. Y. Wang, J. C. Cheng, A. Thakur et al., “YAP transcriptionally regulates ErbB2 to promote liver cell proliferation,” Biochimica et Biophysica Acta. Gene Regulatory Mechanisms, vol. 1861, pp. 54–863, 2018. View at: Publisher Site | Google Scholar
  38. L. Zhao, L.-R. He, R. Zhang et al., “Low expression of IGFBP-3 predicts poor prognosis in patients with esophageal squamous cell carcinoma,” Medical Oncology, vol. 29, no. 4, pp. 2669–2676, 2012. View at: Publisher Site | Google Scholar
  39. H. Cheng, J. Xue, S. Yang et al., “Co-targeting of IGF1R/mTOR pathway by miR-497 and miR-99a impairs hepatocellular carcinoma development,” Oncotarget, vol. 8, no. 29, pp. 47984–47997, 2017. View at: Publisher Site | Google Scholar
  40. Y. Fan, J. J. Li, H. J. Deng et al., “GSTZ1-1 deficiency activates NRF2/IGF 1R axis in HCC via accumulation of oncometabolite succinylacetone,” The EMBO Journal, vol. 38, Article ID e101964, 2019. View at: Publisher Site | Google Scholar
  41. S.-B. Lin, L. Zhou, Z.-Y. Liang, W.-X. Zhou, and Y. Jin, “Expression of GRK2 and IGF1R in hepatocellular carcinoma: clinicopathological and prognostic significance,” Journal of Clinical Pathology, vol. 70, no. 9, pp. 754–759, 2017. View at: Publisher Site | Google Scholar
  42. Y. N. Ye, J. C. Guo, P. Xiao et al., “Macrophages-induced long noncoding RNA H19 up-regulation triggers and activates the miR-193b/MAPK1 axis and promotes cell aggressiveness in hepatocellular carcinoma,” Cancer Letters, vol. 28, pp. 310–322, 2019. View at: Google Scholar
  43. Y. Xiao, Y. Sun, G. D. Liu et al., “Androgen receptor (AR)/miR-520f-3p/SOX9 signaling is involved in the altering hepatocellular carcinoma (HCC) cell sensitivity to the sorafenib therapy under hypoxia via increasing cancer stem cells phenotype,” Cancer Letters, vol. 444, pp. 175–187, 2018. View at: Google Scholar
  44. H. Song, N. Sun, L. Lin et al., “Splicing factor PRPF6 upregulates oncogenic androgen receptor signaling pathway in hepatocellular carcinoma,” Cancer Science, 2020. View at: Publisher Site | Google Scholar
  45. H.-C. Lai, W.-M. Chung, C.-M. Chang et al., “Androgen receptor enhances the efficacy of sorafenib against hepatocellular carcinoma through enriched EpCAM stemness,” Anticancer Research, vol. 40, no. 3, pp. 1285–1295, 2020. View at: Publisher Site | Google Scholar
  46. A. M. Dauki, J. S. Blachly, E. A. Kautto, S. Ezzat, M. H. Abdel-Rahman, and C. C. Coss, “Transcriptionally active androgen receptor splice variants promote hepatocellular carcinoma progression,” Cancer Research, vol. 80, no. 3, pp. 561–575, 2020. View at: Publisher Site | Google Scholar
  47. L. Zhang, J. Chen, J. Yong, L. Qiao, L. Xu, and C. Liu, “An essential role of RNF187 in Notch1 mediated metastasis of hepatocellular carcinoma,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 384, 2019. View at: Publisher Site | Google Scholar
  48. Y. Peng, D. J. Zhou, Y. S. Ye et al., “Epigenetic knockdown of Notch1 inhibits hepatitis B virus X protein-induced hepatocarcinogenesis of L02/HBx cells,” Oncology Reports, vol. 41, pp. 1151–1159, 2019. View at: Google Scholar
  49. W. Zhang, G. Zhangyuan, F. Wang et al., “High preoperative serum globulin in hepatocellular carcinoma is a risk factor for poor survival,” Journal of Cancer, vol. 10, no. 15, pp. 3494–3500, 2019. View at: Publisher Site | Google Scholar
  50. Y. Peng, Q. Wei, Y. He et al., “ALBI versus child-Pugh in predicting outcome of patients with HCC: a systematic review,” Expert Review of Gastroenterology & Hepatology, vol. 14, no. 5, p. 383, 2020. View at: Publisher Site | Google Scholar
  51. L. Wang, Q. Li, J. Zhang, and J. Lu, “A novel prognostic scoring model based on albumin and gamma-glutamyltransferase for hepatocellular carcinoma prognosis,” Cancer Management and Research, vol. 11, pp. 10685–10694, 2019. View at: Publisher Site | Google Scholar
  52. W. Liu, X. S. Li, W. H. Zheng et al., “Preoperative evaluation of the degree of liver fibrosis based on matter-element analysis using serological indicators in patients with hepatocellular carcinoma,” Bioscience Trends, vol. 13, pp. 70–76, 2019. View at: Publisher Site | Google Scholar
  53. B. Kim, H. C. Shin, Y. J. Heo et al., “CCNE1 amplification is associated with liver metastasis in gastric,” Pathology—Research and Practice, vol. 215, Article ID 152434, 2019. View at: Publisher Site | Google Scholar
  54. T. A. Chohan, H. Qian, Y. Pan et al., “Cyclin-dependent kinase-2as a target for cancer therapy:progress in the development of CDK2inhibitors as anti-cancer agents,” Current Medicinal Chemistry, vol. 22, pp. 237–263, 2015. View at: Google Scholar
  55. X. Meng, D. A. Franklin, J. Dong et al., “MDM2-p53 pathway in hepatocellular carcinoma,” Cancer Research, vol. 74, pp. 7161–7167, 2014. View at: Publisher Site | Google Scholar
  56. M. Fan, J. Shen, H. Liu et al., “Downregulation of PRRX1 via the p53-dependent signaling pathway predicts poor prognosis in hepatocellular carcinoma,” Oncology Reports, vol. 38, pp. 1083–1090, 2017. View at: Publisher Site | Google Scholar
  57. E. Coskunpinar, N. Akkaya, P. Yildiz et al., “The significance of HSP90AA1,HSP90AB1 and HSP90B1 gene polymorphisms in a Turkish population with non-small cell lung cancer,” Anticancer Research, vol. 34, pp. 753–757, 2014. View at: Google Scholar
  58. J. He, B. Xiao, X. Y. Li et al., “MiR-486-5p suppresses proliferation and migration of hepatocellular carcinoma cells through downregulation of the E3 ubiquitin ligase CBL,” BioMed Research International, vol. 2019, Article ID 2732057, p. 9, 2019. View at: Publisher Site | Google Scholar
  59. W. Q. Dai, L. Xu, X. N. Yu et al., “OGDHL silencing promotes hepatocellular carcinoma by reprogramming glutamine metabolism,” Journal of Hepatology, vol. 72, pp. 909–923, 2019. View at: Google Scholar
  60. L. Sun, L. Zhang, J. Chen et al., “Activation of tyrosine metabolism in CD13(+) cancer stem cells drives relapse in hepatocellular carcinoma,” Cancer Research and Treatment, vol. 52, pp. 604–621, 2020. View at: Publisher Site | Google Scholar
  61. N. Shang, H. Wang, T. Bank et al., “Focal adhesion kinase and β-catenin cooperate to induce hepatocellular carcinoma,” Hepatology, vol. 70, pp. 1631–1645, 2019. View at: Publisher Site | Google Scholar
  62. Y. Tanaka, R. Tateishi, and K. Koike, “Proteoglycans are attractive biomarkers and therapeutic targets in hepatocellular carcinoma,” International Journal of Molecular Sciences, vol. 19, p. E3070, 2018. View at: Publisher Site | Google Scholar
  63. F. Ma, Y. K. Zhu, X. Liu et al., “Dual-specificity tyrosine phosphorylation-regulated kinase 3 loss activates purine metabolism and promotes hepatocellular carcinoma progression,” Hepatology, vol. 70, pp. 1785–1803, 2019. View at: Publisher Site | Google Scholar

Copyright © 2020 Shanghui 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.

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.