Neurofibromatosis type 1 (NF1) is a cancer predisposition disorder that results from inactivation of the tumor-suppressor Neurofibromin, a negative regulator of RAS signaling. NF1 patients present with a wide range of clinical manifestations and the tumor with highest prevalence is cutaneous neurofibroma (cNF). Most patients harboring cNF suffer greatly from the burden of those tumors, which have no effective medical treatment. Ironically, none of the numerous NF1 mouse models developed so far recapitulate cNF. Here, we discovered that Hoxb7 serves as a lineage marker to trace the developmental origin of cNF neoplastic cells. Ablating Nf1 in the Hoxb7 lineage faithfully recapitulates both human cutaneous and plexiform neurofibroma. In addition, we discovered that modulation of the Hippo pathway acts as a “modifier” for neurofibroma tumorigenesis. This mouse model opens the doors for deciphering the evolution of cNF to identify effective therapies, where none exist today.

Cutaneous neurofibromas (cNF) are a nearly ubiquitous symptom of neurofibromatosis type 1 (NF1), a disorder with a broad phenotypic spectrum caused by germline mutation of the neurofibromatosis type 1 tumour suppressor gene (NF1). Symptoms of NF1 can include learning disabilities, bone abnormalities and predisposition to tumours such as cNFs, plexiform neurofibromas, malignant peripheral nerve sheath tumours and optic nerve tumours. There are no therapies currently approved for cNFs aside from elective surgery, and the molecular aetiology of cNF remains relatively uncharacterised. Furthermore, whereas the biallelic inactivation of NF1 in neoplastic Schwann cells is critical for cNF formation, it is still unclear which additional genetic, transcriptional, epigenetic, microenvironmental or endocrine changes are important. Significant inroads have been made into cNF understanding, including NF1 genotype–phenotype correlations in NF1 microdeletion patients, the identification of recurring somatic mutations, studies of cNF-invading mast cells and macrophages, and clinical trials of putative therapeutic targets such as mTOR, MEK and c-KIT. Despite these advances, several gaps remain in our knowledge of the associated pathogenesis, which is further hampered by a lack of translationally relevant animal models. Some of these questions may be addressed in part by the adoption of genomic analysis techniques. Understanding the aetiology of cNF at the genomic level may assist in the development of new therapies for cNF, and may also contribute to a greater understanding of NF1/RAS signalling in cancers beyond those associated with NF1. Here, we summarise the present understanding of cNF biology, including the pathogenesis, mutational landscape, contribution of the tumour microenvironment and endocrine signalling, and the historical and current state of clinical trials for cNF. We also highlight open access data resources and potential avenues for future research that leverage recently developed genomics-based methods in cancer research.

N-of-1 trials target actionable mutations, yet such approaches do not test genomically-informed therapies in patient tumor models prior to patient treatment. To address this, we developed patient-derived xenograft (PDX) models from fine needle aspiration (FNA) biopsies (FNA-PDX) obtained from primary pancreatic ductal adenocarcinoma (PDAC) at the time of diagnosis. Here, we characterize PDX models established from one primary and two metastatic sites of one patient. We identified an activating KRAS G12R mutation among other mutations in these models. In explant cells derived from these PDX tumor models with a KRAS G12R mutation, treatment with inhibitors of CDKs (including CDK9) reduced phosphorylation of a marker of CDK9 activity (phospho-RNAPII CTD Ser2/5) and reduced viability/growth of explant cells derived from PDAC PDX models. Similarly, a CDK inhibitor reduced phospho-RNAPII CTD Ser2/5, increased apoptosis, and inhibited tumor growth in FNA-PDX and patient-matched metastatic-PDX models. In summary, PDX models can be constructed from FNA biopsies of PDAC which in turn can enable genomic characterization and identification of potential therapies.

Despite decades of intensive search for compounds that modulate the activity of particular protein targets, a large proportion of the human kinome remains as yet undrugged. Effective approaches are therefore required to map the massive space of unexplored compound–kinase interactions for novel and potent activities. Here, we carry out a crowdsourced benchmarking of predictive algorithms for kinase inhibitor potencies across multiple kinase families tested on unpublished bioactivity data. We find the top-performing predictions are based on various models, including kernel learning, gradient boosting and deep learning, and their ensemble leads to a predictive accuracy exceeding that of single-dose kinase activity assays. We design experiments based on the model predictions and identify unexpected activities even for under-studied kinases, thereby accelerating experimental mapping efforts. The open-source prediction algorithms together with the bioactivities between 95 compounds and 295 kinases provide a resource for benchmarking prediction algorithms and for extending the druggable kinome.

Summary The Columbia Cancer Target Discovery and Development (CTD2) Center is developing PANACEA, a resource comprising dose-responses and RNA sequencing (RNA-seq) profiles of 25 cell lines perturbed with ∼400 clinical oncology drugs, to study a tumor-specific drug mechanism of action. Here, this resource serves as the basis for a DREAM Challenge assessing the accuracy and sensitivity of computational algorithms for de novo drug polypharmacology predictions. Dose-response and perturbational profiles for 32 kinase inhibitors are provided to 21 teams who are blind to the identity of the compounds. The teams are asked to predict high-affinity binding targets of each compound among ∼1,300 targets cataloged in DrugBank. The best performing methods leverage gene expression profile similarity analysis as well as deep-learning methodologies trained on individual datasets. This study lays the foundation for future integrative analyses of pharmacogenomic data, reconciliation of polypharmacology effects in different tumor contexts, and insights into network-based assessments of drug mechanisms of action.

