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.     

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.     

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.     

Absolute Proteome and Phosphoproteome Dynamics during the Cell Cycle of Schizosaccharomyces pombe (Fission Yeast)

To quantify cell cycle-dependent fluctuations on a proteome-wide scale, we performed integrative analysis of the proteome and phosphoproteome during the four major phases of the cell cycle in Schizosaccharomyces pombe. In highly synchronized cells, we identified 3753 proteins and 3682 phosphorylation events and relatively quantified 65% of the data across all phases. Quantitative changes during the cell cycle were infrequent and weak in the proteome but prominent in the phosphoproteome. Protein phosphorylation peaked in mitosis, where the median phosphorylation site occupancy was 44%, about 2-fold higher than in other phases. We measured copy numbers of 3178 proteins, which together with phosphorylation site stoichiometry enabled us to estimate the absolute amount of protein-bound phosphate, as well as its change across the cell cycle. Our results indicate that 23% of the average intracellular ATP is utilized by protein kinases to phosphorylate their substrates to drive regulatory processes during cell division. Accordingly, we observe that phosphate transporters and phosphate-metabolizing enzymes are phosphorylated and therefore likely to be regulated in mitosis.Cell replication involves a complex series of highly regulated and evolutionary conserved events, called the “cell cycle.” Aberrations in the cell cycle have severe implications and can cause cancerous growth. A detailed understanding of the cell cycle and its regulation may identify additional targets for cancer therapy (1–3). The cell cycle has been the subject of previous proteomics studies. Olsen et al. (4) measured the dynamics of thousands of proteins and phosphorylation events across cell cycle phases of HeLa cells, providing insights into the underlying regulatory mechanisms and pointing to a general increase in phosphorylation site occupancy during M phase. In a targeted study, Pagliuca et al. (5) investigated interactors of cyclins E1, A2, and B1 in HeLa cells, revealing key mechanistic links between DNA replication and mitosis.Schizosaccharomyces pombe (fission yeast) is a unicellular organism, which can easily be genetically manipulated and carries many cell cycle features similar to metazoan cells. It is an important model organism to study the cell cycle and its checkpoint controls (6). Recent global proteomics studies of yeasts and their cell cycle (7–13) have mainly focused on Saccharomyces cerevisiae (budding yeast), with only a few studies of fission yeast (14, 15), although the fission yeast cell cycle may be more representative of eukaryotic cell cycles (16). However, attention of the proteomics community toward S. pombe is increasing. Recent proteomics studies covered up to 4087 S. pombe proteins (71% of the predicted proteome) and 1544 phosphoproteins in both asynchronous and synchronized cell cultures (17–22); however, a comprehensive analysis of the S. pombe cell cycle is so far missing.Here, we use high resolution mass spectrometry in combination with stable isotope labeling by amino acids in the cell culture (SILAC)¹ method, termed super-SILAC (23), and intensity-based absolute quantification (iBAQ) (24) to measure relative and absolute dynamics of the proteome and phosphoproteome during the cell cycle of fission yeast. We estimate copy numbers for 3178 S. pombe proteins, and we combine these data with calculated phosphorylation site stoichiometry to estimate the total amount of protein-bound phosphate and its dynamics across the cell cycle. Providing the global absolute dynamics and stoichiometry of proteins and their modifications will be a valuable resource for classical and systems biologists alike.     

Phosphoproteomics of Klebsiella pneumoniae NTUH-K2044 Reveals a Tight Link between Tyrosine Phosphorylation and Virulence

