Cancer lineage-specific regulation of YAP responsive elements revealed through large-scale functional epigenomic screens.

YAP is a key transcriptional co-activator of TEADs, it regulates cell growth and is frequently activated in cancer. In Malignant Pleural Mesothelioma (MPM), YAP is activated by loss-of-function mutations in upstream components of the Hippo pathway, while, in Uveal Melanoma (UM), YAP is activated in a Hippo-independent manner. To date, it is unclear if and how the different oncogenic lesions activating YAP impact its oncogenic program, which is particularly relevant for designing selective anti-cancer therapies. Here we show that, despite YAP being essential in both MPM and UM, its interaction with TEAD is unexpectedly dispensable in UM, limiting the applicability of TEAD inhibitors in this cancer type. Systematic functional interrogation of YAP regulatory elements in both cancer types reveals convergent regulation of broad oncogenic drivers in both MPM and UM, but also strikingly selective programs. Our work reveals unanticipated lineage-specific features of the YAP regulatory network that provide important insights to guide the design of tailored therapeutic strategies to inhibit YAP signaling across different cancer types.

Integrated CRISPR screening and drug profiling identifies combination opportunities for EGFR, ALK, and BRAF/MEK inhibitors.

Anti-tumor efficacy of targeted therapies is variable across patients and cancer types. Even in patients with initial deep response, tumors are typically not eradicated and eventually relapse. To address these challenges, we present a systematic screen for targets that limit the anti-tumor efficacy of EGFR and ALK inhibitors in non-small cell lung cancer and BRAF/MEK inhibitors in colorectal cancer. Our approach includes genome-wide CRISPR screens with or without drugs targeting the oncogenic driver ("anchor therapy"), and large-scale pairwise combination screens of anchor therapies with 351 other drugs. Interestingly, targeting of a small number of genes, including MCL1, BCL2L1, and YAP1, sensitizes multiple cell lines to the respective anchor therapy. Data from drug combination screens with EGF816 and ceritinib indicate that dasatinib and agents disrupting microtubules act synergistically across many cell lines. Finally, we show that a higher-order-combination screen with 26 selected drugs in two resistant EGFR-mutant lung cancer cell lines identified active triplet combinations.

A highly multiplexed quantitative phosphosite assay for biology and preclinical studies.

Reliable methods to quantify dynamic signaling changes across diverse pathways are needed to better understand the effects of disease and drug treatment in cells and tissues but are presently lacking. Here, we present SigPath, a targeted mass spectrometry (MS) assay that measures 284 phosphosites in 200 phosphoproteins of biological interest. SigPath probes a broad swath of signaling biology with high throughput and quantitative precision. We applied the assay to investigate changes in phospho-signaling in drug-treated cancer cell lines, breast cancer preclinical models, and human medulloblastoma tumors. In addition to validating previous findings, SigPath detected and quantified a large number of differentially regulated phosphosites newly associated with disease models and human tumors at baseline or with drug perturbation. Our results highlight the potential of SigPath to monitor phosphoproteomic signaling events and to nominate mechanistic hypotheses regarding oncogenesis, response, and resistance to therapy.

Quantitative Proteomics of the Cancer Cell Line Encyclopedia.

Proteins are essential agents of biological processes. To date, large-scale profiling of cell line collections including the Cancer Cell Line Encyclopedia (CCLE) has focused primarily on genetic information whereas deep interrogation of the proteome has remained out of reach. Here, we expand the CCLE through quantitative profiling of thousands of proteins by mass spectrometry across 375 cell lines from diverse lineages to reveal information undiscovered by DNA and RNA methods. We observe unexpected correlations within and between pathways that are largely absent from RNA. An analysis of microsatellite instable (MSI) cell lines reveals the dysregulation of specific protein complexes associated with surveillance of mutation and translation. These and other protein complexes were associated with sensitivity to knockdown of several different genes. These data in conjunction with the wider CCLE are a broad resource to explore cellular behavior and facilitate cancer research.

