Gas-phase rearrangements do not affect site localization reliability in phosphoproteomics data sets

Computational approaches to identify sites of phosphorylation

Due to their oftentimes ambiguous nature, phosphopeptide positional isomers can present challenges in bottom-up mass spectrometry-based workflows as search engine scores alone are often not enough to confidently distinguish them. Additional scoring algorithms can remedy this by providing confidence metrics in addition to these search results, reducing ambiguity. Here we describe challenges to interpreting phosphoproteomics data and review several different approaches to determine sites of phosphorylation for both data-dependent and data-independent acquisition-based workflows. Finally, we discuss open questions regarding neutral losses, gas-phase rearrangement, and false localization rate estimation experienced by both types of acquisition workflows and best practices for managing ambiguity in phosphosite determination.

Phosphate Migration versus the Loss of Phosphoric Acid in Protonated Phosphopeptides: A Computational Study

Residue-specific phosphorylation is a protein post-translational modification that regulates cellular functions. Experimental determination of the exact sites of protein phosphorylation provides an understanding of the signaling and processes at work for a given cellular state. Any experimental artifact that involves migration of the phosphate group during measurement is a concern, as the outcome can lead to erroneous conclusions that may confound studies on cellular signal transduction. Herein, we examine computationally the mechanism by which a phosphate group migrates from one serine residue to another serine in monoprotonated pentapeptides [BA-pSer-Gly-Ser-BB + H]+ → [BA-Ser-Gly-pSer-BB + H]+ (where BA and BB are different combinations of the three basic amino acids, histidine, lysine, and arginine). In addition to moving the phosphate group, the overall mechanism involves transferring a proton from the N-terminal amino acid, BA, to the C-terminal amino acid, BB. This is not a synchronous process, and there is a key high-energy intermediate, structure C, that is zwitterionic with both the basic amino acids protonated and the phosphate group attached to both serine residues and carrying a negative charge. The barriers to moving the phosphate group are calculated to be in the range of 219-274 kJ mol-1 at the B3LYP/6-31G(d) level. These barriers are systematically slightly lower and in good agreement with single-point energy calculations at both M06-2X/6-311++G(d,p) and MP2/6-31++G(d,p) levels. The competitive reaction, loss of phosphoric acid from the protonated pentapeptides, has a barrier in the range of 176-202 kJ mol-1 at the B3LYP/6-31G(d) level. Extension of the theory to M06-2X/6-311++G(d,p)//B3LYP/6-31G(d) and MP2/6-31++G(d,p)// B3LYP/6-31G(d) gives higher values for the loss of phosphoric acid, falling in the range of 196-226 kJ mol-1; these are comparable to the barriers against phosphate migration at the same levels of theory. For larger peptides His-pSer-(Gly)n-Ser-His, where n has values from 2 to 5, the barriers against the loss of phosphoric acid are higher than those against the phosphate group migration. This difference is most pronounced and significant when n = 4 and 5 (the differences are approximately 80 kJ mol-1 under the single-point energy calculations at the M06-2X and MP2 levels). Energy differences using two more recent functionals, M08-HX and MN15, on His-pSer-(Gly)n-Ser-His, where n = 1 and 5, are in good agreement with the M06-2X and MP2 calculations. These results provide the mechanistic rationale for phosphate migration versus other competing reactions in the gas phase under tandem mass spectrometry conditions.

Considerations for defining +80 Da mass shifts in mass spectrometry-based proteomics: phosphorylation and beyond

Post-translational modifications (PTMs) are ubiquitous and key to regulating protein function. Understanding the dynamics of individual PTMs and their biological roles requires robust characterisation. Mass spectrometry (MS) is the method of choice for the identification and quantification of protein modifications. This article focusses on the MS-based analysis of those covalent modifications that induce a mass shift of +80 Da, notably phosphorylation and sulfation, given the challenges associated with their discrimination and pinpointing the sites of modification on a polypeptide chain. Phosphorylation in particular is highly abundant, dynamic and can occur on numerous residues to invoke specific functions, hence robust characterisation is crucial to understanding biological relevance. Showcasing our work in the context of other developments in the field, we highlight approaches for enrichment and site localisation of phosphorylated (canonical and non-canonical) and sulfated peptides, as well as modification analysis in the context of intact proteins (top down proteomics) to explore combinatorial roles. Finally, we discuss the application of native ion-mobility MS to explore the effect of these PTMs on protein structure and ligand binding.

Post-Translational Modifications of Histones Are Versatile Regulators of Fungal Development and Secondary Metabolism

Chromatin structure is a major regulator of DNA-associated processes, such as transcription, DNA repair, and replication. Histone post-translational modifications, or PTMs, play a key role on chromatin dynamics. PTMs are involved in a wide range of biological processes in eukaryotes, including fungal species. Their deposition/removal and their underlying functions have been extensively investigated in yeasts but much less in other fungi. Nonetheless, the major role of histone PTMs in regulating primary and secondary metabolisms of filamentous fungi, including human and plant pathogens, has been pinpointed. In this review, an overview of major identified PTMs and their respective functions in fungi is provided, with a focus on filamentous fungi when knowledge is available. To date, most of these studies investigated histone acetylations and methylations, but the development of new methodologies and technologies increasingly allows the wider exploration of other PTMs, such as phosphorylation, ubiquitylation, sumoylation, and acylation. Considering the increasing number of known PTMs and the full range of their possible interactions, investigations of the subsequent Histone Code, i.e., the biological consequence of the combinatorial language of all histone PTMs, from a functional point of view, are exponentially complex. Better knowledge about histone PTMs would make it possible to efficiently fight plant or human contamination, avoid the production of toxic secondary metabolites, or optimize the industrial biosynthesis of certain beneficial compounds.

