DNA-encoded library-enabled discovery of proximity-inducing small molecules.

Small molecules that induce protein-protein associations represent powerful tools to modulate cell circuitry. We sought to develop a platform for the direct discovery of compounds able to induce association of any two preselected proteins, using the E3 ligase von Hippel-Lindau (VHL) and bromodomains as test systems. Leveraging the screening power of DNA-encoded libraries (DELs), we synthesized ~1 million DNA-encoded compounds that possess a VHL-targeting ligand, a variety of connectors and a diversity element generated by split-and-pool combinatorial chemistry. By screening our DEL against bromodomains in the presence and absence of VHL, we could identify VHL-bound molecules that simultaneously bind bromodomains. For highly barcode-enriched library members, ternary complex formation leading to bromodomain degradation was confirmed in cells. Furthermore, a ternary complex crystal structure was obtained for our most enriched library member with BRD4BD1 and a VHL complex. Our work provides a foundation for adapting DEL screening to the discovery of proximity-inducing small molecules.

Chemogenomic library design strategies for precision oncology, applied to phenotypic profiling of glioblastoma patient cells.

Designing a targeted screening library of bioactive small molecules is a challenging task since most compounds modulate their effects through multiple protein targets with varying degrees of potency and selectivity. We implemented analytic procedures for designing anticancer compound libraries adjusted for library size, cellular activity, chemical diversity and availability, and target selectivity. The resulting compound collections cover a wide range of protein targets and biological pathways implicated in various cancers, making them widely applicable to precision oncology. We characterized the compound and target spaces of the virtual libraries, in comparison with a minimal screening library of 1,211 compounds for targeting 1,386 anticancer proteins. In a pilot screening study, we identified patient-specific vulnerabilities by imaging glioma stem cells from patients with glioblastoma (GBM), using a physical library of 789 compounds that cover 1,320 of the anticancer targets. The cell survival profiling revealed highly heterogeneous phenotypic responses across the patients and GBM subtypes.

Mapping the landscape of genetic dependencies in chordoma.

Identifying the spectrum of genes required for cancer cell survival can reveal essential cancer circuitry and therapeutic targets, but such a map remains incomplete for many cancer types. We apply genome-scale CRISPR-Cas9 loss-of-function screens to map the landscape of selectively essential genes in chordoma, a bone cancer with few validated targets. This approach confirms a known chordoma dependency, TBXT (T; brachyury), and identifies a range of additional dependencies, including PTPN11, ADAR, PRKRA, LUC7L2, SRRM2, SLC2A1, SLC7A5, FANCM, and THAP1. CDK6, SOX9, and EGFR, genes previously implicated in chordoma biology, are also recovered. We find genomic and transcriptomic features that predict specific dependencies, including interferon-stimulated gene expression, which correlates with ADAR dependence and is elevated in chordoma. Validating the therapeutic relevance of dependencies, small-molecule inhibitors of SHP2, encoded by PTPN11, have potent preclinical efficacy against chordoma. Our results generate an emerging map of chordoma dependencies to enable biological and therapeutic hypotheses.

An expanded universe of cancer targets.

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

Targeted brachyury degradation disrupts a highly specific autoregulatory program controlling chordoma cell identity.

Chordomas are rare spinal tumors addicted to expression of the developmental transcription factor brachyury. In chordomas, brachyury is super-enhancer associated and preferentially downregulated by pharmacologic transcriptional CDK inhibition, leading to cell death. To understand the underlying basis of this sensitivity, we dissect the brachyury transcription regulatory network and compare the consequences of brachyury degradation with transcriptional CDK inhibition. Brachyury defines the chordoma super-enhancer landscape and autoregulates through binding its super-enhancer, and its locus forms a transcriptional condensate. Transcriptional CDK inhibition and brachyury degradation disrupt brachyury autoregulation, leading to loss of its transcriptional condensate and transcriptional program. Compared with transcriptional CDK inhibition, which globally downregulates transcription, leading to cell death, brachyury degradation is much more selective, inducing senescence and sensitizing cells to anti-apoptotic inhibition. These data suggest that brachyury downregulation is a core tenet of transcriptional CDK inhibition and motivates developing strategies to target brachyury and its autoregulatory feedback loop.

