Exploring the Antitumor Mechanisms of Zingiberis Rhizoma Combined with Coptidis Rhizoma Using a Network Pharmacology Approach

Wang, Meng; Qi, Youke; Sun, Yongning

doi:https://doi.org/10.1155/2020/8887982

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8887982 | https://doi.org/10.1155/2020/8887982

Exploring the Antitumor Mechanisms of Zingiberis Rhizoma Combined with Coptidis Rhizoma Using a Network Pharmacology Approach

Meng Wang,¹Youke Qi,²and Yongning Sun^1,2

Academic Editor: Rui Liu

Received30 Sept 2020

Revised30 Oct 2020

Accepted25 Nov 2020

Published28 Dec 2020

Abstract

Background. Although the combination of Zingiberis rhizoma (ZR) and Coptidis rhizoma (CR) is a classic traditional Chinese medicine-based herbal pair used for its antitumor effect, the material basis and underlying mechanisms are unclear. Here, a network pharmacology approach was used to elucidate the antitumor mechanisms of ZR-CR. Materials and Methods. To predict the targets of ZR-CR in treating tumors, we constructed protein–protein interactions and hub component-target networks and performed pathway and process enrichment and molecular docking analysis. We used a surface plasmon resonance (SPR) assay to validate the predicted component-target affinities. Hub gene expression and survival analysis in patients with tumors were used to predict the clinical significance. Results. The active components of ZR-CR—shogaol, daucosterol, ginkgetin, berberine, quercetin, chlorogenic acid, and vanillic acid—exhibited antitumor activities via the MAPK, PI3K-AKT, TNF, FOXO, HIF-1, and VEGF signaling pathways. Molecular docking and SPR analyses suggested direct binding of berberine with AKT1 and TP53; quercetin with EGFR and VEGF165; and ginkgetin, isoginkgetin, and daucosterol with VEGF165 with weak affinities. Gene expression levels of the hub targets of ZR-CR were associated with overall survival and disease-free survival in patients with various tumor types. Conclusions. The antitumor components of the ZR-CR herbal pair and the mechanisms underlying their antitumor effects were identified. These antitumor components deserve to be explored further in experimental and clinical studies.

1. Introduction

Although considerable efforts have been made to improve tumor diagnosis and treatment, tumors are still a serious threat to human health. Tumors are recognized as an important cause of morbidity and mortality worldwide; in nearly 100 countries around the world, cancer is the first or second leading cause of premature death [1, 2].

Traditional Chinese medicine (TCM), as an effective treatment method, has been used to treat patients with various clinical diseases, including cancer. TCM has developed over thousands of years because of its curative effects and unique diagnosis and treatment system [3]. The combination of Zingiberis rhizoma (ZR) and Coptidis rhizoma (CR) is a classic TCM-based herb pair, which is recorded in “Treatise on Cold Pathogenic and Miscellaneous Diseases” (“Shang Han Za Bing Lun” in Chinese). This herb pair was mainly used for digestive system diseases according to ancient records but it is currently also used for its antitumor effects. However, the material basis and mechanisms of the antitumor actions of ZR-CR still need to be further explored.

A growing number of theories and practices show that network pharmacology can clarify the potential mechanisms of action of multicomponent and multitarget drugs by integrating multidisciplinary information and analyzing the effects via a network of actions at various levels [4]. Therefore, in our study, a network pharmacology approach was utilized to predict the antitumor mechanisms of ZR-CR. A surface plasmon resonance (SPR) assay was used to validate the predicted component-target affinity.

2. Materials and Methods

2.1. Collection of ZR and CR Information

Using TCMSP [5], TCMIP [6], and TCMID [7], and the keywords “Zingiberis rhizoma (Ganjiang in Chinese)” and “Coptidis rhizoma (Huanglian in Chinese),” compounds in ZR and CR were identified. The composition data of the compounds were obtained, and the datasets of the two Chinese medicine compounds were constructed after aggregation and removal of duplicates.

2.2. Collection of ZR and CR Targets

Similarity Ensemble Approach software (SEAware) [8, 9] was used to reverse target the chemical constituents of ZR and CR by setting the screening condition to human and then calculating and predicting the potential targets of ZR and CR. For targets of compounds that had not been collected by the SEAware, the compound name and SMILES code were queried using PubChem [10]; then, the targets of ingredients predicted using STITCH [11] and the Swiss Target Prediction database [12] were used to supplement the findings. All targets obtained were integrated, and the targets of all components were calibrated to the official name based on the Uniprot database [13].