Mechanistic insights into intramolecular phosphate group transfer during collision induced dissociation of phosphopeptides

We report on the rearrangement chemistry of model phosphorylated peptides during collision-induced dissociation (CID), where intramolecular phosphate group transfers are observed from donor to acceptor residues. Such "scrambling" could result in inaccurate modification localization, potentially leading to misidentifications. Systematic studies presented herein provide mechanistic insights for the unusually high phosphate group rearrangements presented some time ago by Reid and coworkers (Proteomics 2013, 13 [6], 964-973). It is postulated here that a basic residue like histidine can play a key role in mediating the phosphate group transfer by deprotonating the serine acceptor site. The proposed mechanism is consistent with the observation that fast collisional activation by collision-cell CID and higher-energy collisional dissociation (HCD) can shut down rearrangement chemistry. Additionally, the rearrangement chemistry is highly dependent on the charge state of the peptide, mirroring previous studies that less rearrangement is observed under mobile proton conditions.

The Pivotal Role of Protein Phosphorylation in the Control of Yeast Central Metabolism

Protein phosphorylation is the most frequent eukaryotic post-translational modification and can act as either a molecular switch or rheostat for protein functions. The deliberate manipulation of protein phosphorylation has great potential for regulating specific protein functions with surgical precision, rather than the gross effects gained by the over/underexpression or complete deletion of a protein-encoding gene. In order to assess the impact of phosphorylation on central metabolism, and thus its potential for biotechnological and medical exploitation, a compendium of highly confident protein phosphorylation sites (p-sites) for the model organism Saccharomyces cerevisiae has been analyzed together with two more datasets from the fungal pathogen Candida albicans. Our analysis highlights the global properties of the regulation of yeast central metabolism by protein phosphorylation, where almost half of the enzymes involved are subject to this sort of post-translational modification. These phosphorylated enzymes, compared to the nonphosphorylated ones, are more abundant, regulate more reactions, have more protein–protein interactions, and a higher fraction of them are ubiquitinated. The p-sites of metabolic enzymes are also more conserved than the background p-sites, and hundreds of them have the potential for regulating metabolite production. All this integrated information has allowed us to prioritize thousands of p-sites in terms of their potential phenotypic impact. This multi-source compendium should enable the design of future high-throughput (HTP) mutation studies to identify key molecular switches/rheostats for the manipulation of not only the metabolism of yeast, but also that of many other biotechnologically and medically important fungi and eukaryotes.

Neutral Loss Is a Very Common Occurrence in Phosphotyrosine-Containing Peptides Labeled with Isobaric Tags

While developing a multiplexed phosphotyrosine peptide quantification assay, an unexpected observation was made: significant neutral loss from phosphotyrosine (pY) containing peptides. Using a 2000-member peptide library, we sought to systematically investigate this observation by comparing unlabeled peptides with the two highest-plex isobaric tags (iTRAQ8 and TMT10) across CID, HCD, and ETD fragmentation using high resolution high mass accuracy Orbitrap instrumentation. We found pY peptide neutral loss behavior was consistent with reduced proton mobility, and does not occur during ETD. The site of protonation at the peptide N-terminus changes from a primary to a tertiary amine as a result of TMT labeling which would increase the gas phase basicity and reduce proton mobility at this site. This change in fragmentation behavior has implications during instrument method development and interpretation of MS/MS spectra, and therefore ensuing follow-up studies. We show how sites not localized to tyrosine by search and site localization algorithms can be confidently reassigned to tyrosine using neutral loss and phosphotyrosine immonium ions. We believe these findings will be of general interest to those studying pY signal transduction using isobaric tags.

Gas-Phase Rearrangement in Lysine Phosphorylated Peptides During Electron-Transfer Dissociation Tandem Mass Spectrometry

Tandem mass spectrometry (MS/MS) strategies coupled with collision-induced dissociation (CID) or radical-driven fragmentation techniques such as electron-capture dissociation (ECD) or electron-transfer dissociation (ETD) have been successfully used for comprehensive phosphoproteome analysis. However, the unambiguous characterization of the phosphorylation site remains a significant challenge due to phosphate-related neutral losses and gas-phase rearrangements, which have been observed during CID. In particular, for the analysis of labile N-phosphorylated peptides, ECD and ETD are emerging as a complementary method. In contrast to CID, the phosphorylation site of histidine, arginine, and lysine phosphorylated peptides can be characterized by ETD. Here, we present a study on the application of ETD for analysis of phospholysine (pLys) peptides. We show that, depending on the charge state of the precursor ion as well as the presence of basic amino acid side chains, phosphate transfer reactions during the ETD process can be observed leading to ambiguous fragment ion spectra. Basically, pLys is stable under ETD conditions allowing an unambiguous assignment of the site of phosphorylation, but some factors/parameters have to be considered to avoid gas-phase rearrangement which would lead to false positive results in phosphoproteomic studies.

Computational phosphoproteomics: from identification to localization

Analysis of the phosphoproteome by MS has become a key technology for the characterization of dynamic regulatory processes in the cell, since kinase and phosphatase action underlie many major biological functions. However, the addition of a phosphate group to a suitable side chain often confounds informatic analysis by generating product ion spectra that are more difficult to interpret (and consequently identify) relative to unmodified peptides. Collectively, these challenges have motivated bioinformaticians to create novel software tools and pipelines to assist in the identification of phosphopeptides in proteomic mixtures, and help pinpoint or "localize" the most likely site of modification in cases where there is ambiguity. Here we review the challenges to be met and the informatics solutions available to address them for phosphoproteomic analysis, as well as highlighting the difficulties associated with using them and the implications for data standards.

