Massively parallel quantification of phenotypic heterogeneity in single-cell drug responses.

Single-cell analysis tools have made substantial advances in characterizing genomic heterogeneity; however, tools for measuring phenotypic heterogeneity have lagged due to the increased difficulty of handling live biology. Here, we report a single-cell phenotyping tool capable of measuring image-based clonal properties at scales approaching 100,000 clones per experiment. These advances are achieved by exploiting a previously unidentified flow regime in ladder microfluidic networks that, under appropriate conditions, yield a mathematically perfect cell trap. Machine learning and computer vision tools are used to control the imaging hardware and analyze the cellular phenotypic parameters within these images. Using this platform, we quantified the responses of tens of thousands of single cell­derived acute myeloid leukemia (AML) clones to targeted therapy, identifying rare resistance and morphological phenotypes at frequencies down to 0.05%. This approach can be extended to higher-level cellular architectures such as cell pairs and organoids and on-chip live-cell fluorescence assays.

FOXA1 and adaptive response determinants to HER2 targeted therapy in TBCRC 036.

Inhibition of the HER2/ERBB2 receptor is a keystone to treating HER2-positive malignancies, particularly breast cancer, but a significant fraction of HER2-positive (HER2+) breast cancers recur or fail to respond. Anti-HER2 monoclonal antibodies, like trastuzumab or pertuzumab, and ATP active site inhibitors like lapatinib, commonly lack durability because of adaptive changes in the tumor leading to resistance. HER2+ cell line responses to inhibition with lapatinib were analyzed by RNAseq and ChIPseq to characterize transcriptional and epigenetic changes. Motif analysis of lapatinib-responsive genomic regions implicated the pioneer transcription factor FOXA1 as a mediator of adaptive responses. Lapatinib in combination with FOXA1 depletion led to dysregulation of enhancers, impaired adaptive upregulation of HER3, and decreased proliferation. HER2-directed therapy using clinically relevant drugs (trastuzumab with or without lapatinib or pertuzumab) in a 7-day clinical trial designed to examine early pharmacodynamic response to antibody-based anti-HER2 therapy showed reduced FOXA1 expression was coincident with decreased HER2 and HER3 levels, decreased proliferation gene signatures, and increased immune gene signatures. This highlights the importance of the immune response to anti-HER2 antibodies and suggests that inhibiting FOXA1-mediated adaptive responses in combination with HER2 targeting is a potential therapeutic strategy.

Discrete Adaptive Responses to MEK Inhibitor in Subpopulations of Triple-Negative Breast Cancer.

Triple-negative breast cancers contain a spectrum of epithelial and mesenchymal phenotypes. SUM-229PE cells represent a model for this heterogeneity, maintaining both epithelial and mesenchymal subpopulations that are genomically similar but distinct in gene expression profiles. We identified differential regions of open chromatin in epithelial and mesenchymal cells that were strongly correlated with regions of H3K27ac. Motif analysis of these regions identified consensus sequences for transcription factors that regulate cell identity. Treatment with the MEK inhibitor trametinib induced enhancer remodeling that is associated with transcriptional regulation of genes in epithelial and mesenchymal cells. Motif analysis of enhancer peaks downregulated in response to chronic treatment with trametinib identified AP-1 motif enrichment in both epithelial and mesenchymal subpopulations. Chromatin immunoprecipitation sequencing (ChIP-seq) of JUNB identified subpopulation-specific localization, which was significantly enriched at regions of open chromatin. These results indicate that cell identity controls localization of transcription factors and chromatin-modifying enzymes to enhancers for differential control of gene expression. We identified increased H3K27ac at an enhancer region proximal to CXCR7, a G-protein-coupled receptor that increased 15-fold in expression in the epithelial subpopulation during chronic treatment. RNAi knockdown of CXCR7 inhibited proliferation in trametinib-resistant cells. Thus, adaptive resistance to chronic trametinib treatment contributes to proliferation in the presence of the drug. Acquired amplification of KRAS following trametinib dose escalation further contributed to POS cell proliferation. Adaptive followed by acquired gene expression changes contributed to proliferation in trametinib-resistant cells, suggesting inhibition of early transcriptional reprogramming could prevent resistance and the bypass of targeted therapy. IMPLICATIONS: We defined the differential responses to trametinib in subpopulations of a clinically relevant in vitro model of TNBC, and identified both adaptive and acquired elements that contribute to the emergence of drug resistance mediated by increased expression of CXCR7 and amplification of KRAS.

