Although mechanisms of acquired resistance of epidermal growth factor receptor (EGFR)-mutant non-small-cell lung cancers to EGFR inhibitors have been identified, little is known about how resistant clones evolve during drug therapy. Here we observe that acquired resistance caused by the EGFR(T790M) gatekeeper mutation can occur either by selection of pre-existing EGFR(T790M)-positive clones or via genetic evolution of initially EGFR(T790M)-negative drug-tolerant cells. The path to resistance impacts the biology of the resistant clone, as those that evolved from drug-tolerant cells had a diminished apoptotic response to third-generation EGFR inhibitors that target EGFR(T790M); treatment with navitoclax, an inhibitor of the anti-apoptotic factors BCL-xL and BCL-2 restored sensitivity. We corroborated these findings using cultures derived directly from EGFR inhibitor-resistant patient tumors. These findings provide evidence that clinically relevant drug-resistant cancer cells can both pre-exist and evolve from drug-tolerant cells, and they point to therapeutic opportunities to prevent or overcome resistance in the clinic.

The anaplastic lymphoma kinase (ALK) fusion oncogene is observed in 3%-5% of non-small cell lung cancer (NSCLC). Crizotinib and ceritinib, a next-generation ALK tyrosine kinase inhibitor (TKI) active against crizotinib-refractory patients, are clinically available for the treatment of ALK-rearranged NSCLC patients, and multiple next-generation ALK-TKIs are currently under clinical evaluation. These ALK-TKIs exhibit robust clinical activity in ALK-rearranged NSCLC patients; however, the emergence of ALK-TKI resistance restricts the therapeutic effect. To date, various secondary mutations or bypass pathway activation-mediated resistance have been identified, but large parts of the resistance mechanism are yet to be identified. Here, we report the discovery of p-glycoprotein (P-gp/ABCB1) overexpression as a ceritinib resistance mechanism in ALK-rearranged NSCLC patients. P-gp exported ceritinib and its overexpression conferred ceritinib and crizotinib resistance, but not to PF-06463922 or alectinib, which are next-generation ALK inhibitors. Knockdown of ABCB1 or P-gp inhibitors sensitizes the patient-derived cancer cells to ceritinib, in vitro and in vivo. P-gp overexpression was identified in three out of 11 cases with in ALK-rearranged crizotinib or ceritinib resistant NSCLC patients. Our study suggests that alectinib, PF-06463922, or P-gp inhibitor with ceritinib could overcome the ceritinib or crizotinib resistance mediated by P-gp overexpression.

Tumor mutational burden correlates with response to immune checkpoint blockade in multiple solid tumors, although in microsatellite-stable tumors this association is of uncertain clinical utility. Here we uniformly analyzed whole-exome sequencing (WES) of 249 tumors and matched normal tissue from patients with clinically annotated outcomes to immune checkpoint therapy, including radiographic response, across multiple cancer types to examine additional tumor genomic features that contribute to selective response. Our analyses identified genomic correlates of response beyond mutational burden, including somatic events in individual driver genes, certain global mutational signatures, and specific HLA-restricted neoantigens. However, these features were often interrelated, highlighting the complexity of identifying genetic driver events that generate an immunoresponsive tumor environment. This study lays a path forward in analyzing large clinical cohorts in an integrated and multifaceted manner to enhance the ability to discover clinically meaningful predictive features of response to immune checkpoint blockade.

Cancer stem cells (CSCs) correlate with recurrence, metastasis and poor survival in clinical studies. Encouraging results from clinical trials of CSC inhibitors have further validated CSCs as therapeutic targets. ONC201 is a first-in-class small molecule imipridone in Phase I/II clinical trials for advanced cancer. We have previously shown that ONC201 targets self-renewing, chemotherapy-resistant colorectal CSCs via Akt/ERK inhibition and DR5/TRAIL induction. In this study, we demonstrate that the anti-CSC effects of ONC201 involve early changes in stem cell-related gene expression prior to tumor cell death induction. A targeted network analysis of gene expression profiles in colorectal cancer cells revealed that ONC201 downregulates stem cell pathways such as Wnt signaling and modulates genes (ID1, ID2, ID3 and ALDH7A1) known to regulate self-renewal in colorectal, prostate cancer and glioblastoma. ONC201-mediated changes in CSC-related gene expression were validated at the RNA and protein level for each tumor type. Accordingly, we observed inhibition of self-renewal and CSC markers in prostate cancer cell lines and patient-derived glioblastoma cells upon ONC201 treatment. Interestingly, ONC201-mediated CSC depletion does not occur in colorectal cancer cells with acquired resistance to ONC201. Finally, we observed that basal expression of CSC-related genes (ID1, CD44, HES7 and TCF3) significantly correlate with ONC201 efficacy in >1000 cancer cell lines and combining the expression of multiple genes leads to a stronger overall prediction. These proof-of-concept studies provide a rationale for testing CSC expression at the RNA and protein level as a predictive and pharmacodynamic biomarker of ONC201 response in ongoing clinical studies.

