Phosphoproteome Analysis Reveals Regulatory Sites in Major Pathways of Cardiac Mitochondria

Mitochondrial functions are dynamically regulated in the heart. In particular, protein phosphorylation has been shown to be a key mechanism modulating mitochondrial function in diverse cardiovascular phenotypes. However, site-specific phosphorylation information remains scarce for this organ. Accordingly, we performed a comprehensive characterization of murine cardiac mitochondrial phosphoproteome in the context of mitochondrial functional pathways. A platform using the complementary fragmentation technologies of collision-induced dissociation (CID) and electron transfer dissociation (ETD) demonstrated successful identification of a total of 236 phosphorylation sites in the murine heart; 210 of these sites were novel. These 236 sites were mapped to 181 phosphoproteins and 203 phosphopeptides. Among those identified, 45 phosphorylation sites were captured only by CID, whereas 185 phosphorylation sites, including a novel modification on ubiquinol-cytochrome c reductase protein 1 (Ser-212), were identified only by ETD, underscoring the advantage of a combined CID and ETD approach. The biological significance of the cardiac mitochondrial phosphoproteome was evaluated. Our investigations illustrated key regulatory sites in murine cardiac mitochondrial pathways as targets of phosphorylation regulation, including components of the electron transport chain (ETC) complexes and enzymes involved in metabolic pathways (e.g. tricarboxylic acid cycle). Furthermore, calcium overload injured cardiac mitochondrial ETC function, whereas enhanced phosphorylation of ETC via application of phosphatase inhibitors restored calcium-attenuated ETC complex I and complex III activities, demonstrating positive regulation of ETC function by phosphorylation. Moreover, in silico analyses of the identified phosphopeptide motifs illuminated the molecular nature of participating kinases, which included several known mitochondrial kinases (e.g. pyruvate dehydrogenase kinase) as well as kinases whose mitochondrial location was not previously appreciated (e.g. Src). In conclusion, the phosphorylation events defined herein advance our understanding of cardiac mitochondrial biology, facilitating the integration of the still fragmentary knowledge about mitochondrial signaling networks, metabolic pathways, and intrinsic mechanisms of functional regulation in the heart.Mitochondria are the source of energy to sustain life. In addition to their evolutionary origin as an energy-producing organelle, their functionality has integrated into every aspect of life, including the cell cycle, ROS¹ production, apoptosis, and ion balance (1, 2). Our understanding of mitochondrial biology is still growing. Several systems biology approaches have been dedicated to exploring the molecular infrastructure and dynamics of the functional versatility associated with this organelle (3–5).To meet tissue-specific functional demands, mitochondria acquire heterogeneous properties in individual organs, a first statement of their plasticity in function and proteome composition (1, 6). The heterogeneity is evident even in an individual cardiomyocyte (7). A catalogue of the cardiac mitochondrial proteome is emerging via a joint effort (3–5). The dynamics of the mitochondrial proteome manifest at multiple levels, including post-translational modifications, such as phosphorylation. Our investigative goal is to decode this organellar proteome and its post-translational modification in a biological and functional context. In cardiomyocytes, mitochondria are also constantly exposed to fluctuation in energy demands and in ionic conditions. The capacity of mitochondria to cope with such a dynamic environment is essential for the functional role of mitochondria in normal and disease phenotypes (8–10). Unique protein features enabling the mitochondrial proteome to adapt to these biological changes can be interrogated by proteomics tools (10–12). Protein phosphorylation as a rapid and reversible chemical event is an integral component of these protein features (12–14).It has been estimated that one-third of cellular proteins exist in a phosphorylated state at least one time in their lifetime (15). However, only a handful of phosphorylation events have been identified to tune mitochondrial functionality (13, 14, 16) despite the fact that the first demonstration of phosphorylation was reported on a mitochondrial protein more than 5 decades ago (17). Kinases and phosphatases comprise nearly 3% of the human genome (18, 19). In mitochondria, ∼30 kinases and phosphatases have been identified thus far within the expected organellar proteome of a few thousand (3–5, 16). The number of identified mitochondrial phosphoproteins is far below one-third of its proteome size (20). Thus, it appears that the current pool of reported phosphoproteins represents only a small fraction of the anticipated mitochondrial phosphoproteome. The seminal studies from several groups (12–14, 16) demonstrated the prevalence as well as the dynamic nature of phosphorylation in cardiac mitochondria, suggesting that obtaining a comprehensive map of the mitochondrial phosphoproteome is feasible.In this study, we took a systematic approach to tackle the phosphorylation of murine cardiac mitochondrial pathways. We applied the unique strengths of both electron transfer dissociation (ETD) and collision-induced dissociation (CID) LC-MS/MS to screen phosphorylation events in a site-specific fashion. A total of 236 phosphorylation sites in 203 unique phosphopeptides were identified and mapped to 181 phosphoproteins. Novel phosphorylation modifications were discovered in diverse pathways of mitochondrial biology, including ion balance, proteolysis, and apoptosis. Consistent with the role of mitochondria as the major source of energy production under delicate control, metabolic pathways claimed one-third of phosphorylation sites captured in this analysis. To study molecular players steering mitochondrial phosphorylation, we probed the effects of calcium loading on phosphorylation. In addition, a number of kinases with previously unappreciated mitochondrial residence are suggested as potential players modulating mitochondrial pathways. Taken together, the cohort of novel phosphorylation events discovered in this study constitutes an essential step toward the full delineation of the cardiac mitochondrial phosphoproteome.     

Quantitative Site-specific Phosphorylation Dynamics of Human Protein Kinases during Mitotic Progression

Reversible protein phosphorylation is a key regulatory mechanism of mitotic progression. Importantly, protein kinases themselves are also regulated by phosphorylation-dephosphorylation processes; hence, phosphorylation dynamics of kinases hold a wealth of information about phosphorylation networks. Here, we investigated the site-specific phosphorylation dynamics of human kinases during mitosis using synchronization of HeLa suspension cells, kinase enrichment, and high resolution mass spectrometry. In biological triplicate analyses, we identified 206 protein kinases and more than 900 protein kinase phosphorylation sites, including 61 phosphorylation sites on activation segments, and quantified their relative abundances across three specific mitotic stages. Around 25% of the kinase phosphorylation site ratios were found to be changed by at least 50% during mitotic progression. Further network analysis of jointly regulated kinase groups suggested that Cyclin-dependent kinase- and mitogen-activated kinase-centered interaction networks are coordinately down- and up-regulated in late mitosis, respectively. Importantly, our data cover most of the already known mitotic kinases and, moreover, identify attractive candidates for future studies of phosphorylation-based mitotic signaling. Thus, the results of this study provide a valuable resource for cell biologists and provide insight into the system properties of the mitotic phosphokinome.Reversible phosphorylation is a ubiquitous posttranslational protein modification that is involved in the regulation of almost all biological processes (1–3). In human, 518 protein kinases have been identified in the genome that phosphorylate the majority of cellular proteins and increase the diversity of the proteome by severalfold (4). Addition of a phosphate group to a protein can alter its structural, catalytic, and functional properties; hence, kinases require tight regulation to avoid unspecific phosphorylation, which can be deleterious to cells (5–7). As a result, cells use a variety of mechanisms to ensure proper regulation of kinase activities (8). Importantly, most kinases are also in turn regulated through autophosphorylation and phosphorylation by other kinases, thus generating complex phosphorylation networks. In particular, phosphorylation on activation segments is a common mechanism to modulate kinase activities (9–11), but additional phosphorylation sites are also frequently required for fine tuning of kinase localizations and functions (12). Some kinases contain phosphopeptide binding domains that recognize prephosphorylated sites on other kinases, resulting in processive phosphorylation and/or targeting of kinases to distinct cellular locations (13–16). Because such priming phosphorylation events depend on the activities of the priming kinases, these motifs act as conditional docking sites and restrict the interaction with docking kinases to a particular point in time and physiological state. In addition, phosphorylation sites may act through combinatorial mechanisms or through cross-talk with other posttranslational modifications (PTMs)¹ (17, 18), thus further increasing the complexity of kinase regulatory networks.Regulation of kinases is of particular interest in mitosis as most of the mitotic events are regulated by reversible protein phosphorylation (19). During mitosis, error-free segregation of sister chromatids into the two daughter cells is essential to ensure genomic stability. Physically, this process is carried out by the mitotic spindle, a highly dynamic microtubule-based structure. After entry into mitosis, the major microtubule-organizing centers in animal cells, the centrosomes, start to increase microtubule nucleation and move to opposite poles of the cell. Throughout prometaphase, microtubules emanating from centrosomes are captured by kinetochores, protein complexes assembled on centromeric chromosomal DNA. This eventually leads to the alignment of all chromosomes in a metaphase plate. Because proper bipolar attachment of chromosomes to spindle microtubules is essential for the correct segregation of chromosomes, this critical step is monitored by a signaling pathway known as the spindle assembly checkpoint (SAC) (20). This checkpoint is silenced only after all chromosomes have attached to the spindle in a bioriented fashion, resulting in the synchronous segregation of sister chromatids during anaphase. Simultaneously, a so-called central spindle is formed between the separating chromatids, and the formation of a contractile ring initiates cytokinesis. Finally, in telophase, the chromosomes decondense and reassemble into nuclei, whereas remnants of the central spindle form the midbody, marking the site of abscission. Cyclin-dependent kinase 1 (Cdk1), an evolutionarily conserved master mitotic kinase, is activated prior to mitosis and initiates most of the mitotic events. Cdk1 works in close association with other essential mitotic kinases such as Plk1, Aurora A, and Aurora B for the regulation of mitotic progression (19, 21–24). Plk1 and Aurora kinases dynamically localize to different subcellular locations to perform multiple functions during mitosis and are phosphorylated at several conserved sites. Although little is known about the precise roles of these phosphorylation sites, emerging data indicate that they are involved in regulating localization-specific functions (25, 26). Furthermore, the kinases Bub1, BubR1, and TTK (Mps1) and kinases of the Nek family play important roles in maintaining the fidelity and robustness of mitosis (19). Recently, a genome-wide RNA-mediated interference screen identified M phase phenotypes for many kinases that have not previously been implicated in cell cycle functions, indicating that additional kinases have important mitotic functions (27).Although protein phosphorylation plays a pivotal role in the regulation of cellular networks, many phosphorylation events remain undiscovered mainly because of technical limitations (28). The advent of mass spectrometry-based proteomics along with developments in phosphopeptide enrichment methods has enabled large scale global phosphoproteomics studies (29, 30). However, the number of phosphorylation sites identified on kinases is limited compared with other proteins because of their frequently low expression levels. To overcome this problem, small inhibitor-based kinase enrichment strategies were developed, resulting in the identification of more than 200 kinases from HeLa cell lysates (31, 32). This method was also used recently to compare the phosphokinomes during S phase and M phase of the cell cycle, resulting in the identification of several hundreds of M phase-specific kinase phosphorylation sites (31). In the present study, we address the dynamics of the phosphokinome during mitotic progression using large scale cell synchronization at three distinct mitotic stages, small inhibitor-based kinase enrichment, and stable isotope labeling by amino acids in cell culture (SILAC)-based quantitative mass spectrometry. Thus, we determined the mitotic phosphorylation dynamics of more than 900 kinase phosphorylation sites and identified distinctly regulated kinase interaction networks. Our results provide a valuable resource for the dynamics of the kinome during mitotic progression and give insight into the system properties of kinase interaction networks.     

Phosphoproteome Analysis of Functional Mitochondria Isolated from Resting Human Muscle Reveals Extensive Phosphorylation of Inner Membrane Protein Complexes and Enzymes

Mitochondria play a central role in energy metabolism and cellular survival, and consequently mitochondrial dysfunction is associated with a number of human pathologies. Reversible protein phosphorylation emerges as a central mechanism in the regulation of several mitochondrial processes. In skeletal muscle, mitochondrial dysfunction is linked to insulin resistance in humans with obesity and type 2 diabetes. We performed a phosphoproteomics study of functional mitochondria isolated from human muscle biopsies with the aim to obtain a comprehensive overview of mitochondrial phosphoproteins. Combining an efficient mitochondrial isolation protocol with several different phosphopeptide enrichment techniques and LC-MS/MS, we identified 155 distinct phosphorylation sites in 77 mitochondrial phosphoproteins, including 116 phosphoserine, 23 phosphothreonine, and 16 phosphotyrosine residues. The relatively high number of phosphotyrosine residues suggests an important role for tyrosine phosphorylation in mitochondrial signaling. Many of the mitochondrial phosphoproteins are involved in oxidative phosphorylation, tricarboxylic acid cycle, and lipid metabolism, i.e. processes proposed to be involved in insulin resistance. We also assigned phosphorylation sites in mitochondrial proteins involved in amino acid degradation, importers and transporters, calcium homeostasis, and apoptosis. Bioinformatics analysis of kinase motifs revealed that many of these mitochondrial phosphoproteins are substrates for protein kinase A, protein kinase C, casein kinase II, and DNA-dependent protein kinase. Our results demonstrate the feasibility of performing phosphoproteome analysis of organelles isolated from human tissue and provide novel targets for functional studies of reversible phosphorylation in mitochondria. Future comparative phosphoproteome analysis of mitochondria from healthy and diseased individuals will provide insights into the role of abnormal phosphorylation in pathologies, such as type 2 diabetes.Mitochondria are the primary energy-generating systems in eukaryotes. They play a crucial role in oxidative metabolism, including carbohydrate metabolism, fatty acid oxidation, and urea cycle, as well as in calcium signaling and apoptosis (1, 2). Mitochondrial dysfunction is centrally involved in a number of human pathologies, such as type 2 diabetes, Parkinson disease, and cancer (3). The most prevalent form of cellular protein post-translational modifications (PTMs),1 reversible phosphorylation (4–6), is emerging as a central mechanism in the regulation of mitochondrial functions (7, 8). The steadily increasing numbers of reported mitochondrial kinases, phosphatases, and phosphoproteins imply an important role of protein phosphorylation in different mitochondrial processes (9–11).Mass spectrometry (MS)-based proteome analysis is a powerful tool for global profiling of proteins and their PTMs, including protein phosphorylation (12, 13). A variety of proteomics techniques have been developed for specific enrichment of phosphorylated proteins and peptides and for phosphopeptide-specific data acquisition techniques at the MS level (14). Enrichment methods based on affinity chromatography, such as titanium dioxide (TiO₂) (15–17), zwitterionic hydrophilic interaction chromatography (ZIC-HILIC) (18), immobilized metal affinity chromatography (IMAC) (19, 20), and ion exchange chromatography (strong anion exchange and strong cation exchange) (21, 22), have shown high efficiencies for enrichment of phosphopeptides (14). Recently, we demonstrated that calcium phosphate precipitation (CPP) is highly effective for enriching phosphopeptides (23). It is now generally accepted that no single method is comprehensive, but combinations of different enrichment methods produce distinct overlapping phosphopeptide data sets to enhance the overall results in phosphoproteome analysis (24, 25). Phosphopeptide sequencing by mass spectrometry has seen tremendous advances during the last decade (26). For example, MS/MS product ion scanning, multistage activation, and precursor ion scanning are effective methods for identifying serine (Ser), threonine (Thr), and tyrosine (Tyr) phosphorylated peptides (14, 26).A “complete” mammalian mitochondrial proteome was reported by Mootha and co-workers (27) and included 1098 proteins. The mitochondrial phosphoproteome has been characterized in a series of studies, including yeast, mouse and rat liver, porcine heart, and plants (19, 28–31). To date, the largest data set by Deng et al. (30) identified 228 different phosphoproteins and 447 phosphorylation sites in rat liver mitochondria. However, the in vivo phosphoproteome of human mitochondria has not been determined. A comprehensive mitochondrial phosphoproteome is warranted for further elucidation of the largely unknown mechanisms by which protein phosphorylation modulates diverse mitochondrial functions.The percutaneous muscle biopsy technique is an important tool in the diagnosis and management of human muscle disorders and has been widely used to investigate metabolism and various cellular and molecular processes in normal and abnormal human muscle, in particular the molecular mechanism underlying insulin resistance in obesity and type 2 diabetes (32). Skeletal muscle is rich in mitochondria and hence a good source for a comprehensive proteomics and functional analysis of mitochondria (32, 33).The major aim of the present study was to obtain a comprehensive overview of site-specific phosphorylation of mitochondrial proteins in functionally intact mitochondria isolated from human skeletal muscle. Combining an efficient protocol for isolation of skeletal muscle mitochondria with several different state-of-the-art phosphopeptide enrichment methods and high performance LC-MS/MS, we identified 155 distinct phosphorylation sites in 77 mitochondrial phosphoproteins, many of which have not been reported before. We characterized this mitochondrial phosphoproteome by using bioinformatics tools to classify functional groups and functions, including kinase substrate motifs.     

Evolutionary Constraints of Phosphorylation in Eukaryotes,Prokaryotes, and Mitochondria

High accuracy mass spectrometry has proven to be a powerful technology for the large scale identification of serine/threonine/tyrosine phosphorylation in the living cell. However, despite many described phosphoproteomes, there has been no comparative study of the extent of phosphorylation and its evolutionary conservation in all domains of life. Here we analyze the results of phosphoproteomics studies performed with the same technology in a diverse set of organisms. For the most ancient organisms, the prokaryotes, only a few hundred proteins have been found to be phosphorylated. Applying the same technology to eukaryotic species resulted in the detection of thousands of phosphorylation events. Evolutionary analysis shows that prokaryotic phosphoproteins are preferentially conserved in all living organisms, whereas-site specific phosphorylation is not. Eukaryotic phosphosites are generally more conserved than their non-phosphorylated counterparts (with similar structural constraints) throughout the eukaryotic domain. Yeast and Caenorhabditis elegans are two exceptions, indicating that the majority of phosphorylation events evolved after the divergence of higher eukaryotes from yeast and reflecting the unusually large number of nematode-specific kinases. Mitochondria present an interesting intermediate link between the prokaryotic and eukaryotic domains. Applying the same technology to this organelle yielded 174 phosphorylation sites mapped to 74 proteins. Thus, the mitochondrial phosphoproteome is similarly sparse as the prokaryotic phosphoproteomes. As expected from the endosymbiotic theory, phosphorylated as well as non-phosphorylated mitochondrial proteins are significantly conserved in prokaryotes. However, mitochondrial phosphorylation sites are not conserved throughout prokaryotes, consistent with the notion that serine/threonine phosphorylation in prokaryotes occurred relatively recently in evolution. Thus, the phosphoproteome reflects major events in the evolution of life.Reversible protein phosphorylation on serines, threonines, and tyrosines plays a crucial role in regulating processes in all living organisms ranging from prokaryotes to eukaryotes (1). Traditionally, phosphorylation has been detected in single, purified proteins using in vitro assays. Recent advances in mass spectrometry (MS)-based proteomics now allow the identification of in vivo phosphorylation sites with high accuracy (2–7). On-line databases such as PhosphoSite (8), Phospho.ELM (9), and PHOSIDA1 (10) have collected and organized thousands of identified phosphosites. These databases as well as dedicated analysis environments such as NetworKIN (11, 12) offer and use contextual information including structural features, potential kinases, and conservation. They constitute resources that should allow the derivation of general patterns for phosphorylation events. Specifically, the recent availability of data for archaeal, prokaryotic, and diverse eukaryotic phosphoproteomes in these databases should enable investigation of the evolutionary history of this post-translational modification.Prokaryotes have two separate classes of phosphorylation events. Apart from the canonical histidine/aspartate phosphorylation, which has been studied for decades, serine/threonine/tyrosine phosphorylation is also present and has recently become amenable to analysis by MS (13). Bacterial phosphoproteins are involved in protein synthesis, carbohydrate metabolism, and the phosphoenolpyruvate-dependent phosphotransferase system. Recent phosphoproteomics studies of Bacillus subtilis, Escherichia coli, and Lactococcus lactis described around 100 phosphorylation sites on serine, threonine, and tyrosine in each of these species (13–15). Bacterial phosphorylation sites can change in response to environmental conditions (16).Interestingly, even archaea have serine/threonine and tyrosine phosphorylation. A recent study of Halobacterium salinarum described 75 serine/threonine/tyrosine phosphorylation sites on 62 proteins involved in a wide range of cellular processes including a variety of metabolic pathways (17).Although only a few hundred phosphorylation events have been found in prokaryotic species, similar experimental conditions and effort have yielded the detection of thousands of phosphorylation events in eukaryotes ranging from yeast to human (7, 18–21). Before the advent of large scale phosphoproteomics, serine/threonine/tyrosine phosphorylation has been estimated to affect one-third of all proteins (22). Recent large scale phosphoproteomics studies now suggest that more than half of all eukaryotic proteins are phosphorylated (23).A key event in evolution was the endosymbiosis of prokaryotes that enabled the development of a much more complex type of life, the eukaryotic cell. Analyses of mitochondrial genes suggest that the α-proteobacterium Rickettsia prowazekii is the endosymbiotic precursor leading to modern mitochondria (24). Almost all of the mitochondrial genes have migrated to the nuclear genome during subsequent evolution, and it is predicted that 10–15% of eukaryotic nuclear genes of organisms encode mitochondrial proteins (25).Thus, mitochondria with their unique evolutionary position between prokaryotes and eukaryotes form an interesting link for the evolutionary analysis of phosphorylation. Several studies investigated the mitochondrial phosphoproteome in different organisms using gel electrophoresis or specific enrichment methods coupled with mass spectrometry (26–28). Those studies established potential mitochondrial phosphoproteins. Three large scale studies based on affinity enrichment of phosphopeptides and mass spectrometry obtained direct experimental evidence of phosphorylation sites in mitochondria. Lee et al. (29) used a combination of different peptide enrichment strategies and found 80 phosphorylation sites of 48 different proteins from mouse liver. Very recently, a study by Deng et al. (30) characterized the murine cardiac mitochondrial mouse phosphoproteome, covering 236 phosphosites on 181 proteins. Investigating yeast, Reinders et al. (31) assigned 84 phosphorylation sites in 62 proteins.To enable comparative analysis of phosphoproteomes between all domains of life and mitochondria, here we experimentally determined a high accuracy mitochondrial mouse phosphoproteome based on technology conditions similar to those applied to the identification of prokaryotic and eukaryotic phosphoproteomes. We then performed a detailed evolutionary study of the conservation of the identified phosphoproteins and phosphorylation sites in prokaryotes and in eukaryotes. This allowed an initial comparison of the phosphoproteomes of prokaryotes, mitochondria, and eukaryotes.     

A Proteome-wide Domain-centric Perspective on Protein Phosphorylation

Phosphorylation is a widespread post-translational modification that modulates the function of a large number of proteins. Here we show that a significant proportion of all the domains in the human proteome is significantly enriched or depleted in phosphorylation events. A substantial improvement in phosphosites prediction is achieved by leveraging this observation, which has not been tapped by existing methods. Phosphorylation sites are often not shared between multiple occurrences of the same domain in the proteome, even when the phosphoacceptor residue is conserved. This is partly because of different functional constraints acting on the same domain in different protein contexts. Moreover, by augmenting domain alignments with structural information, we were able to provide direct evidence that phosphosites in protein-protein interfaces need not be positionally conserved, likely because they can modulate interactions simply by sitting in the same general surface area.Phosphorylation, the most widespread protein post-translational modification, is an important regulator of protein function. The addition of phosphate groups on serine, threonine, and tyrosine residues can modulate the activity of the target protein by inducing complex conformational changes, by modifying protein electrostatics, and by regulating domain-peptide interactions, as in 14-3-3 or SH2 domains, that specifically recognize phosphorylated residues. The standard experimental technique for the high-throughput identification of phosphorylation sites is mass spectrometry (1).Phosphorylation is catalyzed by protein kinases, a family that in humans comprises ∼540 members (2, 3). It is well understood that these enzymes recognize specific sequence motifs in their substrates (4, 5). Accordingly the sequence around the phosphorylation site is undisputedly the most important feature for phosphosite prediction (6, 7). However the “context,” in a broad sense, where these motifs occur is also important as sequence alone is not enough to achieve the observed specificity of phosphorylation. Therefore, several studies have characterized multiple aspects of phosphosites such as their preference for loops and disordered regions (reviewed in (8)), or the tendency of phosphoserines and phosphothreonines to occur in clusters (9), and these features have been used to improve the performance of phosphosite predictors (6, 7, 10–12). Moreover placing kinases and substrates in the context of protein interaction networks has been shown to improve the prediction of phosphorylation by specific kinases (13).Perhaps one of the most puzzling observations when looking at the phosphoproteome as a whole, is the fact that a large proportion of phosphorylation sites is poorly conserved. This has led to various hypotheses. First some sites may represent nonfunctional, possibly low-stoichiometry, phosphorylation events that are picked up because of the sensitivity of mass-spectrometry (14, 15). Indeed functionally characterized sites and those matching known kinase motifs are more conserved on average (15–17). However, although in biology function often equates with conservation, there could be genuinely functional fast-evolving phosphosites, that are responsible for species-specific differences in signaling and regulation. Moreover in some cases, especially in the regulation of protein-protein interactions, the exact position of the phosphosites may be unimportant (18, 19).Here we explore the issues of “context” and “conservation” of phosphorylation sites from the perspective of protein domains. To this end, we assembled a comprehensive database of phosphosites from publicly available sources and studied their proteome distribution with respect to the location and identity of protein domains. We focus on the human phosphoproteome because it has been very well characterized in a multitude of low- and high-throughput experiments, thus providing the opportunity for a comprehensive, proteome-wide, study. In particular, the issues we want to address are the following:

1. Are specific domain types preferentially phosphorylated? Or conversely are some domains specifically depleted of phosphorylation sites?
2. Can the domain context be used to improve the prediction of phosphorylation sites?
3. What is the conservation pattern of phosphosites when looking at multiple instances of the same domain in the proteome?

     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Phosphorylation Dependence of Hsp27 Multimeric Size and Molecular Chaperone Function

The molecular chaperone Hsp27 exists as a distribution of large oligomers that are disassembled by phosphorylation at Ser-15, -78, and -82. It is controversial whether the unphosphorylated Hsp27 or the widely used triple Ser-to-Asp phospho-mimic mutant is the more active molecular chaperone in vitro. This question was investigated here by correlating chaperone activity, as measured by the aggregation of reduced insulin or α-lactalbumin, with Hsp27 self-association as monitored by analytical ultracentrifugation. Furthermore, because the phospho-mimic is generally assumed to reproduce the phosphorylated molecule, the size and chaperone activity of phosphorylated Hsp27 were compared with that of the phospho-mimic. Hsp27 was triply phosphorylated by MAPKAP-2 kinase, and phosphorylation was tracked by urea-PAGE. An increasing degree of suppression of insulin or α-lactalbumin aggregation correlated with a decreasing Hsp27 self-association, which was the least for phosphorylated Hsp27 followed by the mimic followed by the unphosphorylated protein. It was also found that Hsp27 added to pre-aggregated insulin did not reverse aggregation but did inhibit these aggregates from assembling into even larger aggregates. This chaperone activity appears to be independent of Hsp27 phosphorylation. In conclusion, the most active chaperone of insulin and α-lactalbumin was the Hsp27 (elongated) dimer, the smallest Hsp27 subunit observed under physiological conditions. Next, the Hsp27 phospho-mimic is only a partial mimic of phosphorylated Hsp27, both in self-association and in chaperone function. Finally, the efficient inhibition of insulin aggregation by Hsp27 dimer led to the proposal of two models for this chaperone activity.Oligomeric heat shock protein 27 (Hsp27)² is a ubiquitous mammalian protein with a variety of functions in health and disease (1–8). These functions include ATP-independent chaperone activity in response to environmental stress, e.g. heat shock and oxidative stress, control of apoptosis, and regulation of actin cytoskeleton dynamics. Hsp27 is a member of the α-crystallin small heat shock protein family of which αB-crystallin is the archetype. These proteins are characterized by an α-crystallin domain of 80–90 residues consisting of roughly eight β-strands that form an intermolecular β-sheet interaction interface within a dimer, the basic building subunit of the oligomer (2, 4, 9–11).Hsp27 is in equilibrium between high molecular weight oligomers and much lower molecular weight multimers. It has been reported that unphosphorylated Hsp27 includes predominantly a distribution of high molecular species ranging in size from 12-mer to 35-mer (12–19). Phosphorylation of Hsp27 at serines 15, 78, and 82 by the p38-activated MAPKAP-2 kinase (20–22) or the use of the triple Ser-to-Asp phospho-mimic results in a major shift in the equilibrium toward much smaller multimers (23) and in an alteration of its function (1, 3, 6, 7, 24, 25). The size distribution of the smaller species has been reported to be between monomer and tetramer (12–16, 18, 19).Small heat shock proteins, including Hsp27, behave as ATP-independent molecular chaperones during cellular heat shock. They bind partially unfolded proteins and prevent their aggregation until the proteins can be refolded by larger ATP-dependent chaperones or are digested (7, 8, 26). This function includes the up-regulation and/or phosphorylation of Hsp27.It is not entirely clear what the role of Hsp27 size and phosphorylation state plays in its heat shock function because there are conflicting results in the literature. Some in vitro studies concluded that the unphosphorylated oligomeric Hsp27 (or the murine isoform Hsp25) protects proteins against aggregation better than does the phosphorylation mimic (13, 19, 27), whereas others found no difference (16, 28, 29), and still other studies found that the mimic protects better than does the unphosphorylated wild type (27, 30, 31). In-cell studies found that phosphorylation of Hsp27 was essential for thermo-protection of actin filaments (32), and the Hsp27 phosphorylation mimic decreased inclusion body formation better than did unphosphorylated Hsp27 (33). This study was undertaken to investigate the molecular chaperone function of Hsp27 by correlating chaperone activity with Hsp27 size and by comparing fully phosphorylated Hsp27 with its phospho-mimic.     

