Cancer In Silico Drug Discovery: A Systems Biology Tool for Identifying Candidate Drugs to Target Specific Molecular Tumor Subtypes

Large-scale cancer datasets such as The Cancer Genome Atlas (TCGA) allow researchers to profile tumors based on a wide range of clinical and molecular characteristics. Subsequently, TCGA-derived gene expression profiles can be analyzed with the Connectivity Map (CMap) to find candidate drugs to target tumors with specific clinical phenotypes or molecular characteristics. This represents a powerful computational approach for candidate drug identification, but due to the complexity of TCGA and technology differences between CMap and TCGA experiments, such analyses are challenging to conduct and reproduce. We present Cancer in silico Drug Discovery (CiDD; scheet.org/software), a computational drug discovery platform that addresses these challenges. CiDD integrates data from TCGA, CMap, and Cancer Cell Line Encyclopedia (CCLE) to perform computational drug discovery experiments, generating hypotheses for the following three general problems: (i) determining whether specific clinical phenotypes or molecular characteristics are associated with unique gene expression signatures; (ii) finding candidate drugs to repress these expression signatures; and (iii) identifying cell lines that resemble the tumors being studied for subsequent in vitro experiments. The primary input to CiDD is a clinical or molecular characteristic. The output is a biologically annotated list of candidate drugs and a list of cell lines for in vitro experimentation. We applied CiDD to identify candidate drugs to treat colorectal cancers harboring mutations in BRAF. CiDD identified EGFR and proteasome inhibitors, while proposing five cell lines for in vitro testing. CiDD facilitates phenotype-driven, systematic drug discovery based on clinical and molecular data from TCGA. Mol Cancer Ther; 13(12); 3230–40. ©2014 AACR.

[1]  Andreas Schlicker,et al.  Colorectal cancer intrinsic subtypes predict chemotherapy benefit, deficient mismatch repair and epithelial-to-mesenchymal transition , 2013, International journal of cancer.

[2]  Susumu Goto,et al.  Data, information, knowledge and principle: back to metabolism in KEGG , 2013, Nucleic Acids Res..

[3]  N. Rosen,et al.  Tumor adaptation and resistance to RAF inhibitors , 2013, Nature Medicine.

[4]  Tae-Min Kim,et al.  Clinical applications of next-generation sequencing in colorectal cancers. , 2013, World journal of gastroenterology.

[5]  Philippe Bédard,et al.  Clinical application of high-throughput genomic technologies for treatment selection in breast cancer , 2013, Breast Cancer Research.

[6]  D. Zecchin,et al.  BRAF V600E Is a Determinant of Sensitivity to Proteasome Inhibitors , 2013, Molecular Cancer Therapeutics.

[7]  C. Sander,et al.  Evaluating cell lines as tumour models by comparison of genomic profiles , 2013, Nature Communications.

[8]  M. Salto‐Tellez,et al.  Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data , 2013, PloS one.

[9]  Heidi L. Rehm,et al.  Disease-targeted sequencing: a cornerstone in the clinic , 2013, Nature Reviews Genetics.

[10]  I. Nookaew,et al.  Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods , 2013, Nucleic acids research.

[11]  G. Orphanides,et al.  Subtypes of primary colorectal tumors correlate with response to targeted treatment in colorectal cell lines , 2012, BMC Medical Genomics.

[12]  Michael A. Davies,et al.  Resistance to BRAF Inhibition in BRAF-Mutant Colon Cancer Can Be Overcome with PI3K Inhibition or Demethylating Agents , 2012, Clinical Cancer Research.

[13]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of human colon and rectal cancer , 2012, Nature.

[14]  Chunhua Wang,et al.  Analyzing Gene Expression Profile in K562 Cells Exposed to Sodium Valproate Using Microarray Combined with the Connectivity Map Database , 2012, Journal of biomedicine & biotechnology.

[15]  Gongping Yang,et al.  Finger Vein Recognition Based on (2D)2 PCA and Metric Learning , 2012, Journal of biomedicine & biotechnology.

[16]  Sabine Tejpar,et al.  Identification of a poor-prognosis BRAF-mutant-like population of patients with colon cancer. , 2012, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[17]  Adam A. Margolin,et al.  The Cancer Cell Line Encyclopedia enables predictive modeling of anticancer drug sensitivity , 2012, Nature.

