Integrative Features of the Yeast Phosphoproteome and Protein–Protein Interaction Map
Protein phosphorylation is a reversible , ubiquitous , and fundamentally post-translational modification ( PTM ) that regulates a variety of biological processes; one of its critical roles is the control of protein signaling 1–3 ., Recent advances in mass-spectrometry ( MS ) –based technologies and phosphopeptide enrichment methods have enabled the use of high-throughput in vivo phosphosite mapping 4–7 to identify thousands of phosphoproteins ., To date , around 10 , 000 phosphosites of serine , threonine , or tyrosine residues have been identified in each of many organisms , including human 8–12 , mouse 13 and yeast 14–16 ., Many public databases , such as PHOSIDA 17 , Phospho . ELM 18 , and UniProt 19 , have been developed or expanded to catalog such phosphoproteome data ., Accordingly , the numbers of phosphoproteins that have been identified in various organisms now greatly exceed the numbers known to have roles in protein signaling ., This has raised the question of whether this intracellular phosphorylation , which occurs on such a large scale , has other major roles ., In modern biology , the use of high-throughput screening methods has enabled rapid progress in the disclosure of protein–protein interaction ( PPI ) networks in many organisms 20–27 ., Topological features common to PPI networks ( e . g . , scale-free and small-world properties ) are of prime importance in interpreting intracellular protein behavior and the evolutionary aspects of PPIs 28–31 ., PTM changes the physical characteristics of proteins ., It is therefore probable that reversible PTM has large effects on the dynamic states of intracellular protein-binding patterns and complex formation , and that it controls not only signal transduction but also many other cellular pathways ., However , the impact of PTM on the whole picture of the PPI network has not yet been described ., Here , we describe the intracellular global relationships between protein phosphorylation and physical PPI , as derived from the results of integrative and systematic data-mining of Saccharomyces cerevisiae multi-omics data ( Fig . 1 ) ., New phosphoproteome data on S . cerevisiae were initially obtained by MS–based analysis and unified with data on previously identified phosphoproteomes ., We superimposed the unified phosphoproteome data onto a S . cerevisiae PPI network with other multi-omics data on S . cerevisiae ., From the results , we infer that the tremendous numbers of phosphorylations within a cell have a large impact on PPI diversity , and that intracellular phosphorylation patterns are affected partly by simultaneous phosphorylation of physically bound proteins that is triggered by the action of single kinases ., On the basis of liquid chromatography ( LC ) -MS analysis , we initially identified 1 , 993 S . cerevisiae phosphoproteins containing 6 , 510 phosphosites ., Information on the identified phosphopeptides has been stored in PepBase ( http://pepbase . iab . keio . ac . jp ) ., We unified these new phosphoproteome data with the publicly available phosphoproteome datasets of Holt et al . 16 and UniProt 19 and obtained a total of 3 , 477 phosphoproteins containing 25 , 997 phosphosites ( Fig . 2; Supplementary Table S1 ) ., The pS/pT/pY ratios of this study , the study of Holt et al . , and UniProt were 72%/23%/5% , 72%/23%/5% , and 80%/18%/2% , respectively ., Among the unified phosphoproteome data , 343 phosphoproteins and 2 , 778 phosphosites were not found in the data of Holt et al . or UniProt ., Comparison with S . cerevisiae genomic information 32 revealed that 58 . 5% of the 5 , 815 known and predicted genes were phosphoprotein-encoding genes ( Supplementary Table S2 ) ., Although the use of current high-throughput technologies cannot disclose the entire phosphoproteome picture of a cell , these results imply that most intracellular proteins can be phosphorylated under the appropriate environmental conditions ., The unified phosphoproteome data were superimposed onto the PPI network to generate a “phospho-PPI” network ., PPI data were obtained via DIP ( Database of Interacting Proteins ) 33 and grouped into four categories according to the experimental method used for the PPI assay: all kinds of experimental methods ( “ALL” ) , yeast two-hybrid ( “Y2H” ) , co-immunoprecipitation ( “IMM” ) , and tandem affinity purification ( “TAP” ) ., Among all the protein nodes involved in every category of the phospho-PPI network , the proportion of phosphoproteins was also nearly 60% ( Supplementary Fig . S1 ) ., For example , the phospho-PPI network of the “ALL” category was composed of 4 , 945 proteins , including 2 , 934 phosphoproteins ( 59 . 3% ) and 17 , 215 physical interactions ., To explore specific characteristics of the phospho-PPI network , the number counts of interacting partners of phosphoproteins and nonphosphoproteins were analyzed ( note that throughout this study , the word “nonphosphoprotein” means a protein with no phosphosite identified to date ) ., We found that , in general , phosphoproteins had more interacting partners than nonphosphoproteins ., In each phospho-PPI network of the “ALL” and “Y2H” categories with enough protein nodes for the subsequent statistical analysis , the cumulative percentage distributions of node degrees ( or the number count of interacting partners ) of phosphoproteins and nonphosphoproteins were markedly different ( Fig . 3A and D ) ., For example , in the dataset of “ALL” , 47 . 6% of nonphosphoproteins had three or more interacting partners , but this was true for 67 . 9% of phosphoproteins ., Moreover , in both datasets , about twice as many phosphoproteins as nonphosphoproteins had 10 interacting partners ( Fig . 3B and E ) ., To analyze the statistical significance of this difference in the context of phosphorylation , we prepared randomly generated phospho-PPI networks by “node label shuffling” ( NLS ) , in which the node positions of phosphoproteins and nonphosphoproteins were randomly moved within the phospho-PPI networks ( for details , see Materials and Methods ) ., This demonstrated that the node degree of phosphoproteins was significantly higher than expected from a random distribution ( Fig . 3C and F ) ., Node degree in PPI networks has an exponential relationship with protein expression level 34–36 , perhaps because cellular proteins with more copies have a greater possibility of interacting with others by chance 36 ., Therefore , if the phosphoproteome data are biased by protein abundance and highly abundant proteins tend to be identified as phosphoproteins , there is a strong possibility that the relationship between phosphorylation and node degree is spurious , with no direct causal connection ., In fact , proteome abundance data obtained through a single-cell proteomic analysis combining high-throughput flow cytometry and a library of GFP-tagged yeast strains 37 showed that the number of phosphoproteins in the “ALL” phospho-PPI was skewed , especially among highly abundant proteins ( Fig . 4A and D ) ., However , we demonstrated that in the “ALL” phospho-PPI network there were still significant differences in the node degree levels of phosphoproteins and nonphosphoproteins of similar abundance , and that the differences could be explained independently of protein copy number ( Fig . 4B , C , E and F ) ., Similar results were derived from the phospho-PPI network generated only from the “Y2H” category ( Supplementary Fig . S2 ) ., We further compared the abilities to predict phosphoproteins by using node degree and protein abundance levels above given thresholds ., The predictive power of node degree was markedly higher than that of protein abundance , except in the case of proteins that were extremely abundant ( Supplementary Fig . S3 ) ., If this higher predictive ability were attributable to a spurious relationship associated with the actual intracellular proteome abundance , then the node degree of a protein given by PPI assays would appear to provide a better approximation of the intracellular protein copy number than would single-cell proteomic analysis , which is unlikely ., Protein disorder is also a typical feature of “hub” proteins in PPI networks 38–40 ., Parts of unstructured proteins lack fixed structure , and such disordered regions may have the ability to bind multiple proteins and to diversify PPI networks 38–40 ., Additionally , at the proteome level , phosphorylation occurs at high rates in the disordered regions of proteins 16 , 17 , 41–44 ., Therefore , it is highly likely that protein disorder affects the node degree difference between phosphoproteins and nonphosphoproteins ., For every S . cerevisiae protein registered in UniProt , we calculated the probability of harboring intrinsic disordered regions ( see Materials and Methods ) ., In the “ALL” phospho-PPI network , the ratio of phosphoproteins to nonphosphoproteins increased smoothly with increasing disorder probability level ( Fig . 4G ) ., However , in the same network , the node degree levels of phosphoproteins and nonphosphoproteins of the same disorder probability level were significantly different ( Fig . 4H and I ) ., Even between phosphoproteins that had a low disorder probability of <0 . 1 and nonphosphoproteins that had an extremely high disorder probability of >0 . 9 , the node degree level of the phosphoproteins was significantly higher than that of the nonphosphoproteins ( P\u200a=\u200a0 . 0043 ) ., Similar results were observed in the “Y2H” dataset ( Supplementary Fig . S2 ) ., These results imply that the higher node degree of phosphoproteins than of nonphosphoproteins is at least partly independent of the PPI network diversity produced by unstructured proteins ., Other factors that could influence the relationship between protein phosphorylation and interaction are protein size and protein groups with identical cellular function ., Larger proteins may have a greater chance of being phosphorylated and may provide more binding domains for interactions with other proteins ., However , similar to the results for protein abundance and disorder , statistical significance of the higher node degree of phosphoproteins was observed independently of protein length ( Supplementary Fig . S4 ) ., ( Phosphorylation probability was highly correlated with protein length; Supplementary Fig . S4 . ), In the event that both protein phosphorylation and interaction events occurring in a fraction of proteins confer a particular , identical cellular function , then the global difference in node degree levels of phosphoproteins and nonphosphoproteins would appear to be caused only by differences in function ., However , we found that , for most functional annotations of S . cerevisiae in GO Slim ( a higher level view of Gene Ontology ) , there was a higher node degree level for phosphoproteins than for nonphosphoproteins ( Supplementary Fig . S5 ) ., The average node degree of phosphoproteins is higher than that of nonphosphoproteins 45 , but it was unclear, 1 ) whether this characteristic was observable only in hub proteins or whether it existed broadly at the proteome level; and, 2 ) whether this was a spurious correlation that had emerged because of the presence of some third factor hidden in the complex and intertwining proteomes ., Our results show that , in many cases , this characteristic is present not only in hub proteins but also in proteins that have few interacting partners ., They also imply that these protein interactions or binding patterns are not the result of influence by a third factor but are caused by phosphorylation-dependent cellular activities ., The additive effect of kinase–substrate and phosphatase–substrate reactions is one possible model for interpreting this phenomenon in the phospho-PPI network ., If PPIs include many transient signaling reactions between kinases , phosphatases , and their substrates ( most of which are phosphorylated under certain conditions ) , then the signaling proteins may have interactions additional to the cohesive protein binding interactions in the PPI data ., Indeed , some enzyme–protein substrate interactions are surprisingly stable and can be captured in protein interaction assays 46 ., However , of the 795 yeast phosphorylation and dephosphorylation reactions for which information has previously been published 47 , only 3 . 9% , 1 . 6% , 2 . 4% , and 0 . 8% overlapped with those in our “ALL , ” “Y2H , ” “IMM , ” and “TAP” PPI datasets , respectively ( Supplementary Fig . S6 ) ., Note , however , that these values were significantly higher than those expected from negative controls of the corresponding PPI networks generated by “random edge rewiring” ( RER ) , and similar , significant overlaps between physical PPI and signaling network were obtained by another group 48; for details of RER , see Materials and Methods .,  On the other hand , the node degree levels of at least 600 proteins ( >20% of phosphoproteomes in the “ALL” phospho-PPI network ) might have been related to , and affected by , phosphorylation , as evidenced by the cumulative percentage of phosphoproteins , which was more than 20% higher than that of nonphosphoproteins ( Fig . 3A ) ., In addition to this , many unidentified phosphoproteins are certain to be present in the nonphosphoprotein dataset ., Therefore , it is difficult to interpret such a large difference in the node degree of phosphoproteins and nonphosphoproteins only in terms of the additive effect of signaling reactions , which had such a small overlap with the PPI data ., Furthermore , among the GO Slim ontology groups within the “signal transduction” and “cell cycle” categories , which especially include many signaling proteins , there were no great distinctions between the node degree levels of phosphoproteins and nonphosphoproteins ( although the node degree levels for “cytokinesis” and “response to stress , ” like those for most of the other ontology groups , showed marked differences between phosphoproteins and nonphosphoproteins ) ( Supplementary Fig . S5 ) ., In the phospho-PPI network , phosphoproteins had a greater tendency than nonphosphoproteins to interact with proteins harboring phosphoprotein binding domains ( PPBDs ) ., Out of 10 known PPBDs—14-3-3 , BRCT , C2 , FHA , MH2 , PBD , PTB , SH2 , WD-40 , and WW 49—six ( BRCT , C2 , FHA , SH2 , WD-40 , and WW ) were present in the member proteins of the “ALL” phospho-PPI network , and the average probabilities that phosphoproteins would interact with proteins that had all PPBDs or each type of PPBD were higher than those for nonphosphoproteins ( Fig . 5 ) ., ( The gap between node degree levels of phosphoproteins and nonphosphoproteins was normalized; see Materials and Methods . ), Considering all of these results and perspectives , a reasonable and generalized model that can be used to interpret the higher node degree of phosphoproteins is that reversible and alternative phosphorylation reactions alter the physical characteristics of proteins under various environmental conditions; the interacting or binding partners of phosphoproteins are thereby more diversified than those of nonphosphorylated proteins ., Consistent with this interpretation , phosphoproteins harboring at least two phosphosites had more interacting partners than those with a single phosphosite in the phospho-PPI network ( Supplementary Fig . S7 ) , even though phosphoproteins follow a power-law distribution with regard to phosphosite number counts and only a small fraction of phosphoproteins have multiple phosphosites 50 ., Protein phosphorylation reactions therefore seem to make a large contribution to intracellular PPI diversity ., We further analyzed the phosphorylation patterns of protein pairs forming pair-wise interactions in the phospho-PPI network , and we found that both interacting proteins in each pair tended to be phosphorylated ., For every category of phospho-PPI network , three types of pair-wise interactions were counted , whereby “Both , ” “Either , ” or “Neither” of two interacting proteins were phosphorylated ., The “Both” and “Neither” types of protein interactions were significantly more common in the real phospho-PPI network than was expected from negative controls produced by RER , whereas the “Either” types of protein interactions were significantly less common than expected ( Fig . 6; Supplementary Fig . S8 ) ., Notably , this outcome was independent of whether the node degrees of the phosphoproteins were higher or lower than those of the nonphosphoproteins , because RER does not change the node degree of each protein in a given network 51 ., PPI data contain homodimer and heterodimer information that can be captured by experimental assays such as two-hybrid assays 52 ., Therefore , to check the possibility that the tendency of interacting proteins to have similar phosphorylation patterns was caused by protein interactions between structurally and sequentially homologous proteins with similar phosphosites , we conducted the same analysis as above but using “filtered” phospho-PPI networks , in which interactions between two homologous proteins were eliminated by E-value cut-offs of 1e–10 in the BLASTP program , but no marked change was observed ( Fig . 6; Supplementary Fig . S8 ) ., Proteins involved in signal transduction pathways tend to be phosphorylated , and this is reflected in the PPI data , although the overlaps between such signaling reactions and PPIs are limited ( see above and Supplementary Fig . S6 ) ., Another possible interpretation for the multitude of physical interactions between phosphoproteins is that physically binding proteins that are members of the same protein complex tend to be phosphorylated simultaneously by a single enzyme ., To search for the protein kinases potentially responsible for the co-phosphorylation of proteins forming the same complex , we analyzed a dataset of kinase–substrate relationships with PPI data of the “ALL” category ., In the following analysis , we used 85 and 65 kinases , respectively , from the experimental results of an in vitro kinase–substrate assay 53 and a literature-derived collection of yeast signaling reactions 47 , each having multiple substrates ( Supplementary Table S3 ) ., For each kinase , its multiple substrates were superimposed on the PPI network and the number of “interacting kinate modules” ( IKMs , triangle motifs composed of a kinase and its two physically interacting substrates ) ( Fig . 7A ) 53 was counted and compared with those estimated in negative controls of the PPI network produced by NLS and RER ., This analysis revealed that three kinases from the in vitro assay and 12 from the literature-based collection had significantly higher IKM formability than those expected from both NLS and RER ( P<0 . 05 ) ( Fig . 7B and C; Supplementary Table S3 ) ., Similar results were obtained by using the “filtered” phospho-PPI network ( Supplementary Fig . S9; Supplementary Table S3 ) ., Accordingly , we suggest that , when a protein complex and kinase are in close proximity within the intracellular environment , there is a high chance of simultaneous phosphorylation of member proteins participating in the complex ., This is consistent with the subcellular co-localization of signaling networks recently revealed through the systematic prediction of signaling networks by using phosphoproteome data with an integrated protein network information derived from curated pathway databases , co-occurring terms in abstracts , physical protein interaction assays , mRNA expression profiles , and the genomic context 48 , and by data analysis of time-course phosphoproteome data 54 ., IKMs may enhance the subcellular co-localization of signaling reactions , and/or vice versa ., The literature-derived signaling collection is presumably more enriched with well-investigated reactions and thus may more accurately reflect in vivo signaling ., This may explain why the collection harbored more kinases with high IKM formabilities ( 12 out of 65 ) than the in vitro kinase–substrate relationship data ( three out of 85 ) ., It is plausible that , in living cells , the diversity of protein interactomes ( not only of protein signaling but also of protein complex formation ) is essentially influenced by the large number of phosphorylation events; many reversible phosphorylations might control condition-specific protein binding interactions related to different subcellular processes and molecular machines ., On the other hand , protein phosphorylation patterns also seem to depend largely on intracellular protein interaction diversity ., It is possible that many of the proteins defined as nonphosphoproteins in this study can actually be phosphorylated under appropriate cellular conditions ., Even where this is true , however , the set we defined here as phosphoproteins should be enriched with proteins that are frequently phosphorylated under normal or many different cellular conditions , because the frequently phosphorylated proteins have a higher chance of being identified as phosphoproteins than do the rarely phosphorylated proteins ., Accordingly , the features and models discussed in this study should reflect the overall characteristics of phosphoproteins and nonphosphoproteins among a number of different cellular conditions ., This is supported by the finding that proteins that had two or more phosphosites physically interacted with more proteins than did those with only a single phosphosite ( Supplementary Fig . S7 ) ., Although the quality of current yeast PPI data is also not perfect and the data may include false positives , the observed features with statistical significance should be consequences of the actual behaviors of intracellular proteins , because the effects of such false positives on the statistical tests are supposedly random ., The integrative data-mining of yeast multi-omics data has now shed light on the macroscopic and large-scale relationships between phosphoproteomes and protein interactomes ., Future comprehensive analyses of the in vivo link between protein phosphorylation and physical interaction will yield more insights into the complex and intertwined molecular systems of living cells ., Saccharomyces cerevisiae strain IFO 0233 cells grown continuously on glucose medium 55 were used ., Pelleted cells were vacuum dried and frozen until further analysis ., A Bioruptor UCW-310 ( Cosmo Bio , Tokyo Japan ) was used to disrupt the pellets in 0 . 1 M Tris-HCl ( pH 8 . 0 ) containing 8 M urea , protein phosphatase inhibitor cocktails 1 and 2 ( Sigma ) , and protease inhibitors ( Sigma ) ., The homogenate was centrifuged at 1 , 500g for 10 min and the supernatant was reduced with dithiothreitol , alkylated with iodoacetamide , and digested with Lys-C; this was followed by dilution and trypsin digestion as described 56 ., Digested samples were desalted by using C-18 StageTips 57 ., Phosphopeptide enrichment by hydroxy acid–modified metal oxide chromatography ( HAMMOC ) was performed as reported previously 11 , 58 ., Briefly , digested lysates ( 100 µg each ) were loaded onto a self-packed titania-C8 StageTip in the presence of lactic acid ., After the samples had been washed with 80% acetonitrile containing 0 . 1% TFA , phosphopeptides were eluted by a modified approach using 5% ammonium hydroxide , 5% piperidine , and 5% pyrrolidine in series 59 ., An LTQ-Orbitrap XL ( Thermo Fisher Scientific , Bremen , Germany ) coupled with a Dionex Ultimate 3000 ( Germering , Germany ) and an HTC-PAL autosampler ( CTC Analytics AG , Zwingen , Switzerland ) was used for nanoLC-MS/MS analyses ., An analytical column needle with a “stone-arch” frit 60 was prepared with ReproSil C18 materials ( 3 µm , Dr . Maisch , Ammerbuch , Germany ) ., The injection volume was 5 µL and the flow rate was 500 nL/min ., The mobile phases consisted of ( A ) 0 . 5% acetic acid and ( B ) 0 . 5% acetic acid and 80% acetonitrile ., A three-step linear gradient of 5% to 10% B in 5 min , 10% to 40% B in 60 min , 40% to 100% B in 5 min , and 100% B for 10 min was employed throughout this study ., The MS scan range was m/z 300 to 1500 , and the top 10 precursor ions were selected in MS scans by Orbitrap with R\u200a=\u200a60 , 000 for subsequent MS/MS scans by ion trap in the automated gain control ( AGC ) mode; AGC values of 5 . 00e+05 and 1 . 00e+04 were set for full MS and MS/MS , respectively ., The normalized collision energy was set at 35 . 0 ., A lock mass function was used for the LTQ-Orbitrap to obtain constant mass accuracy during gradient analysis ., Both Mass Navigator v1 . 2 ( Mitsui Knowledge Industry , Tokyo , Japan ) and Mascot Distiller v2 . 2 . 1 . 0 ( Matrix Science , London , UK ) were used to create peak lists based on the recorded fragmentation spectra ., Peptides and proteins were identified by automated database searching using Mascot Server v2 . 2 ( Matrix Science ) against UniProt/SwissProt v56 . 0 with a precursor mass tolerance of 3 ppm , a fragment ion mass tolerance of 0 . 8 Da , and strict trypsin specificity , allowing for up to two missed cleavages ., Carbamidomethylation of cysteine was set as a fixed modification , and oxidation of methionines and phosphorylation of serine , threonine , and tyrosine were allowed as variable modifications ., Phosphopeptide identification and phosphorylated site determination were performed in accordance with a procedure reported previously 11 ., The false discovery rate was estimated to be 1 . 07% using a randomized database ., All annotated MS/MS spectra were stored in PepBase ( http://pepbase . iab . keio . ac . jp ) ., Saccharomyces cerevisiae phosphoproteome data were obtained from Dataset S1 of Holt et al . 16 ., Another collection of formerly identified phosphoproteins and their phosphosites was obtained from UniProt ( release 15 . 14; http://www . uniprot . org/ ) 19 ., All UniProtKB/Swiss-Prot protein entries identified to have at least one phosphosite in high-throughput phosphoproteomics studies were downloaded via the Protein Knowledgebase ( UniProtKB ) in XML format by querying the term scope: “PHOSPHORYLATION LARGE SCALE ANALYSIS AT” ., Some phosphoproteins registered in UniProt had multiple synonyms of UniProt accession ., For integrative analyses and comparisons of yeast multi-omics data , all identities of proteins and genes obtained from different data sources were standardized to UniProt accessions ., If objects ( e . g . gene names , ORF names , and/or locus names ) in a data source did not have UniProt accessions , the objects were standardized to their corresponding UniProt accessions according to the cross-reference list prepared from UniProtKB/Swiss-Prot protein entries obtained from UniProt ( release 15 . 14 ) ., In cases when an object corresponded to multiple synonyms of UniProt accessions , all accessions were used to identify its corresponding objects in other data sources ., The phosphoproteome data newly identified in this study and the former phosphoproteome datasets obtained from Holt et al . and UniProt were unified according to their UniProt accessions ., Positions of phosphosites and their amino acid residues in the unified phosphoproteome data were double-checked by using the proteome sequences obtained from UniProt ( release 15 . 14 ) ., From SGD ( Saccharomyces Genome Database; http://yeastgenome . org ) 32 , annotations of 5 , 815 known and predicted genes were obtained ., ORF names of genes were checked by using the unified phosphoproteome data to determine whether the encoded protein was identified as a phosphoprotein ., The S . cerevisiae PPI network was obtained as XML files ( Scere20081014 ) from DIP ( Database of Interacting Proteins; http://dip . doe-mbi . ucla . edu ) 33 ., We eliminated each interaction entry including three or more “interactors” ( e . g . , in which multiple prey proteins were detected for one bait protein in one experimental assay ) and used only those including two “interactors . ”, Every node in the PPI network was labeled by its corresponding UniProt ID provided in the same XML file ., For the PPI assay , PPI data were further grouped into four categories: all kinds of experimental methods ( “ALL” ) , yeast two-hybrid ( “Y2H” ) , co-immunoprecipitation ( “IMM” ) , and tandem affinity purification ( “TAP” ) ., A “filtered” PPI network was also prepared for each category by eliminating interactions between two similar proteins by using the BLASTP program and an E-value cut-off of 1e–10 ., Unified phosphoproteome data were mapped onto every category of PPI data prepared from DIP according to their UniProt accessions , and a phospho-PPI network was generated ., Throughout this study , proteins that did not correspond to phosphoproteome data were termed “nonphosphoproteins . ”, To prepare negative controls for PPI and phospho-PPI networks , two different processes ( as diagrammed in Fig . 1 ) were appropriately adopted on a case-by-case basis ., “Node label shuffling” ( NLS ) swaps the labels of two randomly selected nodes in a given network; it repeats this operation a sufficient number of times until all pair-wise interactions in the queried network have disappeared or until the number of iterations reaches 1 , 000 times the number of interactions ., “Random edge rewiring” ( RER ) randomly selects two edges in a given network and randomly rewires them ., During this process , each rewiring operation is retried if a pair of nodes redundantly wired by two edges occurs in the network; the iteration termination condition is the same as that of NLS ., Proteome abundance data for S . cerevisiae that were previously acquired through a single-cell proteomics analysis combining high-throughput flow cytometry and a library of GFP-tagged strains 37 were used to analyze the characteristics of protein expression in the phospho-PPI network ., These data were composed of proteome abundance data measured for cells grown in rich ( YEPD ) and synthetic complete ( SD ) medium ., For each cell growth condition , protein names were standardized to UniProt accessions , and protein abundance levels were log-transformed ( base 10 ) and superimposed on each of the phospho-PPI networks of “ALL” and “Y2H . ”, In this case , protein nodes for which the abundance levels were not provided in the abundance data were deleted from the phospho-PPI network ., The protein disorder level of every S . cerevisiae protein registered in UniProt ( release 15 . 14 ) was predicted by the POODLE-W program , which uses the support vector machine–based learning of amino acid sequences of structurally confirmed disordered proteins 61 ., For the analysis , we used the “disorder probability” ( i . e . the probability that a given protein is unstructured ) output by this program ., Saccharomyces cerevisiae gene annotations belonging to “molecular function , ” “biological process , ” or “cellular component” of GO Slim , a higher level view of S . cerevisiae Gene Ontology ( GO ) , were downloaded via the SGD ftp site ., Information on S . cerevisiae proteins , each of which has at least one of 10 known phosphoprotein binding domains ( PPBDs ) , namely 14-3-3 , BRCT , C2 , FHA , MH2 , PBD , PTB , SH2 , WD-40 , and WW 49 , was obtained according to the protein domain annotations of UniProt ( release 15 . 14 ) , which were provided by other protein databases ., To evaluate the tendencies of phosphoproteins and nonphosphoproteins to interact with proteins that had PPBDs , the normalized probabilities of such interactions were defined ., For each protein , the number of interacting protein partners that had PPBDs was divided by the number of all interacting partners ., To find possible IKMs , kinases previously reported to phosphorylate multiple substrates were obtained from data on in vitro substrates recognized by most yeast protein kinases that were measured with the use of proteome chip technology Supplementary Data 2 of Ptacek et al . 53 , as well as from a literature-derived collection of documented yeast signaling reactions Table S3 of Fiedler et al . 47 ., All gene names of substrates in the in vitro kinase–substrate relationship data and ORF name
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 .
To date , high-throughput proteome technologies have revealed that hundreds to thousands of proteins in each of many organisms are phosphorylated under the appropriate environmental conditions ., A critical role of phosphorylation is control of protein signaling ., However , only a fraction of the identified phosphoproteins participate in currently known protein signaling pathways , and the biological relevance of the remainder is unclear ., This has raised the question of whether phosphorylation has other major roles ., In this study , we identified new phosphoproteins in budding yeast by mass spectrometry and unified these new data with publicly available phosphoprotein data ., We then performed an integrative data-mining of large-scale yeast phosphoproteins and protein–protein interactions ( complex formation ) by an exhaustive analysis that incorporated yeast protein information from several other sources ., The phosphoproteome data integration surprisingly showed that nearly 60% of yeast genes encode phosphoproteins , and the subsequent data-mining analysis derived two models interpreting the mutual intracellular effects of large-scale protein phosphorylation and binding interaction ., Biological interpretations of both large-scale intracellular phosphorylation and the topology of protein interaction networks are highly relevant to modern biology ., This study sheds light on how in vivo protein pathways are supported by a combination of protein modification and molecular dynamics .
computational biology/systems biology