Characterizing the Antitumor Effect of Coptis chinensis and Mume Fructus against Colorectal Cancer Based on Pharmacological Analysis
Colorectal cancer (CRC) is the third most diagnosed cancer worldwide and is a significant cause of cancer-related deaths. Previous studies have observed that Coptis chinensis (CC) and Mume Fructus (MF) are effective against CRC, enteritis, and intestinal dysbiosis, but the chemical and pharmacological mechanisms remain poorly understood. In this study, we employed pharmacological network analysis to reveal mechanisms underlying the therapeutic effect of CC and MF against CRC. All compounds and targeted genes were obtained from the traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP). Differentially expressed genes (DEGs) were identified based on GSE146587, GSE156720, and GSE184093 datasets. A protein-protein interaction (PPI) network was constructed to identify putative target genes of CC and MF. Ten key targeted genes were identified, including CCND1, ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, MYC, SERPINE1, and VEGFA. Among these genes, six (ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, and SERPINE1) were positively correlated with levels of effector memory CD4 T cells and natural killer T cells, and three (CCND1, MYC, and VEGFA) were negatively correlated with type 17 T helper cells and CD56dim natural killer cells. Molecular docking analysis showed that four compounds of CC and MF (kaempferol, oleanolic acid, quercetin, and ursolic acid) could affect CRC by interacting with target genes. Our study proved that pharmacological analysis could reliably assess the mechanism of traditional Chinese medicines for treating cancer.
Colorectal cancer (CRC) is the third most common cancer worldwide, accounting for 10% of all cancers . Unhealthy lifestyles and diet are risk factors for CRC , such as heavy drinking, tobacco smoking, and red or processed meat consumption, while supplementing calcium, whole grain, fiber, and milk products can diminish the risk . Common symptoms of CRC include altered bowel habits, hematochezia, iron-deficiency anemia, weight loss, and stomachache . Colonoscopy is considered the optimal method with the highest sensitivity and specificity among screening methods for CRC, which can detect and resect most cases of early CRC. However, colonoscopy is an invasive procedure that can cause patient discomfort and potential complications, limiting its widespread use in physical examinations .
Fortunately, the overall survival of advanced CRC has largely improved due to the advent of novel chemotherapeutic drugs and targeted therapies combining multidisciplinary therapeutics . Nevertheless, these effective therapies have severe side effects, such as stomachache, constipation, and diarrhea. Therefore, new effective and safe therapeutics are still necessary. Traditional Chinese medicine (TCM) is a popular alternative to western medicine that has been applied in Asia for over 2000 years. The efficiency, abundant source, and the low incidence of side effects account for TCM’s increasing popularity in other countries. In clinical cancer therapy, TCM involves treatment via prescriptions or serves as adjuvant therapy with conventional chemotherapy [7–9]. The methanol extracts of lead bioactive isolates of Dita bark are reportedly an effective source of antidepressants, anti-inflammatory drugs, and thrombolytic agents . Besides, the methanol extracts of Hedychium coronarium (Zingiberaceae) have functions against oxidization, enzyme activity, and tumors , while those of Rhamnus triquetra (Wall.) are candidate components to develop new agents against multiple chronic diseases . Additionally, the methanol extracts of Psychotria calocarpa leaves (MEPC) have neuropharmacological, analgesic, antidiarrheal, and antioxidization properties ; Spirulina platensis has great therapeutic potential against diseases involving inflammatory pain . Cannabinoid receptors CB1 and CB2 and type 2 dopaminergic receptor (D2) exhibit antidepressive effects in models of depression and inflammation . Ophiorrhiza rugosa var. prostrata is highly valued at the medicinal level in treating all kinds of pain conditions . Aglaonema hookerianum Schott is recommended as an effective source of neuroprotective agents and as a candidate agent for sexual desire enhancement in cases of neural and sexual disorders . Tabebuia pallida (Lindl.) Miers (T. pallida) have been reported to exert antitumor effects via induction of apoptosis . Gynura possesses multiple pharmacological potentials against diabetes, oxidization, inflammation, bacterial infection, hypertension, and tumors . Andrographis paniculata (Burm. f.) is capable of impeding invasive microorganisms and inhibiting biofilm formation . Besides, Lepidagathis hyalina Nees are potential sources of antioxidants, cytotoxic agents, thrombolytic agents, antianxiety agents, and antidepressants . Moreover, the seeds of Syzygium fruticosum (SF) have health-promoting effects and are the optimal candidate agents used to prevent degenerative diseases . Green Tea Polyphenols are involved in immunoregulation and have antibacterial, antioxidant and anti-inflammatory activities, and it is used against diseases caused by the coronavirus . The methanol extracts of Molineria capitulata have antioxidization, cytotoxicity, thrombolysis, anti-inflammation, and analgesic properties . Last but not least, cannabinoids can be used as a source of new functional food ingredients and nutraceuticals, which are effective in treating and preventing gastrointestinal diseases . Taken together, TCM herbs and their extracts play a huge role in pharmacological studies.
Coptis chinensis (CC) belongs to Ranunculaceae, composed of berberine, coptisine, Coptis chinensis Polysaccharide, quercetin, and so on . CC has been widely used for treating bacillary dysentery, diabetes, pertussis, aphtha, eczema, and sore throat. Mume Fructus (MF) is the fruit of Prunus mume belonging to Rosaceae, generally cultivated in various regions in China, with the greatest output in the provinces south of the Yangtze River. Its major components include eugenol, dodecylic acid, ursolic acid, and quercetin. MF is commonly used for treating cough, bacillary dysentery, and sore throat in China. A combination of CC and MF is generally applied for treating ulcerative colitis, chronic diarrhea, and CRC. Although combined therapy of CC and MF is commonly used in clinical medicine, the active compounds and mechanism underlying their therapeutic effect against CRC remain unclear.
CC and MF combination is a multitarget and multicomponent formulation and is effective against CRC by regulating the intracellular molecular network. To investigate the specific mechanism of CC and MF in treating CRC, we applied bioinformatics and a multi-omics approach involving network pharmacology that has been established as important strategies for understanding TCM [27–30]. For example, Song et al. applied pharmacologic analysis to interpret the mechanism of triptolide in rheumatoid arthritis . Zhang et al. revealed that Bushen Tiansui Formula positively affected cognitive improvement in Alzheimer’s disease by using integrated analysis of pharmacology and serum metabolomics . Yang et al. revealed the pharmacological mechanism of Rhizoma Atractylodis macrocephalae in treating gastritis based on pharmacological analysis . The present study aimed to understand the mechanism of CC and MF treating CRC by using pharmacologic analysis and utilizing a series of CRC expression data from public databases. The workflow of the study is shown in Figure 1.
2. Materials and Methods
2.1. Analysis Based on the Pharmacologic Network
2.1.1. Screening of Active Components in CC and MF
All compounds in CC and MF were obtained from the traditional Chinese medicines integrated database (TCMID, http://www.megabionet.org/tcmid/) . The molecular weight (MW) of compounds is widely acknowledged to affect the efficiency of cell absorption; a high MW is difficult to absorb by cells, and compounds with WM between 180 and 500 are more suitable to be used as medicines, according to the Lipinski principle. Peroral drugs can function only after undergoing absorption, dispersion, metabolism, and excretion (ADME). During the process of ADME, drug-likeness (DL) is an important indicator for determining the active compounds . Therefore, active compounds in CC and MF were included based on the screening criteria of MW ≤ 500 and DL ≥ 0.18.
2.1.2. Constructing Drug-Target Interaction Networks
Protein targets of active compounds in CC and MF were acquired from TCMID and TCMSP databases. Cytoscape (v3.9)  was applied to construct CC/MF-active compounds-protein targets for visualizing the interaction between active compounds of CC/MF and their targets. In the network, drugs, compounds, and protein targets were all indicated as nodes, and their interactions were indicated as lines.
2.1.3. Gene Set Enrichment Analysis on Drug Targets
To analyze gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of drug targets, the R package clusterProfiler (v4.0)  was employed. Functional terms were screened using the criteria count ≥2 and the expression analysis systematic explorer (EASE) score ≤0.05.
2.2. Analyzing Chip Data of CRC
2.2.1. CRC Chip Data
CRC expression profiles of chip data were obtained from gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). Three chips with CRC expression profiles were searched, including GSE156720, GSE146587, and GSE184093. The GSE156720 dataset based on the GPL26963 platform included the expression data of normal samples (n = 3) and CRC samples (n = 3). The GSE146587 dataset based on the GPL17077 platform included expression data of normal samples (n = 6) and stage III CRC samples (n = 6). The GSE184093 dataset based on the GPL20115 platform included data on normal samples (n = 9) and CRC samples (n = 9).
2.2.2. Differential Analysis between Normal and Tumor Samples
The R package limma  was used to perform differential analysis between normal and tumor samples in the three datasets. A paired t-test was conducted to analyze data with normal distribution; otherwise, the Wilcoxon test was conducted to examine the significance of expression profiles between normal and tumor samples. and |fold change (FC)| > 2 were used as screening criteria to identify differentially expressed genes (DEGs). Then, the pheatmap R package (v1.0.12) was employed to output the heatmap of hierarchical clustering and scatter plots to display the expression level.
2.2.3. Identification of Key Drug Targets
The intersection between drug targets and DEGs identified in GSE cohorts was selected and analyzed in the STRING online tool (https://cn.string-db.org/). The results from STRING were visualized by Cytoscape (v3.9). The degree between each interaction was calculated, and top 10 genes with the highest degrees were considered key genes or drug targets.
2.3. Analysis of TCGA Data
2.3.1. Acquisition of TCGA Data and Data Preprocessing
Samples of colorectal cancer were obtained from The Cancer Genome Atlas (TCGA) database (https://tcga-data.nci.nih.gov/). The TCGA-COAD dataset included 471 COAD samples and 41 normal samples, and the TCGA-READ dataset included 167 READ samples and 10 normal samples, where mRNA expression data and gene mutation data were included. The requirement for ethical approval was waivered by the ethics committee. The study met the standards of TCGA for publication. Genes in expression profiles were annotated based on Ensembl gene ID. Gene expression was transferred to log2 expression and normalized through NMM normalization in the limma package. The average expression was selected when multiple expression data of one gene were presented. Genes with averaged expression >1 were retained, and low expressed data was eliminated.
2.3.2. Methodologies for Analyzing TCGA Data
Correlation analysis was assessed by the corrplot package (v0.92). Gene mutation analysis was conducted by maftools (v2.10.05).
2.3.3. Characterizing Tumor Microenvironment (TME)
Single-sample gene set enrichment analysis (ssGSEA)  was used to calculate the proportion of immune cells based on a series of gene sets [40, 41]. To avoid deviation caused by tumor purity, the ESTIMATE algorithm was performed to adjust the enrichment score of immune cells . A total of 24 cell types in TME were analyzed, including endothelial cells, eosinophils, fibroblasts, M0 macrophages, M1 macrophages, M2 macrophages, activated mast cells, resting mast cells, monocytes, neutrophils, memory B cells, naive B cells, activated dendritic cells, resting dendritic cells, activated NK cells, resting NK cells, plasma cells, activated CD4 memory T cells, resting CD4 memory T cells, naive CD4 T cells, CD8 T cells, follicular helper T cells, gamma delta T cells, and regulatory T cells. The expression of four immune checkpoints (PD-L1, CTLA4, PD-1, and PD-L2) was evaluated.
2.4. Molecular Docking for Active Compounds and Drug Targets
To confirm the interaction between key active compounds and protein targets, molecular docking was implemented. Proteins corresponding to drug targets related to CRC were obtained from protein data bank (PDB) (https://www.rcsb.org/pdb/home/home) . Then, the AutoDocking Vina (v1.1.2) software [44, 45] was applied to modify the crystal structure, including adding hydrogen, removing hydrate and ligands, and optimizing and supplementing amino acids. A grid module was employed to implement molecular docking, and binding energy was used as a standard for evaluating the stability of molecular docking. The binding stability was strengthened by increasing the binding energy score.
3.1. Analysis Based on the Pharmacologic Network
We obtained 87 active compounds of CC and MF from TCMSP and TCMID databases, with 48 belonging to CC and 40 to MF. One compound (quercetin) was found in both CC and MF. 56.25% (27/48) and 45.00% (18/40) of active compounds in CC and MF met the criteria of MW ≤ 500 and DL ≥ 0.18. Finally, 44 active compounds were retained for the following analysis. 196 protein targets of 44 active compounds were identified from the TCMSP database. A protein-protein interaction network was constructed based on 196 protein targets and 44 compounds, where 242 nodes were included with 2 drugs, 44 active compounds, and 196 protein targets (Figure 2(a)). Notably, four key compounds of CC and MF were identified with high degrees, including quercetin (Degree = 152), (R)-Canadine (Degree = 32), beta-sitosterol (Degree = 28), and magnoflorine (Degree = 24).
GO analysis revealed the top 10 enriched molecular function terms of CC and MF, including “response to xenobiotic stimulus,” “response to lipopolysaccharide,” “response to oxygen levels,” “response to decreased oxygen levels,” and “reactive oxygen species metabolic process” (Figure 2(b)). KEGG analysis showed that some tumor-related pathways were significantly enriched, such as “chemical carcinogenesis-receptor activation” and “IL-17 signaling pathway” (Figure 2(c)). The above results suggested that CC/MF formulation exhibited multi-target and multicomponent characteristics for treatment via regulating cellular networks.
3.2. Differential Analysis between Normal and CRC Samples
We compared expression data between normal and tumor samples in three datasets (GSE156720, GSE146587, and GSE184093) and identified a series of DEGs by comparing tumor samples to normal samples (Figures 3(a)–3(f)). In the GSE156720 dataset, 832 upregulated and 643 downregulated genes were detected. In the GSE146587 dataset, 1426 upregulated and 1385 downregulated genes were detected. Finally, in the GSE184093 dataset, 1040 upregulated and 1300 downregulated genes were detected. The intersection of DEGs among three datasets yielded 714 upregulated DEGs, and 706 downregulated DEGs in at least two datasets (Figures 3(g) and 3(h)).
Then we assessed the enriched GO terms and KEGG pathways of upregulated and downregulated DEGs, respectively. For upregulated genes, molecular function terms, such as “organelle fission,” “mitotic nuclear division,” “DNA replication,” and “regulation of mitotic cell cycle,” were significantly enriched (, Figure 3(i)). Immune-related pathways such as “cytokine-cytokine receptor interaction” and tumor-related pathways, such as “cell cycle,” “TNF signaling pathway,” and “p53 signaling pathway,” were significantly enriched (, Figure 3(j)), indicating that these upregulated genes were highly involved in tumor proliferation. For downregulated genes, “cellular divalent inorganic cation homeostasis,” “cellular hormone metabolic process,” and “regulation of metal ion transport” were significantly enriched (, Figure 3(k)). Enriched KEGG pathways in downregulated genes were related to cell metabolism, such as “cAMP signaling pathway,” “regulation of lipolysis in adipocytes,” and “pentose and glucuronate interconversions” (, Figure 3(l)).
3.3. Identifying Key Drug Targets Based on PPI Network Analysis
In the previous section, we identified DEGs that were CRC-associated genes. The intersection between potential target genes of CC/MF and CRC-associated genes yielded 28 target genes among upregulated DEGs and 15 target genes in downregulated genes (Figure 4(a)). Then these intersected genes were further analyzed by using the STRING online tool, and the degree of each node was calculated (Figure 4(b)). Top 10 target genes were determined as key drug targets, including CCND1, ICAM1 IL1B, IL-6, MMP1, MMP3, MMP9, MYC, SERPINE1, and VEGFA (Figure 4(c)).
3.4. Analyzing the Features of Key Genes in the TCGA-COAD Dataset
To verify the important role of these 10 key genes, we further characterized their features in two independent datasets (TCGA-COAD and TCGA-READ). We first evaluated them in the TCGA-COAD dataset. Ten key genes (drug targets) were significantly overexpressed in the TCGA-COAD dataset (Figure 4(d)), consistent with the result in GSE cohorts. Correlation analysis revealed positive correlations among these key genes, especially between MMP1 and MMP3, and IL-6 and MMP3 (R > 0.5). By analyzing TCGA mutant gene data, 56 out of 399 samples were found to have mutations in these 10 genes (Figure 4(f)). MMP9 exhibited the highest mutation with a mutation frequency of 5% (Figure 4(f)). Missense mutations accounted for most mutations, with C > T commonly observed.
Then we evaluated the relationship between key genes and TME by analyzing their correlation with immune checkpoints and immune cells. Analysis of the relationship with four important immune checkpoints (PD-L1, CTLA4, PD-1, and PD-L2) in TME showed that ICAM1 expression was strikingly associated with the expression of all four immune checkpoints (, Figure 4(g)). In addition, MMP1 and MMP9 were significantly correlated with CTLA4 and PD-L2 expression (). Furthermore, we calculated the enrichment score of 28 immune-related cells by ssGSEA and analyzed the correlation between the expression of key genes and the enrichment score of these cells (Figure 4(h)). The results showed that ICAM1 was positively correlated with most immune cells, especially natural killer T cells, natural killer cells, and memory CD8 T cells. Notably, only ICAM1 was correlated with activated CD8 T cells, indicating its critical interaction with T cell activation. 7 of 10 genes (ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, and SERPINE1) were positively correlated with memory CD8 T cells, memory CD4 T cells, regulatory T cells, natural killer T cells, and so on. Additionally, CCND1, MYC, and VEGFA were negatively correlated with most immune cells, which suggested that they may play an unfavorable role in immunotherapy.
3.5. Analyzing Molecular Features of Key Genes in the TCGA-READ Dataset
We validated 10 key genes in the TCGA-READ dataset and observed consistent results. Ten key genes were overexpressed in READ samples compared to normal expression, although 3 (IL1B, IL-6, and MMP9) showed no significant difference (Figure 5(a)). Stronger correlations among 10 key genes were observed in the TCGA-READ dataset compared to the TCGA-COAD dataset (Figures 5(b) and 4(e)). Mutational analysis showed that only 6 of 137 samples had mutations of 10 key genes (Figure 5(c)). ICAM1 exhibited the strongest correlation with immune checkpoints, especially PD-L2 and PD-1 (Figure 5(d)). According to the TCGA-COAD dataset results, 7 of 10 key genes were positively correlated with most immune cells, and ICAM1 exhibited the highest correlation (Figure 5(e)). VEGFA, MYC, and CCND1 were negatively correlated with most immune cells, suggesting that they may play a negative role in immunotherapy.
Overall, 10 key genes displayed similar performance in TCGA-COAD and TCGA-READ datasets. Few mutations were found in these 10 genes, indicating that these genes had high specificity in COAD and READ. In immune analysis, ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, and SERPINE1 all showed positive correlations with most immune cells, implying that CC and MF may exert an antitumor activity by interacting with these genes and immune cells, and thus control tumor development.
3.6. Molecular Docking between Active Compounds and Protein Targets
According to the interaction among drugs, compounds, and protein targets, we identified 10 key genes (protein targets) that were putative targets of active compounds. Then we applied the AutoDocking Vina software to conduct molecular docking between active compounds and protein targets (Figures 6(a)–6(j)). The predicted binding energy of 10 protein targets is shown in Table 1. Generally, low binding energy indicates stable binding and a high possibility of interaction. Quercetin was able to interact with all protein targets, and ursolic acid could interact with most protein targets (Figure 6). Among 10 protein targets, ICAM1 could interact with all four active compounds, including kaempferol, oleanolic acid, quercetin, and ursolic acid. These results demonstrated that CC and MF exert an antitumor effect against CRC via the interaction of these active compounds with key protein targets that could modulate TME.
TCM holds an important place in Chinese pharmacological studies and clinical practice. In recent years, overwhelming evidence has substantiated that multiple TCM herbs and their extracts possess antitumor activities. Emodin is a type of anthraquinone derivative distributed in the roots and rhizosphere of diverse plants. It is a natural compound with antitumor activity, predominantly inhibiting cancer cell growth, and proliferation through weakening oncogenic growth signals, such as protein kinase B (AKT), mitogen-activated protein kinase (MAPK), HER-2, Wnt/-catenin, and phosphatidylinositol 3-kinase (PI3K) . Coumarin and related compounds have high bioactivity and low toxicity. They are commonly used in treating prostate cancer, renal cell carcinoma (RCC), and leukemia, with the ability to induce apoptosis, autophagy and cell cycle arrest, antioxidization, and metastasis through activation of various signaling pathways in malignant cells. In the meantime, they can counter the side effects caused by chemotherapy [47, 48]. Multiple natural products have demonstrated antitumor activities by promoting cell apoptosis and inhibiting metastasis to slow cancer spread while reducing cancer cell proliferation and triggering death . Some natural compounds have been shown to have the potential to act as aromatase inhibitors, and thus play a huge role in the treatment of gynecological tumors [50, 51]. TCM extracts also play important roles in studies on colorectal cancer. For example, cruciferous vegetables and their bioactive metabolites play antiviral and antibacterial roles via protecting cells from DNA damage and deactivating oncogenic substances. Additionally, they possess antitumor effects presenting with suppression of metastasis and generation of vessels that supply blood to the tumor through inducing cell apoptosis and inhibiting migration . Flavonoid compounds, with antioxidant, anti-inflammatory, antibacterial, antifungal, and antitumor activities, are highly effective in alleviating clinical symptoms associated with diarrhea, mucositis, neuropathic pain, and chemotherapy for colorectal cancer .
In the present study, we applied pharmacologic network analysis to predict the possible mechanism of CC and MF in treating CRC. GSE cohorts were used to screen key genes, and TCGA datasets were used as independent datasets to validate key genes. Molecular docking was finally performed to confirm the molecular interaction between active compounds and key protein targets.
CC has been used to treat various inflammatory diseases for over a thousand years, to clear away heat and toxic materials and dry dampness, according to Chinese medicine theories. The medicinal value of CC was first documented in ShenMungHerbal, and over 32000 kinds of Chinese medicine formula, including CC in the form of pulvis, pills, decoction, and tablets, can treat diarrhea, emesis, abdominal distension, hyperpyrexia, toothache, diabetes, jaundice, and eczema . Recent evidence suggests that CC has an extensive pharmacological activity against bacteria, viruses, hepatic steatosis, atherosclerosis, reperfusion injury after myocardial ischemia, diabetes, arrhythmia, hypertension, inflammation, and tumors [54, 55]. To date, over 120 compounds have been identified and extracted from CC; most compounds belong to alkaloids, followed by lignans, flavonoids, organic acids, and volatile oil. Interestingly, a study demonstrated that CC has a therapeutic effect against inflammatory bowel disease through inhibiting oxidative stress, protecting intestinal epithelium’s barrier, alleviating pain, regulating T helper cells, and antimicrobial activity . In addition, polysaccharides in CC can regulate immune responses related to the intestinal epithelium and intestinal microenvironment in an effective, dynamic, and dose-independent manner . Moreover, epiberberine in CC is considered a potential drug for treating MKN-45-associated gastric cancer, targeting the p53-independent mitochondrial pathway .
MF has been used in traditional medicine for over 3000 years for treating fatigue, cough, headache, constipation, gastric diseases, and food poisoning. Growing evidence suggests that the major components of MF, including flavonoids, sterols, and organic acids, have multiple pharmacologic functions, such as bacteriostasis, antitumor activity, and regulating gut microbiota. Citric acid, chlorogenic acid, and neochlorogenic acid in MF can significantly alleviate cough and airway inflammation, inhibit the overproduction of mucus, and thus reduce chronic cough resulting from tobacco smoke in a guinea pig model . Xing et al. illustrated that MF extracts could alleviate diarrhea caused by the combined treatment of lapatinib and Xeloda for breast cancer patients . In addition, MF aqueous extracts can inhibit the expression of iNOS and COX-2, and the generation of PGE2, NO, and IL-6 induced by lipopolysaccharide (LPS) through blocking IkappaBalpha degradation and mitogen-activated protein kinase (MAPK) phosphorylation following with NF-kappaB activation .
Functional analysis of dysregulated genes in CRC samples showed that cell cycle-related processes “organelle fission,” “mitotic nuclear division,” “DNA replication,” and “regulation of mitotic cell cycle,” tumor-related pathways, such as “cytokine-cytokine receptor interaction,” and immune-related pathways, such as “Cell cycle,” “TNF signaling pathway,” and “p53 signaling pathway,” were significantly enriched, implying that these dysregulated pathways may be important mechanisms of CRC progression. We identified 10 key genes (CCND1, ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, MYC, SRTPINE1, and VEGFA) as putative targets of CC and MF during treatment against CRC. Some key genes have been reported to be associated with CRC development. In this regard, CCND1 has been documented to regulate the cell cycle in CRC tumorigenesis and development . The MMP protein family promotes the malignancy of CRC cells through epithelial-mesenchymal transition (EMT) and Akt pathways . Moreover, IL-6/IL-11 activates STAT3 in cancer-associated fibroblasts (CAFs) to promote CRC progression and poor prognosis . Finally, a study found that the c-Myc/miR-27b-3p/ATG10 signaling pathway is involved in the drug resistance of chemotherapy in CRC .
Molecular docking demonstrated that the 10 protein targets could interact with four active compounds (kaempferol, oleanolic acid, quercetin, and ursolic acid) of CC and MF, which revealed the potential mechanism of CC and MF for treating CRC. Ten key genes identified in the study may serve as potential therapeutic targets for treating CRC. Furthermore, the study provides new directions for understanding CRC development based on CC and MF pharmacological network analysis. Taken together, our findings lay the groundwork for future pharmacological research and clinical applications, and provide an overview for the development of new drugs for colorectal cancer.
The present study combined network pharmacology, RNA sequencing, gene mutation, immune cell abundance, and molecular docking comprehensively discussed the mechanism underlying the therapeutic efficacy of the combination of CC with MF in CRC treatment, which involved suppression of inflammatory responses, cancer cell metabolism, and cell cycle process. Kaempferol, oleanolic acid, quercetin, and ursolic acid were identified as active compounds in CC and MF, while 10 key genes (including CCND1, ICAM1, IL1B, IL-6, MMP1, MMP3, MMP9, MYC, SRTPINE1, and VEGFA) were discovered as hub genes that played a mediating role in immunotherapy via triggering alterations of immune cells in the body, such as Natural killer T cell and memory CD8 T cell. During clinical practice, CC and MF are commonly used in combination to enhance immunity and inhibit cancer cell growth in cancer patients. The current study provided hitherto undocumented evidence of the mechanism of action of CC and MF in CRC treatment and broadened the therapeutic landscape for this patient population. However, our study still has some limitations. Indeed, no quantitative or experimental analysis of key active compounds was performed. In addition, only in silico research and limited CRC samples were used, warranting further validation in large-scale cell or animal experiments.
Data are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
R-JG and YL contributed to conceptualization, methodology, investigation, and validation; R-JG was responsible for visualization, formal analysis, software, data curation, writing—original draft; supervision, YL was responsible for writing—review and editing. All authors have read and agreed to the published version of the manuscript.
This study was supported by the general project of the Guangdong Bureau of Traditional Chinese Medicine: 20213009. In addition, the authors are indebted to Southern Medical University for providing convenience in the collection of documents.
H. Sung, J. Ferlay, R. L. Siegel et al., “Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 71, no. 3, pp. 209–249, 2021.View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, A. Goding Sauer et al., “Colorectal cancer statistics, 2020,” CA: A Cancer Journal for Clinicians, vol. 70, no. 3, pp. 145–164, 2020.View at: Publisher Site | Google Scholar
J. de Winter, N. P. M. Ezendam, M. J. L. Bours et al., “Is sleep associated with BMI, waist circumference, and diet among long-term colorectal cancer survivors? Results from the population-based PROFILES registry,” Support Care Cancer, vol. 29, no. 12, pp. 7225–7235, 2021.View at: Publisher Site | Google Scholar
K. Thanikachalam and G. Khan, “Colorectal cancer and nutrition,” Nutrients, vol. 11, no. 1, 2019.View at: Publisher Site | Google Scholar
B. A. Alves Martins, G. F. de Bulhões, I. N. Cavalcanti, M. M. Martins, P. G. de Oliveira, and A. M. A. Martins, “Biomarkers in colorectal cancer: the role of translational proteomics research,” Frontiers in Oncology, vol. 9, p. 1284, 2019.View at: Publisher Site | Google Scholar
L. Luo, Y. Ma, Y. Zheng, J. Su, and G. Huang, “Application progress of organoids in colorectal cancer,” Frontiers in Cell and Developmental Biology, vol. 10, Article ID 815067, 2022.View at: Publisher Site | Google Scholar
K.-C. Kim, J.-H. Yook, J. Eisenbraun, B.-S. Kim, and R. Huber, “Quality of life, immunomodulation and safety of adjuvant mistletoe treatment in patients with gastric carcinoma - a randomized, controlled pilot study,” BMC Complementary and Alternative Medicine, vol. 12, no. 1, p. 172, 2012.View at: Publisher Site | Google Scholar
L. Xu, H. Li, Z. Xu et al., “Multi-center randomized double-blind controlled clinical study of chemotherapy combined with or without traditional Chinese medicine on quality of life of postoperative non-small cell lung cancer patients,” BMC Complementary and Alternative Medicine, vol. 12, no. 1, p. 112, 2012.View at: Publisher Site | Google Scholar
K. Gavaraskar, S. Dhulap, and R. R. Hirwani, “Therapeutic and cosmetic applications of Evodiamine and its derivatives-A patent review,” Fitoterapia, vol. 106, pp. 22–35, 2015.View at: Publisher Site | Google Scholar
M. F. Khan, F. B. Kader, M. Arman et al., “Pharmacological insights and prediction of lead bioactive isolates of Dita bark through experimental and computer-aided mechanism,” Biomedicine & Pharmacotherapy, vol. 131, Article ID 110774, 2020.View at: Publisher Site | Google Scholar
K. I. Sinan, U. Akpulat, A. A. Aldahish et al., “LC-MS/HRMS analysis, anti-cancer, anti-enzymatic and anti-oxidant effects of boerhavia diffusa extracts: a potential raw material for functional applications,” Antioxidants, vol. 10, no. 12, 2021.View at: Publisher Site | Google Scholar
J. Iqbal, B. A. Abbasi, R. Ahmad et al., “Phytogenic synthesis of nickel oxide nanoparticles (NiO) using fresh leaves extract of Rhamnus triquetra (wall.) and investigation of its multiple in vitro biological potentials,” Biomedicines, vol. 8, no. 5, 2020.View at: Publisher Site | Google Scholar
T. A. Bristy, N. Barua, A. Montakim Tareq et al., “Deciphering the pharmacological properties of methanol extract of Psychotria calocarpa leaves by in vivo, in vitro and in silico approaches,” Pharmaceuticals (Basel, Switzerland), vol. 13, no. 8, 2020.View at: Publisher Site | Google Scholar
M. A. Freitas, A. Vasconcelos, E. C. D. Gonçalves et al., “Involvement of opioid system and TRPM8/TRPA1 channels in the antinociceptive effect of Spirulina platensis,” Biomolecules, vol. 11, no. 4, 2021.View at: Publisher Site | Google Scholar
G Vieira, J Cavalli, E. C. D. Gonçalves et al., “Antidepressant-like effect of terpineol in an inflammatory model of depression: involvement of the cannabinoid system and D2 dopamine receptor,” Biomolecules, vol. 10, no. 5, 2020.View at: Publisher Site | Google Scholar
M. N. Uddin Chy, M. Adnan, M. R. Chowdhury et al., “Central and peripheral pain intervention by Ophiorrhiza rugosa leaves: Potential underlying mechanisms and insight into the role of pain modulators,” Journal of Ethnopharmacology, vol. 276, Article ID 114182, 2021.View at: Publisher Site | Google Scholar
O. Goni, M. F. Khan, M. M. Rahman et al., “Pharmacological insights on the antidepressant, anxiolytic and aphrodisiac potentials of Aglaonema hookerianum Schott,” Journal of Ethnopharmacology, vol. 268, Article ID 113664, 2021.View at: Publisher Site | Google Scholar
M. M. Rahman, A. S. M. A. Reza, M. A. Khan et al., “Unfolding the apoptotic mechanism of antioxidant enriched-leaves of Tabebuia pallida (Lindl.) Miers in EAC cells and mouse model,” Journal of Ethnopharmacology, vol. 278, Article ID 114297, 2021.View at: Publisher Site | Google Scholar
M. S. Bari, L. Khandokar, E. Haque et al., “Ethnomedicinal uses, phytochemistry, and biological activities of plants of the genus Gynura,” Journal of Ethnopharmacology, vol. 271, Article ID 113834, 2021.View at: Publisher Site | Google Scholar
S. Hossain, Z. Urbi, H. Karuniawati et al., “Andrographis paniculata (burm. F.) wall. Ex nees: an updated review of phytochemistry, antimicrobial pharmacology, and clinical safety and efficacy,” Life (Basel), vol. 11, no. 4, 2021.View at: Google Scholar
F. I. Fahad, N. Barua, M. S. Islam et al., “Investigation of the pharmacological properties of Lepidagathis hyalina nees through experimental approaches,” Life (Basel), vol. 11, no. 3, 2021.View at: Publisher Site | Google Scholar
J. N. R. Moni, M. Adnan, A. M. Tareq et al., “Therapeutic potentials of Syzygium fruticosum fruit (seed) reflected into an array of pharmacological assays and prospective receptors-mediated pathways,” Life (Basel), vol. 11, no. 2, 2021.View at: Publisher Site | Google Scholar
T. E. Tallei, Fatimawali, N. J. Niode et al., “A comprehensive review of the potential use of green tea Polyphenols in the management of COVID-19,” Evidence-Based Complementary And Alternative Medicine: eCAM, vol. 2021, Article ID 7170736, 2021.View at: Publisher Site | Google Scholar
M. A. R. B. Shovo, M. R. Tona, J. Mouah et al., “Computational and pharmacological studies on the antioxidant, thrombolytic, anti-inflammatory, and analgesic activity of Molineria capitulata,” Current Issues in Molecular Biology, vol. 43, no. 2, pp. 434–456, 2021.View at: Publisher Site | Google Scholar
V. Martínez, A. Iriondo De-Hond, F. Borrelli, R. Capasso, M. D. Del Castillo, and R. Abalo, “Cannabidiol and other non-psychoactive cannabinoids for prevention and treatment of gastrointestinal disorders: useful nutraceuticals?” International Journal of Molecular Sciences, vol. 21, no. 9, 2020.View at: Google Scholar
J. Wang, L. Wang, G.-H. Lou et al., “Coptidis Rhizoma: a comprehensive review of its traditional uses, botany, phytochemistry, pharmacology and toxicology,” Pharmaceutical Biology, vol. 57, no. 1, pp. 193–225, 2019.View at: Publisher Site | Google Scholar
T. T. Luo, Y. Lu, S. K. Yan, X. Xiao, X. L. Rong, and J. Guo, “Network pharmacology in research of Chinese medicine formula: methodology, application and prospective,” Chinese Journal of Integrative Medicine, vol. 26, no. 1, pp. 72–80, 2020.View at: Publisher Site | Google Scholar
S. Li and B. Zhang, “Traditional Chinese medicine network pharmacology: theory, methodology and application,” Chinese Journal of Natural Medicines, vol. 11, no. 2, pp. 110–120, 2013.View at: Publisher Site | Google Scholar
R. Zhang, X. Zhu, H. Bai, and K. Ning, “Network pharmacology databases for traditional Chinese medicine: review and assessment,” Frontiers in Pharmacology, vol. 10, p. 123, 2019.View at: Publisher Site | Google Scholar
H. Yuan, Q. Ma, H. Cui et al., “How can synergism of traditional medicines benefit from network pharmacology?” Molecules, vol. 22, no. 7, 2017.View at: Publisher Site | Google Scholar
X. Song, Y. Zhang, E. Dai, L. Wang, and H. Du, “Prediction of triptolide targets in rheumatoid arthritis using network pharmacology and molecular docking,” International Immunopharmacology, vol. 80, Article ID 106179, 2020.View at: Publisher Site | Google Scholar
Z. Zhang, P. Yi, J. Yang et al., “Integrated network pharmacology analysis and serum metabolomics to reveal the cognitive improvement effect of Bushen Tiansui formula on Alzheimer’s disease,” Journal of Ethnopharmacology, vol. 249, Article ID 112371, 2020.View at: Publisher Site | Google Scholar
S. Yang, J. Zhang, Y. Yan et al., “Network pharmacology-based strategy to investigate the pharmacologic mechanisms of atractylodes macrocephala koidz. For the treatment of chronic gastritis,” Frontiers in Pharmacology, vol. 10, p. 1629, 2019.View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014.View at: Publisher Site | Google Scholar
X. Xu, W. Zhang, C. Huang et al., “A novel chemometric method for the prediction of human oral bioavailability,” International Journal of Molecular Sciences, vol. 13, no. 6, pp. 6964–6982, 2012.View at: Publisher Site | Google Scholar
D. Otasek, J. H. Morris, J. Bouças, A. R. Pico, and B. Demchak, “Cytoscape Automation: empowering workflow-based network analysis,” Genome Biology, vol. 20, no. 1, p. 185, 2019.View at: Publisher Site | Google Scholar
T. Wu, E. Hu, S. Xu et al., “clusterProfiler 4.0: a universal enrichment tool for interpreting omics data,” The Innovation, vol. 2, no. 3, Article ID 100141, 2021.View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, 2015.View at: Publisher Site | Google Scholar
D. A. Barbie, P. Tamayo, J. S. Boehm et al., “Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1,” Nature, vol. 462, no. 7269, pp. 108–112, 2009.View at: Publisher Site | Google Scholar
A. M. Newman, C. L. Liu, M. R. Green et al., “Robust enumeration of cell subsets from tissue expression profiles,” Nature Methods, vol. 12, no. 5, pp. 453–457, 2015.View at: Publisher Site | Google Scholar
P. Charoentong, F. Finotello, M. Angelova et al., “Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade,” Cell Reports, vol. 18, no. 1, pp. 248–262, 2017.View at: Publisher Site | Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.View at: Publisher Site | Google Scholar
S. K. Burley, H. M. Berman, G. J. Kleywegt, J. L. Markley, H. Nakamura, and S. Velankar, “Protein Data Bank (PDB): the single global macromolecular structure archive,” Methods in Molecular Biology, vol. 1607, pp. 627–641, 2017.View at: Publisher Site | Google Scholar
O. Trott, A. J. Olson, and A. D. Vina, “Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, pp. 455–461, 2010.View at: Google Scholar
R. Rolta, D. Salaria, V. Kumar et al., “Molecular docking studies of phytocompounds of Rheum emodi Wall with proteins responsible for antibiotic resistance in bacterial and fungal pathogens: in silico approach to enhance the bio-availability of antibiotics,” Journal of Biomolecular Structure and Dynamics, vol. 40, no. 8, pp. 3789–3803, 2022.View at: Publisher Site | Google Scholar
E. K. Akkol, I. I. Tatlı, G. Ş. Karatoprak et al., “Is emodin with anticancer effects completely innocent? Two sides of the coin,” Cancers, vol. 13, no. 11, 2021.View at: Publisher Site | Google Scholar
E. Küpeli Akkol, Y. Genç, B. Karpuz, E. Sobarzo-Sánchez, and R. Capasso, “Coumarins and coumarin-related compounds in pharmacotherapy of cancer,” Cancers, vol. 12, no. 7, 2020.View at: Google Scholar
S. Ahmed, H. Khan, M. Aschner, H. Mirzae, E. Küpeli Akkol, and R. Capasso, “Anticancer potential of furanocoumarins: mechanistic and therapeutic aspects,” International Journal of Molecular Sciences, vol. 21, no. 16, 2020.View at: Publisher Site | Google Scholar
M. R. Islam, F. Islam, M. H. Nafady et al., “Natural small molecules in breast cancer treatment: understandings from a therapeutic viewpoint,” Molecules (Basel, Switzerland), vol. 27, no. 7, 2022.View at: Publisher Site | Google Scholar
M. Balunas, B. Su, R. Brueggemeier, and A. Kinghorn, “Natural products as aromatase inhibitors,” Anti-Cancer Agents in Medicinal Chemistry, vol. 8, no. 6, pp. 646–682, 2008.View at: Publisher Site | Google Scholar
S. Mitra, M. S. Lami, A. Ghosh et al., “Hormonal therapy for gynecological cancers: how far has science progressed toward clinical applications?” Cancers, vol. 14, no. 3, 2022.View at: Publisher Site | Google Scholar
D. Ağagündüz, T. Ö. Şahin, B. Yılmaz, K. D. Ekenci, Ş. Duyar Özer, and R. Capasso, “Cruciferous vegetables and their bioactive metabolites: from prevention to novel therapies of colorectal cancer,” Evidence-Based Complementary And Alternative Medicine: eCAM, vol. 2022, Article ID 1534083, 2022.View at: Publisher Site | Google Scholar
J. Fernández, B. Silván, R. Entrialgo-Cadierno et al., “Antiproliferative and palliative activity of flavonoids in colorectal cancer,” Biomedicine & Pharmacotherapy, vol. 143, Article ID 112241, 2021.View at: Google Scholar
B.-L. Ma and Y.-M. Ma, “Pharmacokinetic properties, potential herb-drug interactions and acute toxicity of oral Rhizoma coptidisalkaloids,” Expert Opinion on Drug Metabolism & Toxicology, vol. 9, no. 1, pp. 51–61, 2013.View at: Publisher Site | Google Scholar
Z. Wang, Y. Yang, M. Liu et al., “Rhizoma coptidis for Alzheimer’s disease and vascular dementia: a literature review,” Current Vascular Pharmacology, vol. 18, no. 4, pp. 358–368, 2020.View at: Publisher Site | Google Scholar
Y. Yang, C. T. Vong, S. Zeng et al., “Tracking evidences of Coptis chinensis for the treatment of inflammatory bowel disease from pharmacological, pharmacokinetic to clinical studies,” Journal of Ethnopharmacology, vol. 268, Article ID 113573, 2021.View at: Publisher Site | Google Scholar
Q. Chen, R. Ren, Q. Zhang et al., “Coptis chinensis Franch polysaccharides provide a dynamically regulation on intestinal microenvironment, based on the intestinal flora and mucosal immunity,” Journal of Ethnopharmacology, vol. 267, Article ID 113542, 2021.View at: Publisher Site | Google Scholar
M. Yu, L. Ren, F. Liang et al., “Effect of epiberberine from Coptis chinensis Franch on inhibition of tumor growth in MKN-45 xenograft mice,” Phytomedicine, vol. 76, Article ID 153216, 2020.View at: Publisher Site | Google Scholar
J. Xiang, X. Liu, S. Zhong et al., “Fructus mume protects against cigarette smoke induced chronic cough Guinea pig,” Journal of Medicinal Food, vol. 23, no. 2, pp. 191–197, 2020.View at: Publisher Site | Google Scholar
H. Xing, L. Zhang, J. Ma, Z. Liu, C. Song, and Y. Liu, “Fructus mume extracts alleviate diarrhea in breast cancer patients receiving the combination therapy of lapatinib and capecitabine,” Frontiers in Pharmacology, vol. 9, p. 516, 2018.View at: Publisher Site | Google Scholar
H.-J. Choi, O.-H. Kang, P.-S. Park et al., “Mume Fructus water extract inhibits pro-inflammatory mediators in lipopolysaccharide-stimulated macrophages,” Journal of Medicinal Food, vol. 10, no. 3, pp. 460–466, 2007.View at: Publisher Site | Google Scholar
C.-Y. Huang, C.-W. Tsai, C.-M. Hsu, W.-S. Chang, H.-A. Shui, and D.-T. Bau, “The significant association of CCND1 genotypes with colorectal cancer in Taiwan,” Tumor Biology, vol. 36, no. 8, pp. 6533–6540, 2015.View at: Publisher Site | Google Scholar
K. Wang, J. Zheng, J. Yu et al., “Knockdown of MMP-1 inhibits the progression of colorectal cancer by suppressing the PI3K/Akt/c-myc signaling pathway and EMT,” Oncology Reports, vol. 43, no. 4, pp. 1103–1112, 2020.View at: Google Scholar
C. Heichler, K. Scheibe, A. Schmied et al., “STAT3 activation through IL-6/IL-11 in cancer-associated fibroblasts promotes colorectal tumour development and correlates with poor prognosis,” Gut, vol. 69, no. 7, pp. 1269–1282, 2020.View at: Publisher Site | Google Scholar
W. Sun, J. Li, L. Zhou et al., “The c-Myc/miR-27b-3p/ATG10 regulatory axis regulates chemoresistance in colorectal cancer,” Theranostics, vol. 10, no. 5, pp. 1981–1996, 2020.View at: Publisher Site | Google Scholar