Deep Learning and Association Rule Mining for Predicting Drug Response in Cancer. A Personalised Medicine Approach

A major challenge in cancer treatment is predicting the clinical response to anticancer drugs for each individual patient. For complex diseases, such as cancer, characterized by high inter-patient variance, the implementation of precision medicine approaches is dependent upon understanding the disease process at the molecular level. While the ‘omics’ era provides unique opportunities to dissect the molecular features of diseases, the ability to apply it to targeted therapeutic efforts is hindered by both the massive size and diverse nature of the ‘omic’ data. Recent advances with Deep Learning Neural Networks (DLNN), suggests that DLNN could be trained on large data sets to efficiently predict therapeutic responses. We present the application of Association Rule Mining (Market Basket Analysis) in combination with Deep Learning to integrate and extract knowledge in the form of easily interpretable rules from the molecular profiles of 689 cancer cell lines and predict pharmacological responses to 139 anti-cancer drugs. The proposed algorithm achieved superior classification and outperformed Random Forests which currently represents the state-of-the-art classification process. Finally, the in silico pipeline presented introduces a novel strategy for identifying drug combinations with high therapeutic potential. For the first time, we demonstrate that DLNN trained on a large pharmacogenomic data set can effectively predict the therapeutic response of specific drugs in specific cancer types, from a large panel of both drugs and cancer cell lines. These findings serve as a proof of concept for the application of DLNN to predict therapeutic responsiveness, a milestone in precision medicine.

[1]  Guigang Zhang,et al.  Deep Learning , 2016, Int. J. Semantic Comput..

[2]  Sergey Plis,et al.  Deep Learning Applications for Predicting Pharmacological Properties of Drugs and Drug Repurposing Using Transcriptomic Data. , 2016, Molecular pharmaceutics.

[3]  A. Anichini,et al.  Towards combinatorial targeted therapy in melanoma: From pre-clinical evidence to clinical application (Review) , 2014, International journal of oncology.

[4]  Nitish Srivastava,et al.  Dropout: a simple way to prevent neural networks from overfitting , 2014, J. Mach. Learn. Res..

[5]  H. Zou,et al.  Regularization and variable selection via the elastic net , 2005 .

[6]  R. Karchin,et al.  Collections of simultaneously altered genes as biomarkers of cancer cell drug response. , 2013, Cancer research.

[7]  S. Mustjoki,et al.  Novel drug candidates for blast phase chronic myeloid leukemia from high-throughput drug sensitivity and resistance testing , 2015, Blood Cancer Journal.

[8]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[9]  Leo Breiman,et al.  Bagging Predictors , 1996, Machine Learning.

[10]  C. Metz Basic principles of ROC analysis. , 1978, Seminars in nuclear medicine.

[11]  D. Haber,et al.  Cell line-based platforms to evaluate the therapeutic efficacy of candidate anticancer agents , 2010, Nature Reviews Cancer.

[12]  R. Bernards,et al.  Enabling personalized cancer medicine through analysis of gene-expression patterns , 2008, Nature.

[13]  A. Schally,et al.  Antagonistic analogs of growth hormone-releasing hormone increase the efficacy of treatment of triple negative breast cancer in nude mice with doxorubicin; A preclinical study , 2014, Oncoscience.

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

[15]  Thomas Beckers,et al.  Histone Deacetylase Inhibitors Suppress the Inducibility of Nuclear Factor-κB by Tumor Necrosis Factor-α Receptor-1 Down-regulation , 2006 .

[16]  Di Chen,et al.  Systematic synergy modeling: understanding drug synergy from a systems biology perspective , 2015, BMC Systems Biology.

[17]  Laura M. Heiser,et al.  A community effort to assess and improve drug sensitivity prediction algorithms , 2014, Nature Biotechnology.

[18]  A. Schally,et al.  Suppression of tumor growth by growth hormone-releasing hormone antagonist JV-1-36 does not involve the inhibition of autocrine production of insulin-like growth factor II in H-69 small cell lung carcinoma. , 2000, Cancer letters.

[19]  S Miyano,et al.  Addiction to the IGF2-ID1-IGF2 circuit for maintenance of the breast cancer stem-like cells , 2016, Oncogene.

