Prediction of triptolide targets in colorectal cancer using network pharmacology and molecular docking

doi:10.21203/rs.3.rs-56695/v1

Download PDF

Research

Prediction of triptolide targets in colorectal cancer using network pharmacology and molecular docking

https://doi.org/10.21203/rs.3.rs-56695/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Network pharmacology is an approach that uses bioinformatics to predict and identify multiple drug targets and interactions in disease. Here, we applied network pharmacology to investigate the potential mechanisms of action of triptolide, an active component in the traditional Chinese medicine Tripterygium wilfordii Hook F, in colorectal cancer (CRC). We first searched public databases for genes and proteins known to be associated with CRC, as well as those predicted to be targets of triptolide, and then used Ingenuity Pathway Analysis (IPA) to identify enriched gene pathways and networks. Networks and pathways that overlapped between CRC-associated proteins and triptolide target proteins were then used to predict candidate protein targets of triptolide in CRC. The following proteins were found to be expressed in both CRC-associated networks and triptolide target networks: JUN, FOS, CASP3, BCL2, IFNG, and VEGFA. Docking studies suggested that triptolide can fit in the binding pocket of the four top candidate triptolide target proteins (CASP3, BCL2, VEGFA and IFNG). The overlapping pathways were activation of neuroinflammation signaling, glucocorticoid receptor signaling, T helper (Th) cell differentiation, Th1/Th2 activation, and colorectal cancer metastasis signaling. These results show that network pharmacology can be used to generate hypotheses about how triptolide exerts therapeutic effects in CRC. Network pharmacology may be a useful method for characterizing multi-target drugs in complex diseases.

Bioinformatics

Triptolide

Colorectal cancer

Ingenuity Pathway Analysis

Network pharmacology

Molecular docking

Colorectal cancer (CRC) is the second leading cause of cancer-related death in the world [1]. Many Asian countries, including China, South Korea, Singapore, and Japan, are experiencing an escalating incidence of CRC [2-9]. A lot of effort has been conducted to enhance the treatment of CRC. Despite improved surgical techniques and advancements in radio- and chemotherapy over the past few decades, the overall survival rate of patients with CRC has not improved substantially [10]. It is therefore imperative to devise novel strategies for safe and effective treatment of CRC.

Triptolide is a diterpene triepoxide purified from Tripterygium wilfordii Hook F, commonly known as 'lei gong teng' or 'thunder god vine', a medicinal plant whose extracts have been used in traditional Chinese medicine for treating rheumatoid arthritis and other inflammatory diseases [11-16]. Recent studies have shown that triptolide kills CRC and other cancer cells in vitro with high potency [17, 18]. Animal studies have shown that triptolide inhibits the growth of CRC cells in a mouse xenograft model [19, 20]. However, the mechanisms underlying the therapeutic effects of triptolide in CRC are unclear.

Network-based drug discovery is a promising, cost-effective drug development approach based on bioinformatics, systems biology and pharmacology. Instead of the current “one target, one drug” approach, network pharmacology utilizes a “network target, multicomponent” strategy. Because network pharmacology can provide a good understanding of the principles of network theory and systems biology, it has been considered the next paradigm in drug discovery [21-24]. Herbal medicines such as triptolide are particularly promising drugs for developing new multicomponent and multitarget synergistic cancer therapeutics [25-27].

In this study, we investigated the potential mechanisms by which triptolide acts in CRC using network pharmacology. Our results revealed potential mechanisms that may underlie the therapeutic effects of triptolide in CRC and showed that network pharmacology is a useful tool to facilitate the development of novel CRC drugs.

2.1. Identification of genes and proteins associated with CRC and potential triptolide target proteins