GSK2801, a BAZ2/BRD9 Bromodomain Inhibitor, Synergizes with BET Inhibitors to Induce Apoptosis in Triple-Negative Breast Cancer.

Screening of an inhibitor library targeting kinases and epigenetic regulators identified several molecules having antiproliferative synergy with extraterminal domain (BET) bromodomain (BD) inhibitors (JQ1, OTX015) in triple-negative breast cancer (TNBC). GSK2801, an inhibitor of BAZ2A/B BDs, of the imitation switch chromatin remodeling complexes, and BRD9, of the SWI/SNF complex, demonstrated synergy independent of BRD4 control of P-TEFb-mediated pause-release of RNA polymerase II. GSK2801 or RNAi knockdown of BAZ2A/B with JQ1 selectively displaced BRD2 at promoters/enhancers of ETS-regulated genes. Additional displacement of BRD2 from rDNA in the nucleolus coincided with decreased 45S rRNA, revealing a function of BRD2 in regulating RNA polymerase I transcription. In 2D cultures, enhanced displacement of BRD2 from chromatin by combination drug treatment induced senescence. In spheroid cultures, combination treatment induced cleaved caspase-3 and cleaved PARP characteristic of apoptosis in tumor cells. Thus, GSK2801 blocks BRD2-driven transcription in combination with BET inhibitor and induces apoptosis of TNBC. IMPLICATIONS: Synergistic inhibition of BDs encoded in BAZ2A/B, BRD9, and BET proteins induces apoptosis of TNBC by a combinatorial suppression of ribosomal DNA transcription and ETS-regulated genes.

Traditional and systems biology based drug discovery for the rare tumor syndrome neurofibromatosis type 2.

Neurofibromatosis 2 (NF2) is a rare tumor suppressor syndrome that manifests with multiple schwannomas and meningiomas. There are no effective drug therapies for these benign tumors and conventional therapies have limited efficacy. Various model systems have been created and several drug targets have been implicated in NF2-driven tumorigenesis based on known effects of the absence of merlin, the product of the NF2 gene. We tested priority compounds based on known biology with traditional dose-concentration studies in meningioma and schwann cell systems. Concurrently, we studied functional kinome and gene expression in these cells pre- and post-treatment to determine merlin deficient molecular phenotypes. Cell viability results showed that three agents (GSK2126458, Panobinostat, CUDC-907) had the greatest activity across schwannoma and meningioma cell systems, but merlin status did not significantly influence response. In vivo, drug effect was tumor specific with meningioma, but not schwannoma, showing response to GSK2126458 and Panobinostat. In culture, changes in both the transcriptome and kinome in response to treatment clustered predominantly based on tumor type. However, there were differences in both gene expression and functional kinome at baseline between meningioma and schwannoma cell systems that may form the basis for future selective therapies. This work has created an openly accessible resource (www.synapse.org/SynodosNF2) of fully characterized isogenic schwannoma and meningioma cell systems as well as a rich data source of kinome and transcriptome data from these assay systems before and after treatment that enables single and combination drug discovery based on molecular phenotype.

Proteomic analysis defines kinase taxonomies specific for subtypes of breast cancer.

