The Phosphoproteome of Bloodstream Form Trypanosoma brucei, Causative Agent of African Sleeping Sickness

The protozoan parasite Trypanosoma brucei is the causative agent of human African sleeping sickness and related animal diseases, and it has over 170 predicted protein kinases. Protein phosphorylation is a key regulatory mechanism for cellular function that, thus far, has been studied in T.brucei principally through putative kinase mRNA knockdown and observation of the resulting phenotype. However, despite the relatively large kinome of this organism and the demonstrated essentiality of several T. brucei kinases, very few specific phosphorylation sites have been determined in this organism. Using a gel-free, phosphopeptide enrichment-based proteomics approach we performed the first large scale phosphorylation site analyses for T.brucei. Serine, threonine, and tyrosine phosphorylation sites were determined for a cytosolic protein fraction of the bloodstream form of the parasite, resulting in the identification of 491 phosphoproteins based on the identification of 852 unique phosphopeptides and 1204 phosphorylation sites. The phosphoproteins detected in this study are predicted from their genome annotations to participate in a wide variety of biological processes, including signal transduction, processing of DNA and RNA, protein synthesis, and degradation and to a minor extent in metabolic pathways. The analysis of phosphopeptides and phosphorylation sites was facilitated by in-house developed software, and this automated approach was validated by manual annotation of spectra of the kinase subset of proteins. Analysis of the cytosolic bloodstream form T. brucei kinome revealed the presence of 44 phosphorylated protein kinases in our data set that could be classified into the major eukaryotic protein kinase groups by applying a multilevel hidden Markov model library of the kinase catalytic domain. Identification of the kinase phosphorylation sites showed conserved phosphorylation sequence motifs in several kinase activation segments, supporting the view that phosphorylation-based signaling is a general and fundamental regulatory process that extends to this highly divergent lower eukaryote.Phosphorylation by eukaryotic protein kinases (PKs)¹ is the most intensively studied posttranslational modification in eukaryotic cells mainly because reversible phosphorylation regulates almost every cellular process (1). PKs make up one of the largest protein superfamilies, comprising the protein-serine/threonine kinases and protein-tyrosine kinases (2).Knowledge of the fundamental role of protein kinases for cell survival in multicellular organisms has led to functional studies of this enzyme class in parasitic protozoa to search for alternative drug targets to treat tropical diseases (3). Trypanosoma brucei causes human African trypanosomiasis (also known as African sleeping sickness) and is responsible for ∼30,000 deaths per annum (4). Initial bioinformatics analysis of the T. brucei genome revealed the presence of 156 conventional protein kinases (ePks) and 20 atypical protein kinases (aPKs) (5). Our characterization of the T. brucei genome using a multilevel hidden Markov model (HMM) library of the kinase catalytic domain (6) followed by manual curation identified 170 ePKs and 12 aPKs.Over the past few years, research on T. brucei protein kinases has largely focused on cell cycle-regulating enzymes because of their potential as targets for African sleeping sickness (3, 7). Indeed several protein kinases have been shown to be essential for the organism (8–12). Most of the above mentioned protein kinases were characterized using techniques such as RNA interference, overexpression of the kinase and/or a kinase-dead version, in vitro and in vivo kinase assays, and immunofluorescence microscopy. However, despite these encouraging investigations, the literature for T. brucei proteins with a predicted role in signal transduction lacks information about their in vivo phosphorylation sites. Knowledge of protein phosphorylation sites, however, is crucial for a complete functional characterization of protein kinases and kinase-substrate relationships (1, 13).At the start of this study, only seven phosphorylation sites had been mapped in any of the trypanosomatids: Mehlert et al. (14) had found that six of seven threonine residues of the Gly-Pro-Glu-Glu-Thr repeat-procyclin are phosphorylated using MALDI mass spectrometry, and da Cunha et al. (15) had determined one phosphorylation site at Ser¹² of Trypanosoma cruzi histone H1 protein by electrospray MS. Over the past few years phosphoproteome studies in other organisms using phosphopeptide enrichment in combination with high accuracy mass spectrometry have revealed thousands of novel protein phosphorylation sites and given new insights into the extent and function of this posttranslational modification (16–20).In this study, we used a gel-free peptide-based phosphoproteomics approach to identify 491 trypanosomal phosphoproteins and, where possible, their precise phosphorylation sites. In addition to phosphorylation on serine and threonine residues, we identified several phosphorylated tyrosine residues in T. brucei proteins despite the lack of conventional tyrosine kinases in the T. brucei genome (5). Furthermore phosphorylation in the activation segment of protein kinases was found to be conserved in T. brucei, suggesting that at least some signaling pathways described in other eukaryotes are conserved in this ancient and highly divergent eukaryote.     