Plasticity of ether lipids promotes ferroptosis susceptibility and evasion.

Ferroptosis-an iron-dependent, non-apoptotic cell death process-is involved in various degenerative diseases and represents a targetable susceptibility in certain cancers1. The ferroptosis-susceptible cell state can either pre-exist in cells that arise from certain lineages or be acquired during cell-state transitions2-5. However, precisely how susceptibility to ferroptosis is dynamically regulated remains poorly understood. Here we use genome-wide CRISPR-Cas9 suppressor screens to identify the oxidative organelles peroxisomes as critical contributors to ferroptosis sensitivity in human renal and ovarian carcinoma cells. Using lipidomic profiling we show that peroxisomes contribute to ferroptosis by synthesizing polyunsaturated ether phospholipids (PUFA-ePLs), which act as substrates for lipid peroxidation that, in turn, results in the induction of ferroptosis. Carcinoma cells that are initially sensitive to ferroptosis can switch to a ferroptosis-resistant state in vivo in mice, which is associated with extensive downregulation of PUFA-ePLs. We further find that the pro-ferroptotic role of PUFA-ePLs can be extended beyond neoplastic cells to other cell types, including neurons and cardiomyocytes. Together, our work reveals roles for the peroxisome-ether-phospholipid axis in driving susceptibility to and evasion from ferroptosis, highlights PUFA-ePL as a distinct functional lipid class that is dynamically regulated during cell-state transitions, and suggests multiple regulatory nodes for therapeutic interventions in diseases that involve ferroptosis.

Selective covalent targeting of GPX4 using masked nitrile-oxide electrophiles.

We recently described glutathione peroxidase 4 (GPX4) as a promising target for killing therapy-resistant cancer cells via ferroptosis. The onset of therapy resistance by multiple types of treatment results in a stable cell state marked by high levels of polyunsaturated lipids and an acquired dependency on GPX4. Unfortunately, all existing inhibitors of GPX4 act covalently via a reactive alkyl chloride moiety that confers poor selectivity and pharmacokinetic properties. Here, we report our discovery that masked nitrile-oxide electrophiles, which have not been explored previously as covalent cellular probes, undergo remarkable chemical transformations in cells and provide an effective strategy for selective targeting of GPX4. The new GPX4-inhibiting compounds we describe exhibit unexpected proteome-wide selectivity and, in some instances, vastly improved physiochemical and pharmacokinetic properties compared to existing chloroacetamide-based GPX4 inhibitors. These features make them superior tool compounds for biological interrogation of ferroptosis and constitute starting points for development of improved inhibitors of GPX4.

Small-Molecule and CRISPR Screening Converge to Reveal Receptor Tyrosine Kinase Dependencies in Pediatric Rhabdoid Tumors.

Cancer is often seen as a disease of mutations and chromosomal abnormalities. However, some cancers, including pediatric rhabdoid tumors (RTs), lack recurrent alterations targetable by current drugs and need alternative, informed therapeutic options. To nominate potential targets, we performed a high-throughput small-molecule screen complemented by a genome-scale CRISPR-Cas9 gene-knockout screen in a large number of RT and control cell lines. These approaches converged to reveal several receptor tyrosine kinases (RTKs) as therapeutic targets, with RTK inhibition effective in suppressing RT cell growth in vitro and against a xenograft model in vivo. RT cell lines highly express and activate (phosphorylate) different RTKs, creating dependency without mutation or amplification. Downstream of RTK signaling, we identified PTPN11, encoding the pro-growth signaling protein SHP2, as a shared dependency across all RT cell lines. This study demonstrates that large-scale perturbational screening can uncover vulnerabilities in cancers with "quiet" genomes.

A GPX4-dependent cancer cell state underlies the clear-cell morphology and confers sensitivity to ferroptosis.