Multiplexed small molecule inhibitors covalently bound to Sepharose beads (MIBs) were used to capture functional kinases in luminal, HER2-enriched and triple negative (basal-like and claudin-low) breast cancer cell lines and tumors. Kinase MIB-binding profiles at baseline without perturbation proteomically distinguished the four breast cancer subtypes. Understudied kinases, whose disease associations and pharmacology are generally unexplored, were highly represented in MIB-binding taxonomies and are integrated into signaling subnetworks with kinases that have been previously well characterized in breast cancer. Computationally it was possible to define subtypes using profiles of less than 50 of the more than 300 kinases bound to MIBs that included understudied as well as metabolic and lipid kinases. Furthermore, analysis of MIB-binding profiles established potential functional annotations for these understudied kinases. Thus, comprehensive MIBs-based capture of kinases provides a unique proteomics-based method for integration of poorly characterized kinases of the understudied kinome into functional subnetworks in breast cancer cells and tumors that is not possible using genomic strategies. The MIB-binding profiles readily defined subtype-selective differential adaptive kinome reprogramming in response to targeted kinase inhibition, demonstrating how MIB profiles can be used in determining dynamic kinome changes that result in subtype selective phenotypic state changes.

Epigenetic Mechanisms Regulating Adaptive Responses to Targeted Kinase Inhibitors in Cancer.

Although targeted inhibition of oncogenic kinase drivers has achieved remarkable patient responses in many cancers, the development of resistance has remained a significant challenge. Numerous mechanisms have been identified, including the acquisition of gatekeeper mutations, activating pathway mutations, and copy number loss or gain of the driver or alternate nodes. These changes have prompted the development of kinase inhibitors with increased selectivity, use of second-line therapeutics to overcome primary resistance, and combination treatment to forestall resistance. In addition to genomic resistance mechanisms, adaptive transcriptional and signaling responses seen in tumors are gaining appreciation as alterations that lead to a phenotypic state change-often observed as an epithelial-to-mesenchymal shift or reversion to a cancer stem cell-like phenotype underpinned by remodeling of the epigenetic landscape. This epigenomic modulation driving cell state change is multifaceted and includes modulation of repressive and activating histone modifications, DNA methylation, enhancer remodeling, and noncoding RNA species. Consequently, the combination of kinase inhibitors with drugs targeting components of the transcriptional machinery and histone-modifying enzymes has shown promise in preclinical and clinical studies. Here, we review mechanisms of resistance to kinase inhibition in cancer, with special emphasis on the rewired kinome and transcriptional signaling networks and the potential vulnerabilities that may be exploited to overcome these adaptive signaling changes.

Enhancer remodeling regulates epigenetic adaptation and resistance to MEK1/2 inhibition in triple-negative breast cancer.

Kinase inhibitors targeting the mitogen/extracellular signal-regulated kinase kinase (MEK)- extracellular signal related kinase (ERK) signaling pathway have limited durability in inhibiting growth of triple-negative breast cancer. We defined genome wide enhancer remodeling following MEK inhibition capable of driving adaptive gene transcription. Targeting positive elongation factor (P-TEFb) transcriptional regulatory complex members can block enhancer remodeling making the response to MEK-ERK inhibition durable.

Kinome and Transcriptome Profiling Reveal Broad and Distinct Activities of Erlotinib, Sunitinib, and Sorafenib in the Mouse Heart and Suggest Cardiotoxicity From Combined Signal Transducer and Activator of Transcription and Epidermal Growth Factor Receptor Inhibition.

BACKGROUND: Most novel cancer therapeutics target kinases that are essential to tumor survival. Some of these kinase inhibitors are associated with cardiotoxicity, whereas others appear to be cardiosafe. The basis for this distinction is unclear, as are the molecular effects of kinase inhibitors in the heart. METHODS AND RESULTS: We administered clinically relevant doses of sorafenib, sunitinib (cardiotoxic multitargeted kinase inhibitors), or erlotinib (a cardiosafe epidermal growth factor receptor inhibitor) to mice daily for 2 weeks. We then compared the effects of these 3 kinase inhibitors on the cardiac transcriptome using RNAseq and the cardiac kinome using multiplexed inhibitor beads coupled with mass spectrometry. We found unexpectedly broad molecular effects of all 3 kinase inhibitors, suggesting that target kinase selectivity does not define either the molecular response or the potential for cardiotoxicity. Using in vivo drug administration and primary cardiomyocyte culture, we also show that the cardiosafety of erlotinib treatment may result from upregulation of the cardioprotective signal transducer and activator of transcription 3 pathway, as co-treatment with erlotinib and a signal transducer and activator of transcription inhibitor decreases cardiac contractile function and cardiomyocyte fatty acid oxidation. CONCLUSIONS: Collectively our findings indicate that preclinical kinome and transcriptome profiling may predict the cardiotoxicity of novel kinase inhibitors, and suggest caution for the proposed therapeutic strategy of combined signal transducer and activator of transcription/epidermal growth factor receptor inhibition for cancer treatment.

MAP3K4 Controls the Chromatin Modifier HDAC6 during Trophoblast Stem Cell Epithelial-to-Mesenchymal Transition.

The first epithelial-to-mesenchymal transition (EMT) occurs in trophoblast stem (TS) cells during implantation. Inactivation of the serine/threonine kinase MAP3K4 in TS cells (TSKI4 cells) induces an intermediate state of EMT, where cells retain stemness, lose epithelial markers, and gain mesenchymal characteristics. Investigation of relationships among MAP3K4 activity, stemness, and EMT in TS cells may reveal key regulators of EMT. Here, we show that MAP3K4 activity controls EMT through the ubiquitination and degradation of HDAC6. Loss of MAP3K4 activity in TSKI4 cells results in elevated HDAC6 expression and the deacetylation of cytoplasmic and nuclear targets. In the nucleus, HDAC6 deacetylates the promoters of tight junction genes, promoting the dissolution of tight junctions. Importantly, HDAC6 knockdown in TSKI4 cells restores epithelial features, including cell-cell adhesion and barrier formation. These data define a role for HDAC6 in regulating gene expression during transitions between epithelial and mesenchymal phenotypes.

Enhancer Remodeling during Adaptive Bypass to MEK Inhibition Is Attenuated by Pharmacologic Targeting of the P-TEFb Complex.

Targeting the dysregulated BRAF-MEK-ERK pathway in cancer has increasingly emerged in clinical trial design. Despite clinical responses in specific cancers using inhibitors targeting BRAF and MEK, resistance develops often involving nongenomic adaptive bypass mechanisms. Inhibition of MEK1/2 by trametinib in patients with triple-negative breast cancer (TNBC) induced dramatic transcriptional responses, including upregulation of receptor tyrosine kinases (RTK) comparing tumor samples before and after one week of treatment. In preclinical models, MEK inhibition induced genome-wide enhancer formation involving the seeding of BRD4, MED1, H3K27 acetylation, and p300 that drives transcriptional adaptation. Inhibition of the P-TEFb-associated proteins BRD4 and CBP/p300 arrested enhancer seeding and RTK upregulation. BRD4 bromodomain inhibitors overcame trametinib resistance, producing sustained growth inhibition in cells, xenografts, and syngeneic mouse TNBC models. Pharmacologic targeting of P-TEFb members in conjunction with MEK inhibition by trametinib is an effective strategy to durably inhibit epigenomic remodeling required for adaptive resistance.Significance: Widespread transcriptional adaptation to pharmacologic MEK inhibition was observed in TNBC patient tumors. In preclinical models, MEK inhibition induces dramatic genome-wide modulation of chromatin, in the form of de novo enhancer formation and enhancer remodeling. Pharmacologic targeting of P-TEFb complex members at enhancers is an effective strategy to durably inhibit such adaptation. Cancer Discov; 7(3); 302-21. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 235.

Inhibition of Lapatinib-Induced Kinome Reprogramming in ERBB2-Positive Breast Cancer by Targeting BET Family Bromodomains.

Therapeutics that target ERBB2, such as lapatinib, often provide initial clinical benefit, but resistance frequently develops. Adaptive responses leading to lapatinib resistance involve reprogramming of the kinome through reactivation of ERBB2/ERBB3 signaling and transcriptional upregulation and activation of multiple tyrosine kinases. The heterogeneity of induced kinases prevents their targeting by a single kinase inhibitor, underscoring the challenge of predicting effective kinase inhibitor combination therapies. We hypothesized that, to make the tumor response to single kinase inhibitors durable, the adaptive kinome response itself must be inhibited. Genetic and chemical inhibition of BET bromodomain chromatin readers suppresses transcription of many lapatinib-induced kinases involved in resistance, including ERBB3, IGF1R, DDR1, MET, and FGFRs, preventing downstream SRC/FAK signaling and AKT reactivation. Combining inhibitors of kinases and chromatin readers prevents kinome adaptation by blocking transcription, generating a durable response to lapatinib, and overcoming the dilemma of heterogeneity in the adaptive response.

