Integrative features of the yeast phosphoproteome and protein-protein interaction map

Following recent advances in high-throughput mass spectrometry (MS)-based proteomics, the numbers of identified phosphoproteins and their phosphosites have greatly increased in a wide variety of organisms. Although a critical role of phosphorylation is control of protein signaling, our understanding of the phosphoproteome remains limited. Here, we report unexpected, large-scale connections revealed between the phosphoproteome and protein interactome by integrative data-mining of yeast multi-omics data. First, new phosphoproteome data on yeast cells were obtained by MS-based proteomics and unified with publicly available yeast phosphoproteome data. This revealed that nearly 60% of ～6,000 yeast genes encode phosphoproteins. We mapped these unified phosphoproteome data on a yeast protein-protein interaction (PPI) network with other yeast multi-omics datasets containing information about proteome abundance, proteome disorders, literature-derived signaling reactomes, and in vitro substratomes of kinases. In the phospho-PPI, phosphoproteins had more interacting partners than nonphosphoproteins, implying that a large fraction of intracellular protein interaction patterns (including those of protein complex formation) is affected by reversible and alternative phosphorylation reactions. Although highly abundant or unstructured proteins have a high chance of both interacting with other proteins and being phosphorylated within cells, the difference between the number counts of interacting partners of phosphoproteins and nonphosphoproteins was significant independently of protein abundance and disorder level. Moreover, analysis of the phospho-PPI and yeast signaling reactome data suggested that co-phosphorylation of interacting proteins by single kinases is common within cells. These multi-omics analyses illuminate how wide-ranging intracellular phosphorylation events and the diversity of physical protein interactions are largely affected by each other.     

Phosphoprotein profiling by PA-GeLC-MS/MS

A significant consequence of protein phosphorylation is to alter protein-protein interactions, leading to dynamic regulation of the components of protein complexes that direct many core biological processes. Recent proteomic studies have populated databases with extensive compilations of cellular phosphoproteins and phosphorylation sites and a similarly deep coverage of the subunit compositions and interactions in multiprotein complexes. However, considerably less data are available on the dynamics of phosphorylation, composition of multiprotein complexes or that define their interdependence. We describe a method to identify candidate phosphoprotein complexes by combining phosphoprotein affinity chromatography, separation by size, denaturing gel electrophoresis, protein identification by tandem mass spectrometry, and informatics analysis. Toward developing phosphoproteome profiling, we have isolated native phosphoproteins using a phosphoprotein affinity matrix, Pro-Q Diamond resin (Molecular Probes-Invitrogen). This resin quantitatively retains phosphoproteins and associated proteins from cell extracts. Pro-Q Diamond purification of a yeast whole cell extract followed by 1-D PAGE separation, proteolysis and ESI LC-MS/MS, a method we term PA-GeLC-MS/MS, yielded 108 proteins, a majority of which were known phosphoproteins. To identify proteins that were purified as parts of phosphoprotein complexes, the Pro-Q eluate was separated into two fractions by size, <100 kDa and >100 kDa, before analysis by PAGE and ESI LC-MS/MS and the component proteins queried against databases to identify protein-protein interactions. The <100 kDa fraction was enriched in phosphoproteins indicating the presence of monomeric phosphoproteins. The >100 kDa fraction contained 171 proteins of 20-80 kDa, nearly all of which participate in known protein-protein interactions. Of these 171, few are known phosphoproteins, consistent with their purification by participation in protein complexes. By comparing the results of our phosphoprotein profiling with the informational databases on phosphoproteomics, protein-protein interactions and protein complexes, we have developed an approach to examining the correlation between protein interactions and protein phosphorylation.     

Proteomic analysis of in vivo phosphorylated synaptic proteins

In the nervous system, protein phosphorylation is an essential feature of synaptic function. Although protein phosphorylation is known to be important for many synaptic processes and in disease, little is known about global phosphorylation of synaptic proteins. Heterogeneity and low abundance make protein phosphorylation analysis difficult, particularly for mammalian tissue samples. Using a new approach, combining both protein and peptide immobilized metal affinity chromatography and mass spectrometry data acquisition strategies, we have produced the first large scale map of the mouse synapse phosphoproteome. We report over 650 phosphorylation events corresponding to 331 sites (289 have been unambiguously assigned), 92% of which are novel. These represent 79 proteins, half of which are novel phosphoproteins, and include several highly phosphorylated proteins such as MAP1B (33 sites) and Bassoon (30 sites). An additional 149 candidate phosphoproteins were identified by profiling the composition of the protein immobilized metal affinity chromatography enrichment. All major synaptic protein classes were observed, including components of important pre- and postsynaptic complexes as well as low abundance signaling proteins. Bioinformatic and in vitro phosphorylation assays of peptide arrays suggest that a small number of kinases phosphorylate many proteins and that each substrate is phosphorylated by many kinases. These data substantially increase existing knowledge of synapse protein phosphorylation and support a model where the synapse phosphoproteome is functionally organized into a highly interconnected signaling network.     

