Skip to main content

An advanced network pharmacology study to explore the novel molecular mechanism of Compound Kushen Injection for treating hepatocellular carcinoma by bioinformatics and experimental verification

Abstract

Background

Compound Kushen Injection (CKI) is a Chinese patent drug that exerts curative effects in the clinical treatment of hepatocellular carcinoma (HCC). This study aimed to explore the targets and potential pharmacological mechanisms of CKI in the treatment of HCC.

Methods

In this study, network pharmacology was used in combination with molecular biology experiments to predict and verify the molecular mechanism of CKI in the treatment of HCC. The constituents of CKI were identified by UHPLC-MS/MS and literature search. The targets corresponding to these compounds and the targets related to HCC were collected based on public databases. To screen out the potential hub targets of CKI in the treatment of HCC, a compound-HCC target network was constructed. The underlying pharmacological mechanism was explored through the subsequent enrichment analysis. Interactive Gene Expression Profiling Analysis and Kaplan-Meier plotter were used to examine the expression and prognostic value of hub genes. Furthermore, the effects of CKI on HCC were verified through molecular docking simulations and cell experiments in vitro.

Results

Network analysis revealed that BCHE, SRD5A2, EPHX2, ADH1C, ADH1A and CDK1 were the key targets of CKI in the treatment of HCC. Among them, only CDK1 was highly expressed in HCC tissues, while the other 5 targets were lowly expressed. Furthermore, the six hub genes were all closely related to the prognosis of HCC patients in survival analysis. Molecular docking revealed that there was an efficient binding potential between the constituents of CKI and BCHE. Experiments in vitro proved that CKI inhibited the proliferation of HepG2 cells and up-regulated SRD5A2 and ADH1A, while down-regulated CDK1 and EPHX2.

Conclusions

This study revealed and verified the targets of CKI on HCC based on network pharmacology and experiments and provided a scientific reference for further mechanism research.

Peer Review reports

Background

Globally, the incidence of liver cancer is increasing annually, and accounts for a large proportion of cancer cases and deaths [1]. According to the statistics in 2020, liver cancer ranks seventh in malignant tumor incidence and third in mortality [2]. Only in China, both the incidence and mortality of liver cancer are higher than the global average, and the survival rate of liver cancer is very low, which is a heavy burden on society and medical care worldwide [3,4,5,6]. Primary liver cancer includes different pathological subtypes, such as hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (ICC), and mixed hepatocarcinoma. HCC accounts for 75–85% of primary liver cancers and is the most common type [2, 7]. Due to its insidious onset, most HCC patients are at an advanced stage at the time of diagnosis. Unfortunately, local (chemoembolization) and surgical treatments are relatively disappointing in improving overall survival (OS) of advanced stage patients. At the same time, traditional chemotherapy methods also do not show promising efficacy in treating HCC due to significant side effects [8,9,10]. Therefore, it is essential to explore new prognostic markers for HCC and develop more effective drugs with less toxic effects.

In recent decades, traditional Chinese medicine (TCM) has been shown to play an important role in the treatment of tumors by inhibiting tumor proliferation, reducing tumor recurrence and metastasis, prolonging survival, and reducing the side effects of conventional treatments, tumors, etc. [11,12,13]. In particular, Chinese medicine injections, which have unique advantages in the treatment of tumors, are widely used in clinical practice because of their characteristics of the convenient application and rapid efficacy without irritating the gastrointestinal tract, as is the case with oral TCM [14,15,16].

Compound Kushen Injection (CKI), a kind of antineoplastic injection commonly used in clinical practice [17,18,19], consists of Kushen (Radix sophorae flavescentis) and Baituling (Rhizoma smilacis glabrae). Early studies have shown that CKI has a significant inhibitory effect on human HCC cells such as SMMC-7721, Hepa1–6, and LPC-H12 [19,20,21]. In addition, CKI significantly improved the clinical symptoms and quality of life of patients with advanced HCC. The combination of CKI and transcatheter arterial chemoembolization (TACE) has been reported to be effective in the treatment of advanced HCC and is therefore worthy of clinical application [22, 23]. The detailed information of CKI can be found in Supplementary file 1.

TCM preparations exert their special therapeutic effects by acting on the biological network in humans. Therefore, it is difficult to elucidate the specific mechanism of TCM in vivo [24, 25]. Until recently, network pharmacology provides a new method to promote the understanding of drug mechanisms [26, 27]. Network pharmacology constructs multicomponent and multitarget models to clarify better the complex interactions between genes, proteins, and metabolites related to diseases and drugs from a network perspective [28, 29]. In addition, molecular docking provides a relatively fast and economical alternative to standard experimental techniques that plays a vital role in the new drug development and discovery projects [30, 31]. AutoDock Vina software is often used to analyze the molecular interactions between protein and ligand [32]. To further explore and predict the molecular mechanism of CKI in the treatment of HCC and to identify targets related to the prognosis of HCC, network pharmacology was used in the present study. In addition, we further verified some of the potential targets of CKI through experiments in vitro. The flowchart of the current study is shown in Fig. 1.

Fig. 1
figure 1

Workflow of exploring the potential pharmacological mechanism of CKI in the treatment of HCC

Materials and methods

Identification of compounds and collection of potential targets for CKI

To identify the constituents of CKI, we conducted a qualitative analysis of CKI using UHPLC-MS/MS. Separations were performed on the Nexera LC-40 system (Shimadzu, Japan) using a Hypersil BDS (150 mm × 4. 6 mm, 5 μm). The mobile phase comprised 0.1% ammonia in water (A) and carbinol (B) at a flow rate of 0.5 ml/min and was eluted with gradient elution: 0–1 min (5–20% B), 1–30 min (20–80% B), 30–60 min (80–60% B). The column temperature was maintained at 25 °C and the injection volume was 5 μl. Analyses were performed with the QE-Orbitrap-MS (Thermo Fisher, USA) using an electrospray ionization (ESI) system. MS was operated in positive mode with the capillary temperature set at 350 °C. The spray voltage in negative mode was 3800 V, while the spray voltage in positive mode was 3200 V. The flow rates for the sheath gas and aux gas were 35arb and 15arb respectively. Three collision energies were used for MS2: low, medium and high. The positive ion mode was 30 eV, 40 eV, 50 eV, while the negative ion mode was 30 eV, 50 eV, 70 eV. The scan mode was set to full scan/ddMS2 and the scan range was 100–1200 Da. The full scan resolution was set to 70,000 FWHM while the resolution of MS2 was 17,500 FWHM. Finally, the retention time, MS fragmentation and UV spectra were used to identify the labeled reference compounds contained in the sample.