Assessing adaptation of the cancer kinome in response to targeted therapies.

Cancer cells are dependent on protein kinase signalling networks to drive proliferation and to promote survival, and, accordingly, kinases continue to represent a major target class for development of anti-cancer therapeutics. Kinase inhibitors nevertheless have yielded only limited success with many different malignancies due to the inability of single agents to sustain a durable clinical response. Cancer cell kinomes are highly resilient and able to bypass targeted kinase inhibition, leading to tumour resistance. A novel platform has been developed to analyse the activity of the expressed kinome using MIBs (multiplexed inhibitor beads), which consist of Sepharose beads with covalently immobilized inhibitors that preferentially bind activated kinases. Coupling MIB capture with MS (MIB-MS) allows simultaneous determination of the activity of over 75% of the expressed kinome, facilitating high-throughput assessment of adaptive kinase responses resulting from deregulated feedback and feedforward regulatory mechanisms. The adaptive response frequently involves transcriptional up-regulation of specific kinases that allow bypass of the targeted kinase. Understanding how the kinome reprogrammes to targeted kinase inhibition will allow novel therapeutic strategies to be developed for durable clinical responses.

Molecular pathways: adaptive kinome reprogramming in response to targeted inhibition of the BRAF-MEK-ERK pathway in cancer.

The central role of the BRAF-MEK-ERK pathway in controlling cell fate has made this pathway a primary target for deregulated activation in cancer. BRaf is activated by Ras proteins allowing Ras oncogenes to constitutively activate the pathway. Activating BRaf mutations are also frequent in several cancers, being the most common oncogenic mutation in thyroid carcinoma and melanoma. There are currently two inhibitors, vemurafenib and dabrafenib, approved for treatment of malignant melanoma having activating BRaf mutations. Concurrent administration of BRAF and MAP-ERK kinase (MEK) inhibitor (trametinib) is significantly more active in patients with BRAF-mutant melanoma than either single agent alone, but progression to resistance ultimately occurs by different mechanisms that increase the activation of extracellular signal-regulated kinase (ERK). Such adaptive changes in tumor cell signaling networks allow bypass of targeted oncoprotein inhibition. This is true with targeted inhibitors for BRaf and MEK as well as specific inhibitors for AKT, mTOR, and many receptor tyrosine kinases such as EGF receptor (EGFR) and HER2. It is this adaptive response to targeted kinase inhibitors that contributes to the failure of single-agent kinase inhibitors to have durable responses. This failure is seen in virtually all cancers treated with single-agent kinase inhibitors, most of which are not as dependent on a single signaling pathway such as BRaf-MEK-ERK in melanoma. Thus, understanding the breadth of adaptive reprogramming responses to specific targeted kinase inhibition will be critical to develop appropriate combination therapies for durable clinical responses.

A RhoC biosensor reveals differences in the activation kinetics of RhoA and RhoC in migrating cells.

RhoA and RhoC GTPases share 92% amino acid sequence identity, yet play different roles in regulating cell motility and morphology. To understand these differences, we developed and validated a biosensor of RhoC activation (RhoC FLARE). This was used together with a RhoA biosensor to compare the spatio-temporal dynamics of RhoA and RhoC activity during cell protrusion/retraction and macropinocytosis. Both GTPases were activated similarly at the cell edge, but in regions more distal from the edge RhoC showed higher activation during protrusion. The two isoforms differed markedly in the kinetics of activation. RhoC was activated concomitantly with RhoA at the cell edge, but distally, RhoC activation preceded RhoA activation, occurring before edge protrusion. During macropinocytosis, differences were observed during vesicle closure and in the area surrounding vesicle formation.