F-Actin Dynamics in Neurospora crassa

This study demonstrates the utility of Lifeact for the investigation of actin dynamics in Neurospora crassa and also represents the first report of simultaneous live-cell imaging of the actin and microtubule cytoskeletons in filamentous fungi. Lifeact is a 17-amino-acid peptide derived from the nonessential Saccharomyces cerevisiae actin-binding protein Abp140p. Fused to green fluorescent protein (GFP) or red fluorescent protein (TagRFP), Lifeact allowed live-cell imaging of actin patches, cables, and rings in N. crassa without interfering with cellular functions. Actin cables and patches localized to sites of active growth during the establishment and maintenance of cell polarity in germ tubes and conidial anastomosis tubes (CATs). Recurrent phases of formation and retrograde movement of complex arrays of actin cables were observed at growing tips of germ tubes and CATs. Two populations of actin patches exhibiting slow and fast movement were distinguished, and rapid (1.2 μm/s) saltatory transport of patches along cables was observed. Actin cables accumulated and subsequently condensed into actin rings associated with septum formation. F-actin organization was markedly different in the tip regions of mature hyphae and in germ tubes. Only mature hyphae displayed a subapical collar of actin patches and a concentration of F-actin within the core of the Spitzenkörper. Coexpression of Lifeact-TagRFP and β-tubulin–GFP revealed distinct but interrelated localization patterns of F-actin and microtubules during the initiation and maintenance of tip growth.Actins are highly conserved proteins found in all eukaryotes and have an enormous variety of cellular roles. The monomeric form (globular actin, or G-actin) can self-assemble, with the aid of numerous actin-binding proteins (ABPs), into microfilaments (filamentous actin, or F-actin), which, together with microtubules, form the two major components of the fungal cytoskeleton. Numerous pharmacological and genetic studies of fungi have demonstrated crucial roles for F-actin in cell polarity, exocytosis, endocytosis, cytokinesis, and organelle movement (6, 7, 20, 34, 35, 51, 52, 59). Phalloidin staining, immunofluorescent labeling, and fluorescent-protein (FP)-based live-cell imaging have revealed three distinct subpopulations of F-actin-containing structures in fungi: patches, cables, and rings (1, 14, 28, 34, 60, 63, 64). Actin patches are associated with the plasma membrane and represent an accumulation of F-actin around endocytic vesicles (3, 26, 57). Actin cables are bundles of actin filaments stabilized with cross-linking proteins, such as tropomyosins and fimbrin, and are assembled by formins at sites of active growth, where they form tracks for myosin V-dependent polarized secretion and organelle transport (10, 16, 17, 27, 38, 47, 48). Cables, unlike patches, are absolutely required for polarized growth in the budding yeast Saccharomyces cerevisiae (34, 38). Contractile actomyosin rings are essential for cytokinesis in budding yeast, whereas in filamentous fungi, actin rings are less well studied but are known to be involved in septum formation (20, 28, 34, 39, 40).Actin cables and patches have been particularly well studied in budding yeast. However, there are likely to be important differences between F-actin architecture and dynamics in budding yeast and those in filamentous fungi, as budding yeasts display only a short period of polarized growth during bud formation, which is followed by isotropic growth over the bud surface (10). Sustained polarized growth during hyphal morphogenesis is a defining feature of filamentous fungi (21), making them attractive models for studying the roles of the actin cytoskeleton in cell polarization, tip growth, and organelle transport.In Neurospora crassa and other filamentous fungi, disruption of the actin cytoskeleton leads to rapid tip swelling, which indicates perturbation of polarized tip growth, demonstrating a critical role for F-actin in targeted secretion to particular sites on the plasma membrane (7, 22, 29, 56). Immunofluorescence studies of N. crassa have shown that F-actin localizes to hyphal tips as “clouds” and “plaques” (7, 54, 59). However, immunolabeling has failed to reveal actin cables in N. crassa and offers limited insights into F-actin dynamics. Live-cell imaging of F-actin architecture and dynamics has not been accomplished in N. crassa, yet it is expected to yield key insights into cell polarization, tip growth, and intracellular transport.We took advantage of a recently developed live-cell imaging probe for F-actin called Lifeact (43). Lifeact is a 17-amino-acid peptide derived from the N terminus of the budding yeast actin-binding protein Abp140 (5, 63) and has recently been demonstrated to be a universal live-cell imaging marker for F-actin in eukaryotes (43). Here, we report the successful application of fluorescent Lifeact fusion constructs for live-cell imaging of F-actin in N. crassa. We constructed two synthetic genes consisting of Lifeact fused to “synthetic” green fluorescent protein (sGFP) (S65T) (henceforth termed GFP) (12) or red fluorescent protein (TagRFP) (33) and expressed these constructs in various N. crassa strains. In all strain backgrounds, fluorescent Lifeact constructs clearly labeled actin patches, cables, and rings and revealed a direct association of F-actin structures with sites of cell polarization and active tip growth. Our results demonstrate the efficacy of Lifeact as a nontoxic live-cell imaging probe in N. crassa.     

High-definition De Novo Sequencing of Crustacean Hyperglycemic Hormone (CHH)-family Neuropeptides

A complete understanding of the biological functions of large signaling peptides (>4 kDa) requires comprehensive characterization of their amino acid sequences and post-translational modifications, which presents significant analytical challenges. In the past decade, there has been great success with mass spectrometry-based de novo sequencing of small neuropeptides. However, these approaches are less applicable to larger neuropeptides because of the inefficient fragmentation of peptides larger than 4 kDa and their lower endogenous abundance. The conventional proteomics approach focuses on large-scale determination of protein identities via database searching, lacking the ability for in-depth elucidation of individual amino acid residues. Here, we present a multifaceted MS approach for identification and characterization of large crustacean hyperglycemic hormone (CHH)-family neuropeptides, a class of peptide hormones that play central roles in the regulation of many important physiological processes of crustaceans. Six crustacean CHH-family neuropeptides (8–9.5 kDa), including two novel peptides with extensive disulfide linkages and PTMs, were fully sequenced without reference to genomic databases. High-definition de novo sequencing was achieved by a combination of bottom-up, off-line top-down, and on-line top-down tandem MS methods. Statistical evaluation indicated that these methods provided complementary information for sequence interpretation and increased the local identification confidence of each amino acid. Further investigations by MALDI imaging MS mapped the spatial distribution and colocalization patterns of various CHH-family neuropeptides in the neuroendocrine organs, revealing that two CHH-subfamilies are involved in distinct signaling pathways.Neuropeptides and hormones comprise a diverse class of signaling molecules involved in numerous essential physiological processes, including analgesia, reward, food intake, learning and memory (1). Disorders of the neurosecretory and neuroendocrine systems influence many pathological processes. For example, obesity results from failure of energy homeostasis in association with endocrine alterations (2, 3). Previous work from our lab used crustaceans as model organisms found that multiple neuropeptides were implicated in control of food intake, including RFamides, tachykinin related peptides, RYamides, and pyrokinins (4–6).Crustacean hyperglycemic hormone (CHH)¹ family neuropeptides play a central role in energy homeostasis of crustaceans (7–17). Hyperglycemic response of the CHHs was first reported after injection of crude eyestalk extract in crustaceans. Based on their preprohormone organization, the CHH family can be grouped into two sub-families: subfamily-I containing CHH, and subfamily-II containing molt-inhibiting hormone (MIH) and mandibular organ-inhibiting hormone (MOIH). The preprohormones of the subfamily-I have a CHH precursor related peptide (CPRP) that is cleaved off during processing; and preprohormones of the subfamily-II lack the CPRP (9). Uncovering their physiological functions will provide new insights into neuroendocrine regulation of energy homeostasis.Characterization of CHH-family neuropeptides is challenging. They are comprised of more than 70 amino acids and often contain multiple post-translational modifications (PTMs) and complex disulfide bridge connections (7). In addition, physiological concentrations of these peptide hormones are typically below picomolar level, and most crustacean species do not have available genome and proteome databases to assist MS-based sequencing.MS-based neuropeptidomics provides a powerful tool for rapid discovery and analysis of a large number of endogenous peptides from the brain and the central nervous system. Our group and others have greatly expanded the peptidomes of many model organisms (3, 18–33). For example, we have discovered more than 200 neuropeptides with several neuropeptide families consisting of as many as 20–40 members in a simple crustacean model system (5, 6, 25–31, 34). However, a majority of these neuropeptides are small peptides with 5–15 amino acid residues long, leaving a gap of identifying larger signaling peptides from organisms without sequenced genome. The observed lack of larger size peptide hormones can be attributed to the lack of effective de novo sequencing strategies for neuropeptides larger than 4 kDa, which are inherently more difficult to fragment using conventional techniques (34–37). Although classical proteomics studies examine larger proteins, these tools are limited to identification based on database searching with one or more peptides matching without complete amino acid sequence coverage (36, 38).Large populations of neuropeptides from 4–10 kDa exist in the nervous systems of both vertebrates and invertebrates (9, 39, 40). Understanding their functional roles requires sufficient molecular knowledge and a unique analytical approach. Therefore, developing effective and reliable methods for de novo sequencing of large neuropeptides at the individual amino acid residue level is an urgent gap to fill in neurobiology. In this study, we present a multifaceted MS strategy aimed at high-definition de novo sequencing and comprehensive characterization of the CHH-family neuropeptides in crustacean central nervous system. The high-definition de novo sequencing was achieved by a combination of three methods: (1) enzymatic digestion and LC-tandem mass spectrometry (MS/MS) bottom-up analysis to generate detailed sequences of proteolytic peptides; (2) off-line LC fractionation and subsequent top-down MS/MS to obtain high-quality fragmentation maps of intact peptides; and (3) on-line LC coupled to top-down MS/MS to allow rapid sequence analysis of low abundance peptides. Combining the three methods overcomes the limitations of each, and thus offers complementary and high-confidence determination of amino acid residues. We report the complete sequence analysis of six CHH-family neuropeptides including the discovery of two novel peptides. With the accurate molecular information, MALDI imaging and ion mobility MS were conducted for the first time to explore their anatomical distribution and biochemical properties.     

