- Open Access
Network pharmacology based research into the effect and mechanism of Yinchenhao Decoction against Cholangiocarcinoma
Chinese Medicine volume 16, Article number: 13 (2021)
Cholangiocarcinoma refers to an epithelial cell malignancy with poor prognosis. Yinchenhao decoction (YCHD) showed positive effects on cancers, and associations between YCHD and cholangiocarcinoma remain unclear. This study aimed to screen out the effective active components of Yinchenhao decoction (YCHD) using network pharmacology, estimate their potential targets, screen out the pathways, as well as delve into the potential mechanisms on treating cholangiocarcinoma.
By the traditional Chinese medicine system pharmacology database and analysis platform (TCMSP) as well as literature review, the major active components and their corresponding targets were estimated and screened out. Using the software Cytoscape 3.6.0, a visual network was established using the active components of YCHD and the targets of cholangiocarcinoma. Based on STRING online database, the protein interaction network of vital targets was built and analyzed. With the Database for Annotation, Visualization, and Integrated Discovery (DAVID) server, the gene ontology (GO) biological processes and the Kyoto encyclopedia of genes and genomes (KEGG) signaling pathways of the targets enrichment were performed. The AutoDock Vina was used to perform molecular docking and calculate the binding affinity. The PyMOL software was utilized to visualize the docking results of active compounds and protein targets. In vivo experiment, the IC50 values and apoptosis rate in PI-A cells were detected using CCK-8 kit and Cell Cycle Detection Kit. The predicted targets were verified by the real-time PCR and western blot methods.
32 effective active components with anti-tumor effects of YCHD were sifted in total, covering 209 targets, 96 of which were associated with cancer. Quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol were identified as the vital active compounds, and AKT1, IL6, MAPK1, TP53 as well as VEGFA were considered as the major targets. The molecular docking revealed that these active compounds and targets showed good binding interactions. These 96 putative targets exerted therapeutic effects on cancer by regulating signaling pathways (e.g., hepatitis B, the MAPK signaling pathway, the PI3K-Akt signaling pathway, and MicroRNAs in cancer). Our in vivo experimental results confirmed that YCHD showed therapeutic effects on cholangiocarcinoma by decreasing IC50 values, down-regulating apoptosis rate of cholangiocarcinoma cells, and lowering protein expressions.
As predicted by network pharmacology strategy and validated by the experimental results, YCHD exerts anti-tumor effectsthrough multiple components, targets, and pathways, thereby providing novel ideas and clues for the development of preparations and the treatment of cholangiocarcinoma.
Cholangiocarcinoma refers to an uncommon biliary adenocarcinoma originating from epithelial cells of the biliary tract. It falls into two subtypes based on anatomical location, namely, intrahepatic cholangiocarcinoma and extrahepatic cholangiocarcinoma . The morbidity and mortality of intrahepatic cholangiocarcinoma have been rising in recent years, which is not consistent with extrahepatic cholangiocarcinoma [2, 3].
The occurrence of cholangiocarcinoma is likely to be associated with a range of risk factors (e.g., cirrhosis, hepatitis B and C infection, obesity, as well as diabetes mellitus), yet the etiology remains unclear [4,5,6]. Existing treatments for cholangiocarcinoma consist of surgery, chemotherapy and targeted radiation. For early stage cholangiocarcinoma patients, surgical treatment is most commonly required for tumor resection, thereby affecting patients’ postoperative recovery with the five-year survival rate (20–40%), and the surgery is not recommended for patients with advanced cholangiocarcinoma [7, 8]. Currently, Regimens based on 5-fluorouracil have been proved to slightly prolong survival time, and cisplatin combined with gemcitabine treatments have achieved the similarresults in comparison with gemcitabine alone [9, 10],whereaswestern medicineis likely to cause adverse reactions . Thus, traditional Chinese medicine, i.e., a complementary and alternative approach, is considered in the treatment of cholangiocarcinoma.
Yinchenhao decoction (YCHD) consists of three herbs, namely, Herba Artemisiae Scopariae (Yinchenhao, YCH), Gardeniae Fructus (Zhizi, ZZ) and Radix Rhei et Rhizoma (Dahuang, DH). YCHD is closely correlated with liver and gallbladderin accordance with the traditional Chinese medicine theory, which can eliminate dampness and heat of liver and gallbladder. Clinically, it has been widely adopted to treat cholestasis, liver disorders and metabolic diseases [12, 13]. Previous researches have also reported that YCHD had positive effects on various cancers including pancreatic cancer and ascites hepatoma [14, 15], however, the associations between YCHD and cholangiocarcinoma haven’t been thoroughly studied and require further researches. Covering numerous chemical compounds acting on multiple targets, the mechanisms of YCHD remain unclear. Network pharmacology refers to a promising methodology integrating pharmacology, molecular biology, electronic technology and bioinformatics to form network relationship among active ingredients of Chinese formulas, relevant targets, pathways, as well as diseases .
The present study was designed to delve into the mechanisms of YCHD on cholangiocarcinoma using network pharmacology methods, as an attempt to be referenced for subsequent pharmacological studies and clinical treatments of cholangiocarcinoma.The flowchart of YCHD in treating cholangiocarcinoma was provided in Fig. 1.
Materials and Methods
Constructing database of candidate compounds
All constituents of YCHD were achieved from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (http://lsp.nwu.edu.cn/tcmsp.php, TCMSP), Traditional Chinese Medicines Integrated Database (http://126.96.36.199:8000/tcmid/, TCMID), and Traditional Chinese Medicine Database@Taiwan (http://tcm.cmu.edu.tw/) [17,18,19]. Oral bioavailability (OB) and drug-likeness (DL) indices recommended by TCMSP were employed to verify druggability of each candidate. OB was referred to as the extent and rate to which the active drug ingredient or active moiety from the drug product was absorbed and available at the site of drug action . High OB appeared to be more probably a drug-like ingredient. DL index was adopted to assess whether the compounds were chemically suitable . Compounds that overcome both OB and drug-like phase screening would be recognized as candidate compounds more possibly. As TCMSP suggests, the molecules with OB ≥ 30% and DL ≥ 0.18 were preserved to display relatively better pharmacologically and then screened out as candidate compounds for subsequent analysis. As a result, thirty-two compounds of YCHD were identified.
Mining cholangiocarcinoma associated targets and target-genes
Protein targets associated with cholangiocarcinoma were provided by the GeneCard databases with “cholangiocarcinoma” as the keywords. All the targets were only limited to “homo sapiens”. Subsequently, protein names of all targets were switched to corresponding gene names in the website of UniProt (https://www.uniprot.org/) or PharmMapper (http://lilab.ecust.edu.cn/pharmmapper/). Furthermore, the same procedure was also performed to extract relevant target genes of herbs of YCHD.
Conducting PPI network
Given that the search of protein interactions and the interaction network is the critical procedure for gaining the insights into cellular organization, bioprocess, and functions, associated targets were input into STRING (Version 10.5, https://string-db.org/) to delve into protein-protein interactions. The network nodes and edges denote proteins and protein-protein associations, respectively. Two PPI interactive networks were built and then visualized by Cytoscape software (version 3.6.0), including estimated YCHD components and cholangiocarcinoma related targets. After merging these two networks as a candidate network following the intersection of PPI data, topological features were analyzed to sift a core PPI network.
Gene Ontology (GO) and pathway enrichment analysis
The DAVID web server (Functional Annotation Result Summary, https://david.ncifcrf.gov/summary.jsp) was adopted to conduct GO enrichment analysis for the candidate target protein obtained after network merging. Subsequently, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted to explore biological pathways where relevant proteins were covered. A P value ≤ 0.05 was considered significant, and enriched GO terms were identified by the hypergeometric test. A bubble chart was plotted via the OmicShare platform, a free online platform to conduct data analysis (http://www.omicshare.com/tools).
Validation of compound‐target interaction
The crystal structures of hub protein targets were obtained from the Protein Data Bank (PDB, https://www.rcsb.org/). The three-dimensional structures of candidate active compounds were downloaded from the PubChem (https://pubchem.ncbi.nlm.nih.gov/), an open archive for chemical information. The downloaded active compounds and hub protein targets were converted to the pdbqt format via the AutoDockTools (version 1.5.6, http://autodock.scripps.edu/). The AutoDock Vina (http://vina.scripps.edu/) was used to perform molecular docking and calculate the binding affinity. The docking results of active compounds and protein targets were visualized with PyMOL software (version 2.2, https://pymol.org/2/).
Cell cultures and cell viability measurements
CCK-8 kit was used to evaluate the effects of quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol on the KKU-M213 cell lines. The KKU-M213 cell lines seeded into 96-wall plates at a density of 8 × 103 cell/100 ul were treated with various concentrations of quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol for 24, 48, and 72 hours, respectively. Then the cells were incubated with 10 µl CCK-8 solution for 90 min. A microculture plate reader was used to measure the optical density at 450 nm. The SPSS 20.0 software was adopted to calculate 50% inhibitory concentration (IC50) value.
The PI-A cells in logarithmic growth of each group containing quercetin and kaempferol treatment were collected and seeded into 96-wall plates at a density of 1 × 106 cell/2 m, and incubated overnight at 37 °C with 5% CO2. The treated cells were washed with cold PBS solution, incubated with Annexin-FITC, and then placed in the water bath for 30 min. The Cell Cycle Detection Kit was applied to detect the cell circle.
Real‐time PCR assay
After an appropriate amount of cells in each group were collected, the supernatant was decanted and then washed with PBS. The TRIzol solution was added to extract total RNA of cells. The ultraviolet spectrophotometer was used to assess the RNA purity and concentration. AKT, P53, MAPK and VEGFA mRNA expression reversed transcription to cDNA. The above cDNA was used as template and amplified according to the real-time PCR instruction. The 20 µL PCR reaction system included 10 µL SYBR Green Mix, 1 µL each for upstream and downstream primer sequences, 2 µL template cDNA, and 2 µL DNase and RNase-free water. They were pre-denaturated at 94 °C for 5 min, denaturated at 94 °C for 30 s, annealed at 60 °C for 30 s, extended at 72 °C for 30 s, totally circulating 40 times. The β-actin acting as internal reference, the relative expression of each target mRNA was calculated by ABI Prism® SDS 2.0.3 using the 2−ΔΔCt method.
Western blot assay
The cells of each group were lysed with RIPA lysate for 30 min, and then transferred to the centrifuge tube. After centrifuging at 12,000 r/min for 10 min, the supernatant was extracted. The quantitative protein concentration was detected by BCA Protein Assay Kit. After separated by SDS-PAGE, 50 µg protein samples were transferred onto the PVDF membrane, sealed with 5% skim milk powder at room temperature for 1 h, and then washed with PBST solution. Rabbit anti AKT, P53, MAPK, VEGFA and IL-6 monoclonal antibodies (1:1000) were added respectively to incubate overnight at 4 °C, then the membrane was washed again, the corresponding secondary antibody were added, and the ECL kit was used to stain. The gray value of each imaging protein band was analyzed by the Gel Imaging System, and compare the relative expression change of each group protein with β-actin as internal reference.
The data were expressed as mean ± SD. The results were analyzed using GraphPad Prism 6.02 and SPSS 20.0 software. Student’s t-tests were developed to compare the between-group quantitative data, and p < 0.05 revealed significant difference.
Compound-compound target network analysis
32 active compounds were screened out in total following the OB ≥ 30% and DL ≥ 0.18 criteria. Table 1 lists the 32 active compounds. With the top 5 degrees, quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol were considered as the vital active compounds. Figure 2 shows that the compound-compound target network consisted of 244 nodes (namely, 3 herbs, 32 active compounds and 209 compound targets) and 640 edges. According to this network, targets in the interior circle exhibited more interactions than those in the exterior circle, and considerable targets were regulated by multiple compounds. For instance, Prostaglandin G/H synthase 2 was modulated by a range of compounds (e.g., crocetin, ammidin and isorhamnetin). It is therefore speculated that active compounds of YCHD might impact multiple targets to effectively treat cholangiocarcinoma. In this network, relationships between active compounds and compound targets as well as potential pharmacological effects of YCHD were visually illustrated.
Identification of intersection target genes
The network scores of genes were mapped from the mentioned calculation to the targets of cholangiocarcinoma disease genes, as well as target genes of herbs of YCHD, respectively, expressing the relationship between the three herbs of YCHD and cholangiocarcinoma disease. 96 target genes were identified affected by cholangiocarcinoma and regulated by YCH, ZZ and DH (Fig. 3; Table 2). As revealed from the results, ZZ were more critical to treat cholangiocarcinoma.
Compound‐compound target‐cholangiocarcinoma Target-Other Human Proteins’ PPI Network analysis
To identify the most highly connected nodes from others, the compound-compound target-cholangiocarcinoma target-other human proteins’ PPI network analysis was conducted. This network is illustrated in Fig. 4, covering 96 nodes and 1496 edges. YCHD exerted its therapeutic effects on cholangiocarcinoma through multiple protein targets. The nodes with top 5 degrees, including AKT1, IL6, MAPK1, TP53, and VEGFA, refer to the major targets in treating cholangiocarcinoma.
GO and pathway enrichment analysis
To clarify the function of estimated protein targets, the GO biological process enrichment analysis was conducted. The top 20 noticeably enriched GO terms are listed in Fig. 5; Table 3. As suggested from the results, the targets of YCHD displayed tight relations to the major biological process, which included positive regulation of transcription from RNA polymerase II promoter, negative regulation of apoptosis process, and positive regulation of transcription, DNA-templated. Thus, compound targets of YCHD exhibited similar functions to the corresponding genes of cholangiocarcinoma disease.
The KEGG pathway enrichment analysis was conducted using DAVID webserver. The target-pathway was built to delve into the mechanisms of potential targets acting on their corresponding signal pathways (Fig. 6; Table 4). YCHD was reported integrating multiple signaling pathways on cancers, immune system, infectious diseases, etc. Furthermore, YCHD probably exerted the therapeutic effects on cholangiocarcinoma by regulating signaling pathways, which included hepatitis B, the MAPK signaling pathway, the PI3K-Akt signaling pathway, and MicroRNAs in cancer.
Herb‐compound target‐cholangiocarcinoma Network analysis
To illustrate the relationship between three herbs of YCHD and their corresponding compound targets and cholangiocarcinoma targets, the herb-compound target-cholangiocarcinoma network was built. Based on this network, 80 nodes(3 herbs, 27 compounds, 46 targets, and 4 pathways) and 247 edges were identified (Fig. 7).
Verification of compound-target interaction
The molecular docking was performed to further investigate interactions between vital active compounds and major targets. The binding affinity lower than − 5.0 kcal/mol indicates that the confirmations have good interactions . In this study, molecular docking results showed that the conformations of active compounds and major protein targets showed good binding interactions, and the interactions were also reliable. The results of binding affinity was shown in Table 5. The conformations of key active compounds and major hub targets were displayed in Fig. 8.
3.7 YCHD inhibited the proliferation of cholangiocarcinoma cells
Five active compounds including quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol were observed to inhibit cholangiocarcinoma cell proliferation in a concentration-dependent manner. The 72 h IC50 values of quercetin, kaempferol, beta-sitosterol, isorhamnetin, stigmasterol were 10.84, 11.54, 48.54, 42.56, and 23.48, respectively, revealing that these five active compounds for the KKU-M213 cell line showed a significant reduction in IC50 values over time. The cell viability curves were shown in Fig. 9.
Quercetin and kaempferol induced apoptosis of cholangiocarcinoma cells
Higher doses of quercetin and kaempferol increased the apoptosis rate in PI-A cells compared with the control group, indicating that cholangiocarcinoma cells treated with quercetin and kaempferol may induce apoptosis (Fig. 10).
Effect of quercetin and kaempferol on AKT, P53, MAPK and VEGFA mRNA exprssion levels in cholangiocarcinoma cells
Several key targets including AKT, P53, MAPK and VEGFA mRNA expression levels were validated by PCR. Quercetin and kaempferol were found to decrease levels of AKT and VEGFA, and increase P53 and MAPK levels. The results were depicted in Fig. 11.
Quercetin combined with kaempferol regulated the protein expression of AKT, P53, MAPK, VEGFA and IL-6 in cell cycle control in cholangiocarcinoma cells
The targets including AKT, P53, MAPK, VEGFA and IL-6 were measured by western blot analysis. As shown in Fig. 12, quercetin and kaempferol decreased the expression of AKT, VEGFA and IL-6, and increased the expression of P53 and MAPK.
Cholangiocarcinoma refers to an epithelial cell malignancy with poor prognosis. Existing therapeutic strategies for patients with cholangiocarcinoma consist of surgery, liver transplantation and chemotherapy. Surgery is a prioritized treatment for all subtypes of cholangiocarcinoma, while lymph nodes and vascular structures should be considered. It was elucidated that for perihilar cholangiocarcinoma patients who have undergone either surgical treatment or liver transplantation, their 5-year survival rates were relatively low . Besides, chemotherapy might cause adverse effects and lower patients’ quality of life [11, 23].
Accordingly, traditional Chinese medicine, an effective and safe complementary and alternative therapy, is considered to be able to treat cholangiocarcinoma. YCHD is a notable traditional Chinese medicine prescription, which consisting of three Chinese herbs, namely, YCH, ZZ and DH; it has been adopted to clear heat, eliminate dampness and remove jaundice . As demonstrated by pharmacological researches, YCHD can modulate inflammatory as well as immune response, ameliorate liver function, and regulate multiple pathways in cancer [13, 25]. YCHD exerts therapeutic effects on liver disorders and metabolic diseases, whereas its mechanisms on cholangiocarcinoma remain unclear. In the present study, network pharmacology approach was adopted to describe the relationship among active compounds, compound targets and signaling pathways, and experimental methods were employed, thereby revealing the potential mechanisms of YCHD.
According to the results achieved here, 32 active compounds in YCHD with 209 compound targets were identified, suggesting that YCHD exerted its pharmacological effects on treating cholangiocarcinoma via multiple targets. Quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol were identified as the vital active compounds with top five degrees. As for quercetin, existing studies have indicated that quercetin is capable of inhibiting cell proliferation of cancer cell lines and regulating cancer metabolisms by modulating the PI3K-Akt-mTOR signaling and MAPK/ERK1/2 pathways [26, 27]. A previous experimental research has proved that quercetin could inhibit inflammatory and tumorigenesis processes and suppress growth and migration of cholangiocarcinoma cell lines . Kaempferol is reported to inhibit the growth and metastasis of cholangiocarcinoma through suppressing the PI3K-AKT pathway and their downstream proteins . As to beta-Sitosterol, it has anticancer properties associated with cell cycle, proliferation, apoptosis, etc. . Isorhamnetin, a flavonoid metabolite, has been demonstrated to be able to produce anti-tumor effects by inhibiting the expression of NF-κB . With regard to stigmasterol, it’s a major phytosterol in herbal plants, which can downregulate levels of inflammatory cytokines involving TNF-α,repress VEGF signaling, produce anti-angiogenic effects, thus inhibit cholangiocarcinoma growth in mice . On the whole, it was speculated that YCHD is a multicomponent formula with multitarget therapeutic effects. Associations between these active compounds and cholangiocarcinoma are supposed to be deeply investigated.
In this study, AKT1, IL6, MAPK1, TP53,and VEGFA were identified as the 5 hub protein targets related to cholangiocarcinoma. Quercetin and kaempferol were proved to decrease the protein expression of AKT, VEGFA and IL-6, and increase P53 and MAPK levels.Previous study has reported that combined mTOR and AKT inhibition may significantly inhibit the tumor growth and proliferation of cholangiocarcinoma cell lines . Similarly, the levels of AKT and IL6 can be decreased by genistein, which is associated with inhibitory effects on growth of cholangiocarcinoma cells . Regarding MAPK1, it can be activated by osteopontin, a phosphorylatedglycoprotein involved in various human cancers, and implicated in activating the MEK/MAPK pathway, therefore promote the growth and metastasis of intrahepatic cholangiocarcinoma . TP53 loss may drive the reprogramming of hepatocytes to biliary cells, which may be related to the formation of intrahepatic cholangiocarcinoma . VEGFA plays an important part in tumor angiogenesis. An experimental research has verified that interactions of LOXL2 with GATA6 can induce the expression of VEGFA, which may promote angiogenesis in cholangiocarcinoma subcutaneous tumorsand tumor growth .
The GO analysis has demonstrated that YCHD is associated with the major biological process (e.g., positive regulation of transcription from RNA polymerase II promoter, negative regulation of apoptosis process, and positive regulation of transcription, DNA-templated). Regarding the RNA polymerase II promoter, a transcriptional regulatory element, was served as a sequence initiating transcription and regulating epigenetics [38, 39]. The RNA polymerase II binding to the tumor suppressor CDC73 contributed to transcriptional repression of oncogenes, conversely, oncogene overexpression could down-regulate CDC73 in tumors and lead to cell proliferation . As for apoptosis, a physiological cell death, can regulate development of organisms, remove cells no longer available or differentiated ones incompatible with the body and maintain organizational homeostasis . Tumor will be induced once cell apoptosis is out of control. Caused by many risk factors, chronicinflammation and cholestasis are driving forces in cholangiocarcinoma development . The accumulation of bile acids from cholestasis lead to pH reduction, increase in apoptosis and activation of ERK1/2,Akt and NF-κB pathways . A previous research  has confirmed that cholangiocarcinoma cell apoptosis was induced by THZ1 inhibiting the synthesis of antiapoptotic protein, and THZ1 could down-regulate the transcriprion of RNA polymerase II promoter in cancer cells and produce anti-tumour effects. In addition, the transcription of DNA-templated possibly related to the anti-tumor property of THZ1 . The existing studies haven’t reported the relationship between DNA-templated transcription and cholangiocarcinoma, and this may become the further research direction. Therefore, we could speculate that the effects of YCHD on cholangiocarcinoma may relate to the above-mentioned biological processes.
As revealed from the KEGG pathway analysis, YCHD produced therapeutic effects on cholangiocarcinoma by regulating pathways (e.g., hepatitis B, the MAPK signaling pathway, the PI3K-Akt signaling pathway, and MicroRNAs in cancer). Recent epidemiological studies reported hepatitis B virus infection as a risk factor of intrahepatic cholangiocarcinoma [45, 46]. Intrahepatic cholangiocarcinoma patients with current and past hepatitis B virus infection were reported to display a better prognosis as compared with those without a history of hepatitis B virus infection .
Experimental results in this study confirmed that the active compounds of YCHD inhibited cell proliferation, induced apoptosis of cholangiocarcinoma cells, thus delayed the progression of cholangiocarcinoma. The MAPK signaling pathway participated in a variety of biological processes such as inflammatory response, cell differentiation, cell apoptosis and tumor invasiveness and so on . A previous study confirmed that the MAPK signaling pathway was activated in intrahepatic cholangiocarcinoma cells by TRIM44, a protein involved in several kinds of cancers, to inhibit cell apoptosis and promote tumour invasion and metastasis . For PI3K-Akt signaling pathway, it has been proved to regulate cell growth and proliferation and be critical to cancers [50, 51]. As one of the most intensively explored signaling pathways in tumorigenesis, PI3K-Akt signaling pathway may participate in the initiation, maintenance and metastasis of cholangiocarcinoma . Existing research has elucidated that theinhibitory of transcription factor 21 mediated by PI3K-Akt signaling pathway could inhibit the progression of cholangiocarcinoma . With respect to microRNAs, a group of small RNAs regulating genes expression, are modulators to suppress or progress tumor . According to one previous research, microRNA-329 had an inhibitory effect on the expression of PTTG1 and inactivated the MAPK signaling pathway to inhibit the cholangiocarcinoma cell proliferation, induce cell cycle arrest as well as promote cell apoptosis, thus prevented the progression of cholangiocarcinoma . Accordingly, it could be assumed that YCHD regulated hepatitis B, the MAPK signaling pathway, the PI3K-Akt signaling pathway, as well as MicroRNAs in cancer, so it could treat cholangiocarcinoma.
The molecular docking validated that vital active compounds and major targets showed good binding interactions. Except for regulation of quercetin on the PI3K-AKT and MAPK pathways described above, quercetin has the ability of inhibiting cell migration and angiogenesis mediated by downregulating VEGFA levels in glioblastoma cells , whereas no study has focused on the relationship of quercetin and VEGFA in treating cholangiocarcinoma. Kaempferol can treat cholangiocarcinoma through the PI3K-AKT pathway . Moreover,kaempferol lowers the levels of TP53, which is involved in attenuating apoptosis mediated by Cisplatin . In terms of beta-sitosterol, it is reported to induce MAPK phosphorylation, downregulate the PI3K/Akt,thus promote cell apoptosis and cell death . Results from KEGG pathway analysis indicated that the action mechanisms of YCHD for cholangiocarcinoma were implicated in hub targets of AKT1 and MAPK1. We could speculate that YCHD exerts therapeutic effects on cholangiocarcinoma through these active compounds, target genes and signaling pathways.
Most of the recent studies have applied the network pharmacology approach to screen active compounds, describe interactions between active compounds and relevant targets, and predict action mechanisms of various diseases [59, 60]. In addition to the network pharmacology methods, this study employed experimental methods to further explore the predictive mechanisms of YCHD for cholangiocarcinoma at molecular and cellular levels. However, there are some limitations in this study. First, some critical targets and active compounds may be ignored due to the incomplete information of databases. Second, multiple signaling pathways of YCHD acting on cholangiocarcinoma were predicted using the network pharmacology strategy, but the contribution of each pathway hasn’t been detected. Third, our current experiment and network pharmacology results provide directions for subsequent researches, whereas mechanisms still need fully investigated in future researches.
In this study, quercetin, kaempferol, beta-sitosterol, isorhamnetin, and stigmasterol were identified as the vital active compounds, and AKT1, IL6, MAPK1, TP53 as well as VEGFA were considered as the major targets. The molecular docking revealed that these active compounds and major targets showed good binding interactions. YCHD may treat cholangiocarcinoma via signaling pathways including hepatitis B, the MAPK signaling pathway, the PI3K-Akt signaling pathway, and MicroRNAs in cancer. Experimental researches provided evidence that YCHD showed therapeutic effects on cholangiocarcinoma by regulating related target protein, inhibiting cell proliferation, and increasing cell apoptosis rate. This study demonstrated potential pharmacological mechanisms of YCHD acting on cholangiocarcinoma; it can be referenced for clinical application of YCHD.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Traditional Chinese Medicine System Pharmacology Database and Analysis Platform
Database for Annotation, Visualization, and Integrated Discovery
Kyoto Encyclopedia of Genes and Genomes
Protein Data Bank
Razumilava N, Gregory J. Gores cholangiocarcinoma. Lancet. 2014;383:2168–79.
Yao D, Kunam V, Li X. A review of the clinical diagnosis and therapy of cholangiocarcinoma. J Int Med Res. 2014;42:3–16.
Patel T. Worldwide trends in mortality from biliary tract malignancies. BMC Cancer. 2002;2:10.
Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet. 2014;383:2168–79.
Palmer WC, Patel T. Are common factors involved in the pathogenesis of primary liver cancers? A meta-analysis of risk factors for intrahepatic cholangiocarcinoma. J Hepatol. 2012;57:69–76.
Charbel H, Al-Kawas F. Cholangiocarcinoma: epidemiology, risk factors, pathogenesis, and diagnosis. Curr Gastroenterol Rep. 2011;13:182–7.
Ruzzenente A, Conci S, Valdegamberi A, et al. Role of surgery in the treatment of intrahepatic cholangiocarcinoma. Eur Rev Med Pharmacol Sci. 2015;19:2892–900.
Huang J, Jiang B, Yang Y, et al. Influencing factors for the prognosis of patients with early-stage intrahepatic cholangiocarcinoma after radical resection. J Clin Hepatol. 2018;34:2157–63.
Lim K, Han S, Oh D, et al. Outcome of infusional 5-fluorouracil, doxorubicin, and mitomycin-C (iFAM) chemotherapy and analysis of prognostic factors in patients with refractory advanced biliary tract cancer. Oncology. 2012;83:57–66.
Valle J, Wasan H, Palmer D, et al. Cisplatin plus gemcitabine versus gemcitabine for biliary tract cancer. N Engl J Med. 2010;362:1273–81.
Liu S, Zhou L, An L. Implementation of comprehensive rehabilitation therapy in postoperative care of patients with cholangiocarcinoma and its impact on patients’ quality of life. Exp Ther Med. 2019;17:2703–7.
Chen Z, Ma X, Zhao Y, et al. Yinchenhao decoction in the treatment of cholestasis: a systematic review and meta-analysis. J Ethnopharmacol. 2015;168:208–16.
Jiang S, Hu X, Liu P. Immunomodulation and liver protection of Yinchenhao decoction against concanavalin A-induced chronic liver injury in mice. J Integr Med. 2015;13:262–8.
Qi J. Clinical Observation on Modified Yinchenhao Decoction Combined with Tegafur Gimeracil Oteracil Potassium Capsule in the Treatment of Advanced Pancreatic Cancer. Guangming J Chin Med. 2018;33:563–4.
Song Z, Song H, Song P. Clinical observation of modified Yinchenhao tang combined with cisplatin Intraperi-toneal perfusion for Ascites Hepatoma. J New Chin Med. 2018;50:139–42.
Li S, Zhang B, Jiang D, et al. Herb network construction and co-module analysis for uncovering the combination rule of traditional Chinese herbal formulae. BMC Bioinformatics. 2010;11(Suppl 11):1–12.
Ru J, Li P, Wang J,et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminformatics. 2014;6:13.
Huang L, Xie D, Yu Y, et al. TCMID 2.0: a comprehensive resource for TCM. Nucleic Acids Res. 2018;46:D1117–20.
Kramer C, Podewitz M, Ertl P, et al. Unique macrocycles in the Taiwan traditional Chinese medicine database. Planta Med. 2015;816:459–66.
Chao W, Lin B. Bioactivities of major constituents isolated from Angelica sinensis (Danggui). Chin Med. 2011;6:29.
Tao W, Xu X, Wang X, et al. Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease. J Ethnopharmacol. 2013;145:1–10.
Li X, Xu X, Wang J, et al. A system-level investigation into the mechanisms of Chinese Traditional Medicine: compound danshen formula for cardiovascular disease treatment. PLoS One. 2012;7:e43918.
He F, Wang M, Li K, et al. Efficacy analysis of hepatic arterial infusion in combination with intravenous gemcitabine chemotherapy for advanced intrahepatic cholangiocarcinoma. Chin J General Surg. 2018;27:219–24.
Guo Y, Li J, Mao T, et al. Effect of Combined Prescription of Linggui Zhugan Tang and Yinchenhao Tang on Nrf2/ARE Signaling Pathway in Rats with Non-alcoholic Steatohepatitis. Chin J Exp Trad Med Formul. 2017;23:108–13.
Huang J, Cheung F, Tan H, et al. Identification of the active compounds and significant pathways of yinchenhao decoction based on network pharmacology. Mol Med Rep. 2017;16:4583–92.
Reyes-Farias M, Carrasco-Pozo C. The anti-cancer effect of Quercetin: molecular implications in cancer metabolism. Int J Mol Sci. 2019;20:3177.
Erdogan S, Turkekul K, Dibirdik I, et al. Midkine downregulation increases the efficacy of quercetin on prostate cancer stem cell survival and migration through PI3K/AKT and MAPK/ERK pathway. Biomed Pharmacother. 2018;107:793–805.
Senggunprai L, Kukongviriyapan V, Prawan A, et al. Quercetin and EGCG exhibit chemopreventive effects in cholangiocarcinoma cells via suppression of JAK/STAT signaling pathway. Phytother Res. 2014;28:841–8.
Qin Y, Cui W, Yang X, et al. Kaempferol inhibits the growth and metastasis of cholangiocarcinoma in vitro and in vivo. Acta Biochim Biophys Sin. 2016;48:238–45.
Bin Sayeed MS, Ameen SS. Beta-Sitosterol. A Promising but Orphan Nutraceutical to Fight Against Cancer. Nutr Cancer. 2015;67:1214–20.
Manu KA, Shanmugam MK, Ramachandran L, et al. Isorhamnetin augments the anti-tumor effect of capecitabine through the negative regulation of NF-κB signaling cascade in gastric cancer. Cancer Lett. 2015;363:28–36.
Kangsamaksin T, Chaithongyot S, Wootthichairangsan C, et al. Lupeol and stigmasterol suppress tumor angiogenesis and inhibit cholangiocarcinoma growth in mice via downregulation of tumor necrosis factor-α. PLoS One. 2017;12:e0189628.
Ewald F, Grabinski N, Grottke A, et al. Combined targeting of AKT and mTOR using MK-2206 and RAD001 is synergistic in the treatment of cholangiocarcinoma. Int J Cancer. 2013;133:2065–76.
Tanjak P, Thiantanawat A, Watcharasit P, et al. Genistein reduces the activation of AKT and EGFR, and the production of IL6 in cholangiocarcinoma cells involving estrogen and estrogen receptors. Int J Oncol. 2018;53:177–88.
Zheng Y, Zhou C, Yu XX, et al. Osteopontin promotes metastasis of intrahepatic cholangiocarcinoma through recruiting MAPK1 and mediating Ser675 phosphorylation of β-Catenin. Cell Death Dis. 2018;9:179.
Hill MA, Alexander WB, Guo B, et al. Kras and Tp53 mutations cause cholangiocyte- and hepatocyte-derived cholangiocarcinoma. Cancer Res. 2018;78:4445–51.
Peng T, Deng X, Tian F, et al. The interaction of LOXL2 with GATA6 induces VEGFA expression and angiogenesis in cholangiocarcinoma. Int J Oncol. 2019;55:657–70.
Kadonaga J. Perspectives on the RNA polymerase II core promoter. Wiley Interdiscip Rev Dev Biol. 2012;1:40–51.
Chujan S, Suriyo T, Satayavivad J. Integrative in silico and in vitro transcriptomics analysis revealed gene expression changes and oncogenic features of normal cholangiocytes after chronic alcohol exposure. Int J Mol Sci. 2019;20:5987.
Rather M, Swamy S, Gopinath K, et al. Transcriptional repression of tumor suppressor CDC73, encoding an RNA polymerase II interactor, by Wilms tumor 1 protein (WT1) promotes cell proliferation: implication for cancer therapeutics. J Biol Chem. 2014;289:968–76.
Celli A, Que F. Dysregulation of apoptosis in the cholangiopathies and cholangiocarcinoma. Semin Liver Dis. 1998;18:177–85.
Andersen JB. Molecular pathogenesis of intrahepatic cholangiocarcinoma. J Hepatobiliary Pancreat Sci. 2015;22:101–13.
Labib PL, Goodchild G, Pereira SP. Molecular pathogenesis of cholangiocarcinoma. BMC Cancer. 2019;19:185.
Huang T, Ding X, Xu G, et al. CDK7 inhibitor THZ1 inhibits MCL1 synthesis and drives cholangiocarcinoma apoptosis in combination with BCL2/BCL-XL inhibitor ABT-263. Cell Death Dis. 2019;10:602.
Zhang H, Yang T, Wu M, et al. Intrahepatic cholangiocarcinoma: epidemiology, risk factors, diagnosis and surgical management. Cancer Lett. 2016;379:198–205.
Jeong S, Tong Y, Sha M, et al. Hepatitis B virus-associated intrahepatic cholangiocarcinoma: a malignancy of distinctive characteristics between hepatocellular carcinoma and intrahepatic cholangiocarcinoma. Oncotarget. 2017;8:17292–300.
Jiang B, Ge R, Sun L, et al. Clinical parameters predicting survival duration after hepatectomy for intrahepatic cholangiocarcinoma. Can J Gastroenterol. 2011;25:603–8.
Chen C, Nelson L, Ávila MA, et al. Mitogen-Activated Protein Kinases (MAPKs) and cholangiocarcinoma: the missing link. Cells. 2019;8:1172.
Peng R, Zhang PF, Zhang C, et al. Elevated TRIM44 promotes intrahepatic cholangiocarcinoma progression by inducing cell EMT via MAPK signaling. Cancer Med. 2018;7:796–808.
Li H, Hu J, Wu S, et al. Auranofin-mediated inhibition of PI3K/AKT/mTOR axis and anticancer activity in non-small cell lung cancer cells. Oncotarget. 2016;7:3548–58.
Burris HA. Overcoming acquired resistance to anticancer therapy: focus on the PI3K/AKT/mTOR pathway. Cancer Chemother Pharmacol. 2013;71:829–42.
Polivka J, Janku F. Molecular targets for cancer therapy in the PI3K/AKT/mTOR pathway. Pharmacol Ther. 2014;142:164–75.
Duan H, Li B, Zhaung X, et al. TCF21 inhibits tumor-associated angiogenesis and suppresses the growth of cholangiocarcinoma by targeting PI3K/Akt and ERK signaling. Am J PhysiolGastrointestLiver Physiol. 2019;316:G763–73.
Chang Y, Yin F, Fan GF, et al. Downregulation of miR-329-3p is associated with worse prognosis in patients withcervical cancer. Eur Rev Med Pharmacol Sci. 2017;21:4045–9.
Hu Z, Zheng C, Su H, et al. MicroRNA-329-mediated PTTG1 downregulation inactivates the MAPK signaling pathway to suppress cell proliferation and tumor growth in cholangiocarcinoma. J Cell Biochem. 2019;120:9964–78.
Liu Y, Tang ZG, Yang JQ, et al. Low concentration of quercetin antagonizes the invasion and angiogenesis of human glioblastoma U251 cells. Onco Targets Ther. 2017;10:4023–8.
Wang Z, Sun W, Sun X, et al. Kaempferol ameliorates Cisplatin induced nephrotoxicity by modulating oxidative stress, inflammation and apoptosis via ERK and NF-κB pathways. AMB Express. 2020;10:58.
Moon D, Lee K, Choi Y, et al. Beta-sitosterol-induced-apoptosis is mediated by the activation of ERK and the downregulation of Akt in MCA-102 murine fibrosarcoma cells. Int Immunopharmacol. 2007;7:1044–53.
Yu G, Zhang Y, Ren W, et al. Network pharmacology-based identification of key pharmacological pathways of Yin-Huang-Qing-Fei capsule acting on chronic bronchitis. Int J Chron Obstruct Pulmon Dis. 2016;12:85–94.
Lee W, Lee C, Kim Y, et al. The methodological trends of traditional herbal medicine employing network pharmacology. Biomolecules. 2019;9:362.
This study was supported by the grant from the National Science Foundation (NSFC) of China (NO.2019XZZX-ZL001).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests regarding the publication of this manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, Z., Lin, T., Liao, X. et al. Network pharmacology based research into the effect and mechanism of Yinchenhao Decoction against Cholangiocarcinoma. Chin Med 16, 13 (2021). https://doi.org/10.1186/s13020-021-00423-4
- Yinchenhao decoction
- Network pharmacology