The intrinsic substrate specificity of the human tyrosine kinome.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Dissecting signaling regulators driving AXL-mediated bypass resistance and associated phenotypes by phosphosite perturbations.

Receptor tyrosine kinase (RTK)-targeted therapies are often effective but invariably limited by drug resistance. A major mechanism of acquired resistance involves "bypass" switching to alternative pathways driven by non-targeted RTKs that restore proliferation. One such RTK is AXL whose overexpression, frequently observed in bypass resistant tumors, drives both cell survival and associated malignant phenotypes such as epithelial-to-mesenchymal (EMT) transition and migration. However, the signaling molecules and pathways eliciting these responses have remained elusive. To explore these coordinated effects, we generated a panel of mutant lung adenocarcinoma PC9 cell lines in which each AXL intracellular tyrosine residue was mutated to phenylalanine. By integrating measurements of phosphorylation signaling and other phenotypic changes associated with resistance through multivariate modeling, we mapped signaling perturbations to specific resistant phenotypes. Our results suggest that AXL signaling can be summarized into two clusters associated with progressive disease and poor clinical outcomes in lung cancer patients. These clusters displayed favorable Abl1 and SFK motifs and their phosphorylation was consistently decreased by dasatinib. High-throughput kinase specificity profiling showed that AXL likely activates the SFK cluster through FAK1 which is known to complex with Src. Moreover, the SFK cluster overlapped with a previously established focal adhesion kinase (FAK1) signature conferring EMT-mediated erlotinib resistance in lung cancer cells. Finally, we show that downstream of this kinase signaling, AXL and YAP form a positive feedback loop that sustains drug tolerant persister cells. Altogether, this work demonstrates an approach for dissecting signaling regulators by which AXL drives erlotinib resistance-associated phenotypic changes.

The Fanconi anemia core complex promotes CtIP-dependent end-resection to drive homologous recombination at DNA double-strand breaks.

Homologous Recombination (HR) is a high-fidelity repair mechanism of DNA Double-Strand Breaks (DSBs), which are induced by irradiation, genotoxic chemicals or physiological DNA damaging processes. DSBs are also generated as intermediates during the repair of interstrand crosslinks (ICLs). In this context, the Fanconi anemia (FA) core complex, which is effectively recruited to ICLs, promotes HR-mediated DSB-repair. However, whether the FA core complex also promotes HR at ICL-independent DSBs remains controversial. Here, we identified the FA core complex members FANCL and Ube2T as HR-promoting factors in a CRISPR/Cas9-based screen with cells carrying the DSB-repair reporter DSB-Spectrum. Using isogenic cell-line models, we validated the HR-function of FANCL and Ube2T, and demonstrated a similar function for their ubiquitination-substrate FANCD2. We further show that FANCL and Ube2T are directly recruited to DSBs and are required for the accumulation of FANCD2 at these break sites. Mechanistically, we demonstrate that FANCL ubiquitin ligase activity is required for the accumulation of the nuclease CtIP at DSBs, and consequently for optimal end-resection and Rad51 loading. CtIP overexpression rescues HR in FANCL-deficient cells, validating that FANCL primarily regulates HR by promoting CtIP recruitment. Together, these data demonstrate that the FA core complex and FANCD2 have a dual genome maintenance function by promoting repair of DSBs as well as the repair of ICLs.

Global view of the RAF-MEK-ERK module and its immediate downstream effectors.

Small molecule inhibitors of BRAF and MEK have proven effective at inhibiting tumor growth in melanoma patients, however this efficacy is limited due to the almost universal development of drug resistance. To provide advanced insight into the signaling responses that occur following kinase inhibition we have performed quantitative (phospho)-proteomics of human melanoma cells treated with either dabrafenib, a BRAF inhibitor; trametinib, a MEK inhibitor or SCH772984, an ERK inhibitor. Over nine experiments we identified 7827 class I phosphorylation sites on 4960 proteins. This included 54 phosphorylation sites that were significantly down-modulated after exposure to all three inhibitors, 34 of which have not been previously reported. Functional analysis of these novel ERK targets identified roles for them in GTPase activity and regulation, apoptosis and cell-cell adhesion. Comparison of the results presented here with previously reported phosphorylation sites downstream of ERK showed a limited degree of overlap suggesting that ERK signaling responses may be highly cell line and cue specific. In addition we identified 26 phosphorylation sites that were only responsive to dabrafenib. We provide further orthogonal experimental evidence for 3 of these sites in human embryonic kidney cells over-expressing BRAF as well as further computational insights using KinomeXplorer. The validated phosphorylation sites were found to be involved in actin regulation, which has been proposed as a novel mechanism for inhibiting resistance development. These results would suggest that the linearity of the BRAF-MEK-ERK module is at least context dependent.