[20]  R. Maxwell,et al.  The synergistic interaction of MEK and PI3K inhibitors is modulated by mTOR inhibition , 2012, British Journal of Cancer.

[21]  Mehmet Gönen,et al.  Drug susceptibility prediction against a panel of drugs using kernelized Bayesian multitask learning , 2014, Bioinform..

[22]  X. Chen,et al.  Random forests for genomic data analysis. , 2012, Genomics.

[23]  J. Falkenburg,et al.  Differential Effects of Selective Inhibitors Targeting the PI3K/AKT/mTOR Pathway in Acute Lymphoblastic Leukemia , 2013, PloS one.

[24]  Jürgen Schmidhuber,et al.  Deep learning in neural networks: An overview , 2014, Neural Networks.

[25]  山田 英貴 Detailed characterization of a homozygously deleted region corresponding to a candidate tumor suppressor locus at 21q11-21 in human lung cancer , 2008 .

[26]  Dai Shimizu,et al.  Prognostic relevance of SAMSN1 expression in gastric cancer. , 2016, Oncology letters.

[27]  R. Geyer,et al.  In Human Leukemia Cells Ephrin-B–Induced Invasive Activity Is Supported by Lck and Is Associated with Reassembling of Lipid Raft Signaling Complexes , 2008, Molecular Cancer Research.

[28]  L. Augenlicht,et al.  HDAC2 deficiency sensitizes colon cancer cells to TNFalpha-induced apoptosis through inhibition of NF-kappaB activity. , 2008, Experimental cell research.

[29]  M. Zloh,et al.  Natural chalcones as dual inhibitors of HDACs and NF-κB , 2012, Oncology reports.

[30]  A. Sami,et al.  Targeting the PI3K/AKT/mTOR signaling pathway in glioblastoma: novel therapeutic agents and advances in understanding , 2013, Tumor Biology.

[31]  A. Villa,et al.  The POF1B candidate gene for premature ovarian failure regulates epithelial polarity , 2011, Journal of Cell Science.

[32]  Thomas Beckers,et al.  Histone deacetylase inhibitors suppress the inducibility of nuclear factor-kappaB by tumor necrosis factor-alpha receptor-1 down-regulation. , 2006, Cancer research.

[33]  A. Martin,et al.  Assessment of epidermal growth factor receptor (EGFR, ErbB1) and HER2 (ErbB2) protein expression levels and response to lapatinib (Tykerb®, GW572016) in an expanded panel of human normal and tumour cell lines , 2007, Cell proliferation.

[34]  S. Ramaswamy,et al.  Systematic identification of genomic markers of drug sensitivity in cancer cells , 2012, Nature.

[35]  Y. Zhong,et al.  Overexpression of a transcription factor LYL1 induces T- and B-cell lymphoma in mice , 2007, Oncogene.

[36]  Xiaohua Ma,et al.  Mechanisms of drug combinations: interaction and network perspectives , 2009, Nature Reviews Drug Discovery.

[37]  Stefano Volinia,et al.  Mutational Landscape and Gene Expression Patterns in Adult Acute Myeloid Leukemias with Monosomy 7 as a Sole Abnormality. , 2017, Cancer research.

[38]  L. Lou,et al.  NL-103, a novel dual-targeted inhibitor of histone deacetylases and hedgehog pathway, effectively overcomes vismodegib resistance conferred by Smo mutations , 2014, Pharmacology research & perspectives.

[39]  Fangfang Shi,et al.  Ribavirin Inhibits the Activity of mTOR/eIF4E, ERK/Mnk1/eIF4E Signaling Pathway and Synergizes with Tyrosine Kinase Inhibitor Imatinib to Impair Bcr-Abl Mediated Proliferation and Apoptosis in Ph+ Leukemia , 2015, PloS one.

[40]  Jeffrey Heer,et al.  SpanningAspectRatioBank Easing FunctionS ArrayIn ColorIn Date Interpolator MatrixInterpola NumObjecPointI Rectang ISchedu Parallel Pause Scheduler Sequen Transition Transitioner Transiti Tween Co DelimGraphMLCon IData JSONCon DataField DataSc Dat DataSource Data DataUtil DirtySprite LineS RectSprite , 2011 .