[18]  Mari Mino-Kenudson,et al.  EGFR-mediated re-activation of MAPK signaling contributes to insensitivity of BRAF mutant colorectal cancers to RAF inhibition with vemurafenib. , 2012, Cancer discovery.

[19]  Yu Shyr,et al.  Survival in BRAF V600-mutant advanced melanoma treated with vemurafenib. , 2012, The New England journal of medicine.

[20]  S. Kopetz,et al.  Antitumor activity of BRAF inhibitor vemurafenib in preclinical models of BRAF-mutant colorectal cancer. , 2012, Cancer research.

[21]  R. Bernards,et al.  Unresponsiveness of colon cancer to BRAF(V600E) inhibition through feedback activation of EGFR , 2012, Nature.

[22]  P. Gibbs,et al.  Impact of BRAF mutation and microsatellite instability on the pattern of metastatic spread and prognosis in metastatic colorectal cancer , 2011, Cancer.

[23]  Yuan Hu,et al.  Protection of SH-SY5Y Neuronal Cells from Glutamate-Induced Apoptosis by 3,6′-Disinapoyl Sucrose, a Bioactive Compound Isolated from Radix Polygala , 2011, Journal of biomedicine & biotechnology.

[24]  A. Hauschild,et al.  Improved survival with vemurafenib in melanoma with BRAF V600E mutation. , 2011, The New England journal of medicine.

[25]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[26]  J. O’Leary,et al.  BRAFV600E: Implications for Carcinogenesis and Molecular Therapy , 2011, Molecular Cancer Therapeutics.

[27]  David S. Wishart,et al.  DrugBank 3.0: a comprehensive resource for ‘Omics’ research on drugs , 2010, Nucleic Acids Res..

[28]  S. Sivapalaratnam,et al.  More outreach for young scientists. , 2010, Nature.

[29]  K. Baggerly Disclose all data in publications. , 2010, Nature.

[30]  J. Desai,et al.  PLX4032 in metastatic colorectal cancer patients with mutant BRAF tumors. , 2010 .

[31]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[32]  Ultan McDermott,et al.  Personalized cancer therapy with selective kinase inhibitors: an emerging paradigm in medical oncology. , 2009, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[33]  Leyla Isik,et al.  Cancer-specific high-throughput annotation of somatic mutations: computational prediction of driver missense mutations. , 2009, Cancer research.

[34]  Sampsa Hautaniemi,et al.  Deciphering downstream gene targets of PI3K/mTOR/p70S6K pathway in breast cancer , 2008, BMC Genomics.

[35]  W. Gullick,et al.  Development of breast cancer immunotherapy using MUC1-retargeted T lymphocytes , 2006, Breast Cancer Research.

[36]  A. Lièvre,et al.  KRAS mutations as an independent prognostic factor in patients with advanced colorectal cancer treated with cetuximab. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[37]  Robert B. Russell,et al.  SuperTarget and Matador: resources for exploring drug-target relationships , 2007, Nucleic Acids Res..

[38]  J. Davis Bioinformatics and Computational Biology Solutions Using R and Bioconductor , 2007 .

[39]  Paul A Clemons,et al.  The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease , 2006, Science.

[40]  S. Kakar,et al.  BRAF mutation, CpG island methylator phenotype and microsatellite instability occur more frequently and concordantly in mucinous than non‐mucinous colorectal cancer , 2006, International journal of cancer.

[41]  D. Geman,et al.  Simple decision rules for classifying human cancers from gene expression profiles , 2005, Bioinform..

[42]  P. Hall,et al.  An expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Rafael A. Irizarry,et al.  Bioinformatics and Computational Biology Solutions using R and Bioconductor , 2005 .

[44]  Carlo Rago,et al.  Mutant PIK3CA promotes cell growth and invasion of human cancer cells. , 2005, Cancer cell.

[45]  H. Varmus,et al.  KRAS Mutations and Primary Resistance of Lung Adenocarcinomas to Gefitinib or Erlotinib , 2005, PLoS medicine.

[46]  Steven Henikoff,et al.  SIFT: predicting amino acid changes that affect protein function , 2003, Nucleic Acids Res..

[47]  A. Nicholson,et al.  Mutations of the BRAF gene in human cancer , 2002, Nature.

[48]  Robert Gentleman,et al.  gCMAP: user-friendly connectivity mapping with R , 2014, Bioinform..

[49]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .