Phosphoproteomic analysis reveals interconnected system-wide responses to perturbations of kinases and phosphatases in yeast

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.

Regulation of pol III transcription by nutrient and stress signaling pathways

Transcription by RNA polymerase III (pol III) is responsible for ~15% of total cellular transcription through the generation of small structured RNAs such as tRNA and 5S RNA. The coordinate synthesis of these molecules with ribosomal protein mRNAs and rRNA couples the production of ribosomes and their tRNA substrates and balances protein synthetic capacity with the growth requirements of the cell. Ribosome biogenesis in general and pol III transcription in particular is known to be regulated by nutrient availability, cell stress and cell cycle stage and is perturbed in pathological states. High throughput proteomic studies have catalogued modifications to pol III subunits, assembly, initiation and accessory factors but most of these modifications have yet to be linked to functional consequences. Here we review our current understanding of the major points of regulation in the pol III transcription apparatus, the targets of regulation and the signaling pathways known to regulate their function. This article is part of a Special Issue entitled: Transcription by Odd Pols.

Profiling the kinome: current capabilities and future challenges

Protein kinases are the second largest human protein family, but in terms of research interest, both basic and applied, they are surely the most popular. Over the past decade, many techniques and approaches for studying the kinome have been described and the pace of development is ever increasing. Presently, a molecular biologist can approach the kinome from many different angles: what kinases are active during a specific cell state of interest or become activated in response to a specific stimulus? What are the effects of controlling the activation status of an individual kinase? What substrates are targeted by a particular kinase, either in general or under particular conditions? And what kinase is responsible for targeting a specific phosphorylation site of interest? These are some of the more commonly asked questions during any kinase-centric research project and different strategies have been devised for answering such queries. In this review, we outline the most promising of these approaches, particularly those with a capacity for high-throughput studies. This article is part of a Special Issue entitled: From protein structures to clinical applications.

Response to hyperosmotic stress

An appropriate response and adaptation to hyperosmolarity, i.e., an external osmolarity that is higher than the physiological range, can be a matter of life or death for all cells. It is especially important for free-living organisms such as the yeast Saccharomyces cerevisiae. When exposed to hyperosmotic stress, the yeast initiates a complex adaptive program that includes temporary arrest of cell-cycle progression, adjustment of transcription and translation patterns, and the synthesis and retention of the compatible osmolyte glycerol. These adaptive responses are mostly governed by the high osmolarity glycerol (HOG) pathway, which is composed of membrane-associated osmosensors, an intracellular signaling pathway whose core is the Hog1 MAP kinase (MAPK) cascade, and cytoplasmic and nuclear effector functions. The entire pathway is conserved in diverse fungal species, while the Hog1 MAPK cascade is conserved even in higher eukaryotes including humans. This conservation is illustrated by the fact that the mammalian stress-responsive p38 MAPK can rescue the osmosensitivity of hog1Δ mutations in response to hyperosmotic challenge. As the HOG pathway is one of the best-understood eukaryotic signal transduction pathways, it is useful not only as a model for analysis of osmostress responses, but also as a model for mathematical analysis of signal transduction pathways. In this review, we have summarized the current understanding of both the upstream signaling mechanism and the downstream adaptive responses to hyperosmotic stress in yeast.

Validation of regulated protein phosphorylation events in yeast by quantitative mass spectrometry analysis of purified proteins

Global phosphoproteomic studies based on MS have generated qualitative and quantitative data describing protein phosphorylation events in various biological systems. Since high-throughput data for protein modifications are inherently incomplete, we developed a strategy to extend and validate such primary datasets. We selected interesting protein candidates from a global screen in yeast and employed a modified histidine biotin tag that allows tandem affinity purifications of our targets under denaturing conditions. Products in question can be digested directly from affinity resins and phosphopeptides can be further enriched via TiO(2) before MS analysis. Our robust protocol can be amended for SILAC as well as iTRAQ quantifications or label-free approaches based on selective reaction monitoring, allowing completion of the phosphorylation pattern in a first step, followed by a detailed analysis of the phosphorylation kinetics. We exemplify the value of such a strategy by an in-depth analysis of Pan1, a highly phosphorylated factor involved in early steps of endocytosis. The study of Pan1 under osmotic stress conditions in different mutant backgrounds allowed us to differentiate between mitogen-activated protein kinase Hog1 driven and Hog1 independent stress responses.

Mass spectrometry-based proteomics and network biology