Distinct Transcriptional Programming Drive Response to MAPK Inhibition in BRAF V600-Mutant Melanoma Patient-Derived Xenografts.

Inhibitors targeting BRAF and its downstream kinase MEK produce robust response in patients with advanced BRAF V600-mutant melanoma. However, the duration and depth of response vary significantly between patients; therefore, predicting response a priori remains a significant challenge. Here, we utilized the Novartis collection of patient-derived xenografts to characterize transcriptional alterations elicited by BRAF and MEK inhibitors in vivo, in an effort to identify mechanisms governing differential response to MAPK inhibition. We show that the expression of an MITF-high, "epithelial-like" transcriptional program is associated with reduced sensitivity and adaptive response to BRAF and MEK inhibitor treatment. On the other hand, xenograft models that express an MAPK-driven "mesenchymal-like" transcriptional program are preferentially sensitive to MAPK inhibition. These gene-expression programs are somewhat similar to the MITF-high and -low phenotypes described in cancer cell lines, but demonstrate an inverse relationship with drug response. This suggests a discrepancy between in vitro and in vivo experimental systems that warrants future investigations. Finally, BRAF V600-mutant melanoma relies on either MAPK or alternative pathways for survival under BRAF and MEK inhibition in vivo, which in turn predicts their response to further pathway suppression using a combination of BRAF, MEK, and ERK inhibitors. Our findings highlight the intertumor heterogeneity in BRAF V600-mutant melanoma, and the need for precision medicine strategies to target this aggressive cancer.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Tumor-derived IFN triggers chronic pathway agonism and sensitivity to ADAR loss.

Interferons (IFNs) are cytokines that play a critical role in limiting infectious and malignant diseases 1-4 . Emerging data suggest that the strength and duration of IFN signaling can differentially impact cancer therapies, including immune checkpoint blockade 5-7 . Here, we characterize the output of IFN signaling, specifically IFN-stimulated gene (ISG) signatures, in primary tumors from The Cancer Genome Atlas. While immune infiltration correlates with the ISG signature in some primary tumors, the existence of ISG signature-positive tumors without evident infiltration of IFN-producing immune cells suggests that cancer cells per se can be a source of IFN production. Consistent with this hypothesis, analysis of patient-derived tumor xenografts propagated in immune-deficient mice shows evidence of ISG-positive tumors that correlates with expression of human type I and III IFNs derived from the cancer cells. Mechanistic studies using cell line models from the Cancer Cell Line Encyclopedia that harbor ISG signatures demonstrate that this is a by-product of a STING-dependent pathway resulting in chronic tumor-derived IFN production. This imposes a transcriptional state on the tumor, poising it to respond to the aberrant accumulation of double-stranded RNA (dsRNA) due to increased sensor levels (MDA5, RIG-I and PKR). By interrogating our functional short-hairpin RNA screen dataset across 398 cancer cell lines, we show that this ISG transcriptional state creates a novel genetic vulnerability. ISG signature-positive cancer cells are sensitive to the loss of ADAR, a dsRNA-editing enzyme that is also an ISG. A genome-wide CRISPR genetic suppressor screen reveals that the entire type I IFN pathway and the dsRNA-activated kinase, PKR, are required for the lethality induced by ADAR depletion. Therefore, tumor-derived IFN resulting in chronic signaling creates a cellular state primed to respond to dsRNA accumulation, rendering ISG-positive tumors susceptible to ADAR loss.

TRPS1 Is a Lineage-Specific Transcriptional Dependency in Breast Cancer.

Perturbed epigenomic programs play key roles in tumorigenesis, and chromatin modulators are candidate therapeutic targets in various human cancer types. To define singular and shared dependencies on DNA and histone modifiers and transcription factors in poorly differentiated adult and pediatric cancers, we conducted a targeted shRNA screen across 59 cell lines of 6 cancer types. Here, we describe the TRPS1 transcription factor as a strong breast cancer-specific hit, owing largely to lineage-restricted expression. Knockdown of TRPS1 resulted in perturbed mitosis, apoptosis, and reduced tumor growth. Integrated analysis of TRPS1 transcriptional targets, chromatin binding, and protein interactions revealed that TRPS1 is associated with the NuRD repressor complex. These findings uncover a transcriptional network that is essential for breast cancer cell survival and propagation.