Analyzing yeast protein-protein interaction data obtained from different sources

High-throughput methods for detecting protein interactions, such as mass spectrometry and yeast two-hybrid assays, continue to produce vast amounts of data that may be exploited to infer protein function and regulation. As this article went to press, the pool of all published interaction information on Saccharomyces cerevisiae was 15,143 interactions among 4,825 proteins, and power-law scaling supports an estimate of 20,000 specific protein interactions. To investigate the biases, overlaps, and complementarities among these data, we have carried out an analysis of two high-throughput mass spectrometry (HMS)-based protein interaction data sets from budding yeast, comparing them to each other and to other interaction data sets. Our analysis reveals 198 interactions among 222 proteins common to both data sets, many of which reflect large multiprotein complexes. It also indicates that a "spoke" model that directly pairs bait proteins with associated proteins is roughly threefold more accurate than a "matrix" model that connects all proteins. In addition, we identify a large, previously unsuspected nucleolar complex of 148 proteins, including 39 proteins of unknown function. Our results indicate that existing large-scale protein interaction data sets are nonsaturating and that integrating many different experimental data sets yields a clearer biological view than any single method alone.     

Analysis of desiccation-induced candidate phosphoproteins from Craterostigma plantagineum isolated with a modified metal oxide affinity chromatography procedure

Reversible protein phosphorylation/dephosphorylation is crucial for regulation of many cellular events, and increasing evidence indicates that this post-translational modification is also involved in the complex process of acquisition of desiccation tolerance. To analyze the phosphoproteome of the desiccation tolerant resurrection plant Craterostigma plantagineum, MOAC-enriched proteins from leaves at different stages of a de-/rehydration cycle were separated by 2-D PAGE and detected by phosphoprotein-specific staining. Using this strategy 20 putative phosphoproteins were identified by MALDI-TOF MS and MS/MS, which were not detected when total proteins were analyzed. The characterized desiccation-related phosphoproteins CDeT11-24 and CDeT6-19 were used as internal markers to validate the specificity of the analyses. For 16 of the identified proteins published evidence suggests that they are phosphoproteins. Comparative analysis of the 2-D gels showed that spot intensities of most identified putative phosphoproteins change during the de-/rehydration cycle. This suggests an involvement of these proteins in desiccation tolerance. Nearly all changes in the phosphoproteome of C. plantagineum, which are triggered by dehydration, are reversed within 4 days of rehydration, which is in agreement with physiological observations. Possible functions of selected proteins are discussed in the context of the de-/rehydration cycle.     

Evolution of networks and sequences in eukaryotic cell cycle control

The molecular networks regulating the G1-S transition in budding yeast and mammals are strikingly similar in network structure. However, many of the individual proteins performing similar network roles appear to have unrelated amino acid sequences, suggesting either extremely rapid sequence evolution, or true polyphyly of proteins carrying out identical network roles. A yeast/mammal comparison suggests that network topology, and its associated dynamic properties, rather than regulatory proteins themselves may be the most important elements conserved through evolution. However, recent deep phylogenetic studies show that fungal and animal lineages are relatively closely related in the opisthokont branch of eukaryotes. The presence in plants of cell cycle regulators such as Rb, E2F and cyclins A and D, that appear lost in yeast, suggests cell cycle control in the last common ancestor of the eukaryotes was implemented with this set of regulatory proteins. Forward genetics in non-opisthokonts, such as plants or their green algal relatives, will provide direct information on cell cycle control in these organisms, and may elucidate the potentially more complex cell cycle control network of the last common eukaryotic ancestor.     

Profiling phosphoproteins of yeast mitochondria reveals a role of phosphorylation in assembly of the ATP synthase