Increased Lipid Accumulation in the Chlamydomonas reinhardtii sta7-10 Starchless Isoamylase Mutant and Increased Carbohydrate Synthesis in Complemented Strains

The accumulation of bioenergy carriers was assessed in two starchless mutants of Chlamydomonas reinhardtii (the sta6 [ADP-glucose pyrophosphorylase] and sta7-10 [isoamylase] mutants), a control strain (CC124), and two complemented strains of the sta7-10 mutant. The results indicate that the genetic blockage of starch synthesis in the sta6 and sta7-10 mutants increases the accumulation of lipids on a cellular basis during nitrogen deprivation relative to that in the CC124 control as determined by conversion to fatty acid methyl esters. However, this increased level of lipid accumulation is energetically insufficient to completely offset the loss of cellular starch that is synthesized by CC124 during nitrogen deprivation. We therefore investigated acetate utilization and O₂ evolution to obtain further insights into the physiological adjustments utilized by the two starchless mutants in the absence of starch synthesis. The results demonstrate that both starchless mutants metabolize less acetate and have more severely attenuated levels of photosynthetic O₂ evolution than CC124, indicating that a decrease in overall anabolic processes is a significant physiological response in the starchless mutants during nitrogen deprivation. Interestingly, two independent sta7-10:STA7 complemented strains exhibited significantly greater quantities of cellular starch and lipid than CC124 during acclimation to nitrogen deprivation. Moreover, the complemented strains synthesized significant quantities of starch even when cultured in nutrient-replete medium.Microalgae are able to efficiently convert sunlight, water, and CO₂ into a variety of products suitable for renewable energy applications, including H₂, carbohydrates, and lipids (11, 12, 16, 21, 38, 41, 44). The unicellular green alga Chlamydomonas reinhardtii has emerged as a model organism for studying algal physiology, photosynthesis, metabolism, nutrient stress, and the synthesis of bioenergy carriers (12, 15, 19, 24, 32). During acclimation to nitrogen deprivation, C. reinhardtii cells accumulate significant quantities of starch and form lipid bodies (4, 5, 8, 26, 28, 30, 34, 43, 46, 48). Despite the significance of these products in algal physiology and in biofuels applications, the metabolic, enzymatic, and regulatory mechanisms controlling the partitioning of metabolites into these distinct carbon stores in algae are poorly understood. Several C. reinhardtii starch mutants with various phenotypic changes in starch content and structure have been isolated (2,–4). Two of these, the sta6 and sta7 mutants, contain single-gene disruptions that result in “starchless” phenotypes with severely attenuated levels of starch granule accumulation (2, 4, 34, 39, 40, 48).The disrupted loci in the two isolated starchless mutants are distinct and each mutant has a unique phenotype (7, 40). In the sta6 mutant, the small, catalytic subunit of ADP-glucose pyrophosphorylase (AGPase-SS) is disrupted (2, 4, 48), and this mutant accumulates less than 1% of the starch observed in wild-type (WT) cells under conditions of nitrogen deprivation. The sta7 mutant contains a disrupted isoamylase gene (7, 8, 10, 39, 40) and also has severely attenuated levels of starch, but it accumulates a soluble glycogen-like product (4, 9). In this study, we conducted an examination of the unique physiological acclimations that are utilized by these mutants to adapt to the loss of starch synthesis. As the genetic lesions in these two mutants are distinct and block starch synthesis via two very different mechanisms, we investigated the physiological consequences of starch inhibition in both of these mutants from a holistic bioenergy perspective, which included photosynthetic parameters and the overall yields of lipids and carbohydrates, the two primary bioenergy carriers in C. reinhardtii. Specifically, we examined whether the inability to synthesize starch would result in the accumulation of additional lipid, alter cellular growth or cell size, affect acetate utilization, and/or influence photosynthetic O₂ evolution. Our data indicate that both the sta6 (BAFJ5) and sta7 (sta7-10) mutants accumulate more lipid than the CC124 control during nitrogen deprivation. However, the additional lipid does not completely offset the loss of starch synthesis from a complete energetic perspective. Increased lipid accumulation during nitrogen stress has also been reported for a variety of starch mutants in recent papers (26, 27, 46). A significant feature in both of the starchless mutants studied here is that O₂ evolution and acetate utilization are diminished during nitrogen stress, which is undesirable from an overall bioenergy perspective. Remarkably, complementation of sta7-10 with genomic DNA encoding the wild-type isoamylase gene resulted in cells that were larger than those of the sta6, sta7-10, and CC124 strains, exhibited the highest total lipid levels during nitrogen deprivation, and overaccumulated starch even in nutrient-replete medium.     

