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.

An atlas of substrate specificities for the human serine/threonine kinome.

Protein phosphorylation is one of the most widespread post-translational modifications in biology1,2. With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes3,4. For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible3. Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways.

Host protein kinases required for SARS-CoV-2 nucleocapsid phosphorylation and viral replication.

Multiple coronaviruses have emerged independently in the past 20 years that cause lethal human diseases. Although vaccine development targeting these viruses has been accelerated substantially, there remain patients requiring treatment who cannot be vaccinated or who experience breakthrough infections. Understanding the common host factors necessary for the life cycles of coronaviruses may reveal conserved therapeutic targets. Here, we used the known substrate specificities of mammalian protein kinases to deconvolute the sequence of phosphorylation events mediated by three host protein kinase families (SRPK, GSK-3, and CK1) that coordinately phosphorylate a cluster of serine and threonine residues in the viral N protein, which is required for viral replication. We also showed that loss or inhibition of SRPK1/2, which we propose initiates the N protein phosphorylation cascade, compromised the viral replication cycle. Because these phosphorylation sites are highly conserved across coronaviruses, inhibitors of these protein kinases not only may have therapeutic potential against COVID-19 but also may be broadly useful against coronavirus-mediated diseases.

The FDA-approved drug Alectinib compromises SARS-CoV-2 nucleocapsid phosphorylation and inhibits viral infection in vitro.

While vaccines are vital for preventing COVID-19 infections, it is critical to develop new therapies to treat patients who become infected. Pharmacological targeting of a host factor required for viral replication can suppress viral spread with a low probability of viral mutation leading to resistance. In particular, host kinases are highly druggable targets and a number of conserved coronavirus proteins, notably the nucleoprotein (N), require phosphorylation for full functionality. In order to understand how targeting kinases could be used to compromise viral replication, we used a combination of phosphoproteomics and bioinformatics as well as genetic and pharmacological kinase inhibition to define the enzymes important for SARS-CoV-2 N protein phosphorylation and viral replication. From these data, we propose a model whereby SRPK1/2 initiates phosphorylation of the N protein, which primes for further phosphorylation by GSK-3a/b and CK1 to achieve extensive phosphorylation of the N protein SR-rich domain. Importantly, we were able to leverage our data to identify an FDA-approved kinase inhibitor, Alectinib, that suppresses N phosphorylation by SRPK1/2 and limits SARS-CoV-2 replication. Together, these data suggest that repurposing or developing novel host-kinase directed therapies may be an efficacious strategy to prevent or treat COVID-19 and other coronavirus-mediated diseases.

Systematic analysis of the intersection of disease mutations with protein modifications.

BACKGROUND: Perturbed posttranslational modification (PTM) landscapes commonly cause pathological phenotypes. The Cancer Genome Atlas (TCGA) project profiles thousands of tumors allowing the identification of spontaneous cancer-driving mutations, while Uniprot and dbSNP manage genetic disease-associated variants in the human population. PhosphoSitePlus (PSP) is the most comprehensive resource for studying experimentally observed PTM sites and the only repository with daily updates on functional annotations for many of these sites. To elucidate altered PTM landscapes on a large scale, we integrated disease-associated mutations from TCGA, Uniprot, and dbSNP with PTM sites from PhosphoSitePlus. We characterized each dataset individually, compared somatic with germline mutations, and analyzed PTM sites intersecting directly with disease variants. To assess the impact of mutations in the flanking regions of phosphosites, we developed DeltaScansite, a pipeline that compares Scansite predictions on wild type versus mutated sequences. Disease mutations are also visualized in PhosphoSitePlus. RESULTS: Characterization of somatic variants revealed oncoprotein-like mutation profiles of U2AF1, PGM5, and several other proteins, showing alteration patterns similar to germline mutations. The union of all datasets uncovered previously unknown losses and gains of PTM events in diseases unevenly distributed across different PTM types. Focusing on phosphorylation, our DeltaScansite workflow predicted perturbed signaling networks consistent with calculations by the machine learning method MIMP. CONCLUSIONS: We discovered oncoprotein-like profiles in TCGA and mutations that presumably modify protein function by impacting PTM sites directly or by rewiring upstream regulation. The resulting datasets are enriched with functional annotations from PhosphoSitePlus and present a unique resource for potential biomarkers or disease drivers.

15 years of PhosphoSitePlus®: integrating post-translationally modified sites, disease variants and isoforms.