Mitochondria are crucial for numerous cellular processes, yet the regulation of mitochondrial functions is only understood in part. Recent studies indicated that the number of mitochondrial phosphoproteins is higher than expected; however, the effect of reversible phosphorylation on mitochondrial structure and function has only been defined in a few cases. It is thus crucial to determine authentic protein phosphorylation sites from highly purified mitochondria in a genetically tractable organism. The yeast Saccharomyces cerevisiae is a major model organism for the analysis of mitochondrial functions. We isolated highly pure yeast mitochondria and performed a systematic analysis of phosphorylation sites by a combination of different enrichment strategies and mass spectrometry. We identified 80 phosphorylation sites in 48 different proteins. These mitochondrial phosphoproteins are involved in critical mitochondrial functions, including energy metabolism, protein biogenesis, fatty acid metabolism, metabolite transport, and redox regulation. By combining yeast genetics and in vitro biochemical analysis, we found that phosphorylation of a serine residue in subunit g (Atp20) regulates dimerization of the mitochondrial ATP synthase. The authentic phosphoproteome of yeast mitochondria will represent a rich source to uncover novel roles of reversible protein phosphorylation.     

Phosphoproteomic analysis of the mouse brain cytosol reveals a predominance of protein phosphorylation in regions of intrinsic sequence disorder

We analyzed the mouse forebrain cytosolic phosphoproteome using sequential (protein and peptide) IMAC purifications, enzymatic dephosphorylation, and targeted tandem mass spectrometry analysis strategies. In total, using complementary phosphoenrichment and LC-MS/MS strategies, 512 phosphorylation sites on 540 non-redundant phosphopeptides from 162 cytosolic phosphoproteins were characterized. Analysis of protein domains and amino acid sequence composition of this data set of cytosolic phosphoproteins revealed that it is significantly enriched in intrinsic sequence disorder, and this enrichment is associated with both cellular location and phosphorylation status. The majority of phosphorylation sites found by MS were located outside of structural protein domains (97%) but were mostly located in regions of intrinsic sequence disorder (86%). 368 phosphorylation sites were located in long regions of disorder (over 40 amino acids long), and 94% of proteins contained at least one such long region of disorder. In addition, we found that 58 phosphorylation sites in this data set occur in 14-3-3 binding consensus motifs, linear motifs that are associated with unstructured regions in proteins. These results demonstrate that in this data set protein phosphorylation is significantly depleted in protein domains and significantly enriched in disordered protein sequences and that enrichment of intrinsic sequence disorder may be a common feature of phosphoproteomes. This supports the hypothesis that disordered regions in proteins allow kinases, phosphatases, and phosphorylation-dependent binding proteins to gain access to target sequences to regulate local protein conformation and activity.     

Motif-specific sampling of phosphoproteomes

Phosphoproteomics, the targeted study of a subfraction of the proteome which is modified by phosphorylation, has become an indispensable tool to study cell signaling dynamics. We described a methodology that linked phosphoproteome and proteome analysis based on Ba2+ binding properties of amino acids. This technology selected motif-specific phosphopeptides independent of the system under analysis. MudPIT (Multidimensional Identification Technology) identified 1037 precipitated phosphopeptides from as little as 250 microg of proteins. To extend coverage of the phosphoproteome, we sampled the nuclear extract of HeLa cells with three values of Ba2+ ions molarity. The presence of more than 70% of identified phosphoproteins was further substantiated by their nonmodified peptides. Upon isoproterenol stimulation of HEK cells, we identified an increasing number of phosphoproteins from MAPK cascades and AKAP signaling hubs. We quantified changes in both protein and phosphorylation levels of 197 phosphoproteins including a critical kinase, MAPK1. Integration of differential phosphorylation of MAPK1 with knowledge bases constructed modules that correlated well with its role as node in cross-talk of canonical pathways.     

Phosphoproteome mapping of cardiomyocyte mitochondria in a rat model of heart failure