Clear-cell carcinomas (CCCs) are a histological group of highly aggressive malignancies commonly originating in the kidney and ovary. CCCs are distinguished by aberrant lipid and glycogen accumulation and are refractory to a broad range of anti-cancer therapies. Here we identify an intrinsic vulnerability to ferroptosis associated with the unique metabolic state in CCCs. This vulnerability transcends lineage and genetic landscape, and can be exploited by inhibiting glutathione peroxidase 4 (GPX4) with small-molecules. Using CRISPR screening and lipidomic profiling, we identify the hypoxia-inducible factor (HIF) pathway as a driver of this vulnerability. In renal CCCs, HIF-2α selectively enriches polyunsaturated lipids, the rate-limiting substrates for lipid peroxidation, by activating the expression of hypoxia-inducible, lipid droplet-associated protein (HILPDA). Our study suggests targeting GPX4 as a therapeutic opportunity in CCCs, and highlights that therapeutic approaches can be identified on the basis of cell states manifested by morphological and metabolic features in hard-to-treat cancers.

Renal medullary carcinomas depend upon SMARCB1 loss and are sensitive to proteasome inhibition.

Renal medullary carcinoma (RMC) is a rare and deadly kidney cancer in patients of African descent with sickle cell trait. We have developed faithful patient-derived RMC models and using whole-genome sequencing, we identified loss-of-function intronic fusion events in one SMARCB1 allele with concurrent loss of the other allele. Biochemical and functional characterization of these models revealed that RMC requires the loss of SMARCB1 for survival. Through integration of RNAi and CRISPR-Cas9 loss-of-function genetic screens and a small-molecule screen, we found that the ubiquitin-proteasome system (UPS) was essential in RMC. Inhibition of the UPS caused a G2/M arrest due to constitutive accumulation of cyclin B1. These observations extend across cancers that harbor SMARCB1 loss, which also require expression of the E2 ubiquitin-conjugating enzyme, UBE2C. Our studies identify a synthetic lethal relationship between SMARCB1-deficient cancers and reliance on the UPS which provides the foundation for a mechanism-informed clinical trial with proteasome inhibitors.

Stabilization of the Max Homodimer with a Small Molecule Attenuates Myc-Driven Transcription.

The transcription factor Max is a basic-helix-loop-helix leucine zipper (bHLHLZ) protein that forms homodimers or interacts with other bHLHLZ proteins, including Myc and Mxd proteins. Among this dynamic network of interactions, the Myc/Max heterodimer has crucial roles in regulating normal cellular processes, but its transcriptional activity is deregulated in a majority of human cancers. Despite this significance, the arsenal of high-quality chemical probes to interrogate these proteins remains limited. We used small molecule microarrays to identify compounds that bind Max in a mechanistically unbiased manner. We discovered the asymmetric polycyclic lactam, KI-MS2-008, which stabilizes the Max homodimer while reducing Myc protein and Myc-regulated transcript levels. KI-MS2-008 also decreases viable cancer cell growth in a Myc-dependent manner and suppresses tumor growth in vivo. This approach demonstrates the feasibility of modulating Max with small molecules and supports altering Max dimerization as an alternative approach to targeting Myc.

Small-molecule targeting of brachyury transcription factor addiction in chordoma.

Chordoma is a primary bone cancer with no approved therapy1. The identification of therapeutic targets in this disease has been challenging due to the infrequent occurrence of clinically actionable somatic mutations in chordoma tumors2,3. Here we describe the discovery of therapeutically targetable chordoma dependencies via genome-scale CRISPR-Cas9 screening and focused small-molecule sensitivity profiling. These systematic approaches reveal that the developmental transcription factor T (brachyury; TBXT) is the top selectively essential gene in chordoma, and that transcriptional cyclin-dependent kinase (CDK) inhibitors targeting CDK7/12/13 and CDK9 potently suppress chordoma cell proliferation. In other cancer types, transcriptional CDK inhibitors have been observed to downregulate highly expressed, enhancer-associated oncogenic transcription factors4,5. In chordoma, we find that T is associated with a 1.5-Mb region containing 'super-enhancers' and is the most highly expressed super-enhancer-associated transcription factor. Notably, transcriptional CDK inhibition leads to preferential and concentration-dependent downregulation of cellular brachyury protein levels in all models tested. In vivo, CDK7/12/13-inhibitor treatment substantially reduces tumor growth. Together, these data demonstrate small-molecule targeting of brachyury transcription factor addiction in chordoma, identify a mechanism of T gene regulation that underlies this therapeutic strategy, and provide a blueprint for applying systematic genetic and chemical screening approaches to discover vulnerabilities in genomically quiet cancers.