Semaphorin3A Signaling Mediated by Fyn-dependent Tyrosine Phosphorylation of Collapsin Response Mediator Protein 2 at Tyrosine 32

Collapsin response mediator protein 2 (CRMP2) is an intracellular protein that mediates signaling of Semaphorin3A (Sema3A), a repulsive axon guidance molecule. Fyn, a Src-type tyrosine kinase, is involved in the Sema3A signaling. However, the relationship between CRMP2 and Fyn in this signaling pathway is still unknown. In our research, we demonstrated that Fyn phosphorylated CRMP2 at Tyr³² residues in HEK293T cells. Immunohistochemical analysis using a phospho-specific antibody at Tyr³² of CRMP showed that Tyr³²-phosphorylated CRMP was abundant in the nervous system, including dorsal root ganglion neurons, the molecular and Purkinje cell layer of adult cerebellum, and hippocampal fimbria. Overexpression of a nonphosphorylated mutant (Tyr³² to Phe³²) of CRMP2 in dorsal root ganglion neurons interfered with Sema3A-induced growth cone collapse response. These results suggest that Fyn-dependent phosphorylation of CRMP2 at Tyr³² is involved in Sema3A signaling.Collapsin response mediator proteins (CRMPs)⁴ have been identified as intracellular proteins that mediate Semaphorin3A (Sema3A) signaling in the nervous system (1). CRMP2 is one of the five members of the CRMP family. CRMPs also mediate signal transduction of NT3, Ephrin, and Reelin (2–4). CRMPs interact with several intracellular molecules, including tubulin, Numb, kinesin1, and Sra1 (5–8). CRMPs are involved in axon guidance, axonal elongation, cell migration, synapse maturation, and the generation of neuronal polarity (1, 2, 4, 5).CRMP family proteins are known to be the major phosphoproteins in the developing brain (1, 9). CRMP2 is phosphorylated by several Ser/Thr kinases, such as Rho kinase, cyclin-dependent kinase 5 (Cdk5), and glycogen synthase kinase 3β (GSK3β) (2, 10–13). The phosphorylation sites of CRMP2 by these kinases are clustered in the C terminus and have already been identified. Rho kinase phosphorylates CRMP2 at Thr⁵⁵⁵ (10). Cdk5 phosphorylates CRMP2 at Ser⁵²², and this phosphorylation is essential for sequential phosphorylations by GSK3β at Ser⁵¹⁸, Thr⁵¹⁴, and Thr⁵⁰⁹ (2, 11–13). These phosphorylations disrupt the interaction of CRMP2 with tubulin or Numb (2, 3, 13). The sequential phosphorylation of CRMP2 by Cdk5 and GSK3β is an essential step in Sema3A signaling (11, 13). Furthermore, the neurofibrillary tangles in the brains of people with Alzheimer disease contain hyperphosphorylated CRMP2 at Thr⁵⁰⁹, Ser⁵¹⁸, and Ser⁵²² (14, 15).CRMPs are also substrates of several tyrosine kinases. The phosphorylation of CRMP2 by Fes/Fps and Fer has been shown to be involved in Sema3A signaling (16, 17). Phosphorylation of CRMP2 at Tyr⁴⁷⁹ by a Src family tyrosine kinase Yes regulates CXCL12-induced T lymphocyte migration (18). We reported previously that Fyn is involved in Sema3A signaling (19). Fyn associates with PlexinA2, one of the components of the Sema3A receptor complex. Fyn also activates Cdk5 through the phosphorylation at Tyr¹⁵ of Cdk5 (19). In dorsal root ganglion (DRG) neurons from fyn-deficient mice, Sema3A-induced growth cone collapse response is attenuated compared with control mice (19). Furthermore, we recently found that Fyn phosphorylates CRMP1 and that this phosphorylation is involved in Reelin signaling (4). Although it has been shown that CRMP2 is involved in Sema3A signaling (1, 11, 13), the relationship between Fyn and CRMP2 in Sema3A signaling and the tyrosine phosphorylation site(s) of CRMPs remain unknown.Here, we show that Fyn phosphorylates CRMP2 at Tyr³². Using a phospho-specific antibody against Tyr³², we determined that the residue is phosphorylated in vivo. A nonphosphorylated mutant CRMP2Y32F inhibits Sema3A-induced growth cone collapse. These results indicate that tyrosine phosphorylation by Fyn at Tyr³² is involved in Sema3A signaling.     

Systematic Analysis of the Phosphoproteome and Kinase-substrate Networks in the Mouse Testis

Spermatogenesis is a complex process closely associated with the phosphorylation-orchestrated cell cycle. Elucidating the phosphorylation-based regulations should advance our understanding of the underlying molecular mechanisms. Here we present an integrative study of phosphorylation events in the testis. Large-scale phosphoproteome profiling in the adult mouse testis identified 17,829 phosphorylation sites in 3955 phosphoproteins. Although only approximately half of the phosphorylation sites enriched by IMAC were also captured by TiO₂, both the phosphoprotein data sets identified by the two methods significantly enriched the functional annotation of spermatogenesis. Thus, the phosphoproteome profiled in this study is a highly useful snapshot of the phosphorylation events in spermatogenesis. To further understand phosphoregulation in the testis, the site-specific kinase-substrate relations were computationally predicted for reconstructing kinase-substrate phosphorylation networks. A core sub-kinase-substrate phosphorylation networks among the spermatogenesis-related proteins was retrieved and analyzed to explore the phosphoregulation during spermatogenesis. Moreover, network-based analyses demonstrated that a number of protein kinases such as MAPKs, CDK2, and CDC2 with statistically more site-specific kinase-substrate relations might have significantly higher activities and play an essential role in spermatogenesis, and the predictions were consistent with previous studies on the regulatory roles of these kinases. In particular, the analyses proposed that the activities of POLO-like kinases (PLKs) might be dramatically higher, while the prediction was experimentally validated by detecting and comparing the phosphorylation levels of pT210, an indicator of PLK1 activation, in testis and other tissues. Further experiments showed that the inhibition of POLO-like kinases decreases cell proliferation by inducing G2/M cell cycle arrest. Taken together, this systematic study provides a global landscape of phosphoregulation in the testis, and should prove to be of value in future studies of spermatogenesis.Spermatogenesis is a complex sperm-generating process involving the mitosis of spermatogonia, meiosis of spermatocytes, and spermiogenesis of spermatids. Sperms are produced in the male testis at the speed of ∼1000 sperm per heart beat (1), which indicate that spermatogenesis is an extremely dynamic process in the testis. The protein expression levels during spermatogenesis have been well studied by high-throughput proteomic studies, and over 7000 proteins have been identified in the mammalian testis (2–4). However, the dynamic regulatory events that orchestrate this complex process have yet to be elucidated. Because phosphorylation, an important and ubiquitous post-translational modification (PTM)¹, is one of the most critical regulatory mechanisms of the cell cycle (5), which is particularly active during spermatogenesis, a number of pioneering studies have contributed to our understanding of phosphoregulation in spermatogenesis. For example, mitogen-activated protein kinases (MAPKs) such as ERK1/2, were found to play an important role in ectoplasmic specialization dynamics during spermatogenesis (6). As important regulators of the cell cycle (7, 8), the POLO-like kinases (PLKs) especially PLK1, were found to be required at multiple stages of spermatogenesis (9–12). Thus, a systematic analysis of phosphorylation in the testis is of great importance for advancing the current understanding of the molecular mechanisms of spermatogenesis.In order to elucidate the phosphorylation-mediated regulation of spermatogenesis, the characterization of the testicular phosphoproteome could serve as a straightforward start. Recently, rapid progress in mass spectrometry based proteomic technologies has greatly advanced to a state-of-art stage at which thousands of PTM sites can be identified in a single run (13). Although a large proportion of these studies were carried out in cell lines, only a handful of studies have contributed to the identification of phosphoproteome in testis and sperm (14–19). For example, Huttlin et al. identified ∼36,000 phosphorylation sites in 6296 proteins from nine tissues, including the 3-week testis of immature mice (16). Because neither elongated spermatids nor sperms exist in such immature testes (16), it might be impossible to identify of phosphorylation sites across all stages of spermatogenesis from the samples. Moreover, a recent study characterized the testicular phosphoproteome in the perfluorododecanoic acid (PFDoA)-exposed rats, and demonstrated the importance of MAPK pathway and CDC2 protein phosphorylation in the toxicity of PFDoA (19). Taken together, despite the fact that a number of studies have been carried out (18), our understanding of the testicular phosphoproteome is still limited, and more effort needs to be expended on this area.In coordination with the exploration of the phosphoproteome, the technology for analyzing kinase-substrate relations (KSRs) has also greatly advanced. In addition to conventional experimental approaches, a number of computational studies have been carried out (8, 20–23), whereas network approaches have attracted growing attention (20, 22). In 2007, Linding et al. first constructed a human kinase-substrate phosphorylation network (KSPN) through the prediction of site-specific kinase-substrate relations (ssKSRs) with a novel algorithm NetworKIN (24, 25). Combined with sequence-based predictions using the Group-based Prediction System (GPS) algorithm and protein–protein interactions (PPIs), we also developed iGPS (in vivo GPS) software to reconstruct ssKSR-based KSPNs in eukaryotes, and achieved a superior performance compared with NetworKIN (26). With these computational tools, network-based analyses can be performed for mining phospho-signatures from the phosphoproteomic data. For example, based on a hypothesis that a kinase with higher activity will phosphorylate more sites, we designed a novel computational method of kinase activity analysis (KAA) to statistically identify kinases with significantly more or less phosphorylation sites (20, 26). Using the human whole phosphoproteome as a background, we totally detected 60 kinases with higher activities (i.e. with more sites) from a human liver phosphoproteome (26). Our hypothesis and methodology was successfully supported by following studies, which used phospho-specific antibodies to validate the modification levels of the activity-associated autophosphorylation sites in the predicted kinases (27, 28). In particular, the two studies demonstrated that kinases predicted with significantly higher activities can act as important regulators in distinct biological processes by regulating the KSPNs (27, 28). However, such an analysis of potentially differential kinase activities in spermatogenesis still remains to be performed.In this study, we systematically profiled the phosphoproteome in the adult mouse testes. Using phosphopeptide enrichment methods, including immobilized metal affinity chromatography (IMAC) and Titanium dioxide (TiO₂), high-throughput mass spectrometry identified 17,829 phosphorylation sites in 3955 proteins in the adult mouse testis. Although only approximately half of the phosphorylation sites enriched by IMAC were also enriched by TiO₂, statistical analyses of the gene ontology (GO) terms consistently found the GO term “spermatogenesis” to be significantly over-represented. Thus, as the first comprehensive phosphoproteome in mature testis, these results provide an in-depth picture of phosphorylation in spermatogenesis. To further investigate phosphoregulation, the ssKSRs were predicted and employed to re-construct the KSPNs in the testis. Based on the working concept that kinases with a higher level of activity phosphorylate more sites (26), the predicted ssKSRs were used to predict the kinase activity profiles. Although the overlap of different phosphoproteome data sets is limited, the kinase activity profiles indicate a pattern of consistently high activity for a number of kinases, including the MAPKs, CDKs, and especially the POLO-like kinases (PLKs). Through Western blot detection of the phosphorylation levels of T210, which is positively correlated with PLK1 activation (29–32), it was observed that PLK1 was highly activated in testis. The PLKs inhibition assay results showed that PLKs activities are critical for cell proliferation in the spermatocyte GC2 cell line, whereas PLKs inhibition generated G2/M arrest. Taken together, this study of the testicular phosphoproteome provides a systematic understanding of the phosphorylation that occurs during spermatogenesis, with the results able to serve as a resource for future investigation.     

Use of Quantitative Mass Spectrometric Analysis to Elucidate the Mechanisms of Phospho-priming and Auto-activation of the Checkpoint Kinase Rad53 in Vivo

