Dynamic Phosphoproteome Profiling of Zebrafish Embryonic Fibroblasts during Cold Acclimation

Temperature affects almost all aspects of the fish life. To cope with low temperature, fish have evolved the ability of cold acclimation for survival. However, intracellular signaling events underlying cold acclimation in fish remain largely unknown. Here, the formation of cold acclimation in zebrafish embryonic fibroblasts (ZF4) is monitored and the phosphorylation events during the process are investigated through a large‐scale quantitative phosphoproteomic approach. In total, 11 474 phosphorylation sites are identified on 4066 proteins and quantified 5772 phosphosites on 2519 proteins. Serine, threonine, and tyrosine (Ser/Thr/Tyr) phosphorylation accounted for 85.5%, 13.3%, and 1.2% of total phosphosites, respectively. Among all phosphosites, 702 phosphosites on 510 proteins show differential regulation during cold acclimation of ZF4 cells. These phosphosites are divided into six clusters according to their dynamic changes during cold exposure. Kinase–substrate prediction reveals that mitogen‐activated protein kinase (MAPK) among the kinase groups is predominantly responsible for phosphorylation of these phosphosites. The differentially regulated phosphoproteins are functionally associated with various cellular processes such as regulation of actin cytoskeleton and MAPK signaling pathway. These data enrich the database of protein phosphorylation sites in zebrafish and provide key clues for the elucidation of intracellular signaling networks during cold acclimation of fish.     

Evolution of protein phosphorylation for distinct functional modules in vertebrate genomes

Recent publications have revealed that the evolution of phosphosites is influenced by the local protein structures and whether the phosphosites have characterized functions or not. With knowledge of the wide functional range of phosphorylation, we attempted to clarify whether the evolutionary conservation of phosphosites is different among distinct functional modules. We grouped the phosphosites in the human genome into the modules according to the functional categories of KEGG (Kyoto Encyclopedia of Genes and Genomes) and investigated their evolutionary conservation in vertebrate genomes from mouse to zebrafish. We have found that the phosphosites in the vertebrate-specific functional modules (VFMs), such as cellular signaling processes and responses to stimuli, are evolutionarily more conserved than those in the basic functional modules (BFMs), such as metabolic and genetic processes. The phosphosites in the VFMs are also significantly more conserved than their flanking regions, whereas those in the BFMs are not. These results hold for both serine/threonine and tyrosine residues, although the fraction of phosphorylated tyrosine residues is increased in the VFMs. Moreover, the difference in the evolutionary conservation of the phosphosites between the VFMs and BFMs could not be explained by the difference in the local protein structures. There is also a higher fraction of phosphosites with known functions in the VFMs than BFMs. Based on these findings, we have concluded that protein phosphorylation may play more dominant roles for the VFMs than BFMs during the vertebrate evolution. As phosphorylation is a quite rapid biological reaction, the VFMs that quickly respond to outer stimuli and inner signals might heavily depend on this regulatory mechanism. Our results imply that phosphorylation may have an essential role in the evolution of vertebrates.     

CPhos: A program to calculate and visualize evolutionarily conserved functional phosphorylation sites

Profiling using high‐throughput MS has discovered an overwhelming number of novel protein phosphorylation sites (“phosphosites”). However, the functional relevance of these sites is not always clear. In light of recent studies on the evolutionary mechanism of phosphorylation, we have developed CPhos, a Java program that can assess the conservation of phosphosites among species using an information theory‐based approach. The degree of conservation established using CPhos can be used to assess the functional significance of phosphosites. CPhos has a user friendly graphical user interface and is available both as a web service and as a standalone Java application to assist phosphoproteomic researchers in analyzing and prioritizing lists of phosphosites for further experimental validation. CPhos can be accessed or downloaded at http://helixweb.nih.gov/CPhos/ .     

Quantitative proteomics indicate a strong correlation of mitotic phospho-/dephosphorylation with non-structured regions of substrates