Identification of Extracellular Signal-regulated Kinase 1 (ERK1) Direct Substrates using Stable Isotope Labeled Kinase Assay-Linked Phosphoproteomics

Kinase mediated phosphorylation signaling is extensively involved in cellular functions and human diseases, and unraveling phosphorylation networks requires the identification of substrates targeted by kinases, which has remained challenging. We report here a novel proteomic strategy to identify the specificity and direct substrates of kinases by coupling phosphoproteomics with a sensitive stable isotope labeled kinase reaction. A whole cell extract was moderately dephosphorylated and subjected to in vitro kinase reaction under the condition in which ¹⁸O-ATP is the phosphate donor. The phosphorylated proteins are then isolated and identified by mass spectrometry, in which the heavy phosphate (+85.979 Da) labeled phosphopeptides reveal the kinase specificity. The in vitro phosphorylated proteins with heavy phosphates are further overlapped with in vivo kinase-dependent phosphoproteins for the identification of direct substrates with high confidence. The strategy allowed us to identify 46 phosphorylation sites on 38 direct substrates of extracellular signal-regulated kinase 1, including multiple known substrates and novel substrates, highlighting the ability of this high throughput method for direct kinase substrate screening.Protein phosphorylation regulates almost all aspects of cell life, such as cell cycle, migration, and apoptosis (1), and deregulation of protein phosphorylation is one of the most frequent causes or consequences of human diseases including cancers, diabetes, and immune disorders (2). Up till now, however, known substrates are far from saturation for the majority of protein kinases (3); thus, mapping comprehensive kinase-substrate relationships is essential to understanding biological mechanisms and uncovering new drug targets (4).Accompanied with advances of high-speed and high-resolution mass spectrometry, the technique of kinase substrate screening using proteomic strategy is quickly evolving (5–7). Mass spectrometry has been extensively used for kinase-substrate interaction mapping (8) and global phosphorylation profiling (9). Although thousands of phosphorylation sites have been detected, complex phosphorylation cascade and crosstalk between pathways make it difficult for large-scale phosphoproteomics to reveal direct relationships between protein kinases and their substrates (10, 11). Extensive statistics, bioinformatics, and downstream biochemical assays are mandatory for the substrate verification (12, 13). Another strategy uses purified, active kinases to phosphorylate cell extracts in vitro, followed by mass spectrometric analysis to identify phosphoproteins. This approach inevitably faces the major challenge of separating real sites phosphorylated by target kinase and the phosphorylation triggered by endogenous kinases from cell lysates (14). Analog-sensitive kinase allele (15) overcomes the issue by utilizing the engineered kinase that can exclusively take a bulky-ATP analog under the reaction condition. Analog-sensitive kinase allele has been coupled with γ-thiophosphate analog ATP to facilitate the mass spectrometric analysis (16–18).We have introduced kinase assay-linked phosphoproteomics (KALIP)¹ to link the in vitro substrate identification and physiological phosphorylation events together in a high throughput manner (19, 20). The strategy, however, has only been applied to identify direct substrates of tyrosine kinases. In this study, we expanded the application of KALIP to serine/threonine kinases by introducing a quantitative strategy termed Stable Isotope Labeled Kinase Assay-Linked Phosphoproteomics (siKALIP). The method was applied to identify direct substrates of extracellular signal-regulated kinase 1 (ERK1), a serine/threonine kinase acting as an essential component of the Mitogen-activated protein kinase (MAPK) signal transduction pathway (21). A defect in the MAP/ERK pathway causes uncontrolled growth, which likely leads to cancer (22) and other diseases (23–25). ERK1 can be activated by growth factors such as platelet-derived growth factor (PDGF), epidermal growth factor (EGF), and nerve growth factor (NGF) (26). Upon stimulation, ERK1 phosphorylates hundreds of substrates in various cellular compartments including cytoplasm, nucleus, and membrane (27). Among 38 ERK1 direct substrates identified by siKALIP, more than one third are previously discovered by classical molecular biology approaches, highlighting high specificity and sensitivity of the strategy. The results also support the hypothesis that ERK1 plays complex roles in multiple pathways that are essential for the cell growth regulation.     

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.     

A Trk/HKT-Type K+ Transporter from Trypanosoma brucei

The molecular mechanisms of K⁺ homeostasis are only poorly understood for protozoan parasites. Trypanosoma brucei subsp. parasites, the causative agents of human sleeping sickness and nagana, are strictly extracellular and need to actively concentrate K⁺ from their hosts’ body fluids. The T. brucei genome contains two putative K⁺ channel genes, yet the trypanosomes are insensitive to K⁺ antagonists and K⁺ channel-blocking agents, and they do not spontaneously depolarize in response to high extracellular K⁺ concentrations. However, the trypanosomes are extremely sensitive to K⁺ ionophores such as valinomycin. Surprisingly, T. brucei possesses a member of the Trk/HKT superfamily of monovalent cation permeases which so far had only been known from bacteria, archaea, fungi, and plants. The protein was named TbHKT1 and functions as a Na⁺-independent K⁺ transporter when expressed in Escherichia coli, Saccharomyces cerevisiae, or Xenopus laevis oocytes. In trypanosomes, TbHKT1 is expressed in both the mammalian bloodstream stage and the Tsetse fly midgut stage; however, RNA interference (RNAi)-mediated silencing of TbHKT1 expression did not produce a growth phenotype in either stage. The presence of HKT genes in trypanosomatids adds a further piece to the enigmatic phylogeny of the Trk/HKT superfamily of K⁺ transporters. Parsimonial analysis suggests that the transporters were present in the first eukaryotes but subsequently lost in several of the major eukaryotic lineages, in at least four independent events.Potassium (K⁺) is the most abundant cation in the cytosol of any cell and hence an essential macronutrient for life on earth. Concentrative K⁺ uptake across the plasma membrane is energized directly by ATPases and indirectly by the negative membrane potential or by coupling, via symport or antiport, to other transport processes such as H⁺ flux. The ancestry of K⁺ transporters renders them ideal subjects for phylogenetic comparisons. Indeed, the different kinds of known K⁺ transporters—pumps, channels, permeases, symporters, and antiporters—are all found in bacteria (43). Eukaryotes do not appear to have invented further mechanisms of K⁺ transport; on the contrary, some families of K⁺ transporters were lost over the course of eukaryote evolution, particularly among the metazoa (53).The Trk/HKT superfamily (TC transporter classification 2.A.38 [43]) consists of bacterial TrkH and KtrB, plant HKT, and fungal Trk transporters (15). These proteins share a topology with 8 transmembrane (TM) domains and, sandwiched between odd- and even-numbered TM domains, 4 shorter hydrophobic helices that resemble the P-loops of K⁺ channels (14, 27, 55). In the K⁺ channel, these pore-forming loops end in the filter residues glycine-tyrosine-glycine, which coordinate K⁺ by means of their backbones’ carbonyl oxygens (13). The P-loop-like helices of Trk/HKT transporters end in a single conserved glycine (48), and these glycines have been shown to determine K⁺ selectivity over Na⁺ of the transporters (34, 49). Thus, a Trk/HKT monomer with 8 TM domains and 4 P-loops is thought to have a similar pore architecture to a K⁺ channel tetramer with two TM domains and one P-loop per subunit. The Trk/HKT transporters are important for cellular K⁺ acquisition in microorganisms, since trk null mutant yeast or bacteria exhibit growth phenotypes on media containing low K⁺ concentrations (20, 46). The roles of the Trk/HKT transporters in plants are more diverse, including Na⁺ distribution (10, 33, 47), osmoregulation (32), and salt tolerance (39). So far, no HKT/Trk transporter has been described from the metazoa or protista.Trypanosoma brucei subsp. parasites comprise the causative agents of human and livestock trypanosomosis: sleeping sickness and nagana, respectively. The distribution of the parasites is restricted by that of their vector, the blood-sucking tsetse fly (Glossina spp.), to the so-called tsetse belt comprising 36 countries between the Sahara desert and the Kalahari (3, 8). African trypanosomes proliferate extracellularly in the blood, evading the mammalian immune response by antigenic variation. Untreated sleeping sickness is fatal. There is an urgent need for new and better drugs since the current ones, the arsenical melarsoprol in particular, suffer from severe side effects (31). In the mammalian bloodstream, the parasites encounter a rich and steady supply of nutrients, readily imported by specific permeases or endocytosed via receptors (7). Research on trypanosomal nutrient uptake has so far concentrated on transporters of organic substrates: nucleobases, nucleosides, sugars, and amino acids (4, 12, 26, 30, 35, 56). Little is known about how the parasites import inorganic nutrients. The malaria parasite Plasmodium falciparum possesses two putative K⁺ channel subunits with 6 TM domains and one P-loop (19, 52). Disruption of an orthologous gene in Plasmodium berghei strongly impaired the development of the malaria parasites in the mosquito (18). However, these putative channels have not yet been proven to be permeable to K⁺. The T. brucei genome (6) is annotated to contain two putative K⁺ channels; in addition, a putative ATPase has been identified resembling fungal Na⁺/K⁺ efflux ATPases (5, 45). None of these has been addressed experimentally. Here we present the identification and characterization of TbHKT1 (Tb10.70.2940), a Trk/HKT-type K⁺ transporter from Trypanosoma brucei and representative of a new clade of Trk/HKT genes from kinetoplastid parasites.     