Mitochondria are complex organelles essential to cardiomyocyte survival. Protein phosphorylation is emerging as a key regulator of mitochondrial function. In the study reported here, we analyzed subsarcolemmal (SSM) mitochondria harvested from rats who have received 4 weeks of aldosterone/salt treatment to simulate the neurohormonal profile of human congestive heart failure. Our objective was to obtain an initial qualitative inventory of the phosphoproteins in this biologic system. SSM mitochondria were harvested, and the phosphoproteome was analyzed with a gel-free bioanalytical platform. Mitochondrial proteins were digested with trypsin, and the digests were enriched for phosphopeptides with immobilized metal ion affinity chromatography. The phosphopeptides were analyzed by ion trap liquid chromatography–tandem mass spectrometry, and the phosphoproteins identified via database searches. Based on MS/MS and MS³ data, we characterized a set of 42 phosphopeptides that encompassed 39 phosphorylation sites. These peptides mapped to 26 proteins, for example, long-chain specific acyl-CoA dehydrogenase, Complex III subunit 6, and mitochondrial import receptor TOM70. Collectively, the characterized phosphoproteins belong to diverse functional modules, including bioenergetic pathways, protein import machinery, and calcium handling. The phosphoprotein panel discovered in this study provides a foundation for future differential phosphoproteome profiling toward an integrated understanding of the role of mitochondrial phosphorylation in heart failure.     

A brief history of TOR

The TOR (target of rapamycin) serine/threonine kinases are fascinating in that they influence many different aspects of eukaryote physiology including processes often dysregulated in disease. Beginning with the initial characterization of rapamycin as an antifungal agent, studies with yeast have contributed greatly to our understanding of the molecular pathways in which TORs operate. Recently, building on advances in quantitative MS, the rapamycin-dependent phosphoproteome in the budding yeast Saccharomyces cerevisiae was elucidated. These studies emphasize the central importance of TOR and highlight its many previously unrecognized functions. One of these, the regulation of intermediary metabolism, is discussed.     

First the CDKs, now the DDKs.

In budding yeast, Dbf4p and Cdc7p control initiation of DNA synthesis. They form a protein kinase - Cdc7p being the catalytic subunit and Dbf4p a cyclin-like molecule that activates the kinase in late G1 phase. Dbf4p also targets Cdc7p to origins of replication, where probable substrates include certain Mcm proteins. Recent studies have identified Dbf4p- and Cdc7p-related proteins in fission yeast and metazoans. These homologues also phosphorylate Mcm proteins and could have a similar function to that of Dbf4p-Cdc7p in budding yeast. Thus, it seems likely that, like the cyclin-dependent kinases (CDKs), the Dbf4p-Cdc7p activity is conserved in all eukaryotes.     

Analysis of the splicing machinery in fission yeast: a comparison with budding yeast and mammals

Based on genetic and bioinformatic analysis, 80 proteins from the newly sequenced Schizosaccharomyces pombe genome appear to be splicing factors. The fission yeast splicing factors were compared to those of Homo sapiens and Saccharomyces cerevisiae in order to determine the extent of conservation or divergence that has occurred over the billion years of evolution that separate these organisms. Our results indicate that many of the factors present in all three organisms have been well conserved throughout evolution. It is calculated that 38% of the fission yeast splicing factors are more similar to the human proteins than to the budding yeast proteins (>10% more similar or similar over a greater region). Many of the factors in this category are required for recognition of the 3′ splice site. Ten fission yeast splicing factors, including putative regulatory factors, have human homologs, but no apparent budding yeast homologs based on sequence data alone. Many of the budding yeast factors that are absent in fission yeast are associated with the U1 and U4/U6.U5 snRNP. Collectively the data presented in this survey indicate that of the two yeasts, S.pombe contains a splicing machinery more closely reflecting the archetype of a spliceosome.     

Modular complexes that regulate actin assembly in budding yeast

The actin cytoskeleton of budding yeast contains an extensive set of actin-associated proteins with conserved mammalian counterparts. For more than 20 years, yeast has been used as a model organism to dissect the in vivo functions of these factors, revealing an intricate web of genetic interactions in the cell. Now, a surge of biochemical reports is defining the physical interactions and activities of these proteins and providing mechanistic insights into their cellular roles. The emerging view is that most actin-associated proteins do not act alone but, rather, associate to form modular protein complexes that regulate actin assembly and organization.     

Topological structure analysis of the protein-protein interaction network in budding yeast

Interaction detection methods have led to the discovery of thousands of interactions between proteins, and discerning relevance within large-scale data sets is important to present-day biology. Here, a spectral method derived from graph theory was introduced to uncover hidden topological structures (i.e. quasi-cliques and quasi-bipartites) of complicated protein-protein interaction networks. Our analyses suggest that these hidden topological structures consist of biologically relevant functional groups. This result motivates a new method to predict the function of uncharacterized proteins based on the classification of known proteins within topological structures. Using this spectral analysis method, 48 quasi-cliques and six quasi-bipartites were isolated from a network involving 11,855 interactions among 2617 proteins in budding yeast, and 76 uncharacterized proteins were assigned functions.     

