Evidence-Based Complementary and Alternative Medicine
Volume 2020 (2020), Article ID 4863015, 12 pages
https://doi.org/10.1155/2020/4863015

Research Article

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

Shanghui Liu,¹ Run Wang,² Yan Lou,¹ and Jia Liu¹

¹School of Fundamental Sciences, China Medical University, Shenyang 110122, Liaoning, China
²Department of Neurosurgery, The 1st Hospital of China Medical University, Shenyang 110001, Liaoning, China

Correspondence should be addressed to Jia Liu; jliu03@cmu.edu.cn

Received 8 June 2020; Revised 20 August 2020; Accepted 26 August 2020; Published 7 September 2020

Academic Editor: Victor Kuete

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.

Abstract.

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 [15–17]. 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⁻¹ [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⁻¹ 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.

Table 1: Information of active ingredients in PTH.

Traditional Chinese medicine (TCM)
No.
Active ingredient
Chemical Molecular Formula

Rhizoma Notoginseng
1
Notoginsenoside R1
C47H80O18

2
Ginsenoside Re
C48H82O18

3
Ginsenoside Rg1
C42H72O14

4
Ginsenoside Rb1
C54H92O23

5
Astragaloside A
C41H68O14

6
Ginsenoside Rd
C48H82O18

Calculus Bovis
1
Taurine
C2H7NO3S

2
Taurocholic acid
C26H45NO7S

3
Glycocholic acid hydrate
C26H43NO6

4
Cholic acid
C24H40O5

5
Glycodeoxycholic acid hydrate
C26H43NO5

6
Chenodeoxycholic acid
C24H40O4

7
Deoxycholic acid
C24H40O4

Moschus
1
Muscone
C16H30O

Snake galls
1
Taurocholic acid
C26H45NO7S

2
Glycocholic acid hydrate
C26H43NO6

3
Sodium taurodeoxylate
C26H44NO6S

4
Taurochenodeoxycholic acid
C26H45NO6S

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).

Figure 1: Venn of GEO DEmRNA and target genes.

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.

Figure 2: Heatmap plot of molecular docking between active ingredients and target genes. Red represents higher binding activity, green represents lower binding activity levels, and gray represents no connection between target gene and active ingredient.

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).

Figure 3: Enriched gene ontology terms for biological processes from active ingredients of PTH.

Figure 4: Enriched gene ontology terms for cell component from active ingredients of PTH.

Figure 5: Enriched gene ontology terms for molecular function from active ingredients of PTH.

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).

Figure 6: Enriched KEGG pathways of potential targets from active ingredients of PTH.

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.

Figure 7: Component-target-pathway network. Red circles represent components, blue ellipses represent targets, and green squares represent pathways.

Table 2: The degree value of 16 compounds in PTH.

No.
Chemical compound
Degree

1
Sodium taurodeoxylate
48

2
Taurochenodeoxycholic acid
46

3
Glycocholic acid hydrate
43

4
Glycodeoxycholic acid hydrate
41

5
Ginsenoside Re
40

6
Ginsenoside Rg1
40

7
Notoginsenoside R1
40

8
Chenodeoxycholic acid
39

9
Deoxycholic acid
38

10
Astragaloside A
34

11
Cholic acid
37

12
Ginsenoside Rb1
37

13
Ginsenoside Rd
37

14
Muscone
37

15
Taurocholic acid
34

16
Taurine
0

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).

Figure 8: PPI network of target proteins.

Figure 9: Key subnetwork with core genes.

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.

Table 3: Molecular docking of core genes and active components.

Core gene
Active ingredient
Combined energy

ESR1
Glycocholic acid hydrate
−6.4

AR
Chenodeoxycholic acid
−7.3

ALB
Chenodeoxycholic acid
−9.5

NOTCH1
Ginsenoside Rb1
−7.4

ERBB2
Taurochenodeoxycholic acid
−7.4

IGF1R
Sodium taurodeoxylate
−7.1

(a)

(b)

(c)

(d)

(e)

(f)

Figure 10: Molecular docking charts of core genes and active ingredients. (a) ESR1 and glycocholic acid hydrate. (b) AR and chenodeoxycholic acid. (c) ALB and chenodeoxycholic acid. (d) IGF1 and sodium taurodeoxylate. (e) ERBB2 and taurochenodeoxycholic acid. (f) NOTCH1 and ginsenoside Rb1.

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.

References

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.

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.

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.

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.

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.

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.

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.

Pharmacopoeia Commission, Pharmacopoeia of the People’s Republic of China, Chinese Medical Science and Technology Press, Beijing, China, 2015.

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.

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.

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.

J. Y. Zhao, D. X. Chen, S. Lin et al., “Hepatoprotective activity of Pien Tze Huang Gan Bao in CCl₄-induced chronic liver injury models,” International Journal of Clinical and Experimental Medicine, vol. 12, pp. 11324–11332, 2019.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

H. Song, N. Sun, L. Lin et al., “Splicing factor PRPF6 upregulates oncogenic androgen receptor signaling pathway in hepatocellular carcinoma,” Cancer Science, 2020.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

X. Meng, D. A. Franklin, J. Dong et al., “MDM2-p53 pathway in hepatocellular carcinoma,” Cancer Research, vol. 74, pp. 7161–7167, 2014.

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.

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.

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.

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.

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.

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.

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.

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.