SWI/SNF chromatin-remodeling factor Smarcd3/Baf60c controls epithelial-mesenchymal transition by inducing Wnt5a signaling.

We previously identified a gene signature predicted to regulate the epithelial-mesenchymal transition (EMT) in both epithelial tissue stem cells and breast cancer cells. A phenotypic RNA interference (RNAi) screen identified the genes within this 140-gene signature that promoted the conversion of mesenchymal epithelial cell adhesion molecule-negative (EpCAM-) breast cancer cells to an epithelial EpCAM+/high phenotype. The screen identified 10 of the 140 genes whose individual knockdown was sufficient to promote EpCAM and E-cadherin expression. Among these 10 genes, RNAi silencing of the SWI/SNF chromatin-remodeling factor Smarcd3/Baf60c in EpCAM- breast cancer cells gave the most robust transition from the mesenchymal to epithelial phenotype. Conversely, expression of Smarcd3/Baf60c in immortalized human mammary epithelial cells induced an EMT. The mesenchymal-like phenotype promoted by Smarcd3/Baf60c expression resulted in gene expression changes in human mammary epithelial cells similar to that of claudin-low triple-negative breast cancer cells. These mammary epithelial cells expressing Smarcd3/Baf60c had upregulated Wnt5a expression. Inhibition of Wnt5a by either RNAi knockdown or blocking antibody reversed Smarcd3/Baf60c-induced EMT. Thus, Smarcd3/Baf60c epigenetically regulates EMT by activating WNT signaling pathways.

MicroRNA 9-3p targets ß1 integrin to sensitize claudin-low breast cancer cells to MEK inhibition.

MEK1/2 inhibitors such as AZD6244 are in clinical trials for the treatment of multiple cancers, including breast cancer. Targeted kinase inhibition can induce compensatory kinome changes, rendering single therapeutic agents ineffective. To identify target proteins to be used in a combinatorial approach to inhibit tumor cell growth, we used a novel strategy that identified microRNAs (miRNAs) that synergized with AZD6244 to inhibit the viability of the claudin-low breast cancer cell line MDA-MB-231. Screening of a miRNA mimic library revealed the ability of miR-9-3p to significantly enhance AZD6244-induced extracellular signal-regulated kinase inhibition and growth arrest, while miR-9-3p had little effect on growth alone. Promoter methylation of mir-9 genes correlated with low expression of miR-9-3p in different breast cancer cell lines. Consistent with miR-9-3p having synthetic enhancer tumor suppressor characteristics, miR-9-3p expression in combination with MEK inhibitor caused a sustained loss of c-MYC expression and growth inhibition. The ß1 integrin gene (ITGB1) was identified as a new miR-9-3p target, and the growth inhibition seen with small interfering RNA knockdown or antibody blocking of ITGB1 in combination with MEK inhibitor phenocopied the growth inhibition seen with miR-9-3p plus AZD6244. The miRNA screen led to identification of a druggable protein, ITGB1, whose functional inhibition synergizes with MEK inhibitor.

Dynamic reprogramming of the kinome in response to targeted MEK inhibition in triple-negative breast cancer.

Kinase inhibitors have limited success in cancer treatment because tumors circumvent their action. Using a quantitative proteomics approach, we assessed kinome activity in response to MEK inhibition in triple-negative breast cancer (TNBC) cells and genetically engineered mice (GEMMs). MEK inhibition caused acute ERK activity loss, resulting in rapid c-Myc degradation that induced expression and activation of several receptor tyrosine kinases (RTKs). RNAi knockdown of ERK or c-Myc mimicked RTK induction by MEK inhibitors, and prevention of proteasomal c-Myc degradation blocked kinome reprogramming. MEK inhibitor-induced RTK stimulation overcame MEK2 inhibition, but not MEK1 inhibition, reactivating ERK and producing drug resistance. The C3Tag GEMM for TNBC similarly induced RTKs in response to MEK inhibition. The inhibitor-induced RTK profile suggested a kinase inhibitor combination therapy that produced GEMM tumor apoptosis and regression where single agents were ineffective. This approach defines mechanisms of drug resistance, allowing rational design of combination therapies for cancer.