Correction of copy number induced false positives in CRISPR screens.

Cell autonomous cancer dependencies are now routinely identified using CRISPR loss-of-function viability screens. However, a bias exists that makes it difficult to assess the true essentiality of genes located in amplicons, since the entire amplified region can exhibit lethal scores. These false-positive hits can either be discarded from further analysis, which in cancer models can represent a significant number of hits, or methods can be developed to rescue the true-positives within amplified regions. We propose two methods to rescue true positive hits in amplified regions by correcting for this copy number artefact. The Local Drop Out (LDO) method uses the relative lethality scores within genomic regions to assess true essentiality and does not require additional orthogonal data (e.g. copy number value). LDO is meant to be used in screens covering a dense region of the genome (e.g. a whole chromosome or the whole genome). The General Additive Model (GAM) method models the screening data as a function of the known copy number values and removes the systematic effect from the measured lethality. GAM does not require the same density as LDO, but does require prior knowledge of the copy number values. Both methods have been developed with single sample experiments in mind so that the correction can be applied even in smaller screens. Here we demonstrate the efficacy of both methods at removing the copy number effect and rescuing hits from some of the amplified regions. We estimate a 70-80% decrease of false positive hits with either method in regions of high copy number compared to no correction.

Project DRIVE: A Compendium of Cancer Dependencies and Synthetic Lethal Relationships Uncovered by Large-Scale, Deep RNAi Screening.

Elucidation of the mutational landscape of human cancer has progressed rapidly and been accompanied by the development of therapeutics targeting mutant oncogenes. However, a comprehensive mapping of cancer dependencies has lagged behind and the discovery of therapeutic targets for counteracting tumor suppressor gene loss is needed. To identify vulnerabilities relevant to specific cancer subtypes, we conducted a large-scale RNAi screen in which viability effects of mRNA knockdown were assessed for 7,837 genes using an average of 20 shRNAs per gene in 398 cancer cell lines. We describe findings of this screen, outlining the classes of cancer dependency genes and their relationships to genetic, expression, and lineage features. In addition, we describe robust gene-interaction networks recapitulating both protein complexes and functional cooperation among complexes and pathways. This dataset along with a web portal is provided to the community to assist in the discovery and translation of new therapeutic approaches for cancer.

CRISPR Screens Provide a Comprehensive Assessment of Cancer Vulnerabilities but Generate False-Positive Hits for Highly Amplified Genomic Regions.

UNLABELLED: CRISPR/Cas9 has emerged as a powerful new tool to systematically probe gene function. We compared the performance of CRISPR to RNAi-based loss-of-function screens for the identification of cancer dependencies across multiple cancer cell lines. CRISPR dropout screens consistently identified more lethal genes than RNAi, implying that the identification of many cellular dependencies may require full gene inactivation. However, in two aneuploid cancer models, we found that all genes within highly amplified regions, including nonexpressed genes, scored as lethal by CRISPR, revealing an unanticipated class of false-positive hits. In addition, using a CRISPR tiling screen, we found that sgRNAs targeting essential domains generate the strongest lethality phenotypes and thus provide a strategy to rapidly define the protein domains required for cancer dependence. Collectively, these findings not only demonstrate the utility of CRISPR screens in the identification of cancer-essential genes, but also reveal the need to carefully control for false-positive results in chromosomally unstable cancer lines. SIGNIFICANCE: We show in this study that CRISPR-based screens have a significantly lower false-negative rate compared with RNAi-based screens, but have specific liabilities particularly in the interrogation of regions of genome amplification. Therefore, this study provides critical insights for applying CRISPR-based screens toward the systematic identification of new cancer targets. Cancer Discov; 6(8); 900-13. ©2016 AACR.See related commentary by Sheel and Xue, p. 824See related article by Aguirre et al., p. 914This article is highlighted in the In This Issue feature, p. 803.