[41]  P. Kranz,et al.  Depletion of the thiol oxidoreductase ERp57 in tumor cells inhibits proliferation and increases sensitivity to ionizing radiation and chemotherapeutics , 2015, Oncotarget.

[42]  Thomas Lengauer,et al.  ROCR: visualizing classifier performance in R , 2005, Bioinform..

[43]  Julio Saez-Rodriguez,et al.  Machine Learning Prediction of Cancer Cell Sensitivity to Drugs Based on Genomic and Chemical Properties , 2012, PloS one.

[44]  Pascal Vincent,et al.  Representation Learning: A Review and New Perspectives , 2012, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[45]  James S Blachly,et al.  Targeting PI3‐kinase (PI3K), AKT and mTOR axis in lymphoma , 2014, British journal of haematology.

[46]  Atsushi Imiya,et al.  Machine Learning and Data Mining in Pattern Recognition , 2013, Lecture Notes in Computer Science.

[47]  O. Myklebost,et al.  Potential for treatment of liposarcomas with the MDM2 antagonist Nutlin‐3A , 2007, International journal of cancer.

[48]  K. Kiura,et al.  Synergistic Effect of Olaparib with Combination of Cisplatin on PTEN-Deficient Lung Cancer Cells , 2012, Molecular Cancer Research.

[49]  C. Baldus,et al.  Outlook on PI3K/AKT/mTOR inhibition in acute leukemia , 2015, Molecular and Cellular Therapies.

[50]  William Stafford Noble,et al.  Machine learning applications in genetics and genomics , 2015, Nature Reviews Genetics.

[51]  K. Vandyke,et al.  SAMSN1 Is a Tumor Suppressor Gene in Multiple Myeloma12 , 2014, Neoplasia.

[52]  David Craft,et al.  The value of prior knowledge in machine learning of complex network systems , 2016 .

[53]  Tomasz Imielinski,et al.  Mining association rules between sets of items in large databases , 1993, SIGMOD Conference.

[54]  M. A. King,et al.  Inhibition of cholesterol metabolism underlies synergy between mTOR pathway inhibition and chloroquine in bladder cancer cells , 2016, Oncogene.

[55]  Paula T Hammond,et al.  Layer-by-layer nanoparticles for systemic codelivery of an anticancer drug and siRNA for potential triple-negative breast cancer treatment. , 2013, ACS nano.

[56]  T G Myers,et al.  DT-Diaphorase expression and tumor cell sensitivity to 17-allylamino, 17-demethoxygeldanamycin, an inhibitor of heat shock protein 90. , 1999, Journal of the National Cancer Institute.

[57]  C. Dickey,et al.  Potential synergy between tau aggregation inhibitors and tau chaperone modulators , 2013, Alzheimer's Research & Therapy.

[58]  R. Lo Receptor tyrosine kinases in cancer escape from BRAF inhibitors , 2012, Cell Research.

[59]  Nci Dream Community A community effort to assess and improve drug sensitivity prediction algorithms , 2014 .

[60]  J. Weinstein Drug discovery: Cell lines battle cancer , 2012, Nature.

[61]  R. Place,et al.  HDAC inhibition prevents NF-kappa B activation by suppressing proteasome activity: down-regulation of proteasome subunit expression stabilizes I kappa B alpha. , 2005, Biochemical pharmacology.

[62]  A. Giassi,et al.  Immunohistochemistry in the Diagnosis of Mucinous Neoplasms Involving the Ovary: The Added Value of SATB2 and Biomarker Discovery Through Protein Expression Database Mining , 2016, International journal of gynecological pathology : official journal of the International Society of Gynecological Pathologists.

[63]  King-Sun Fu,et al.  IEEE Transactions on Pattern Analysis and Machine Intelligence Publication Information , 2004, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[64]  Nazanin Heidari,et al.  Targeting PI3K/AKT/mTOR network for treatment of leukemia , 2015, Cellular and Molecular Life Sciences.