The cell cycle checkpoint kinases play central roles in the genome maintenance of eukaryotes. Activation of the yeast checkpoint kinase Rad53 involves Rad9 or Mrc1 adaptor-mediated phospho-priming by Mec1 kinase, followed by auto-activating phosphorylation within its activation loop. However, the mechanisms by which these adaptors regulate priming phosphorylation of specific sites and how this then leads to Rad53 activation remain poorly understood. Here we used quantitative mass spectrometry to delineate the stepwise phosphorylation events in the activation of endogenous Rad53 in response to S phase alkylation DNA damage, and we show that the two Rad9 and Mrc1 adaptors, the four N-terminal Mec1-target TQ sites of Rad53 (Rad53-SCD1), and Rad53-FHA2 coordinate intimately for optimal priming phosphorylation to support substantial Rad53 auto-activation. Rad9 or Mrc1 alone can mediate surprisingly similar Mec1 target site phosphorylation patterns of Rad53, including previously undetected tri- and tetraphosphorylation of Rad53-SCD1. Reducing the number of TQ motifs turns the SCD1 into a proportionally poorer Mec1 target, which then requires the presence of both Mrc1 and Rad9 for sufficient priming and auto-activation. The phosphothreonine-interacting Rad53-FHA domains, particularly FHA2, regulate phospho-priming by interacting with the checkpoint mediators but do not seem to play a major role in the phospho-SCD1-dependent auto-activation step. Finally, mutation of all four SCD1 TQ motifs greatly reduces Rad53 activation but does not eliminate it, and residual Rad53 activity in this mutant is dependent on Rad9 but not Mrc1. Altogether, our results provide a paradigm for how phosphorylation site clusters and checkpoint mediators can be involved in the regulation of signaling relay in protein kinase cascades in vivo and elucidate an SCD1-independent Rad53 auto-activation mechanism through the Rad9 pathway. The work also demonstrates the power of mass spectrometry for in-depth analyses of molecular mechanisms in cellular signaling in vivo.Eukaryotic cells are most vulnerable to exogenous DNA-damaging agents during the S phase of the cell cycle, when unprogrammed DNA lesions interfere with the tightly choreographed DNA replication process. DNA damage during this phase leads to the activation of two overlapping checkpoint pathways in Saccharomyces cerevisiae, the DNA replication checkpoint and the intra-S-phase DNA damage checkpoint (1, 2). Phospho-priming for auto-activation of the central checkpoint kinase Rad53 by the upstream kinase Mec1/Tel1 depends on Mrc1 as an adaptor in the DNA replication checkpoint pathway and Rad9 as an adaptor in the DNA damage checkpoint pathway (3–10). Rad53, a well-accepted model system for studying the function and regulation of Chk2-like kinases, contains two forkhead-associated (FHA)¹ domains (FHA1 and -2) and two SQ/TQ cluster domains (SCD1 and -2) enriched in Mec1/Tel1-target phosphorylation sites (11–13).Mrc1 normally is a replisome component that functionally couples DNA Pol ε with Cdc45 and MCM helicase during replication fork progression (14, 15). As the replication forks are stalled by replication stress, the recruited checkpoint sensor kinase Mec1 phosphorylates the SCD of Mrc1, which abolishes its N-terminal interaction with Pol ε and enables Mrc1 to recruit Rad53 and promote Rad53 phosphorylation by Mec1 as an initial step in the activation of Rad53 in the Mrc1 branch (6, 14, 16). Alanine substitution of all Mec1 target sites of Mrc1 (designated the mrc1-AQ allele) has been shown to selectively disable its checkpoint function for Rad53 activation without affecting its DNA replication functions (4). In response to DNA damage, Rad9 is able to associate with damaged chromatin via its BRCT and Tudor domains, which tether it to Ser129-phosphorylated histone H2A (γH2A) and Lys79-methylated histone H3, respectively (17, 18). Alternatively, the recruitment of Rad9 onto damaged DNA could also be facilitated by its phosphorylation by CDK1, which enables the specific interaction of Rad9 with Dpb11, allowing the formation of the ternary complex of Dpb11, Mec1, and Rad9 (19, 20). Similar to Mrc1, Mec1 activates the adaptor function of Rad9 by phosphorylation of its SCD, which then binds to the Rad53-FHA domains to promote Rad53 phosphorylation by Mec1 (3, 5, 10).Beyond serving as scaffolds to recruit Rad53, Mrc1 and Rad9 have been shown to promote Rad53 phosphorylation by Mec1 in a dose-dependent manner in vitro (3, 16), underlining their adaptor role to enhance the enzyme–substrate (Mec1–Rad53) interaction. However, how they can specifically regulate the priming phosphorylation at specific sites and how this then leads to Rad53 activation remains poorly understood. Finally, hyperphosphorylated Rad9 has also been shown to catalyze the auto-phosphorylation of recombinant Rad53 (21), but it remains to be examined whether and how this occurs in vivo.The activation of SCD-FHA containing kinases such as human Chk2 and fission yeast Cds1 has been suggested to involve a two-step phosphorylation process: first, SCD phosphorylation by an ATM/ATR-like kinase leads to intermolecular binding to the FHA domain of another Chk2/Cds1 monomer, which then results in dimerization/oligomerization-dependent auto-phosphorylation within the kinase activation loop (22–26). In addition to the characteristic N-terminal SCD-FHA module of Chk2-like kinases, Rad53 contains another SCD2-FHA2 module C-terminal to its kinase domain. Similar to its orthologues, Rad53 activation has been proposed to depend on SCD1 phosphorylation (but not SCD2 phosphorylation) and partially redundant functions of the two FHA domains (9, 27–29). However, although Rad53-FHA1 can interact with SCD1 in a phospho-threonine (pT)-dependent manner in vitro (9, 28), it appears to be required for Rad53 activation only in G2/M-arrested cells (27, 29). In contrast, the FHA2 domain, which seems to be more important overall for Rad53 activation, does not appreciably bind phospho-SCD1 peptides in vitro (27, 28). Thus, the mechanisms by which Mrc1, Rad9, SCD1 phosphorylation, and FHA domains interact during checkpoint-dependent Rad53 priming and auto-activation remain to be elucidated.Quantitative mass spectrometric analysis has revolutionized the functional analysis of cellular signaling pathways, including site-specific phosphorylation events of key signaling molecules (30–33), but an important caveat is that MS studies often involve protein tags or nonphysiological expression levels that can interfere with normal protein functions. For example, the integration of a triple HA tag into the endogenous RAD53 gene locus has been shown to reduce Rad53 protein levels, resulting in significantly altered checkpoint activity (34). In this study we used quantitative MS analyses to dissect the stepwise phosphorylation events of endogenous, untagged Rad53 in response to MMS-induced alkylation DNA damage and replication stress during the S phase. Together with functional analyses, our results delineate how the two Mec1 adaptors Rad9 and Mrc1 can coordinate with the four SCD1 priming sites (T5, T8, T12, and T15) to regulate the phospho-priming of Rad53 by Mec1. In addition, an SCD1-priming independent Rad53 auto-activation mechanism and the specific roles of the FHA domains during Rad53 hyperphosphorylation are also elucidated in this work.     

Tyrosine Phosphorylation of Growth Factor Receptor-bound Protein-7 by Focal Adhesion Kinase in the Regulation of Cell Migration, Proliferation, and Tumorigenesis

We have previously reported that growth factor receptor-bound protein-7 (Grb7), an Src-homology 2 (SH2)-containing adaptor protein, enables interaction with focal adhesion kinase (FAK) to regulate cell migration in response to integrin activation. To further elucidate the signaling events mediated by FAK·Grb7 complexes in promoting cell migration and other cellular functions, we firstly examined the phos pho ryl a ted tyrosine site(s) of Grb7 by FAK using an in vivo mutagenesis. We found that FAK was capable of phos pho rylating at least 2 of 12 tyrosine residues within Grb7, Tyr-188 and Tyr-338. Moreover, mutations converting the identified Tyr to Phe inhibited integrin-dependent cell migration as well as impaired cell proliferation but not survival compared with the wild-type control. Interestingly, the above inhibitory effects caused by the tyrosine phos pho ryl a tion-deficient mutants are probably attributed to their down-regulation of phospho-Tyr-397 of FAK, thereby implying a mechanism by competing with wild-type Grb7 for binding to FAK. Consequently, these tyrosine phos pho ryl a tion-deficient mutants evidently altered the phospho-Tyr-118 of paxillin and phos pho ryl a tion of ERK1/2 but less on phospho-Ser-473 of AKT, implying their involvement in the FAK·Grb7-mediated cellular functions. Additionally, we also illustrated that the formation of FAK·Grb7 complexes and Grb7 phos pho ryl a tion by FAK in an integrin-dependent manner were essential for cell migration, proliferation and anchorage-independent growth in A431 epidermal carcinoma cells, indicating the importance of FAK·Grb7 complexes in tumorigenesis. Our data provide a better understanding on the signal transduction event for FAK·Grb7-mediated cellular functions as well as to shed light on a potential therapeutic in cancers.Growth factor receptor bound protein-7 (Grb7)² is initially identified as a SH2 domain-containing adaptor protein bound to the activated EGF receptor (1). Grb7 is composed of an N-terminal proline-rich region, following a putative RA (Ras-associating) domain and a central PH (pleckstrin homology) domain and a BPS motif (between PH and SH2 domains), and a C-terminal SH2 domain (2–6). Despite the lack of enzymatic activity, the presence of multiple protein-protein interaction domains allows Grb7 family adaptor proteins to participate in versatile signal transduction pathways and, therefore, to regulate many cellular functions (4–6). A number of signaling molecules has been reported to interact with these featured domains, although most of the identified Grb7 binding partners are mediated through its SH2 domain. For example, the SH2 domain of Grb7 has been demonstrated to be capable of binding to the phospho-tyrosine sites of EGF receptor (1), ErbB2 (7), ErbB3 and ErbB4 (8), Ret (9), platelet-derived growth factor receptor (10), insulin receptor (11), SHPTP2 (12), Tek/Tie2 (13), caveolin (14), c-Kit (15), EphB1 (16), G6f immunoreceptor protein (17), Rnd1 (18), Shc (7), FAK (19), and so on. The proceeding α-helix of the PH domain of Grb7 is the calmodulin-binding domain responsible for recruiting Grb7 to plasma membrane in a Ca²⁺-dependent manner (20), and the association between the PH domain of Grb7 and phosphoinositides is required for the phosphorylation by FAK (21). Two additional proteins, NIK (nuclear factor κB-inducing kinase) and FHL2 (four and half lim domains isoform 2), in association with the GM region (Grb and Mig homology region) of Grb7 are also reported, although the physiological functions for these interactions remain unknown (22, 23). Recently, other novel roles in translational controls and stress responses through the N terminus of Grb7 are implicated for the findings of Grb7 interacting with the 5′-untranslated region of capped targeted KOR (kappa opioid receptor) mRNA and the Hu antigen R of stress granules in an FAK-mediated phosphorylation manner (24, 25).Unlike its member proteins Grb10 and Grb14, the role of Grb7 in cell migration is unambiguous and well documented. This is supported by a series of studies. Firstly, Grb7 family members share a significantly conserved molecular architecture with the Caenorhabditis elegans Mig-10 protein, which is involved in neuronal cell migration during embryonic development (4, 5, 26), suggesting that Grb7 may play a role in cell migration. Moreover, Grb7 is often co-amplified with Her2/ErbB2 in certain human cancers and tumor cell lines (7, 27, 28), and its overexpression resulted in invasive and metastatic consequences of various cancers and tumor cells (23, 29–33). On the contrary, knocking down Grb7 by RNA interference conferred to an inhibitory outcome of the breast cancer motility (34). Furthermore, interaction of Grb7 with autophosphorylated FAK at Tyr-397 could promote integrin-mediated cell migration in NIH 3T3 and CHO cells, whereas overexpression of its SH2 domain, an dominant negative mutant of Grb7, inhibited cell migration (19, 35). Recruitment and phosphorylation of Grb7 by EphB1 receptors enhanced cell migration in an ephrin-dependent manner (16). Recently, G7–18NATE, a selective Grb7-SH2 domain affinity cyclic peptide, was demonstrated to efficiently block cell migration of tumor cells (32, 36). In addition to cell migration, Grb7 has been shown to play a role in a variety of physiological and pathological events, for instance, kidney development (37), tumorigenesis (7, 14, 38–41), angiogenic activity (20), proliferation (34, 42, 43), anti-apoptosis (44), gene expression regulation (24), Silver-Russell syndrome (45), rheumatoid arthritis (46), atopic dermatitis (47), and T-cell activation (17, 48). Nevertheless, it remains largely unknown regarding the downstream signaling events of Grb7-mediated various functions. In particular, given the role of Grb7 as an adaptor molecule and its SH2 domain mainly interacting with upstream regulators, it will be interesting to identify potential downstream effectors through interacting with the functional GM region or N-terminal proline-rich region.In this report, we identified two tyrosine phosphorylated sites of Grb7 by FAK and deciphered the signaling targets downstream through these phosphorylated tyrosine sites to regulate various cellular functions such as cell migration, proliferation, and survival. In addition, our study sheds light on tyrosine phosphorylation of Grb7 by FAK involved in tumorigenesis.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Dbf4-Cdc7 Phosphorylation of Mcm2 Is Required for Cell Growth

The Dbf4-Cdc7 kinase (DDK) is required for the activation of the origins of replication, and DDK phosphorylates Mcm2 in vitro. We find that budding yeast Cdc7 alone exists in solution as a weakly active multimer. Dbf4 forms a likely heterodimer with Cdc7, and this species phosphorylates Mcm2 with substantially higher specific activity. Dbf4 alone binds tightly to Mcm2, whereas Cdc7 alone binds weakly to Mcm2, suggesting that Dbf4 recruits Cdc7 to phosphorylate Mcm2. DDK phosphorylates two serine residues of Mcm2 near the N terminus of the protein, Ser-164 and Ser-170. Expression of mcm2-S170A is lethal to yeast cells that lack endogenous MCM2 (mcm2Δ); however, this lethality is rescued in cells harboring the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Mcm2 is required for cell growth.The Cdc7 protein kinase is required throughout the yeast S phase to activate origins (1, 2). The S phase cyclin-dependent kinase also activates yeast origins of replication (3–5). It has been proposed that Dbf4 activates Cdc7 kinase in S phase, and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6). However, it is not known how Dbf4-Cdc7 (DDK)² acts during S phase to trigger the initiation of DNA replication. DDK has homologs in other eukaryotic species, and the role of Cdc7 in activation of replication origins during S phase may be conserved (7–10).The Mcm2-7 complex functions with Cdc45 and GINS to unwind DNA at a replication fork (11–15). A mutation of MCM5 (mcm5-bob1) bypasses the cellular requirements for DBF4 and CDC7 (16), suggesting a critical physiologic interaction between Dbf4-Cdc7 and Mcm proteins. DDK phosphorylates Mcm2 in vitro with proteins purified from budding yeast (17, 18) or human cells (19). Furthermore, there are mutants of MCM2 that show synthetic lethality with DBF4 mutants (6, 17), suggesting a biologically relevant interaction between DBF4 and MCM2. Nevertheless, the physiologic role of DDK phosphorylation of Mcm2 is a matter of dispute. In human cells, replacement of MCM2 DDK-phosphoacceptor residues with alanines inhibits DNA replication, suggesting that Dbf4-Cdc7 phosphorylation of Mcm2 in humans is important for DNA replication (20). In contrast, mutation of putative DDK phosphorylation sites at the N terminus of Schizosaccharomyces pombe Mcm2 results in viable cells, suggesting that phosphorylation of S. pombe Mcm2 by DDK is not critical for cell growth (10).In budding yeast, Cdc7 is present at high levels in G₁ and S phase, whereas Dbf4 levels peak in S phase (18, 21, 22). Furthermore, budding yeast DDK binds to chromatin during S phase (6), and it has been shown that Dbf4 is required for Cdc7 binding to chromatin in budding yeast (23, 24), fission yeast (25), and Xenopus (9). Human and fission yeast Cdc7 are inert on their own (7, 8), but Dbf4-Cdc7 is active in phosphorylating Mcm proteins in budding yeast (6, 26), fission yeast (7), and human (8, 10). Based on these data, it has been proposed that Dbf4 activates Cdc7 kinase in S phase and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6, 9, 18, 21–24). However, a mechanistic analysis of how Dbf4 activates Cdc7 has not yet been accomplished. For example, the multimeric state of the active Dbf4-Cdc7 complex is currently disputed. A heterodimer of fission yeast Cdc7 (Hsk1) in complex with fission yeast Dbf4 (Dfp1) can phosphorylate Mcm2 (7). However, in budding yeast, oligomers of Cdc7 exist in the cell (27), and Dbf4-Cdc7 exists as oligomers of 180 and 300 kDa (27).DDK phosphorylates the N termini of human Mcm2 (19, 20, 28), human Mcm4 (10), budding yeast Mcm4 (26), and fission yeast Mcm6 (10). Although the sequences of the Mcm N termini are poorly conserved, the DDK sites identified in each study have neighboring acidic residues. The residues of budding yeast Mcm2 that are phosphorylated by DDK have not yet been identified.In this study, we find that budding yeast Cdc7 is weakly active as a multimer in phosphorylating Mcm2. However, a low molecular weight form of Dbf4-Cdc7, likely a heterodimer, has a higher specific activity for phosphorylation of Mcm2. Dbf4 or DDK, but not Cdc7, binds tightly to Mcm2, suggesting that Dbf4 recruits Cdc7 to Mcm2. DDK phosphorylates two serine residues of Mcm2, Ser-164 and Ser-170, in an acidic region of the protein. Mutation of Ser-170 is lethal to yeast cells, but this phenotype is rescued by the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Ser-170 of Mcm2 is required for budding yeast growth.