Skip to main content
  • Research article
  • Open access
  • Published:

Study on the mechanisms of compound Kushen injection for the treatment of gastric cancer based on network pharmacology



As an effective prescription for gastric cancer (GC), Compound Kushen Injection (CKI) has been widely used even though few molecular mechanism analyses have been carried out.


In this study, we identified 16 active ingredients and 60 GC target proteins. Then, we established a compound-predicted target network and a GC target protein-protein interaction (PPI) network by Cytoscape 3.5.1 and systematically analyzed the potential targets of CKI for the treatment of GC. Finally, molecular docking was applied to verify the key targets. In addition, we analyzed the mechanism of action of the predicted targets by Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses.


The results showed that the potential targets, including CCND1, PIK3CA, AKT1, MAPK1, ERBB2, and MMP2, are the therapeutic targets of CKI for the treatment of GC. Functional enrichment analysis indicated that CKI has a therapeutic effect on GC by synergistically regulating some biological pathways, such as the cell cycle, pathways in cancer, the PI3K-AKT signaling pathway, the mTOR signaling pathway, and the FoxO signaling pathway. Moreover, molecular docking simulation indicated that the compounds had good binding activity to PIK3CA, AKT1, MAPK1, ERBB2, and MMP2 in vivo.


This research partially highlighted the molecular mechanism of CKI for the treatment of GC, which has great potential in the identification of the effective compounds in CKI and biomarkers to treat GC.

Peer Review reports


Gastric cancer (GC) is one of the most common malignant tumors, and its mortality rate ranks third in the world, making it the fifth most frequently diagnosed cancer according to the latest report of Global Cancer Statistics (GLOBOCAN 2018). The global incidence rate varies widely, with the highest incidence rates in East Asia and Eastern Europe, while North America, parts of Africa and Northern Europe have the lowest incidence [1]. Environmental factors, including Helicobacter pylori (H. pylori) infection, smoking, high salt intake, and other dietary factors, can increase the risk of GC [2]. GC is often diagnosed at an advanced stage, and surgery can possibly cure patients with GC; moreover, a multimodality approach can increase survival [3]. Combined in different forms and with different drugs, chemotherapy and radiotherapy are common methods for the treatment of GC, and chemotherapy shows the best supportive care in advanced GC [4]. However, the adverse effects of radiotherapy and chemotherapy frequently cause physical damage to patients, such as upper limb swelling on the affected side, decreased white blood cells and immunity, gastrointestinal reactions, bone marrow suppression, immune dysfunction and organ damage, causing tremendous pain to the patients [5, 6]. As an important part of complementary and alternative medicine, traditional Chinese medicine has been widely used in the clinical treatment of cancer in Asian countries, especially in China and Japan [7]. Compound Kushen Injection (CKI) consists of two herbs, Kushen (Radix Sophorae Flavescentis) and Baituling (Rhizoma Smilacis Glabrae). CKI is often used to treat many types of solid tumors, such as hepatic carcinoma, lung cancer and GC, and it is increasingly being used for cancer-related pain and has been used clinically in China for more than 15 years [8, 9]. Kushen has a long history of being used for the treatment of solid tumors, inflammation and other diseases, such as gynecological diseases [10]. Studies have found that Baituling can not only enhance the cytotoxicity of Kushen in CKI, but also activate immune-related pathways [11]. CKI mainly contains various anticancer ingredients such as matrine and oxymatrine, which can inhibit tumor cell growth, induce the apoptosis of cancer cells, and resist cancer cell metastasis, and is unsusceptible to multidrug resistance and protective against overactive immune responses in the body [12]. Research has shown that the matrine in CKI can inhibit cell proliferation and induce apoptosis in gastric carcinoma cells (SGC-7901) [13]. CKI can limit mouse sarcoma growth by reducing the phosphorylation of ERK and AKT kinases and BAD. CKI can also block TRPV1 signaling by inhibiting TRPV1-mediated capsaicin-induced ERK phosphorylation and reducing tumor-induced proinflammatory cytokine production, indirectly limiting cancer pain by reducing tumor growth [14]. The chemical composition of traditional Chinese medicine (TCM) compound is complex, and they are usually prepared by the synergy of multiple Chinese medicines and multiple targets to exert pharmacodynamic effects [15]. The network pharmacology “disease-gene-target-drug” research model is very consistent with the synergistic effect of the TCM “multicomponent, multichannel, multitarget” model, which more effectively reveals the relationship between TCM compound preparation and disease [16]. As a new field, network pharmacology provides a deeper insight into the basic mechanisms of TCM theory by using systems biology, highly integrated data analysis strategies and visualizations of the interpretations [17]. TCM network pharmacology has the role of discovering biologically active compounds and elucidating the mechanisms of action of TCM formulas, with the significance of accelerating the discovery of TCM and improving current drug discovery strategies [18]. Besides, molecular docking methods were used to assess the binding of components in CKI to key targets to confirm the drug’s effect on the target [19].