Dynamic Regulation of α-Actinin's Calponin Homology Domains on F-Actin.

α-Actinin is an essential actin cross-linker involved in cytoskeletal organization and dynamics. The molecular conformation of α-actinin's actin-binding domain (ABD) regulates its association with actin and thus mutations in this domain can lead to severe pathogenic conditions. A point mutation at lysine 255 in human α-actinin-4 to glutamate increases the binding affinity resulting in stiffer cytoskeletal structures. The role of different ABD conformations and the effect of K255E mutation on ABD conformations remain elusive. To evaluate the impact of K255E mutation on ABD binding to actin we use all-atom molecular dynamics and free energy calculation methods and study the molecular mechanism of actin association in both wild-type α-actinin and in the K225E mutant. Our models illustrate that the strength of actin association is indeed sensitive to the ABD conformation, predict the effect of K255E mutation--based on simulations with the K237E mutant chicken α-actinin--and evaluate the mechanism of α-actinin binding to actin. Furthermore, our simulations showed that the calmodulin domain binding to the linker region was important for regulating the distance between actin and ABD. Our results provide valuable insights into the molecular details of this critical cellular phenomenon and further contribute to an understanding of cytoskeletal dynamics in health and disease.

Integrative analyses of proteomics and RNA transcriptomics implicate mitochondrial processes, protein folding pathways and GWAS loci in Parkinson disease.

BACKGROUND: Parkinson disease (PD) is a neurodegenerative disease characterized by the accumulation of alpha-synuclein (SNCA) and other proteins in aggregates termed "Lewy Bodies" within neurons. PD has both genetic and environmental risk factors, and while processes leading to aberrant protein aggregation are unknown, past work points to abnormal levels of SNCA and other proteins. Although several genome-wide studies have been performed for PD, these have focused on DNA sequence variants by genome-wide association studies (GWAS) and on RNA levels (microarray transcriptomics), while genome-wide proteomics analysis has been lacking. METHODS: This study employed two state-of-the-art technologies, three-stage Mass Spectrometry Tandem Mass Tag Proteomics (12 PD, 12 controls) and RNA-sequencing transcriptomics (29 PD, 44 controls), evaluated in the context of PD GWAS implicated loci and microarray transcriptomics (19 PD, 24 controls). The technologies applied for this study were performed in a set of overlapping prefrontal cortex (Brodmann area 9) samples obtained from PD patients and sex and age similar neurologically healthy controls. RESULTS: After appropriate filters, proteomics robustly identified 3558 unique proteins, with 283 of these (7.9 %) significantly different between PD and controls (q-value < 0.05). RNA-sequencing identified 17,580 protein-coding genes, with 1095 of these (6.2 %) significantly different (FDR p-value < 0.05); only 166 of the FDR significant protein-coding genes (0.94 %) were present among the 3558 proteins characterized. Of these 166, eight genes (4.8 %) were significant in both studies, with the same direction of effect. Functional enrichment analysis of the proteomics results strongly supports mitochondrial-related pathways, while comparable analysis of the RNA-sequencing results implicates protein folding pathways and metallothioneins. Ten of the implicated genes or proteins co-localized to GWAS loci. Evidence implicating SNCA was stronger in proteomics than in RNA-sequencing analyses. CONCLUSIONS: We report the largest analysis of proteomics in PD to date, and the first to combine this technology with RNA-sequencing to investigate GWAS implicated loci. Notably, differentially expressed protein-coding genes were more likely to not be characterized in the proteomics analysis, which lessens the ability to compare across platforms. Combining multiple genome-wide platforms offers novel insights into the pathological processes responsible for this disease by identifying pathways implicated across methodologies.