The human embryonic stem cell proteome revealed by multidimensional fractionation followed by tandem mass spectrometry

Human embryonic stem cells (hESCs) have received considerable attention due to their therapeutic potential and usefulness in understanding early development and cell fate commitment. In order to appreciate the unique properties of these pluripotent, self-renewing cells, we have performed an in-depth multidimensional fractionation followed by LC-MS/MS analysis of the hESCs harvested from defined media to elucidate expressed, phosphorylated, O-linked β-N-acetylglucosamine (O-GlcNAc) modified, and secreted proteins. From the triplicate analysis, we were able to assign more than 3000 proteins with less than 1% false-discovery rate. This analysis also allowed us to identify nearly 500 phosphorylation sites and 68 sites of O-GlcNAc modification with the same high confidence. Investigation of the phosphorylation sites allowed us to deduce the set of kinases that are likely active in these cells. We also identified more than 100 secreted proteins of hESCs that likely play a role in extracellular matrix formation and remodeling, as well as autocrine signaling for self-renewal and maintenance of the undifferentiated state. Finally, by performing in-depth analysis in triplicate, spectral counts were obtained for these proteins and posttranslationally modified peptides, which will allow us to perform relative quantitative analysis between these cells and any derived cell type in the future.

Fragmentation of integral membrane proteins in the gas phase

Integral membrane proteins (IMPs) are of great biophysical and clinical interest because of the key role they play in many cellular processes. Here, a comprehensive top down study of 152 IMPs and 277 soluble proteins from human H1299 cells including 11 087 fragments obtained from collisionally activated dissociation (CAD), 6452 from higher-energy collisional dissociation (HCD), and 2981 from electron transfer dissociation (ETD) shows their great utility and complementarity for the identification and characterization of IMPs. A central finding is that ETD is ∼2-fold more likely to cleave in soluble regions than threshold fragmentation methods, whereas the reverse is observed in transmembrane domains with an observed ∼4-fold bias toward CAD and HCD. The location of charges just prior to dissociation is consistent with this directed fragmentation: protons remain localized on basic residues during ETD but easily mobilize along the backbone during collisional activation. The fragmentation driven by these protons, which is most often observed in transmembrane domains, both is of higher yield and occurs over a greater number of backbone cleavage sites. Further, while threshold dissociation events in transmembrane domains are on average 10.1 (CAD) and 9.2 (HCD) residues distant from the nearest charge site (R, K, H, N-terminus), fragmentation is strongly influenced by the N- or C-terminal position relative to that site: the ratio of observed b- to y-fragments is ∼1:3 if the cleavage occurs >7 residues N-terminal and ∼3:1 if it occurs >7 residues C-terminal to the nearest basic site. Threshold dissociation products driven by a mobilized proton appear to be strongly dependent on not only relative position of a charge site but also N- or C-terminal directionality of proton movement.

Confident and sensitive phosphoproteomics using combinations of collision induced dissociation and electron transfer dissociation

We present a workflow using an ETD-optimised version of Mascot Percolator and a modified version of SLoMo (turbo-SLoMo) for analysis of phosphoproteomic data. We have benchmarked this against several database searching algorithms and phosphorylation site localisation tools and show that it offers highly sensitive and confident phosphopeptide identification and site assignment with PSM-level statistics, enabling rigorous comparison of data acquisition methods. We analysed the Plasmodium falciparum schizont phosphoproteome using for the first time, a data-dependent neutral loss-triggered-ETD (DDNL) strategy and a conventional decision-tree method. At a posterior error probability threshold of 0.01, similar numbers of PSMs were identified using both methods with a 73% overlap in phosphopeptide identifications. The false discovery rate associated with spectral pairs where DDNL CID/ETD identified the same phosphopeptide was < 1%. 72% of phosphorylation site assignments using turbo-SLoMo without any score filtering, were identical and 99.8% of these cases are associated with a false localisation rate of < 5%. We show that DDNL acquisition is a useful approach for phosphoproteomics and results in an increased confidence in phosphopeptide identification without compromising sensitivity or duty cycle. Furthermore, the combination of Mascot Percolator and turbo-SLoMo represents a robust workflow for phosphoproteomic data analysis using CID and ETD fragmentation.

Biological significance

Protein phosphorylation is a ubiquitous post-translational modification that regulates protein function. Mass spectrometry-based approaches have revolutionised its analysis on a large-scale but phosphorylation sites are often identified by single phosphopeptides and therefore require more rigorous data analysis to unsure that sites are identified with high confidence for follow-up experiments to investigate their biological significance. The coverage and confidence of phosphoproteomic experiments can be enhanced by the use of multiple complementary fragmentation methods. Here we have benchmarked a data analysis pipeline for analysis of phosphoproteomic data generated using CID and ETD fragmentation and used it to demonstrate the utility of a data-dependent neutral loss triggered ETD fragmentation strategy for high confidence phosphopeptide identification and phosphorylation site localisation.

•

We report and benchmark a data analysis pipeline for phosphoproteomic data analysis.

•

Combined use of Mascot Percolator and turbo-SLoMo to compare fragmentation methods

•

CID and ETD fragmentation for phosphorylation site identification

•

Demonstrate the utility of data-dependent neutral loss triggered ETD fragmentation

•

High confidence of phosphoproteomic analysis using ETD/CID spectral pairs

Quantification of competing H3PO4 versus HPO3 + H2O neutral losses from regioselective 18O-labeled phosphopeptides