Predicting the response of a specific cancer to a therapy is a major goal in modern oncology that should ultimately lead to a personalised treatment. High-throughput screenings of potentially active compounds against a panel of genomically heterogeneous cancer cell lines have unveiled multiple relationships between genomic alterations and drug responses. Various computational approaches have been proposed to predict sensitivity based on genomic features, while others have used the chemical properties of the drugs to ascertain their effect. In an effort to integrate these complementary approaches, we developed machine learning models to predict the response of cancer cell lines to drug treatment, quantified through IC₅₀ values, based on both the genomic features of the cell lines and the chemical properties of the considered drugs. Models predicted IC₅₀ values in a 8-fold cross-validation and an independent blind test with coefficient of determination R² of 0.72 and 0.64 respectively. Furthermore, models were able to predict with comparable accuracy (R² of 0.61) IC50s of cell lines from a tissue not used in the training stage. Our in silico models can be used to optimise the experimental design of drug-cell screenings by estimating a large proportion of missing IC₅₀ values rather than experimentally measuring them. The implications of our results go beyond virtual drug screening design: potentially thousands of drugs could be probed in silico to systematically test their potential efficacy as anti-tumour agents based on their structure, thus providing a computational framework to identify new drug repositioning opportunities as well as ultimately be useful for personalized medicine by linking the genomic traits of patients to drug sensitivity.

Inhibitors that target the receptor tyrosine kinase (RTK)/Ras/mitogen-activated protein kinase (MAPK) pathway have led to clinical responses in lung and other cancers, but some patients fail to respond and in those that do resistance inevitably occurs (Balak et al., 2006; Kosaka et al., 2006; Rudin et al., 2013; Wagle et al., 2011). To understand intrinsic and acquired resistance to inhibition of MAPK signaling, we performed CRISPR-Cas9 gene deletion screens in the setting of BRAF, MEK, EGFR, and ALK inhibition. Loss of KEAP1, a negative regulator of NFE2L2/NRF2, modulated the response to BRAF, MEK, EGFR, and ALK inhibition in BRAF-, NRAS-, KRAS-, EGFR-, and ALK-mutant lung cancer cells. Treatment with inhibitors targeting the RTK/MAPK pathway increased reactive oxygen species (ROS) in cells with intact KEAP1, and loss of KEAP1 abrogated this increase. In addition, loss of KEAP1 altered cell metabolism to allow cells to proliferate in the absence of MAPK signaling. These observations suggest that alterations in the KEAP1/NRF2 pathway may promote survival in the presence of multiple inhibitors targeting the RTK/Ras/MAPK pathway.

Evolved resistance to tyrosine kinase inhibitor (TKI)-targeted therapies remains a major clinical challenge. In epidermal growth factor receptor (EGFR) mutant non-small-cell lung cancer (NSCLC), failure of EGFR TKIs can result from both genetic and epigenetic mechanisms of acquired drug resistance. Widespread reports of histologic and gene expression changes consistent with an epithelial-to-mesenchymal transition (EMT) have been associated with initially surviving drug-tolerant persister cells, which can seed bona fide genetic mechanisms of resistance to EGFR TKIs. While therapeutic approaches targeting fully resistant cells, such as those harboring an EGFRT790M mutation, have been developed, a clinical strategy for preventing the emergence of persister cells remains elusive. Using mesenchymal cell lines derived from biopsies of patients who progressed on EGFR TKI as surrogates for persister populations, we performed whole-genome CRISPR screening and identified fibroblast growth factor receptor 1 (FGFR1) as the top target promoting survival of mesenchymal EGFR mutant cancers. Although numerous previous reports of FGFR signaling contributing to EGFR TKI resistance in vitro exist, the data have not yet been sufficiently compelling to instigate a clinical trial testing this hypothesis, nor has the role of FGFR in promoting the survival of persister cells been elucidated. In this study, we find that combining EGFR and FGFR inhibitors inhibited the survival and expansion of EGFR mutant drug-tolerant cells over long time periods, preventing the development of fully resistant cancers in multiple vitro models and in vivo. These results suggest that dual EGFR and FGFR blockade may be a promising clinical strategy for both preventing and overcoming EMT-associated acquired drug resistance and provide motivation for the clinical study of combined EGFR and FGFR inhibition in EGFR-mutated NSCLCs.