Next, we investigated the other compounds in CKI that had been reported in the literature to supplement our list of CKI-ingredients, and 23 compounds were selected for further study [20, 33, 34]. The canonical simplified molecular input line entry specification (SMILES) of 16 compounds were retrieved from the PubChem database [35] (https://pubchem.ncbi.nlm.nih.gov/) and exported. To search for the targets corresponding to these compounds, we imported the above data into the Search Tool for Interactions of Chemicals (STITCH) [36] (http://stitch.embl.de/), SuperPred [37] (http://prediction.charite.de/), SwissTargetPrediction [38] (http://www.swisstargetprediction.ch/), and Traditional Chinese Medicine Database and Analysis Platform (TCMSP) [39] (http://tcmspw.com/tcmsp.php).

Collection of HCC-related targets

Therapeutic Target Database (TTD), Gene Expression Omnibus (GEO), and Cancer Genome Atlas (TCGA) were the main sources of HCC-related targets. In the TTD [40] (http://bidd.nus.edu.sg/group/ttd/ttd.asp), HCC related targets were identified by searching for “hepatocellular carcinoma”. For the dataset in GEO [41] (http://www.ncbi.nlm.nih.gov/GEO/), the differentially expressed genes (DEG) in each microarray were first screened using the limma package [42]. Then the RobustRankAggreg package [43] was used to integrate genes that were generally identified as differentially expressed in different datasets. All datasets used for the analysis met the following two criteria: (1) The samples used in the dataset were tissues obtained from human HCC and corresponding adjacent or normal tissues. (2) The data set containing at least 40 samples. Data from the TCGA-LIHC project [44] were obtained from UCSC Xena (https://xenabrowser.net/datapages/) and analyzed in R using the edgeR package [45].

Network establishment

In the current study, three undirected networks were constructed: (1) a compound-putative target network, which included compounds of CKI and their corresponding targets; (2) a compound-HCC target network, which contained shared targets of CKI-ingredients and HCC; and (3) a drug-compound-target-pathway network. The above networks were visualized using Cytoscape 3.7.1 [46] (http://cytoscape.org/). Using the “Analyze Network” tool in Cytoscape3.7.1 software, the topological characteristics of each node in the HCC-related CKI compound-putative target network, namely Degree, Betweenness, and Closeness, were calculated to evaluate the importance of each node in the network.

Enrichment analysis

To illustrate the role of putative targets of CKI in the treatment of HCC in biological processes and signaling pathways, a Gene Ontology (GO) enrichment analysis and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using the Database for Annotation, Visualization, and Integrated Discovery [47] (DAVID, https://david.ncifcrf.gov/). These analyses were performed for the targets in the compound-HCC target network, and the results were visualized by the ggplot2 package (http://docs.ggplot2.org/current).

Survival analysis and correlation analysis of hub genes

Kaplan-Meier plotter (KM plotter, http://kmplot.com/analysis/), a database of clinical and gene expression data used to study the molecular basis of disease and identify biomarkers related to survival [48, 49]. The database contains disease-free survival and overall survival data based on the GEO, EGA (European Genome-Phenome Archive), and TCGA databases to calculate the hazard ratio (HR) with a 95% confidence interval and the P value of the log-rank test to assess the association between gene expression and survival [50].

Visualization of the results of expression level analysis and correlation analysis was carried out through Gene Expression Profiling Interactive Analysis [50] (GEPIA, http://GEPIA.cancerpku.cn/index.html), a web-based tool that allows the analysis of different tumor data in TCGA and Genotype Tissue Expression (GTEx). It contains 9736 tumor samples and 8587 normal tissue samples from 33 types of malignant tumors [51]. In addition, GEPIA provides customizable features such as differential expression analysis of tumors and normal. Using GEPIA, we have demonstrated the expression of key targets in HCC and normal tissues and the correlation between these targets.

Molecular docking

The molecular docking includes the following three steps:

  1. (1)

    The preparation of the receptors. The three-dimensional (3D) crystal structures of the key targets were extracted from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Database [52] (PDB, https://www. rcsb.org/), and the protein structures of these key targets were processed using AutoDock Tools (ADT) [53], including dehydration, hydrogenation, removal of small molecule ligands, and calculation of Gasteiger charge. The position of the active pocket was determined by the ligand in the crystal structure, and the structures were saved in *pdbqt format.

  2. (2)

    The preparation of ligands. The *SDF files of the 2D structure of all bioactives were downloaded from PubChem database (https://pubchem.ncbi.nlm.nih.gov/). Next, they were converted into the corresponding 3D structures with the help of ChemOffice software (https://www.chemdraw.com.cn/product.html), and their energy were minimized and saved as *mol2 files. Finally, the small ligand molecules in *mol2 format were exported to *pdbqt format through ADT.

  3. (3)

    Molecular docking simulation. Through AutoDock Vina [54], the molecular docking simulation of the bioactives and the target protein were carried out in turn, and the Affinity were extracted. The results were visualized and analyzed by PyMOL (http://www.PyMOL.org).

In vitro cell experimental verification

Reagents and cell line

CKI (Batch No. 20200329, total alkaloid concentration of 25 mg/ml) was provided by Shanxi Zhendong Pharmaceutical Co., Ltd. (China). Dulbecco modified eagle medium (DMEM), fetal bovine serum (FBS), penicillin-streptomycin, trypsin-EDTA, and phosphate buffered saline (PBS) were purchased from GIBCO (NY, USA). Cell Counting Kit-8 (CCK8) and BeyoClick™ EdU Cell Proliferation Kit with Alexa Fluor 488 were purchased from Beyotime (Beijing, China). RIPA lysis buffer, phenylmethylsulfonyl fluoride (PMSF) and BCA Protein Assay Kit were obtained from Biorigin (Beijing, China). Primary antibodies steroid 5-α-reductase 2 (SRD5A2) were purchased from Bioss (Beijing, China). Alcohol dehydrogenase 1A (ADH1A) was purchased from Solarbio (Beijing, China). Soluble epoxide hydrolase 2 (EPHX2), cyclin-dependent kinase 1 (CDK1), β-actin, and secondary antibody (anti-rabbit IgG and anti-mouse IgG) were procured from Proteintech (Wuhan, China). Nitrocellulose (NC) membranes and Immobilion Western Chemilum HRP Substrate were purchased from Merck Millipore (USA). Trizol reagent was obtained from Invitrogen (USA). SYBR Green Realtime PCR Master Mix was purchased from TOYOBO (Japan).

The human hepatocellular carcinoma cell line HepG2 was obtained from Procell (Wuhan, China) and cultured in DMEM supplemented with 10% FBS and 1% penicillin-streptomycin. Cells were maintained in a cell incubator at 37 °C and 5% CO2.

Cell proliferation assays

HepG2 cells in the logarithmic growth phase were processed, into a cell suspension, and seeded in 96-well plates at a density of 1 × 104 cells / well. After 24 h, cells were treated with different concentrations of CKI (0, 0.125, 0.25, 0.5, 1, 2, 4, 8, and 16 mg/ml) diluted with DMEM for 24 h, 48 h, and 72 h, respectively. Then, 10 μl of CCK8 reagent was added to each well to detect cell proliferation. The viability of the CKI treated cells was calculated using the optical density (OD) at 450 nm.

$$\mathrm{Cell}\ \mathrm{viability}=\left({\mathrm{OD}}_{\mathrm{CKI}}-{\mathrm{OD}}_{\mathrm{blank}}\right)/\left({\mathrm{OD}}_{\mathrm{control}}-{\mathrm{OD}}_{\mathrm{blank}}\right)\times 100\%$$

The EdU staining method was used to verify the effect of CKI on the proliferation of HepG2 cells. Cells were seeded in a confocal dish and treated with different concentrations of CKI (0, 1, 2, 4 mg/ml) after 24 h. According to the instructions of the EdU kit, the HepG2 cells treated with CKI for 48 h were fluorescently stained. Then the cells were observed with a fluorescence microscope (Olympus FV3000, Japan).

Real-time polymerase chain reaction (RT-PCR)

The Trizol reagent was used to extract total RNA from the cell samples under the same treatment conditions as in the WB experiment. The SYBR Green Realtime PCR Master Mix was used to determine mRNA expression in each sample according to the manufacturer’s protocols. Then the 2-ΔΔCT method was applied to calculate the relative expression level of the target gene. The primers used in the PCR were synthesized by Shenggong Biotech (Table 1).

Table 1 Primer sequences for RT-PCR

Western blot (WB)

HepG2 cells treated with CKI (0, 1, 2, 4 mg/ml) for 48 h were collected and total protein was extracted using RIPA lysis buffer supplement with 1% PMSF. The BCA method was used to calculate the protein concentration. After that, 5 × loading buffer was added to the protein samples and heated at 99 °C for 15 min to denature the protein. Then, the protein samples (20 μg/lane) were separated by 10% SDS-PAGE gel electrophoresis and transferred to NC membranes. Subsequently, the membranes were blocked in 5% skim milk for 1 h at room temperature. Afterward, the specific primary antibodies were added and incubated overnight at 4 °C, followed by incubation with the secondary antibodies at room temperature for 1 h. Finally, Immobilion Western Chemilum HRP Substrate was used to visualize the protein bands, and the gray values of the bands were quantified by ImageJ software (https://imagej.nih.gov/ij/).

Statistical analysis

Statistical analysis of experimental data was performed by GraphPad Prism8 (Inc. La Jolla, CA, United States). All data were expressed as mean ± standard deviation (SD) and were evaluated by Dunnett’s multiple comparisons test of one-way ANOVA. P-values of < 0.05 were considered to indicate a statistically significant difference. The calculation of the topological characteristics (Degree, Betweenness and Closeness) of each node in the network was performed using the “Analyze Network” tool in Cytoscape3.7.1.

Results

Construction of compound-HCC target network

Through UHPLC-MS/MS, 10 compounds were identified in CKI (Table 2). Through a literature search, additional 13 compounds were added in CKI. Finally, the 3D-chemical structures of 14 chemical components and 217 corresponding targets were obtained from public databases such as PubChem, STITCH, SuperPred and TCMSP. Table 2 shows the basic information of the 14 constituents in CKI. As shown in Fig. 2, the distributions of the mixed reference standards and the samples in the chromatogram were approximately the same. Based on the above data, the compound-putative target network (Fig. 3A) was established, including 231 nodes (14 compound nodes and 217 target nodes) and 555 edges.

Table 2 Basic information of the 14 compounds in CKI
Fig. 2
figure 2

The result of UHPLC-MS/MS. (A) CKI. (B) Mixed reference standards (MRS)

Fig. 3
figure 3

Construction of compound-HCC target network. A Compound-putative target network of CKI (red round rectangles represent compounds, and yellow octagons represent targets). B Identification of DEGs in TCGA and GEO. C Compound-HCC target network (red diamonds represent compounds, and blue hexagons represent targets of CKI for HCC treatment, the size of the node was proportional to its degree)

By consulting TTD, GEO, and TCGA, we identified HCC-related targets, of which 380 DEGs were obtained from GEO and 4276 DEGs were obtained from TCGA. After taking the intersection of the above two DEG sets, 343 DEGs were obtained (Fig. 3B). Integrating the 16 targets related to the treatment of HCC from TTD and the DEGs as above-mentioned, a total of 358 disease targets were obtained. Detailed information on the 358 disease targets that can be found in Supplementary file 2.

In order to explore the targets of CKI in the treatment of HCC, the CKI compound-putative target network and the targets related to HCC were merged, and the targets that did not overlap were removed. Then, 14 putative targets of CKI in the treatment of HCC were intuitively identified. As shown in Fig. 3C, the HCC-related CKI compound-putative target network consisted of 25 nodes (11 compound nodes and 14 target nodes) and 41 edges.

GO and KEGG enrichment analysis

To clarify the pharmacological mechanism of CKI in the treatment of HCC from the system level, GO and KEGG enrichment analyses were performed for the 14 key targets. In the results of GO analysis, 22 entries were selected based on P < 0.01, of which 12 were biological processes mainly related to drug metabolic process, drug response and xenobiotic metabolic process, etc., 9 were molecular functions mainly involving enzyme binding, chromatin binding and oxidoreductase activity, etc., and 1 was cell components, i.e. organelle membrane (Fig. 4A). KEGG enrichment analysis showed that the key genes were significantly enriched in 5 pathways, namely Drug metabolism-cytochrome P450 (hsa00982), Chemical carcinogenesis (hsa05204), Arachidonic acid metabolism (hsa00590), Retinol metabolism (hsa00830) and Metabolism of xenobiotics by cytochrome P450 (hsa00980) (Fig. 4B). Figure 4C shows two important metabolic pathways, namely drug metabolism-cytochrome P450 and arachidonic acid metabolism pathway.

Fig. 4
figure 4

Enrichment analysis of the compound-HCC target network. A The result of GO enrichment analysis. B The result of KEGG enrichment analysis. C Arachidonic acid metabolism pathway and the drug metabolism-cytochrome P450 pathway. D Drug-compound-target-pathway network (orange ellipses represent CKI, red diamonds represent compounds, blue hexagons represent shared targets of CKI and HCC, and green Vs represent pathways)

To systematically and holistically explain the mechanism of CKI in the treatment of HCC, Cytoscape software was used to construct a drug-compound-target-pathway network. As shown in Fig. 4D, there were a total of 31 nodes (1 CKI node, 11 compound nodes, 14 target nodes, and 5 pathway nodes) and 69 edges.

Survival analysis and correlation analysis of hub genes

According to the criterion of degree ≥2 (median degree), betweenness ≥0.0062 (average betweenness), closeness ≥0.3774 (average closeness), and edge count ≥2 (median edge count), six targets, including BCHE, SRD5A2, EPHX2, ADH1C, ADH1A, and CDK1, were identified as potential hub targets of CKI in the treatment of HCC. Figure 5A shows the expression levels of the six hub genes in tumor and normal tissues, respectively. BCHE, SRD5A2, EPHX2, ADH1C and ADH1A were lowly expressed in HCC, while CDK1 was highly expressed in HCC.

Fig. 5
figure 5

Expression and correlation analysis of the six hub genes of CKI in the treatment of HCC. A Expression levels of the six hub genes in the tumor and normal groups. B Correlation analysis between CDK1 and the other five important targets

Correlation analysis was performed for CDK1 and the above 5 low-expressed genes (BCHE, SRD5A2, EPHX2, ADH1C, and ADH1A). As shown in Fig. 5B, the expression of SRD5A2, EPHX2, ADH1C, and ADH1A in HCC showed a strong negative correlation with the expression of CDK1.

Survival analysis of these six hub targets was performed by KM plotter. The results showed that low-expressed SRD5A2 [HR = 0.59 (0.41–0.83), P = 0.0026], EPHX2 [HR = 0.51 (0.36–0.73), P = 0.00013], ADH1C [HR = 0.46 (0.3–0.7), P = 0.00018], and ADH1A [HR = 0.52 (0.36–0.74), P = 0.00028] and high-expressed CDK1 were all related to a poor prognosis in HCC patients. Only BCHE was not independent of the OS of HCC (P > 0.05) (Fig. 6). However, the level of BCHE was found to be an important survival factor for patients with prostate cancer [55].

Fig. 6
figure 6

Survival analysis of the six hub genes of CKI in the treatment of HCC

Molecular docking verification

Molecular docking simulation was used to analyze the interaction between the eight active compounds of CKI (9α-hydroxymatrine, isomatrine, lamprolobine, matrine, oxymatrine, sophocarpine, sophoranol, and sophoridine) and BCHE, the target with the highest degree. The crystal structure of BCHE (PDB ID: 5K5E) was retrieved from the PDB, and the 3D structures of the eight compounds were downloaded from the PubChem database. The docking results are listed in Table 3. As shown in Fig. 7, 9α-hydroxymatrine formed a hydrogen bond with Tyr332; isomatrine formed a hydrogen bond with Trp82; lamprolobine formed hydrogen bonds with Ser198 and Gly117; matrine formed a hydrogen bond with Thr120; oxymatrine formed a hydrogen bond with Tyr332; sophocarpine formed a hydrogen bond with Glu197, and sophoranol formed a hydrogen bond with Thr120.

Table 3 Information on the docking of BCHE-related compounds with BCHE
Fig. 7
figure 7

Detailed CKI ingredients-BCHE active site docking simulation

In vitro cell experimental verification

CKI inhibits the proliferation of hepatocellular carcinoma cells

To determine the effects of CKI on HepG2 cell viability, a CCK-8 assay was performed. As shown in Fig. 8A, CKI reduced the viability of HepG2 cells in a dose-dependent manner. In addition, the IC50 values of CKI after 24 h, 48 h, and 72 h were 1.616, 1.572 and 1.498 mg/ml, respectively. Also, the results of EdU showed f that the proliferation of HepG2 cells decreased with increasing dosage of CKI (Fig. 8B). Hoechst 33,342 staining showed that the number of HepG2 cells was significantly reduced under the intervention of 2 mg/ml and 4 mg/ml CKI for 48 h. Overall, these results suggest that CKI inhibits the growth and proliferation of HepG2 cells in a dose-dependent manner.

Fig. 8
figure 8

Effects of CKI on cell viability and proliferation of HepG2 cells. A Cell viability curve of HepG2 cells with CKI treated for 24 h, 48 h and 72 h. B The effect of CKI on the proliferation of HepG2 cells. Data are presented as mean ± SD of three independent experiments. (**** P < 0.0001)

Regulations of CKI on ADH1A, CDK1, EPHX2, and SRD5A2

The regulation of CKI on the potential targets (ADH1A, CDK1, EPHX2 and SRD5A2) in HepG2 cells were investigated by WB and RT-qPCR (Fig. 9A-B). The results showed that, CKI significantly increased the protein and mRNA expression of ADH1A and SRD5A2 in HepG2 cells compared with the control group. In contrast, the protein and mRNA expression of CDK1 and EPHX2 were significantly down-regulated by CKI. Notably, the prominent regulation of CKI on CDK1 was concentration-dependent.

Fig. 9
figure 9

Effects of CKI on the expressions of ADH1A, CDK1, EPHX2 and SRD5A2 in HepG2 cells. A ADH1A, CDK1, EPHX2 and SRD5A2 protein in HepG2 cells with CKI treated. B ADH1A, CDK1, EPHX2 and SRD5A2 mRNA in HepG2 cells with CKI treated. Data are presented as mean ± SD of three independent experiments. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001)

Discussion

HCC is the most common malignant tumor in the clinic and one of the leading causes of cancer-related deaths worldwide. In China, CKI is widely used in the treatment of cancer pain and is often combined with chemotherapy and radiotherapy in clinical practice [56]. Studies have shown that the combination of CKI with TACE has significant clinical efficacy in the treatment of HCC [22, 57]. In this study, the network pharmacology method was used in combination with molecular biology experiments to explore the potential targets and pharmacological mechanisms of CKI in the treatment of HCC.

Topological characteristics are important parameters for evaluating the importance of nodes in the network, including Degree, Betweenness and Closeness. “Degree” is defined as the number of edges associated with node i. “Betweenness” means the number of shortest paths between pairs of nodes passing through node i. “Closeness” represents the reciprocal of the sum of the distances between node i and the other nodes. The higher the above three node values, the greater the importance of the node in the network [58,59,60,61]. Based on the topological characteristics and edge counts of each node, we screened out 6 potential key targets (BCHE, SRD5A2, EPHX2, ADH1C, ADH1A, and CDK1) for CKI in the treatment of HCC. Among them, CDK1 was the only one highly expressed in HCC, and there was a strong negative correlation between its expression and the expression of the other five targets (BCHE, SRD5A2, EPHX2, ADH1C, and ADH1A).

Furthermore, the survival analysis results showed that high expression of CDK1 was closely related to the poor prognosis of HCC patients. In contrast, the high expression of SRD5A2, EPHX2, ADH1C and ADH1A indicated that HCC patients had a longer overall survival time. In addition, this study proved that CKI had a significant inhibitory effect on the proliferation of HepG2 cells.

Cyclin-dependent kinase 1 (CDK1) is a serine/threonine protein kinase and belongs to the cyclin-dependent kinase family (CDKs). The activated CDKs are necessary for cancer cell proliferation [62], and CDK1 is essential for the phase transition of G1/S and G2/M [63]. Zhao et al. found that CDK1 interacts with apoptin in HCC and is one of the key players in the activation of apoptin-induced tumor-specific apoptosis [64]. Clinical studies have shown that anti-CDK1 in HCC was one of the most effective ways to improve clinical efficacy [65]. In this study, we found that CKI caused a dose-dependent down-regulation of CDK1 in HepG2 cells. Mechanism studies have shown that matrine, one of the main ingredients of CKI, exhibits time-dependent inhibition of CDK1 expression in prostate cancer cells, resulting in o cell cycle arrest in the G0/G1 phase, thereby inhibiting cancer cell proliferation [66].

The alcohol dehydrogenase (ADH) family includes seven enzymes (ADH1–7) [67]. ADH1, including ADH1A, ADH1B, and ADH1C, are mainly expressed in the liver and play an important role in the conversion of ethanol to acetaldehyde, a kind of carcinogenic metabolite, especially in the excretory phase. Furthermore, these enzymes are essential for oral alcohol metabolism [68, 69]. Studies have shown that ADH1C polymorphisms were associated with various cancer risks such as gastric cancer and oral squamous cell carcinoma [70, 71]. It has been confirmed that low expression of ADH1C is associated with poor prognosis in patients with non-small cell lung cancer [67]. In HCC, both ADH1C and ADH1A have been shown to be downregulated [72, 73]. Inhibition of ADH1A expression would promote the occurrence and progression of HCC. However, high expression of ADH1A was related to a reasonable survival rate of HCC patients [74, 75]. Our experimental results showed that CKI increased the expression of ADH1A in HepG2 cells.

Butyrylcholinesterase (BCHE) is a glycoprotein synthesized by the liver and secreted into the bloodstream. It is a nonspecific cholinesterase and widely distributed in the nervous system, small intestine, and adipose tissue [76]. Recently, it has been reported to be a solid biochemical marker indicating liver organ damage. Functional destruction of the liver is associated with a decrease in plasma BCHE activity [77]. Therefore, BCHE has been found to be involved in the pathogenesis of chronic liver disease in patients with advanced HCC [78,79,80,81]. Besides, a lower level of BCHE in the serum is closely associated with the advanced stage and poor prognosis of various cancers, such as gastric cancer, renal cancer, bladder cancer, prostate cancer, and cervical cancer [76, 82,83,84,85]. Here, we carried out a molecular docking verification of the highest degree target BCHE and its related CKI-compounds. We found that all eight active compounds of CKI (9α-hydroxymatrine, isomatrine, lamprolobine, matrine, oxymatrine, sophocarpine, sophoranol, and sophoridine) had relatively high binding potential with the active site of BECE, indicating that there is likely an interaction between CKI and BCHE.

The EPHX2 gene encodes soluble epoxide hydrolase (sEH), a bifunctional enzyme of the epoxide hydrolase family, an enzyme that promotes increased apoptosis by inducing oxidative stress and inflammation [86]. Thomassen et al. found that sEH degrades biologically active epoxy fatty acids derived from arachidonic acid through the metabolism of P450 cytochromes (CYP) in the arachidonic acid pathway [87]. In addition, silencing of EPHX2 has been proved to reduce tumor cell viability, induce apoptosis, and inhibit androgen receptor signaling in prostate cancer [88]. For HCC, EPHX2 has been identified as a prognostic biomarker in previous studies [89]. In our research, CKI down-regulated the expression of EPHX2 in HCC.

SRD5A2 encodes steroid 5-α-reductase 2, an important enzyme in androgen metabolism. This gene has been found to expressed in cells sensitive to androgens, such as prostate, breast glands, and liver cells. It has been considered a risk factor for breast cancer and has been closely related to prostate cancer prognosis [90,91,92,93]. Previous studies have shown that SRD5A2 was down-regulated in tumors and negatively correlated with biochemical recurrence [91, 94]. In addition, epigenetic abnormal hypermethylation in the SRD5A2 promoter region can be used as an important molecular target for the detection of HCC [95], and the analysis of methylated SRD5A2 may help to accurately diagnose HCC, especially the early diagnosis of HCC [96, 97]. Our WB and PCR experimental results both showed that the expression level of SRD5A2 was significantly up-regulated in HepG2 cells after CKI treatment.

To further elucidate the multiple mechanisms of CKI in HCC from a systemic perspective, GO and KEGG enrichment analyses were performed. The results of GO analysis showed that CKI may be involved in drug response, metabolic drug process, and enzyme binding to affect HCC. In the metabolic pathway enrichment analysis, the potential targets of CKI were significantly enriched in drug metabolism-cytochrome P450, arachidonic acid metabolism and other metabolic pathways. Most of the drugs are metabolized in the liver by drug-metabolizing enzymes, and the primary drug-metabolizing enzymes are severely dysregulated in HCC, leading to chemotherapy failur e[98]. Moreover, dysfunction of the drug metabolism-cytochrome P450 pathway (hsa00982) has been reported to induce drug resistance or adverse reactions during chemotherapy in cancer by interrupting drug metabolism and promoting drug excretion [99]. Arachidonic acid metabolism has been proved to significantly impact the occurrence and development of various malignant diseases [100]. It has been reported that activation of arachidonic acid metabolism induces liver inflammation [101]. Studies have shown that berberine may induce apoptosis of HCC cells by inhibiting the metabolic pathway of arachidonic acid [102]. Therefore, we speculated that CKI may play an important therapeutic role in HCC by regulating drug metabolism-cytochrome P450 and arachidonic acid metabolic pathway.

Conclusions

In this study, we systematically investigated the potential targets and pharmacological mechanisms of CKI in the treatment of HCC through network pharmacology analysis combined with enrichment analysis, survival analysis and other experimental methods. The network analysis showed that BCHE, SRD5A2, EPHX2, ADH1C, ADH1A and CDK1 were the key targets of CKI in treating HCC. The survival analysis results showed that SRD5A2, EPHX2, ADH1C, ADH1A and CDK1 were closely related to the prognosis of HCC patients. GO and KEGG enrichment analysis revealed that CKI could exert its therapeutic effect by regulating drug metabolism-cytochrome P450 and arachidonic acid metabolism. In addition, based on the cell experiments in vitro, we confirmed that CKI had a significant inhibitory effect on the proliferation of HepG2 cells. The results of WB and PCR experiments suggest that the anti-HCC effect of CKI may be related to the down-regulation of CDK1 and EPHX2 and the up-regulation of SRD5A2 and ADH1A. In future experiments, we would investigate the roles and functions of these key targets in more depth and details.

In conclusion, this study systematically analyzed the pharmacological mechanism of CKI in the treatment of HCC. Further, it confirmed the anti-cancer effects and potential targets of CKI in the clinical treatment of HCC.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ADH:

Alcohol dehydrogenase

ADH1A:

Alcohol dehydrogenase 1A

ADH1C:

Alcohol dehydrogenase 1C

ADT:

AutoDock Tools

BCHE:

Butylcholinesterase

CDK1:

Cyclin-dependent kinase 1

CDKs:

Cyclin-dependent kinase

CCK8:

Cell Counting Kit-8

CKI:

Compound Kushen Injection

DAVID:

Database for Annotation, Visualization, and Integrated Discovery

DEGs:

Differentially expressed genes

DMEM:

Dulbecco modified eagle medium

EGA:

European Genome-Phenome Archive

EPHX2:

Epoxide hydrolase 2

FBS:

Fetal bovine serum

GEO:

Gene Expression Omnibus

GEPIA:

Gene Expression Profiling Interactive Analysis

GO:

Gene Ontology

GTEx:

Genotype Tissue Expression

HCC:

Hepatocellular carcinoma

HR:

Hazard ratio

ICC:

Intrahepatic cholangiocarcinoma

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KM:

Kaplan-Meier

NC membranes:

Nitrocellulose membranes

OD:

Optical density

OS:

Overall survival

PBS:

Phosphate buffered saline

PCR:

Polymerase Chain Reaction

PDB:

Protein Database

RCSB:

Research Collaboratory for Structural Bioinformatics

SD:

Standard deviation

sEH:

Soluble epoxide hydrolase

SMILES:

Simplified molecular input line entry specification

SRD5A2:

Steroid 5-α-reductase 2

STITCH:

Search Tool for Interactions of Chemicals

TACE:

Transcatheter arterial chemoembolization

TCGA:

The Cancer Genome Atlas

TCM:

Traditional Chinese medicine

TCMSP:

Traditional Chinese Medicine Database and Analysis Platform.

TTD:

Therapeutic Target Database

WB:

Western blot

References

  1. Liu Z, Jiang Y, Yuan H, Fang Q, Cai N, Suo C, et al. The trends in incidence of primary liver cancer caused by specific etiologies: results from the global burden of disease study 2016 and implications for liver cancer prevention. J Hepatol. 2019;70(4):674–83.

    Article  PubMed  Google Scholar 

  2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  3. Torre LA, Siegel RL, Ward EM, Jemal A. Global Cancer Incidence and Mortality Rates and Trends--An Update. Cancer Epidemiol Biomark Prev. 2016;25(1):16–27.

    Article  Google Scholar 

  4. Liu Z, Mao X, Jiang Y, Cai N, Jin L, Zhang T, et al. Changing trends in the disease burden of primary liver cancer caused by specific etiologies in China. Cancer Med. 2019;8(12):5787–99.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.

    Article  PubMed  Google Scholar 

  6. Feng RM, Zong YN, Cao SM, Xu RH. Current cancer situation in China: good or bad news from the 2018 global Cancer statistics? Cancer Commun (Lond). 2019;39(1):22.

    Article  Google Scholar 

  7. Gelband H, Chen C, Chen W. Liver Cancer. In: Cancer: disease control priorities, 3rd. Volume 3, edn. Edited by Gelband H, Jha P, Sankaranarayanan R, Horton S. Washington, DC: The International Bank for Reconstruction and Development/The World Bank press; 2015.

    Chapter  Google Scholar 

  8. Ma L, Chua MS, Andrisani O, So S. Epigenetics in hepatocellular carcinoma: an update and future therapy perspectives. World J Gastroenterol. 2014;20(2):333–45.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Chen S, Cao Q, Wen W, Wang H. Targeted therapy for hepatocellular carcinoma: challenges and opportunities. Cancer Lett. 2019;460:1–9.

    Article  CAS  PubMed  Google Scholar 

  10. Eatrides J, Wang E, Kothari N, Kim R. Role of systemic therapy and future directions for hepatocellular carcinoma. Cancer Control. 2017;24(3):1073274817729243.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Nie J, Zhao C, Deng LI, Chen J, Yu B, Wu X, et al. Efficacy of traditional Chinese medicine in treating cancer. Biomed Rep. 2016;4(1):3–14.

    Article  CAS  PubMed  Google Scholar 

  12. Hong M, Wang N, Tan HY, Tsao SW, Feng Y. MicroRNAs and Chinese medicinal herbs: New possibilities in Cancer therapy. Cancers (Basel). 2015;7(3):1643–57.

    Article  CAS  Google Scholar 

  13. Xiang Y, Guo Z, Zhu P, Chen J, Huang Y. Traditional Chinese medicine as a cancer treatment: modern perspectives of ancient but advanced science. Cancer Med. 2019;8(5):1958–75.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Wang H, Hu H, Rong H, Zhao X. Effects of compound Kushen injection on pathology and angiogenesis of tumor tissues. Oncol Lett. 2019;17(2):2278–82.

    CAS  PubMed  Google Scholar 

  15. Wang Z, Qi F, Cui Y, Zhao L, Sun X, Tang W, et al. An update on Chinese herbal medicines as adjuvant treatment of anticancer therapeutics. Biosci Trends. 2018;12(3):220–39.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang D, Wu J, Duan X, Wang K, Ni M, Liu S, et al. Network Meta-analysis of Chinese herbal injections plus the FOLFOX regimen for the treatment of colorectal Cancer in China. Integr Cancer Ther. 2019;18:1534735419827098.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Yanju B, Yang L, Hua B, Hou W, Shi Z, Li W, et al. A systematic review and meta-analysis on the use of traditional Chinese medicine compound kushen injection for bone cancer pain. Support Care Cancer. 2014;22(3):825–36.

    Article  PubMed  Google Scholar 

  18. Chen H, Yao X, Li T, Lam CW, Zhang H, Wang J, et al. Compound Kushen injection plus platinum-based chemotherapy for stage IIIB/IV non-small cell lung cancer: a protocol for meta-analysis of randomized clinical trials following the PRISMA guidelines. Medicine (Baltimore). 2019;98(52):e18552.

    Article  Google Scholar 

  19. Yang Y, Sun M, Yao W, Wang F, Li X, Wang W, et al. Compound kushen injection relieves tumor-associated macrophage-mediated immunosuppression through TNFR1 and sensitizes hepatocellular carcinoma to sorafenib. J Immunother Cancer. 2020;8(1).

  20. Gao L, Wang KX, Zhou YZ, Fang JS, Qin XM, Du GH. Uncovering the anticancer mechanism of compound Kushen injection against HCC by integrating quantitative analysis, network analysis and experimental validation. Sci Rep. 2018;8(1):624.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Wang KX, Chen YP, Lu AP, Du GH, Qin XM, Guan DG, et al. A metabolic data-driven systems pharmacology strategy for decoding and validating the mechanism of compound Kushen injection against HCC. J Ethnopharmacol. 2021;274:114043.

    Article  CAS  PubMed  Google Scholar 

  22. Zhang Y, Hui F, Yang Y, Chu H, Qin X, Zhao M, et al. Can Kushen injection combined with TACE improve therapeutic efficacy and safety in patients with advanced HCC? A systematic review and network meta-analysis. Oncotarget. 2017;8(63):107258–72.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ma X, Li RS, Wang J, Huang YQ, Li PY, Wang J, et al. The therapeutic efficacy and safety of compound Kushen injection combined with Transarterial chemoembolization in Unresectable hepatocellular carcinoma: an update systematic review and Meta-analysis. Front Pharmacol. 2016;7:70.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zhang L, Wu C, Zhang Y, Liu F, Wang X, Zhao M, et al. Comparison of efficacy and toxicity of traditional Chinese medicine (TCM) herbal mixture LQ and conventional chemotherapy on lung cancer metastasis and survival in mouse models. PLoS One. 2014;9(10):e109814.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Xu HY, Liu ZM, Fu Y, Zhang YQ, Yu JJ, Guo FF, et al. Exploiture and application of an internet-based computation platform for integrative pharmacology of traditional Chinese medicine. China J Chin Mater Med. 2017;42(18):3633–8.

    Google Scholar 

  26. Li S, Zhang B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11(2):110–20.

    Article  PubMed  Google Scholar 

  27. Zhang R, Zhu X, Bai H, Ning K. Network pharmacology databases for traditional Chinese medicine: review and assessment. Front Pharmacol. 2019;10:123.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhang S, Shan L, Li Q, Wang X, Li S, Zhang Y, et al. Systematic analysis of the multiple bioactivities of green tea through a network pharmacology approach. Evid Based Complement Alternat Med. 2014;2014:512081.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhang Y, Mao X, Su J, Geng Y, Guo R, Tang S, et al. A network pharmacology-based strategy deciphers the underlying molecular mechanisms of Qixuehe capsule in the treatment of menstrual disorders. Chin Med. 2017;12:23.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Yuriev E, Holien J, Ramsland PA. Improvements, trends, and new ideas in molecular docking: 2012-2013 in review. J Mol Recognit. 2015;28(10):581–604.

    Article  CAS  PubMed  Google Scholar 

  31. Pinzi L, Rastelli G. Molecular docking: shifting paradigms in drug discovery. Int J Mol Sci. 2019;20(18).

  32. Di Muzio E, Toti D, Polticelli F. DockingApp: a user friendly interface for facilitated docking simulations with AutoDock Vina. J Comput Aided Mol Des. 2017;31(2):213–8.

    Article  PubMed  Google Scholar 

  33. Aung TN, Nourmohammadi S, Qu Z, Harata-Lee Y, Cui J, Shen HY, et al. Fractional deletion of compound Kushen injection indicates cytokine signaling pathways are critical for its perturbation of the cell cycle. Sci Rep. 2019;9(1):14200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang W, You RL, Qin WJ, Hai LN, Fang MJ, Huang GH, et al. Anti-tumor activities of active ingredients in compound Kushen injection. Acta Pharmacol Sin. 2015;36(6):676–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kim S, Thiessen PA, Bolton EE, Chen J, Fu G, Gindulyte A, et al. PubChem substance and compound databases. Nucleic Acids Res. 2016;44(D1):D1202–13.

    Article  CAS  PubMed  Google Scholar 

  36. Szklarczyk D, Santos A, von Mering C, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2016;44(D1):D380–4.

    Article  CAS  PubMed  Google Scholar 

  37. Nickel J, Gohlke BO, Erehman J, Banerjee P, Rong WW, Goede A, et al. SuperPred: update on drug classification and target prediction. Nucleic Acids Res. 2014;42(Web Server issue):W26–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Gfeller D, Grosdidier A, Wirth M, Daina A, Michielin O, Zoete V. SwissTargetPrediction: a web server for target prediction of bioactive small molecules. Nucleic Acids Res. 2014;42(Web Server issue):W32–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wang Y, Zhang S, Li F, Zhou Y, Zhang Y, Wang Z, et al. Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. 2020;48(D1):D1031–d1041.

    CAS  PubMed  Google Scholar 

  41. Clough E, Barrett T. The gene expression omnibus database. Methods Mol Biol. 2016;1418:93–110.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28(4):573–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Zhang Z, Li H, Jiang S, Li R, Li W, Chen H, et al. A survey and evaluation of web-based tools/databases for variant analysis of TCGA data. Brief Bioinform. 2019;20(4):1524–41.

    Article  CAS  PubMed  Google Scholar 

  45. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  46. Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.

    Article  CAS  PubMed  Google Scholar 

  47. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  PubMed  Google Scholar 

  48. Hou GX, Liu P, Yang J, Wen S. Mining expression and prognosis of topoisomerase isoforms in non-small-cell lung cancer by using Oncomine and Kaplan-Meier plotter. PLoS One. 2017;12(3):e0174515.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Lacny S, Wilson T, Clement F, Roberts DJ, Faris P, Ghali WA, et al. Kaplan-Meier survival analysis overestimates cumulative incidence of health-related events in competing risk settings: a meta-analysis. J Clin Epidemiol. 2018;93:25–35.

    Article  PubMed  Google Scholar 

  50. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–w102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Sun C, Yuan Q, Wu D, Meng X, Wang B. Identification of core genes and outcome in gastric cancer using bioinformatics analysis. Oncotarget. 2017;8(41):70271–80.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. 2009;30(16):2785–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Gu Y, Chow MJ, Kapoor A, Mei W, Jiang Y, Yan J, et al. Biphasic alteration of Butyrylcholinesterase (BChE) during prostate Cancer development. Transl Oncol. 2018;11(4):1012–22.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Guo YM, Huang YX, Shen HH, Sang XX, Ma X, Zhao YL, et al. Efficacy of compound Kushen injection in relieving Cancer-related pain: a systematic review and Meta-analysis. Evid Based Complement Alternat Med. 2015;2015:840742.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Cho WC, Chen HY. Transcatheter arterial chemoembolization combined with or without Chinese herbal therapy for hepatocellular carcinoma: meta-analysis. Expert Opin Investig Drugs. 2009;18(5):617–35.

    Article  CAS  PubMed  Google Scholar 

  58. Liu H, Zeng L, Yang K, Zhang G. A network pharmacology approach to explore the pharmacological mechanism of Xiaoyao powder on Anovulatory infertility. Evid Based Complement Alternat Med. 2016;2016:2960372.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Jeong H, Mason SP, Barabási AL, Oltvai ZN. Lethality and centrality in protein networks. Nature. 2001;411(6833):41–2.

    Article  CAS  PubMed  Google Scholar 

  60. Kivimäki I, Lebichot B, Saramäki J, Saerens M. Two betweenness centrality measures based on randomized shortest paths. Sci Rep. 2016;6:19668.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Sabidussi G. The centrality of a graph. Psychometrika. 1966;31(4):581–603.

    Article  CAS  PubMed  Google Scholar 

  62. Malumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm. Nat Rev Cancer. 2009;9(3):153–66.

    Article  CAS  PubMed  Google Scholar 

  63. Dorée M, Hunt T. From Cdc2 to Cdk1: when did the cell cycle kinase join its cyclin partner? J Cell Sci. 2002;115(Pt 12):2461–4.

    Article  PubMed  Google Scholar 

  64. Wu CX, Wang XQ, Chok SH, Man K, Tsang SHY, Chan ACY, et al. Blocking CDK1/PDK1/β-catenin signaling by CDK1 inhibitor RO3306 increased the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma. Theranostics. 2018;8(14):3737–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Zhao J, Han SX, Ma JL, Ying X, Liu P, Li J, et al. The role of CDK1 in apoptin-induced apoptosis in hepatocellular carcinoma cells. Oncol Rep. 2013;30(1):253–9.

    Article  CAS  PubMed  Google Scholar 

  66. Chang J, Hu S, Wang W, Li Y, Zhi W, Lu S, et al. Matrine inhibits prostate cancer via activation of the unfolded protein response/endoplasmic reticulum stress signaling and reversal of epithelial to mesenchymal transition. Mol Med Rep. 2018;18(1):945–57.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Wang P, Zhang L, Huang C, Huang P, Zhang J. Distinct prognostic values of alcohol dehydrogenase family members for non-small cell lung Cancer. Med Sci Monit. 2018;24:3578–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Wang HL, Zhou PY, Liu P, Zhang Y. ALDH2 and ADH1 genetic polymorphisms may contribute to the risk of gastric cancer: a meta-analysis. PLoS One. 2014;9(3):e88779.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Birley AJ, James MR, Dickson PA, Montgomery GW, Heath AC, Martin NG, et al. ADH single nucleotide polymorphism associations with alcohol metabolism in vivo. Hum Mol Genet. 2009;18(8):1533–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Duell EJ, Sala N, Travier N, Muñoz X, Boutron-Ruault MC, Clavel-Chapelon F, et al. Genetic variation in alcohol dehydrogenase (ADH1A, ADH1B, ADH1C, ADH7) and aldehyde dehydrogenase (ALDH2), alcohol consumption and gastric cancer risk in the European prospective investigation into Cancer and nutrition (EPIC) cohort. Carcinogenesis. 2012;33(2):361–7.

    Article  CAS  PubMed  Google Scholar 

  71. Brocic M, Supic G, Zeljic K, Jovic N, Kozomara R, Zagorac S, et al. Genetic polymorphisms of ADH1C and CYP2E1 and risk of oral squamous cell carcinoma. Otolaryngol Head Neck Surg. 2011;145(4):586–93.

    Article  PubMed  Google Scholar 

  72. Liu Y, Zhu X, Zhu J, Liao S, Tang Q, Liu K, et al. Identification of differential expression of genes in hepatocellular carcinoma by suppression subtractive hybridization combined cDNA microarray. Oncol Rep. 2007;18(4):943–51.

    PubMed  Google Scholar 

  73. Gaviria-Calle M, Duque-Jaramillo A, Aranzazu M, Di Filippo D, Montoya M, Roldán I, et al. Polymorphisms in alcohol dehydrogenase (ADH1) and cytochrome p450 2E1 (CYP2E1) genes in patients with cirrhosis and/or hepatocellular carcinoma. Biomedica. 2018;38(4):555–68.

    PubMed  Google Scholar 

  74. Zahid KR, Yao S, Khan ARR, Raza U, Gou D. mTOR/HDAC1 crosstalk mediated suppression of ADH1A and ALDH2 links alcohol metabolism to hepatocellular carcinoma onset and progression in silico. Front Oncol. 2019;9:1000.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Gao Q, Zhu H, Dong L, Shi W, Chen R, Song Z, et al. Integrated Proteogenomic characterization of HBV-related hepatocellular carcinoma. Cell. 2019;179(2):561–577.e522.

    Article  CAS  PubMed  Google Scholar 

  76. Poetsch N, Sturdza A, Aust S, Polterauer S, Grimm C, Schwameis R, et al. The value of pretreatment serum butyrylcholinesterase level as a novel prognostic biomarker in patients with cervical cancer treated with primary (chemo-)radiation therapy. Strahlenther Onkol. 2019;195(5):430–40.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Pavo N, Raderer M, Goliasch G, Wurm R, Strunk G, Cho A, et al. Subclinical involvement of the liver is associated with prognosis in treatment naïve cancer patients. Oncotarget. 2017;8(46):81250–60.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Santarpia L, Alfonsi L, Pasanisi F, De Caprio C, Scalfi L, Contaldo F. Predictive factors of survival in patients with peritoneal carcinomatosis on home parenteral nutrition. Nutrition. 2006;22(4):355–60.

    Article  PubMed  Google Scholar 

  79. Donadon M, Cimino M, Procopio F, Morenghi E, Montorsi M, Torzilli G. Potential role of cholinesterases to predict short-term outcome after hepatic resection for hepatocellular carcinoma. Updat Surg. 2013;65(1):11–8.

    Article  Google Scholar 

  80. Ohashi N, Tsuji N, Naito Y, Iwakura T, Isobe S, Ono M, et al. Relationship between urinary fractional excretion of sodium and life prognosis in liver cirrhosis patients. Hepatol Res. 2013;43(11):1156–62.

    Article  CAS  PubMed  Google Scholar 

  81. Takeda H, Nishikawa H, Iguchi E, Ohara Y, Sakamoto A, Hatamaru K, et al. Impact of pretreatment serum cholinesterase level in unresectable advanced hepatocellular carcinoma patients treated with sorafenib. Mol Clin Oncol. 2013;1(2):241–8.

    Article  PubMed  Google Scholar 

  82. Gu SZ, Zhao XH, Quan P, Li SB, Pan BR. Alterations of serum cholinesterase in patients with gastric cancer. World J Gastroenterol. 2005;11(29):4604–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Koie T, Ohyama C, Mikami J, Iwamura H, Fujita N, Sato T, et al. Preoperative butyrylcholinesterase level as an independent predictor of overall survival in clear cell renal cell carcinoma patients treated with nephrectomy. Sci World J. 2014;2014:948305.

    Article  Google Scholar 

  84. Koie T, Ohyama C, Yamamoto H, Hatakeyama S, Imai A, Yoneyama T, et al. Significance of preoperative butyrylcholinesterase as an independent predictor of survival in patients with muscle-invasive bladder cancer treated with radical cystectomy. Urol Oncol. 2014;32(6):820–5.

    Article  CAS  PubMed  Google Scholar 

  85. Koie T, Ohyama C, Hatakeyama S, Imai A, Yoneyama T, Hashimoto Y, et al. Significance of preoperative butyrylcholinesterase as an independent predictor of biochemical recurrence-free survival in patients with prostate cancer treated with radical prostatectomy. Int J Clin Oncol. 2016;21(2):379–83.

    Article  CAS  PubMed  Google Scholar 

  86. Roumeliotis S, Roumeliotis A, Stamou A, Panagoutsos S, Manolopoulos VG, Tsetsos F, et al. Association of rs11780592 polymorphism in the human soluble epoxide hydrolase gene (EPHX2) with oxidized LDL and mortality in patients with diabetic chronic kidney disease. Oxidative Med Cell Longev. 2021;2021:8817502.

    Article  Google Scholar 

  87. Thomassen M, Tan Q, Kruse TA. Gene expression meta-analysis identifies chromosomal regions and candidate genes involved in breast cancer metastasis. Breast Cancer Res Treat. 2009;113(2):239–49.

    Article  PubMed  Google Scholar 

  88. Vainio P, Gupta S, Ketola K, Mirtti T, Mpindi JP, Kohonen P, et al. Arachidonic acid pathway members PLA2G7, HPGD, EPHX2, and CYP4F8 identified as putative novel therapeutic targets in prostate cancer. Am J Pathol. 2011;178(2):525–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Zhang R, Ye J, Huang H, Du X. Mining featured biomarkers associated with vascular invasion in HCC by bioinformatics analysis with TCGA RNA sequencing data. Biomed Pharmacother. 2019;118:109274.

    Article  CAS  PubMed  Google Scholar 

  90. Wigley WC, Prihoda JS, Mowszowicz I, Mendonca BB, New MI, Wilson JD, et al. Natural mutagenesis study of the human steroid 5 alpha-reductase 2 isozyme. Biochemistry. 1994;33(5):1265–70.

    Article  CAS  PubMed  Google Scholar 

  91. Johnston WL, Catton CN, Swallow CJ. Unbiased data mining identifies cell cycle transcripts that predict non-indolent Gleason score 7 prostate cancer. BMC Urol. 2019;19(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  92. Li J, Coates RJ, Gwinn M, Khoury MJ. Steroid 5-{alpha}-reductase type 2 (SRD5a2) gene polymorphisms and risk of prostate cancer: a HuGE review. Am J Epidemiol. 2010;171(1):1–13.

    Article  PubMed  Google Scholar 

  93. Rubicz R, Zhao S, Wright JL, Coleman I, Grasso C, Geybels MS, et al. Gene expression panel predicts metastatic-lethal prostate cancer outcomes in men diagnosed with clinically localized prostate cancer. Mol Oncol. 2017;11(2):140–50.

    Article  CAS  PubMed  Google Scholar 

  94. Zhang D, Li Q, Qu HC, Yu T, Liu YR. Associations between the SRD5A2 gene V89L and TA repeat polymorphisms and breast cancer risk: a meta-analysis. Genet Mol Res. 2015;14(3):9004–12.

    Article  CAS  PubMed  Google Scholar 

  95. Yu SJ, Kim H, Min H, Sohn A, Cho YY, Yoo JJ, et al. Targeted proteomics predicts a sustained complete-response after Transarterial chemoembolization and clinical outcomes in patients with hepatocellular carcinoma: a prospective cohort study. J Proteome Res. 2017;16(3):1239–48.

    Article  CAS  PubMed  Google Scholar 

  96. Moribe T, Iizuka N, Miura T, Stark M, Tamatsukuri S, Ishitsuka H, et al. Identification of novel aberrant methylation of BASP1 and SRD5A2 for early diagnosis of hepatocellular carcinoma by genome-wide search. Int J Oncol. 2008;33(5):949–58.

    CAS  PubMed  Google Scholar 

  97. Tsunedomi R, Ogawa Y, Iizuka N, Sakamoto K, Tamesa T, Moribe T, et al. The assessment of methylated BASP1 and SRD5A2 levels in the detection of early hepatocellular carcinoma. Int J Oncol. 2010;36(1):205–12.

    CAS  PubMed  Google Scholar 

  98. Ul-Islam S, Ahmed MB, Shehzad A, Ul-Islam M, Lee YS. Failure of chemotherapy in hepatocellular carcinoma due to impaired and Dysregulated primary liver drug metabolizing enzymes and drug transport proteins: what to do? Curr Drug Metab. 2018;19(10):819–29.

    Article  CAS  PubMed  Google Scholar 

  99. Zhao L, Lei H, Shen L, Tang J, Wang Z, Bai W, et al. Prognosis genes in gastric adenocarcinoma identified by cross talk genes in disease-related pathways. Mol Med Rep. 2017;16(2):1232–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Yang P, Cartwright CA, Li J, Wen S, Prokhorova IN, Shureiqi I, et al. Arachidonic acid metabolism in human prostate cancer. Int J Oncol. 2012;41(4):1495–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Chiu AP, Tschida BR, Sham TT, Lo LH, Moriarity BS, Li XX, et al. HBx-K130M/V131I promotes liver Cancer in transgenic mice via AKT/FOXO1 signaling pathway and Arachidonic acid metabolism. Mol Cancer Res. 2019;17(7):1582–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  102. Yarla NS, Bishayee A, Sethi G, Reddanna P, Kalle AM, Dhananjaya BL, et al. Targeting arachidonic acid pathway by natural products for cancer prevention and therapy. Semin Cancer Biol. 2016;40-41:48–81.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This study was supported by the Young Scientists Training Program of Beijing University of Chinese Medicine and the National Nature Science Foundation of China [Grant numbers 81673829]. The funding body approved the study’s design, analysis and interpretation of data, and publication of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

SL and ZQM conceived and designed the study. SL, YYT, CW, ZHH and JQH performed the in vitro experiments. SL and SYG analysed and interpreted the research data. XKL, LMY and JYZ provided significant suggestions on the methodology. SL and ZQM majorly contributed to writing the manuscript. CGF, XJYL and WZ helped revise the manuscript. AS helped touch up the manuscript. XMZ, MMW, JRW were responsible for funding acquisition. All authors have read and approved the final draft.

Corresponding author

Correspondence to Jiarui Wu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Detailed information of Compound Kushen injection.

Additional file 2: Table S2.

Detailed information of 358 disease targets.

Additional file 3: Figure S1.

Original images of WB.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, S., Meng, Z., Tan, Y. et al. An advanced network pharmacology study to explore the novel molecular mechanism of Compound Kushen Injection for treating hepatocellular carcinoma by bioinformatics and experimental verification. BMC Complement Med Ther 22, 54 (2022). https://doi.org/10.1186/s12906-022-03530-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12906-022-03530-3

Keywords