The talin dimer structure orientation is mechanically regulated.

Formation of a stable cell-substrate contact can be regulated by mechanical force, especially at the focal adhesion. Individual proteins that make up the focal adhesions, such as talin, can exhibit mechanosensing. We previously described one mode of talin mechanosensing in which the vinculin-binding site of talin is exposed after force-induced stretch of a single talin rod domain. Here, we describe a second mode of talin mechanosensing in which the talin dimer itself can adopt different orientations in response to mechanical stimulation. Using molecular dynamics models, we demonstrate that the C-terminus region of the talin dimer is flexible mainly at the linker between the dimerization helices and the nearby actin-binding helical bundle. Our molecular dynamics simulations reveal two possible orientations of the talin dimer at its C-terminus. The extracellular matrix (ECM)-bound integrins cross-linked by talin can be forced apart leading to an elongated orientation of the talin dimer, and the ECM-bound integrins can be forced together by the ECM producing a collapsed orientation of the talin dimer. Formation of the elongated orientation is shown to be more favorable. Switching between the two talin dimer orientations constitutes a mode of mechanosensing.

The interaction of vinculin with actin.

Vinculin can interact with F-actin both in recruitment of actin filaments to the growing focal adhesions and also in capping of actin filaments to regulate actin dynamics. Using molecular dynamics, both interactions are simulated using different vinculin conformations. Vinculin is simulated either with only its vinculin tail domain (Vt), with all residues in its closed conformation, with all residues in an open I conformation, and with all residues in an open II conformation. The open I conformation results from movement of domain 1 away from Vt; the open II conformation results from complete dissociation of Vt from the vinculin head domains. Simulation of vinculin binding along the actin filament showed that Vt alone can bind along the actin filaments, that vinculin in its closed conformation cannot bind along the actin filaments, and that vinculin in its open I conformation can bind along the actin filaments. The simulations confirm that movement of domain 1 away from Vt in formation of vinculin 1 is sufficient for allowing Vt to bind along the actin filament. Simulation of Vt capping actin filaments probe six possible bound structures and suggest that vinculin would cap actin filaments by interacting with both S1 and S3 of the barbed-end, using the surface of Vt normally occluded by D4 and nearby vinculin head domain residues. Simulation of D4 separation from Vt after D1 separation formed the open II conformation. Binding of open II vinculin to the barbed-end suggests this conformation allows for vinculin capping. Three binding sites on F-actin are suggested as regions that could link to vinculin. Vinculin is suggested to function as a variable switch at the focal adhesions. The conformation of vinculin and the precise F-actin binding conformation is dependent on the level of mechanical load on the focal adhesion.

A molecular trajectory of α-actinin activation.

The mechanisms by which living cells respond to mechanical stimuli are not yet fully understood. It has been suggested that mechanosensing proteins play an important role in mechanotransduction because their binding affinities are directly affected by the external stress. α-Actinin is an actin cross-linker and may act as a mechanosensor in adhesion sites. Its interaction with vinculin is suggested to be mechanically regulated. In this study, the free energy of activation is explored using the umbrella sampling method. An activation trajectory is generated in which α-actinin's vinculin-binding site swings out of the rod domain, leading to approximately an 8 kcal/mol free energy release. The activation trajectory reveals several local and global conformational changes along the activation pathway accompanied by the breakage of a number of key interactions stabilizing the inhibited structure. These results may shed light on the role of α-actinin in cellular mechanotransduction and focal adhesion formation.

Phosphorylation primes vinculin for activation.

Vinculin phosphorylation has been implicated as a potential mechanism for focal adhesion growth and maturation. Four vinculin residues-Y100, S1033, S1045, and Y1065-are phosphorylated by kinases during focal adhesion maturation. In this study, phosphorylation at each of these residues is simulated using molecular dynamics models. The simulations demonstrate that once each phosphorylated vinculin structure is at equilibrium, significant local conformational changes result that may impact either vinculin activation or vinculin binding to actin and PIP2. Simulation of vinculin activation after phosphorylation shows that the added phosphoryl groups can prime vinculin for activation. It remains to be seen if vinculin can be phosphorylated at S1033 in vivo, but these simulations highlight that in the event of a S1033 phophorylation vinculin will likely be primed for activation.