Previous studies have confirmed the effectiveness of CKI in the treatment of GC [20, 21] but did not explore the molecular biology. To better analyze and predict the molecular mechanism of CKI in the treatment of GC, this study adopted the network pharmacology and molecular docking method. The therapeutic target prediction of active ingredients and drug target network analysis were systematically carried out and may provide a reference for the basic experimental research and clinical rational application of CKI in the treatment of GC. The detailed workflow of the network pharmacology-based study of CKI is shown in Fig. 1.

Fig. 1
figure 1

Workflow for CKI in treating GC



Collection of target proteins related to GC

The Therapeutic Target Database (TTD) [22] ( is a database that allows users to search for targets, biomarkers, and drugs related to various disease conditions. Online Mendelian Inheritance in Man (OMIM) [23] ( is an online database of continuously updated human genes and genetic diseases. Pharmacogenomics Knowledge Base (PharmGKB) [24] ( is a database of high-quality pharmacogenomic data that provides information on genetic variation, pharmacokinetics, and pharmacodynamic pathways, as well as data on the effects of mutations on drug-related phenotypes. The disease gene search engine with evidence sentences (DigSEE) [25] ( database can directly query the links between diseases, genes and biological processes. DisGeNET [26] ( is one of the most comprehensive genetic disease association databases. The above databases were searched for “gastric cancer” as a keyword to collect proteins related to GC.

Collection of CKI active ingredients

Following a literature search [27], the 23 active ingredients contained in CKI were selected for research (Additional file 1: Table S1), and the three-dimensional structure data of 16 active ingredients were searched and derived from the PubChem database [28].

Prediction of potential targets of CKI

The 3D chemical structure data of 16 active ingredients were imported into the Search Tool for Interactions of Chemicals (STITCH), SuperPred, and SwissTargetPrediction databases for retrieval. STITCH shows the link between proteins and small molecules in biological processes, visualizing them in a network, and provides information such as the strength of their association; it can also be used to search for potential targets for compounds [29]. SuperPred uses chemical similarity principles to predict possible biological targets of drugs and demonstrates the interaction between compounds and predicted targets [30]. SwissTargetPrediction performs a combination of similarity measurements based on known 2D and 3D chemical structures to predict the corresponding potential bioactive targets [31]. The predicted multiple target information of the compounds and the obtained information were introduced into Cytoscape 3.5.1 ( to obtain a compound-predicted target network map. Cytoscape is bioinformatics analysis software that visualizes biological pathways and intermolecular interaction networks. It provides a basic set of data integration, analysis and visualization capabilities for complex network analyses [32].

Construction of related protein interaction networks

Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) 10.5 ( is a database of known and predicted protein interactions that contains direct and indirect associations of proteins. The number and quality of interacting proteins can be set according to their confidence settings, and it has a score for each protein interaction. The higher the score is, the higher the confidence of the protein interaction [33]. The proteins related to GC collected by the databases and the proteins predicted according to the 3D structure of the compound were input into the STRING 10.5 database, and the species selection “Homo sapiens” was selected as the confidence data with a scoring value greater than 0.7. Then, the data were introduced into Cytoscape to construct a protein-protein interaction (PPI) network related to GC and a PPI network related to the predicted target of the compound.

Network merge

The active components of CKI and the predicted target PPI networks were combined with the GC-PPI networks. The overlapping proteins in the two networks are likely to be potential targets for the treatment of GC by the active ingredients of CKI. A potential key target network for the CKI treatment of GC was constructed by Cytoscape, and the potential targets of the network were systematically analyzed.

GO functional and KEGG pathway enrichment analysis

To illustrate the role of the key potential targets in gene function and signaling pathways, this study used the Database for Annotation, Visualization and Integrated Discovery (DAVID) v 6.8 [34] ( to perform GO functional enrichment analyses and KEGG pathway enrichment analyses on the potential key targets in the CKI network for the treatment of GC.

Molecular docking simulation

The molecular docking of systemsDock [19] ( is a web server based on the docking program designed by AutoDock VINA with a high-precision scoring function. It has the characteristics of high precision and fast speed. The docking score for systemsDock is the negative logarithm of the experimental dissociation/inhibition constant (pKd/pKi), which directly indicates the binding strength.


PPI network of GC targets

A total of 60 target proteins (after removal of duplicates) related to GC were retrieved from the TTD, OMIM, PharmGKB, DigSEE, and DisGeNET databases (Additional file 2: Table S2). The PPI network was constructed in the STRING 10.5 database and included 154 nodes (60 disease-related nodes and 94 other related human target nodes), which constitute 1484 protein interactions (Fig. 2).

Fig. 2
figure 2

PPI network of GC targets (red nodes are disease-related targets and yellow nodes are 94 other related human targets)

Compound-predicted target network

Basic information on the 16 active ingredients in CKI is shown in Table 1. The active compound-predicted target network (Fig. 3) consists of 190 nodes (16 compound nodes and 174 target nodes) that constitute 331 active compound-predicted target linkages (Additional file 3: Table S3).

Table 1 Information on the active ingredients of CKI
Fig. 3
figure 3

Compound-predicted target network of CKI (16 compound nodes are red, and 174 target nodes are blue-green)

Potential key target network for the treatment of GC with CKI

The CKI active compound-predicted target PPI network was combined with the GC-PPI network to remove proteins that did not intersect, and the potential targets of CKI for the treatment of GC were intuitively obtained. Essentially, the node connections in the network identify the same targets in the CKI active component prediction target PPI network as those found in the GC-PPI network, that is, the potential targets of the CKI active ingredient for treating GC. Using Cytoscape software to obtain the potential target degree, the greater the degree is, the greater the possibility that the CKI active ingredients will act on GC through the targets. The potential target network of CKI for the treatment of GC is shown in the Fig. 4 and includes 43 targets (6 targets for disease and active ingredients, 4 disease-related targets, 7 active-related targets, and 26 other related human targets) and 358 potential links. The targets of the GC-CKI active component potential target network have an average degree of 16.65. Among the targets, the degree of 22 targets is greater than or equal to 17, indicating that they are likely to be vital potential targets for the treatment of GC with the active ingredients in CKI.

Fig. 4
figure 4

Potential key target network of CKI for the treatment of GC (red diamonds represent the same targets for GC and CKI components. Pink squares indicate the targets of the disease. Yellow hexagons represent the targets of the CKI active ingredients. Blue circles represent PPI-related human targets)

Molecular docking verification

Fourteen of the potential targets were directly related to the CKI active components, and these 14 targets were imported into the Protein Data Bank (PDB) ( [35] database to find their 3D structure. The 3D structure of CCND1 could not be obtained, but the remaining 13 protein structures were available. The 3D structures of the remaining 13 targets were introduced into systemsDock for molecular docking verification. Most of the docking scores were larger than 5, which showed that the proteins possessed good binding activity [36]. The docking results are shown in Fig. 5. The docking details are shown in Table 2. The protein structures were showed by Pymol [37].

Fig. 5
figure 5

Detailed target-compound interactions of the docking simulation

Table 2 Information on molecular docking

GO functional and KEGG pathway enrichment analysis

Using the MCODE algorithm in Cytoscape for cluster analysis [38], two modules were obtained after clustering. Module 1 contained 28 nodes with 157 edges, and Module 2 contained 7 nodes with 20 edges (Fig. 6). Using the DAVID platform for GO functional enrichment analysis, the roles of the 43 target proteins involved in the Module 1 and Module 2 networks in gene function were studied. For the respective modules, 22 and 13 GO entries were determined based on the false discovery rate (FDR) and P values (FDR < 0.01 and P < 0.01). KEGG function enrichment analysis was performed by using the DAVID platform, and the roles of the 43 target proteins involved in the Module 1 and Module 2 networks in signaling pathways were studied. A total of 76 KEGG entries were determined based on the FDR and P values (FDR < 0.01 and P < 0.01) (Fig. 7). The overall mechanism prediction of the effect of CKI on the treatment of gastric cancer is shown in Fig. 8.

Fig. 6
figure 6

Module analysis infographic

Fig. 7
figure 7

Graph of the enrichment analysis results (A. Graph of Module 1 GO analysis results; B. Graph of Module 2 GO analysis results; and C. Chart of KEGG enrichment analysis results)

Fig. 8
figure 8

Drug-ingredient-target-pathway network. This network intuitively illustrates the interactions of the chemical components, predictive targets and pathways of CKI for the treatment of GC. (The yellow node represents CKI. The pink nodes represent the three active components in CKI. The red nodes represent the corresponding key targets. The purple nodes represent the target KEGG enrichment pathways)


The incidence of GC is extremely high in many countries worldwide, and the incidence of GC is much higher in China than in many other countries [39]. CKI has been clinically used in China for over 15 years to treat various types of solid tumors, including GC. This study constructed a PPI network related to GC, a CKI active component-predicted target network, and a potential target network of CKI for the treatment of GC, and then systematically analyzed the mechanism of action of CKI in the treatment of GC. The molecular docking method was used to verify the correspondence between the target and the component, and the good docking score reflects the effectiveness of the component.

The potential target network of CKI for the treatment of GC contained a total of 43 targets, sorted into CCND1, PIK3CA, PIK3CG, PIK3CD, AKT1, PIK3CB, CDKN1B, EGF, CDKN1A, CDK2, MAPK1, JUN, CDK4, PTEN, MTOR, RB1, MAPK8, FoxO1, CCNB1, FOS, PCNA, EGFR, FoxO3, NFKB1, CDK6, RPTOR, IKBKB, CCNA1, RPS6KB1, CCNA2, ERBB2, CDC6, CCND3, CCNE1, RAF1, MCM5, RICTOR, CCNB2, MMP2, MLST8, MMP14, TIMP2 and MSH6 according to the their degree value, and the degree of the first 21 targets is higher than the average. In addition, among the 43 targets, the following were targets for both GC and CKI: CCND1, PIK3CA, AKT1, MAPK1, ERBB2 and MMP2. These targets may be key targets of CKI for the treatment of GC.

Cyclin D1 (CCND1) is an important cell cycle regulatory protein that primarily controls the transition of proteins from G1 to S in the cell cycle. The overexpression of CCND1 can cause a variety of tumors, increasing tumor cell proliferation, and affecting the prognosis of cancer patients [40]. Kaplan-Meier analysis indicated that an overexpression of CCND1 affected the survival of GC patients. Existing studies have shown that forcing the expression of miR-33a can inhibit the overexpression of CCND1 and inhibit the proliferation of GC cells [41, 42]. A study showed that CCND1 in MCF-7 breast cancer cells using CKI was significantly down-regulated compared with untreated cells [43].Therefore, CCND1 is predicted to be a significant target of CKI for the treatment of GC. Phosphatidylinositol-3-kinase catalytic subunit (PIK3CA) is the third most frequently mutated gene in GC [44, 45]. Clinical trial studies have shown that PIK3CA mutations are associated with increased tumor aggressiveness, especially in locoregional disease, and Akt activation in GC [46]. AKT1, also known as protein kinase B (PKB), is a serine/threonine protein kinase. Existing experiments have demonstrated that the activation of AKT inhibits the abundance of p53, leading to a decrease in the transcription of miR-365, an upregulation of CCND1 and cdc25A, and the promotion of gastric cell proliferation [47]. Therefore, it is speculated that both PIK3CA and AKT1 may be targets of GC treated with CKI. Furthermore, GC can be treated by inhibiting the expression of PIK3CA and AKT1, leading to an inhibition of the growth and proliferation of GC cells. Mitogen-activated protein kinases (MAPKs) are a group of silk/threonine kinases that are activated by a variety of signals. Mitogen-activated protein kinase 1 (MAPK1), also known as extracellular signal-regulated kinase 2 (ERK2), is involved in the regulation of cell differentiation and proliferation. MAPK1 is the core of the ras-MAPK signaling pathway. Activated MAPK1 transmits signals to the nucleus, phosphorylates and activates various nuclear transcription factors, and promotes and regulates some gene expression related to cell growth and differentiation. Abnormalities in MAPK1 can lead to the development of multiple cancers [48]. It has been demonstrated that MAPK1 can be mediated by multiple microRNAs to reduce the proliferation and metastasis of GC cells [49, 50]. Transcriptome analysis of CKI-treated MDA-MB-231 cells revealed that PIK3CA, AKT1, and MAPK1 all affect cell migration, which may be related to the role of CKI in the treatment of cancer [51]. The epidermal growth factor-related protein (ErbB) family of receptor tyrosine kinases plays an important role in epithelial cell development. The overexpression of human epidermal growth factor receptor 2 (ERBB2) promotes the differentiation and proliferation of cancer cells, which is recognized as one of the key causes of GC and seriously affects the prognosis of GC patients [52, 53]. At present, there is no in vivo or in vitro experiment to confirm that CKI achieves the therapeutic effect of GC by interfering with ERBB2, which is an important potential direction for future research. MMPs are a family of highly conserved endopeptidases that participate in a variety of physiological and pathological processes, such as inflammatory response, tissue ischemia and hypoxia injury, tissue fibrosis, neovascularization, and tumor adhesion, invasion, and metastasis [54]. MMP2 plays an important role in the invasion and metastasis of GC and can promote tumor invasion and metastasis by degrading the extracellular matrix, promoting neovascularization, and regulating cell adhesion [55]. Li et al. [9] found that various components of CKI can act on MMP2 and confirmed by Western blotting that CKI significantly downregulated MMP2.

The clustering analysis was carried out by using the MCODE algorithm in Cytoscape, and two modules were obtained by clustering. Among them, Module 1 had a higher number of nodes, edges, and scores, which may be the main functional module. To further understand the overall process of CKI in the treatment of GC, the GO and KEGG enrichment functions of the two modules obtained in this study were analyzed.

In the GO functional enrichment analysis, the target proteins were divided into different functional modules. Module 1 had the most biological process-related entries, including the regulation of cell mitosis, cell proliferation and transcription, protein phosphorylation, and protein signaling. Module 2 contained 13 related items, showing a strong correlation with phosphatidylinositol and its kinase-mediated signaling.

In the KEGG pathway enrichment analysis, the different modules obtained 76 related items by FDR < 0.01 screening analysis. The 76 related pathways included cancer pathways and disease pathways, most of which were related to signaling pathways, such as the PI3K-AKT signaling pathway, FoxO signaling pathway, mTOR signaling pathway, ErbB signaling pathway, AMPK signaling pathway, and HIF-1 signaling pathway. Therefore, the predicted mechanism of CKI for the treatment of GC is related to these pathways.

PI3K signaling is a critical regulator of many important cellular processes, including cell growth, metabolism, survival, metastasis, and resistance to chemotherapy. PI3K and AKT proteins are significantly overexpressed in tumor tissues. AKT is a downstream effector of PTEN/PI3K, and the inactivation of PTEN will result in an abnormal activation of AKT. p-AKT is an activated form of AKT that plays an important role in cell proliferation, differentiation and survival. Studies have found that it is abnormally activated in various tumors, such as lung cancer, breast cancer, and colorectal cancer. mTOR is located downstream of the PI3K/AKT/mTOR signaling pathway. After activation, it mediates important downstream signaling molecules and affects cell proliferation and inhibits apoptosis. High mTOR expression is found in many malignant tumors [38, 56, 57]. The epidermal growth factor receptor (EGFR) family plays an important role in growth and development, and its alteration is associated with a variety of conditions, including cancer. This receptor family mediates cell proliferation and survival through the MAPK and PI3K/KT signaling pathways. Distortion of the ErbB3/PI3 kinase pathway is common in GC, especially in diffuse tumors, leading to an excessive activation of the PI3K-AKT signaling pathway [58]. Based on a specific genotype, the FoxO family is a transcription factor that plays an important role in apoptosis, cell cycle arrest, DNA repair, and oxidative stress by regulating downstream target gene transcription levels. The most important pathway for the interaction with FoxO in different types of cancer is the PI3K-AKT pathway [59, 60]. In a current experiment, transcriptomic analysis of CKI-treated non-small cell lung cancer cell lines revealed that CKI does regulate ErbB, MAPK, PI3K/Akt and other pathways. Flow cytometry analysis confirmed that mTOR is located downstream of PI3K/Akt and is the main pathway for CKI to regulate cancer cells [61]. Therefore, these potential targets and pathways may be key for CKI treatment of GC.

The limitations of our study were as follows: Firstly, our study has some limitations because it focuses on genes that were obtained data in the network databases, without regard to gender, age, tumor classification, and staging. Secondly, our data research results need to be verified by corresponding experimental studies.


In summary, this study applied the network pharmacology method to study complex network relationships regarding the multicomponent, multitarget and multichannel aspects of CKI for the treatment of GC. We predicted that CKI controls the growth and metastasis of GC cells by regulating vital targets such as CCND1, PIK3CA, and AKT1, as well as important related pathways. The results of the study preliminarily verified and predicted the molecular mechanism of CKI against GC but still need further experimental verification. The results provide a basis for further exploration of the mechanism of action and of CKI for the treatment of GC and provides a reference for the study of the more complex mechanism of action of this Chinese herbal compound.

Availability of data and materials

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



Cyclin D1


Compound Kushen Injection


Database for Annotation, Visualization and Integrated Discovery


Extracellular signal-regulated kinase 2


False discovery rate


Gastric cancer


Genome Ontology


Kyoto Encyclopedia of Genes and Genomes


Mitogen-activated protein kinases


Online Mendelian Inheritance in Man


Protein Data Bank


Pharmacogenomics Knowledge Base


Protein kinase B


Negative logarithm of the experimental dissociation/suppression constant

PPI network:

Protein-protein interaction network


Search Tool for the Retrieval of Interacting Genes


Traditional Chinese medicine


Therapeutic Target Database


  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;0:1–31.

    Google Scholar 

  2. Marquéslespier JM, Gonzálezpons M, Cruzcorrea M. Current perspectives on gastric Cancer. Gastroenterol Clin N Am. 2016;45(3):413–28.

    Article  Google Scholar 

  3. Charalampakis N, Economopoulou P, Kotsantis I, Tolia M, Schizas D, Liakakos T, Elimova E, Ajani JA, Psyrri A. Medical management of gastric cancer: a 2017 update. Cancer Med. 2018;7(1):123–33.

    Article  PubMed  Google Scholar 

  4. Coccolini F, Montori G, Ceresoli M, Cima S, Valli MC, Nita GE, Heyer A, Catena F, Ansaloni L. Advanced gastric cancer: what we know and what we still have to learn. World J Gastroenterol. 2016;22(3):1139.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Xue JX, Zhu ZY, Bian WH, Yao C. The Traditional Chinese Medicine Kangai Injection as an Adjuvant Method in Combination with Chemotherapy for the Treatment of Breast Cancer in Chinese Patients: A Meta-Analysis. Evid Based Complement Alternat Med. 2018;1:1–16.

    Article  Google Scholar 

  6. Tao JJ, Visvanathan K, Wolff AC. Long term side effects of adjuvant chemotherapy in patients with early breast cancer. Breast. 2015;24(supplement 2):S149–53.

    Article  PubMed  Google Scholar 

  7. Qi F, Li A, Inagaki Y, Gao J, Li J, Kokudo N, Li XK, Tang W. Chinese herbal medicines as adjuvant treatment during chemo- or radio-therapy for cancer. Biosci Trends. 2010;4(6):297–307.

    PubMed  Google Scholar 

  8. Guo Y-m, Huang Y-x, Shen H-h, Sang X-x, Ma X, Zhao Y-l, Xiao X-h. Efficacy of compound Kushen injection in relieving cancer-related pain: a systematic review and meta-analysis. Evid Based Complement Alternat Med. 2015;2015:840742.

    PubMed  PubMed Central  Google Scholar 

  9. 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  CAS  Google Scholar 

  10. Sun FG, Sun XY. “Sophora flavescens”, in Shen Nong Ben Cao Jing. Scientific and Technical Documents Publishing House; 1999. p. 59.

    Google Scholar 

  11. Shen H, Qu Z, Harata-Lee Y, Aung TN, Cui J, Wang W, Kortschak RD, Adelson DL. Understanding the mechanistic contribution of herbal extracts in compound Kushen injection with Transcriptome analysis. Front Oncol. 2019;9:632.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Wang W, You RL, Qin WJ, Hai LN, Fang MJ, Huang GH, Kang RX, Li MH, Qiao YF, Li JW, Li AP. Anti-tumor activities of active ingredients in compound Kushen injection. Acta Pharmacol Sin. 2015;36(6):676.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Dai ZJ, Gao J, Ji ZZ, Wang XJ, Ren HT, Liu XX, Wu WY, Kang HF, Guan HT. Matrine induces apoptosis in gastric carcinoma cells via alteration of Fas/FasL and activation of caspase-3. J Ethnopharmacol. 2009;123:91–6.

    Article  CAS  PubMed  Google Scholar 

  14. Zhao Z, Fan H, Higgins T, Qi J, Haines D, Trivett A, Oppenheim JJ, Wei H, Li J, Lin H, Howard OM. Fufang Kushen injection inhibits sarcoma growth and tumor-induced hyperalgesia via TRPV1 signaling pathways. Cancer Lett. 2014;355(2):232–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Gao L, Wang XD, Niu YY, Duan DD, Yang X, Hao J, Zhu CH, Chen D, Wang KX, Qin XM, Wu XZ. Molecular targets of Chinese herbs: a clinical study of hepatoma based on network pharmacology. Sci Rep. 2016;6:24944.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Yang M, Chen JL, Xu LW, Ji G. Navigating Traditional Chinese Medicine Network Pharmacology and Computational Tools. Evid Based Complement Alternat Med. 2013;12:731969.

    Google Scholar 

  17. Fang J, Wang L, Wu T, Yang C, Gao L, Cai H, Liu J, Fang S, Chen Y, Tan W, Wang Q. Network pharmacology-based study on the mechanism of action for herbal medicines in Alzheimer treatment. J Ethnopharmacol. 2017;196:281–92.

    Article  CAS  PubMed  Google Scholar 

  18. Zhang Y, Mao X, Su J, Geng Y, Guo R, Tang S, Li J, Xiao X, Xu H, Yang H. A network pharmacology-based strategy deciphers the underlying molecular mechanisms of Qixuehe capsule in the treatment of menstrual disorders. Chin Med. 2017;12(1):23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Hsin K-Y, Matsuoka Y, Asai Y, Kamiyoshi K, Watanabe T, Kawaoka Y, Kitano H. Systemsdock: a web server for network pharmacology-based prediction and analysis. Nucleic Acids Res. 2016;44(Web Server issue:W507–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhang D, Wu J, Wang K, Duan X, Liu S, Zhang B. Which are the best Chinese herbal injections combined with XELOX regimen for gastric cancer? Medicine. 2018;97(12):e0127.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhang D, Zheng J, Ni M, Wu J, Wang K, Duan X, Zhang X, Zhang B. Comparative efficacy and safety of Chinese herbal injections combined with the FOLFOX regimen for treating gastric cancer in China a network meta-analysis. Oncotarget. 2017;8(40):68873–89.

    PubMed  PubMed Central  Google Scholar 

  22. Qin C, Zhang C, Zhu F, Xu F, Chen SY, Zhang P, Li YH, Yang SY, Wei YQ, Tao L, Chen YZ. Therapeutic target database update 2014: a resource for targeted therapeutics. Nucleic Acids Res. 2014;42(D1):1118–23.

    Article  CAS  Google Scholar 

  23. Amberger JS, Hamosh A. Searching online Mendelian inheritance in man (OMIM): a knowledgebase of human genes and genetic phenotypes. Curr Protoc Bioinformatics. 2017;58(1):1–2.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Thorn CF, Klein TE, Altman RB. PharmGKB: the pharmacogenomics Knowledge Base. Methods Mol Biol. 2013;1015:311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kim J, So S, Lee HJ, Park JC, Kim JJ, Lee H. DigSee: disease gene search engine with evidence sentences (version cancer). Nucleic Acids Res. 2013;41(W1):W510–7.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Piñero J, Bravo À, QueraltRosinach N. DisGeNET a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45(D1):D833–9.

    Article  PubMed  CAS  Google Scholar 

  27. Ma Y, Gao HM, Liu J, Chen LM, Zhang QW, Wang ZM. Identifcation and determination of the chemical constituents in a herbal preparation, compound Kushen injection, by Hplc and Lc-dad-Ms/Ms. J LIQ Chromatography R T. 2014;37:207–20.

    Article  CAS  Google Scholar 

  28. Kim S, Thiessen PA, Bolton EE, Chen J, Gang F, Gindulyte A, Han L, He J, He S, Shoemaker BA, Wang J, Yu B, Zhang J, Bryant SH. PubChem substance and compound databases. Nucleic Acids Res. 2016;44(D1):D1202–13.

    Article  CAS  PubMed  Google Scholar 

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

  30. Nickel J, Gohlke B-O, Erehman J, Banerjee P, Rong WW, Goede A, Dunkel M, Preissner R. SuperPred: update on drug classification and target prediction. Nucleic Acids Res. 2014;42(W1):W26–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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):32–8.

    Article  CAS  Google Scholar 

  32. Franz M, Lopes CT, Huck G, Dong Y, Sumer O, Bader GD. Cytoscape.Js: a graph theory library for visualisation and analysis. Bioinformatics. 2016;32(2):309–11.

    CAS  PubMed  Google Scholar 

  33. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, Jensen LJ, von Mering C. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(Database issue):D362–8.

    Article  CAS  PubMed  Google Scholar 

  34. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.

    Article  Google Scholar 

  35. Rose PW, Prlić A, Altunkaya A, Bi C, Bradley AR, Christie CH, Di Costanzo L, Duarte JM, Dutta S, Feng Z, Green RK, Goodsell DS, Hudson B, Kalro T, Lowe R, Peisach E, Randle C, Rose AS, Shao C, Tao Y-P, Valasatava Y, Voigt M, Westbrook JD, Woo J, Yang H, Young JY, Zardecki C, Berman HM, Burley SK. The RCSB protein data bank: integrative view of protein, gene and 3D structural information. Nucleic Acids Res. 2017;45(Database issue):D271–81.

    CAS  PubMed  Google Scholar 

  36. Hsin KY, Ghosh S, Kitano H. Combining machine learning systems and multiple docking simulation packages to improve docking prediction reliability for network pharmacology. PLoS ONE. 2013;8:e83922.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. DeLano WL (2002) The PyMOL molecular graphics system.

    Google Scholar 

  38. He P, Zhang Z, Liao W, Xu D, Ming F, Kang Y. Screening of gene signatures for rheumatoid arthritis and osteoarthritis based on bioinformatics analysis. Mol Med Rep. 2016;14(2):1587–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ajani JA, D’Amico TA, Almhanna K, Bentrem DJ, Chao J, Das P, Denlinger CS, Fanta P, Farjah F, Fuchs CS, Gerdes H, Gibson M, Glasgow RE, Hayman JA, Hochwald S, Hofstetter WL, Ilson DH, Jaroszewski D, Johung KL, Keswani RN, Kleinberg LR, Korn WM, Leong S, Linn C, Lockhart AC, Ly QP, Mulcahy MF, Orringer MB, Perry KA, Poultsides GA, Scott WJ, Strong VE, Washington MK, Weksler B, Willett CG, Wright CD, Zelman D, McMillian N, Sundar H. Gastric Cancer, Version 3.2016, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw. 2016;14(10):1286.

    Article  PubMed  Google Scholar 

  40. Hu Z, Zhou Z, Xiong G, Wang Y, Lai Y, Deng L, Yang J. Cyclin D1 G870A polymorphism and the risk of hepatocellular carcinoma in a Chinese population. Tumour Biol. 2014;35(6):5607–12.

    Article  CAS  PubMed  Google Scholar 

  41. Shan YS, Hsu HP, Lai MD, Hung YH, Wang CY, Yen MC, Chen YL. Cyclin D1 overexpression correlates with poor tumor differentiation and prognosis in gastric cancer. Oncol Lett. 2017;14(4):4517.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Wang Y, Zhou X, Shan B, Han J, Wang F, Fan X, Lv Y, Chang L, Liu W. Downregulation of microRNA-33a promotes cyclin-dependent kinase 6, cyclin D1 and PIM1 expression and gastric cancer cell proliferation. Mol Med Rep. 2015;12(5):6491.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Qu Z, Cui J, Harata-Lee Y, Aung TN, Feng Q, Raison JM, Kortschak RD, Adelson DL. Identifcation of candidate anti-cancer molecular mechanisms of compound kushen injection using functional genomics. Oncotarget. 2016;7(40):66003–19.

    PubMed  PubMed Central  Google Scholar 

  44. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of gastric adenocarcinoma. Nature. 2014;513:202–9.

    Article  CAS  Google Scholar 

  45. Cristescu R, Lee J, Nebozhyn M, Kim KM, Ting JC, Wong SS, Liu J, Yue YG, Wang J, Yu K, Ye XS, Do IG, Liu S, Gong L, Fu J, Jin JG, Choi MG, Sohn TS, Lee JH, Bae JM, Kim ST, Park SH, Sohn I, Jung SH, Tan P, Chen R, Hardwick J, Kang WK, Ayers M, Hongyue D, Reinhard C, Loboda A, Kim S, Aggarwal A. Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes. Nat Med. 2015;21:449–56.

    Article  CAS  PubMed  Google Scholar 

  46. Kim JW, Lee HS, Nam KH, Ahn S, Kim JW, Ahn SH, Park DJ, Kim HH, Lee KW. PIK3CA mutations are associated with increased tumor aggressiveness and Akt activation in gastric cancer. Oncotarget. 2017;8(53):90948–58.

    PubMed  PubMed Central  Google Scholar 

  47. Guo SL, Ye H, Teng Y, Wang YL, Yang G, Li XB, Zhang C, Yang X, Yang ZZ, Yang X. Akt-p53-miR-365-cyclin D1/cdc25A axis contributes to gastric tumorigenesis induced by PTEN deficiency. Nat Commun. 2013;4:1.

    Google Scholar 

  48. Albanell J, Rojo F, Baselga J. Pharmacodynamic studies with the epidermal growth factor receptor tyrosine kinase inhibitor ZD1839. Semin Oncol. 2001;28(5 Suppl 16):56–66.

    Article  CAS  PubMed  Google Scholar 

  49. Bojian F, Haorong W. MiR-378 inhibits progression of human gastric cancer MGC-803 cells by targeting MAPK1 in vitro. Oncol Res. 2012;20(12):557.

    Google Scholar 

  50. Lingyun D, Shengying W, Zhiguang S. Long noncoding RNA GAPLINC promotes gastric cancer cell proliferation by acting as a molecular sponge of miR-378 to modulate MAPK1 expression. Onco Targets Ther. 2018;11:2797–804.

    Article  Google Scholar 

  51. Nourmohammadi S, Aung TN, Cui J, Pei JV, De Ieso ML, Harata-Lee Y, Qu Z, Adelson DL, Yool AJ. Effect of compound Kushen injection, a natural compound mixture, and its identified chemical components on migration and invasion of Colon, brain, and breast Cancer cell lines. Front Oncol. 2019;9:314.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Mahlknecht G, Maron R, Mancini M, Schechter B, Sela M, Yarden Y. Aptamer to ErbB-2/HER2 enhances degradation of the target and inhibits tumorigenic growth. Proc Natl Acad Sci U S A. 2013;110(20):8170–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Duarte HO, Balmaña M, Mereiter S, Osório H, Gomes J, Reis CA. Gastric Cancer cell glycosylation as a modulator of the ErbB2 oncogenic receptor. Int J Mol Sci. 2017;18(11):2262.

    Article  PubMed Central  CAS  Google Scholar 

  54. Jabå Oå Ska-Trypuä A, Matejczyk M, Rosochacki S. Matrix metalloproteinases (MMPs), the main extracellular matrix (ECM) enzymes in collagen degradation, as a target for anticancer drugs. J Enzyme Inhib Med Chem. 2016;23(sup1):177–83.

    Google Scholar 

  55. Burlaka AP, Ganusevich II, Gafurov MR, Lukin SM, Sidorik EP. Stomach Cancer: interconnection between the redox state, activity of MMP-2, MMP-9 and stage of tumor growth. Cancer Microenviron. 2016;9(1):27–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Fang WL, Huang KH, Lan YT, Lin CH, Chang SC, Chen MH, Chao Y, Lin WC, Lo SS, Li AF, Wu CW, Chiou SH, Shyr YM. Mutations in PI3K/AKT pathway genes and amplifications of, PIK3CA, are associated with patterns of recurrence in gastric cancers. Oncotarget. 2016;7(5):6201–20.

    Article  PubMed  Google Scholar 

  57. Tapia O, Riquelme I, Leal P, Sandoval A, Aedo S, Weber H, Letelier P, Bellolio E, Villaseca M, Garcia P, Roa JC. The PI3K/AKT/mTOR pathway is activated in gastric cancer with potential prognostic and predictive significance. Virchows Arch. 2014;465(1):25–33.

    Article  CAS  PubMed  Google Scholar 

  58. Ying J, Xu Q, Liu B, Zhang G, Chen L, Pan H. The expression of the PI3K/AKT/mTOR pathway in gastric cancer and its role in gastric cancer prognosis. Onco Targets. 2015;8(default:2427–33.

    Article  CAS  Google Scholar 

  59. Zhang XL, Yang YS, Xu DP, Qu JH, Guo MZ, Gong Y, Huang J. Comparative study on overexpression of HER2/neu and HER3 in gastric cancer. World J Surg. 2009;33:2112–8.

    Article  PubMed  Google Scholar 

  60. Brunet A, Greer EL. FoxO transcription factors at the interface between longevity and tumor suppression. Oncogene. 2005;24(50):7410–25.

    Article  PubMed  CAS  Google Scholar 

  61. Zhang J, Qu ZP, Yao H, Sun LL, Harata-Lee Y, Cui J, Aung NT, Liu XM, You GL, Wang W, Hai L, Adelson DL. An effective drug sensitizing agent increases gefitinib treatment by down regulating PI3K/Akt/mTOR pathway and up regulating autophagy in nonsmall cell lung cancer. Biomed Pharmacother. 2019;118:109169.

    Article  PubMed  CAS  Google Scholar 

Download references


No applicable.


The study was supported by the National Natural Science Foundation of China to Jiarui Wu (Grant No. 81473547; Grant No. 81673829) and Young Scientists Training Program of Beijing University of Chinese Medicine. The funding body approved the design of the study, analysis, and interpretation of data, and publication of the manuscript.

Author information

Authors and Affiliations



ZW and WJR conceived and designed the study. ZYL, JSS and MZQ collected the ingredients of CKI, targets of CKI and GC. LXK, LSY, GSY and NMW performed the network pharmacology analysis and molecular docking verification, and ZW and ZJY wrote the manuscript. All authors were responsible for reviewing data. All authors read and approved the final manuscript.

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 that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, W., Wu, J., Zhu, Y. et al. Study on the mechanisms of compound Kushen injection for the treatment of gastric cancer based on network pharmacology. BMC Complement Med Ther 20, 6 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: