Identification of drugs including a dopamine receptor antagonist that selectively target cancer stem cells.

Selective targeting of cancer stem cells (CSCs) offers promise for a new generation of therapeutics. However, assays for both human CSCs and normal stem cells that are amenable to robust biological screens are limited. Using a discovery platform that reveals differences between neoplastic and normal human pluripotent stem cells (hPSC), we identify small molecules from libraries of known compounds that induce differentiation to overcome neoplastic self-renewal. Surprisingly, thioridazine, an antipsychotic drug, selectively targets the neoplastic cells, and impairs human somatic CSCs capable of in vivo leukemic disease initiation while having no effect on normal blood SCs. The drug antagonizes dopamine receptors that are expressed on CSCs and on breast cancer cells as well. These results suggest that dopamine receptors may serve as a biomarker for diverse malignancies, demonstrate the utility of using neoplastic hPSCs for identifying CSC-targeting drugs, and provide support for the use of differentiation as a therapeutic strategy.

RPL22 is a tumor suppressor in MSI-high cancers and a splicing regulator of MDM4.

Microsatellite instability-high (MSI-H) tumors are malignant tumors that, despite harboring a high mutational burden, often have intact TP53. One of the most frequent mutations in MSI-H tumors is a frameshift mutation in RPL22, a ribosomal protein. Here, we identified RPL22 as a modulator of MDM4 splicing through an alternative splicing switch in exon 6. RPL22 loss increases MDM4 exon 6 inclusion and cell proliferation and augments resistance to the MDM inhibitor Nutlin-3a. RPL22 represses the expression of its paralog, RPL22L1, by mediating the splicing of a cryptic exon corresponding to a truncated transcript. Therefore, damaging mutations in RPL22 drive oncogenic MDM4 induction and reveal a common splicing circuit in MSI-H tumors that may inform therapeutic targeting of the MDM4-p53 axis and oncogenic RPL22L1 induction.

Comprehensive structure-function analysis reveals gain- and loss-of-function mechanisms impacting oncogenic KRAS activity.

To dissect variant-function relationships in the KRAS oncoprotein, we performed deep mutational scanning (DMS) screens for both wild-type and KRAS G12D mutant alleles. We defined the spectrum of oncogenic potential for nearly all possible KRAS variants, identifying several novel transforming alleles and elucidating a model to describe the frequency of KRAS mutations in human cancer as a function of transforming potential, mutational probability, and tissue-specific mutational signatures. Biochemical and structural analyses of variants identified in a KRAS G12D second-site suppressor DMS screen revealed that attenuation of oncogenic KRAS can be mediated by protein instability and conformational rigidity, resulting in reduced binding affinity to effector proteins, such as RAF and PI3-kinases, or reduced SOS-mediated nucleotide exchange activity. These studies define the landscape of single amino acid alterations that modulate the function of KRAS, providing a resource for the clinical interpretation of KRAS variants and elucidating mechanisms of oncogenic KRAS inactivation for therapeutic exploitation.

RPL22 is a tumor suppressor in MSI-high cancers and a key splicing regulator of MDM4.

Microsatellite instability high (MSI-H) tumors are malignant tumors that, despite harboring a high mutational burden, often have intact TP53. One of the most frequent mutations in MSI-H tumors is a frameshift mutation in RPL22, a ribosomal protein. Here, we identified RPL22 as a modulator of MDM4 splicing through an alternative splicing switch in exon 6. RPL22 loss increases MDM4 exon 6 inclusion, cell proliferation, and augments resistance to the MDM inhibitor Nutlin-3a. RPL22 represses expression of its paralog, RPL22L1, by mediating the splicing of a cryptic exon corresponding to a truncated transcript. Therefore, damaging mutations in RPL22 drive oncogenic MDM4 induction and reveal a common splicing circuit in MSI-H tumors that may inform therapeutic targeting of the MDM4-p53 axis and oncogenic RPL22L1 induction.

Massively parallel phenotyping of coding variants in cancer with Perturb-seq.

