Abstract Background Zuojinwan (ZJW), a famous Chinese medicine formula, has been widely used to treat colorectal cancer (CRC). However, its bioactive compounds, potential targets, and molecular mechanism remain largely elusive. Aim A network pharmacology-based strategy combined with molecular docking studies and in vitro validation were employed to investigate bioactive compounds, potential targets, and molecular mechanism of ZJW against CRC. Materials and Methods Bioactive compounds and potential targets of ZJW, as well as related genes of CRC, were acquired from public databases. Important ingredients, potential targets, and signaling pathways were determined through bioinformatics analysis, including protein–protein interaction (PPI), the Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Subsequently, molecular docking and cell experiments were performed to further verify the findings. Results A total of 36 bioactive ingredients of ZJW and 163 gene targets of ZJW were identified. The network analysis revealed that quercetin, baicalein, wogonin, beta-sitosterol, and isorhamnetin may be candidate agents. The AKT1, JUN, CDKN1A, BCL2L1, and NCOA1 could become potential drug targets. The KEGG indicated that PI3K-AKT signaling pathway may play an important role in the effect of ZJW against CRC. Molecular docking suggested that quercetin, baicalein, and wogonin combined well with AKT1 and JUN. The in vitro experiment showed that quercetin, the most important ingredient of ZJW, could induce apoptosis of HCT116 cells through PI3K-Akt signaling pathway. This finding was congruent with the prediction obtained through the network pharmacology approach. Conclusion This study comprehensively illuminated the active ingredients, potential targets, and molecular mechanism of ZJW against CRC. It also provided a promising approach to uncover the scientific basis and therapeutic mechanism of traditional Chinese medicine (TCM) formula treating for disease. Keywords: Zuojinwan, colorectal cancer, network pharmacology, PI3K/AKT signal pathway Introduction Colorectal cancer (CRC) is considered the third most common cancer among men and second among women, and the second most common cause of cancer death worldwide.[42]1 The incidence and mortality of CRC is rapidly rising particularly in low- and middle-income countries, and it is estimated that the global burden of CRC increases by 60% over 2.2 million new case and 1.1 million cancer death by 2030.[43]2,[44]3 Despite recent advances in therapy and multidisciplinary care, patients with CRC continue to suffer from severe adverse actions, which impair prognosis and reduce the quality of life. More than 50% of CRC patients were diagnosed at an advanced stage, and the 5-year survival rate is <14%.[45]4–[46]6 Hence, the development of alternative therapies with low toxicity is important for providing more effective clinical therapy and alleviating adverse reactions. Traditional Chinese medicine (TCM), as a crucial complementary and alternative medicine, offers beneficial effects to patients with cancer and has been widely accepted worldwide.[47]7 Previous studies demonstrated that traditional herbal medicine, as adjuvant therapy, combined with chemotherapy or radiotherapy, could enhance treatment effectiveness, alleviate adverse actions, improve quality of life, and prolong survival time.[48]8–[49]12 Owing to its effectiveness and fewer side effects, it is increasingly popular in Western countries. Zuojinwan (ZJW), a classical TCM formula, was first described in Danxi Xinfa for a long history.[50]13 Besides, ZJW was also officially listed in the Chinese Pharmacopoeia.[51]14 ZJW, consists of Coptis chinensis Franch (R. coptidis) and Evodia ruticarpa (A.Juss.) Hook.f. and Thomson (E. rutaecarpa) at a ratio of 6:1 (w/w), was initially employed in treating gastro-intestinal disorders in China.[52]15–[53]17 From bench to bedside, an increasing body of evidence suggests that ZJW exerts favorable anti-cancer effects. It has been demonstrated that ZJW significantly suppresses the growth of eight kinds of human cancer cell lines by inducing mitochondrial apoptosis.[54]14 Besides, oral administration of ZJW notably inhibited tumor growth in HepG2 xenograft-bearing immunocompetent mice and breast cancer ZHENG model.[55]18,[56]19 Furthermore, ZJW extracts inhibit CRC cells growth and invasion, and enhance the sensitivity of chemotherapy of CRC cells and gastric cancer cells to chemotherapy.[57]20–[58]22 Similar to other TCM formulae, ZJW involves multiple components, targets, and pathways against tumors. However, its major bioactive ingredients and pharmacological molecular mechanism against CRC remain poorly understood. Network pharmacology is a promising approach to uncover the scientific basis and therapeutic mechanism of TCM formulae, and identify new drugs.[59]23 With the rapid development of bioinformatics and pharmacology, network pharmacology comprehensively investigates interrelationships of drugs, targets and diseases, and shows the network of drug-targets-disease in a visual way. This approach clearly observes the effects of drugs on disease,[60]24 is consistent with the TCM theory that emphasizes synergy between Chinese medicine.[61]25 Besides, new important bioactive compound of TCM formula could be found based on network pharmacology. In the present study, we aim to explore and verify the molecular mechanisms and pathways of the main bioactive ingredients of ZJW in the treatment of CRC based on a network pharmacology with molecular docking and in vitro validation. The workflow is shown in [62]Figure 1. Figure 1. [63]Figure 1 [64]Open in a new tab Workflow for network pharmacology-based prediction and validation of anti-CRC mechanism. Materials and Methods Data Preparation Composition of ZJW All ingredients of ZJW were acquired from three public databases including Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform (TCMSP, available online: [65]http://tcmspw.com/tcmsp.php), TCM Database@Taiwan ([66]http://tcm.cmu.edu.tw/), and A Bioinformation Analysis Tool for Molecular mechANism of Traditional Chinese Medicine (BATMAN-TCM, available online: [67]http://bionet.ncpsb.org/batman-tcm/index.php/home/index/result/jobI d/batman-I2019-04-16-17745-1555420409.html) and related studies. Pharmacokinetic ADME Prediction The TCMSP database was used to select all compounds based on absorption, distribution, metabolism, and excretion (ADME) criteria. In this study, oral bioavailability (OB) and drug-likeness (DL) were adopted to identify the bioactive ingredients of ZJW. OB described the percentage of an oral dose capable of producing pharmacological activity.[68]26 DL was applied to determine the similarity of a compound to conventional drugs and its physicochemical properties.[69]27 Only compounds that had related targets and met the ADME criteria (ie, OB threshold≥30% and DL threshold≥0.18) were included in the study for subsequent research. Prediction Targets of Compounds The TCMSP and DRUGBANK ([70]https://www.drugbank.ca/) databases, a unique bioinformatics resource which combines detailed drug data with comprehensive drug target in formation, were utilized to predict the targets of all ingredients.[71]28 The Uniprot database ([72]https://www.uniprot.org/), which provides a comprehensive, high-quality, and free source of protein sequence and functional information, was applied to convert protein target names to their corresponding gene symbols. The predicted targets of CRC were collected from the GeneCards ([73]https://www.genecards.org/, version 4.9.0) and OMIM ([74]http://www.omim.org/, updated April 23, 2019) databases. Bioinformatics Analysis Protein–Protein Interaction (PPI) Data The overlapping gene of compounds and CRC were considered as hub genes and analyzed using online STRING ([75]https://string db.org/cgi/network.pl?taskId=wJDcOUv6JCzz, version 11.0, updated January 19, 2019) to obtain the PPI, with the species limited to “Homo sapiens” and a confidence score >0.990. The TSV format file, which was downloaded from the STRING database, was imported into the Cytoscape3.7.1 for analysis. Cytoscape 3.7.1, which is widely applied to network pharmacology research, was used to construct and visualize the network. Furthermore, it provides a basic set of characters for data integration, analysis, and visualization for complicated network analysis. Two key topological parameters (ie, degree and betweenness centrality) characterize the most important nodes in the network; higher quantitative values of topological parameters indicate greater importance of the node. Enrichment Analysis Based on the R-package clusterProfiler, the Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to the target proteins. A p<0.001 denoted statistical significance. The KEGG database, which was established in 1995 by the Kanehisa laboratory in the Bioinformatics Center of Kyoto University, is typically used in network pharmacology to annotate the pathways of target genes. Network Construction Construction of the network was performed as follows: (1) herbs-bioactive compound-target network of ZJW (C-T network); (2) herbs-bioactive compound-disease-target network of treatment with ZJW for CRC (C-D-T); (3) PPI network of treatment with ZJW for CRC; and (4) bioactive compound-pathway-target (C-P-T) network. Validation Molecular Docking Molecular docking, which combines known protein with small compounds, is a crucial technology of network pharmacology. Molecular docking was performed using the Molecular Operating Environment (MOE) (v2015.10) software to validate compound-target interaction.[76]29 The 3D structure of protein was obtained from the Protein Data Bank (PDB) ([77]http://www.rcsb.org/) and then imputed it into MOE to construct mating pocket of molecular docking after deleting water molecules, preparing protein structure and minimizing energy. The mating pocket was performed molecular docking with compounds whose structure was downloaded from PubChem ([78]https://pubchem.ncbi.nlm.nih.gov/). Cell Experiments Cell Culture HCT116 and HT29 cells, a human colon cancer cell line, were purchased from GeneChem Co., Ltd. (Shanghai, China). The cells were cultured in RPMI-1640 medium containing 10% fetal bovine serum, 100 µg/mL ampicillin, and 100 µg/mL streptomycin, at 37°C with 5% CO2 and 95% air. Cell Viability Assay In brief, HCT116 and HT29 cells were seeded in 96-well plates at 4×10^3 cells/well and treated with different concentrations of quercetin (0, 100, 150, 200, and 250 μΜ) for 24, 48, and 72 h. The Cell Counting Kit-8 (CCK-8) assay was utilized to detected cell proliferation at 450 nm wavelength. The Graphpad prism 7.0 software was used to analyze and plot the data. Flow Cytometry Analysis for Apoptosis HCT116 and HT29 cells were incubated in 6-well plates at 37°C and treated with quercetin at different concentrations (0, 100, 150, and 200μΜ) for 72 h. Subsequently, cell apoptosis was measured through flow cytometry with the Annexin V-fluorescein isothiocyanate/propidium iodide (Annexin V-FITC/PI) apoptosis detection kit (40302ES20; Yeasen Biotechnology Co. Ltd., Shanghai, China) according to the manufacturer instruction. Briefly, treated cells were collected and washed twice with cold phosphate-buffered saline, resuspended in 400 µL 1×binding buffer with 5 µL Annexin V-FITC and 5 µL PI staining solution, and incubated for 15 min in the dark at room temperature. Next, 400 µL 1×binding buffer was added and mixed. Finally, cell apoptosis was determined using flow cytometry (Guava PCA; Millipore, USA), and data were analyzed using the software FlowJo 7.6 (De Novo Software, Los Angeles, CA, USA). Western Blotting Assay Total protein was obtained from treated cells in radioimmunoprecipitation assay buffer (P0013B; Beyotime, Shanghai, China) and centrifuged at 12,000 × g and 4°C for 15 min. Proteins were separated via sodium dodecyl-sulfate polyacrylamide gel electrophoresis and transferred onto polyvinylidene fluoride membranes (Millipore, Bedford, MA, USA). After blocking in 5% bovine serum albumin (Gen-view Scientific Inc., USA) for 1.5 h at room temperature, the membranes were incubated with primary antibodies against glyceraldehyde-3-phosphate dehydrogenase (Abclonal), PI3K (cell signal technique), AKT (cell signal technique), p-AKT (cell signal technique), BCL2 (Abclonal), and BAX (cell signal technique) at 4°C overnight. Next, the membranes were incubated with horseradish-peroxidase-conjugated secondary antibodies (proteintech) at room temperature for 1 h. Subsequently, the membranes were incubated with horseradish-peroxidase-conjugated secondary antibodies (proteintech) at room temperature for 1 h. The Enhanced Chemiluminescence Detection Kit ([79]P10100; New Cell & Molecular Biotech Co., Ltd.) was used to detect and visualize protein bands. Results Data Statistics Ingredients of ZJW ZJW consists of R. coptidis and E. rutaecarpa. A total of 245 ingredients were recognized in ZJW, including 62 in R. coptidis and 183 in E. rutaecarpa. Among the 245 ingredients, eight were duplicated and thus eliminated, while the remaining 237 ingredients were eventually included in the analysis. According to a reference, three ingredients were included: alpha-bulnesene, aspartate, and nerolidol. The details are described in [80]Table S1 (supplementary material). Acquiring Compounds by ADME Screening and C-T Network Analysis Two ADME-related models (ie, OB and DL) were utilized to screen all compounds with related targets. In total, after removing two duplicated compounds (ie, berberine and quercetin), 36 compounds (OB≥30% and DL≥0.18) were included in the study (13 in R. coptidis and 25 in E. rutaecarpa). The results are listed in [81]Table 1. The C-T network includes 226 nodes and 587 edges shown in [82]Figure 2. The pink ellipse represents 188 potential target genes of compounds. The blue triangle, red v and yellow diamond stand for compounds of E. rutaecarpa, R. coptidis and common compounds of E. rutaecarpa and R. coptidis, respectively. Furthermore, higher degrees indicated larger sizes, which present important nodes in the network. Quercetin, baicalein, wogonin, beta-sitosterol, and isorhamnetin, which exhibited higher degrees, play major roles in the effect of ZJW against CRC. Table 1. Bioactive Compounds of ZJW Mol ID Molecule Name OB (%) DL MOL001454 berberine 36.86 0.78 MOL000098 quercetin 46.43 0.28 MOL002894 berberrubine 35.74 0.73 MOL002897 epiberberine 43.09 0.78 MOL002903 (R)-Canadine 55.37 0.77 MOL002904 Berlambine 36.68 0.82 MOL002907 Corchoroside A_qt 104.95 0.78 MOL000622 Magnograndiolide 63.71 0.19 MOL000785 palmatine 64.6 0.65 MOL001458 coptisine 30.67 0.86 MOL002668 Worenine 45.83 0.87 MOL002714 baicalein 33.52 0.21 MOL000173 wogonin 30.68 0.23 MOL002662 rutaecarpine 40.3 0.6 MOL000354 isorhamnetin 49.6 0.31 MOL000358 beta-sitosterol 36.91 0.75 MOL000359 sitosterol 36.91 0.75 MOL003943 Rutalinidine 40.89 0.22 MOL003947 1-methyl-2-[(Z)-pentadec-10-enyl]-4-quinolone 48.45 0.46 MOL003956 dihydrorutaecarpine 42.27 0.6 MOL003957 1-methyl-2-pentadecyl-4-quinolone 44.52 0.46 MOL003958 Evodiamine 86.02 0.64 MOL003960 1-(5,7,8-trimethoxy-2,2-dimethylchromen-6-yl)ethanone 30.39 0.18 MOL003963 hydroxyevodiamine 72.11 0.71 MOL003964 1-methyl-2-undecyl-4-quinolone 47.59 0.27 MOL003972 1-methyl-2-nonyl-4-quinolone 48.42 0.2 MOL003974 Evocarpine 48.66 0.36 MOL003988 2-Hydroxy-3-formyl-7-methoxycarbazole 83.08 0.18 MOL003994 24-methyl-31-norlanost-9(11)-enol 38 0.75 MOL004004 6-OH-Luteolin 46.93 0.28 MOL004014 Evodiamide 73.77 0.28 MOL004018 Goshuyuamide I 83.19 0.39 MOL004019 GoshuyuamideII 69.11 0.43 MOL004020 gossypetin 35 0.31 MOL004021 Gravacridoneshlirine 63.73 0.54 MOL004025 N-(2-Methylaminobenzoyl) tryptamine 56.96 0.26 [83]Open in a new tab Abbreviations: OB, oral bioavailability; DL, drug-likeness. Figure 2. [84]Figure 2 [85]Open in a new tab C-T network of ZJW. Pink ellipse represents genes. The blue triangle, red v and yellow diamond stand for compounds ZJW. The green diamond stand for E. rutaecarpa and R. coptidis. Potential Target Genes and C-D-T Network Analysis After removing duplication values and transferring gene symbol, 188 potential target genes of compounds were collected from the TCMSP and DRUGBANK databases. CRC was used as a keyword for the search, and a total of 7481 related target genes for CRC were collected from the GeneCards and OMIM databases. Among 188 of compounds and 7481 CRC-related target genes, there were 163 overlapping related target genes, considered as the hub genes for subsequent research. The results are shown in [86]Table S2 (see supplementary material). The C-D-T network contains 202 nodes and 786 edges ([87]Figure 3). The red circle stands for gene, while the blue, pink, and navy-blue triangles stand for compounds of ZJW. The green diamond stands for CRC, E. rutaecarpa, and R. coptidis. The size of the shape is proportional to the degree. Higher degrees indicated that the compounds or genes play more important anti-CRC roles. Figure 3. [88]Figure 3 [89]Open in a new tab C-D-T network of ZJW. The red circle stands for gene, blue, pink and navy-blue triangle stand for compounds of ZJW. The green diamond stand for CRC, E. rutaecarpa and R. coptidis. Bioinformatics Analysis Analysis on CRC-Related PPI The 163 CRC-related hub genes were submitted to STRING to obtain PPI data. The network was subsequently constructed and visualized using Cytoscape 3.7.1. As shown in [90]Figure 4, the PPI network included 69 nodes and 76 edges. Higher degrees indicated larger node sizes; red indicates the highest degree, and yellow represents the lowest degree. The nodes with the top five degrees were AKT1, JUN, CDKN1A, BCL2L1, and NCOA1. These nodes may play a pivotal role in the PPI network for the anti-CRC effect. Figure 4. [91]Figure 4 [92]Open in a new tab PPI network of anti-CRC-related protein. Higher degrees indicated larger node sizes, and red indicates higher degree, and yellow represents the lowest degree. Enrichment Analyses GO Analysis The R-package clusterProfiler was used to elucidate the biological process (BP), cell composition (CC), and molecular function (MF) annotation of the 163 hub target proteins. There were 902 GO terms, including 828 of BP, 24 of CC, and 50 of MF. Most GO terms of BP were related to response to lipopolysaccharide, response to molecule of bacterial origin, response to oxygen levels, response to oxidative stress, response to reactive oxygen species, response to steroid hormone, reactive oxygen species metabolic process, response to hypoxia, response to antibiotic, response to decreased oxygen levels, etc. The main terms of CC were associated with membrane raft, membrane microdomain, membrane region, caveola, transcription factor complex, etc. MF enrichment was mainly involved in nuclear receptor activity, transcription factor activity, direct ligand-regulated sequence-specific DNA binding, steroid hormone receptor activity, cytokine receptor binding, DNA-binding transcription activator activity, RNA polymerase II-specific, etc. The top 20 terms of BP, CC, and MF, which were ranked according to their adjusted p-value, are shown in [93]Figure 5. Low-adjusted p-values and red color indicated greater enrichment of the GO terms. Figure 5. [94]Figure 5 [95]Open in a new tab BP, CC, MF of GO analysis. KEGG Analysis and C-P-T Network Analysis The KEGG pathway annotation indicated that 149 of the 163 target genes were enriched and yielded 97 pathways (p< 0.001). It verified CRC-related pathways including the PI3K-AKT signaling pathway, MAPK signaling pathway, TNF signaling pathway, HIF-1 signaling pathway, etc. Furthermore, the top 20 KEGG pathways with high counts are shown in [96]Figure 6. The results distinctly suggested that the PI3K-AKT signaling pathway, which enriched a larger count among CRC-related pathways, played a key role in the treatment of CRC. The C-P-T network, containing 150 nodes and 762 edges, was constructed ([97]Figure 7). The red circle nodes denote genes, pink triangle nodes denote compounds of ZJW, blue diamond nodes denote pathways, and green octagon nodes denote E. rutaecarpa and R. coptidis. Figure 6. [98]Figure 6 [99]Open in a new tab Top 20 KEGG pathway with high count. Figure 7. [100]Figure 7 [101]Open in a new tab C-P-T network of top 20 pathway. The red circle nodes stand for gene, pink triangle nodes stand for compounds of ZJW, blue diamond stand for pathway, green octagon nodes stand for E. rutaecarpa and R. coptidis. Validation Molecular Docking Molecular docking was applied to validate the binding action mode of AKT1 and JUN with quercetin, wogonin, and baicalein, respectively, according to the network analysis. The results revealed that AKT1 interacted with quercetin, baicalein, and wogonin, as shown in [102]Figure 8A–[103]C, respectively. As shown in [104]Figure 8A, the structure of quercetin could form two hydrogen bonds and one π-π bond with Ile290, Gln79, and Trp80 in AKT1. In [105]Figure 8B, the structure of baicalein could form one hydrogen bond, one H-π bond, and two π-H bonds with Gln79, Trp80, Leu210, and Asp292 in AKT1. As shown in [106]Figure 8C, one hydrogen bond was formed between wogonin and AKT1. Figure 8. [107]Figure 8 [108]Open in a new tab (A–C) represent AKT1 interacted action mode with quercetin, baicalein and wogonin respectively. The results showed that JUN interacted with quercetin, wogonin, and baicalein, as shown in [109]Figure 9A–[110]C, respectively. As shown in [111]Figure 9A, the structure of quercetin could form a hydrogen bond and a π-H bond with Ser144 and Val150 in JUN, respectively. Similarly, in [112]Figure 9B, a hydrogen bond and a π-H bond were formed to bind baicalein with Glu243 and Gly143 in JUN, respectively. Wogonin interacted with Gly143 and Val150 in JUN through two π-H bonds ([113]Figure 9C). Figure 9. [114]Figure 9 [115]Open in a new tab (A–C) represent JUN interacted action mode with quercetin, baicalein and wogonin respectively. In conclusion, there is an interplay between ingredients and protein targets with different bonds. The binding energy was counted to assess the matching degree of ingredients with proteins. The detail binding energy score is shown in [116]Table 2. Lower binding energy indicated greater stability. The results suggested that the ingredients could bind to the active site of protein targets. Table 2. The Binding Energy Score of Ingredients in ZJW Binding with AKT1 and JUN. AKT1, RAC-Alpha Serine/Threonine-Protein Kinase; JUN, Transcription Factor AP-1 Ligands Receptors Quercetin Baicalein Wogonin AKT1 −6.9308 −6.5591 −6.9289 JUN −6.9 −7.0652 −7.3975 [117]Open in a new tab Quercetin Suppressed the Proliferation of CRC Cells CCK-8 assays were performed to determine the anti-cancer effect of treatment with quercetin at different concentrations (0, 100, 150, 200, and 250 µM) for 24, 48, and 72 h in HCT116 and HT29 cells. As shown in [118]Figure 10A and [119]B, the viability of HCT116 and HT29 cells was markedly reduced by quercetin in a dose- and time-dependent manner. The IC50 values for quercetin were 225.7, 186.3, and138.6 µM in HCT116 cells; 221.6, 185.4, 157.4 µM in HT29 cells for 24, 48, and 72 h, respectively. Figure 10. [120]Figure 10 [121]Open in a new tab Quercetin inhibits the viability of CRC cells as measured by the CCK8 assay. (A) HCT116, (B) HT29 cells were treated with various concentrations of quercetin for 24, 48, and 72 h, respectively. Data are presented as the mean ± SD from at least three independent experiments. Quercetin Promoted Apoptosis of CRC Cells The effect of quercetin on the apoptosis of HCT116 and HT29 cells was detected through flow cytometric analysis. After treatment with 0, 100, 150, and 200 µM of quercetin for 72 h, CRC cells were stained with Annexin V-FITC and PI to determine the degree of apoptosis. [122]Figure 11 shows that the percentage of apoptotic cells was notably increased in a dose-dependent manner after treatment with quercetin, strongly indicating that quercetin induced apoptosis in HCT116 and HT29 cells. Figure 11. [123]Figure 11 [124]Open in a new tab HCT116 (A) and HT29 (B) cells were treated with various concentrations quercetin. Apoptosis was analyzed by flow cytometry after Annexin V-FITC/PI staining. (C) The percentage of apoptotic cells was presented as the mean ± SD of three independent experiments, **P < 0.01 vs control; one-way ANOVA, followed by Dunnett’s post-hoc test. Quercetin Inhibits the PI3K/AKT Pathway in CRC Cells We further explored the molecular mechanism of quercetin involved in inducing apoptosis in HCT116 and HT29 cells. The PI3K/AKT pathway, which is a key pathway associated with the effect of ZJW against CRC, was analyzed through Western blotting. As shown in [125]Figure 12, the protein levels of PI3K, AKT, p-AKT, and BCL2 were significantly decreased in a dose-dependent manner. In contrast, the protein levels of BAX were significantly increased. The results demonstrated that quercetin induced apoptosis in HCT116 and HT29 cells via the PI3K/AKT pathway. Figure 12. [126]Figure 12 [127]Open in a new tab Quercetin inhibits the PI3K/AKT pathway in CRC cells. Representative Western blots showing the status of PI3K, AKT, p-AKT, BCL2, BAX in HCT116 (A1) and HT29 (B1) cells. GAPDH was used as an internal control. The protein levels of PI3K, AKT, p-AKT, BCL2, BAX in HCT116 (A2) and HT29 (B2) cells. Values are presented as a mean ± SEM. * P < 0.05, ** P < 0.01 vs control. Discussion In the present study, a network pharmacology-based approach, molecular docking and in vitro experimental validation were performed to reveal bioactive compounds and the molecular mechanism of ZJW for treating CRC. A total of 36 bioactive compounds and 163 protein targets were selected from public databases. The PPI network showed the interaction of each protein and indicated that AKT1, JUN, CDKN1A, BCL2L1, and NCOA1 may play key roles in the effect of ZJW against CRC. The GO enrichment analysis annotated the function of protein targets from BP, CC, and MF. The KEGG pathway enrichment analysis revealed that 149 of the 163 target proteins significantly (p<0.001) enriched 97 cancer-related signal pathways, including the PI3K-AKT signaling pathway, MAPK signaling pathway, TNF signaling pathway, HIF-1 signaling pathway, etc. Molecular docking suggested that important ingredients combined well with target proteins. The cell experiments verified that the key ingredient quercetin inhibited proliferation and induced apoptosis by suppressing the PI3K-AKT signaling pathway. In recent years, naturally bioactive compounds are attracting considerable attention as anti-cancer agents, owing to their high therapeutic value and low systemic toxicity.[128]30 Previous studies demonstrated that berberine and evodiamine, the major active compounds of ZJW, exhibited anti-cancer effects (colon cancer, breast cancer, liver cancer, lung cancer, CRC, etc.) by inhibiting proliferation, and inducing apoptosis and autophagic cell death.[129]31–[130]35 However, to the best of our knowledge, this is the first study to report that quercetin, which shows the highest degree among all compounds, is the most important compound of ZJW in the treatment of CRC. Quercetin belongs to the group of polyphenolic compounds termed flavonoids and is abundantly distributed in plants, foods, and beverages.[131]36,[132]37 Previous evidence showed that quercetin exerted its anti-cancer effects by suppressing proliferation and metastasis, and inducing apoptosis and autophagy.[133]38–[134]41 Chuang et al demonstrated that quercetin enhanced the anti-cancer activity of trichostatin A by upregulating the p300 protein in lung cancer.[135]42 Lu et al showed that quercetin inhibited the migration and invasion of HCCLM3 cells by downregulating the expression of p-AKT1, Matrix Metalloproteinase (MMP) MMP-2, and MMP-9.[136]43 Refolo et al revealed that quercetin suppressed proliferation and apoptosis was mediated by CB1 receptor.[137]44 Moreover, quercetin enhances hypoxia-mediated apoptosis in HCT116 cells through AMPK-related signaling pathway.[138]45,[139]46 Furthermore, quercetin increased the sensitivity of MSI CRC cells to 5-FU.[140]47 The PPI results indicated that AKT1, JUN, BCL2L, CDKN1A, and NCOA1, especially AKT1, may play a key role in the anti-CRC process. The KEGG results indicated that the PI3K-AKT signaling pathway, which exhibited the highest gene count enrichment among the cancer-related signal pathways, plays a critical role in the effect of ZJW against CRC. AKT1, belonging to the AKT family, regulates biological functions, such as cell survival, proliferation, metabolism, and growth to affect cancer growth.[141]48,[142]49 AKT is involved in the treatment of numerous types of cancer via various signaling pathways, especially the PI3K-AKT signaling pathway. Accumulating evidence indicates that the PI3K/AKT signaling pathway, which is one of the most frequently activated signal-transduction pathways in cancer, is involved in cell cycle, cell proliferation, growth, migration, angiogenesis, and apoptosis.[143]50–[144]53 Liu et al demonstrated that aloperine induced apoptosis and arrested the cell cycle in hepatocellular carcinoma cells through the PI3K-AKT signaling pathway.[145]54 Hossan et al verified that cerberin inhibited cell proliferation, migration, and survival in pancreatic, triple-negative breast, and non-small cell lung cancer through the PI3K/AKT signaling pathway.[146]55 Hence, molecular docking and in vitro experiments were performed to further verify the molecular mechanism of quercetin which possesses the highest degree and may be the most important bioactive compound of ZJW for treating for CRC. The results of molecular docking showed that quercetin interacted with AKT1 through three bonds, indicating that AKT1 may be a key action target for the effect of quercetin against CRC. Moreover, the cell experiments demonstrated that quercetin for the treatment of CRC suppressed the proliferation and induced apoptosis in HCT116 cell by inhibiting the PI3K-AKT signaling pathway with the downregulation of PI3K, AKT, and BCL2. The present study has several limitations. Firstly, the public databases investigated in the study are constantly updated; thus, some other bioactive ingredients and target genes may not have been included in our analysis. In addition, other signaling pathways (eg, MAPK, TNF, HIF-1 signaling pathways) may also be involved in the anti-tumor effect of ZJW. Further studies are warranted to examine the potential involvement of these pathways. Moreover, the quercetin could not really represent for ZJW, even though it was found the most important bioactive ingredient of ZJW in the treatment of CRC. Therefore, additional research is required to further explore the molecular mechanism of ZJW in the treatment of CRC in vitro and in vivo. Conclusion In the present study, we combined network pharmacology-based prediction, molecular docking, and in vitro experiments to verify the targets of ZJW and the mechanisms through which its main bioactive ingredients exert their effectively anti-CRC. We demonstrated that ZJW treated CRC by activating its targets, which play key roles in inhibiting cell proliferation and promoting apoptosis. Based on this multidisciplinary strategy, the present study provided a promising approach for the treatment of disease using TCM. Acknowledgments