Encapsulated Klebsiella pneumoniae is the predominant causative agent of pyogenic liver abscess, an emerging infectious disease that often complicates metastatic meningitis or endophthalmitis. The capsular polysaccharide on K. pneumoniae surface was determined as the key to virulence. Although the regulation of capsular polysaccharide biosynthesis is largely unclear, it was found that protein-tyrosine kinases and phosphatases are involved. Therefore, the identification and characterization of such kinases, phosphatases, and their substrates would advance our knowledge of the underlying mechanism in capsule formation and could contribute to the development of new therapeutic strategies. Here, we analyzed the phosphoproteome of K. pneumoniae NTUH-K2044 with a shotgun approach and identified 117 unique phosphopeptides along with 93 in vivo phosphorylated sites corresponding to 81 proteins. Interestingly, three of the identified tyrosine phosphorylated proteins, namely protein-tyrosine kinase (Wzc), phosphomannomutase (ManB), and undecaprenyl-phosphate glycosyltransferase (WcaJ), were found to be distributed in the cps locus and thus were speculated to be involved in the converging signal transduction of capsule biosynthesis. Consequently, we decided to focus on the lesser studied ManB and WcaJ for mutation analysis. The capsular polysaccharides of WcaJ mutant (WcaJY5F) were dramatically reduced quantitatively, and the LD₅₀ increased by 200-fold in a mouse peritonitis model compared with the wild-type strain. However, the capsular polysaccharides of ManB mutant (ManBY26F) showed no difference in quantity, and the LD₅₀ increased by merely 6-fold in mice test. Our study provided a clear trend that WcaJ tyrosine phosphorylation can regulate the biosynthesis of capsular polysaccharides and result in the pathogenicity of K. pneumoniae NTUH-K2044.Protein phosphorylation is one of the most biologically relevant and ubiquitous post-translational modifications in both eukaryotic and prokaryotic organisms. It is best known that protein phosphorylation is a reversible enzyme-catalyzed process that is controlled by various kinases and phosphatases. The aberrant functions often result in irregular protein phosphorylation and ultimately lead to serious disease states such as malignant transformation, immune disorders, and pathogenic infections in mammals (1, 2). Recently, accumulating evidences suggest that Ser/Thr/Tyr phosphorylations also contribute to regulate a diverse range of cellular responses and physiological processes in prokaryotes (1). Among them, tyrosine phosphorylation in encapsulated bacteria has been discovered to play key roles in capsular polysaccharide (CPS1; K antigen) biosynthesis, which leads to virulence (3, 4). This thick layer of exopolysaccharide on many pathogenic bacteria can act as a physical boundary to evade phagocytosis and complement-mediated killing and further inhibit complement activation of the host (1, 5, 6).In 1996, Acinetobacter johnsonii protein-tyrosine kinase (Ptk) was first discovered and categorized under the bacterial protein-tyrosine kinase (BY-kinase) family (1, 7, 8). Shortly after, its function in bacterial exopolysaccharide production and transport was characterized (1, 7, 8). From then on, many more bacterial tyrosine kinases such as Wzc of Escherichia coli (1, 9) and EpsB of Pseudomonas solanacearum (10, 11) were found to possess this conserved property; deletion of such tyrosine kinases will result in the loss of exopolysaccharide production (12). Therefore, several experiments were conducted to investigate the role of the downstream substrates of the tyrosine kinases in different strains of bacteria, and some targeted proteins were found to participate in the exopolysaccharide anabolism (13, 14). These findings demonstrated a direct relationship between bacterial tyrosine phosphorylation and exopolysaccharide biosynthesis that was directly reflected in the strain virulence.In the past, the functional roles of the critical components involved in protein phosphorylation were defined by basic biochemical and genetic approaches (1). However, there exists a salient gap between the growing number of identified protein-tyrosine kinases/phosphatases and the relative paucity of protein substrates characterized to date. Genomic sequence analyses and advanced high resolution/high accuracy MS systems with vastly improved phosphopeptide enrichment strategies are among the two key enabling technologies that allow a high efficiency identification of the scarcely detectable site-specific phosphorylations in bacterial systems (15). Mann et al. (16) were the first to initiate a systematic study of the phosphoproteome of B. subtilis in 2007 followed by similar site-specific phosphoproteomics analyses of E. coli (17), Lactococcus lactis (18), and Halobacterium salinarum (19). These pioneering works have since set the foundation in bacterial phosphoproteomics but have not been specifically carried out to address a particular biological issue of causal relevance to virulence or pathogenesis.Klebsiella pneumoniae is a Gram-negative, non-motile, facultative anaerobic, and rod-shaped bacterium. It is commonly found in water and soil (20) as well as on plants (21) and mucosal surfaces of mammals, such as human, horse, and swine (22, 23). It was demonstrated that CPS on the surface of K. pneumoniae is the prime factor of virulence and toxicity in causing pyogenic liver abscess (PLA), a common intra-abdominal infection with a high 10–30% mortality rate worldwide (24–29). There are also variations in virulence in regard to different capsular serotypes; K1 and K2 were found to be especially pathogenic in causing PLA in a mouse model (30) compared with other serotypes, which show little or no effect (31–34). The K. pneumoniae NTUH-K2044 (K2044) strain, encapsulated with K1 antigen (35), was isolated from clinical K. pneumoniae liver abscess patients. It has become an important emerging pathogen (36) because it usually complicates metastatic septic endophthalmitis and irreversible central nervous system infections independent of host underlying diseases (30, 34). The transmission rate is high (37), and it often rapidly leads to outbreaks of community-acquired infections, such as bacteremia, nosocomial pneumonia, and sepsis, common in immunocompromised individuals (38).In this study, we wanted to prove that the biosynthesis of CPS is mediated through tyrosine phosphorylation of a subset of proteins. An MS-based systematic phosphoproteomics analysis was conducted on K2044 to identify tyrosine phosphorylated proteins that are also associated with CPS biosynthesis. We further validated the relationship between tyrosine phosphorylation on those proteins and virulence of K2044 by site-directed mutagenesis, CPS quantification, serum killing, and mouse lethality assay.     

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]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

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.     

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]     

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.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[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]     