Genome sequencing studies have identified millions of somatic variants in cancer, but it remains challenging to predict the phenotypic impact of most. Experimental approaches to distinguish impactful variants often use phenotypic assays that report on predefined gene-specific functional effects in bulk cell populations. Here, we develop an approach to functionally assess variant impact in single cells by pooled Perturb-seq. We measured the impact of 200 TP53 and KRAS variants on RNA profiles in over 300,000 single lung cancer cells, and used the profiles to categorize variants into phenotypic subsets to distinguish gain-of-function, loss-of-function and dominant negative variants, which we validated by comparison with orthogonal assays. We discovered that KRAS variants did not merely fit into discrete functional categories, but spanned a continuum of gain-of-function phenotypes, and that their functional impact could not have been predicted solely by their frequency in patient cohorts. Our work provides a scalable, gene-agnostic method for coding variant impact phenotyping, with potential applications in multiple disease settings.

Allosteric inhibition of PPM1D serine/threonine phosphatase via an altered conformational state.

PPM1D encodes a serine/threonine phosphatase that regulates numerous pathways including the DNA damage response and p53. Activating mutations and amplification of PPM1D are found across numerous cancer types. GSK2830371 is a potent and selective allosteric inhibitor of PPM1D, but its mechanism of binding and inhibition of catalytic activity are unknown. Here we use computational, biochemical and functional genetic studies to elucidate the molecular basis of GSK2830371 activity. These data confirm that GSK2830371 binds an allosteric site of PPM1D with high affinity. By further incorporating data from hydrogen deuterium exchange mass spectrometry and sedimentation velocity analytical ultracentrifugation, we demonstrate that PPM1D exists in an equilibrium between two conformations that are defined by the movement of the flap domain, which is required for substrate recognition. A hinge region was identified that is critical for switching between the two conformations and was directly implicated in the high-affinity binding of GSK2830371 to PPM1D. We propose that the two conformations represent active and inactive forms of the protein reflected by the position of the flap, and that binding of GSK2830371 shifts the equilibrium to the inactive form. Finally, we found that C-terminal truncating mutations proximal to residue 400 result in destabilization of the protein via loss of a stabilizing N- and C-terminal interaction, consistent with the observation from human genetic data that nearly all PPM1D mutations in cancer are truncating and occur distal to residue 400. Taken together, our findings elucidate the mechanism by which binding of a small molecule to an allosteric site of PPM1D inhibits its activity and provides insights into the biology of PPM1D.

SMAD4 represses FOSL1 expression and pancreatic cancer metastatic colonization.

Metastasis is a complex and poorly understood process. In pancreatic cancer, loss of the transforming growth factor (TGF)-ß/BMP effector SMAD4 is correlated with changes in altered histopathological transitions, metastatic disease, and poor prognosis. In this study, we use isogenic cancer cell lines to identify SMAD4 regulated genes that contribute to the development of metastatic colonization. We perform an in vivo screen identifying FOSL1 as both a SMAD4 target and sufficient to drive colonization to the lung. The targeting of these genes early in treatment may provide a therapeutic benefit.

The Pea3 Ets transcription factor regulates differentiation of multipotent progenitor cells during mammary gland development.

The Pea3 Ets transcription factor is overexpressed in breast tumors suggesting that it plays a role in mammary oncogenesis. However, the normal biological function of Pea3 in the mammary gland is not known. Here we report that Pea3 was expressed in the epithelium of the mouse mammary anlagen commensurate with their genesis, and at later times in the nipple and mammary ducts of female embryos. In adult mice Pea3 transcripts peaked at the onset of puberty and early pregnancy, times of active epithelial cell proliferation and differentiation. Pea3 was expressed in all progenitor cap cells and rare body cells of terminal end buds, and in the myoepithelial cells of ducts and alveoli. Analyses of the mammary glands of Pea3-null mice during puberty revealed an increased number of terminal end buds and an increased fraction of proliferating progenitor cells within these structures compared to their wild type littermates. Tissue transplant experiments demonstrated that these phenotypes were intrinsic to the Pea3-null mammary epithelium. During pregnancy, mammary glands isolated from Pea3-null females had impaired alveolar development as revealed by a decreased fraction of alveolar structures. We performed in vitro colony forming assays of mammary epithelial cells and discovered that loss of Pea3 altered the distribution of specific multipotent progenitor cells. Double-immunofluorescence confirmed that multipotential progenitors co-expressing markers of the myoepithelial and luminal epithelial lineage were amplified in the mammary glands of Pea3-null mice by comparison to their wild type counterparts. We propose that Pea3 functions in multipotential progenitors to regulate their lineage-specific differentiation potential.

Rhabdoid Tumors Are Sensitive to the Protein-Translation Inhibitor Homoharringtonine.

PURPOSE: Rhabdoid tumors are devastating pediatric cancers in need of improved therapies. We sought to identify small molecules that exhibit in vitro and in vivo efficacy against preclinical models of rhabdoid tumor. EXPERIMENTAL DESIGN: We screened eight rhabdoid tumor cell lines with 481 small molecules and compared their sensitivity with that of 879 other cancer cell lines. Genome-scale CRISPR-Cas9 inactivation screens in rhabdoid tumors were analyzed to confirm target vulnerabilities. Gene expression and CRISPR-Cas9 data were queried across cell lines and primary rhabdoid tumors to discover biomarkers of small-molecule sensitivity. Molecular correlates were validated by manipulating gene expression. Subcutaneous rhabdoid tumor xenografts were treated with the most effective drug to confirm in vitro results. RESULTS: Small-molecule screening identified the protein-translation inhibitor homoharringtonine (HHT), an FDA-approved treatment for chronic myelogenous leukemia (CML), as the sole drug to which all rhabdoid tumor cell lines were selectively sensitive. Validation studies confirmed the sensitivity of rhabdoid tumor to HHT was comparable with that of CML cell lines. Low expression of the antiapoptotic gene BCL2L1, which encodes Bcl-XL, was the strongest predictor of HHT sensitivity, and HHT treatment consistently depleted Mcl-1, the synthetic-lethal antiapoptotic partner of Bcl-XL. Rhabdoid tumor cell lines and primary-tumor samples expressed low BCL2L1, and overexpression of BCL2L1 induced resistance to HHT in rhabdoid tumor cells. Furthermore, HHT treatment inhibited rhabdoid tumor cell line and patient-derived xenograft growth in vivo. CONCLUSIONS: Rhabdoid tumor cell lines and xenografts are highly sensitive to HHT, at least partially due to their low expression of BCL2L1. HHT may have therapeutic potential against rhabdoid tumors.

Selective USP7 inhibition elicits cancer cell killing through a p53-dependent mechanism.

Ubiquitin specific peptidase 7 (USP7) is a deubiquitinating enzyme (DUB) that removes ubiquitin tags from specific protein substrates in order to alter their degradation rate and sub-cellular localization. USP7 has been proposed as a therapeutic target in several cancers because it has many reported substrates with a role in cancer progression, including FOXO4, MDM2, N-Myc, and PTEN. The multi-substrate nature of USP7, combined with the modest potency and selectivity of early generation USP7 inhibitors, has presented a challenge in defining predictors of response to USP7 and potential patient populations that would benefit most from USP7-targeted drugs. Here, we describe the structure-guided development of XL177A, which irreversibly inhibits USP7 with sub-nM potency and selectivity across the human proteome. Evaluation of the cellular effects of XL177A reveals that selective USP7 inhibition suppresses cancer cell growth predominantly through a p53-dependent mechanism: XL177A specifically upregulates p53 transcriptional targets transcriptome-wide, hotspot mutations in TP53 but not any other genes predict response to XL177A across a panel of ~500 cancer cell lines, and TP53 knockout rescues XL177A-mediated growth suppression of TP53 wild-type (WT) cells. Together, these findings suggest TP53 mutational status as a biomarker for response to USP7 inhibition. We find that Ewing sarcoma and malignant rhabdoid tumor (MRT), two pediatric cancers that are sensitive to other p53-dependent cytotoxic drugs, also display increased sensitivity to XL177A.

MDM2 and MDM4 Are Therapeutic Vulnerabilities in Malignant Rhabdoid Tumors.