2.3. Collection of Tumor Targets

GeneCards [14], NCBI-gene (https://www.ncbi.nlm.nih.gov/), and DisGeNET [15] databases were searched to identify tumor-related targets using keywords such as “tumor,” “cancer,” and “carcinoma.” The target names were unified using the Uniprot database [13] by integrating the calculated potential targets of ZR and CR with a set relevance . This process enabled us to obtain the potential antitumor targets of the active components of ZR and CR. The R programming language was used to match the targets of the components with those of the disease and construct a Venn diagram.

2.4. Network Construction and Analysis

The target interactions of ZR and CR were introduced into the STRING database [16] to obtain protein–protein interactions, and the active ingredients and targets of ZR and CR were imported into the Cytoscape 3.6.1 software [17] to construct a comprehensive network diagram. The top 20 hub target networks were obtained using the Cytohubba module, and the NetworkAnalyzer module was used to analyze the degree and related parameters of the node. The potential transcription factors (TFs) based on the predicted antitumor targets of ZR and CR were investigated using the Metascape database (TRRUST module) [18, 19].

2.5. Pathway and Process Enrichment Analyses

Pathway and process enrichment analyses of the predicted antitumor targets of ZR and CR were conducted and visualized using the Metascape database [18]. The hub antitumor targets of ZR and CR were introduced into the DAVID database for GO analysis and KEGG pathway analysis [20, 21]. GO analysis was used to annotate and classify targets according to biological process (BP), molecular function (MF), and cellular component (CC). The top 20 terms were imported into the R programming language software for visual display.

2.6. Tumor Types Targeted by ZR and CR

The tumor-related targets of ZR and CR were imported into the Enrichr database [22] to identify the relationship between targets and diseases. The results were visualized using the R language software.

2.7. Molecular Docking

The hub target network was constructed by molecular docking between key targets and their corresponding components. The interaction between small molecules and targets was simulated using professional software to calculate their binding strengths as follows. The crystal structures of four hub targets (AKT1, EGFR, TP53, and VEGFA) were downloaded from the Protein Data Bank database [23] for molecular docking studies. The small molecules corresponding to the key target were generated in the lowest energy conformation. Then, the protein crystals and small molecules were introduced into the Molecular Operating Environment (MOE) software for structural preprocessing, and the dock module in the MOE and AutoDock Vina software was used to complete the molecular docking. The small molecule corresponding to the hub target albumin was oleic acid, which has a highly flexible long chain that might reduce the docking study accuracy, and therefore, we did not conduct a related experiment.

2.8. Materials and Buffer Preparation

Human AKT1, EGFR, and VEGF165 proteins were obtained from ACROBiosystems (Beijing, China), and human TP53 protein was purchased from Active Motif (Shanghai, China). Berberine (CAS: 2086-83-1), quercetin (CAS: 117-39-5), and daucosterol (CAS: 474-58-8) were purchased from Solarbio (Beijing, China), and ginkgetin (CAS: 481-46-9) and isoginkgetin (CAS: 548-19-6) were purchased from MedChemExpress (Shanghai, China). KH₂PO₄, NaCl, Na₂HPO₄, and KCl were obtained from Sinopharm Chemical Reagent (Shanghai, China), and dimethyl sulfoxide (DMSO) was purchased from Sangon Biotech (Shanghai, China).

The Series S Sensor Chip CM5 and Amine Coupling kit were obtained from GE Healthcare (Little Chalfont, UK). The prepared buffers were the running buffer (pH 7.4): (2 mM KH₂PO₄, 10 mM Na₂HPO₄, 137 mM NaCl, and 2.7 mM KCl) with 0.05% Tween 20 and 5% DMSO and immobilization buffer (pH 4.0): 10 mM sodium acetate.

2.9. Affinity Test

The SPR assay was conducted using a Biacore 8K instrument (GE Healthcare) with the running buffer. Briefly, 50 μg/mL of the target proteins in 10 mM NaAc (pH 4.5; Sinopharm Chemical Reagent, Shanghai, China) was then injected into the sample channel of the CM5 chips at a flow rate of 10 μL/min to reach an immobilization level. The detection was performed according to the protocol provided by GE Healthcare. Gradient concentrations of the components (0, 1, 10, and 100 μM) were dissolved in the running buffer and then injected into the channel at a flow rate of 30 μL/min, with an association phase of 60 s, followed by a 90-second dissociation. The affinity was analyzed using the Biacore 8K software.

2.10. Hub Gene Expression and Survival Analysis in Patients

The relative gene expression level of 20 hub antitumor targets of ZR and CR was analyzed using the GEPIA database [24]. The differential expression between tumor and adjacent normal tissues for hub genes across all The Cancer Genome Atlas (TCGA) tumors was analyzed using the TIMER database [25]. The GEPIA database was also used to conduct overall survival and disease-free survival analyses based on gene expression. Tumor names as in TCGA were adopted in the analysis process (Table S1).

3. Results

3.1. Potential Tumor Targets of ZR and CR

We identified 157 and 55 chemical components of ZR and CR, respectively, for a total of 212 ingredients (Table S2). Then, 416 corresponding potential targets were selected, consisting of 249 and 288 targets in ZR and CR, respectively, which included 121 shared targets between ZR and CR. In addition, a total of 31490 tumor-related targets were identified. After intersecting the component and tumor targets, 239 common targets were obtained, including 79 targets that were shared by ZR and CR. The 239 common targets were considered potential targets for ZR and CR against tumors (Figure 1).

3.2. Protein–Protein Interaction and Component-Target Networks

The protein–protein interaction network of potential targets of ZR and CR against tumors contained 239 nodes and 2814 edges (Figure 2). The hub target network shows the top 20 targets of ZR and CR against tumors (Figure 3 and Table 1). Nodes changed from orange to red and small to large, indicating that the degree gradually increased, and the importance of the targets was determined by the degree value. Upon analysis of 239 predicted antitumor targets of ZR and CR, 20 key potential TFs were identified (Figure 4).

The network of the components in ZR and CR corresponding to tumor targets contained 380 nodes and 1244 edges (Figure 5), and the network of hub components and targets contained 72 nodes and 117 edges (Figure 6), where the nodes represent the components of ZR and CR or the targets corresponding to the components, whereas the edges indicate interactions between components and targets.

3.3. Pathway and Process Enrichment

The top 20 terms in the pathway and process enrichment analyses of the predicted antitumor targets of ZR and CR via the Metascape database are shown in Figure 7. The GO enrichment analysis indicated that the hub antitumor targets of ZR and CR are involved in 144, 18, and 42 BP-, CC-, and MF-related entries, respectively. The top-ranking entries based on value are shown in Figure 8; the larger the dots, the more the number of enriched targets and the redder the dots, the smaller the value. The top 20 KEGG enrichment entries of the hub antitumor targets of ZR and CR are shown in Figure 9. The top 20 disease types of Enrichr enrichment entries of the hub targets are shown in Figure 10, which contains multiple tumor types.

(a)

(b)

(c)

3.4. Molecular Docking between Components and Target Proteins

The MOE molecular docking scoring value was used to evaluate the interactions between the small molecules and proteins (Table 2); the smaller the scoring value, the stronger the interaction. The results showed that berberine and AKT1, berberine and TP53, and quercetin and EGFR may have stronger binding effects than observed with other interactions. The binding free energy of the first molecule docking posture of these three interactions is shown in Figure 11. Because numerous components corresponded to the VEGFA protein, AutoDock Vina was used for molecular docking. The top 10 terms identified by AutoDock Vina are shown in Table 3; the smaller the binding energy, the stronger the binding capacity of the target protein and component.

(a)

(b)

(c)

3.5. Affinity of Components for Target Proteins

The relative response values of the SPR affinity tests suggested that several components of ZR and CR can directly bind to its tumor target protein (Table 4 and Figure 12), and the affinity between ginkgetin, isoginkgetin, and daucosterol and their common target VEGF165 was weak.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

3.6. Clinical Significance of Hub Targets

The relative gene expression levels of 20 hub antitumor targets of ZR and CR in different tumor tissues were compared. The profiles showed that different hub targets are differentially expressed in multiple tumor tissues (Figure 13), and AKT1, TP53, VEGFA, and EGFR are differentially expressed between multiple tumors and adjacent normal tissues (Figure 14), such as breast invasive carcinoma, cholangiocarcinoma, prostate adenocarcinoma, and thyroid carcinoma. The results also showed that hub gene expression levels are associated with unfavorable overall survival and disease-free survival in patients with multiple tumor types (Figures 15 and 16), such as patients with kidney renal clear cell carcinoma, pancreatic adenocarcinoma, and brain lower grade glioma.

4. Discussion

Presently, the ZR-CR herbal pair is widely used in the field of antitumor TCM. In this study, the mechanisms of antitumor action of the ZR-CR herbal pair were investigated using a network pharmacology approach. The network pharmacology analysis suggested that the active components of ZR and CR, such as gingerdione, shogaol, zingerone, sitosterol, daucosterol, ginkgetin, isoginkgetin, berberine, quercetin, noroxyhydrastinine, chlorogenic acid, vanillic acid, and ethyl caffeate, identified using the TCMSP, TCMIP, and TCMID platforms, could serve as effective therapeutic agents for the treatment of tumors via multiple mechanisms.

The component-target network showed that, although multiple components corresponded to a single target, a single component also corresponded to multiple targets, indicating that the effects of the ZR-CR herbal pair on tumors are mediated through multiple components and multiple targets. Among the numerous targets, hub targets such as AKT1, TP53, VEGFA, EGFR, and MAPK1 may play important roles. Furthermore, hub targets such as TP53, JUN, ESR1, EP300, PPARG, and AR are also TFs that may be involved in the regulation of their target genes.

Pathway and process enrichment analyses suggested that the ZR-CR herbal pair is involved in mechanisms such as the metabolism of lipids, response to oxidative stress, and regulation of hormone levels, cell migration, MAPK cascade, and programmed cell death to exert antitumor effects. The enrichment analyses of hub targets suggested that ZR and CR are involved in estrogen, MAPK, PI3K-AKT, TNF, FOXO, HIF-1, and VEGF signaling pathways, which are also closely related to antitumor activities. For example, abnormalities in MAPK signaling play a crucial part in the progression of cancer [26]. FOXO is a subfamily of the forkhead transcription factor family that is involved in cell fate decisions [27]. HIF-1 functions as a signaling hub that coordinates the activities of many TFs and signaling molecules that affect tumorigenesis [28].

Disease type enrichment analysis showed that the ZR-CR herbal pair plays a role in tumors of the digestive, urinary, and gynecological systems. This finding is also consistent with the extensive application of ZR and CR in the treatment of digestive system diseases and tumors in TCM [29, 30].

Molecular docking analysis and SPR assays were conducted based on the hub targets of ZR and CR against tumors to preliminarily screen and validate our predictions. The results suggested that berberine directly binds to AKT1 and TP53 and quercetin to EGFR and VEGF165. Ginkgetin, isoginkgetin, and daucosterol directly bind to VEGF165 but with weak affinities. Although affinity is a crucial functional factor of pharmacological activity, its strength does not necessarily determine the potency of the effect. For example, one study investigated the effect of affinity and antigen internalization on the uptake and penetration of anti-HER2 antibodies in solid tumors. The findings suggested that antibodies with lower affinity penetrate tumors more effectively when rates of antibody-antigen dissociation are higher than those of antigen internalization [31]. Some studies have shown that these identified compounds also indirectly affect target proteins and exhibit antitumor activities. For example, berberine induces p53 expression and thereby decreases the mitochondrial membrane potential and induces cytochrome C release [32]. The AKT-related mitochondrial pathway may be partly involved in berberine-induced apoptosis of gastric cancer cells [33]. Quercetin prevents prostate cancer progression in an in vivo model by inhibiting EGFR signaling [34]. In addition, the combination of quercetin and metformin exerts synergistic antitumor effects by inhibiting the VEGF/PI3K/AKT pathway [35]. Ginkgetin suppresses VEGF-mediated angiogenesis during cancer development [36]. Furthermore, isoginkgetin, an isomer of ginkgetin, markedly decreases MMP9 expression and invasion by inhibiting PI3K/AKT [37]. Daucosterol linoleate suppresses VEGF, MMP2, and MMP9 expression in breast cancer [38].

These previous studies mostly focused on evaluating pharmacological parameters, functional phenotype changes, and molecular mechanisms limited to some classical signaling pathways. Consequently, direct interactions between active molecules and target proteins are often not considered [39]. The current results show the direct binding of small-molecule compounds to tumor-related proteins, which confirmed that ZR and CR components act on multiple tumors.

Furthermore, hub targets, such as AKT1, TP53, VEGFA, and EGFR, were shown to be differentially expressed between multiple tumors and adjacent normal tissues, and their expression levels can impact overall survival and disease-free survival in patients with tumors. These pan-cancer analyses are expected to benefit the clinical application of the ZR-CR herbal pair. The hub targets identified in our study are considered targets of currently available antitumor drugs such as EGFR inhibitors (gefitinib, erlotinib, and lapatinib) [40]. Several anti-VEGF drugs have also been approved for certain advanced cancers, but they have only exhibited limited benefits on the overall survival of patients with tumors and have rarely resulted in durable responses [41]. TCM components often exhibit less potent effects than modern drugs with clear targets, but the overall efficacy advantage of TCM preparations for tumor treatment is not only based on a simple additive effect of multiple “weak effects.” In vitro and in vivo EGFR mutant non-small-cell lung cancer models have shown that 20% of individual effective drug doses can completely block MAPK signaling when used in RAF+MEK+ERK or EGFR+RAF+MEK+ERK inhibitor combinations [42]. Our findings suggest that the synergistic antitumor effect of the compounds in ZR and CR deserves more detailed exploration.

Although some components of ZR-CR with activity against tumors in the hub network were not analyzed further, they are not unimportant. Related studies have been reported, such as one that showed that shogaol causes cancer cell death by inducing cell cycle arrest and apoptosis [43]. Additionally, the antitumor activity of chlorogenic acid is mediated through the induction of cancer cell differentiation [44], and vanillic acid inhibits hypoxia-induced HIF-1α expression in various human cancer cell lines [45]. Therefore, these components of the ZR-CR herbal pair also warrant further experimental and clinical studies.

5. Conclusions

The active ingredients of the ZR-CR herbal pair and their underlying antitumor mechanisms were partly identified. The network pharmacology approach revealed that the components of ZR-CR, such as shogaol, daucosterol, ginkgetin, berberine, quercetin, chlorogenic acid, and vanillic acid, have potential efficacy as antitumor agents; they mediate their activity via multiple mechanisms, such as MAPK, PI3K-AKT, TNF, FOXO, HIF-1, and VEGF signaling pathways. The ZR-CR herbal pair could be useful in the treatment of tumors, such as those of the digestive, urinary, and gynecological systems. This network-based investigation could facilitate the elucidation of the mechanisms of action of the ZR-CR herbal pair against tumors. The molecular docking approach and SPR assay identified some components of ZR-CR that directly bind to hub target proteins to exert their antitumor effects. The gene expression levels of hub targets of the ZR-CR herbal pair are associated with overall survival and disease-free survival in patients with multiple tumor types. Finally, the antitumor components of the ZR-CR warrant further experimental and clinical investigations, as they have not been discussed adequately in our study.

Data Availability

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

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Meng Wang and Youke Qi contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant number 81874377) and the internal fund of the Shanghai Jiao Tong University Affiliated Sixth People’s Hospital (grant number LYZY-0297).