Abundant neutral losses of 98 Da are often observed upon ion trap CID-MS/MS of protonated phosphopeptide ions. Two competing fragmentation pathways are involved in this process, namely, the direct loss of H3PO4 from the phosphorylated residue and the combined losses of HPO3 and H2O from the phosphorylation site and from an additional site within the peptide, respectively. These competing pathways produce product ions with different structures but the same m/z values, potentially limiting the utility of CID-MS(3) for phosphorylation site localization. To quantify the relative contributions of these pathways and to determine the conditions under which each pathway predominates, we have examined the ion trap CID-MS/MS fragmentation of a series of regioselective (18)O-phosphate ester labeled phosphopeptides prepared using novel solution-phase amino acid synthesis and solid-phase peptide synthesis methodologies. By comparing the intensity of the -100 Da (-H3PO3 (18)O) versus -98 Da (-[HPO3 + H2O]) neutral loss product ions formed upon MS/MS, quantification of the two pathways was achieved. Factors that affect the extent of formation of the competing neutral losses were investigated, with the combined loss pathway predominantly occurring under conditions of limited proton mobility, and with increased combined losses observed for phosphothreonine compared with phosphoserine-containing peptides. The combined loss pathway was found to be less dominant under ion activation conditions associated with HCD-MS/MS. Finally, the contribution of carboxylic acid functional groups and backbone amide bonds to the water loss in the combined loss fragmentation pathway was determined via methyl esterification and by examination of a phosphopeptide lacking side-chain hydroxyl groups.

Examining factors that influence erroneous phosphorylation site localization via competing fragmentation and rearrangement reactions during ion trap CID-MS/MS and -MS(3.)

Factors influencing the magnitude of competing fragmentation and intramolecular phosphate group rearrangement reactions during CID-MS/MS and CID-MS(3) of protonated phosphopeptide ions in ion trap mass spectrometers, and their effect on phosphorylation site localization using automated search algorithms, have been examined by systematically varying the peptide composition, the identity, number, and position of the phosphorylated "donor" and nonphosphorylated "acceptor" residues, and the proton mobility of the precursor ion charge states for a synthetic phosphopeptide library. CID-MS(3) of product ions formed via combined neutral losses of HPO3 and H2 O, rather than direct loss of H3 PO4 from phosphotyrosine containing peptides yielded incorrect phosphorylation site assignments, while correct phosphorylation site assignments for phosphothreonine and phosphoserine containing peptides were highly dependant on the relative abundance of these competing fragmentation pathways. Abundant phosphate group rearrangement product ions were observed from CID-MS/MS of multiply protonated phosphopeptide ions, with increased rearrangement under nonmobile or partially mobile protonation conditions, and as a function of the identity and number of the donor and acceptor residues. A clear inverse trend was observed between the amplitude of these rearrangement reactions and the confidence for phosphorylation site localization, and rearrangement played a contributing role in erroneous phosphorylation site assignment for several peptides.

Quantitative global phosphoproteomics of human umbilical vein endothelial cells after activation of the Rap signaling pathway

The small GTPase Rap1 is required for proper cell-cell junction formation and also plays a key role in mediating cAMP-induced tightening of adherens junctions and subsequent increased barrier function of endothelial cells. To further study how Rap1 controls barrier function, we performed quantitative global phosphoproteomics in human umbilical vein endothelial cells (HUVECs) prior to and after Rap1 activation by the Epac-selective cAMP analog 8-pCPT-2'-O-Me-cAMP-AM (007-AM). Tryptic digests were labeled using stable isotope dimethyl labeling, enriched with phosphopeptides by strong cation exchange (SCX), followed by titanium(iv) immobilized metal affinity chromatography (Ti(4+)-IMAC) and analyzed by high resolution mass spectrometry. We identified 19 859 unique phosphopeptides containing 17 278 unique phosphosites on 4594 phosphoproteins, providing the largest HUVEC phosphoproteome to date. Of all identified phosphosites, 220 (∼1%) were more than 1.5-fold up- or downregulated upon Rap activation, in two independent experiments. Compatible with the function of Rap1, these alterations were found predominantly in proteins regulating the actin cytoskeleton, cell-cell junctions and cell adhesion.

Moving away from the reference genome: evaluating a peptide sequencing tagging approach for single amino acid polymorphism identifications in the genus Populus

The genetic diversity across natural populations of the model organism, Populus, is extensive, containing a single nucleotide polymorphism roughly every 200 base pairs. When deviations from the reference genome occur in coding regions, they can impact protein sequences. Rather than relying on a static reference database to profile protein expression, we employed a peptide sequence tagging (PST) approach capable of decoding the plasticity of the Populus proteome. Using shotgun proteomics data from two genotypes of P. trichocarpa, a tag-based approach enabled the detection of 6653 unexpected sequence variants. Through manual validation, our study investigated how the most abundant chemical modification (methionine oxidation) could masquerade as a sequence variant (Ala→Ser) when few site-determining ions existed. In fact, precise localization of an oxidation site for peptides with more than one potential placement was indeterminate for 70% of the MS/MS spectra. We demonstrate that additional fragment ions made available by high energy collisional dissociation enhances the robustness of the peptide sequence tagging approach (81% of oxidation events could be exclusively localized to a methionine). We are confident that augmenting fragmentation processes for a PST approach will further improve the identification of single amino acid polymorphism in Populus and potentially other species as well.

CTR1 phosphorylates the central regulator EIN2 to control ethylene hormone signaling from the ER membrane to the nucleus in Arabidopsis

The gaseous phytohormone ethylene C(2)H(4) mediates numerous aspects of growth and development. Genetic analysis has identified a number of critical elements in ethylene signaling, but how these elements interact biochemically to transduce the signal from the ethylene receptor complex at the endoplasmic reticulum (ER) membrane to transcription factors in the nucleus is unknown. To close this gap in our understanding of the ethylene signaling pathway, the challenge has been to identify the target of the CONSTITUTIVE TRIPLE RESPONSE1 (CTR1) Raf-like protein kinase, as well as the molecular events surrounding ETHYLENE-INSENSITIVE2 (EIN2), an ER membrane-localized Nramp homolog that positively regulates ethylene responses. Here we demonstrate that CTR1 interacts with and directly phosphorylates the cytosolic C-terminal domain of EIN2. Mutations that block the EIN2 phosphorylation sites result in constitutive nuclear localization of the EIN2 C terminus, concomitant with constitutive activation of ethylene responses in Arabidopsis. Our results suggest that phosphorylation of EIN2 by CTR1 prevents EIN2 from signaling in the absence of ethylene, whereas inhibition of CTR1 upon ethylene perception is a signal for cleavage and nuclear localization of the EIN2 C terminus, allowing the ethylene signal to reach the downstream transcription factors. These findings significantly advance our understanding of the mechanisms underlying ethylene signal transduction.

Automated and high confidence protein phosphorylation site localization using complementary collision-activated dissociation and electron transfer dissociation tandem mass spectrometry

Reversible protein phosphorylation plays a critical role in cell signaling and is responsible for the regulation of many biological processes in most living organisms. The low stoichiometry of protein phosphorylation requires sensitive analysis by tandem mass spectrometry. However, incomplete peptide fragmentation and the loss of labile phosphate groups complicate identification of the site of the phosphate motif. Here, we have implemented and evaluated a novel approach for phospho-site localization by the combined use of peptide tandem mass spectrometry data obtained using both collision-activated dissociation and electron transfer dissociation, an approach termed the Cscore. The scoring algorithm used in the Cscore was adapted from the widely used Ascore method. The analytical benefit of integrating the product ion information of both ETD and CAD data are evident by increased confidence in phospho-site localization and the number of assigned phospho-sites at a fixed false-localization rate. The average calculated Cscore from a large data set (>7000 phosphopeptide MS/MS spectra) was ∼32 compared to ∼23 and ∼17 for the Ascore using collision-activated dissociation (CAD) or electron transfer dissociation (ETD), respectively. Compared with the Ascore using either CAD or ETD, the Cscore identified up to 88% more phosphorylation sites. Using a phosphopeptide library revealed that the score threshold for obtaining a false-localization rate of 0.5% was lower for the Cscore than either the Ascore (CAD) or the Ascore (ETD).

Peptide identification by tandem mass spectrometry with alternate fragmentation modes

The high-throughput nature of proteomics mass spectrometry is enabled by a productive combination of data acquisition protocols and the computational tools used to interpret the resulting spectra. One of the key components in mainstream protocols is the generation of tandem mass (MS/MS) spectra by peptide fragmentation using collision induced dissociation, the approach currently used in the large majority of proteomics experiments to routinely identify hundreds to thousands of proteins from single mass spectrometry runs. Complementary to these, alternative peptide fragmentation methods such as electron capture/transfer dissociation and higher-energy collision dissociation have consistently achieved significant improvements in the identification of certain classes of peptides, proteins, and post-translational modifications. Recognizing these advantages, mass spectrometry instruments now conveniently support fine-tuned methods that automatically alternate between peptide fragmentation modes for either different types of peptides or for acquisition of multiple MS/MS spectra from each peptide. But although these developments have the potential to substantially improve peptide identification, their routine application requires corresponding adjustments to the software tools and procedures used for automated downstream processing. This review discusses the computational implications of alternative and alternate modes of MS/MS peptide fragmentation and addresses some practical aspects of using such protocols for identification of peptides and post-translational modifications.

Characterization of O-sulfopeptides by negative ion mode tandem mass spectrometry: superior performance of negative ion electron capture dissociation

Positive ion mode collision-activated dissociation tandem mass spectrometry (CAD MS/MS) of O-sulfopeptides precludes determination of sulfonated sites due to facile proton-driven loss of the highly labile sulfonate groups. A previously proposed method for localizing peptide and protein O-sulfonation involves derivatization of nonsulfonated tyrosines followed by positive ion CAD MS/MS of the corresponding modified sulfopeptides for diagnostic sulfonate loss. This indirect method relies upon specific and complete derivatization of nonsulfonated tyrosines. Alternative MS/MS activation methods, including positive ion metastable atom-activated dissociation (MAD) and metal-assisted electron transfer dissociation (ETD) or electron capture dissociation (ECD) provide varying degrees of sulfonate retention. Sulfonate retention has also been reported following negative ion MAD and electron detachment dissociation (EDD), which also operates in negative ion mode in which sulfonate groups are less labile than in positive ion mode. However, an MS/MS activation technique that can effectively preserve sulfonate groups while providing extensive backbone fragmentation (translating to sequence information, including sulfonated sites) with little to no noninformative small molecule neutral loss has not previously been realized. Here, we report that negative ion CAD, EDD, and negative ETD (NETD) result in sulfonate retention mainly at higher charge states with varying degrees of fragmentation efficiency and sequence coverage. Similar to previous observations from CAD of sulfonated glycosaminoglycan anions, higher charge states translate to a higher probability of deprotonation at the sulfonate groups thus yielding charge-localized fragmentation without loss of the sulfonate groups. However, consequently, higher sulfonate retention comes at the price of lower sequence coverage in negative ion CAD. Fragmentation efficiency/sequence coverage averaged 19/6% and 33/20% in EDD and NETD, respectively, both of which are only applicable to multiply-charged anions. In contrast, the recently introduced negative ion ECD showed an average fragmentation efficiency of 69% and an average sequence coverage of 82% with complete sulfonate retention from singly- and doubly-deprotonated sulfopeptide anions.

Comparison of CID, ETD and metastable atom-activated dissociation (MAD) of doubly and triply charged phosphorylated tau peptides

The fragmentation behavior of the 2+ and 3+ charge states of eleven different phosphorylated tau peptides was studied using collision-induced dissociation (CID), electron transfer dissociation (ETD) and metastable atom-activated dissociation (MAD). The synthetic peptides studied contain up to two known phosphorylation sites on serine or threonine residues, at least two basic residues, and between four and eight potential sites of phosphorylation. CID produced mainly b-/y-type ions with abundant neutral losses of the phosphorylation modification. ETD produced c-/z-type ions in highest abundance but also showed numerous y-type ions at a frequency about 50% that of the z-type ions. The major peaks observed in the ETD spectra correspond to the charge-reduced product ions and small neutral losses from the charge-reduced peaks. ETD of the 2+ charge state of each peptide generally produced fewer backbone cleavages than the 3+ charge state, consistent with previous reports. Regardless of charge state, MAD achieved more extensive backbone cleavage than CID or ETD, while retaining the modification(s) in most cases. In all but one case, unambiguous modification site determination was achieved with MAD. MAD produced 15-20% better sequence coverage than CID and ETD for both the 2+ and 3+ charge states and very different fragmentation products indicating that the mechanism of fragmentation in MAD is unique and complementary to CID and ETD.

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

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

193 nm ultraviolet photodissociation of imidazolinylated Lys-N peptides for de novo sequencing

The goal of many MS/MS de novo sequencing strategies is to generate a single product ion series that can be used to determine the precursor ion sequence. Most methods fall short of achieving such simplified spectra, and the presence of additional ion series impede peptide identification. The present study aims to solve the problem of confounding ion series by enhancing the formation of "golden" sets of a, b, and c ions for sequencing. Taking advantage of the characteristic mass differences between the golden ions allows N-terminal fragments to be readily identified while other ion series are excluded. By combining the use of Lys-N, an alternate protease, to produce peptides with lysine residues at each N-terminus with subsequent imidazolinylation of the ε-amino group of each lysine, peptides with highly basic sites localized at each N-terminus are generated. Subsequent MS/MS analysis by using 193 nm ultraviolet photodissociation (UVPD) results in enhanced formation of the diagnostic golden pairs and golden triplets that are ideal for de novo sequencing.

Combined electron transfer dissociation-collision-induced dissociation fragmentation in the mass spectrometric distinction of leucine, isoleucine, and hydroxyproline residues in Peptide natural products

Distinctions between isobaric residues have been a major challenge in mass spectrometric peptide sequencing. Here, we propose a methodology for distinction among isobaric leucine, isoleucine, and hydroxyproline, a commonly found post-translationally modified amino acid with a nominal mass of 113 Da, through a combined electron transfer dissociation-collision-induced dissociation approach. While the absence of c and z(•) ions, corresponding to the Yyy-Xxx (Xxx = Leu, Ile, or Hyp) segment, is indicative of the presence of hydroxyproline, loss of isopropyl (Δm = 43 Da) or ethyl radicals (Δm = 29 Da), through collisional activation of z radical ions, are characteristic of leucine or isoleucine, respectively. Radical migration processes permit distinctions even in cases where the specific z(•) ions, corresponding to the Yyy-Leu or -Ile segments, are absent or of low intensity. This tandem mass spectrometric (MS(n)) method has been successfully implemented in a liquid chromatography-MS(n) platform to determine the identity of 23 different isobaric residues from a mixture of five different peptides. The approach is convenient for distinction of isobaric residues from any crude peptide mixture, typically encountered in natural peptide libraries or proteomic analysis.

Detection of two minor phosphorylation sites for bovine κ-casein macropeptide by reversed-phase liquid chromatography-tandem mass spectrometry

This work addresses the characterization of phosphopeptides in bovine κ-casein macropeptide by reversed-phase liquid chromatography-electrospray ionization-tandem mass spectrometry (RPLC-ESI-MS(2)). Two different mass spectrometers, equipped with an ion trap (IT) or a quadrupole time-of-flight (Q-TOF) analyzer, were used to perform an accurate phosphorylation site assignment. A total of 8 phosphopeptides from 26 identified peptides were characterized. MS(2) spectra of phosphopeptides were dominated by the neutral loss of a phosphoric acid molecule (H(3)PO(4)) and sufficient informative fragment ions resulting from peptide backbone cleavages enabling the elucidation of the phosphopeptide sequence. A higher number of sequence informative b and y ions were detected using a Q-TOF instead of an IT analyzer. In addition to the well-established phosphorylation sites at Ser(149) and Ser(127), this study also revealed the presence of two minor phosphorylation sites at Thr(145) and Ser(166). These findings indicate that RPLC-ESI-MS(2) on a Q-TOF analyzer is a useful technique for identifying low-abundance phosphorylation sites in caseins.

Phosphorylation of human eukaryotic initiation factor 2γ: novel site identification and targeted PKC involvement

Eukaryotic translation requires a suite of proteins known as eukaryotic initiation factors (eIFs). These molecular effectors oversee the highly regulated initiation phase of translation. Essential to eukaryotic translation initiation is the protein eIF2, a heterotrimeric protein composed of the individually distinct subunits eIF2α, eIF2β, and eIF2γ. The ternary complex, formed when eIF2 binds to GTP and Met-tRNA(i), is responsible for shuttling Met-tRNA(i) onto the awaiting 40S ribosome. As a necessary component for translation initiation, much attention has been given to the phosphorylation of eIF2α. Despite several previous investigations into eIF2 phosphorylation, most have centered on α- or β-subunit phosphorylation and little is known regarding γ-subunit phosphorylation. Herein, we report eight sites of phosphorylation on the largest eIF2 subunit with seven novel phosphosite identifications via high resolution mass spectrometry. Of the eight sites identified, three are located in either the switch regions or nucleotide binding pocket domain. In addition, we have identified a possible kinase of eIF2, protein kinase C (PKC), which is capable of phosphorylating threonine 66 (thr-66) on the intact heterotrimer. These findings may shed new light on the regulation of ternary complex formation and alternate molecular effectors involved in this process prior to 80S ribosome formation and subsequent translation elongation and termination.