Comprehensive substrate specificity profiling of the human Nek kinome reveals unexpected signaling outputs.

Human NimA-related kinases (Neks) have multiple mitotic and non-mitotic functions, but few substrates are known. We systematically determined the phosphorylation-site motifs for the entire Nek kinase family, except for Nek11. While all Nek kinases strongly select for hydrophobic residues in the -3 position, the family separates into four distinct groups based on specificity for a serine versus threonine phospho-acceptor, and preference for basic or acidic residues in other positions. Unlike Nek1-Nek9, Nek10 is a dual-specificity kinase that efficiently phosphorylates itself and peptide substrates on serine and tyrosine, and its activity is enhanced by tyrosine auto-phosphorylation. Nek10 dual-specificity depends on residues in the HRD+2 and APE-4 positions that are uncommon in either serine/threonine or tyrosine kinases. Finally, we show that the phosphorylation-site motifs for the mitotic kinases Nek6, Nek7 and Nek9 are essentially identical to that of their upstream activator Plk1, suggesting that Nek6/7/9 function as phospho-motif amplifiers of Plk1 signaling.

Engineering allosteric regulation in protein kinases.

Phosphoregulation, in which the addition of a negatively charged phosphate group modulates protein activity, enables dynamic cellular responses. To understand how new phosphoregulation might be acquired, we mutationally scanned the surface of a prototypical yeast kinase (Kss1) to identify potential regulatory sites. The data revealed a set of spatially distributed "hotspots" that might have coevolved with the active site and preferentially modulated kinase activity. By engineering simple consensus phosphorylation sites at these hotspots, we rewired cell signaling in yeast. Using the same approach with a homolog yeast mitogen-activated protein kinase, Hog1, we introduced new phosphoregulation that modified its localization and signaling dynamics. Beyond revealing potential use in synthetic biology, our findings suggest that the identified hotspots contribute to the diversity of natural allosteric regulatory mechanisms in the eukaryotic kinome and, given that some are mutated in cancers, understanding these hotspots may have clinical relevance to human disease.

Hierarchical Organization Endows the Kinase Domain with Regulatory Plasticity.

The functional diversity of kinases enables specificity in cellular signal transduction. Yet how more than 500 members of the human kinome specifically receive regulatory inputs and convey information to appropriate substrates-all while using the common signaling output of phosphorylation-remains enigmatic. Here, we perform statistical co-evolution analysis, mutational scanning, and quantitative live-cell assays to reveal a hierarchical organization of the kinase domain that facilitates the orthogonal evolution of regulatory inputs and substrate outputs while maintaining catalytic function. We find that three quasi-independent "sectors"-groups of evolutionarily coupled residues-represent functional units in the kinase domain that encode for catalytic activity, substrate specificity, and regulation. Sector positions impact both disease and pharmacology: the catalytic sector is significantly enriched for somatic cancer mutations, and residues in the regulatory sector interact with allosteric kinase inhibitors. We propose that this functional architecture endows the kinase domain with inherent regulatory plasticity.

Exploiting Temporal Collateral Sensitivity in Tumor Clonal Evolution.

The prevailing approach to addressing secondary drug resistance in cancer focuses on treating the resistance mechanisms at relapse. However, the dynamic nature of clonal evolution, along with potential fitness costs and cost compensations, may present exploitable vulnerabilities-a notion that we term "temporal collateral sensitivity." Using a combined pharmacological screen and drug resistance selection approach in a murine model of Ph(+) acute lymphoblastic leukemia, we indeed find that temporal and/or persistent collateral sensitivity to non-classical BCR-ABL1 drugs arises in emergent tumor subpopulations during the evolution of resistance toward initial treatment with BCR-ABL1-targeted inhibitors. We determined the sensitization mechanism via genotypic, phenotypic, signaling, and binding measurements in combination with computational models and demonstrated significant overall survival extension in mice. Additional stochastic mathematical models and small-molecule screens extended our insights, indicating the value of focusing on evolutionary trajectories and pharmacological profiles to identify new strategies to treat dynamic tumor vulnerabilities.