[65]  G. Zaman,et al.  Selective Targeting of CTNNB1-, KRAS- or MYC-Driven Cell Growth by Combinations of Existing Drugs , 2015, PloS one.

[66]  P. Jaccard THE DISTRIBUTION OF THE FLORA IN THE ALPINE ZONE.1 , 1912 .

[67]  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.

[68]  S. Ramaswamy,et al.  Assessment of ABT-263 activity across a cancer cell line collection leads to a potent combination therapy for small-cell lung cancer , 2015, Proceedings of the National Academy of Sciences.

[69]  R. Place,et al.  HDAC inhibition prevents NF-κB activation by suppressing proteasome activity: Down-regulation of proteasome subunit expression stabilizes IκBα , 2005 .

[70]  J. Dick,et al.  Oncogenic potential of the transcription factor LYL1 in acute myeloblastic leukemia , 2005, Leukemia.

[71]  G. Hannon,et al.  Patient-derived tumor xenografts: transforming clinical samples into mouse models. , 2013, Cancer research.

[72]  Peter Dalgaard,et al.  R Development Core Team (2010): R: A language and environment for statistical computing , 2010 .

[73]  W. Youden,et al.  Index for rating diagnostic tests , 1950, Cancer.

[74]  Dimitris Kletsas,et al.  Chronic p53-independent p21 expression causes genomic instability by deregulating replication licensing , 2016, Nature Cell Biology.

[75]  Jaideep Srivastava,et al.  Selecting the right objective measure for association analysis , 2004, Inf. Syst..

[76]  Takeshi Inoue,et al.  The Use of Olaparib (AZD2281) Potentiates SN-38 Cytotoxicity in Colon Cancer Cells by Indirect Inhibition of Rad51-Mediated Repair of DNA Double-Strand Breaks , 2014, Molecular Cancer Therapeutics.

[77]  Hidetoshi Shimodaira,et al.  Pvclust: an R package for assessing the uncertainty in hierarchical clustering , 2006, Bioinform..

[78]  Emanuel J. V. Gonçalves,et al.  A Landscape of Pharmacogenomic Interactions in Cancer , 2016, Cell.

[79]  Mingming Jia,et al.  COSMIC: exploring the world's knowledge of somatic mutations in human cancer , 2014, Nucleic Acids Res..

[80]  R. Shoemaker The NCI60 human tumour cell line anticancer drug screen , 2006, Nature Reviews Cancer.

[81]  A. Knudson Mutation and cancer: statistical study of retinoblastoma. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[82]  Shumei Zhai,et al.  P-glycoprotein-evading anti-tumor activity of a novel tubulin and HSP90 dual inhibitor in a non-small-cell lung cancer model. , 2014, Journal of pharmacological sciences.

[83]  Michael Hahsler,et al.  Visualizing association rules in hierarchical groups , 2016, Journal of Business Economics.

[84]  Luhua Lai,et al.  Deep Learning for Drug-Induced Liver Injury , 2015, J. Chem. Inf. Model..

[85]  D. Cho,et al.  PI3K inhibition potentiates Bcl-2-dependent apoptosis in renal carcinoma cells , 2013, Journal of cellular and molecular medicine.

[86]  J. Vázquez-Prado,et al.  Regulation of mTORC1 Complex Assembly and Signaling by GRp58/ERp57 , 2011, Molecular and Cellular Biology.

[87]  Petr Konečný,et al.  Pathophysiologically relevant in vitro tumor models for drug screening. , 2015, Drug discovery today.

[88]  M. Mandal,et al.  PI3K and Akt as molecular targets for cancer therapy: current clinical outcomes , 2012, Acta Pharmacologica Sinica.

[89]  Jiri Bartek,et al.  An Oncogene-Induced DNA Damage Model for Cancer Development , 2008, Science.

[90]  Y. Hirasawa,et al.  Actin organization associated with the expression of multidrug resistant phenotype in osteosarcoma cells and the effect of actin depolymerization on drug resistance. , 1998, Cancer letters.

[91]  B. Matthews Comparison of the predicted and observed secondary structure of T4 phage lysozyme. , 1975, Biochimica et biophysica acta.