Protein phosphorylation plays a critical role in the regulation and progression of mitosis. >40,000 phosphorylated residues and the associated kinases have been identified to date via proteomic analyses. Although some of these phosphosites are associated with regulation of either protein-protein interactions or the catalytic activity of the substrate protein, the roles of most mitotic phosphosites remain unclear. In this study, we examined structural properties of mitotic phosphosites and neighboring residues to understand the role of heavy phosphorylation in non-structured domains. Quantitative mass spectrometry analysis of mitosis-arrested and non-arrested HeLa cells revealed >4100 and > 2200 residues either significantly phosphorylated or dephosphorylated, respectively, at mitotic entry. The calculated disorder scores of amino acid sequences of neighboring individual phosphosites revealed that >70% of dephosphorylated phosphosites exist in disordered regions, whereas 50% of phosphorylated sites exist in non-structured domains. A clear inverse correlation was observed between probability of phosphorylation in non-structured domain and increment of phosphorylation in mitosis. These results indicate that at entry to mitosis, a significant number of phosphate groups are removed from non-structured domains and transferred to more-structured domains. Gene ontology term analysis revealed that mitosis-related proteins are heavily phosphorylated, whereas RNA-related proteins are both dephosphorylated and phosphorylated, suggesting that heavy phosphorylation/dephosphorylation in non-structured domains of RNA-binding proteins plays a role in dynamic rearrangement of RNA-containing organelles, as well as other intracellular environments.     

PhoPepMass: A database and search tool assisting human phosphorylation peptide identification from mass spectrometry data

Protein phosphorylation, one of the most important protein post-translational modifications, is involved in various biological processes, and the identification of phosphorylation peptides (phosphopeptides) and their corresponding phosphorylation sites (phosphosites) will facilitate the understanding of the molecular mechanism and function of phosphorylation. Mass spectrometry (MS) provides a high-throughput technology that enables the identification of large numbers of phosphosites. PhoPepMass is designed to assist human phosphopeptide identification from MS data based on a specific database of phophopeptide masses and a multivariate hypergeometric matching algorithm. It contains 244,915 phosphosites from several public sources. Moreover, the accurate masses of peptides and fragments with phosphosites were calculated. It is the first database that provides a systematic resource for the query of phosphosites on peptides and their corresponding masses. This allows researchers to search certain proteins of which phosphosites have been reported, to browse detailed phosphopeptide and fragment information, to match masses from MS analyses with defined threshold to the corresponding phosphopeptide, and to compare proprietary phosphopeptide discovery results with results from previous studies. Additionally, a database search software is created and a “two-stage search strategy” is suggested to identify phosphopeptides from tandem mass spectra of proteomics data. We expect PhoPepMass to be a useful tool and a source of reference for proteomics researchers. PhoPepMass is available at https://www.scbit.org/phopepmass/index.html.     

Protein abundance is key to distinguish promiscuous from functional phosphorylation based on evolutionary information

In eukaryotic cells, protein phosphorylation is an important and widespread mechanism used to regulate protein function. Yet, of the thousands of phosphosites identified to date, only a few hundred at best have a characterized function. It was recently shown that these functional sites are significantly more conserved than phosphosites of unknown function, stressing the importance of considering evolutionary conservation in assessing the global functional landscape of phosphosites. This leads us to review studies that examined the impact of phosphorylation on evolutionary conservation. While all these studies have shown that conservation is greater among phosphorylated sites compared with non-phosphorylated ones, the magnitude of this difference varies greatly. Further, not all studies have considered key factors that may influence the rate of phosphosite evolution. Such key factors are their localization in ordered or disordered regions, their stoichiometry or the abundance of their corresponding protein. Here we take into account all of these factors simultaneously, which reveals remarkable evolutionary patterns. First, while it is well established that protein conservation increases with abundance, we show that phosphosites partly follow an opposite trend. More precisely, Saccharomyces cerevisiae phosphosites present among abundant proteins are 1.5 times more likely to diverge in the closely related species Saccharomyces bayanus when compared with phosphosites present in the 5 per cent least abundant proteins. Second, we show that conservation is coupled to stoichiometry, whereby sites frequently phosphorylated are more conserved than those rarely phosphorylated. Finally, we provide a model of functional and noisy or 'accidental' phosphorylation that explains these observations.     

Protein kinases phosphorylate long disordered regions in intrinsically disordered proteins

Phosphorylation is a major post‐translational modification that plays a central role in signaling pathways. Protein kinases phosphorylate substrates (phosphoproteins) by adding phosphate at Ser/Thr or Tyr residues (phosphosites). A large amount of data identifying and describing phosphosites in phosphoproteins has been reported but the specificity of phosphorylation is not fully resolved. In this report, data of kinase‐substrate pairs identified by the Kinase‐Interacting Substrate Screening (KISS) method were used to analyze phosphosites in intrinsically disordered regions (IDRs) of intrinsically disordered proteins. We compared phosphorylated and nonphosphorylated IDRs and found that the phosphorylated IDRs were significantly longer than nonphosphorylated IDRs. The phosphorylated IDR is often the longest IDR (71%) in a phosphoprotein when only a single phosphosite exists in the IDR, and when the phosphoprotein has multiple phosphosites in an IDR(s), the phosphosites are primarily localized in a single IDR (78%) and this IDR is usually the longest one (81%). We constructed a stochastic model of phosphorylation to estimate the effect of IDR length. The model that accounted for IDR length produced more realistic results when compared with a model that excluded the IDR length. We propose that the IDR length is a significant determinant for locating kinase phosphorylation sites in phosphoproteins.     

Cooperativity within proximal phosphorylation sites is revealed from large-scale proteomics data

Background  

Phosphorylation is the most prevalent post-translational modification on eukaryotic proteins. Multisite phosphorylation enables a specific combination of phosphosites to determine the speed, specificity and duration of biological response. Until recent years, the lack of high quality data limited the possibility for analyzing the properties of phosphorylation at the proteome scale and in the context of a wide range of conditions. Thanks to advances of mass spectrometry technologies, thousands of phosphosites from in-vivo experiments were identified and archived in the public domain. Such resource is appropriate to derive an unbiased view on the phosphosites properties in eukaryotes and on their functional relevance.     

PHOSIDA (phosphorylation site database): management,structural and evolutionary investigation,and prediction of phosphosites

PHOSIDA , a phosphorylation site database, integrates thousands of high-confidence in vivo phosphosites identified by mass spectrometry-based proteomics in various species. For each phosphosite, PHOSIDA lists matching kinase motifs, predicted secondary structures, conservation patterns, and its dynamic regulation upon stimulus. Using support vector machines, PHOSIDA also predicts phosphosites.     

Phosphorylation network rewiring by gene duplication

Elucidating how complex regulatory networks have assembled during evolution requires a detailed understanding of the evolutionary dynamics that follow gene duplication events, including changes in post‐translational modifications. We compared the phosphorylation profiles of paralogous proteins in the budding yeast Saccharomyces cerevisiae to that of a species that diverged from the budding yeast before the duplication of those genes. We found that 100 million years of post‐duplication divergence are sufficient for the majority of phosphorylation sites to be lost or gained in one paralog or the other, with a strong bias toward losses. However, some losses may be partly compensated for by the evolution of other phosphosites, as paralogous proteins tend to preserve similar numbers of phosphosites over time. We also found that up to 50% of kinase–substrate relationships may have been rewired during this period. Our results suggest that after gene duplication, proteins tend to subfunctionalize at the level of post‐translational regulation and that even when phosphosites are preserved, there is a turnover of the kinases that phosphorylate them.     

Untargeted metabolomics unravels functionalities of phosphorylation sites in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

Background

Coordinated through a complex network of kinases and phosphatases, protein phosphorylation regulates essentially all cellular processes in eukaryotes. Recent advances in proteomics enable detection of thousands of phosphorylation sites (phosphosites) in single experiments. However, functionality of the vast majority of these sites remains unclear and we lack suitable approaches to evaluate functional relevance at a pace that matches their detection.

Results

Here, we assess functionality of 26 phosphosites by introducing phosphodeletion and phosphomimic mutations in 25 metabolic enzymes and regulators from the TOR and HOG signaling pathway in Saccharomyces cerevisiae by phenotypic analysis and untargeted metabolomics. We show that metabolomics largely outperforms growth analysis and recovers 10 out of the 13 previously characterized phosphosites and suggests functionality for several novel sites, including S79 on the TOR regulatory protein Tip41. We analyze metabolic profiles to identify consequences underlying regulatory phosphorylation events and detecting glycerol metabolism to have a so far unknown influence on arginine metabolism via phosphoregulation of the glycerol dehydrogenases. Further, we also find S508 in the MAPKK Pbs2 as a potential link for cross-talking between HOG signaling and the cell wall integrity pathway.

Conclusions

