Publications

Summary

Traditional deep mutational scanning (DMS) encodes every single amino-acid substitution from a wild-type sequence. We hypothesize that combinatorial DMS (CDMS) libraries, incorporating all mutations in all combinations, can enable the discovery of high-affinity protein (super)binders by capturing epistatic, non-linear amino-acid interactions. Here, we introduce origin-independent and context-exhaustive high-throughput integration of combinatorial DMS libraries (ORCHID), which systematically maps regions of wild-type-independent epistasis across all mutational contexts and trajectories. For benchmarking, we build a high-throughput peptide display assay measuring PIN1^WW-domain affinity for a CDMS peptide library containing phosphoserine via amber codon suppression. ORCHID raises prediction accuracy by 45% over non-epistatic models. We identify and structurally characterize two epistatic superbinders, SPY-tide and LYR-tide, binding 3- to 5-fold tighter than current optimal PIN1^WW-domain binders through “fold-and-turn” and register-shifted conformational changes. We also identify two molecular determinants of epistasis, PIN1^F25 and PIN1^R14, that natively encode non-linear binding and, when mutated, abolish it. Hence, natural proteins recognize peptides non-linearly, offering opportunities for improved binder design.

Abstract

Targeted therapies for receptor tyrosine kinases are effective but invariably limited by drug resistance. In EGFR-mutant lung cancer, AXL activation drives resistance to erlotinib by restoring cell survival and migration. To map these signaling mechanisms, we generated a panel of lung adenocarcinoma PC9 cell lines with phenylalanine substitutions at each intracellular AXL tyrosine residue. By integrating phosphorylation data with phenotypic changes via multivariate modeling, we found that AXL signaling organizes into two clusters enriched for Abl1 and SFK substrate motifs. A peptide specificity screen identified FAK1 as a top proximal substrate of AXL. Downstream, AXL-mediated YAP1 activation was found to sustain drug tolerance, while combined inhibition eliminated persister cells in vitro. These AXL and YAP pathways correlate with disease progression and poor clinical outcomes in erlotinib-treated patients. Collectively, this study dissects the specific signaling regulators by which AXL drives erlotinib resistance.

Abstract: Dynamic regulation of protein tyrosine phosphorylation (pTyr) by kinases and phosphatases enables cells to sense and respond to environmental changes. The widely used chemical pervanadate induces the accumulation of pTyr in mammalian cell lines. This effect is primarily attributed to its inhibition of protein tyrosine phosphatases (PTPs), leading to the assertion that PTPs are master gatekeepers of intracellular pTyr homeostasis. Here, we used several approaches to reveal that pervanadate disrupted cellular redox homeostasis and directly activated tyrosine kinases of the SRC family through the oxidation of specific cysteine residues. Mass spectrometry and biophysical approaches showed that pervanadate-induced oxidation of cysteine-188 and cysteine-280 activated SRC by disrupting autoinhibitory intramolecular interactions between the catalytic domain and the SH2/SH3 domains and by impairing SH2 domain binding to phosphopeptides, including the regulatory carboxyl-terminal tail phosphotyrosine-530. Redox-sensitive cysteine residues were essential for SRC to promote the overgrowth of mouse fibroblasts. Our findings call for a reevaluation of pervanadate-based experiments and demonstrate that SRC cysteines control its oncogenic properties.

Abstract: Animal germ cells deploy a specialized small RNA-based silencing system, called the PIWI-interacting RNA (piRNA) pathway, to prevent unwanted expression of transposable elements (TEs) and maintain genome integrity. In Drosophila melanogaster germ cells, the majority of piRNA populations originate from dual-strand piRNA clusters, genomic regions highly enriched in TE fragments, via an elaborate machinery centered on the Heterochromatin Protein 1 homolog, Rhino. Although Rhino binds to peptides carrying tri-methylated H3K9 in vitro, it is not fully understood why in vivo only a fraction of H3K9me3-decorated heterochromatin is occupied by Rhino. Recent work revealed that Rhino is recruited to a subset of piRNA clusters by Kipferl. Here we identify a Kipferl-independent mode of Rhino recruitment that, in addition to the previously established role of H3K9me3, also depends on the histone H3 lysine 27 methyltransferase Enhancer of Zeste. At Kipferl-independent sites, we find that Rhino specifically binds to loci marked by both H3K9me3 and H3K27me3 via its chromodomain. Although the exact mechanism of how Rhino binding is influenced by dual histone modifications remains unclear from a structural and biochemical perspective, our work suggests that combinatorial modifications may regulate the specificity of chromatin-binding protein interactions. These findings provide an enhanced understanding of how Rhino targets piRNA source loci, highlighting the sophisticated epigenetic landscape governing TE silencing in Drosophila germ cells.

Abstract: During the repair of interstrand crosslinks (ICLs) a DNA double-strand break (DSB) is generated. The Fanconi anemia (FA) core complex, which is recruited to ICLs, promotes high-fidelity repair of this DSB by homologous recombination (HR). However, whether the FA core complex also promotes HR at ICL-independent DSBs, for example induced by ionizing irradiation or nucleases, remains controversial. Here, we identified the FA core complex members FANCL and Ube2T as HR-promoting factors in a CRISPR/Cas9-based screen. Using isogenic cell line models, we further demonstrated an HR-promoting function of FANCL and Ube2T, and of their ubiquitination substrate FANCD2. We show that FANCL and Ube2T localize at DSBs in a FANCM-dependent manner, and are required for the DSB accumulation of FANCD2. Mechanistically, we demonstrate that FANCL ubiquitin ligase activity is required for the accumulation of CtIP at DSBs, thereby promoting end resection and Rad51 loading. Together, these data demonstrate a dual genome maintenance function of the FA core complex and FANCD2 in promoting repair of both ICLs and DSBs. 

Abstract: Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood. 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.

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.

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.

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.

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.