Phosphoproteome profile of human lung cancer cell line A549

As an in vitro model for type II human lung cancer, A549 cells resist cytotoxicity via phosphorylation of proteins as demonstrated by many studies. However, to date, no large-scale phosphoproteome investigation has been conducted on A549. Here, we performed a systematical analysis of the phosphoproteome of A549 by using mass spectrometry (MS)-based strategies. This investigation led to the identification of 337 phosphorylation sites on 181 phosphoproteins. Among them, 67 phosphoproteins and 230 phosphorylation sites identified appeared to be novel with no previous characterization in lung cancer. Based on their known functions as reported in the literature, these phosphoproteins were functionally organized into highly interconnected networks. Western blotting and immunohistochemistry analyses were performed to validate the expression of a bottleneck phosphoprotein YAP1 in cancer cell lines and tissues. This dataset provides a valuable resource for further studies on phosphorylation and lung carcinogenesis.     

Enrichment techniques employed in phosphoproteomics

Rapid changes of protein phosphorylation play a crucial role in the regulation of many cellular processes. Being post-translationally modified, phosphoproteins are often present in quite low abundance and tend to co-exist with their unphosphorylated isoform within the cell. To make their identification more practicable, the use of enrichment protocols is often required. The enrichment strategies can be performed either at the level of phosphoproteins or at the level of phosphopeptides. Both approaches have their advantages and disadvantages. Most enriching strategies are based on chemical modifications, affinity chromatography to capture peptides and proteins containing negatively charged phosphate groups onto a positively charged matrix, or immunoprecipitation by phospho-specific antibodies. In this article, the most up-to-date enrichment techniques are discussed, taking into account their optimization, and highlighting their advantages and disadvantages. Moreover, these methods are compared to each other, revealing their complementary nature in providing comprehensive coverage of the phosphoproteome.     

Phosphoproteome analysis of the human Chang liver cells using SCX and a complementary mass spectrometric strategy

The liver is the largest organ in the body, with many complex, essential functions, such as metabolism, deintoxication, and secretion, often regulated via post-translational modifications, especially phosphorylation. Thus, the detection of phosphoproteins and phosphorylation sites is important to comprehensively explore human liver biological function. The human Chang liver cell line is among the first derived from non-malignant tissue, and its phosphoproteome profile has never been globally analyzed. To develop the complete phosphoproteome and probe the roles of protein phosphorylation in normal human liver, we adopted a shotgun strategy based on strong cation exchange chromatograph, titanium dioxide and LC-MS/MS to isolate and identify phosphorylated proteins. Two types of MS approach, Q-TOF and IT, were used and compared to identify phosphosites from complex protein mixtures of these cells. A total of 1035 phosphorylation sites and 686 phosphorylated peptides were identified from 607 phosphoproteins. A search using the public database of PhosphoSite showed that approximately 344 phosphoproteins and 760 phosphorylation sites appeared to be novel. In addition, N-terminal phosphorylated peptides were a greater fraction of all identified phosphopeptides. With GOfact analysis, we found that most of the identified phosphoproteins are involved in regulating metabolism, consistent with the liver's role as a key metabolic organ.     

Computational framework for analysis of prey-prey associations in interaction proteomics identifies novel human protein-protein interactions and networks

Large-scale protein-protein interaction data sets have been generated for several species including yeast and human and have enabled the identification, quantification, and prediction of cellular molecular networks. Affinity purification-mass spectrometry (AP-MS) is the preeminent methodology for large-scale analysis of protein complexes, performed by immunopurifying a specific "bait" protein and its associated "prey" proteins. The analysis and interpretation of AP-MS data sets is, however, not straightforward. In addition, although yeast AP-MS data sets are relatively comprehensive, current human AP-MS data sets only sparsely cover the human interactome. Here we develop a framework for analysis of AP-MS data sets that addresses the issues of noise, missing data, and sparsity of coverage in the context of a current, real world human AP-MS data set. Our goal is to extend and increase the density of the known human interactome by integrating bait-prey and cocomplexed preys (prey-prey associations) into networks. Our framework incorporates a score for each identified protein, as well as elements of signal processing to improve the confidence of identified protein-protein interactions. We identify many protein networks enriched in known biological processes and functions. In addition, we show that integrated bait-prey and prey-prey interactions can be used to refine network topology and extend known protein networks.