Environmental Stress Affects the Activity of Metabolic and Growth Factor Signaling Networks and Induces Autophagy Markers in MCF7 Breast Cancer Cells

Phosphoproteomic techniques are contributing to our understanding of how signaling pathways interact and regulate biological processes. This technology is also being used to characterize how signaling networks are remodeled during disease progression and to identify biomarkers of signaling pathway activity and of responses to cancer therapy. A potential caveat in these studies is that phosphorylation is a very dynamic modification that can substantially change during the course of an experiment or the retrieval and processing of cellular samples. Here, we investigated how exposure of cells to ambient conditions modulates phosphorylation and signaling pathway activity in the MCF7 breast cancer cell line. About 1.5% of 3,500 sites measured showed a significant change in phosphorylation extent upon exposure of cells to ambient conditions for 15 min. The effects of this perturbation in modifying phosphorylation patterns did not involve random changes due to stochastic activation of kinases and phosphatases. Instead, exposure of cells to ambient conditions elicited an environmental stress reaction that involved a coordinated response to a metabolic stress situation, which included: (1) the activation of AMPK; (2) the inhibition of PI3K, AKT, and ERK; (3) an increase in markers of protein synthesis inhibition at the level of translation elongation; and (4) an increase in autophagy markers. We also observed that maintaining cells in ice modified but did not completely abolish this metabolic stress response. In summary, exposure of cells to ambient conditions affects the activity of signaling networks previously implicated in metabolic and growth factor signaling. Mass spectrometry data have been deposited to the ProteomeXchange with identifier PXD000472.Phosphorylation is a posttranslational modification involving the addition of phosphate groups to serine, threonine and tyrosine residues on target proteins. This modification, regulated by kinases and phosphatases that phosphorylate and dephosphorylate these amino acids respectively, controls many aspects of protein biochemistry including stability, localization, ability to interact with other molecules and enzymatic activity (1, 2). In addition to playing a pivotal role in regulating most biological processes, alterations in biochemical pathways regulated by protein phosphorylation contribute to the pathophysiology of various diseases including cancer, diabetes and neurodegeneration (2–6).In recent years the development of MS techniques has allowed the study of protein phosphorylation on an untargeted and global scale. As a consequence, signaling processes can now be studied with unprecedented depth and coverage (7–10). Phosphoproteomics has also been applied to investigate how signaling networks are modulated during disease progression and for the identification of biomarkers that classify patients according to prognosis or treatment response (11–15). A potential caveat in the interpretation of such experiments is that protein phosphorylation is a dynamic modification that can be affected by variables difficult to control including cell confluence, circadian rhythms, shear stress and other types of environmental stresses including exposure to ambient conditions (16–22). Thus, during the course of an experiment variations or delays in sample retrieval and processing can potentially alter the quantitative characteristics of the phosphoproteome (17, 18, 22). Similar problems could in principle occur in a clinical environment where several hours may elapse from patient sample collection to processing or preservation (16, 17, 23). Delays because of ethical and practical considerations may also affect collection and preservation of post-mortem samples (24, 25). As a consequence, it can in principle be introduced variability and artifacts that may potentially confound the interpretation of data obtained from large-scale as well as targeted phosphoproteomics experiments (16).To our knowledge, there are no reports that systematically evaluate, in an untargeted manner, how exposure to environmental stress modulates the phosphoproteome of human cells in culture. Here, we used the MCF7 breast cancer cell line to investigate how ambient conditions alter phosphorylation and to evaluate signaling pathways that may be modulated by environmental stress. We found several phosphorylation events that increased or decreased after 15 min exposure of cells to ambient conditions at room temperature (RT)¹. We then studied whether these changes in phosphorylation were a random effect due to stochastic inactivation of kinases and phosphatases or whether these were the consequence of actual responses involving specific signaling pathways. Our data indicate that the phosphorylations regulated by environmental conditions at RT are the initial steps of a complex adaptive response to a metabolic stress. Data supporting these conclusions include the observation that ambient conditions at RT activated catabolic pathways regulated by AMPK and GSK3β and inactivated anabolic pathways involving the AKT, ERK and mTOR signaling nodes. When we compared the responses to ambient conditions at RT or on ice, we found that maintaining cells on ice induced a different adaptive response rather than an attenuated one. We also found that the adaptation response to ambient conditions at RT triggered a functional biological process that involved the initiation of macroautophagy (hereafter referred as autophagy) and the activation of a pathway known to inhibit protein synthesis at the level of translation elongation. Thus our study also defines experimental conditions that can be used to study the mechanisms involved in the process of autophagy.