For 15 years the mission of PhosphoSitePlus® (PSP, https://www.phosphosite.org) has been to provide comprehensive information and tools for the study of mammalian post-translational modifications (PTMs). The number of unique PTMs in PSP is now more than 450 000 from over 22 000 articles and thousands of MS datasets. The most important areas of growth in PSP are in disease and isoform informatics. Germline mutations associated with inherited diseases and somatic cancer mutations have been added to the database and can now be viewed along with PTMs and associated quantitative information on novel 'lollipop' plots. These plots enable researchers to interactively visualize the overlap between disease variants and PTMs, and to identify mutations that may alter phenotypes by rewiring signaling networks. We are expanding the sequence space to include over 30 000 human and mouse isoforms to enable researchers to explore the important but understudied biology of isoforms. This represents a necessary expansion of sequence space to accommodate the growing precision and depth of coverage enabled by ongoing advances in mass spectrometry. Isoforms are aligned using a new algorithm. Exploring the worlds of PTMs and disease mutations in the entire isoform space will hopefully lead to new biomarkers, therapeutic targets, and insights into isoform biology.

A Curated Resource for Phosphosite-specific Signature Analysis.

Signaling pathways are orchestrated by post-translational modifications (PTMs) such as phosphorylation. However, pathway analysis of PTM data sets generated by mass spectrometry (MS)-based proteomics is typically performed at a gene-centric level because of the lack of appropriately curated PTM signature databases and bioinformatic tools that leverage PTM site-specific information. Here we present the first version of PTMsigDB, a database of modification site-specific signatures of perturbations, kinase activities and signaling pathways curated from more than 2,500 publications. We adapted the widely used single sample Gene Set Enrichment Analysis approach to utilize PTMsigDB, enabling PTMSignature Enrichment Analysis (PTM-SEA) of quantitative MS data. We used a well-characterized data set of epidermal growth factor (EGF)-perturbed cancer cells to evaluate our approach and demonstrated better representation of signaling events compared with gene-centric methods. We then applied PTM-SEA to analyze the phosphoproteomes of cancer cells treated with cell-cycle inhibitors and detected mechanism-of-action specific signatures of cell cycle kinases. We also applied our methods to analyze the phosphoproteomes of PI3K-inhibited human breast cancer cells and detected signatures of compounds inhibiting PI3K as well as targets downstream of PI3K (AKT, MAPK/ERK) covering a substantial fraction of the PI3K pathway. PTMsigDB and PTM-SEA can be freely accessed at https://github.com/broadinstitute/ssGSEA2.0.

Integration of protein phosphorylation, acetylation, and methylation data sets to outline lung cancer signaling networks.

Protein posttranslational modifications (PTMs) have typically been studied independently, yet many proteins are modified by more than one PTM type, and cell signaling pathways somehow integrate this information. We coupled immunoprecipitation using PTM-specific antibodies with tandem mass tag (TMT) mass spectrometry to simultaneously examine phosphorylation, methylation, and acetylation in 45 lung cancer cell lines compared to normal lung tissue and to cell lines treated with anticancer drugs. This simultaneous, large-scale, integrative analysis of these PTMs using a cluster-filtered network (CFN) approach revealed that cell signaling pathways were outlined by clustering patterns in PTMs. We used the t-distributed stochastic neighbor embedding (t-SNE) method to identify PTM clusters and then integrated each with known protein-protein interactions (PPIs) to elucidate functional cell signaling pathways. The CFN identified known and previously unknown cell signaling pathways in lung cancer cells that were not present in normal lung epithelial tissue. In various proteins modified by more than one type of PTM, the incidence of those PTMs exhibited inverse relationships, suggesting that molecular exclusive "OR" gates determine a large number of signal transduction events. We also showed that the acetyltransferase EP300 appears to be a hub in the network of pathways involving different PTMs. In addition, the data shed light on the mechanism of action of geldanamycin, an HSP90 inhibitor. Together, the findings reveal that cell signaling pathways mediated by acetylation, methylation, and phosphorylation regulate the cytoskeleton, membrane traffic, and RNA binding protein-mediated control of gene expression.

Splice-Aware Multiple Sequence Alignment of Protein Isoforms.

Multiple sequence alignment (MSA) is a classic problem in computational genomics. In typical use, MSA software is expected to align a collection of homologous genes, such as orthologs from multiple species or duplication-induced paralogs within a species. Recent focus on the importance of alternatively-spliced isoforms in disease and cell biology has highlighted the need to create MSAs that more effectively accommodate isoforms. MSAs are traditionally constructed using scoring criteria that prefer alignments with occasional mismatches over alignments with long gaps. Alternatively spliced protein isoforms effectively contain exon-length insertions or deletions (indels) relative to each other, and demand an alternative approach. Some improvements can be achieved by making indel penalties much smaller, but this is merely a patchwork solution. In this work we present Mirage, a novel MSA software package for the alignment of alternatively spliced protein isoforms. Mirage aligns isoforms to each other by first mapping each protein sequence to its encoding genomic sequence, and then aligning isoforms to one another based on the relative genomic coordinates of their constitutive codons. Mirage is highly effective at mapping proteins back to their encoding exons, and these protein-genome mappings lead to extremely accurate intra-species alignments; splice site information in these alignments is used to improve the accuracy of inter-species alignments of isoforms. Mirage alignments have also revealed the ubiquity of dual-coding exons, in which an exon conditionally encodes multiple open reading frames as overlapping spliced segments of frame-shifted genomic sequence.

Crk Tyrosine Phosphorylation Regulates PDGF-BB-inducible Src Activation and Breast Tumorigenicity and Metastasis.

The activity of Src family kinases (Src being the prototypical member) is tightly regulated by differential phosphorylation on Tyr416 (positive) and Tyr527 (negative), a duet that reciprocally regulates kinase activity. The latter negative regulation of Src on Tyr527 is mediated by C-terminal Src kinase (CSK) that phosphorylates Tyr527 and maintains Src in a clamped negative regulated state by promoting an intramolecular association. Here it is demonstrated that the SH2- and SH3-domain containing adaptor protein CrkII, by virtue of its phosphorylation on Tyr239, regulates the Csk/Src signaling axis to control Src activation. Once phosphorylated, the motif (PIpYARVIQ) forms a consensus sequence for the SH2 domain of CSK to form a pTyr239-CSK complex. Functionally, when expressed in Crk-/- MEFs or in Crk+/+ HS683 cells, Crk Y239F delayed PDGF-BB-inducible Src Tyr416 phosphorylation. Moreover, expression of Crk Y239F in HS683 cells delayed Src kinase activation and suppressed the cell-invasive and -transforming phenotypes. Finally, through loss-of-function and epistasis experiments using CRISPR-Cas9-engineered 4T1 murine breast cancer cells, Crk Tyr239 is implicated in breast cancer tumor growth and metastasis in orthotopic immunocompetent 4T1 mice model of breast adenocarcinoma. These findings delineate a novel role for Crk Tyr239 phosphorylation in the regulation of Src kinases, as well as a potential molecular explanation for a long-standing question as to how Crk regulates the activation of Src kinases.Implications: These findings provide new perspectives on the versatility of Crk in cancer by demonstrating how Crk mechanistically drives, through a tyrosine phosphorylation-dependent manner, tumor growth, and metastasis. Mol Cancer Res; 16(1); 173-83. ©2017 AACR.

Clustergrammer, a web-based heatmap visualization and analysis tool for high-dimensional biological data.

Most tools developed to visualize hierarchically clustered heatmaps generate static images. Clustergrammer is a web-based visualization tool with interactive features such as: zooming, panning, filtering, reordering, sharing, performing enrichment analysis, and providing dynamic gene annotations. Clustergrammer can be used to generate shareable interactive visualizations by uploading a data table to a web-site, or by embedding Clustergrammer in Jupyter Notebooks. The Clustergrammer core libraries can also be used as a toolkit by developers to generate visualizations within their own applications. Clustergrammer is demonstrated using gene expression data from the cancer cell line encyclopedia (CCLE), original post-translational modification data collected from lung cancer cells lines by a mass spectrometry approach, and original cytometry by time of flight (CyTOF) single-cell proteomics data from blood. Clustergrammer enables producing interactive web based visualizations for the analysis of diverse biological data.

Isotype Determination of Antibodies.

On identifying a new monoclonal antibody, or in characterizing antibodies in sera evoked by disease or immunization, it is particularly informative to determine the serological class of the antibodies. The serological class of the protein is determined by the structure of the antibody constant region. Several methods of class or isotype determination are outlined in this unit: sandwich ELISA, electrophoresis, and immunofixation, or use of a variety of commercially available kits. Different purification schemes and approaches for enzymatic fragmentation of the antibodies depend on the class or isotype of the antibody, so this information streamlines these processes. © 2017 by John Wiley & Sons, Inc.

Double-Immunodiffusion Assay for Detecting Specific Antibodies (Ouchterlony).

The method first described by Ouchterlony in 1948 is a classic and simple technique that permits evaluation and comparison of antibodies in animal or human sera directed against protein or complex carbohydrate antigens. It is a low-tech procedure that may provide information on the relative quantity of antibody activity and the nature of the antigenic epitopes in different preparations. © 2017 by John Wiley & Sons, Inc.

Reciprocal regulation of Abl kinase by Crk Y251 and Abi1 controls invasive phenotypes in glioblastoma.

Crk is the prototypical member of a class of Src homology 2 (SH2) and Src homology 3 (SH3) domain-containing adaptor proteins that positively regulate cell motility via the activation of Rac1 and, in certain tumor types such as GBM, can promote cell invasion and metastasis by mechanisms that are not well understood. Here we demonstrate that Crk, via its phosphorylation at Tyr251, promotes invasive behavior of tumor cells, is a prominent feature in GBM, and correlating with aggressive glioma grade IV staging and overall poor survival outcomes. At the molecular level, Tyr251 phosphorylation of Crk is negatively regulated by Abi1, which competes for Crk binding to Abl and attenuates Abl transactivation. Together, these results show that Crk and Abi1 have reciprocal biological effects and act as a molecular rheostat to control Abl activation and cell invasion. Finally, these data suggest that Crk Tyr251 phosphorylation regulate invasive cell phenotypes and may serve as a biomarker for aggressive GBM.

Enzyme-Linked Immunosorbent Assays.

This unit describes six different ELISA systems for the detection of specific antibodies, soluble antigens, or cell-surface antigens. In all six systems, soluble reactants are removed from solution after specifically binding to solid-phase reactants. In the first four protocols, solid-phase reactants are prepared by adsorbing an antigen or antibody onto plastic microtiter plates; in the next two protocols, the solid-phase reactants are cell-associated molecules. In all protocols, the solid-phase reagents are incubated with secondary or tertiary reactants covalently coupled to an enzyme. Unbound conjugates are washed out and a chromogenic or fluorogenic substrate is added. As the substrate is hydrolyzed by the bound enzyme conjugate, a colored or fluorescent product is generated. Finally, the product is detected visually or with a microtiter plate reader. The amount of product generated is proportional to the amount of analysate in the test mixture. One of the support protocols can be used to optimize the different ELISAs. A second support protocol presents a method for preparing alkaline phosphatase conjugates.

PhosphoSitePlus, 2014: mutations, PTMs and recalibrations.

PhosphoSitePlus(®) (PSP, http://www.phosphosite.org/), a knowledgebase dedicated to mammalian post-translational modifications (PTMs), contains over 330,000 non-redundant PTMs, including phospho, acetyl, ubiquityl and methyl groups. Over 95% of the sites are from mass spectrometry (MS) experiments. In order to improve data reliability, early MS data have been reanalyzed, applying a common standard of analysis across over 1,000,000 spectra. Site assignments with P > 0.05 were filtered out. Two new downloads are available from PSP. The 'Regulatory sites' dataset includes curated information about modification sites that regulate downstream cellular processes, molecular functions and protein-protein interactions. The 'PTMVar' dataset, an intersect of missense mutations and PTMs from PSP, identifies over 25,000 PTMVars (PTMs Impacted by Variants) that can rewire signaling pathways. The PTMVar data include missense mutations from UniPROTKB, TCGA and other sources that cause over 2000 diseases or syndromes (MIM) and polymorphisms, or are associated with hundreds of cancers. PTMVars include 18 548 phosphorlyation sites, 3412 ubiquitylation sites, 2316 acetylation sites, 685 methylation sites and 245 succinylation sites.

Signatures of natural selection on mutations of residues with multiple posttranslational modifications.

Posttranslational modifications (PTMs) regulate molecular structures and functions of proteins by covalently binding to amino acids. Hundreds of thousands of PTMs have been reported for the human proteome, with multiple PTMs known to affect tens of thousands of lysine (K) residues. Our molecular evolutionary analyses show that K residues with multiple PTMs exhibit greater conservation than those with a single PTM, but the difference is rather small. In contrast, short-term evolutionary trends revealed in an analysis of human population variation exhibited a much larger difference. Lysine residues with three PTMs show 1.8-fold enrichment of Mendelian disease-associated variants when compared with K residues with two PTMs, with the latter showing 1.7-fold enrichment of these variants when compared with the K residues with one PTM. Rare polymorphisms in humans show a similar trend, which suggests much greater negative selection against mutations of K residues with multiple PTMs within population. Conversely, common polymorphisms are overabundant at unmodified K residues and at K residues with fewer PTMs. The observed difference between inter- and intraspecies patterns of purifying selection on residues with PTMs suggests extensive species-specific drifting of PTM positions. These results suggest that the functionality of a protein is likely conserved, without necessarily conserving the PTM positions over evolutionary time.

PhosphoSitePlus: a comprehensive resource for investigating the structure and function of experimentally determined post-translational modifications in man and mouse.

PhosphoSitePlus (http://www.phosphosite.org) is an open, comprehensive, manually curated and interactive resource for studying experimentally observed post-translational modifications, primarily of human and mouse proteins. It encompasses 1,30,000 non-redundant modification sites, primarily phosphorylation, ubiquitinylation and acetylation. The interface is designed for clarity and ease of navigation. From the home page, users can launch simple or complex searches and browse high-throughput data sets by disease, tissue or cell line. Searches can be restricted by specific treatments, protein types, domains, cellular components, disease, cell types, cell lines, tissue and sequences or motifs. A few clicks of the mouse will take users to substrate pages or protein pages with sites, sequences, domain diagrams and molecular visualization of side-chains known to be modified; to site pages with information about how the modified site relates to the functions of specific proteins and cellular processes and to curated information pages summarizing the details from one record. PyMOL and Chimera scripts that colorize reactive groups on residues that are modified can be downloaded. Features designed to facilitate proteomic analyses include downloads of modification sites, kinase-substrate data sets, sequence logo generators, a Cytoscape plugin and BioPAX download to enable pathway visualization of the kinase-substrate interactions in PhosphoSitePlus®.

The BioPAX community standard for pathway data sharing.

Biological Pathway Exchange (BioPAX) is a standard language to represent biological pathways at the molecular and cellular level and to facilitate the exchange of pathway data. The rapid growth of the volume of pathway data has spurred the development of databases and computational tools to aid interpretation; however, use of these data is hampered by the current fragmentation of pathway information across many databases with incompatible formats. BioPAX, which was created through a community process, solves this problem by making pathway data substantially easier to collect, index, interpret and share. BioPAX can represent metabolic and signaling pathways, molecular and genetic interactions and gene regulation networks. Using BioPAX, millions of interactions, organized into thousands of pathways, from many organisms are available from a growing number of databases. This large amount of pathway data in a computable form will support visualization, analysis and biological discovery.

Akt-RSK-S6 kinase signaling networks activated by oncogenic receptor tyrosine kinases.

Receptor tyrosine kinases (RTKs) activate pathways mediated by serine-threonine kinases, such as the PI3K (phosphatidylinositol 3-kinase)-Akt pathway, the Ras-MAPK (mitogen-activated protein kinase)-RSK (ribosomal S6 kinase) pathway, and the mTOR (mammalian target of rapamycin)-p70 S6 pathway, that control important aspects of cell growth, proliferation, and survival. The Akt, RSK, and p70 S6 family of protein kinases transmits signals by phosphorylating substrates on an RxRxxS/T motif (R, arginine; S, serine; T, threonine; and x, any amino acid). We developed a large-scale proteomic approach to identify more than 300 substrates of this kinase family in cancer cell lines driven by the c-Met, epidermal growth factor receptor (EGFR), or platelet-derived growth factor receptor alpha (PDGFRalpha) RTKs. We identified a subset of proteins with RxRxxS/T sites for which phosphorylation was decreased by RTK inhibitors (RTKIs), as well as by inhibitors of the PI3K, mTOR, and MAPK pathways, and we determined the effects of small interfering RNA directed against these substrates on cell viability. Phosphorylation of the protein chaperone SGTA (small glutamine-rich tetratricopeptide repeat-containing protein alpha) at serine-305 was essential for PDGFRalpha stabilization and cell survival in PDGFRalpha-dependent cancer cells. Our approach provides a new view of RTK and Akt-RSK-S6 kinase signaling, revealing previously unidentified Akt-RSK-S6 kinase substrates that merit further consideration as targets for combination therapy with RTKIs.