SHP2 mediates RAS activation downstream of multiple receptor tyrosine kinases (RTKs) and cancer cell lines dependent on RTKs are in general dependent on SHP2. Profiling of the allosteric SHP2 inhibitor SHP099 across cancer cell lines harboring various RTK dependencies reveals that FGFR-dependent cells are often insensitive to SHP099 when compared to EGFR-dependent cells. We find that FGFR-driven cells depend on SHP2 but exhibit resistance to SHP2 inhibitors in vitro and in vivo. Treatment of such models with SHP2 inhibitors results in an initial decrease in phosphorylated ERK1/2 (p-ERK) levels, however p-ERK levels rapidly rebound within two hours. This p-ERK rebound is blocked by FGFR inhibitors or high doses of SHP2 inhibitors. Mechanistically, compared with EGFR-driven cells, FGFR-driven cells tend to express high levels of RTK negative regulators such as the SPRY family proteins, which are rapidly downregulated upon ERK inhibition. Moreover, over-expression of SPRY4 in FGFR-driven cells prevents MAPK pathway reactivation and sensitizes them to SHP2 inhibitors. We also identified two novel combination approaches to enhance the efficacy of SHP2 inhibitors, either with a distinct site 2 allosteric SHP2 inhibitor or with a RAS-SOS1 interaction inhibitor. Our findings suggest the rapid FGFR feedback activation following initial pathway inhibition by SHP2 inhibitors may promote the open conformation of SHP2 and lead to resistance to SHP2 inhibitors. These findings may assist to refine patient selection and predict resistance mechanisms in the clinical development of SHP2 inhibitors and to suggest strategies for discovering SHP2 inhibitors that are more effective against upstream feedback activation.

Phosphoinositide-3- kinase (PI3K) signaling regulates cellular proliferation, survival and metabolism, and its aberrant activation is one of the most frequent oncogenic events across human cancers. In the last few decades, research focused on the development of PI3K inhibitors, from preclinical tool compounds to the highly specific medicines approved to treat patients with cancer. Herein we discuss current paradigms for PI3K inhibitors in cancer therapy, focusing on clinical data and mechanisms of action. We also discuss current limitations in the use of PI3K inhibitors including toxicities and mechanisms of resistance, with specific emphasis on approaches aimed to improve their efficacy.

Combinations of anti-cancer drugs can overcome resistance and provide new treatments1,2. The number of possible drug combinations vastly exceeds what could be tested clinically. Efforts to systematically identify active combinations and the tissues and molecular contexts in which they are most effective could accelerate the development of combination treatments. Here we evaluate the potency and efficacy of 2,025 clinically relevant two-drug combinations, generating a dataset encompassing 125 molecularly characterized breast, colorectal and pancreatic cancer cell lines. We show that synergy between drugs is rare and highly context-dependent, and that combinations of targeted agents are most likely to be synergistic. We incorporate multi-omic molecular features to identify combination biomarkers and specify synergistic drug combinations and their active contexts, including in basal-like breast cancer, and microsatellite-stable or KRAS-mutant colon cancer. Our results show that irinotecan and CHEK1 inhibition have synergistic effects in microsatellite-stable or KRAS-TP53 double-mutant colon cancer cells, leading to apoptosis and suppression of tumour xenograft growth. This study identifies clinically relevant effective drug combinations in distinct molecular subpopulations and is a resource to guide rational efforts to develop combinatorial drug treatments.

Many gene fusions are reported in tumours and for most their role remains unknown. As fusions are used for diagnostic and prognostic purposes, and are targets for treatment, it is crucial to assess their function in cancer. To systematically investigate the role of fusions in tumour cell fitness, we utilized RNA-sequencing data from 1011 human cancer cell lines to functionally link 8354 fusion events with genomic data, sensitivity to >350 anti-cancer drugs and CRISPR-Cas9 loss-of-fitness effects. Established clinically-relevant fusions were identified. Overall, detection of functional fusions was rare, including those involving cancer driver genes, suggesting that many fusions are dispensable for tumour fitness. Therapeutically actionable fusions involving RAF1, BRD4 and ROS1 were verified in new histologies. In addition, recurrent YAP1-MAML2 fusions were identified as activators of Hippo-pathway signaling in multiple cancer types. Our approach discriminates functional fusions, identifying new drivers of carcinogenesis and fusions that could have clinical implications.

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.