Malignant rhabdoid tumors (MRT) are highly aggressive pediatric cancers that respond poorly to current therapies. In this study, we screened several MRT cell lines with large-scale RNAi, CRISPR-Cas9, and small-molecule libraries to identify potential drug targets specific for these cancers. We discovered MDM2 and MDM4, the canonical negative regulators of p53, as significant vulnerabilities. Using two compounds currently in clinical development, idasanutlin (MDM2-specific) and ATSP-7041 (MDM2/4-dual), we show that MRT cells were more sensitive than other p53 wild-type cancer cell lines to inhibition of MDM2 alone as well as dual inhibition of MDM2/4. These compounds caused significant upregulation of the p53 pathway in MRT cells, and sensitivity was ablated by CRISPR-Cas9-mediated inactivation of TP53. We show that loss of SMARCB1, a subunit of the SWI/SNF (BAF) complex mutated in nearly all MRTs, sensitized cells to MDM2 and MDM2/4 inhibition by enhancing p53-mediated apoptosis. Both MDM2 and MDM2/4 inhibition slowed MRT xenograft growth in vivo, with a 5-day idasanutlin pulse causing marked regression of all xenografts, including durable complete responses in 50% of mice. Together, these studies identify a genetic connection between mutations in the SWI/SNF chromatin-remodeling complex and the tumor suppressor gene TP53 and provide preclinical evidence to support the targeting of MDM2 and MDM4 in this often-fatal pediatric cancer. SIGNIFICANCE: This study identifies two targets, MDM2 and MDM4, as vulnerabilities in a deadly pediatric cancer and provides preclinical evidence that compounds inhibiting these proteins have therapeutic potential.

A dominant-negative effect drives selection of TP53 missense mutations in myeloid malignancies.

TP53, which encodes the tumor suppressor p53, is the most frequently mutated gene in human cancer. The selective pressures shaping its mutational spectrum, dominated by missense mutations, are enigmatic, and neomorphic gain-of-function (GOF) activities have been implicated. We used CRISPR-Cas9 to generate isogenic human leukemia cell lines of the most common TP53 missense mutations. Functional, DNA-binding, and transcriptional analyses revealed loss of function but no GOF effects. Comprehensive mutational scanning of p53 single-amino acid variants demonstrated that missense variants in the DNA-binding domain exert a dominant-negative effect (DNE). In mice, the DNE of p53 missense variants confers a selective advantage to hematopoietic cells on DNA damage. Analysis of clinical outcomes in patients with acute myeloid leukemia showed no evidence of GOF for TP53 missense mutations. Thus, a DNE is the primary unit of selection for TP53 missense mutations in myeloid malignancies.

CREB5 Promotes Resistance to Androgen-Receptor Antagonists and Androgen Deprivation in Prostate Cancer.

Androgen-receptor (AR) inhibitors, including enzalutamide, are used for treatment of all metastatic castration-resistant prostate cancers (mCRPCs). However, some patients develop resistance or never respond. We find that the transcription factor CREB5 confers enzalutamide resistance in an open reading frame (ORF) expression screen and in tumor xenografts. CREB5 overexpression is essential for an enzalutamide-resistant patient-derived organoid. In AR-expressing prostate cancer cells, CREB5 interactions enhance AR activity at a subset of promoters and enhancers upon enzalutamide treatment, including MYC and genes involved in the cell cycle. In mCRPC, we found recurrent amplification and overexpression of CREB5. Our observations identify CREB5 as one mechanism that drives resistance to AR antagonists in prostate cancers.

An alternative splicing switch in FLNB promotes the mesenchymal cell state in human breast cancer.

Alternative splicing of mRNA precursors represents a key gene expression regulatory step and permits the generation of distinct protein products with diverse functions. In a genome-scale expression screen for inducers of the epithelial-to-mesenchymal transition (EMT), we found a striking enrichment of RNA-binding proteins. We validated that QKI and RBFOX1 were necessary and sufficient to induce an intermediate mesenchymal cell state and increased tumorigenicity. Using RNA-seq and eCLIP analysis, we found that QKI and RBFOX1 coordinately regulated the splicing and function of the actin-binding protein FLNB, which plays a causal role in the regulation of EMT. Specifically, the skipping of FLNB exon 30 induced EMT by releasing the FOXC1 transcription factor. Moreover, skipping of FLNB exon 30 is strongly associated with EMT gene signatures in basal-like breast cancer patient samples. These observations identify a specific dysregulation of splicing, which regulates tumor cell plasticity and is frequently observed in human cancer.

Mutational processes shape the landscape of TP53 mutations in human cancer.

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations1,2. Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models3-8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

PRMT1-Mediated Translation Regulation Is a Crucial Vulnerability of Cancer.