[92]  Rong Qin,et al.  SAMSN1 Is Highly Expressed and Associated with a Poor Survival in Glioblastoma Multiforme , 2013, PloS one.

[93]  Dai Shimizu,et al.  Suppression of SAMSN1 Expression is Associated with the Malignant Phenotype of Hepatocellular Carcinoma , 2015, Annals of Surgical Oncology.

[94]  E. Kohn,et al.  The MAPK pathway across different malignancies: A new perspective , 2014, Cancer.

[95]  Bonnie E. Shook-Sa,et al.  . CC-BY-NC-ND 4 . 0 International licenseIt is made available under a is the author / funder , who has granted medRxiv a license to display the preprint in perpetuity , 2021 .

[96]  F. Finkernagel,et al.  Histone Deacetylase 6 Represents a Novel Drug Target in the Oncogenic Hedgehog Signaling Pathway , 2014, Molecular Cancer Therapeutics.

[97]  J. McCubrey,et al.  Cytotoxic activity of the novel Akt inhibitor, MK-2206, in T-cell acute lymphoblastic leukemia , 2012, Leukemia.

[98]  Yan Wu,et al.  Interaction of the Tumor Suppressor PTEN/MMAC with a PDZ Domain of MAGI3, a Novel Membrane-associated Guanylate Kinase* , 2000, The Journal of Biological Chemistry.

[99]  Pier Paolo Pandolfi,et al.  The PTEN–PI3K pathway: of feedbacks and cross-talks , 2008, Oncogene.

[100]  Ji-Eun Jung,et al.  The ID1-CULLIN3 Axis Regulates Intracellular SHH and WNT Signaling in Glioblastoma Stem Cells. , 2016, Cell reports.

[101]  Víctor J Cid,et al.  A comprehensive functional analysis of PTEN mutations: implications in tumor- and autism-related syndromes. , 2011, Human molecular genetics.

[102]  Lionel Feigenbaum,et al.  Smad3 reduces susceptibility to hepatocarcinoma by sensitizing hepatocytes to apoptosis through downregulation of Bcl-2. , 2006, Cancer cell.

[103]  Wei Li,et al.  An essential role for the Id1/PI3K/Akt/NFkB/survivin signalling pathway in promoting the proliferation of endothelial progenitor cells in vitro , 2011, Molecular and Cellular Biochemistry.

[104]  C. Datz,et al.  Synergistic inhibition of growth of lung carcinomas by antagonists of growth hormone-releasing hormone in combination with docetaxel , 2006, Proceedings of the National Academy of Sciences.

[105]  Dimitris Kletsas,et al.  Cytotoxicity and estrogenicity of Invisalign appliances. , 2009, American journal of orthodontics and dentofacial orthopedics : official publication of the American Association of Orthodontists, its constituent societies, and the American Board of Orthodontics.

[106]  J. Utikal,et al.  Combined BRAF and MEK inhibition versus BRAF inhibition alone in melanoma. , 2014, The New England journal of medicine.

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

[108]  Hiroyuki Hamakawa,et al.  Suppression of invasion and metastasis in aggressive salivary cancer cells through targeted inhibition of ID1 gene expression. , 2016, Cancer letters.

[109]  Kurt Hornik,et al.  Multilayer feedforward networks are universal approximators , 1989, Neural Networks.

[110]  Naren L. Banik,et al.  Drug Resistance in Glioblastoma: A Mini Review , 2012, Neurochemical Research.

[111]  W. Hsu,et al.  Comparison of Chemotherapy Response with P-Glycoprotein, Multidrug Resistance-Related Protein-1, and Lung Resistance-Related Protein Expression in Untreated Small Cell Lung Cancer , 2005, Lung.

[112]  J. Downing,et al.  Loss of oncogenic Notch1 with resistance to a PI3K inhibitor in T-cell leukaemia , 2014, Nature.

[113]  Masato Hareyama,et al.  The combination of olaparib and camptothecin for effective radiosensitization , 2012, Radiation oncology.

[114]  Thomas G. Dietterich What is machine learning? , 2020, Archives of Disease in Childhood.