Vinculin activation is necessary for complete talin binding.

Focal adhesions are critical to a number of cellular processes that involve mechanotransduction and mechanical interaction with the cellular environment. The growth and strengthening of these focal adhesions is dependent on the interaction between talin and vinculin. This study investigates said interaction and how vinculin activation influences it. Using molecular dynamics, the interaction between talin's vinculin binding site (VBS) and vinculin's domain 1 (D1) is simulated both before and after vinculin activation. The simulations of VBS binding to vinculin before activation suggest the proximity of the vinculin tail to D1 prevents helical movement in D1 and thus prevents binding of VBS. In contrast, interaction of VBS with activated vinculin shows the possibility of complete VBS insertion into D1. In the simulations of both activated and autoinhibited vinculin where VBS fails to fully bind, VBS demonstrates significant hydrophobic interaction with surface residues in D1. These interactions link VBS to D1 even without its proper insertion into the hydrophobic core. Together these simulations suggest VBS binds to vinculin with the following mechanism: 1), VBS links to D1 via surface hydrophobic interactions; 2), vinculin undergoes activation and D1 is moved away from the vinculin tail; 3), helices in D1 undergo conformational change to allow VBS binding; and 4), VBS inserts itself into the hydrophobic core of D1.

A molecular dynamics investigation of vinculin activation.

Vinculin activation plays a critical role in focal adhesion initiation and formation. In its native state, vinculin is in an autoinhibitory conformation in which domain 1 prevents interaction of the vinculin tail domain with actin by steric hindrance. Once activated, vinculin is able to interact with both actin and talin. Several hypotheses have been put forth addressing the mechanisms of vinculin activation. One set of studies suggests that vinculin interaction with talin is sufficient to cause activation, whereas another set of studies suggests that a simultaneous interaction with several binding partners is necessary to achieve vinculin activation. Using molecular-dynamics (MD) simulations, we investigate the mechanisms of vinculin activation and suggest both a trajectory of conformational changes leading to vinculin activation, and key structural features that are likely involved in stabilizing the autoinhibited conformation. Assuming that the simultaneous interaction of vinculin with both actin and talin causes a stretching force on vinculin, and that vinculin activation results from a removal of steric hindrance blocking the actin-binding sites, we simulate with MD the stretching and activation of vinculin. The MD simulations are further confirmed by normal-mode analysis and simulation after residue modification. Taken together, the results of these simulations suggest that bending of the vinculin-binding-site region in vinculin away from the vinculin tail is the likely trajectory of vinculin activation.

Molecular mechanics of the alpha-actinin rod domain: bending, torsional, and extensional behavior.

alpha-Actinin is an actin crosslinking molecule that can serve as a scaffold and maintain dynamic actin filament networks. As a crosslinker in the stressed cytoskeleton, alpha-actinin can retain conformation, function, and strength. alpha-Actinin has an actin binding domain and a calmodulin homology domain separated by a long rod domain. Using molecular dynamics and normal mode analysis, we suggest that the alpha-actinin rod domain has flexible terminal regions which can twist and extend under mechanical stress, yet has a highly rigid interior region stabilized by aromatic packing within each spectrin repeat, by electrostatic interactions between the spectrin repeats, and by strong salt bridges between its two anti-parallel monomers. By exploring the natural vibrations of the alpha-actinin rod domain and by conducting bending molecular dynamics simulations we also predict that bending of the rod domain is possible with minimal force. We introduce computational methods for analyzing the torsional strain of molecules using rotating constraints. Molecular dynamics extension of the alpha-actinin rod is also performed, demonstrating transduction of the unfolding forces across salt bridges to the associated monomer of the alpha-actinin rod domain.