Supplementary Materials

Table S1: tumor name in TCGA. Table S2: components of Zingiberis rhizoma and Coptidis rhizome. (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
S. Gupta, A. Harper, Y. Ruan et al., “International trends in the incidence of cancer among adolescents and young adults,” JNCI: Journal of the National Cancer Institut, vol. 112, no. 11, pp. 1105–1117, 2020.
View at: Publisher Site | Google Scholar
X. Liu, L. J. Xiu, J. P. Jiao et al., “Traditional Chinese medicine integrated with chemotherapy for stage IV non-surgical gastric cancer: a retrospective clinical analysis,” Journal of integrative medicine, vol. 15, no. 6, pp. 469–475, 2017.
View at: Publisher Site | Google Scholar
A. L. Hopkins, “Network pharmacology: the next paradigm in drug discovery,” Nature Chemical Biology, vol. 4, no. 11, pp. 682–690, 2008.
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
H. Y. Xu, Y. Q. Zhang, Z. M. Liu et al., “ETCM: an encyclopaedia of traditional Chinese medicine,” Nucleic Acids Research, vol. 47, no. D1, pp. D976–D982, 2019.
View at: Publisher Site | Google Scholar
L. Huang, D. Xie, Y. Yu et al., “TCMID 2.0: a comprehensive resource for TCM,” Nucleic Acids Research, vol. 46, no. D1, pp. D1117–D1120, 2018.
View at: Publisher Site | Google Scholar
M. J. Keiser, B. L. Roth, B. N. Armbruster, P. Ernsberger, J. J. Irwin, and B. K. Shoichet, “Relating protein pharmacology by ligand chemistry,” Nature Biotechnology, vol. 25, no. 2, pp. 197–206, 2007.
View at: Publisher Site | Google Scholar
M. J. Keiser, V. Setola, J. J. Irwin et al., “Predicting new molecular targets for known drugs,” Nature, vol. 462, no. 7270, pp. 175–181, 2009.
View at: Publisher Site | Google Scholar
S. Kim, J. Chen, T. Cheng et al., “PubChem 2019 update: improved access to chemical data,” Nucleic Acids Research, vol. 47, no. D1, pp. D1102–D1109, 2019.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Santos, C. von Mering, L. J. Jensen, P. Bork, and M. Kuhn, “STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data,” Nucleic Acids Research, vol. 44, no. D1, pp. D380–D384, 2016.
View at: Publisher Site | Google Scholar
A. Daina, O. Michielin, and V. Zoete, “SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules,” Nucleic Acids Research, vol. 47, no. W1, pp. W357–W364, 2019.
View at: Publisher Site | Google Scholar
UniProt Consortium, “UniProt: a worldwide hub of protein knowledge,” Nucleic Acids Research, vol. 47, no. D1, pp. D506–D515, 2019.
View at: Publisher Site | Google Scholar
G. Stelzer, N. Rosen, I. Plaschkes et al., “The GeneCards suite: from gene data mining to disease genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, pp. 1.30.1–1.30.33, 2016.
View at: Google Scholar
J. Piñero, J. M. Ramírez-Anguita, J. Saüch-Pitarch et al., “The DisGeNET knowledge platform for disease genomics: 2019 update,” Nucleic Acids Research, vol. 48, no. D1, pp. D845–D855, 2020.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, p. 1523, 2019.
View at: Publisher Site | Google Scholar
H. Han, J. W. Cho, S. Lee et al., “TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions,” Nucleic Acids Research, vol. 46, no. D1, pp. D380–D386, 2018.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, pp. 44–57, 2009.
View at: Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists,” Nucleic Acids Research, vol. 37, pp. 1–13, 2009.
View at: Google Scholar
M. V. Kuleshov, M. R. Jones, A. D. Rouillard et al., “Enrichr: a comprehensive gene set enrichment analysis web server 2016 update,” Nucleic Acids Research, vol. 44, no. W1, pp. W90–W97, 2016.
View at: Publisher Site | Google Scholar
S. K. Burley, H. M. Berman, C. Bhikadiya et al., “RCSB Protein Data Bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy,” Nucleic Acids Research, vol. 47, no. D1, pp. D464–D474, 2019.
View at: Publisher Site | Google Scholar
Z. Tang, B. Kang, C. Li, T. Chen, and Z. Zhang, “GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis,” Nucleic Acids Research, vol. 47, no. W1, pp. W556–W560, 2019.
View at: Publisher Site | Google Scholar
T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
View at: Publisher Site | Google Scholar
A. S. Dhillon, S. Hagan, O. Rath, and W. Kolch, “MAP kinase signalling pathways in cancer,” Oncogene, vol. 26, no. 22, pp. 3279–3290, 2007.
View at: Publisher Site | Google Scholar
M. Farhan, H. Wang, U. Gaur, P. J. Little, J. Xu, and W. Zheng, “FOXO signaling pathways as therapeutic targets in cancer,” International Journal of Biological Sciences, vol. 13, no. 7, pp. 815–827, 2017.
View at: Publisher Site | Google Scholar
K. Balamurugan, “HIF-1 at the crossroads of hypoxia, inflammation, and cancer,” International Journal of Cancer, vol. 138, no. 5, pp. 1058–1066, 2016.
View at: Publisher Site | Google Scholar
S. Chrubasik, M. H. Pittler, and B. D. Roufogalis, “Zingiberis rhizoma: a comprehensive review on the ginger effect and efficacy profiles,” Phytomedicine, vol. 12, no. 9, pp. 684–701, 2005.
View at: Publisher Site | 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
S. I. Rudnick, J. Lou, C. C. Shaller et al., “Influence of affinity and antigen internalization on the uptake and penetration of anti-HER2 antibodies in solid tumors,” Cancer Research, vol. 71, no. 6, pp. 2250–2259, 2011.
View at: Publisher Site | Google Scholar
J. P. Lin, J. S. Yang, J. H. Lee, W. T. Hsieh, and J. G. Chung, “Berberine induces cell cycle arrest and apoptosis in human gastric carcinoma SNU-5 cell line,” World Journal of Gastroenterology, vol. 12, no. 1, pp. 21–28, 2006.
View at: Publisher Site | Google Scholar
T. Yi, L. Zhuang, G. Song, B. Zhang, G. Li, and T. Hu, “Akt signaling is associated with the berberine-induced apoptosis of human gastric cancer cells,” Nutrition and Cancer, vol. 67, no. 3, pp. 523–531, 2015.
View at: Publisher Site | Google Scholar
A. B. Firdous, G. Sharmila, S. Balakrishnan et al., “Quercetin, a natural dietary flavonoid, acts as a chemopreventive agent against prostate cancer in an in vivo model by inhibiting the EGFR signaling pathway,” Food & Function, vol. 5, no. 10, pp. 2632–2645, 2014.
View at: Publisher Site | Google Scholar
S. Sun, F. Gong, P. Liu, and Q. Miao, “Metformin combined with quercetin synergistically repressed prostate cancer cells via inhibition of VEGF/PI3K/Akt signaling pathway,” Gene, vol. 664, pp. 50–57, 2018.
View at: Publisher Site | Google Scholar
W. H. Hu, G. K. Chan, R. Duan et al., “Synergy of ginkgetin and resveratrol in suppressing VEGF-induced angiogenesis: a therapy in treating colorectal cancer,” Cancers, vol. 11, article 1828, 2019.
View at: Google Scholar
S. O. Yoon, S. Shin, H. J. Lee, H. K. Chun, and A. S. Chung, “Isoginkgetin inhibits tumor cell invasion by regulating phosphatidylinositol 3-kinase/Akt-dependent matrix metalloproteinase-9 expression,” Molecular Cancer Therapeutics, vol. 5, no. 11, pp. 2666–2675, 2006.
View at: Publisher Site | Google Scholar
B. Han, P. Jiang, W. Liu et al., “Role of daucosterol linoleate on breast cancer: studies on apoptosis and metastasis,” Journal of Agricultural and Food Chemistry, vol. 66, no. 24, pp. 6031–6041, 2018.
View at: Publisher Site | Google Scholar
P. F. Tu, K. W. Zeng, L. X. Liao, and X. M. Song, “Progress of methodology for identifying target protein of natural active small molecules,” Zhongguo Zhong Yao Za Zhi, vol. 41, no. 1, pp. 6–13, 2016.
View at: Publisher Site | Google Scholar
M. Ammad-Ud-Din, S. A. Khan, K. Wennerberg, and T. Aittokallio, “Systematic identification of feature combinations for predicting drug response with Bayesian multi-view multi-task linear regression,” Bioinformatics, vol. 33, no. 14, pp. i359–i368, 2017.
View at: Publisher Site | Google Scholar
R. R. Ramjiawan, A. W. Griffioen, and D. G. Duda, “Anti-angiogenesis for cancer revisited: is there a role for combinations with immunotherapy?” Angiogenesis, vol. 20, no. 2, pp. 185–204, 2017.
View at: Publisher Site | Google Scholar
J. M. Neto, E. Nadal, E. Bosdriesz et al., “Multiple low dose therapy as an effective strategy to treat EGFR inhibitor-resistant NSCLC tumours,” Nature Communications, vol. 11, no. 1, p. 3157, 2020.
View at: Publisher Site | Google Scholar
L. W. Qi, Z. Zhang, C. F. Zhang et al., “Anti-colon cancer effects of 6-shogaol through G2/M cell cycle arrest by p53/p21-cdc2/cdc25A crosstalk,” The American Journal of Chinese Medicine, vol. 43, no. 4, pp. 743–756, 2015.
View at: Publisher Site | Google Scholar
S. Huang, L. L. Wang, N. N. Xue et al., “Chlorogenic acid effectively treats cancers through induction of cancer cell differentiation,” Theranostics, vol. 9, no. 23, pp. 6745–6763, 2019.
View at: Publisher Site | Google Scholar
J. Gong, S. Zhou, and S. Yang, “Vanillic acid suppresses HIF-1α expression via inhibition of mTOR/p70s6K/4E-BP1 and Raf/MEK/ERK pathways in human colon cancer HCT116 cells,” International Journal of Molecular Sciences, vol. 20, no. 3, p. 465, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Meng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

597

Downloads

1342

Citations