Cancer-associated fibroblasts (CAFs) are highly heterogeneous. With the lack of a comprehensive understanding of CAFs' functional distinctions, it remains unclear how cancer treatments could be personalized based on CAFs in a patient's tumor. We have established a living biobank of CAFs derived from biopsies of patients' non-small lung cancer (NSCLC) that encompasses a broad molecular spectrum of CAFs in clinical NSCLC. By functionally interrogating CAF heterogeneity using the same therapeutics received by patients, we identify three functional subtypes: (1) robustly protective of cancers and highly expressing HGF and FGF7; (2) moderately protective of cancers and highly expressing FGF7; and (3) those providing minimal protection. These functional differences among CAFs are governed by their intrinsic TGF-β signaling, which suppresses HGF and FGF7 expression. This CAF functional classification correlates with patients' clinical response to targeted therapies and also associates with the tumor immune microenvironment, therefore providing an avenue to guide personalized treatment.

BH3 mimetics such as ABT-263 induce apoptosis in a subset of cancer models. However, these drugs have shown limited clinical efficacy as single agents in small-cell lung cancer (SCLC) and other solid tumor malignancies, and rational combination strategies remain underexplored. To develop a novel therapeutic approach, we examined the efficacy of ABT-263 across >500 cancer cell lines, including 311 for which we had matched expression data for select genes. We found that high expression of the proapoptotic gene Bcl2-interacting mediator of cell death (BIM) predicts sensitivity to ABT-263. In particular, SCLC cell lines possessed greater BIM transcript levels than most other solid tumors and are among the most sensitive to ABT-263. However, a subset of relatively resistant SCLC cell lines has concomitant high expression of the antiapoptotic myeloid cell leukemia 1 (MCL-1). Whereas ABT-263 released BIM from complexes with BCL-2 and BCL-XL, high expression of MCL-1 sequestered BIM released from BCL-2 and BCL-XL, thereby abrogating apoptosis. We found that SCLCs were sensitized to ABT-263 via TORC1/2 inhibition, which led to reduced MCL-1 protein levels, thereby facilitating BIM-mediated apoptosis. AZD8055 and ABT-263 together induced marked apoptosis in vitro, as well as tumor regressions in multiple SCLC xenograft models. In a Tp53; Rb1 deletion genetically engineered mouse model of SCLC, the combination of ABT-263 and AZD8055 significantly repressed tumor growth and induced tumor regressions compared with either drug alone. Furthermore, in a SCLC patient-derived xenograft model that was resistant to ABT-263 alone, the addition of AZD8055 induced potent tumor regression. Therefore, addition of a TORC1/2 inhibitor offers a therapeutic strategy to markedly improve ABT-263 activity in SCLC.

Fewer than half of children with high-risk neuroblastoma survive. Many of these tumors harbor high-level amplification of MYCN, which correlates with poor disease outcome. Using data from our large drug screen we predicted, and subsequently demonstrated, that MYCN-amplified neuroblastomas are sensitive to the BCL-2 inhibitor ABT-199. This sensitivity occurs in part through low anti-apoptotic BCL-xL expression, high pro-apoptotic NOXA expression, and paradoxical, MYCN-driven upregulation of NOXA. Screening for enhancers of ABT-199 sensitivity in MYCN-amplified neuroblastomas, we demonstrate that the Aurora Kinase A inhibitor MLN8237 combines with ABT-199 to induce widespread apoptosis. In diverse models of MYCN-amplified neuroblastoma, including a patient-derived xenograft model, this combination uniformly induced tumor shrinkage, and in multiple instances led to complete tumor regression.

In eukaryotic cells, the SH2 and PTB domains mediate protein-protein interactions by recognizing phosphotyrosine residues on target proteins. Here we make the unexpected finding that the C2 domain of PKCdelta directly binds to phosphotyrosine peptides in a sequence-specific manner. We provide evidence that this domain mediates PKCdelta interaction with a Src binding glycoprotein, CDCP1. The crystal structure of the PKCdelta C2 domain in complex with an optimal phosphopeptide reveals a new mode of phosphotyrosine binding in which the phosphotyrosine moiety forms a ring-stacking interaction with a histidine residue of the C2 domain. This is also the first example of a protein Ser/Thr kinase containing a domain that binds phosphotyrosine.