We demonstrate that metabolic profiles can be exploited for gaining insight into regulatory consequences and biological roles of phosphosites. Altogether, untargeted metabolomics is a fast, sensitive and informative approach appropriate for future large-scale functional analyses of phosphosites.

     

The in vivo phosphorylation sites in multiple isoforms of amphiphysin I from rat brain nerve terminals

Amphiphysin I (amphI) is dephosphorylated by calcineurin during nerve terminal depolarization and synaptic vesicle endocytosis (SVE). Some amphI phosphorylation sites (phosphosites) have been identified with in vitro studies or phosphoproteomics screens. We used a multifaceted strategy including 32P tracking to identify all in vivo amphI phosphosites and determine their relative abundance and potential relevance to SVE. AmphI was extracted from 32P-labeled synaptosomes, phosphopeptides were isolated from proteolytic digests using TiO2 chromatography, and mass spectrometry revealed 13 sites: serines 250, 252, 262, 268, 272, 276, 285, 293, 496, 514, 539, and 626 and Thr-310. These were distributed into two clusters around the proline-rich domain and the C-terminal Src homology 3 domain. Hierarchical phosphorylation of Ser-262 preceded phosphorylation of Ser-268, -272, -276, and -285. Off-line HPLC separation and two-dimensional tryptic mapping of 32P-labeled amphI revealed that Thr-310, Ser-293, Ser-285, Ser-272, Ser-276, and Ser-268 contained the highest 32P incorporation and were the most stimulus-sensitive. Individually Thr-310 and Ser-293 were the most abundant phosphosites, incorporating 16 and 23% of the 32P. The multiple phosphopeptides containing Ser-268, Ser-276, Ser-272, and Ser-285 had 27% of the 32P. Evidence for a role for at least one proline-directed protein kinase and one non-proline-directed kinase was obtained. Four phosphosites predicted for non-proline-directed kinases, Ser-626, -250, -252, and -539, contained low amounts of 32P and were not depolarization-responsive. At least one alternatively spliced amphI isoform was identified in synaptosomes as being constitutively phosphorylated because it did not incorporate 32P during the 1-h labeling period. Multiple phosphosites from amphI-co-migrating synaptosomal proteins were also identified, including SGIP (Src homology 3 domain growth factor receptor-bound 2 (Grb2)-like (endophilin)-interacting protein 1), AAK1, eps15R, MAP6, alpha/beta-adducin, and HCN1. The results reveal two sets of amphI phosphosites that are either dynamically turning over or constitutively phosphorylated in nerve terminals and improve understanding of the role of individual amphI sites or phosphosite clusters in synaptic SVE.     

Phosphorylation states greatly regulate the activity and gating properties of Cav3.1 T-type Ca2+ channels

Ca_v3.1 T-type Ca²⁺ channels play pivotal roles in neuronal low-threshold spikes, visceral pain, and pacemaker activity. Phosphorylation has been reported to potently regulate the activity and gating properties of Ca_v3.1 channels. However, systematic identification of phosphorylation sites (phosphosites) in Ca_v3.1 channel has been poorly investigated. In this work, we analyzed rat Ca_v3.1 protein expressed in HEK-293 cells by mass spectrometry, identified 30 phosphosites located at the cytoplasmic regions, and illustrated them as a Ca_v3.1 phosphorylation map which includes the reported mouse Ca_v3.1 phosphosites. Site-directed mutagenesis of the phosphosites to Ala residues and functional analysis of the phospho-silent Ca_v3.1 mutants expressed in Xenopus oocytes showed that the phospho-silent mutation of the N-terminal Ser18 reduced its current amplitude with accelerated current kinetics and negatively shifted channel availability. Remarkably, the phospho-silent mutations of the C-terminal Ser residues (Ser1924, Ser2001, Ser2163, Ser2166, or Ser2189) greatly reduced their current amplitude without altering the voltage-dependent gating properties. In contrast, the phosphomimetic Asp mutations of Ca_v3.1 on the N- and C-terminal Ser residues reversed the effects of the phospho-silent mutations. Collectively, these findings demonstrate that the multiple phosphosites of Ca_v3.1 at the N- and C-terminal regions play crucial roles in the regulation of the channel activity and voltage-dependent gating properties.     

A systematic comparative and structural analysis of protein phosphorylation sites based on the mtcPTM database