We searched the National Center for Biotechnology Information (NCBI) GenBank Database (https://www.ncbi.nlm.nih.gov/) for genes related to the search term “colorectal cancer”. Search hits were filtered to retain only studies performed in Homo sapiens. Potential human protein targets of triptolide were obtained from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (http://lsp.nwu.edu.cn/tcmsp.php) [32] and PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [33].

2.2. Prediction of pathways and networks associated with CRC and triptolide targets

The protein symbols of candidate triptolide targets and CRC-associated genes were uploaded into Ingenuity Pathway Analysis (IPA) software version 2019 (https://www.qiagenbioinformatics.com, Redwood City, CA, US) for network and pathway analysis. IPA was used to construct pathways and networks based on known interactions between genes and proteins. Enrichment analysis of CRC target Gene Ontology (GO) enrichment and network was performed using R (version 3.6.0 for Windows) and Cytoscape version 3.6.1 (http://www.cytoscape.org). Pathways and networks were ranked according to the number of molecules in pathways and networks, and cut-off p value < 0.05 was used to identify significantly enriched pathways/networks. Pathways and networks involving CRC-associated genes and candidate triptolide targets were identified using the “Compare” module in IPA.

2.3. Construction of common networks and core target screening

Proteins that were both previously associated with CRC and predicted to be targets of triptolide were collected using Venny 2.1.0 software (http://bioinfogp.cnb.csic.es/tools/venny/), and potential protein-protein interactions were analyzed in STRING 10.5 (https://string-db.org/) after filtering by Homo sapiens. Target proteins whose values for the topological attributes of node degree distribution and betweenness centrality were above the means were defined as "core targets". The degree of a node is the number of nodes to which it is linked, while betweenness reflects the extent to which nodes lie between one another. Finally, core targets were analyzed using IPA.

2.4. Prediction of binding between triptolide and candidate target proteins

The crystal structures of candidate proteins bound by triptolide were downloaded from the RCSB Protein Data Bank (http://www.pdb.org/) and modified using YASARA software (http://www.modekeji.cn) to remove ligands, add hydrogen, remove water, and optimize and patch amino acids. Before docking, ChemBioDraw 3D was used to generate 3D chemical structures and minimize binding energies for all candidate triptolide targets. We used YASARA software to test the precision of docking between triptolide and candidate target proteins because it had the highest accuracy and consistency (data not shown). The best docking poses were identified as those showing the smallest root mean square deviation (RMSD) between the predicted conformation and the observed X-ray crystallographic conformation. Models with an RMSD ≤ 4 Å were considered reliable and those with an RMSD ≤ 2 Å were considered accurate [34]. The network pharmacology approach in this study is summarized in Fig. 1.

3.1. Networks and enriched functions in CRC-associated genes

A total of 3,298 CRC-associated genes were identified in GenBank. IPA identified a total of 596 pathways and 25 networks associated with these genes. The top/majority of pathways were involved in molecular mechanisms of cancer, CRC metastasis signaling and Wnt/β-catenin signaling. The majority of networks were involved in cancer, cellular movement, organismal injury and abnormalities, cellular development, embryonic development, organismal development, cell-to-cell signaling and interaction, protein synthesis, and RNA damage and repair (Fig. 2).

Gene Ontology (GO) enrichment and network analysis of CRC-associated proteins showed that the top three functions were epithelial cell proliferation, ameboidal-type cell migration and regulation of vasculature development (Fig. 3).

3.2. Networks and enriched functions in triptolide target genes and proteins

A total of 33 proteins were identified as candidate triptolide targets, and IPA revealed a total of 294 enriched pathways and 10 networks. The most significantly enriched pathways were neuroinflammation signaling, glucocorticoid receptor signaling, T helper (Th) cell differentiation, Th1 and Th2 activation, and CRC metastasis signaling. The top networks identified were involved in gene expression, cellular function and maintenance, cell cycle, inflammatory response, organismal injury and abnormalities, cell-to-cell signaling and interaction, cell death and survival, dermatological diseases and conditions, and infectious diseases (Fig. 4).

3.3. Networks of shared proteins and special proteins targeted by triptolide

The intersections between the set of potential triptolide targets and CRC-related proteins were analyzed using Venny software, which identified 29 shared proteins (Fig. 5A). STRING suggested that the 29 proteins can interact with one another via 269 interactions (edges) (Fig. 5B).

Proteins linked to CRC and potentially targeted by triptolide participate in several canonical pathways involved in a range of biological activities. To demonstrate the ability of our integrative bioinformatic approach to propose specific protein targets for further mechanistic studies, we selected the top pathway in the IPA categories “molecular mechanisms of cancer” and “colorectal cancer metastasis signaling” that were linked to CRC and targeted by triptolide. Several nodes in this pathway emerged as potential direct targets of triptolide in CRC: JUN, FOS, CASP3, BCL2, IFNG, and VEGFA (Fig. 5C and D). Combining these results with STRING analysis, we identified JUN, FOS, CASP3, BCL2, IFNG, and VEGFA as particularly likely targets of triptolide in CRC.

3.4. Predicted binding of triptolide to target proteins in CRC

To further validate candidate triptolide targets in CRC, we tested the precision of docking between triptolide and the following potential target proteins by YASARA software (Fig. 6): (A) CASP3 (PDB: 5IAN), (B) BCL2 (PDB: 6GL8), (C) VEGFA (PDB: 4KZN) and (D) IFNG (PDB: 1EKU). As shown in Fig. 6, triptolide binds to the active sites of these target proteins and interacts with several amino acid residues, with most interactions being hydrophobic.

For instance, in the combination of triptolide with VEGFA, there are different hydrophobic interactions between triptolide and residues of VEGFA such as Gln-22, Tyr-25, His-27 and Pro-28. In addition, triptolide can form a hydrogen bond with a length of 1.8 Å and a bond energy of 8.40 kJ/mol with the Arg-207 residue of CASP3, and its aromatic ring forms a π-π interaction with the aromatic ring in Phe-256 of CASP3. Overall, these results provide further evidence that these four proteins may act as triptolide targets in CRC.

The network pharmacology approach is relatively new and was first proposed by Li et al. in 2014 [28]. Because it provides a more complete understanding of network theory and systems biology, it has been considered the next paradigm in drug discovery [29, 30]. Network pharmacology has been used to study pathways of interaction between drugs and proteins or genes and diseases, and it is capable of describing complexities among biological systems, drugs and diseases from a network perspective [31-34]. Therefore, the development of network pharmacology techniques that can predict multiple drug-target interactions may hold the key to future drug discoveries in complex diseases such as CRC. In this work, we integrated information from publicly available databases to predict interactions between triptolide and its potential targets in CRC, as well as the signaling pathways and networks involved. Pathway analysis suggested that triptolide regulates the activation of neuroinflammation signaling, glucocorticoid receptor signaling, Th cell differentiation, Th1 and Th2 activation, and metastasis signaling in CRC.

These results are consistent with several in vitro and in vivo studies. Previous work [35] suggested that triptolide is able to induce G1 cell cycle arrest by inhibiting transcriptional activation of E2F1. Triptolide was found to reduce both tumor number and tumor size in mice carrying mutations to promote growth of adenomatous polyposis coli as well as in mice treated with azoxymethane/dextran sodium sulfate to induce cancer[36]. Triptolide effectively inhibits CRC cell proliferation, colony formation, and organoid growth in vitro, and these effects are associated with down-regulation of target genes transcribed by RNA polymerase III. Moreover, triptolide inhibits cyclooxygenase-2 and inducible nitric oxide synthase expression in human colon cancer [37].

Our study presents several limitations. Although we predicted and verified potential triptolide target proteins by molecular docking, our results are based on bioinformatic predictions and require experimental confirmation. As next steps, we intend to use triptolide to treat CRC cells, identify its targets using transcriptomics and proteomics, and thereby elucidate its mechanisms of action.

In summary, our results predict that the therapeutic effects of triptolide in CRC are mediated, at least in part, via JUN, FOS, CASP3, BCL2, IFNG, and VEGFA. These results may be useful in guiding further research to determine the molecular targets of triptolide in CRC, and they may inspire further application of network pharmacology to drug discovery.

Funding

The present study was supported by The National Natural Science Foundation of China (grant no. U1804179), The Key Scientific and Technological Projects in Henan Province (grant no. 202102310190), The Henan Science and Technology Innovation Team, Investigation on Plant Resources in Dabie Mountains and The Study and Utilization of Active Components of Special Plants (grant no. 2017083), and The Nanhu Scholars Program for Young Scholars of Xinyang Normal University (grant no. 2018001).

Authors' contributions

Xinqiang SONG conceived and designed the study, and wrote the manuscript.

Yu Zhang, Erqin Dai, Qiyue Zhang, Nongyi Zheng, Lei Wang, Hongtao Du analyzed the data.

All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Weitz, J., et al., Colorectal cancer. Lancet, 2005. 365(9454): p. 153-65.
Brody, H., Colorectal cancer. Nature, 2015. 521(7551): p. S1.
Dekker, E., et al., Colorectal cancer. Lancet, 2019. 394(10207): p. 1467-1480.
Kwon, W.A. and H.K. Seo, Usefulness of the mass screening program for colorectal cancer in China: further long-term validation is needed to confirm its value. Ann Transl Med, 2020. 8(7): p. 427.
Luu, H.N., et al., Association Between Leukocyte Telomere Length and Colorectal Cancer Risk in the Singapore Chinese Health Study. Clin Transl Gastroenterol, 2019. 10(5): p. 1-9.
Kobayashi, K., et al., Predictors of invasive cancer of large laterally spreading colorectal tumors: A multicenter study in Japan. JGH Open, 2020. 4(1): p. 83-89.
Okada, T., et al., Colorectal cancer risk factors in asymptomatic Chilean population: a survey of international collaboration between Japan and Chile. Eur J Cancer Prev, 2020. 29(2): p. 127-133.
Shin, D.W., et al., Disparities in the Participation Rate of Colorectal Cancer Screening by Fecal Occult Blood Test among People with Disabilities: A National Database Study in South Korea. Cancer Res Treat, 2020. 52(1): p. 60-73.
Yamaji, T., et al., Fat mass and obesity-associated gene polymorphisms, pre-diagnostic plasma adipokine levels and the risk of colorectal cancer: The Japan Public Health Center-based Prospective Study. PLoS One, 2020. 15(2): p. e0229005.
Sao Juliao, G.P., et al., New Strategies in Rectal Cancer. Surg Clin North Am, 2017. 97(3): p. 587-604.
Xi, C., et al., Toxicity of triptolide and the molecular mechanisms involved. Biomed Pharmacother, 2017. 90: p. 531-541.
Hou, W., B. Liu, and H. Xu, Triptolide: Medicinal chemistry, chemical biology and clinical progress. Eur J Med Chem, 2019. 176: p. 378-392.
Fan, D., et al., The Effect of Triptolide in Rheumatoid Arthritis: From Basic Research towards Clinical Translation. Int J Mol Sci, 2018. 19(2).
Fan, D., et al., Triptolide Modulates TREM-1 Signal Pathway to Inhibit the Inflammatory Response in Rheumatoid Arthritis. Int J Mol Sci, 2016. 17(4): p. 498.
Song, X., et al., Prediction of triptolide targets in rheumatoid arthritis using network pharmacology and molecular docking. Int Immunopharmacol, 2020. 80: p. 106179.
Zhang, L., et al., Novel Carboxylated Chitosan-Based Triptolide Conjugate for the Treatment of Rheumatoid Arthritis. Pharmaceutics, 2020. 12(3).
Noel, P., et al., Triptolide and Its Derivatives as Cancer Therapies. Trends Pharmacol Sci, 2019. 40(5): p. 327-341.
Song, W., et al., Preclinical Pharmacokinetics of Triptolide: A Potential Antitumor Drug. Curr Drug Metab, 2019. 20(2): p. 147-154.
Wei, Y.M., et al., Triptolide, A Potential Autophagy Modulator. Chin J Integr Med, 2019. 25(3): p. 233-240.
Xu, H. and B. Liu, Triptolide-targeted delivery methods. Eur J Med Chem, 2019. 164: p. 342-351.
Danhof, M., Systems pharmacology - Towards the modeling of network interactions. Eur J Pharm Sci, 2016. 94: p. 4-14.
Ye, H., et al., Drug Repositioning Through Network Pharmacology. Curr Top Med Chem, 2016. 16(30): p. 3646-3656.
Boezio, B., et al., Network-based Approaches in Pharmacology. Mol Inform, 2017. 36(10).
Zhang, W., et al., Systems Pharmacology for Investigation of the Mechanisms of Action of Traditional Chinese Medicine in Drug Discovery. Front Pharmacol, 2019. 10: p. 743.
Engin, H.B., et al., Network-based strategies can help mono- and poly-pharmacology drug discovery: a systems biology view. Curr Pharm Des, 2014. 20(8): p. 1201-7.
Hao da, C. and P.G. Xiao, Network pharmacology: a Rosetta Stone for traditional Chinese medicine. Drug Dev Res, 2014. 75(5): p. 299-312.
Kibble, M., et al., Network pharmacology applications to map the unexplored target space and therapeutic potential of natural products. Nat Prod Rep, 2015. 32(8): p. 1249-66.
Li, X.J., Z.Z. Jiang, and L.Y. Zhang, Triptolide: progress on research in pharmacodynamics and toxicology. J Ethnopharmacol, 2014. 155(1): p. 67-79.
Hopkins, A.L., Network pharmacology. Nat Biotechnol, 2007. 25(10): p. 1110-1.
Hopkins, A.L., Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol, 2008. 4(11): p. 682-90.
Li, S. and B. Zhang, Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med, 2013. 11(2): p. 110-20.
Zhang, B., X. Wang, and S. Li, An Integrative Platform of TCM Network Pharmacology and Its Application on a Herbal Formula, Qing-Luo-Yin. Evid Based Complement Alternat Med, 2013. 2013: p. 456747.
Zheng, C., et al., System-level multi-target drug discovery from natural products with applications to cardiovascular diseases. Mol Divers, 2014. 18(3): p. 621-35.
Zhou, W., et al., Systems Pharmacology in Small Molecular Drug Discovery. Int J Mol Sci, 2016. 17(2): p. 246.
Oliveira, A., et al., Triptolide abrogates growth of colon cancer and induces cell cycle arrest by inhibiting transcriptional activation of E2F. Lab Invest, 2015. 95(6): p. 648-659.
Liang, X., et al., Inhibition of RNA polymerase III transcription by Triptolide attenuates colorectal tumorigenesis. J Exp Clin Cancer Res, 2019. 38(1): p. 217.
Tong, X., et al., Triptolide inhibits cyclooxygenase-2 and inducible nitric oxide synthase expression in human colon cancer and leukemia cells. Acta Biochim Biophys Sin (Shanghai), 2007. 39(2): p. 89-95.

Download PDF

Version 1

posted

You are reading this latest preprint version

Prediction of triptolide targets in colorectal cancer using network pharmacology and molecular docking

Status:

Version 1

Abstract

Figures

1. Introduction

2. Materials And Methods

2.1. Identification of genes and proteins associated with CRC and potential triptolide target proteins

2.2. Prediction of pathways and networks associated with CRC and triptolide targets

2.3. Construction of common networks and core target screening

2.4. Prediction of binding between triptolide and candidate target proteins

3. Results

3.1. Networks and enriched functions in CRC-associated genes

3.2. Networks and enriched functions in triptolide target genes and proteins

3.3. Networks of shared proteins and special proteins targeted by triptolide

3.4. Predicted binding of triptolide to target proteins in CRC

4. Discussion

Declarations

Funding

Authors' contributions

Competing interests

References

Status:

Version 1