Modelling of the gas-phase phosphate group loss and rearrangement in phosphorylated peptides

The gas-phase dissociation of phosphorylated peptides was modelled using a combination of quantum mechanics and the Rice-Ramsperger-Kassel-Marcus theory. Potential energy surfaces and unimolecular reaction rates for several low-energy fragmentation and rearrangement pathways were estimated, and a general mechanism was proposed. The neutral loss of the phosphoric acid was mainly an outcome of the intramolecular nucleophilic substitution mechanism. The mechanism involves a nucleophilic attack of the phosphorylated amino acid N-terminal carbonyl oxygen on β-carbon, yielding a cyclic five-membered oxazoline product ion. Regardless of the proton mobility, the pathway was charge directed either by a mobile proton or by a positively charged side chain of some basic residue. Although the mechanistic aspects of the phosphate loss are not influenced by the proton mobility environment, it does affect ion abundances. Results suggest that under the mobile proton environment, the interplay between phosphoric acid neutral loss product ion and backbone cleavage fragments should occur. On the other hand, when proton mobility is limited, neutral loss product ion may predominate. The fragmentation dynamics of phosphoserine versus phosphothreonine containing peptides suggests that H(3)PO(4) neutral loss from phosphothreonine containing peptides is less abundant than that from their phosphoserine containing analogs. During the low-energy CID of phosphorylated peptides in the millisecond time range, typical for ion trap instruments, a phosphate group rearrangement may happen, resulting in an interchange between the phosphorylated and the hydroxylated residues. Unimolecular dissociation rate constants imply the low abundance of such scrambled product ions.

Proteomic responses in Arabidopsis thaliana seedlings treated with ethylene

Ethylene (ET) is a volatile hormone that modulates fruit ripening, plant growth, development and stress responses. Key components of the ET-signaling pathway identified by genetic dissection in Arabidopsis thaliana include five ET receptors, the negative regulator CTR1 and the positive regulator EIN2, all of which localize to the endoplasmic reticulum. Mechanisms of signaling among these proteins are still unresolved and targets of ET responses are not fully known. So, we used mass spectrometry to identify proteins in microsomal membrane preparations from etiolated A. thaliana seedlings maintained in ambient air or treated with ET for 3 h. We compared 3814 proteins from ET-exposed seedlings and controls and identified 304 proteins with significant accumulation changes. The proteins with increased accumulation were involved in ET biosynthesis, cell morphogenesis, oxidative stress and vesicle secretion while those with decreased accumulation were ribosomal proteins and proteins positively regulated by brassinosteroid, another hormone involved in cell elongation. Several proteins, including EIN2, appeared to be differentially phosphorylated upon ET treatment, which suggests that the activity or stability of these proteins may be controlled by phosphorylation. TUA3, a component of microtubules that contributes to cellular morphological change, exhibited both increased accumulation and differential phosphorylation upon ET treatment. To verify the role of TUA3 in the ET response, tua3 mutants were evaluated. Mutant seedlings had altered ET-associated growth movements. The data indicate that ET perception leads to rapid proteomic change and that these changes are an important part of signaling and development. The data serve as a foundation for exploring ET signaling through systems biology.

Enhancing the identification of phosphopeptides from putative basophilic kinase substrates using Ti (IV) based IMAC enrichment

Metal and metal oxide chelating-based phosphopeptide enrichment technologies provide powerful tools for the in-depth profiling of phosphoproteomes. One weakness inherent to current enrichment strategies is poor binding of phosphopeptides containing multiple basic residues. The problem is exacerbated when strong cation exchange (SCX) is used for pre-fractionation, as under low pH SCX conditions phosphorylated peptides with multiple basic residues elute with the bulk of the tryptic digest and therefore require more stringent enrichment. Here, we report a systematic evaluation of the characteristics of a novel phosphopeptide enrichment approach based on a combination of low pH SCX and Ti(4+)-immobilized metal ion affinity chromatography (IMAC) comparing it one-to-one with the well established low pH SCX-TiO(2) enrichment method. We also examined the effect of 1,1,1,3,3,3-hexafluoroisopropanol (HFP), trifluoroacetic acid (TFA), or 2,5-dihydroxybenzoic acid (DHB) in the loading buffer, as it has been hypothesized that high levels of TFA and the perfluorinated solvent HFP improve the enrichment of phosphopeptides containing multiple basic residues. We found that Ti(4+)-IMAC in combination with TFA in the loading buffer, outperformed all other methods tested, enabling the identification of around 5000 unique phosphopeptides containing multiple basic residues from 400 μg of a HeLa cell lysate digest. In comparison, ∼ 2000 unique phosphopeptides could be identified by Ti(4+)-IMAC with HFP and close to 3000 by TiO(2). We confirmed, by motif analysis, the basic phosphopeptides enrich the number of putative basophilic kinases substrates. In addition, we performed an experiment using the SCX/Ti(4+)-IMAC methodology alongside the use of collision-induced dissociation (CID), higher energy collision induced dissociation (HCD) and electron transfer dissociation with supplementary activation (ETD) on considerably more complex sample, consisting of a total of 400 μg of triple dimethyl labeled MCF-7 digest. This analysis led to the identification of over 9,000 unique phosphorylation sites. The use of three peptide activation methods confirmed that ETD is best capable of sequencing multiply charged peptides. Collectively, our data show that the combination of SCX and Ti(4+)-IMAC is particularly advantageous for phosphopeptides with multiple basic residues.