mtcPTM is an online repository of human and mouse phosphosites in which data are hierarchically organized to preserve biologically relevant experimental information, thus allowing straightforward comparisons of phosphorylation patterns found under different conditions. The database also contains the largest available collection of atomic models of phosphorylatable proteins. Detailed analysis of this structural dataset reveals that phosphorylation sites are found in a heterogeneous range of structural and sequence contexts. mtcPTM is available on the web .     

Novel Binding Partners and Differentially Regulated Phosphorylation Sites Clarify Eps8 as a Multi-Functional Adaptor

Eps8 is involved in both cell signalling and receptor trafficking. It is a known phosphorylation substrate for two proteins involved in the fibroblast growth factor receptor (FGFR) signalling pathway: the receptor itself and Src. Here we report a differential proteomic analysis of Eps8 aimed to identify specific FGFR and Src family kinase dependent phosphosites and co-associated phosphodependent binding partners. This study reveals a total of 22 Eps8 pTyr and pSer/Thr phosphorylation sites, including those that are dependent on Src family and FGFR kinase activity. Peptide affinity purification of proteins that bind to a selection of the pTyr phosphosites has identified a range of novel Eps8 binding partners including members of the intracellular vesicle trafficking machinery (clathrin and AP-2), proteins which have been shown to regulate activated receptor trafficking (NBR1 and Vav2), and proteins involved in receptor signalling (IRS4 and Shp2). Collectively this study significantly extends the understanding of Eps8 post-translational modification by regulated phosphorylation, identifies novel Eps8 binding partners implicated in receptor trafficking and signalling, and confirms the functions of Eps8 at the nexus of receptor signalling and vesicular trafficking.     

Identification of the PLK2-Dependent Phosphopeptidome by Quantitative Proteomics

Polo-like kinase 2 (PLK2) has been recently recognized as the major enzyme responsible for phosphorylation of α-synuclein at S129 in vitro and in vivo, suggesting that this kinase may play a key role in the pathogenesis of Parkinson''s disease and other synucleinopathies. Moreover PLK2 seems to be implicated in cell division, oncogenesis, and synaptic regulation of the brain. However little is known about the phosphoproteome generated by PLK2 and, consequently the overall impact of PLK2 on cellular signaling. To fill this gap we exploited an approach based on in vitro kinase assay and quantitative phosphoproteomics. A proteome-derived peptide library obtained by digestion of undifferentiated human neuroblastoma cell line was exhaustively dephosphorylated by lambda phosphatase followed by incubation with or without PLK2 recombinant kinase. Stable isotope labeling based quantitative phosphoproteomics was applied to identify the phosphosites generated by PLK2. A total of 98 unique PLK2-dependent phosphosites from 89 proteins were identified by LC-MS/MS. Analysis of the primary structure of the identified phosphosites allowed the detailed definition of the kinase specificity and the compilation of a list of potential PLK2 targets among those retrieved in PhosphositePlus, a curated database of in cell/vivo phosphorylation sites.     

Spatial features of proteins related to their phosphorylation and associated structural changes

Protein phosphorylation is widely used in biological regulatory processes. The study of spatial features related to phosphorylation sites is necessary to increase the efficacy of recognition of phosphorylation patterns in protein sequences. Using the data on phosphosites found in amino acid sequences, we mapped these sites onto 3D structures and studied the structural variability of the same sites in different PDB entries related to the same proteins. Solvent accessibility was calculated for the residues known to be phosphorylated. A significant change in accessibility was shown for many sites, but several ones were determined as buried in all the structures considered. Most phosphosites were found in coil regions. However, a significant portion was located in the structurally stable ordered regions. Comparison of structures with the same sites in modified and unmodified states showed that the region surrounding a site could be significantly shifted due to phosphorylation. Comparison between non‐modified structures (as well as between the modified ones) suggested that phosphorylation stabilizes one of the possible conformations. The local structure around the site could be changed due to phosphorylation, but often the initial conformation of the site surrounding is not altered within bounds of a rather large substructure. In this case, we can observe an extensive displacement within a protein domain. Phosphorylation without structural alteration seems to provide the interface for domain‐domain or protein‐protein interactions. Accounting for structural features is important for revealing more specific patterns of phosphorylation. It is also necessary for explaining structural changes as a basis for regulatory processes.