SUMMARYThe Columbia Cancer Target Discovery and Development (CTD2) Center has developed PANACEA (PANcancer Analysis of Chemical Entity Activity), a collection of dose-response curves and perturbational profiles for 400 clinical oncology drugs in cell lines selected to optimally represent 19 cancer subtypes. This resource, developed to study tumor-specific drug mechanism of action, was instrumental in hosting a DREAM Challenge to assess computational models for de novo drug polypharmacology prediction. Dose-response and perturbational profiles for 32 kinase inhibitors were provided to 21 participating teams, who did not know the identity or nature of the compounds, and they were asked to predict high-affinity binding among ~1,300 possible protein targets. Best performing methods leveraged both gene expression profile similarity analysis, and deep-learning methodologies trained on individual datasets. This study lays the foundation for future integrative analyses of pharmacogenomic data, reconciliation of polypharmacology effects in different tumor contexts, and insights into network-based assessment of context-specific drug mechanism of action.

BackgroundWe have identified molecules that exhibit synthetic lethality in cells with loss of the neurofibromin 1 (NF1) tumor suppressor gene. However, recognizing tumors that have inactivation of the NF1 tumor suppressor function is challenging because the loss may occur via mechanisms that do not involve mutation of the genomic locus. Degradation of the NF1 protein, independent of NF1 mutation status, phenocopies inactivating mutations to drive tumors in human glioma cell lines. NF1 inactivation may alter the transcriptional landscape of a tumor and allow a machine learning classifier to detect which tumors will benefit from synthetic lethal molecules.ResultsWe developed a strategy to predict tumors with low NF1 activity and hence tumors that may respond to treatments that target cells lacking NF1. Using RNAseq data from The Cancer Genome Atlas (TCGA), we trained an ensemble of 500 logistic regression classifiers that integrates mutation status with whole transcriptomes to predict NF1 inactivation in glioblastoma (GBM). On TCGA data, the classifier detected NF1 mutated tumors (test set area under the receiver operating characteristic curve (AUROC) mean = 0.77, 95% quantile = 0.53 – 0.95) over 50 random initializations. On RNA-Seq data transformed into the space of gene expression microarrays, this method produced a classifier with similar performance (test set AUROC mean = 0.77, 95% quantile = 0.53 – 0.96). We applied our ensemble classifier trained on the transformed TCGA data to a microarray validation set of 12 samples with matched RNA and NF1 protein-level measurements. The classifier’s NF1 score was associated with NF1 protein concentration in these samples.ConclusionsWe demonstrate that TCGA can be used to train accurate predictors of NF1 inactivation in GBM. The ensemble classifier performed well for samples with very high or very low NF1 protein concentrations but had mixed performance in samples with intermediate NF1 concentrations. Nevertheless, high-performing and validated predictors have the potential to be paired with targeted therapies and personalized medicine.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3519-7) contains supplementary material, which is available to authorized users.

Despite decades of intensive search for compounds that modulate the activity of particular proteins, there are currently small-molecule probes available only for a small proportion of the human proteome. Effective approaches are therefore required to map the massive space of unexplored compound-target interactions for novel and potent activities. Here, we carried out a crowdsourced benchmarking of the accuracy of machine learning (ML) algorithms at predicting kinase inhibitor potencies across multiple kinase families. A total of 268 ML predictions were scored in unpublished bioactivity data sets. Top-performing algorithms used kernel learning, gradient boosting and deep learning, with predictive accuracy exceeding that of target activity assays. Subsequent experiments carried out based on the the top-performing model predictions demonstrated that these models and their ensemble can improve the accuracy of experimental mapping efforts, especially for so far under-studied kinases. The open-source ML algorithms together with the novel dose-response data for 905 bioactivities between 95 compounds and 295 kinases provide a unique resource for extending the druggable kinome.AI approaches, that can leverage the information extracted from similar kinases and compounds to predict the activity of so far unexplored interactions.The Challenge was implemented in a screening-based, pre-competitive drug discovery project in collaboration with the NIH-supported Illuminating the Druggable Genome (IDG) program ( https://commonfund.nih.gov/idg ), with the common aim to establish kinome-wide target profiles of small-molecule agents, and thereby to extend the druggability of the human kinome space by providing activity information on under-studied proteins. The specific questions this Challenge sought to address were: (i) What are the best computational modelling approaches for predicting quantitative compound-target activity profiles?; (ii) What are the optimal molecular and chemical descriptors for maximal prediction accuracy?; and (iii) What are the most predictive bioactivity assays and publicly available datasets? The Challenge attracted 212 active participants, and a total of 268 predictions were scored, covering a wide range of ML approaches, including deep and kernel learning and gradient boosting decision trees. Here, we describe the benchmarking results from the Challenge, and the use of top-performing models for identifying novel kinase inhibitor activities.