Regulatory evolution in proteins by turnover and lineage-specific changes of cyclin-dependent kinase consensus sites

Evolution of SLiM-mediated hijack functions in intrinsically disordered viral proteins

Viruses and their hosts are involved in an ‘arms race’ where they continually evolve mechanisms to overcome each other. It has long been proposed that intrinsic disorder provides a substrate for the evolution of viral hijack functions and that short linear motifs (SLiMs) are important players in this process. Here, we review evidence in support of this tenet from two model systems: the papillomavirus E7 protein and the adenovirus E1A protein. Phylogenetic reconstructions reveal that SLiMs appear and disappear multiple times across evolution, providing evidence of convergent evolution within individual viral phylogenies. Multiple functionally related SLiMs show strong coevolution signals that persist across long distances in the primary sequence and occur in unrelated viral proteins. Moreover, changes in SLiMs are associated with changes in phenotypic traits such as host range and tropism. Tracking viral evolutionary events reveals that host switch events are associated with the loss of several SLiMs, suggesting that SLiMs are under functional selection and that changes in SLiMs support viral adaptation. Fine-tuning of viral SLiM sequences can improve affinity, allowing them to outcompete host counterparts. However, viral SLiMs are not always competitive by themselves, and tethering of two suboptimal SLiMs by a disordered linker may instead enable viral hijack. Coevolution between the SLiMs and the linker indicates that the evolution of disordered regions may be more constrained than previously thought. In summary, experimental and computational studies support a role for SLiMs and intrinsic disorder in viral hijack functions and in viral adaptive evolution.

A functionally divergent intrinsically disordered region underlying the conservation of stochastic signaling

Stochastic signaling dynamics expand living cells' information processing capabilities. An increasing number of studies report that regulators encode information in their pulsatile dynamics. The evolutionary mechanisms that lead to complex signaling dynamics remain uncharacterized, perhaps because key interactions of signaling proteins are encoded in intrinsically disordered regions (IDRs), whose evolution is difficult to analyze. Here we focused on the IDR that controls the stochastic pulsing dynamics of Crz1, a transcription factor in fungi downstream of the widely conserved calcium signaling pathway. We find that Crz1 IDRs from anciently diverged fungi can all respond transiently to calcium stress; however, only Crz1 IDRs from the Saccharomyces clade support pulsatility, encode extra information, and rescue fitness in competition assays, while the Crz1 IDRs from distantly related fungi do none of the three. On the other hand, we find that Crz1 pulsing is conserved in the distantly related fungi, consistent with the evolutionary model of stabilizing selection on the signaling phenotype. Further, we show that a calcineurin docking site in a specific part of the IDRs appears to be sufficient for pulsing and show evidence for a beneficial increase in the relative calcineurin affinity of this docking site. We propose that evolutionary flexibility of functionally divergent IDRs underlies the conservation of stochastic signaling by stabilizing selection.

The CRK2-CYC13 complex functions as an S-phase cyclin-dependent kinase to promote DNA replication in Trypanosoma brucei

BACKGROUND

Faithful DNA replication is essential to maintain genomic stability in all living organisms, and the regulatory pathway for DNA replication initiation is conserved from yeast to humans. The evolutionarily ancient human parasite Trypanosoma brucei, however, lacks many of the conserved DNA replication factors and may employ unusual mechanisms for DNA replication. Neither the S-phase cyclin-dependent kinase (CDK) nor the regulatory pathway governing DNA replication has been previously identified in T. brucei.

RESULTS

Here we report that CRK2 (Cdc2-related kinase 2) complexes with CYC13 (Cyclin13) and functions as an S-phase CDK to promote DNA replication in T. brucei. We further show that CRK2 phosphorylates Mcm3, a subunit of the Mcm2-7 sub-complex of the Cdc45-Mcm2-7-GINS complex, and demonstrate that Mcm3 phosphorylation by CRK2 facilitates interaction with Sld5, a subunit of the GINS sub-complex of the Cdc45-Mcm2-7-GINS complex.

CONCLUSIONS

These results identify the CRK2-CYC13 complex as an S-phase regulator in T. brucei and reveal its role in regulating DNA replication through promoting the assembly of the Cdc45-Mcm2-7-GINS complex.

The evolutionary scaling of cellular traits imposed by the drift barrier

Owing to internal homeostatic mechanisms, cellular traits may experience long periods of stable selective pressures, during which the stochastic forces of drift and mutation conspire to generate variation. However, even in the face of invariant selection, the drift barrier defined by the genetic effective population size, which is negatively associated with organism size, can have a substantial influence on the location and dispersion of the long-term steady-state distribution of mean phenotypes. In addition, for multilocus traits, the multiplicity of alternative, functionally equivalent states can draw mean phenotypes away from selective optima, even in the absence of mutation bias. Using a framework for traits with an additive genetic basis, it is shown that 1) optimal phenotypic states may be only rarely achieved; 2) gradients of mean phenotypes with respect to organism size (i.e., allometric relationships) are likely to be molded by differences in the power of random genetic drift across the tree of life; and 3) for any particular set of population-genetic conditions, significant variation in mean phenotypes may exist among lineages exposed to identical selection pressures. These results provide a potentially useful framework for understanding numerous aspects of cellular diversification and illustrate the risks of interpreting such variation in a purely adaptive framework.

A Theoretical Framework for Evolutionary Cell Biology

One of the last uncharted territories in evolutionary biology concerns the link with cell biology. Because all phenotypes ultimately derive from events at the cellular level, this connection is essential to building a mechanism-based theory of evolution. Given the impressive developments in cell biological methodologies at the structural and functional levels, the potential for rapid progress is great. The primary challenge for theory development is the establishment of a quantitative framework that transcends species boundaries. Two approaches to the problem are presented here: establishing the long-term steady-state distribution of mean phenotypes under specific regimes of mutation, selection, and drift and evaluating the energetic costs of cellular structures and functions. Although not meant to be the final word, these theoretical platforms harbor potential for generating insight into a diversity of unsolved problems, ranging from genome structure to cellular architecture to aspects of motility in organisms across the Tree of Life.

Evidence for Stabilizing Selection Driving Mutational Turnover of Short Motifs in the Eukaryotic Complementary Sex Determiner (Csd) Protein