In the life sciences, a new paradigm is emerging that places networks of interacting molecules between genotype and phenotype. These networks are dynamically modulated by a multitude of factors, and the properties emerging from the network as a whole determine observable phenotypes. This paradigm is usually referred to as systems biology, network biology, or integrative biology. Mass spectrometry (MS)-based proteomics is a central life science technology that has realized great progress toward the identification, quantification, and characterization of the proteins that constitute a proteome. Here, we review how MS-based proteomics has been applied to network biology to identify the nodes and edges of biological networks, to detect and quantify perturbation-induced network changes, and to correlate dynamic network rewiring with the cellular phenotype. We discuss future directions for MS-based proteomics within the network biology paradigm.

Regulation of yeast pyruvate kinase by ultrasensitive allostery independent of phosphorylation

Allostery and covalent modification are major means of fast-acting metabolic regulation. Their relative roles in responding to environmental changes remain, however, unclear. Here we examine this issue, using as a case study the rapid decrease in pyruvate kinase flux in yeast upon glucose removal. The main pyruvate kinase isozyme (Cdc19) is phosphorylated in response to environmental cues. It also exhibits positively cooperative (ultrasensitive) allosteric activation by fructose-1,6-bisphosphate (FBP). Glucose removal causes accumulation of Cdc19's substrate, phosphoenolpyruvate. This response is retained in strains with altered protein-kinase-A or AMP-activated-protein-kinase activity or with CDC19 carrying mutated phosphorylation sites. In contrast, yeast engineered with a CDC19 point mutation that ablates FBP-based regulation fail to accumulate phosphoenolpyruvate. They also fail to grow on ethanol and slowly resume growth upon glucose upshift. Thus, while yeast pyruvate kinase is covalently modified in response to glucose availability, its activity is controlled almost exclusively by ultrasensitive allostery.

Mathematical modelling of the MAP kinase pathway using proteomic datasets

The advances in proteomics technologies offer an unprecedented opportunity and valuable resources to understand how living organisms execute necessary functions at systems levels. However, little work has been done up to date to utilize the highly accurate spatio-temporal dynamic proteome data generated by phosphoprotemics for mathematical modeling of complex cell signaling pathways. This work proposed a novel computational framework to develop mathematical models based on proteomic datasets. Using the MAP kinase pathway as the test system, we developed a mathematical model including the cytosolic and nuclear subsystems; and applied the genetic algorithm to infer unknown model parameters. Robustness property of the mathematical model was used as a criterion to select the appropriate rate constants from the estimated candidates. Quantitative information regarding the absolute protein concentrations was used to refine the mathematical model. We have demonstrated that the incorporation of more experimental data could significantly enhance both the simulation accuracy and robustness property of the proposed model. In addition, we used the MAP kinase pathway inhibited by phosphatases with different concentrations to predict the signal output influenced by different cellular conditions. Our predictions are in good agreement with the experimental observations when the MAP kinase pathway was inhibited by phosphatase PP2A and MKP3. The successful application of the proposed modeling framework to the MAP kinase pathway suggests that our method is very promising for developing accurate mathematical models and yielding insights into the regulatory mechanisms of complex cell signaling pathways.

Proteome-wide analysis of lysine acetylation suggests its broad regulatory scope in Saccharomyces cerevisiae

Post-translational modification of proteins by lysine acetylation plays important regulatory roles in living cells. The budding yeast Saccharomyces cerevisiae is a widely used unicellular eukaryotic model organism in biomedical research. S. cerevisiae contains several evolutionary conserved lysine acetyltransferases and deacetylases. However, only a few dozen acetylation sites in S. cerevisiae are known, presenting a major obstacle for further understanding the regulatory roles of acetylation in this organism. Here we use high resolution mass spectrometry to identify about 4000 lysine acetylation sites in S. cerevisiae. Acetylated proteins are implicated in the regulation of diverse cytoplasmic and nuclear processes including chromatin organization, mitochondrial metabolism, and protein synthesis. Bioinformatic analysis of yeast acetylation sites shows that acetylated lysines are significantly more conserved compared with nonacetylated lysines. A large fraction of the conserved acetylation sites are present on proteins involved in cellular metabolism, protein synthesis, and protein folding. Furthermore, quantification of the Rpd3-regulated acetylation sites identified several previously known, as well as new putative substrates of this deacetylase. Rpd3 deficiency increased acetylation of the SAGA (Spt-Ada-Gcn5-Acetyltransferase) complex subunit Sgf73 on K33. This acetylation site is located within a critical regulatory domain in Sgf73 that interacts with Ubp8 and is involved in the activation of the Ubp8-containing histone H2B deubiquitylase complex. Our data provides the first global survey of acetylation in budding yeast, and suggests a wide-ranging regulatory scope of this modification. The provided dataset may serve as an important resource for the functional analysis of lysine acetylation in eukaryotes.

Analysing signalling networks by mass spectrometry

Sequence analysis of the human genome and the association of genetic aberrations with diseases have provided a rough framework whereby the impact of individual genotypes can be assessed. To fully understand the effect of individual and co-occurring genetic aberrations, as well as their individual and collected contribution to the development of diseases, it is critical to analyse the matching proteome and to determine how the organisation, expression level and function of protein networks are affected. Sensitive mass spectrometric platforms in combination with innovative workflows allow qualitative and quantitative analyses of the cellular as well as the extracellular proteome. Importantly, in addition to specifically identifying the content of the proteome, several aspects of the proteomic organisation can be analysed including protein complexes, protein modifications, enzymatic activities and subcellular/organelle localisation. Together, these measurements will provide novel insight into the biological effect of disease-causing mutations ultimately coupling genotype and phenotype.

An N-terminal acidic region of Sgs1 interacts with Rpa70 and recruits Rad53 kinase to stalled forks

DNA replication fork stalling poses a major threat to genome stability. This is counteracted in part by the intra-S phase checkpoint, which stabilizes arrested replication machinery, prevents cell-cycle progression and promotes DNA repair. The checkpoint kinase Mec1/ATR and RecQ helicase Sgs1/BLM contribute synergistically to fork maintenance on hydroxyurea (HU). Both enzymes interact with replication protein A (RPA). We identified and deleted the major interaction sites on Sgs1 for Rpa70, generating a mutant called sgs1-r1. In contrast to a helicase-dead mutant of Sgs1, sgs1-r1 did not significantly reduce recovery of DNA polymerase α at HU-arrested replication forks. However, the Sgs1 R1 domain is a target of Mec1 kinase, deletion of which compromises Rad53 activation on HU. Full activation of Rad53 is achieved through phosphorylation of the Sgs1 R1 domain by Mec1, which promotes Sgs1 binding to the FHA1 domain of Rad53 with high affinity. We propose that the recruitment of Rad53 by phosphorylated Sgs1 promotes the replication checkpoint response on HU. Loss of the R1 domain increases lethality selectively in cells lacking Mus81, Slx4, Slx5 or Slx8.

Identification and validation of inhibitor-responsive kinase substrates using a new paradigm to measure kinase-specific protein phosphorylation index

Regulation of all cellular processes requires dynamic regulation of protein phosphorylation. We have developed an unbiased system to globally quantify the phosphorylation index for substrates of a specific kinase by independently quantifying phosphorylated and total substrate molecules in a reverse in-gel kinase assay. Non-phosphorylated substrate molecules are first quantified in the presence and absence of a specific stimulus. Total substrate molecules are then measured after complete chemical dephosphorylation, and a ratio of phosphorylated to total substrate is derived. To demonstrate the utility of this approach, we profiled and quantified changes in phosphorylation index for Protein Kinase CK2 substrates that respond to a small-molecule inhibitor. A broad range of inhibitor-induced changes in phosphorylation was observed in cultured cells. Differences among substrates in the kinetics of phosphorylation change were also revealed. Comparison of CK2 inhibitor-induced changes in phosphorylation in cultured cells and in mouse peripheral blood lymphocytes in vivo revealed distinct kinetic and depth-of-response profiles. This technology provides a new approach to facilitate functional analyses of kinase-specific phosphorylation events. This strategy can be used to dissect the role of phosphorylation in cellular events, to facilitate kinase inhibitor target validation studies, and to inform in vivo analyses of kinase inhibitor drug efficacy.

M-Track: detecting short-lived protein-protein interactions in vivo

We developed a protein-proximity assay in yeast based on fusing a histone lysine methyltransferase onto a bait and its substrate onto a prey. Upon binding, the prey is stably methylated and detected by methylation-specific antibodies. We applied this approach to detect varying interaction affinities among proteins in a mitogen-activated protein kinase pathway and to detect short-lived interactions between protein phosphatase 2A and its substrates that have so far escaped direct detection.

Research resource: identification of novel growth hormone-regulated phosphorylation sites by quantitative phosphoproteomics

GH and GH receptors are expressed throughout life, and GH elicits a diverse range of responses, including growth and altered metabolism. It is therefore important to understand the full spectrum of GH signaling pathways and cellular responses. We applied mass spectrometry-based phosphoproteomics combined with stable isotope labeling with amino acids in cell culture to identify proteins rapidly phosphorylated in response to GH in 3T3-F442A preadipocytes. We identified 132 phosphosites in 95 proteins that exhibited rapid (5 or 15 min) GH-dependent statistically significant increases in phosphorylation by more than or equal to 50% and 96 phosphosites in 46 proteins that were down-regulated by GH by more than or equal to 30%. Several of the GH-stimulated phosphorylation sites were known (e.g. regulatory Thr/Tyr in Erks 1 and 2, Tyr in signal transducers and activators of transcription (Stat) 5a and 5b, Ser939 in tuberous sclerosis protein (TSC) 2 or tuberin). The remaining 126 GH-stimulated sites were not previously associated with GH. Kyoto Encyclopedia of Genes and Genomes pathway analysis of GH-stimulated sites indicated enrichment in proteins associated with the insulin and mammalian target of rapamycin (mTOR) pathways, regulation of the actin cytoskeleton, and focal adhesions. Akt/protein kinase A consensus sites (RXRXXS/T) were the most commonly phosphorylated consensus sites. Immunoblotting confirmed GH-stimulated phosphorylation of all seven novel GH-dependent sites tested [regulatory sites in proline-rich Akt substrate, 40 kDA (PRAS40), regulatory associated protein of mTOR, ATP-citrate lyase, Na⁺/H⁺ exchanger-1, N-myc downstream regulated gene 1, and Shc]). The immunoblot results suggest that many, if not most, of the GH-stimulated phosphosites identified in this large-scale quantitative phosphoproteomics analysis, including sites in multiple proteins in the Akt/ mTOR complex 1 pathway, are phosphorylated in response to GH. Their identification significantly broadens our thinking of GH-regulated cell functions.

Phosphoproteomic analysis of leukemia cells under basal and drug-treated conditions identifies markers of kinase pathway activation and mechanisms of resistance

Protein kinase signaling is fundamental to cell homeostasis and is deregulated in all cancers but varies between patients. Understanding the mechanisms underlying this heterogeneity is critical for personalized targeted therapies. Here, we used a recently established LC-MS/MS platform to profile protein phosphorylation in acute myeloid leukemia cell lines with different sensitivities to kinase inhibitors. The compounds used in this study were originally developed to target Janus kinase, phosphatidylinositol 3-kinase, and MEK. After further validation of the technique, we identified several phosphorylation sites that were inhibited by these compounds but whose intensities did not always correlate with growth inhibition sensitivity. In contrast, several hundred phosphorylation sites that correlated with sensitivity/resistance were not in general inhibited by the compounds. These results indicate that markers of pathway activity may not always be reliable indicators of sensitivity of cancer cells to inhibitors that target such pathways, because the activity of parallel kinases can contribute to resistance. By mining our data we identified protein kinase C isoforms as one of such parallel pathways being more active in resistant cells. Consistent with the view that several parallel kinase pathways were contributing to resistance, inhibitors that target protein kinase C, MEK, and Janus kinase potentiated each other in arresting the proliferation of multidrug-resistant cells. Untargeted/unbiased approaches, such as the one described here, to quantify the activity of the intended target kinase pathway in concert with the activities of parallel kinase pathways will be invaluable to personalize therapies based on kinase inhibitors.

Computational approaches for analyzing information flow in biological networks

The advancements in "omics" (proteomics, genomics, lipidomics, and metabolomics) technologies have yielded large inventories of genes, transcripts, proteins, and metabolites. The challenge is to find out how these entities work together to regulate the processes by which cells respond to external and internal signals. Mathematical and computational modeling of signaling networks has a key role in this task, and network analysis provides insights into biological systems and has applications for medicine. Here, we review experimental and theoretical progress and future challenges toward this goal. We focus on how networks are reconstructed from data, how these networks are structured to control the flow of biological information, and how the design features of the networks specify biological decisions.

Global detection of protein kinase D-dependent phosphorylation events in nocodazole-treated human cells

Protein kinase D (PKD) is a cytosolic serine/threonine kinase implicated in regulation of several cellular processes such as response to oxidative stress, directed cell migration, invasion, differentiation, and fission of the vesicles at the trans-Golgi network. Its variety of functions must be mediated by numerous substrates; however, only a couple of PKD substrates have been identified so far. Here we perform stable isotope labeling of amino acids in cell culture-based quantitative phosphoproteomic analysis to detect phosphorylation events dependent on PKD1 activity in human cells. We compare relative phosphorylation levels between constitutively active and kinase dead PKD1 strains of HEK293 cells, both treated with nocodazole, a microtubule-depolymerizing reagent that disrupts the Golgi complex and activates PKD1. We identify 124 phosphorylation sites that are significantly down-regulated upon decrease of PKD1 activity and show that the PKD target motif is significantly enriched among down-regulated phosphorylation events, pointing to the presence of direct PKD1 substrates. We further perform PKD1 target motif analysis, showing that a proline residue at position +1 relative to the phosphorylation site serves as an inhibitory cue for PKD1 activity. Among PKD1-dependent phosphorylation events, we detect predominantly proteins with localization at Golgi membranes and function in protein sorting, among them several sorting nexins and members of the insulin-like growth factor 2 receptor pathway. This study presents the first global detection of PKD1-dependent phosphorylation events and provides a wealth of information for functional follow-up of PKD1 activity upon disruption of the Golgi network in human cells.

Dissecting phosphorylation networks: lessons learned from yeast

Protein phosphorylation continues to be regarded as one of the most important post-translational modifications found in eukaryotes and has been implicated in key roles in the development of a number of human diseases. In order to elucidate roles for the 518 human kinases, phosphorylation has routinely been studied using the budding yeast Saccharomyces cerevisiae as a model system. In recent years, a number of technologies have emerged to globally map phosphorylation in yeast. In this article, we review these technologies and discuss how these phosphorylation mapping efforts have shed light on our understanding of kinase signaling pathways and eukaryotic proteomic networks in general.

Phosphoproteomic analysis: an emerging role in deciphering cellular signaling in human embryonic stem cells and their differentiated derivatives

Cellular signaling is largely controlled by protein phosphorylation. This post-translational modification (PTM) has been extensively analyzed when examining one or a few protein phosphorylation events that effect cell signaling. However, protein kinase-driven signaling networks, comprising total (phospho)proteomes, largely control cell fate. Therefore, large-scale analysis of differentially regulated protein phosphorylation is central to elucidating complex cellular events, including maintenance of pluripotency and differentiation of embryonic stem cells (ESCs). The current technology of choice for total phosphoproteome and combined total proteome plus total phosphoproteome (termed (phospho)proteome) analyses is multidimensional liquid chromatography-(MDLC) tandem mass spectrometry (MS/MS). Advances in the use of MDLC for separation of peptides comprising total (phospho)proteomes, phosphopeptide enrichment, separation of enriched fractions, and quantitative peptide identification by MS/MS have been rapid in recent years, as have improvements in the sensitivity, speed, and accuracy of mass spectrometers. Increasingly deep coverage of (phospho)proteomes is allowing an improved understanding of changes in protein phosphorylation networks as cells respond to stimuli and progress from one undifferentiated or differentiated state to another. Although MDLC-MS/MS studies are powerful, understanding the interpretation of the data is important, and targeted experimental pursuit of biological predictions provided by total (phospho)proteome analyses is needed. (Phospho)proteomic analyses of pluripotent stem cells are in their infancy at this time. However, such studies have already begun to contribute to an improved and accelerated understanding of basic pluripotent stem cell signaling and fate control, especially at the systems-biology level.

Cross-talk between phosphorylation and lysine acetylation in a genome-reduced bacterium

The effect of kinase, phosphatase and N-acetyltransferase deletions on proteome phosphorylation and acetylation was investigated in Mycoplasma pneumoniae. Bi-directional cross-talk between post-transcriptional modifications suggests an underlying regulatory molecular code in prokaryotes.

Post-translational modifications (PTMs) change the chemical properties of proteins, conferring diversity beyond the amino-acid sequence. Proteins are often modified on multiple sites. A PTM code has been proposed, whereby modifications at specific positions influence further modifications. These regulatory circuits though have rarely been studied on a large-scale; conservation in prokaryotes remains elusive.

Here, we studied two important PTMs– phosphorylation and lysine acetylation in the small bacterium Mycoplasma pneumoniae. We combined genetics and quantitative mass spectrometry to measure the effect of systematic kinase, phosphatase and N-acetyltransferase deletions on proteome abundance, phosphorylation and lysine acetylation.

The data set represents a comprehensive analysis of both phosphorylation and lysine acetylation in a single prokaryote. It reveals (1) proteins often carry multiple modifications and multiple types of PTMs, reminiscent of the PTM code proposed in eukaryotes, (2) phosphorylation exerts pleiotropic effect on proteins abundances, phosphorylation, but also lysine acetylation, (3) the cross-talk between the two PTMs is bi-directional and (4) PTMs are frequently located at interaction interfaces and in multifunctional proteins, illustrating how PTMs could modulate protein functions affecting the way they interact.

The study provides an unbiased and quantitative view on cross-talk between phosphorylation and lysine acetylation. It suggests that these regulatory circuits are a fundamental principle of regulation that might have evolved before the divergence of prokaryotes and eukaryotes.