Kinome-wide decoding of network-attacking mutations rewiring cancer signaling.

Cancer cells acquire pathological phenotypes through accumulation of mutations that perturb signaling networks. However, global analysis of these events is currently limited. Here, we identify six types of network-attacking mutations (NAMs), including changes in kinase and SH2 modulation, network rewiring, and the genesis and extinction of phosphorylation sites. We developed a computational platform (ReKINect) to identify NAMs and systematically interpreted the exomes and quantitative (phospho-)proteomes of five ovarian cancer cell lines and the global cancer genome repository. We identified and experimentally validated several NAMs, including PKCγ M501I and PKD1 D665N, which encode specificity switches analogous to the appearance of kinases de novo within the kinome. We discover mutant molecular logic gates, a drift toward phospho-threonine signaling, weakening of phosphorylation motifs, and kinase-inactivating hotspots in cancer. Our method pinpoints functional NAMs, scales with the complexity of cancer genomes and cell signaling, and may enhance our capability to therapeutically target tumor-specific networks.

Unmasking determinants of specificity in the human kinome.

Protein kinases control cellular responses to environmental cues by swift and accurate signal processing. Breakdowns in this high-fidelity capability are a driving force in cancer and other diseases. Thus, our limited understanding of which amino acids in the kinase domain encode substrate specificity, the so-called determinants of specificity (DoS), constitutes a major obstacle in cancer signaling. Here, we systematically discover several DoS and experimentally validate three of them, named the αC1, αC3, and APE-7 residues. We demonstrate that DoS form sparse networks of non-conserved residues spanning distant regions. Our results reveal a likely role for inter-residue allostery in specificity and an evolutionary decoupling of kinase activity and specificity, which appear loaded on independent groups of residues. Finally, we uncover similar properties driving SH2 domain specificity and demonstrate how the identification of DoS can be utilized to elucidate a greater understanding of the role of signaling networks in cancer (Creixell et al., 2015 [this issue of Cell]).

Pathway and network analysis of cancer genomes.

Genomic information on tumors from 50 cancer types cataloged by the International Cancer Genome Consortium (ICGC) shows that only a few well-studied driver genes are frequently mutated, in contrast to many infrequently mutated genes that may also contribute to tumor biology. Hence there has been large interest in developing pathway and network analysis methods that group genes and illuminate the processes involved. We provide an overview of these analysis techniques and show where they guide mechanistic and translational investigations.

Integrative analysis of kinase networks in TRAIL-induced apoptosis provides a source of potential targets for combination therapy.

Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) is an endogenous secreted peptide and, in preclinical studies, preferentially induces apoptosis in tumor cells rather than in normal cells. The acquisition of resistance in cells exposed to TRAIL or its mimics limits their clinical efficacy. Because kinases are intimately involved in the regulation of apoptosis, we systematically characterized kinases involved in TRAIL signaling. Using RNA interference (RNAi) loss-of-function and cDNA overexpression screens, we identified 169 protein kinases that influenced the dynamics of TRAIL-induced apoptosis in the colon adenocarcinoma cell line DLD-1. We classified the kinases as sensitizers or resistors or modulators, depending on the effect that knockdown and overexpression had on TRAIL-induced apoptosis. Two of these kinases that were classified as resistors were PX domain-containing serine/threonine kinase (PXK) and AP2-associated kinase 1 (AAK1), which promote receptor endocytosis and may enable cells to resist TRAIL-induced apoptosis by enhancing endocytosis of the TRAIL receptors. We assembled protein interaction maps using mass spectrometry-based protein interaction analysis and quantitative phosphoproteomics. With these protein interaction maps, we modeled information flow through the networks and identified apoptosis-modifying kinases that are highly connected to regulated substrates downstream of TRAIL. The results of this analysis provide a resource of potential targets for the development of TRAIL combination therapies to selectively kill cancer cells.

CoreFlow: a computational platform for integration, analysis and modeling of complex biological data.

A major challenge in mass spectrometry and other large-scale applications is how to handle, integrate, and model the data that is produced. Given the speed at which technology advances and the need to keep pace with biological experiments, we designed a computational platform, CoreFlow, which provides programmers with a framework to manage data in real-time. It allows users to upload data into a relational database (MySQL), and to create custom scripts in high-level languages such as R, Python, or Perl for processing, correcting and modeling this data. CoreFlow organizes these scripts into project-specific pipelines, tracks interdependencies between related tasks, and enables the generation of summary reports as well as publication-quality images. As a result, the gap between experimental and computational components of a typical large-scale biology project is reduced, decreasing the time between data generation, analysis and manuscript writing. CoreFlow is being released to the scientific community as an open-sourced software package complete with proteomics-specific examples, which include corrections for incomplete isotopic labeling of peptides (SILAC) or arginine-to-proline conversion, and modeling of multiple/selected reaction monitoring (MRM/SRM) results. BIOLOGICAL SIGNIFICANCE: CoreFlow was purposely designed as an environment for programmers to rapidly perform data analysis. These analyses are assembled into project-specific workflows that are readily shared with biologists to guide the next stages of experimentation. Its simple yet powerful interface provides a structure where scripts can be written and tested virtually simultaneously to shorten the life cycle of code development for a particular task. The scripts are exposed at every step so that a user can quickly see the relationships between the data, the assumptions that have been made, and the manipulations that have been performed. Since the scripts use commonly available programming languages, they can easily be transferred to and from other computational environments for debugging or faster processing. This focus on 'on the fly' analysis sets CoreFlow apart from other workflow applications that require wrapping of scripts into particular formats and development of specific user interfaces. Importantly, current and future releases of data analysis scripts in CoreFlow format will be of widespread benefit to the proteomics community, not only for uptake and use in individual labs, but to enable full scrutiny of all analysis steps, thus increasing experimental reproducibility and decreasing errors. This article is part of a Special Issue entitled: Can Proteomics Fill the Gap Between Genomics and Phenotypes?

Computational approaches to identify functional genetic variants in cancer genomes.

The International Cancer Genome Consortium (ICGC) aims to catalog genomic abnormalities in tumors from 50 different cancer types. Genome sequencing reveals hundreds to thousands of somatic mutations in each tumor but only a minority of these drive tumor progression. We present the result of discussions within the ICGC on how to address the challenge of identifying mutations that contribute to oncogenesis, tumor maintenance or response to therapy, and recommend computational techniques to annotate somatic variants and predict their impact on cancer phenotype.

Navigating cancer network attractors for tumor-specific therapy.

Cells employ highly dynamic signaling networks to drive biological decision processes. Perturbations to these signaling networks may attract cells to new malignant signaling and phenotypic states, termed cancer network attractors, that result in cancer development. As different cancer cells reach these malignant states by accumulating different molecular alterations, uncovering these mechanisms represents a grand challenge in cancer biology. Addressing this challenge will require new systems-based strategies that capture the intrinsic properties of cancer signaling networks and provide deeper understanding of the processes by which genetic lesions perturb these networks and lead to disease phenotypes. Network biology will help circumvent fundamental obstacles in cancer treatment, such as drug resistance and metastasis, empowering personalized and tumor-specific cancer therapies.

Mutational properties of amino acid residues: implications for evolvability of phosphorylatable residues.

As François Jacob pointed out over 30 years ago, evolution is a tinkering process, and, as such, relies on the genetic diversity produced by mutation subsequently shaped by Darwinian selection. However, there is one implicit assumption that is made when studying this tinkering process; it is typically assumed that all amino acid residues are equally likely to mutate or to result from a mutation. Here, by reconstructing ancestral sequences and computing mutational probabilities for all the amino acid residues, we refute this assumption and show extensive inequalities between different residues in terms of their mutational activity. Moreover, we highlight the importance of the genetic code and physico-chemical properties of the amino acid residues as likely causes of these inequalities and uncover serine as a mutational hot spot. Finally, we explore the consequences that these different mutational properties have on phosphorylation site evolution, showing that a higher degree of evolvability exists for phosphorylated threonine and, to a lesser extent, serine in comparison with tyrosine residues. As exemplified by the suppression of serine's mutational activity in phosphorylation sites, our results suggest that the cell can fine-tune the mutational activities of amino acid residues when they reside in functional protein regions.