Computational Analyses Connect Small-Molecule Sensitivity to Cellular Features Using Large Panels of Cancer Cell Lines.

We recently pioneered several analyses of small-molecule sensitivity data collected from large-scale perturbation of hundreds of cancer cell lines with hundreds of small molecules, with cell viability measured as a readout of compound sensitivity. We performed these studies using cancer cell lines previously annotated with cellular, genomic, and basal gene-expression features. By combining small-molecule sensitivity data with these other datasets, we identified new candidate biomarkers of sensitivity, gained insights into small-molecule mechanisms of action, and proposed candidate hypotheses for cancer dependencies (including candidate combination therapies). Nevertheless, given the size of these datasets, we expect that many connections between cellular features and small-molecule sensitivity remain underexplored. In this chapter, we provide a step-by-step account of foundational data-analysis methods underlying our published studies, including working MATLAB code applied to our own public datasets. These procedures will allow others to repeat analyses of our data with new parameters, in additional contexts, and to adapt our procedures to their own datasets.

Synergistic Effects of Stereochemistry and Appendages on the Performance Diversity of a Collection of Synthetic Compounds.

Target- and phenotype-agnostic assessments of biological activity have emerged as viable strategies for prioritizing scaffolds, structural features, and synthetic pathways in screening sets, with the goal of increasing performance diversity. Here, we describe the synthesis of a small library of functionalized stereoisomeric azetidines and its biological annotation by "cell painting," a multiplexed, high-content imaging assay capable of measuring many hundreds of compound-induced changes in cell morphology in a quantitative and unbiased fashion. Using this approach, we systematically compare the degrees to which a core scaffold's biological activity, inferred from its effects on cell morphology, is affected by variations in stereochemistry and appendages. We show that stereoisomerism and appendage diversification can produce effects of similar magnitude, and that the concurrent use of these strategies results in a broader sampling of biological activity.

Modeling the impact of drug interactions on therapeutic selectivity.

Combination therapies that produce synergistic growth inhibition are widely sought after for the pharmacotherapy of many pathological conditions. Therapeutic selectivity, however, depends on the difference between potency on disease-causing cells and potency on non-target cell types that cause toxic side effects. Here, we examine a model system of antimicrobial compound combinations applied to two highly diverged yeast species. We find that even though the drug interactions correlate between the two species, cell-type-specific differences in drug interactions are common and can dramatically alter the selectivity of compounds when applied in combination vs. single-drug activity-enhancing, diminishing, or inverting therapeutic windows. This study identifies drug combinations with enhanced cell-type-selectivity with a range of interaction types, which we experimentally validate using multiplexed drug-interaction assays for heterogeneous cell cultures. This analysis presents a model framework for evaluating drug combinations with increased efficacy and selectivity against pathogens or tumors.

A precision oncology approach to the pharmacological targeting of mechanistic dependencies in neuroendocrine tumors.