Protein post-translational modifications (PTMs) represent important regulatory states that when combined have been hypothesized to act as molecular codes and to generate a functional diversity beyond genome and transcriptome. We systematically investigate the interplay of protein phosphorylation with other post-transcriptional regulatory mechanisms in the genome-reduced bacterium Mycoplasma pneumoniae. Systematic perturbations by deletion of its only two protein kinases and its unique protein phosphatase identified not only the protein-specific effect on the phosphorylation network, but also a modulation of proteome abundance and lysine acetylation patterns, mostly in the absence of transcriptional changes. Reciprocally, deletion of the two putative N-acetyltransferases affects protein phosphorylation, confirming cross-talk between the two PTMs. The measured M. pneumoniae phosphoproteome and lysine acetylome revealed that both PTMs are very common, that (as in Eukaryotes) they often co-occur within the same protein and that they are frequently observed at interaction interfaces and in multifunctional proteins. The results imply previously unreported hidden layers of post-transcriptional regulation intertwining phosphorylation with lysine acetylation and other mechanisms that define the functional state of a cell.

Functional wiring of the yeast kinome revealed by global analysis of genetic network motifs

A combinatorial genetic perturbation strategy was applied to interrogate the yeast kinome on a genome-wide scale. We assessed the global effects of gene overexpression or gene deletion to map an integrated genetic interaction network of synthetic dosage lethal (SDL) and loss-of-function genetic interactions (GIs) for 92 kinases, producing a meta-network of 8700 GIs enriched for pathways known to be regulated by cognate kinases. Kinases most sensitive to dosage perturbations had constitutive cell cycle or cell polarity functions under standard growth conditions. Condition-specific screens confirmed that the spectrum of kinase dosage interactions can be expanded substantially in activating conditions. An integrated network composed of systematic SDL, negative and positive loss-of-function GIs, and literature-curated kinase-substrate interactions revealed kinase-dependent regulatory motifs predictive of novel gene-specific phenotypes. Our study provides a valuable resource to unravel novel functional relationships and pathways regulated by kinases and outlines a general strategy for deciphering mutant phenotypes from large-scale GI networks.

Phosphosite mapping of P-type plasma membrane H+-ATPase in homologous and heterologous environments

Phosphorylation is an important posttranslational modification of proteins in living cells and primarily serves regulatory purposes. Several methods were employed for isolating phosphopeptides from proteolytically digested plasma membranes of Arabidopsis thaliana. After a mass spectrometric analysis of the resulting peptides we could identify 10 different phosphorylation sites in plasma membrane H(+)-ATPases AHA1, AHA2, AHA3, and AHA4/11, five of which have not been reported before, bringing the total number of phosphosites up to 11, which is substantially higher than reported so far for any other P-type ATPase. Phosphosites were almost exclusively (9 of 10) in the terminal regulatory domains of the pumps. The AHA2 isoform was subsequently expressed in the yeast Saccharomyces cerevisiae. The plant protein was phosphorylated at multiple sites in yeast, and surprisingly, seven of nine of the phosphosites identified in AHA2 were identical in the plant and fungal systems even though none of the target sequences in AHA2 show homology to proteins of the fungal host. These findings suggest an unexpected accessibility of the terminal regulatory domain of plasma membrane H(+)-ATPase to protein kinase action.

Mass spectrometry-based proteomics for systems biology

Mass spectrometry (MS)-based proteomics has significantly contributed to the development of systems biology, a new paradigm for the life sciences in which biological processes are addressed in terms of dynamic networks of interacting molecules. Because of its advanced analytical capabilities, MS-based proteomics has been used extensively to identify the components of biological systems, and it is the method of choice to consistently quantify the effects of network perturbation in time and space. Herein, we review recent contributions of MS to systems biology and discuss several examples that illustrate the importance of mass spectrometry to elucidate the components and interactions of molecular networks.

Dependence of Chs2 ER export on dephosphorylation by cytoplasmic Cdc14 ensures that septum formation follows mitosis

Sequestration of Cdc14 from the cytoplasm ensures Chs2 ER retention after MEN activation. The interdependence of chromosome segregation, MEN activation, decrease in mitotic CDK activity, and Cdc14 dispersal provides an effective mechanism for cells to order late mitotic events.

Cytokinesis, which leads to the physical separation of two dividing cells, is normally restrained until after nuclear division. In Saccharomyces cerevisiae, chitin synthase 2 (Chs2), which lays down the primary septum at the mother–daughter neck, also ensures proper actomyosin ring constriction during cytokinesis. During the metaphase-to-anaphase transition, phosphorylation of Chs2 by the mitotic cyclin-dependent kinase (Cdk1) retains Chs2 at the endoplasmic reticulum (ER), thereby preventing its translocation to the neck. Upon Cdk1 inactivation at the end of mitosis, Chs2 is exported from the ER and targeted to the neck. The mechanism for triggering Chs2 ER export thus far is unknown. We show here that Chs2 ER export requires the direct reversal of the inhibitory Cdk1 phosphorylation sites by Cdc14 phosphatase, the ultimate effector of the mitotic exit network (MEN). We further show that only Cdc14 liberated by the MEN after completion of chromosome segregation, and not Cdc14 released in early anaphase by the Cdc fourteen early anaphase release pathway, triggers Chs2 ER exit. Presumably, the reduced Cdk1 activity in late mitosis further favors dephosphorylation of Chs2 by Cdc14. Thus, by requiring declining Cdk1 activity and Cdc14 nuclear release for Chs2 ER export, cells ensure that septum formation is contingent upon chromosome separation and exit from mitosis.