Quantitative phosphoproteomics strategies for understanding protein kinase-mediated signal transduction pathways

Protein phosphorylation is a central regulatory mechanism of cell signaling pathways. This highly controlled biochemical process is involved in most cellular functions, and defects in protein kinases and phosphatases have been implicated in many diseases, highlighting the importance of understanding phosphorylation-mediated signaling networks. However, phosphorylation is a transient modification, and phosphorylated proteins are often less abundant. Therefore, the large-scale identification and quantification of phosphoproteins and their phosphorylation sites under different conditions are one of the most interesting and challenging tasks in the field of proteomics. Both 2D gel electrophoresis and liquid chromatography-tandem mass spectrometry serve as key phosphoproteomic technologies in combination with prefractionation, such as enrichment of phosphorylated proteins/peptides. Recently, new possibilities for quantitative phosphoproteomic analysis have been offered by technical advances in sample preparation, enrichment, separation, instrumentation, quantification and informatics. In this article, we present an overview of several strategies for quantitative phosphoproteomics and discuss how phosphoproteomic analysis can help to elucidate signaling pathways that regulate various cellular processes.

Catch me if you can: mass spectrometry-based phosphoproteomics and quantification strategies

Phosphorylation of proteins is one of the most prominent PTMs and for instance a key regulator of signal transduction. In order to improve our understanding of cellular phosphorylation events, considerable effort has been devoted to improving the analysis of phosphorylation by MS-based proteomics. Different enrichment strategies for phosphorylated peptides/proteins, such as immunoaffinity chromatography (IMAC) or titanium dioxide, have been established and constantly optimized for subsequent MS analysis. Concurrently, specific MS techniques were developed for more confident identification and phosphorylation site localization. In addition, more attention is paid to the LC-MS instrumentation to avoid premature loss of phosphorylated peptides within the analytical system. Despite major advances in all of these fields, the analysis of phosphopeptides still remains far from being routine in proteomics. However, to reveal cellular regulation by phosphorylation events, not only qualitative information about the phosphorylation status of proteins but also, in particular, quantitative information about distinct changes in phosphorylation patterns upon specific stimulation is mandatory. Thus, yielded insights are of outstanding importance for the emerging field of systems biology. In this review, we will give an insight into the historical development of phosphoproteome analysis and discuss its recent progress particularly regarding phosphopeptide quantification and assessment of phosphorylation stoichiometry.

Enhanced characterization of singly protonated phosphopeptide ions by femtosecond laser-induced ionization/dissociation tandem mass spectrometry (fs-LID-MS/MS)

To develop an improved understanding of the regulatory role that post-translational modifications (PTMs) involving phosphorylation play in the maintenance of normal cellular function, tandem mass spectrometry (MS/MS) strategies coupled with ion activation techniques such as collision-induced dissociation (CID) and electron-transfer dissociation (ETD) are typically employed to identify the presence and site-specific locations of the phosphate moieties within a given phosphoprotein of interest. However, the ability of these techniques to obtain sufficient structural information for unambiguous phosphopeptide identification and characterization is highly dependent on the ion activation method employed and the properties of the precursor ion that is subjected to dissociation. Herein, we describe the application of a recently developed alternative ion activation technique for phosphopeptide analysis, termed femtosecond laser-induced ionization/dissociation (fs-LID). In contrast to CID and ETD, fs-LID is shown to be particularly suited to the analysis of singly protonated phosphopeptide ions, yielding a wide range of product ions including a, b, c, x, y, and z sequence ions, as well as ions that are potentially diagnostic of the positions of phosphorylation (e.g., 'a(n)+1-98'). Importantly, the lack of phosphate moiety losses or phosphate group 'scrambling' provides unambiguous information for sequence identification and phosphorylation site characterization. Therefore, fs-LID-MS/MS can serve as a complementary technique to established methodologies for phosphoproteomic analysis.

Confident phosphorylation site localization using the Mascot Delta Score

Large scale phosphorylation analysis is more and more getting into focus of proteomic research. Although it is now possible to identify thousands of phosphorylated peptides in a biological system, confident site localization remains challenging. Here we validate the Mascot Delta Score (MD-score) as a simple method that achieves similar sensitivity and specificity for phosphosite localization as the published Ascore, which is mainly used in conjunction with Sequest. The MD-score was evaluated using liquid chromatography-tandem MS data of 180 individually synthesized phosphopeptides with precisely known phosphorylation sites. We tested the MD-score for a wide range of commonly available fragmentation methods and found it to be applicable throughout with high statistical significance. However, the different fragmentation techniques differ strongly in their ability to localize phosphorylation sites. At 1% false localization rate, the highest number of correctly assigned phosphopeptides was achieved by higher energy collision induced dissociation in combination with an Orbitrap mass analyzer followed very closely by low resolution ion trap spectra obtained after electron transfer dissociation. Both these methods are significantly better than low resolution spectra acquired after collision induced dissociation and multi stage activation. Score thresholds determined from simple calibration functions for each fragmentation method were stable over replicate analyses of the phosphopeptide set. The MD-score outperforms the Ascore for tyrosine phosphorylated peptides and we further show that the ability to call sites correctly increases with increasing distance of two candidate sites within a peptide sequence. The MD-score does not require complex computational steps which makes it attractive in terms of practical utility. We provide all mass spectra and the synthetic peptides to the community so that the development of present and future localization software can be benchmarked and any laboratory can determine MD-scores and localization probabilities for their individual analytical set up.