Through an shRNA screen, we identified the protein arginine methyltransferase Prmt1 as a vulnerable intervention point in murine p53/Rb-null osteosarcomas, the human counterpart of which lacks effective therapeutic options. Depletion of Prmt1 in p53-deficient cells impaired tumor initiation and maintenance in vitro and in vivo Mechanistic studies reveal that translation-associated pathways were enriched for Prmt1 downstream targets, implicating Prmt1 in translation control. In particular, loss of Prmt1 led to a decrease in arginine methylation of the translation initiation complex, thereby disrupting its assembly and inhibiting translation. p53/Rb-null cells were sensitive to p53-induced translation stress, and analysis of human cancer cell line data from Project Achilles further revealed that Prmt1 and translation-associated pathways converged on the same functional networks. We propose that targeted therapy against Prmt1 and its associated translation-related pathways offer a mechanistic rationale for treatment of osteosarcomas and other cancers that exhibit dependencies on translation stress response. Cancer Res; 77(17); 4613-25. ©2017 AACR.

Serotonin transporter antagonists target tumor-initiating cells in a transgenic mouse model of breast cancer.

Accumulating data suggests that the initiation and progression of human breast tumors is fueled by a rare subpopulation of tumor cells, termed breast tumor-initiating cells (BTIC), which resist radiotherapy and chemotherapy. Consequently, therapies that abrogate BTIC activity are needed to achieve durable cures for breast cancer patients. To identify such therapies we used a sensitive assay to complete a high-throughput screen of small molecules, including approved drugs, with BTIC-rich mouse mammary tumor cell populations. We found that inhibitors of the serotonin reuptake transporter (SERT) and serotonin receptors, which include approved drugs used to treat mood disorders, were potent inhibitors of mouse BTIC activity as determined by functional sphere-forming assays and the initiation of tumor formation by transplant of drug-exposed tumor cells into syngeneic mice. Moreover, sertraline (Zoloft), a selective serotonin reuptake inhibitor (SSRI), synergized with docetaxel (Taxotere) to shrink mouse breast tumors in vivo. Hence drugs targeting the serotonergic system might be repurposed to treat breast cancer patients to afford more durable breast cancer remissions.

Integrative radiogenomic profiling of squamous cell lung cancer.

Radiotherapy is one of the mainstays of anticancer treatment, but the relationship between the radiosensitivity of cancer cells and their genomic characteristics is still not well defined. Here, we report the development of a high-throughput platform for measuring radiation survival in vitro and its validation in comparison with conventional clonogenic radiation survival analysis. We combined results from this high-throughput assay with genomic parameters in cell lines from squamous cell lung carcinoma, which is standardly treated by radiotherapy, to identify parameters that predict radiation sensitivity. We showed that activation of NFE2L2, a frequent event in lung squamous cancers, confers radiation resistance. An expression-based, in silico screen nominated inhibitors of phosphoinositide 3-kinase (PI3K) as NFE2L2 antagonists. We showed that the selective PI3K inhibitor, NVP-BKM120, both decreased NRF2 protein levels and sensitized NFE2L2 or KEAP1-mutant cells to radiation. We then combined results from this high-throughput assay with single-sample gene set enrichment analysis of gene expression data. The resulting analysis identified pathways implicated in cell survival, genotoxic stress, detoxification, and innate and adaptive immunity as key correlates of radiation sensitivity. The integrative and high-throughput methods shown here for large-scale profiling of radiation survival and genomic features of solid-tumor-derived cell lines should facilitate tumor radiogenomics and the discovery of genotype-selective radiation sensitizers and protective agents.

Systematic interrogation of 3q26 identifies TLOC1 and SKIL as cancer drivers.

UNLABELLED: 3q26 is frequently amplified in several cancer types with a common amplified region containing 20 genes. To identify cancer driver genes in this region, we interrogated the function of each of these genes by loss- and gain-of-function genetic screens. Specifically, we found that TLOC1 (SEC62) was selectively required for the proliferation of cell lines with 3q26 amplification. Increased TLOC1 expression induced anchorage-independent growth, and a second 3q26 gene, SKIL (SNON), facilitated cell invasion in immortalized human mammary epithelial cells. Expression of both TLOC1 and SKIL induced subcutaneous tumor growth. Proteomic studies showed that TLOC1 binds to DDX3X, which is essential for TLOC1-induced transformation and affected protein translation. SKIL induced invasion through upregulation of SLUG (SNAI2) expression. Together, these studies identify TLOC1 and SKIL as driver genes at 3q26 and more broadly suggest that cooperating genes may be coamplified in other regions with somatic copy number gain. SIGNIFICANCE: These studies identify TLOC1 and SKIL as driver genes in 3q26. These observations provide evidence that regions of somatic copy number gain may harbor cooperating genes of different but complementary functions.