Systematic phosphorylation analysis of human mitotic protein complexes

Progression through mitosis depends on a large number of protein complexes that regulate the major structural and physiological changes necessary for faithful chromosome segregation. Most, if not all, of the mitotic processes are regulated by a set of mitotic protein kinases that control protein activity by phosphorylation. Although many mitotic phosphorylation events have been identified in proteome-scale mass spectrometry studies, information on how these phosphorylation sites are distributed within mitotic protein complexes and which kinases generate these phosphorylation sites is largely lacking. We used systematic protein-affinity purification combined with mass spectrometry to identify 1818 phosphorylation sites in more than 100 mitotic protein complexes. In many complexes, the phosphorylation sites were concentrated on a few subunits, suggesting that these subunits serve as "switchboards" to relay the kinase-regulatory signals within the complexes. Consequent bioinformatic analyses identified potential kinase-substrate relationships for most of these sites. In a subsequent in-depth analysis of key mitotic regulatory complexes with the Aurora kinase B (AURKB) inhibitor Hesperadin and a new Polo-like kinase (PLK1) inhibitor, BI 4834, we determined the kinase dependency for 172 phosphorylation sites on 41 proteins. Combination of the results of the cellular studies with Scansite motif prediction enabled us to identify 14 sites on six proteins as direct candidate substrates of AURKB or PLK1.

The quantitative proteome of a human cell line

The majority of all proteins expressed in the human osteosarcoma cell line U2OS were absolutely quantified by mass spectrometry. The quantified proteins span a concentration range of seven orders of magnitude up to 20 000 000 copies per cell.

The generation of mathematical models of biological processes, the simulation of these processes under different conditions, and the comparison and integration of multiple data sets are explicit goals of systems biology that require the knowledge of the absolute quantity of the system's components. To date, systematic estimates of cellular protein concentrations have been exceptionally scarce. Here, we provide a quantitative description of the proteome of a commonly used human cell line in two functional states, interphase and mitosis. We show that these human cultured cells express at least ∼10 000 proteins and that the quantified proteins span a concentration range of seven orders of magnitude up to 20 000 000 copies per cell. We discuss how protein abundance is linked to function and evolution.

Phosphorylation provides a negative mode of regulation for the yeast Rab GTPase Sec4p

The Rab family of Ras-related GTPases are part of a complex signaling circuitry in eukaryotic cells, yet we understand little about the mechanisms that underlie Rab protein participation in such signal transduction networks, or how these networks are integrated at the physiological level. Reversible protein phosphorylation is widely used by cells as a signaling mechanism. Several phospho-Rabs have been identified, however the functional consequences of the modification appear to be diverse and need to be evaluated on an individual basis. In this study we demonstrate a role for phosphorylation as a negative regulatory event for the action of the yeast Rab GTPase Sec4p in regulating polarized growth. Our data suggest that the phosphorylation of the Rab Sec4p prevents interactions with its effector, the exocyst component Sec15p, and that the inhibition may be relieved by a PP2A phosphatase complex containing the regulatory subunit Cdc55p.

The Saccharomyces cerevisiae fermentation stress response protein Igd1p/Yfr017p regulates glycogen levels by inhibiting the glycogen debranching enzyme

Wine fermentation imposes a number of stresses on Saccharomyces cerevisiae, and wine yeasts respond to this harsh environment by altering their transcriptional profile (Marks et al., 2008). We have labeled this change in gene expression patterns the fermentation stress response (FSR). An important component of the FSR is the increased expression of 62 genes for which no function has been identified for their protein products. We hypothesize that a function for these proteins may only be revealed late in grape must fermentation, when the yeast cells are facing conditions much more extreme than those normally encountered in laboratory media. We used affinity copurification to identify interaction partners for the FSR protein Yfr017p, and found that it interacts specifically with the glycogen debranching enzyme (Gdb1p). The expression of both of these proteins is strongly induced during wine fermentation. Therefore, we investigated the role of Yfr017p in glycogen metabolism by constructing wine yeast strains that lack this protein. These YFR017C null cells displayed a significant reduction in their ability to accumulate glycogen during aerobic growth and fermentation. Moreover, Yfr017p inhibits Gdb1p activity in vitro. These results suggest that Yfr017p functions as an inhibitor of Gdb1p, enhancing the ability of yeast cells to store glucose as glycogen. Therefore, we propose IGD1 (for inhibitor of glycogen debranching) as a gene name for the YFR017C ORF.