Recent clinical data have linked KRAS/TP53 comutation (mut) to resistance to radiotherapy (RT), but supporting laboratory in vivo evidence is lacking. In addition, the ability of different radiation doses, with/without epidermal growth factor receptor (EGFR)-directed treatment, to achieve local tumor control as a function of KRAS status is unknown. Here, we assessed clonogenic radiation survival of a panel of annotated lung cancer cell lines. KRASmut/TP53mut was associated with the highest radioresistance in nonisogenic and isogenic comparisons. To validate these findings, isogenic TP53mut NCI-H1703 models, KRASmut or wild-type (wt), were grown as heterotopic xenografts in nude mice. A clinical RT schedule of 30 fractions over 6 weeks was employed. The dose that controlled 50% of tumors (TCD50 ) was calculated. The TCD50 for KRASwt/TP53mut xenografts was 43.1 Gy whereas KRASmut/TP53mut tumors required a 1.9-fold higher TCD50 of 81.4 Gy. The EGFR inhibitor erlotinib radiosensitized KRASmut but not KRASwt cells and xenografts. The TCD50 associated with adding erlotinib to RT was 58.8 Gy for KRASmut, that is, a ~1.4-fold dose enhancement. However, the EGFR antibody cetuximab did not have a radiosensitizing effect. In conclusion, we demonstrate for the first time that KRASmut in a TP53mut background confers radioresistance when studying a clinical RT schedule and local control rather than tumor growth delay. Despite the known unresponsiveness of KRASmut tumors to EGFR inhibitors, erlotinib radiosensitized KRASmut tumors. Our data highlight KRAS/TP53 comutation as a candidate biomarker of radioresistance that can be at least partially reversed by dose escalation or the addition of a targeted agent.

α-Ketoglutarate (αKG) is a key node in many important metabolic pathways. The αKG analog N-oxalylglycine (NOG) and its cell-permeable prodrug dimethyloxalylglycine (DMOG) are extensively used to inhibit αKG-dependent dioxygenases. However, whether NOG interference with other αKG-dependent processes contributes to its mode of action remains poorly understood. Here we show that, in aqueous solutions, DMOG is rapidly hydrolyzed, yielding methyloxalylglycine (MOG). MOG elicits cytotoxicity in a manner that depends on its transport by monocarboxylate transporter 2 (MCT2) and is associated with decreased glutamine-derived tricarboxylic acid-cycle flux, suppressed mitochondrial respiration and decreased ATP production. MCT2-facilitated entry of MOG into cells leads to sufficiently high concentrations of NOG to inhibit multiple enzymes in glutamine metabolism, including glutamate dehydrogenase. These findings reveal that MCT2 dictates the mode of action of NOG by determining its intracellular concentration and have important implications for the use of (D)MOG in studying αKG-dependent signaling and metabolism.

Identifying robust, patient-specific, and predictive biomarkers presents a major obstacle in precision oncology. To optimize patient-specific therapeutic strategies, here we couple pathway knowledge with large-scale drug sensitivity, RNAi, and CRISPR-Cas9 screening data from 460 cell lines. Pathway activity levels are found to be strong predictive biomarkers for the essentiality of 15 proteins, including the essentiality of MAD2L1 in breast cancer patients with high BRCA-pathway activity. We also find strong predictive biomarkers for the sensitivity to 31 compounds, including BCL2 and microtubule inhibitors (MTIs). Lastly, we show that Bcl-xL inhibition can modulate the activity of a predictive biomarker pathway and re-sensitize lung cancer cells and tumors to MTI therapy. Overall, our results support the use of pathways in helping to achieve the goal of precision medicine by uncovering dozens of predictive biomarkers.

In patients with metastatic cancer, spatial heterogeneity of somatic alterations may lead to incomplete assessment of a cancer's mutational profile when analyzing a single tumor biopsy. In this study, we perform sequencing of cell-free DNA (cfDNA) and distinct metastatic tissue samples from ten rapid autopsy cases with pre-treated metastatic cancer. We show that levels of heterogeneity in genetic biomarkers vary between patients but that gene expression signatures representative of the tumor microenvironment are more consistent. Across nine patients with plasma samples available, we are able to detect 62/62 truncal and 47/121 non-truncal point mutations in cfDNA. We observe that mutation clonality in cfDNA is correlated with the number of metastatic lesions in which the mutation is detected and use this result to derive a clonality threshold to classify truncal and non-truncal driver alterations with reasonable specificity. In contrast, mutation truncality is more often incorrectly assigned when studying single tissue samples. Our results demonstrate the utility of a single cfDNA sample relative to that of single tissue samples when treating patients with metastatic cancer.