Short linear motifs (SLiMs) can play pivotal functional roles in proteins, such as targeting proteins to specific subcellular localizations, modulating the efficiency of translation and tagging proteins for degradation. Until recently we had little knowledge about SLiM evolution. Only a few amino acids in these motifs are functionally important, making them likely to evolve ex nihilo and suggesting that they can play key roles in protein evolution. Several reports now suggest that these motifs can appear and disappear while their function in the protein is preserved, a process sometimes referred to as “turnover”. However, there has been a lack of specific experiments to determine whether independently evolved motifs do indeed have the same function, which would conclusively determine whether the process of turnover actually occurs. In this study, we experimentally detected evidence for such a mutational turnover process for nuclear localization signals (NLS) during the post-duplication divergence of the Complementary sex determiner (Csd) and Feminizer (Fem) proteins in the honeybee (Apis mellifera) lineage. Experiments on the nuclear transport activity of protein segments and those of the most recent common ancestor (MRCA) sequences revealed that three new NLS motifs evolved in the Csd protein during the post-duplication divergence while other NLS motifs were lost that existed before duplication. A screen for essential and newly evolved amino acids revealed that new motifs in the Csd protein evolved by one or two missense mutations coding for lysine. Amino acids that were predating the duplication were also essential in the acquisition of the C1 motif suggesting that the ex nihilo origin was constrained by preexisting amino acids in the physical proximity. Our data support a model in which stabilizing selection maintains the constancy of nuclear transport function but allowed mutational turnover of the encoding NLS motifs.

Conservation of kinase-phosphorylation site pairings: Evidence for an evolutionarily dynamic phosphoproteome

Kinase-mediated protein phosphorylation is a central mechanism for regulation of cellular responses and phenotypes. While considerable information is available regarding the evolutionary relationships within the kinase family, as well as the evolutionary conservation of phosphorylation sites, each aspect of this partnership is typically considered in isolation, despite their clear functional relationship. Here, to offer a more holistic perspective on the evolution of protein phosphorylation, the conservation of protein phosphorylation sites is considered in the context of the conservation of the corresponding modifying kinases. Specifically, conservation of defined kinase-phosphorylation site pairings (KPSPs), as well as of each of the component parts (the kinase and the phosphorylation site), were examined across a range of species. As expected, greater evolutionary distance between species was generally associated with lower probability of KPSP conservation, and only a small fraction of KPSPs were maintained across all species, with the vast majority of KPSP losses due to the absence of the phosphorylation site. This supports a model in which a relatively stable kinome promotes the emergence of functional substrates from an evolutionarily malleable phosphoproteome.

Loss of conserved ubiquitylation sites in conserved proteins during human evolution

Ubiquitylation of lysine residues in proteins serves a pivotal role in the efficient removal of misfolded or unused proteins and in the control of various regulatory pathways by monitoring protein activity that may lead to protein degradation. The loss of ubiquitylated lysines may affect the ubiquitin‑mediated regulatory network and result in the emergence of novel phenotypes. The present study analyzed mouse ubiquitylation data and orthologous proteins from 62 mammals to identify 193 conserved ubiquitylation sites from 169 proteins that were lost in the Euarchonta lineage leading to humans. A total of 8 proteins, including betaine homocysteine S‑methyltransferase, clin and CBS domain divalent metal cation transport mediator 3, ribosome‑binding protein 1 and solute carrier family 37 member 4, lost 1 conserved lysine residue, which was ubiquitylated in the mouse ortholog, following the human‑chimpanzee divergence. A total of 17 of the lost ubiquitylated lysines are also known to be modified by acetylation and/or succinylation in mice. In 8 cases, a novel lysine evolved at positions flanking the lost conserved lysine residues, potentially as a method of compensation. We hypothesize that the loss of ubiquitylation sites during evolution may lead to the development of advantageous phenotypes, which are then fixed by selection. The ancestral ubiquitylation sites identified in the present study may be a useful resource for investigating the association between loss of ubiquitylation sites and the emergence of novel phenotypes during evolution towards modern humans.

Most m6A RNA Modifications in Protein-Coding Regions Are Evolutionarily Unconserved and Likely Nonfunctional

Methylation of the adenosine base at the nitrogen-6 position (m6A) is the most prevalent internal posttranscriptional modification of mRNAs in many eukaryotes. Despite the rapid progress in the transcriptome-wide mapping of m6As, identification of proteins responsible for writing, reading, and erasing m6As, and elucidation of m6A functions in splicing, RNA stability, translation, and other processes, it is unknown whether most observed m6A modifications are functional. To address this question, we respectively analyze the evolutionary conservation of yeast and human m6As in protein-coding regions. Relative to comparable unmethylated As, m6As are overall no more conserved in yeasts and only slightly more conserved in mammals. Furthermore, yeast m6As and comparable unmethylated As have no significant difference in single nucleotide polymorphism (SNP) density or SNP site frequency spectrum. The same is true in human. The methylation status of a gene, not necessarily the specific sites methylated in the gene, is subject to purifying selection for no more than ∼20% of m6A-modified genes. These observations suggest that most m6A modifications in protein-coding regions are nonfunctional and nonadaptive, probably resulting from off-target activities of m6A methyltransferases. In addition, our reanalysis invalidates the recent claim of positive selection for newly acquired m6A modifications in human evolution. Regarding the small number of evolutionarily conserved m6As, evidence suggests that a large proportion of them are likely functional; they should be prioritized in future functional characterizations of m6As. Together, these findings have important implications for understanding the biological significance of m6A and other posttranscriptional modifications.

Multivalent Interactions with Fbw7 and Pin1 Facilitate Recognition of c-Jun by the SCFFbw7 Ubiquitin Ligase

Many regulatory proteins, including the transcription factor c-Jun, are highly enriched in disordered protein regions that govern growth, division, survival, differentiation, and response to signals. The stability of c-Jun is controlled by poorly understood regulatory interactions of its disordered region with both the E3 ubiquitin ligase SCF^Fbw7 and prolyl cis-trans isomerase Pin1. We use nuclear magnetic resonance and fluorescence studies of c-Jun to demonstrate that multisite c-Jun phosphorylation is required for high-affinity interaction with Fbw7. We show that the Pin1 WW and PPIase domains interact in a dynamic complex with multiply phosphorylated c-Jun. Importantly, Pin1 isomerizes a pSer-Pro peptide bond at the c-Jun N terminus that affects binding to Fbw7 and thus modulates the ubiquitin-mediated degradation of c-Jun. Our findings support the general principle that multiple weak binding motifs within disordered regions can synergize to yield high-affinity interactions and provide rapidly evolvable means to build and fine-tune regulatory events.

Co-occurring protein phosphorylation are functionally associated

Post-translational modifications (PTMs) add a further layer of complexity to the proteome and regulate a wide range of cellular protein functions. With the increasing number of known PTM sites, it becomes imperative to understand their functional interplays. In this study, we proposed a novel analytical strategy to explore functional relationships between PTM sites by testing their tendency to be modified together (co-occurrence) under the same condition, and applied it to proteome-wide human phosphorylation data collected under 88 different laboratory or physiological conditions. Co-occurring phosphorylation occurs significantly more frequently than randomly expected and include many known examples of cross-talk or functional connections. Such pairs, either within the same phosphoprotein or between interacting partners, are more likely to be in sequence or structural proximity, be phosphorylated by the same kinases, participate in similar biological processes, and show residue co-evolution across vertebrates. In addition, we also found that their co-occurrence states tend to be conserved in orthologous phosphosites in the mouse proteome. Together, our results support that the co-occurring phosphorylation are functionally associated. Comparison with existing methods further suggests that co-occurrence analysis can be a useful complement to uncover novel functional associations between PTM sites.

In addition to gene expression and translation control, post-translational modifications (PTMs) represent another level to regulate proteins functions. Different PTM sites within a protein usually co-operate to fulfill their functional roles. Recent advances in high-throughput mass spectrometry (MS) technologies have facilitated the proteome-wide identification of PTM sites, giving rise to both challenge and opportunity to understand their functional relationships. Previously, several data mining approaches have been developed to explore the global PTM interplays. In this study, we proposed to infer functional associations between PTM sites from the correlation of their modification status across many biological conditions, which was not exploited before. In practice, we tested if a pair of sites are modified together under the same condition significantly more often than expected (co-occurrence). As a proof of principle, we applied this analytical strategy to human phosphorylation because we could collect data sets of proteome-wide coverage under 88 different conditions. We demonstrated that sites with co-occurring phosphorylation status are functionally associated from several lines of evidence. The co-occurrence analysis can also uncover functionally connected phosphosites with clear biological evidence which are missed by other approaches. With increasing proteome-wide data for other types of PTMs under different conditions, the co-occurrence analysis can be integrated with other methods to identify novel PTM associations.

Selection maintains signaling function of a highly diverged intrinsically disordered region

Intrinsically disordered regions (IDRs) are characterized by their lack of stable secondary or tertiary structure and comprise a large part of the eukaryotic proteome. Although these regions play a variety of signaling and regulatory roles, they appear to be rapidly evolving at the primary sequence level. To understand the functional implications of this rapid evolution, we focused on a highly diverged IDR in Saccharomyces cerevisiae that is involved in regulating multiple conserved MAPK pathways. We hypothesized that under stabilizing selection, the functional output of orthologous IDRs could be maintained, such that diverse genotypes could lead to similar function and fitness. Consistent with the stabilizing selection hypothesis, we find that diverged, orthologous IDRs can mostly recapitulate wild-type function and fitness in S. cerevisiae We also find that the electrostatic charge of the IDR is correlated with signaling output and, using phylogenetic comparative methods, find evidence for selection maintaining this quantitative molecular trait despite underlying genotypic divergence.

Computational Analysis of the Predicted Evolutionary Conservation of Human Phosphorylation Sites

Protein kinase-mediated phosphorylation is among the most important post-translational modifications. However, few phosphorylation sites have been experimentally identified for most species, making it difficult to determine the degree to which phosphorylation sites are conserved. The goal of this study was to use computational methods to characterize the conservation of human phosphorylation sites in a wide variety of eukaryotes. Using experimentally-determined human sites as input, homologous phosphorylation sites were predicted in all 432 eukaryotes for which complete proteomes were available. For each pair of species, we calculated phosphorylation site conservation as the number of phosphorylation sites found in both species divided by the number found in at least one of the two species. A clustering of the species based on this conservation measure was concordant with phylogenies based on traditional genomic measures. For a subset of the 432 species, phosphorylation site conservation was compared to conservation of both protein kinases and proteins in general. Protein kinases exhibited the highest degree of conservation, while general proteins were less conserved and phosphorylation sites were least conserved. Although preliminary, these data tentatively suggest that variation in phosphorylation sites may play a larger role in explaining phenotypic differences among organisms than differences in the complements of protein kinases or general proteins.

Motif co-regulation and co-operativity are common mechanisms in transcriptional, post-transcriptional and post-translational regulation

A substantial portion of the regulatory interactions in the higher eukaryotic cell are mediated by simple sequence motifs in the regulatory segments of genes and (pre-)mRNAs, and in the intrinsically disordered regions of proteins. Although these regulatory modules are physicochemically distinct, they share an evolutionary plasticity that has facilitated a rapid growth of their use and resulted in their ubiquity in complex organisms. The ease of motif acquisition simplifies access to basal housekeeping functions, facilitates the co-regulation of multiple biomolecules allowing them to respond in a coordinated manner to changes in the cell state, and supports the integration of multiple signals for combinatorial decision-making. Consequently, motifs are indispensable for temporal, spatial, conditional and basal regulation at the transcriptional, post-transcriptional and post-translational level. In this review, we highlight that many of the key regulatory pathways of the cell are recruited by motifs and that the ease of motif acquisition has resulted in large networks of co-regulated biomolecules. We discuss how co-operativity allows simple static motifs to perform the conditional regulation that underlies decision-making in higher eukaryotic biological systems. We observe that each gene and its products have a unique set of DNA, RNA or protein motifs that encode a regulatory program to define the logical circuitry that guides the life cycle of these biomolecules, from transcription to degradation. Finally, we contrast the regulatory properties of protein motifs and the regulatory elements of DNA and (pre-)mRNAs, advocating that co-regulation, co-operativity, and motif-driven regulatory programs are common mechanisms that emerge from the use of simple, evolutionarily plastic regulatory modules.

Short linear motifs - ex nihilo evolution of protein regulation

Short sequence motifs are ubiquitous across the three major types of biomolecules: hundreds of classes and thousands of instances of DNA regulatory elements, RNA motifs and protein short linear motifs (SLiMs) have been characterised. The increase in complexity of transcriptional, post-transcriptional and post-translational regulation in higher Eukaryotes has coincided with a significant expansion of motif use. But how did the eukaryotic cell acquire such a vast repertoire of motifs? In this review, we curate the available literature on protein motif evolution and discuss the evidence that suggests SLiMs can be acquired by mutations, insertions and deletions in disordered regions. We propose a mechanism of ex nihilo SLiM evolution – the evolution of a novel SLiM from “nothing” – adding a functional module to a previously non-functional region of protein sequence. In our model, hundreds of motif-binding domains in higher eukaryotic proteins connect simple motif specificities with useful functions to create a large functional motif space. Accessible peptides that match the specificity of these motif-binding domains are continuously created and destroyed by mutations in rapidly evolving disordered regions, creating a dynamic supply of new interactions that may have advantageous phenotypic novelty. This provides a reservoir of diversity to modify existing interaction networks. Evolutionary pressures will act on these motifs to retain beneficial instances. However, most will be lost on an evolutionary timescale as negative selection and genetic drift act on deleterious and neutral motifs respectively. In light of the parallels between the presented model and the evolution of motifs in the regulatory segments of genes and (pre-)mRNAs, we suggest our understanding of regulatory networks would benefit from the creation of a shared model describing the evolution of transcriptional, post-transcriptional and post-translational regulation.

Prediction of Functionally Important Phospho-Regulatory Events in Xenopus laevis Oocytes

The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.

Bora downregulation results in radioresistance by promoting repair of double strand breaks

Following DNA double-strand breaks cells activate several DNA-damage response protein kinases, which then trigger histone H2AX phosphorylation and the accumulation of proteins such as MDC1, p53-binding protein 1, and breast cancer gene 1 at the damage site to promote DNA double-strand breaks repair. We identified a novel biomarker, Bora (previously called C13orf34), that is associated with radiosensitivity. In the current study, we set out to investigate how Bora might be involved in response to irradiation. We found a novel function of Bora in DNA damage repair response. Bora down-regulation increased colony formation in cells exposed to irradiation. This increased resistance to irradiation in Bora-deficient cells is likely due to a faster rate of double-strand breaks repair. After irradiation, Bora-knockdown cells displayed increased G2-M cell cycle arrest and increased Chk2 phosphorylation. Furthermore, Bora specifically interacted with the tandem breast cancer gene 1 C-terminal domain of MDC1 in a phosphorylation dependent manner, and overexpression of Bora could abolish irradiation induced MDC1 foci formation. In summary, Bora may play a significant role in radiosensitivity through the regulation of MDC1 and DNA repair.

Evolution of bacterial protein-tyrosine kinases and their relaxed specificity toward substrates

It has often been speculated that bacterial protein-tyrosine kinases (BY-kinases) evolve rapidly and maintain relaxed substrate specificity to quickly adopt new substrates when evolutionary pressure in that direction arises. Here, we report a phylogenomic and biochemical analysis of BY-kinases, and their relationship to substrates aimed to validate this hypothesis. Our results suggest that BY-kinases are ubiquitously distributed in bacterial phyla and underwent a complex evolutionary history, affected considerably by gene duplications and horizontal gene transfer events. This is consistent with the fact that the BY-kinase sequences represent a high level of substitution saturation and have a higher evolutionary rate compared with other bacterial genes. On the basis of similarity networks, we could classify BY kinases into three main groups with 14 subgroups. Extensive sequence conservation was observed only around the three canonical Walker motifs, whereas unique signatures proposed the functional speciation and diversification within some subgroups. The relationship between BY-kinases and their substrates was analyzed using a ubiquitous substrate (Ugd) and some Firmicute-specific substrates (YvyG and YjoA) from Bacillus subtilis. No evidence of coevolution between kinases and substrates at the sequence level was found. Seven BY-kinases, including well-characterized and previously uncharacterized ones, were used for experimental studies. Most of the tested kinases were able to phosphorylate substrates from B. subtilis (Ugd, YvyG, and YjoA), despite originating from very distant bacteria. Our results are consistent with the hypothesis that BY-kinases have evolved relaxed substrate specificity and are probably maintained as rapidly evolving platforms for adopting new substrates.

Evolutionary constraint and disease associations of post-translational modification sites in human genomes

Interpreting the impact of human genome variation on phenotype is challenging. The functional effect of protein-coding variants is often predicted using sequence conservation and population frequency data, however other factors are likely relevant. We hypothesized that variants in protein post-translational modification (PTM) sites contribute to phenotype variation and disease. We analyzed fraction of rare variants and non-synonymous to synonymous variant ratio (Ka/Ks) in 7,500 human genomes and found a significant negative selection signal in PTM regions independent of six factors, including conservation, codon usage, and GC-content, that is widely distributed across tissue-specific genes and function classes. PTM regions are also enriched in known disease mutations, suggesting that PTM variation is more likely deleterious. PTM constraint also affects flanking sequence around modified residues and increases around clustered sites, indicating presence of functionally important short linear motifs. Using target site motifs of 124 kinases, we predict that at least ∼180,000 motif-breaker amino acid residues that disrupt PTM sites when substituted, and highlight kinase motifs that show specific negative selection and enrichment of disease mutations. We provide this dataset with corresponding hypothesized mechanisms as a community resource. As an example of our integrative approach, we propose that PTPN11 variants in Noonan syndrome aberrantly activate the protein by disrupting an uncharacterized cluster of phosphorylation sites. Further, as PTMs are molecular switches that are modulated by drugs, we study mutated binding sites of PTM enzymes in disease genes and define a drug-disease network containing 413 novel predicted disease-gene links.

Evolutionary conservation of the CDK targets in eukaryotic DNA replication initiation

A fundamental requirement for all organisms is the faithful duplication and transmission of the genetic material. Failure to accurately copy and segregate the genome during cell division leads to loss of genetic information and chromosomal abnormalities. Such genome instability is the hallmark of the earliest stages of tumour formation. Cyclin-dependent kinase (CDK) plays a vital role in regulating the duplication of the genome within the eukaryotic cell cycle. Importantly, this kinase is deregulated in many cancer types and is an emerging target of chemotherapeutics. In this review, I will consider recent advances concerning the role of CDK in replication initiation across eukaryotes. The implications for strict CDK-dependent regulation of genome duplication in the context of the cell cycle will be discussed.

Detecting functional divergence after gene duplication through evolutionary changes in posttranslational regulatory sequences

Gene duplication is an important evolutionary mechanism that can result in functional divergence in paralogs due to neo-functionalization or sub-functionalization. Consistent with functional divergence after gene duplication, recent studies have shown accelerated evolution in retained paralogs. However, little is known in general about the impact of this accelerated evolution on the molecular functions of retained paralogs. For example, do new functions typically involve changes in enzymatic activities, or changes in protein regulation? Here we study the evolution of posttranslational regulation by examining the evolution of important regulatory sequences (short linear motifs) in retained duplicates created by the whole-genome duplication in budding yeast. To do so, we identified short linear motifs whose evolutionary constraint has relaxed after gene duplication with a likelihood-ratio test that can account for heterogeneity in the evolutionary process by using a non-central chi-squared null distribution. We find that short linear motifs are more likely to show changes in evolutionary constraints in retained duplicates compared to single-copy genes. We examine changes in constraints on known regulatory sequences and show that for the Rck1/Rck2, Fkh1/Fkh2, Ace2/Swi5 paralogs, they are associated with previously characterized differences in posttranslational regulation. Finally, we experimentally confirm our prediction that for the Ace2/Swi5 paralogs, Cbk1 regulated localization was lost along the lineage leading to SWI5 after gene duplication. Our analysis suggests that changes in posttranslational regulation mediated by short regulatory motifs systematically contribute to functional divergence after gene duplication.

How a protein is controlled is intimately linked to its function. Therefore, evolution can drive the functional divergence of proteins by tweaking their regulation, even if enzymatic capacities are preserved. Changes in posttranslational regulation (protein phosphorylation, degradation, subcellular localization, etc.) could therefore represent key mechanisms in functional divergence and lead to different phenotypic outcomes. Since disordered protein regions contain sites of protein modification and interaction (known as short linear motifs) and evolve rapidly relative to domains encoding enzymatic functions, these regions are good candidates to harbour sequence changes that underlie changes in function. In this study, we develop a statistical framework to identify changes in rate of evolution specific to protein regulatory sequences and identify hundreds of short linear motifs in disordered regions that are likely to have diverged after the whole-genome duplication in budding yeast. We show that these divergent motifs are much more frequent in paralogs than in single-copy proteins, and that they are more frequent in duplicate pairs that have functionally diverged. Our analysis suggests that changes in short linear motifs in disordered protein regions could be important molecular mechanisms of functional divergence after gene duplication.

Turnover of protein phosphorylation evolving under stabilizing selection

Most proteins are regulated by posttranslational modifications and changes in these modifications contribute to evolutionary changes as well as to human diseases. Phosphorylation of serines, threonines, and tyrosines are the most common modifications identified to date in eukaryotic proteomes. While the mode of action and the function of most phosphorylation sites remain unknown, functional studies have shown that phosphorylation affects protein stability, localization and ability to interact. Two broad modes of action have been described for protein phosphorylation. The first mode corresponds to the canonical and qualitative view whereby single phosphorylation sites act as molecular switches that either turn on or off specific protein functions through direct or allosteric effects. The second mode is more akin to a rheostat than a switch. In this case, a group of phosphorylation sites in a given protein region contributes collectively to the modification of the protein, irrespective of the precise position of individual sites, through an aggregate property. Here we discuss these two types of regulation and examine how they affect the rate and patterns of protein phosphorylation evolution. We describe how the evolution of clusters of phosphorylation sites can be studied under the framework of complex traits evolution and stabilizing selection.

CDK phosphorylation of SLD-2 is required for replication initiation and germline development in C. elegans

Timely phosphorylation of SLD-2 by CDK is essential for proper replication initiation and cell proliferation in the germline of C. elegans.

Cyclin-dependent kinase (CDK) plays a vital role in proliferation control across eukaryotes. Despite this, how CDK mediates cell cycle and developmental transitions in metazoa is poorly understood. In this paper, we identify orthologues of Sld2, a CDK target that is important for DNA replication in yeast, and characterize SLD-2 in the nematode worm Caenorhabditis elegans. We demonstrate that SLD-2 is required for replication initiation and the nuclear retention of a critical component of the replicative helicase CDC-45 in embryos. SLD-2 is a CDK target in vivo, and phosphorylation regulates the interaction with another replication factor, MUS-101. By mutation of the CDK sites in sld-2, we show that CDK phosphorylation of SLD-2 is essential in C. elegans. Finally, using a phosphomimicking sld-2 mutant, we demonstrate that timely CDK phosphorylation of SLD-2 is an important control mechanism to allow normal proliferation in the germline. These results determine an essential function of CDK in metazoa and identify a developmental role for regulated SLD-2 phosphorylation.

Functional divergence and evolutionary turnover in mammalian phosphoproteomes

Protein phosphorylation is a key mechanism to regulate protein functions. However, the contribution of this protein modification to species divergence is still largely unknown. Here, we studied the evolution of mammalian phosphoregulation by comparing the human and mouse phosphoproteomes. We found that 84% of the positions that are phosphorylated in one species or the other are conserved at the residue level. Twenty percent of these conserved sites are phosphorylated in both species. This proportion is 2.5 times more than expected by chance alone, suggesting that purifying selection is preserving phosphoregulation. However, we show that the majority of the sites that are conserved at the residue level are differentially phosphorylated between species. These sites likely result from false-negative identifications due to incomplete experimental coverage, false-positive identifications and non-functional sites. In addition, our results suggest that at least 5% of them are likely to be true differentially phosphorylated sites and may thus contribute to the divergence in phosphorylation networks between mouse and humans and this, despite residue conservation between orthologous proteins. We also showed that evolutionary turnover of phosphosites at adjacent positions (in a distance range of up to 40 amino acids) in human or mouse leads to an over estimation of the divergence in phosphoregulation between these two species. These sites tend to be phosphorylated by the same kinases, supporting the hypothesis that they are functionally redundant. Our results support the hypothesis that the evolutionary turnover of phosphorylation sites contributes to the divergence in phosphorylation profiles while preserving phosphoregulation. Overall, our study provides advanced analyses of mammalian phosphoproteomes and a framework for the study of their contribution to phenotypic evolution.

Understanding how differences in cellular regulation lead to phenotypic differences between species remains an open challenge in evolutionary genetics. The extensive phosphorylation data currently available allows to compare the human and mouse phosphoproteomes and to measure changes in their phosphoregulation. We found a general conservation of phosphorylation sites between these two species. However, a fraction of sites are conserved at the sequence level (the same amino acid is present in both species) but differ in their phosphorylation status. These sites represent candidate sites that have the potential to explain differences between human and mouse signalling networks that do not depend on the divergence of orthologous residues. Furthermore, we identified several sites where to a phosphorylation site in one species corresponds a non-phosphorylatable residue in the other one. These cases represent clear differences in protein regulation. Recent studies suggest that phosphorylation sites can shift position during evolution, leading to configurations in which pairs of divergent phosphorylation sites are functionally redundant. We identified more than 100 putative such cases, suggesting that divergence in amino acid does not necessarily imply functional divergence when comparing phosphoproteomes. Overall, our study provides new key concepts and data for the study of how regulatory differences may be linked to phenotypic ones at the network level.

Evolution and functional cross-talk of protein post-translational modifications

Protein post-translational modifications (PTMs) allow the cell to regulate protein activity and play a crucial role in the response to changes in external conditions or internal states. Advances in mass spectrometry now enable proteome wide characterization of PTMs and have revealed a broad functional role for a range of different types of modifications. Here we review advances in the study of the evolution and function of PTMs that were spurred by these technological improvements. We provide an overview of studies focusing on the origin and evolution of regulatory enzymes as well as the evolutionary dynamics of modification sites. Finally, we discuss different mechanisms of altering protein activity via post-translational regulation and progress made in the large-scale functional characterization of PTM function.

Mitotic phosphorylation of Exo84 disrupts exocyst assembly and arrests cell growth

Mitotic phosphorylation of Exo84 disrupts the assembly of the exocyst complex, thereby inhibiting exocytosis of select secreted cargoes and cell surface expansion.

The rate of eukaryotic cell growth is tightly controlled for proper progression through each cell cycle stage and is important for cell size homeostasis. It was previously shown that cell growth is inhibited during mitosis when cells are preparing for division. However, the mechanism for growth arrest at this stage is unknown. Here we demonstrate that exocytosis of a select group of cargoes was inhibited before the metaphase–anaphase transition in the budding yeast Saccharomyces cerevisiae. The cyclin-dependent kinase, Cdk1, when bound to the mitotic cyclin Clb2, directly phosphorylated Exo84, a component of the exocyst complex essential for exocytosis. Mitotic phosphorylation of Exo84 disrupted the assembly of the exocyst complex, thereby affecting exocytosis and cell surface expansion. Our study demonstrates the coordination between membrane trafficking and cell cycle progression and provides a molecular mechanism by which cell growth is controlled during the cell division cycle.

Cell cycle control across the eukaryotic kingdom

Almost two billion years of evolution have generated a vast and amazing variety of eukaryotic life with approximately 8.7 million extant species. Growth and reproduction of all of these organisms depend on faithful duplication and distribution of their chromosomes to the newly forming daughter cells in a process called the cell cycle. However, most of what is known today about cell cycle control comes from a few model species that belong to the unikonts; that is, to only one of five 'supergroups' that comprise the eukaryotic kingdom. Recently, analyzing species from distantly related clades is providing insights into general principles of cell cycle regulation and shedding light on its evolution. Here, referring to animal and fungal as opposed to non-unikont systems, especially flowering plants from the archaeplastid supergroup, we compare the conservation of central cell cycle regulator functions, the structure of network topologies, and the evolutionary dynamics of substrates of core cell cycle kinases.

SH3 interactome conserves general function over specific form

Src homology 3 (SH3) domains bind peptides to mediate protein-protein interactions that assemble and regulate dynamic biological processes. We surveyed the repertoire of SH3 binding specificity using peptide phage display in a metazoan, the worm Caenorhabditis elegans, and discovered that it structurally mirrors that of the budding yeast Saccharomyces cerevisiae. We then mapped the worm SH3 interactome using stringent yeast two-hybrid and compared it with the equivalent map for yeast. We found that the worm SH3 interactome resembles the analogous yeast network because it is significantly enriched for proteins with roles in endocytosis. Nevertheless, orthologous SH3 domain-mediated interactions are highly rewired. Our results suggest a model of network evolution where general function of the SH3 domain network is conserved over its specific form.

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.

Gains of ubiquitylation sites in highly conserved proteins in the human lineage

Background

Post-translational modification of lysine residues of specific proteins by ubiquitin modulates the degradation, localization, and activity of these target proteins. Here, we identified gains of ubiquitylation sites in highly conserved regions of human proteins that occurred during human evolution.

Results

We analyzed human ubiquitylation site data and multiple alignments of orthologous mammalian proteins including those from humans, primates, other placental mammals, opossum, and platypus. In our analysis, we identified 281 ubiquitylation sites in 252 proteins that first appeared along the human lineage during primate evolution: one protein had four novel sites; four proteins had three sites each; 18 proteins had two sites each; and the remaining 229 proteins had one site each. PML, which is involved in neurodevelopment and neurodegeneration, acquired three sites, two of which have been reported to be involved in the degradation of PML. Thirteen human proteins, including ERCC2 (also known as XPD) and NBR1, gained human-specific ubiquitylated lysines after the human-chimpanzee divergence. ERCC2 has a Lys/Gln polymorphism, the derived (major) allele of which confers enhanced DNA repair capacity and reduced cancer risk compared with the ancestral (minor) allele. NBR1 and eight other proteins that are involved in the human autophagy protein interaction network gained a novel ubiquitylation site.

Conclusions

The gain of novel ubiquitylation sites could be involved in the evolution of protein degradation and other regulatory networks. Although gains of ubiquitylation sites do not necessarily equate to adaptive evolution, they are useful candidates for molecular functional analyses to identify novel advantageous genetic modifications and innovative phenotypes acquired during human evolution.

Systematic functional prioritization of protein posttranslational modifications

Protein function is often regulated by posttranslational modifications (PTMs), and recent advances in mass spectrometry have resulted in an exponential increase in PTM identification. However, the functional significance of the vast majority of these modifications remains unknown. To address this problem, we compiled nearly 200,000 phosphorylation, acetylation, and ubiquitination sites from 11 eukaryotic species, including 2,500 newly identified ubiquitylation sites for Saccharomyces cerevisiae. We developed methods to prioritize the functional relevance of these PTMs by predicting those that likely participate in cross-regulatory events, regulate domain activity, or mediate protein-protein interactions. PTM conservation within domain families identifies regulatory "hot spots" that overlap with functionally important regions, a concept that we experimentally validated on the HSP70 domain family. Finally, our analysis of the evolution of PTM regulation highlights potential routes for neutral drift in regulatory interactions and suggests that only a fraction of modification sites are likely to have a significant biological role.

Predicting kinase substrates using conservation of local motif density

MOTIVATION

Protein kinases represent critical links in cell signaling. A central problem in computational biology is to systematically identify their substrates.

RESULTS

This study introduces a new method to predict kinase substrates by extracting evolutionary information from multiple sequence alignments in a manner that is tolerant to degenerate motif positioning. Given a known consensus, the new method (ConDens) compares the observed density of matches to a null model of evolution and does not require labeled training data. We confirmed that ConDens has improved performance compared with several existing methods in the field. Further, we show that it is generalizable and can predict interesting substrates for several important eukaryotic kinases where training data is not available.

AVAILABILITY AND IMPLEMENTATION

ConDens can be found at http://www.moseslab.csb.utoronto.ca/andyl/.

CONTACT

alan.moses@utoronto.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluation and properties of the budding yeast phosphoproteome

We have assembled a reliable phosphoproteomic data set for budding yeast Saccharomyces cerevisiae and have investigated its properties. Twelve publicly available phosphoproteome data sets were triaged to obtain a subset of high-confidence phosphorylation sites (p-sites), free of “noisy” phosphorylations. Analysis of this combined data set suggests that the inventory of phosphoproteins in yeast is close to completion, but that these proteins may have many undiscovered p-sites. Proteins involved in budding and protein kinase activity have high numbers of p-sites and are highly over-represented in the vast majority of the yeast phosphoproteome data sets. The yeast phosphoproteome is characterized by a few proteins with many p-sites and many proteins with a few p-sites. We confirm a tendency for p-sites to cluster together and find evidence that kinases may phosphorylate off-target amino acids that are within one or two residues of their cognate target. This suggests that the precise position of the phosphorylated amino acid is not a stringent requirement for regulatory fidelity. Compared with nonphosphorylated proteins, phosphoproteins are more ancient, more abundant, have longer unstructured regions, have more genetic interactions, more protein interactions, and are under tighter post-translational regulation. It appears that phosphoproteins constitute the raw material for pathway rewiring and adaptation at various evolutionary rates.

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.

Dynamic regulation of Emi2 by Emi2-bound Cdk1/Plk1/CK1 and PP2A-B56 in meiotic arrest of Xenopus eggs

In vertebrates, unfertilized eggs are arrested at metaphase of meiosis II by Mos and Emi2, an inhibitor of the APC/C ubiquitin ligase. In Xenopus, Cdk1 phosphorylates Emi2 and both destabilizes and inactivates it, whereas Mos recruits PP2A phosphatase to antagonize the Cdk1 phosphorylation. However, how Cdk1 phosphorylation inhibits Emi2 is largely unknown. Here we show that multiple N-terminal Cdk1 phosphorylation motifs bind cyclin B1-Cdk1 itself, Plk1, and CK1δ/ε to inhibit Emi2. Plk1, after rebinding to other sites by self-priming phosphorylation, partially destabilizes Emi2. Cdk1 and CK1δ/ε sequentially phosphorylate the C-terminal APC/C-docking site, thereby cooperatively inhibiting Emi2 from binding the APC/C. In the presence of Mos, however, PP2A-B56β/ε bind to Emi2 and keep dephosphorylating it, particularly at the APC/C-docking site. Thus, Emi2 stability and activity are dynamically regulated by Emi2-bound multiple kinases and PP2A phosphatase. Our data also suggest a general role for Cdk1 substrate phosphorylation motifs in M phase regulation.

Identification of novel phosphorylation modification sites in human proteins that originated after the human-chimpanzee divergence

MOTIVATION

Phosphorylation modifications of specific protein residues are involved in a wide range of biological processes such as modulation of intracellular signal networks. Here, we present the development and application of a bioinformatics procedure for systematic identification of human-specific phosphorylation sites in proteins that may have occurred after the human-chimpanzee divergence.

RESULTS

We collected annotated human phosphorylation sites and compared each site to orthologous mammalian proteins across taxa including chimpanzee, orangutan, rhesus macaque, marmoset, mouse, dog, cow, elephant, opossum and platypus. We identified 37 human-specific gains of annotated phosphorylation sites in 35 proteins: 22 serines, 12 threonines and 3 tyrosines. The novel phosphorylation sites are situated in highly conserved segments of the protein. Proteins with novel phosphorylation sites are involved in crucial biological processes such as cell division (AURKB, CASC5, MKI67 and PDCD4) and chromatin remodeling (HIRA, HIRIP3, HIST1H1T, NAP1L4 and LRWD1). Modified phosphorylatable residues produce novel target sites for protein kinases such as cyclin-dependent kinases and casein kinases, possibly resulting in rewiring and fine-tuning of phosphorylation regulatory networks. The potential human-specific phosphorylation sites identified in this study are useful as candidates for functional analysis to identify novel phenotypes in humans.

CONTACT

hahny@cau.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Moving from transcriptional to phospho-evolution: generalizing regulatory evolution?

Much of biological diversity is thought to arise from changes in regulatory networks. Although the role of transcriptional regulation has been well established, the contribution to evolution of changes at other levels of regulation has yet to be addressed. Using examples from the literature and recent studies on the evolution of protein phosphorylation, we argue that protein regulatory networks also play a prime role in generating diversity within and between species. Because there are several analogies between the regulation of protein functions by kinases and the regulation of gene expression by transcription factors, the principles that guide transcriptional regulatory evolution can also be explored in kinase-substrate networks. These comparisons will allow us to generalize existing models of evolution across the complex layers of the cell's regulatory links.

The Mcm complex: unwinding the mechanism of a replicative helicase

The Mcm2-7 complex serves as the eukaryotic replicative helicase, the molecular motor that both unwinds duplex DNA and powers fork progression during DNA replication. Consistent with its central role in this process, much prior work has illustrated that Mcm2-7 loading and activation are landmark events in the regulation of DNA replication. Unlike any other hexameric helicase, Mcm2-7 is composed of six unique and essential subunits. Although the unusual oligomeric nature of this complex has long hampered biochemical investigations, recent advances with both the eukaryotic as well as the simpler archaeal Mcm complexes provide mechanistic insight into their function. In contrast to better-studied homohexameric helicases, evidence suggests that the six Mcm2-7 complex ATPase active sites are functionally distinct and are likely specialized to accommodate the regulatory constraints of the eukaryotic process.

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

Background

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

Results

We present statistically rigorous tests on the spatial and functional properties of a collection of ~70,000 reported phosphosites. We show that the distribution of phosphosites positioning along the protein tends to occur as dense clusters of Serine/Threonines (pS/pT) and between Serine/Threonines and Tyrosines, but generally not as much between Tyrosines (pY) only. This phenomenon is more ubiquitous than anticipated and is pertinent for most eukaryotic proteins: for proteins with ≥ 2 phosphosites, 54% of all pS/pT sites are within 4 amino acids of another site. We found a strong tendency for clustered pS/pT to be activated by the same kinase. Large-scale analyses of phosphopeptides are thus consistent with a cooperative function within the cluster.

Conclusions

We present evidence supporting the notion that clusters of pS/pT but generally not pY should be considered as the elementary building blocks in phosphorylation regulation. Indeed, closely positioned sites tend to be activated by the same kinase, a signal that overrides the tendency of a protein to be activated by a single or only few kinases. Within these clusters, coordination and positional dependency is evident. We postulate that cellular regulation takes advantage of such design. Specifically, phosphosite clusters may increase the robustness of the effectiveness of phosphorylation-dependent response.

Reviewers

Reviewed by Joel Bader, Frank Eisenhaber, Emmanuel Levy (nominated by Sarah Teichmann). For the full reviews, please go to the Reviewers' comments section.

Global analysis of Cdk1 substrate phosphorylation sites provides insights into evolution

To explore the mechanisms and evolution of cell-cycle control, we analyzed the position and conservation of large numbers of phosphorylation sites for the cyclin-dependent kinase Cdk1 in the budding yeast Saccharomyces cerevisiae. We combined specific chemical inhibition of Cdk1 with quantitative mass spectrometry to identify the positions of 547 phosphorylation sites on 308 Cdk1 substrates in vivo. Comparisons of these substrates with orthologs throughout the ascomycete lineage revealed that the position of most phosphorylation sites is not conserved in evolution; instead, clusters of sites shift position in rapidly evolving disordered regions. We propose that the regulation of protein function by phosphorylation often depends on simple nonspecific mechanisms that disrupt or enhance protein-protein interactions. The gain or loss of phosphorylation sites in rapidly evolving regions could facilitate the evolution of kinase-signaling circuits.

Comparative analysis reveals conserved protein phosphorylation networks implicated in multiple diseases

Protein kinases enable cellular information processing. Although numerous human phosphorylation sites and their dynamics have been characterized, the evolutionary history and physiological importance of many signaling events remain unknown. Using target phosphoproteomes determined with a similar experimental and computational pipeline, we investigated the conservation of human phosphorylation events in distantly related model organisms (fly, worm, and yeast). With a sequence-alignment approach, we identified 479 phosphorylation events in 344 human proteins that appear to be positionally conserved over approximately 600 million years of evolution and hence are likely to be involved in fundamental cellular processes. This sequence-alignment analysis suggested that many phosphorylation sites evolve rapidly and therefore do not display strong evolutionary conservation in terms of sequence position in distantly related organisms. Thus, we devised a network-alignment approach to reconstruct conserved kinase-substrate networks, which identified 778 phosphorylation events in 698 human proteins. Both methods identified proteins tightly regulated by phosphorylation as well as signal integration hubs, and both types of phosphoproteins were enriched in proteins encoded by disease-associated genes. We analyzed the cellular functions and structural relationships for these conserved signaling events, noting the incomplete nature of current phosphoproteomes. Assessing phosphorylation conservation at both site and network levels proved useful for exploring both fast-evolving and ancient signaling events. We reveal that multiple complex diseases seem to converge within the conserved networks, suggesting that disease development might rely on common molecular networks.

Evolution of phosphoregulation: comparison of phosphorylation patterns across yeast species

Analysis of the phosphoproteomes and the gene interaction networks of divergent yeast species defines the relative contribution of changes in protein phosphorylation pathways to the generation of phenotypic diversity.

The extent by which different cellular components generate phenotypic diversity is an ongoing debate in evolutionary biology that is yet to be addressed by quantitative comparative studies. We conducted an in vivo mass-spectrometry study of the phosphoproteomes of three yeast species (Saccharomyces cerevisiae, Candida albicans, and Schizosaccharomyces pombe) in order to quantify the evolutionary rate of change of phosphorylation. We estimate that kinase–substrate interactions change, at most, two orders of magnitude more slowly than transcription factor (TF)–promoter interactions. Our computational analysis linking kinases to putative substrates recapitulates known phosphoregulation events and provides putative evolutionary histories for the kinase regulation of protein complexes across 11 yeast species. To validate these trends, we used the E-MAP approach to analyze over 2,000 quantitative genetic interactions in S. cerevisiae and Sc. pombe, which demonstrated that protein kinases, and to a greater extent TFs, show lower than average conservation of genetic interactions. We propose therefore that protein kinases are an important source of phenotypic diversity.

Natural selection at a population level requires phenotypic diversity, which at the molecular level arises by mutation of the genome of each individual. What kinds of changes at the level of the DNA are most important for the generation of phenotypic differences remains a fundamental question in evolutionary biology. One well-studied source of phenotypic diversity is mutation in gene regulatory regions that results in changes in gene expression, but what proportion of phenotypic diversity is due to such mutations is not entirely clear. We investigated the relative contribution to phenotypic diversity of mutations in protein-coding regions compared to mutations in gene regulatory sequences. Given the important regulatory role played by phosphorylation across biological systems, we focused on mutations in protein-coding regions that alter protein–protein interactions involved in the binding of kinases to their substrate proteins. We studied the evolution of this “phosphoregulation” by analyzing the in vivo complement of phosphorylated proteins (the “phosphoproteome”) in three highly diverged yeast species—the budding yeast Saccharomyces cerevisiae, the pathogenic yeast Candida albicans, and the fission yeast Schizosaccharomyces pombe—and integrating those data with existing data on thousands of known genetic interactions from S. cerevisiae and Sc. pombe. We show that kinase–substrate interactions are altered at a rate that is at most two orders of magnitude slower than the alteration of transcription factor (TF)–promoter interactions, whereas TFs and kinases both show a faster than average rate of functional divergence estimated by the cross-species analysis of genetic interactions. Our data provide a quantitative estimate of the relative frequencies of different kinds of functionally relevant mutations and demonstrate that, like mutations in gene regulatory regions, mutations that result in changes in kinase–substrate interactions are an important source of phenotypic diversity.

Factors affecting the diversity of DNA replication licensing control in eukaryotes

Replication of eukaryotic genomes is limited to once per cell cycle, by a two-step mechanism. DNA replication origins are first "licensed" during G1 phase by loading of an inactive DNA helicase (Mcm2-7) into pre-replicative complexes (pre-RCs). Initiation then occurs during S phase, triggered by cyclin-dependent kinases (CDKs), which promote recruitment of proteins required for helicase activation and replisome assembly. CDKs and the anaphase promoting complex/cyclosome (APC/C) restrict licensing to G1 phase by directly and indirectly regulating pre-RC components, including ORC, Cdc6, Cdt1, and Mcm2-7. Despite the fundamental importance of licensing regulation, the mechanisms by which pre-RC components are regulated differ widely across Eukarya. Here we show that even within the genus Saccharomyces, Cdc6 is regulated differently in different species. We propose that two factors contribute to the rapid evolution of licensing regulation. The first is redundancy: eliminating any single pre-RC-regulatory mechanism has very little affect on viability. The second is interchangeability: we show that regulatory mechanisms can be swapped between pre-RC components without compromising the block to re-replication. These experiments provide a framework for understanding the diversity of licensing regulation in eukaryotes and provide new tools for manipulating the chromosome-replication cycle.

Genotype-phenotype mapping developmental biology confronts the toolkit paradox

The quest to understand the relationship between an organism's DNA sequence and three-dimensional form is an interdisciplinary task, integrating diverse fields of the life sciences. The relevance of the metaphor of a genotype-phenotype map is explored from a developmental perspective, in light of the recent concept of a "molecular toolkit" of protein-coding genes, and the widespread view that analyzing the logic and mechanics of gene regulation at multiple levels is key to explaining how morphology is genetically encoded. We discuss the challenges of decoding genomes despite variable genetic backgrounds, the dynamically changing physical and molecular contexts of the internal environment during development, and the impact of external forces on morphogenesis.