Network pharmacology based research into the effect and mechanism of Yinchenhao Decoction against Cholangiocarcinoma

Background 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. Methods 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. Results 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. Conclusions 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.


Background
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 [1]. 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 [11]. 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 Conclusions: 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. Keywords: Yinchenhao decoction, Network pharmacology, Mechanism, Cholangiocarcinoma Fig. 1 The whole framework of this study based on network pharmacology for investigating pharmacological mechanisms of YCHD acting on cholangiocarcinoma pharmacology, molecular biology, electronic technology and bioinformatics to form network relationship among active ingredients of Chinese formulas, relevant targets, pathways, as well as diseases [16].
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.  [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 [20]. High OB appeared to be more probably a drug-like ingredient. DL index was adopted to assess whether the compounds were chemically suitable [21]. 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.unipr ot.org/) or PharmMapper (http://lilab .ecust .edu.cn/pharm mappe r/). 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 ://strin g-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 .ncifc rf.gov/summa ry.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.omics hare.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 ://pubch em.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://autod ock.scrip ps.edu/). The AutoDock Vina (http://vina.scrip ps.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 × 10 3 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.

Apoptosis detection
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 × 10 6 cell/2 m, and incubated overnight at 37 °C with 5% CO 2 . 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 predenaturated 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.

Statistical analysis
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 [22]. 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.

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

Discussion
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 [1]. 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 [24]. 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 [28]. Kaempferol is reported to inhibit the growth and metastasis of cholangiocarcinoma through suppressing the PI3K-AKT pathway and their downstream proteins [29]. As to beta-Sitosterol, it has anticancer properties associated with cell cycle, proliferation, apoptosis, etc. [30]. Isorhamnetin, a flavonoid metabolite, has been demonstrated to be able to produce anti-tumor effects by inhibiting the expression of NF-κB [31]. 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 [32]. 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 [33]. Similarly, the levels of AKT and IL6 can be decreased by genistein, which is associated with inhibitory effects on growth of cholangiocarcinoma cells [34]. Regarding , 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 [35]. TP53 loss may drive the reprogramming of hepatocytes to biliary cells, which may be related to the formation of intrahepatic cholangiocarcinoma [36]. 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 [37].
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 [40]. 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 [41]. 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 [42]. The accumulation of bile acids from cholestasis lead to pH reduction, increase in apoptosis and activation of ERK1/2,Akt and NF-κB pathways [43]. A previous research [44] 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 [44]. 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 [47].
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 [48]. 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 [49]. 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 [52]. Existing research has elucidated that theinhibitory of transcription factor 21 mediated by PI3K-Akt signaling pathway could inhibit the progression of cholangiocarcinoma [53]. With respect to microRNAs, a group of small RNAs regulating genes expression, are modulators to suppress or progress tumor [54]. According to one previous research, micro-RNA-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 [55]. Accordingly, it could be assumed that YCHD regulated hepatitis B, the MAPK signaling pathway, the PI3K-Akt Fig. 7 The herb-compound Target-cholangiocarcinoma network 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 [56], whereas no study has focused on the relationship of quercetin and VEGFA in treating cholangiocarcinoma. Kaempferol can treat cholangiocarcinoma through the PI3K-AKT pathway [29]. Moreover,kaempferol lowers the levels of TP53, which is involved in attenuating apoptosis mediated by Cisplatin [57]. In terms of beta-sitosterol, it is reported to induce MAPK phosphorylation, downregulate the PI3K/Akt,thus promote cell apoptosis and cell death [58]. 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   All data were expressed as mean ± SD 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.

Conclusions
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.  , d). All data were expressed as mean ± SD ( ** P < 0.01; *** P < 0.001)