Regulation and control of metabolic fluxes in microbes

After about ten years of research renaissance in metabolism, the present challenge is to understand how metabolic fluxes are controlled by a complex interplay of overlapping regulatory mechanisms. Reconstruction of various regulatory network topologies is steaming, illustrating that we underestimated the broad importance of post-translational modifications such as enzyme phosphorylation or acetylation for microbial metabolism. With the growing topological knowledge, the functional relevance of these regulatory events becomes an even more pressing need. A major knowledge gap resides in the regulatory network of protein-metabolite interactions, simply because we lacked pertinent methods for systematic analyses - but a start has now been made. Perhaps most dramatic was the conceptual shift in our perception of metabolism from an engine of cellular operation to a generator of input and feedback signals for regulatory circuits that govern many important decisions on cell proliferation, differentiation, death, and naturally metabolism.

Network-based tools for the identification of novel drug targets

In the past few years, network-based tools have become increasingly important in the identification of novel molecular targets for drug development. Systems-based approaches to predict signal transduction-related drug targets have developed into an especially promising field. Here, we summarize our studies, which indicate that modular bridges and overlaps of protein-protein interaction and signaling networks may be of key importance in future drug design. Intermodular nodes are very efficient in mediating the transmission of perturbations between signaling modules and are important in network cooperation. The analysis of stress-induced rearrangements of the yeast interactome by the ModuLand modularization algorithm indicated that components of modular overlap are key players in cellular adaptation to stress. Signaling crosstalk was much more pronounced in humans than in Caenorhabditis elegans or Drosophila melanogaster in the SignaLink (http://www.SignaLink.org) database, a uniformly curated database of eight major signaling pathways. We also showed that signaling proteins that participate in multiple pathways included multiple established drug targets and drug target candidates. Lastly, we caution that the pervasive overlap of cellular network modules implies that wider use of multitarget drugs to partially inhibit multiple individual proteins will be necessary to modify specific cellular functions, because targeting single proteins for complete disruption usually affects multiple cellular functions with little specificity for a particular process. Tools for analyzing network topology and especially network dynamics have great potential to identify alternative sets of targets for developing multitarget drugs.

ResponseNet: revealing signaling and regulatory networks linking genetic and transcriptomic screening data

Cellular response to stimuli is typically complex and involves both regulatory and metabolic processes. Large-scale experimental efforts to identify components of these processes often comprise of genetic screening and transcriptomic profiling assays. We previously established that in yeast genetic screens tend to identify response regulators, while transcriptomic profiling assays tend to identify components of metabolic processes. ResponseNet is a network-optimization approach that integrates the results from these assays with data of known molecular interactions. Specifically, ResponseNet identifies a high-probability sub-network, composed of signaling and regulatory molecular interaction paths, through which putative response regulators may lead to the measured transcriptomic changes. Computationally, this is achieved by formulating a minimum-cost flow optimization problem and solving it efficiently using linear programming tools. The ResponseNet web server offers a simple interface for applying ResponseNet. Users can upload weighted lists of proteins and genes and obtain a sparse, weighted, molecular interaction sub-network connecting their data. The predicted sub-network and its gene ontology enrichment analysis are presented graphically or as text. Consequently, the ResponseNet web server enables researchers that were previously limited to separate analysis of their distinct, large-scale experiments, to meaningfully integrate their data and substantially expand their understanding of the underlying cellular response. ResponseNet is available at http://bioinfo.bgu.ac.il/respnet.

Setting the standards for signal transduction research

Major advances in high-throughput technology platforms, coupled with increasingly sophisticated computational methods for systematic data analysis, have provided scientists with tools to better understand the complexity of signaling networks. In this era of massive and diverse data collection, standardization efforts that streamline data gathering, analysis, storage, and sharing are becoming a necessity. Here, we give an overview of current technologies to study signal transduction. We argue that along with the opportunities the new technologies open, their heterogeneous nature poses critical challenges for data handling that are further increased when data are to be integrated in mathematical models. Efficient standardization through markup languages and data annotation is a sine qua non condition for a systems-level analysis of signaling processes. It remains to be seen the extent to which and the speed at which the emerging standardization efforts will be embraced by the signaling community.

Science Signaling Podcast: 21 December 2010

Targeted removal of individual enzymes elicits changes throughout the entire network of kinases and phosphatases in yeast.