We introduce and validate a new precision oncology framework for the systematic prioritization of drugs targeting mechanistic tumor dependencies in individual patients. Compounds are prioritized on the basis of their ability to invert the concerted activity of master regulator proteins that mechanistically regulate tumor cell state, as assessed from systematic drug perturbation assays. We validated the approach on a cohort of 212 gastroenteropancreatic neuroendocrine tumors (GEP-NETs), a rare malignancy originating in the pancreas and gastrointestinal tract. The analysis identified several master regulator proteins, including key regulators of neuroendocrine lineage progenitor state and immunoevasion, whose role as critical tumor dependencies was experimentally confirmed. Transcriptome analysis of GEP-NET-derived cells, perturbed with a library of 107 compounds, identified the HDAC class I inhibitor entinostat as a potent inhibitor of master regulator activity for 42% of metastatic GEP-NET patients, abrogating tumor growth in vivo. This approach may thus complement current efforts in precision oncology.

RWEN: response-weighted elastic net for prediction of chemosensitivity of cancer cell lines.

Motivation: In recent years there have been several efforts to generate sensitivity profiles of collections of genomically characterized cell lines to panels of candidate therapeutic compounds. These data provide the basis for the development of in silico models of sensitivity based on cellular, genetic, or expression biomarkers of cancer cells. However, a remaining challenge is an efficient way to identify accurate sets of biomarkers to validate. To address this challenge, we developed methodology using gene-expression profiles of human cancer cell lines to predict the responses of these cell lines to a panel of compounds. Results: We developed an iterative weighting scheme which, when applied to elastic net, a regularized regression method, significantly improves the overall accuracy of predictions, particularly in the highly sensitive response region. In addition to application of these methods to actual chemical sensitivity data, we investigated the effects of sample size, number of features, model sparsity, signal-to-noise ratio, and feature correlation on predictive performance using a simulation framework, particularly for situations where the number of covariates is much larger than sample size. While our method aims to be useful in therapeutic discovery and understanding of the basic mechanisms of action of drugs and their targets, it is generally applicable in any domain where predictions of extreme responses are of highest importance. Availability and implementation: The iterative and other weighting algorithms were implemented in R. The code is available at https://github.com/kiwtir/RWEN. The CTRP data are available at ftp://caftpd.nci.nih.gov/pub/OCG-DCC/CTD2/Broad/CTRPv2.1_2016_pub_NatChemBiol_12_109/ and the Sanger data at ftp://ftp.sanger.ac.uk/pub/project/cancerrxgene/releases/release-6.0/. Supplementary information: Supplementary data are available at Bioinformatics online.

A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles.

We previously piloted the concept of a Connectivity Map (CMap), whereby genes, drugs, and disease states are connected by virtue of common gene-expression signatures. Here, we report more than a 1,000-fold scale-up of the CMap as part of the NIH LINCS Consortium, made possible by a new, low-cost, high-throughput reduced representation expression profiling method that we term L1000. We show that L1000 is highly reproducible, comparable to RNA sequencing, and suitable for computational inference of the expression levels of 81% of non-measured transcripts. We further show that the expanded CMap can be used to discover mechanism of action of small molecules, functionally annotate genetic variants of disease genes, and inform clinical trials. The 1.3 million L1000 profiles described here, as well as tools for their analysis, are available at https://clue.io.

CTD2 Dashboard: a searchable web interface to connect validated results from the Cancer Target Discovery and Development Network.

Database URL: https://ctd2-dashboard.nci.nih.gov/.

Decomposing Oncogenic Transcriptional Signatures to Generate Maps of Divergent Cellular States.

The systematic sequencing of the cancer genome has led to the identification of numerous genetic alterations in cancer. However, a deeper understanding of the functional consequences of these alterations is necessary to guide appropriate therapeutic strategies. Here, we describe Onco-GPS (OncoGenic Positioning System), a data-driven analysis framework to organize individual tumor samples with shared oncogenic alterations onto a reference map defined by their underlying cellular states. We applied the methodology to the RAS pathway and identified nine distinct components that reflect transcriptional activities downstream of RAS and defined several functional states associated with patterns of transcriptional component activation that associates with genomic hallmarks and response to genetic and pharmacological perturbations. These results show that the Onco-GPS is an effective approach to explore the complex landscape of oncogenic cellular states across cancers, and an analytic framework to summarize knowledge, establish relationships, and generate more effective disease models for research or as part of individualized precision medicine paradigms.