Diversity Within the O-linked Protein Glycosylation Systems of Acinetobacter Species

The opportunistic human pathogen Acinetobacter baumannii is a concern to health care systems worldwide because of its persistence in clinical settings and the growing frequency of multiple drug resistant infections. To combat this threat, it is necessary to understand factors associated with disease and environmental persistence of A. baumannii. Recently, it was shown that a single biosynthetic pathway was responsible for the generation of capsule polysaccharide and O-linked protein glycosylation. Because of the requirement of these carbohydrates for virulence and the non-template driven nature of glycan biogenesis we investigated the composition, diversity, and properties of the Acinetobacter glycoproteome. Utilizing global and targeted mass spectrometry methods, we examined 15 strains and found extensive glycan diversity in the O-linked glycoproteome of Acinetobacter. Comparison of the 26 glycoproteins identified revealed that different A. baumannii strains target similar protein substrates, both in characteristics of the sites of O-glycosylation and protein identity. Surprisingly, glycan micro-heterogeneity was also observed within nearly all isolates examined demonstrating glycan heterogeneity is a widespread phenomena in Acinetobacter O-linked glycosylation. By comparing the 11 main glycoforms and over 20 alternative glycoforms characterized within the 15 strains, trends within the glycan utilized for O-linked glycosylation could be observed. These trends reveal Acinetobacter O-linked glycosylation favors short (three to five residue) glycans with limited branching containing negatively charged sugars such as GlcNAc3NAcA4OAc or legionaminic/pseudaminic acid derivatives. These observations suggest that although highly diverse, the capsule/O-linked glycan biosynthetic pathways generate glycans with similar characteristics across all A. baumannii.Acinetobacter baumannii is an emerging opportunistic pathogen of increasing significance to health care institutions worldwide (1–3). The growing number of identified multiple drug resistant (MDR)¹ strains (2–4), the ability of isolates to rapidly acquire resistance (3, 4), and the propensity of this agent to survive harsh environmental conditions (5) account for the increasing number of outbreaks in intensive care, burn, or high dependence health care units since the 1970s (2–5). The burden on the global health care system of MDR A. baumannii is further exacerbated by standard infection control measures often being insufficient to quell the spread of A. baumannii to high risk individuals and generally failing to remove A. baumannii from health care institutions (5). Because of these concerns, there is an urgent need to identify strategies to control A. baumannii as well as understand the mechanisms that enable its persistence in health care environments.Surface glycans have been identified as key virulence factors related to persistence and virulence within the clinical setting (6–8). Acinetobacter surface carbohydrates were first identified and studied in A. venetianus strain RAG-1, leading to the identification of a gene locus required for synthesis and export of the surface carbohydrates (9, 10). These carbohydrate synthesis loci are variable yet ubiquitous in A. baumannii (11, 12). Comparison of 12 known capsule structures from A. baumannii with the sequences of their carbohydrate synthesis loci has provided strong evidence that these loci are responsible for capsule synthesis with as many as 77 distinct serotypes identified by molecular serotyping (11). Because of the non-template driven nature of glycan synthesis, the identification and characterization of the glycans themselves are required to confirm the true diversity. This diversity has widespread implications for Acinetobacter biology as the resulting carbohydrate structures are not solely used for capsule biosynthesis but can be incorporated and utilized by other ubiquitous systems, such as O-linked protein glycosylation (13, 14).Although originally thought to be restricted to species such as Campylobacter jejuni (15, 16) and Neisseria meningitidis (17), bacterial protein glycosylation is now recognized as a common phenomenon within numerous pathogens and commensal bacteria (18, 19). Unlike eukaryotic glycosylation where robust and high-throughput technologies now exist to enrich (20–22) and characterize both the glycan and peptide component of glycopeptides (23–25), the diversity (glycan composition and linkage) within bacterial glycosylation systems makes few technologies broadly applicable to all bacterial glycoproteins. Because of this challenge a deeper understanding of the glycan diversity and substrates of glycosylation has been largely unachievable for the majority of known bacterial glycosylation systems. The recent implementation of selective glycopeptide enrichment methods (26, 27) and the use of multiple fragmentation approaches (28, 29) has facilitated identification of an increasing number of glycosylation substrates independent of prior knowledge of the glycan structure (30–33). These developments have facilitated the undertaking of comparative glycosylation studies, revealing glycosylation is widespread in diverse genera and far more diverse then initially thought. For example, Nothaft et al. were able to show N-linked glycosylation was widespread in the Campylobacter genus and that two broad groupings of the N-glycans existed (34).During the initial characterization of A. baumannii O-linked glycosylation the use of selective enrichment of glycopeptides followed by mass spectrometry analysis with multiple fragmentation technologies was found to be an effective means to identify multiple glycosylated substrates in the strain ATCC 17978 (14). Interestingly in this strain, the glycan utilized for protein modification was identical to a single subunit of the capsule (13) and the loss of either protein glycosylation or glycan synthesis lead to decreases in biofilm formation and virulence (13, 14). Because of the diversity in the capsule carbohydrate synthesis loci and the ubiquitous distribution of the PglL O-oligosaccharyltransferase required for protein glycosylation, we hypothesized that the glycan variability might be also extended to O-linked glycosylation. This diversity, although common in surface carbohydrates such as the lipopolysaccharide of numerous Gram-negative pathogens (35), has only recently been observed within bacterial proteins glycosylation system that are typically conserved within species (36) and loosely across genus (34, 37).In this study, we explored the diversity within the O-linked protein glycosylation systems of Acinetobacter species. Our analysis complements the recent in silico studies of A. baumannii showing extensive glycan diversity exists in the carbohydrate synthesis loci (11, 12). Employing global strategies for the analysis of glycosylation, we experimentally demonstrate that the variation in O-glycan structure extends beyond the genetic diversity predicted by the carbohydrate loci alone and targets proteins of similar properties and identity. Using this knowledge, we developed a targeted approach for the detection of protein glycosylation, enabling streamlined analysis of glycosylation within a range of genetic backgrounds. We determined that; O-linked glycosylation is widespread in clinically relevant Acinetobacter species; inter- and intra-strain heterogeneity exist within glycan structures; glycan diversity, although extensive results in the generation of glycans with similar properties and that the utilization of a single glycan for capsule and O-linked glycosylation is a general feature of A. baumannii but may not be a general characteristic of all Acinetobacter species such as A. baylyi.     

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?