Mining a tandem mass spectrometry database to determine the trends and global factors influencing peptide fragmentation

Basophile: accurate fragment charge state prediction improves peptide identification rates

In shotgun proteomics, database search algorithms rely on fragmentation models to predict fragment ions that should be observed for a given peptide sequence. The most widely used strategy (Naive model) is oversimplified, cleaving all peptide bonds with equal probability to produce fragments of all charges below that of the precursor ion. More accurate models, based on fragmentation simulation, are too computationally intensive for on-the-fly use in database search algorithms. We have created an ordinal-regression-based model called Basophile that takes fragment size and basic residue distribution into account when determining the charge retention during CID/higher-energy collision induced dissociation (HCD) of charged peptides. This model improves the accuracy of predictions by reducing the number of unnecessary fragments that are routinely predicted for highly-charged precursors. Basophile increased the identification rates by 26% (on average) over the Naive model, when analyzing triply-charged precursors from ion trap data. Basophile achieves simplicity and speed by solving the prediction problem with an ordinal regression equation, which can be incorporated into any database search software for shotgun proteomic identification.

Influence of N-terminal residue composition on the structure of proline-containing b2+ ions

To probe the structural implications of the proline residue on its characteristic peptide fragmentation patterns, in particular its unusual cleavage at its C-terminus in formation of a b(2) ion in XxxProZzz sequences, the structures of a series of proline-containing b(2)(+) ions were studied by using action infrared multiphoton dissociation (IRMPD) spectroscopy and fragment ion hydrogen-deuterium exchange (HDX). Five different Xxx-Pro b(2)(+) ions were studied, with glycine, alanine, isoleucine, valine, or histidine in the N-terminal position. The residues selected feature different sizes, chain lengths, and gas phase basicities to explore whether the structure of the N-terminal residue influences the Xxx-Pro b(2)(+) ion structure. In proteins, the proline side chain-to-backbone attachment causes its peptide bonds to be in the cis conformation more than any other amino acid, although trans is still favored over cis. However, HP is the only b(2)(+) ion studied here that forms the diketopiperazine exclusively. The GP, AP, IP, and VP b(2)(+) ions formed from protonated tripeptide precursors predominantly featured oxazolone structures with small diketopiperazine contributions. In contrast to the b(2)(+) ions generated from tripeptides, synthetic cyclic dipeptides VP and HP were confirmed to have exclusive diketopiperazine structures.

Improving Peptide identification using empirical scoring systems

Peptides and proteins are routinely identified from peptide fragmentation spectra acquired in a mass spectrometer, analyzed by database search engines. The types of fragments that can be formed are known, and it is also well appreciated that certain fragment types are more common or more informative than others. However, most search engines do not use detailed knowledge of peptide fragmentation, but rather consider a limited range of fragments, giving each an equivalent weighting in their scoring system that decides which results are likely to be correct. This chapter discusses efforts to make use of information about the frequency of observation of different fragment ion types in order to produce more sophisticated and sensitive scoring systems and demonstrates how these new scoring systems are particularly powerful for analysis of electron capture or electron transfer dissociation data.

Interpretation of tandem mass spectra of posttranslationally modified peptides

Tandem mass spectrometry provides a sensitive means of analyzing the amino acid sequence of peptides and modified peptides by providing accurate mass measurements of precursor and fragment ions. Modern mass spectrometry instrumentation is capable of rapidly generating many thousands of tandem mass spectra and protein database search engines have been developed to match the experimental data to peptide candidates. In most studies there is a schism between discarding perfectly valid data and including nonsensical peptide identifications-this is currently a major bottleneck in data-analysis and it calls for an understanding of tandem mass spectrometry data. Manual evaluation of the data and perhaps experimental cross-checking of the MS data can save many months of experimental work trying to do biological follow-ups based on erroneous identifications. Especially for posttranslationally modified peptides there is a need for manual validation of the data because search algorithms seldom have been optimized for the identification of modified peptides and because there are many pitfalls for the unwary. This chapter describes some of the issues that should be considered when interpreting and validating tandem mass spectra and gives some useful tables to aid this process.

Algorithms for database-dependent search of MS/MS data

The frequent used bottom-up strategy for identification of proteins and their associated modifications generate nowadays typically thousands of MS/MS spectra that normally are matched automatically against a protein sequence database. Search engines that take as input MS/MS spectra and a protein sequence database are referred as database-dependent search engines. Many programs both commercial and freely available exist for database-dependent search of MS/MS spectra and most of the programs have excellent user documentation. The aim here is therefore to outline the algorithm strategy behind different search engines rather than providing software user manuals. The process of database-dependent search can be divided into search strategy, peptide scoring, protein scoring, and finally protein inference. Most efforts in the literature have been put in to comparing results from different software rather than discussing the underlining algorithms. Such practical comparisons can be cluttered by suboptimal implementation and the observed differences are frequently caused by software parameters settings which have not been set proper to allow even comparison. In other words an algorithmic idea can still be worth considering even if the software implementation has been demonstrated to be suboptimal. The aim in this chapter is therefore to split the algorithms for database-dependent searching of MS/MS data into the above steps so that the different algorithmic ideas become more transparent and comparable. Most search engines provide good implementations of the first three data analysis steps mentioned above, whereas the final step of protein inference are much less developed for most search engines and is in many cases performed by an external software. The final part of this chapter illustrates how protein inference is built into the VEMS search engine and discusses a stand-alone program SIR for protein inference that can import a Mascot search result.

Learning from decoys to improve the sensitivity and specificity of proteomics database search results

The statistical validation of database search results is a complex issue in bottom-up proteomics. The correct and incorrect peptide spectrum match (PSM) scores overlap significantly, making an accurate assessment of true peptide matches challenging. Since the complete separation between the true and false hits is practically never achieved, there is need for better methods and rescoring algorithms to improve upon the primary database search results. Here we describe the calibration and False Discovery Rate (FDR) estimation of database search scores through a dynamic FDR calculation method, FlexiFDR, which increases both the sensitivity and specificity of search results. Modelling a simple linear regression on the decoy hits for different charge states, the method maximized the number of true positives and reduced the number of false negatives in several standard datasets of varying complexity (18-mix, 49-mix, 200-mix) and few complex datasets (E. coli and Yeast) obtained from a wide variety of MS platforms. The net positive gain for correct spectral and peptide identifications was up to 14.81% and 6.2% respectively. The approach is applicable to different search methodologies- separate as well as concatenated database search, high mass accuracy, and semi-tryptic and modification searches. FlexiFDR was also applied to Mascot results and showed better performance than before. We have shown that appropriate threshold learnt from decoys, can be very effective in improving the database search results. FlexiFDR adapts itself to different instruments, data types and MS platforms. It learns from the decoy hits and sets a flexible threshold that automatically aligns itself to the underlying variables of data quality and size.

The use of selected reaction monitoring in quantitative proteomics

Selected reaction monitoring (SRM) has a long history of use in the area of quantitative MS. In recent years, the approach has seen increased application to quantitative proteomics, facilitating multiplexed relative and absolute quantification studies in a variety of organisms. This article discusses SRM, after introducing the context of quantitative proteomics (specifically primarily absolute quantification) where it finds most application, and considers topics such as the theory and advantages of SRM, the selection of peptide surrogates for protein quantification, the design of optimal SRM co-ordinates and the handling of SRM data. A number of published studies are also discussed to demonstrate the impact that SRM has had on the field of quantitative proteomics.

TUTORIAL ON DE NOVO PEPTIDE SEQUENCING USING MS/MS MASS SPECTROMETRY

This paper is a self-contained introductory tutorial on the problem in proteomics known as peptide sequencing using tandem mass spectrometry. This tutorial deals specifically with de novo sequencing methods (as opposed to database search methods). We first give an introduction to peptide sequencing, its importance and history and some background on proteins. Next we show the relationship between a peptide and the final spectrum produced from a tandem mass spectrometer, together with a description of the various sources of complications that arise during the process of generating the mass spectrum. From there we model the computational problem of de novo peptide sequencing, which is basically the reverse problem of identifying the peptide which produced the spectrum. We then present several major approaches to solve it (including reviewing some of the current algorithms in each approach), and also discuss related problems and post-processing approaches.

Investigation of scrambled ions in tandem mass spectra. Part 1. Statistical characterization

Scrambled ions have become the focus of recent investigations of peptide fragmentation. Here, an investigation of more than 390,000 high quality CID mass spectra is presented to explore the extent of scrambled ions in mass spectra and the possible fragmentation rules during scramble reactions. For the former, scrambled ions generally make up more than 10 % of mass spectra in number, although the abundances are less than 0.1 of the base peak. For the latter, relatively preferential re-opening sites were found for aliphatic residues Ala, Ile, Leu, and other residues such as Met, Gln, Ser, Phe, and Thr, whereas disfavored sites were found for basic residues Arg, Lys, and His, and Trp for both scrambled b and a ions. Similar preferential order in re-opening reaction was found in the reaction of losing internal residues when cleavage occurs at C-terminal side of 20 residues. However, when cleavage occurs at N-terminal side, Glu, Phe, and Trp become the most preferential sites. These results provide a deep insight into cleavage rules during scramble reactions for prediction of peptide mass spectra. Also, an additional investigation of whether scrambled ions could help discriminate false identifications from correct identifications was performed. Probing the number fraction of scrambled ions in falsely and correctly interpreted spectra and analyzing the correlation between scrambled ions and SEQUEST scores XCorr and Sp showed scrambled ions could at some extent help improve the discrimination in singly charged identifications, whereas no improvement was found for multiply charged results.

Formation of y + 10 and y + 11 ions in the collision-induced dissociation of peptide ions

Tandem mass spectra of peptide ions, acquired in shotgun proteomic studies of selected proteins, tissues, and organisms, commonly include prominent peaks that cannot be assigned to the known fragmentation product ions (y, b, a, neutral losses). In many cases these persist even when creating consensus spectra for inclusion in spectral libraries, where it is important to determine whether these peaks represent new fragmentation paths or arise from impurities. Using spectra from libraries and synthesized peptides, we investigate a class of fragment ions corresponding to y(n-1) + 10 and y(n-1) + 11, where n is the number of amino acid residues in the peptide. These 10 and 11 Da differences in mass of the y ion were ascribed before to the masses of [+ CO - H(2)O] and [+ CO - NH(3)], respectively. The mechanism is suggested to involve dissociation of the N-terminal residue at the CH-CO bond following loss of H(2)O or NH(3). MS(3) spectra of these ions show that the location of the additional 10 or 11 Da is at the N-terminal residue. The y(n-1) + 10 ion is most often found in peptides with N-terminal proline, asparagine, and histidine, and also with serine and threonine in the adjacent position. The y(n-1) + 11 ion is observed predominantly with histidine and asparagine at the N-terminus, but also occurs with asparagine in positions two through four. The intensities of the y(n-1) + 10 ions decrease with increasing peptide length. These data for y(n-1) + 10 and y(n-1) + 11 ion formation may be used to improve peptide identification from tandem mass spectra.

Occurrence of C-terminal residue exclusion in peptide fragmentation by ESI and MALDI tandem mass spectrometry

By screening a data set of 392 synthetic peptides MS/MS spectra, we found that a known C-terminal rearrangement was unexpectedly frequently occurring from monoprotonated molecular ions in both ESI and MALDI tandem mass spectrometry upon low and high energy collision activated dissociations with QqTOF and TOF/TOF mass analyzer configuration, respectively. Any residue localized at the C-terminal carboxylic acid end, even a basic one, was lost, provided that a basic amino acid such arginine and to a lesser extent histidine and lysine was present in the sequence leading to a fragment ion, usually depicted as (b(n-1) + H(2)O) ion, corresponding to a shortened non-scrambled peptide chain. Far from being an epiphenomenon, such a residue exclusion from the peptide chain C-terminal extremity gave a fragment ion that was the base peak of the MS/MS spectrum in certain cases. Within the frame of the mobile proton model, the ionizing proton being sequestered onto the basic amino acid side chain, it is known that the charge directed fragmentation mechanism involved the C-terminal carboxylic acid function forming an anhydride intermediate structure. The same mechanism was also demonstrated from cationized peptides. To confirm such assessment, we have prepared some of the peptides that displayed such C-terminal residue exclusion as a C-terminal backbone amide. As expected in this peptide amide series, the production of truncated chains was completely suppressed. Besides, multiply charged molecular ions of all peptides recorded in ESI mass spectrometry did not undergo such fragmentation validating that any mobile ionizing proton will prevent such a competitive C-terminal backbone rearrangement. Among all well-known nondirect sequence fragment ions issued from non specific loss of neutral molecules (mainly H(2)O and NH(3)) and multiple backbone amide ruptures (b-type internal ions), the described C-terminal residue exclusion is highly identifiable giving raise to a single fragment ion in the high mass range of the MS/MS spectra. The mass difference between this signal and the protonated molecular ion corresponds to the mass of the C-terminal residue. It allowed a straightforward identification of the amino acid positioned at this extremity. It must be emphasized that a neutral residue loss can be misattributed to the formation of a y(m-1) ion, i.e., to the loss of the N-terminal residue following the a(1)-y(m-1) fragmentation channel. Extreme caution must be adopted when reading the direct sequence ion on the positive ion MS/MS spectra of singly charged peptides not to mix up the attribution of the N- and C-terminal amino acids. Although such peculiar fragmentation behavior is of obvious interest for de novo peptide sequencing, it can also be exploited in proteomics, especially for studies involving digestion protocols carried out with proteolytic enzymes other than trypsin (Lys-N, Glu-C, and Asp-N) that produce arginine-containing peptides.

Plasma proteins present in osteoarthritic synovial fluid can stimulate cytokine production via Toll-like receptor 4

Introduction

Osteoarthritis (OA) is a degenerative disease characterized by cartilage breakdown in the synovial joints. The presence of low-grade inflammation in OA joints is receiving increasing attention, with synovitis shown to be present even in the early stages of the disease. How the synovial inflammation arises is unclear, but proteins in the synovial fluid of affected joints could conceivably contribute. We therefore surveyed the proteins present in OA synovial fluid and assessed their immunostimulatory properties.

Methods

We used mass spectrometry to survey the proteins present in the synovial fluid of patients with knee OA. We used a multiplex bead-based immunoassay to measure levels of inflammatory cytokines in serum and synovial fluid from patients with knee OA and from patients with rheumatoid arthritis (RA), as well as in sera from healthy individuals. Significant differences in cytokine levels between groups were determined by significance analysis of microarrays, and relations were determined by unsupervised hierarchic clustering. To assess the immunostimulatory properties of a subset of the identified proteins, we tested the proteins' ability to induce the production of inflammatory cytokines by macrophages. For proteins found to be stimulatory, the macrophage stimulation assays were repeated by using Toll-like receptor 4 (TLR4)-deficient macrophages.

Results

We identified 108 proteins in OA synovial fluid, including plasma proteins, serine protease inhibitors, proteins indicative of cartilage turnover, and proteins involved in inflammation and immunity. Multiplex cytokine analysis revealed that levels of several inflammatory cytokines were significantly higher in OA sera than in normal sera, and levels of inflammatory cytokines in synovial fluid and serum were, as expected, higher in RA samples than in OA samples. As much as 36% of the proteins identified in OA synovial fluid were plasma proteins. Testing a subset of these plasma proteins in macrophage stimulation assays, we found that Gc-globulin, α₁-microglobulin, and α₂-macroglobulin can signal via TLR4 to induce macrophage production of inflammatory cytokines implicated in OA.

Conclusions

Our findings suggest that plasma proteins present in OA synovial fluid, whether through exudation from plasma or production by synovial tissues, could contribute to low-grade inflammation in OA by functioning as so-called damage-associated molecular patterns in the synovial joint.

Mechanism elucidation for nonstochastic femtosecond laser-induced ionization/dissociation: from amino acids to peptides

Femtosecond laser-induced ionization/dissociation (fs-LID) has been demonstrated as a novel ion activation method for use in tandem mass spectrometry. The technique opens the door to unique structural information about biomolecular samples that is not easily accessed by traditional means. fs-LID is able to cleave strong bonds while keeping weaker bonds intact. This feature has been found to be particularly useful for the mapping of post-translational modifications such as phosphorylation, which is difficult to achieve by conventional proteomic studies. Here we investigate the laser-ion interaction on a fundamental level through the characterization of fs-LID spectra for the protonated amino acids and two series of derivatized samples. The findings are used to better understand the fs-LID spectra of synthetic peptides. This is accomplished by exploring the effects of several single-residue substitutions.

A high-confidence human plasma proteome reference set with estimated concentrations in PeptideAtlas

Human blood plasma can be obtained relatively noninvasively and contains proteins from most, if not all, tissues of the body. Therefore, an extensive, quantitative catalog of plasma proteins is an important starting point for the discovery of disease biomarkers. In 2005, we showed that different proteomics measurements using different sample preparation and analysis techniques identify significantly different sets of proteins, and that a comprehensive plasma proteome can be compiled only by combining data from many different experiments. Applying advanced computational methods developed for the analysis and integration of very large and diverse data sets generated by tandem MS measurements of tryptic peptides, we have now compiled a high-confidence human plasma proteome reference set with well over twice the identified proteins of previous high-confidence sets. It includes a hierarchy of protein identifications at different levels of redundancy following a clearly defined scheme, which we propose as a standard that can be applied to any proteomics data set to facilitate cross-proteome analyses. Further, to aid in development of blood-based diagnostics using techniques such as selected reaction monitoring, we provide a rough estimate of protein concentrations using spectral counting. We identified 20,433 distinct peptides, from which we inferred a highly nonredundant set of 1929 protein sequences at a false discovery rate of 1%. We have made this resource available via PeptideAtlas, a large, multiorganism, publicly accessible compendium of peptides identified in tandem MS experiments conducted by laboratories around the world.

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.

Tandem mass spectrometry strategies for phosphoproteome analysis

Protein phosphorylation is involved in nearly all essential biochemical pathways and the deregulation of phosphorylation events has been associated with the onset of numerous diseases. A multitude of tandem mass spectrometry (MS/MS) and multistage MS/MS (i.e., MS(n) ) strategies have been developed in recent years and have been applied toward comprehensive phosphoproteomic analysis, based on the interrogation of proteolytically derived phosphopeptides. However, the utility of each of these MS/MS and MS(n) approaches for phosphopeptide identification and characterization, including phosphorylation site localization, is critically dependant on the properties of the precursor ion (e.g., polarity and charge state), the specific ion activation method that is employed, and the underlying gas-phase ion chemistries, mechanisms and other factors that influence the gas-phase fragmentation behavior of phosphopeptide ions. This review therefore provides an overview of recent studies aimed at developing an improved understanding of these issues, and highlights the advantages and limitations of both established (e.g., CID) and newly maturing (e.g., ECD, ETD, photodissociation, etc.) yet complementary, ion activation techniques. This understanding is expected to facilitate the continued refinement of existing MS/MS strategies, and the development of novel MS/MS techniques for phosphopeptide analysis, with great promise in providing new insights into the role of protein phosphorylation on normal biological function, and in the onset and progression of disease. © 2011 Wiley Periodicals, Inc., Mass Spec Rev 30:600-625, 2011.

CrossWork: software-assisted identification of cross-linked peptides

The increased interest in chemical cross-linking for probing protein structure and interaction has led to a large increase in literature describing new cross-linkers and search programs. However, this has not led to a corresponding increase in the analysis of large and complex proteins. A major obstacle is that the new cross-linkers are either not readily available and/or have a low reactivity. In combination with aging search programs that are slow and have low sensitivity, or new search programs that are described but not released, these efforts do little to advance the field of cross-linking. Here we present a method pipeline for chemical cross-linking, using two standard cross-linkers, BS3 and BS2G, combined with our freely available CrossWork search program. By this approach we generate cross-link data sufficient to derive structural information for large and complex proteins. CrossWork searches batches of tandem mass-spectrometric data, and identifies cross-linked and non-cross-linked peptides using a standard PC. We tested CrossWork by searching mass-spectrometric datasets of cross-linked complement factor C3 against small (1 protein) and large (1000 proteins) search spaces, and show that the resulting distance constraints agree with the established structures. We further investigated the structure of the multi-domain ERp72, and combined the individual domains of ERp72 into a single structure.

Charge states of y ions in the collision-induced dissociation of doubly charged tryptic peptide ions

Bonds that break in collision-induced dissociation (CID) are often weakened by a nearby proton, which can, in principle, be carried away by either of the product fragments. Since peptide backbone dissociation is commonly charge-directed, relative intensities of charge states of product y- and b-ions depend on the final location of that proton. This study examines y-ion charge distributions for dissociation of doubly charged peptide ions, using a large reference library of peptide ion fragmentation generated from ion-trap CID of peptide ions from tryptic digests. Trends in relative intensities of y(2+) and y(1+) ions are examined as a function of bond cleavage position, peptide length (n), residues on either side of the bond and effects of residues remote from the bond. It is found that y(n-2)/b(2) dissociation is the most sensitive to adjacent amino acids, that y(2+)/y(1+) steadily increase with increasing peptide length, that the N-terminal amino acid can have a major influence in all dissociations, and in some cases other residues remote from the bond cleavage exert significant effects. Good correlation is found between the values of y(2+)/y(1+) for the peptide and the proton affinities of the amino acids present at the dissociating peptide bond. A few deviations from this correlation are rationalized by specific effects of the amino acid residues. These correlations can be used to estimate trends in y(2+)/y(1+) ratios for peptide ions from amino acid proton affinities.

Utility of characteristic QTOF MS/MS fragmentation for MHC class I peptides

Systematic investigation of cellular process by mass spectrometric detection of peptides obtained from proteins digestion or directly from immuno-purification can be a powerful tool when used appropriately. The true sequence of these peptides is defined by the interpretation of spectral data using a variety of available algorithms. However peptide match algorithm scoring is typically based on some, but not all, of the mechanisms of peptide fragmentation. Although algorithm rules for soft ionization techniques generally fit very well to tryptic peptides, manual validation of spectra is often required for endogenous peptides such as MHC class I molecules where traditional trypsin digest techniques are not used. This study summarizes data mining and manual validation of hundreds of peptide sequences from MHC class I molecules in publically available data files. We herein describe several important features to improve and quantify manual validation for these endogenous peptides--post automated algorithm searching. Important fragmentation patterns are discussed for the studied MHC Class I peptides. These findings lead to practical rules that are helpful when performing manual validation. Furthermore, these observations may be useful to improve current peptide search algorithms or development of novel software tools.

A global analysis of peptide fragmentation variability

Understanding the fragmentation process in MS/MS experiments is vital when trying to validate the results of such experiments, and one way of improving our understanding is to analyze existing data. We here present our findings from an analysis of a large and diverse data set of MS/MS-based peptide identifications, in which each peptide has been identified from multiple spectra, recorded on two commonly used types of electrospray instruments. By analyzing these data we were able to study fragmentation variability on three levels: (i) variation in detection rates and intensities for fragment ions from the same peptide sequence measured multiple times on a single instrument; (ii) consistency of rank-based fragmentation patterns; and (iii) a set of general observations on fragment ion occurrence in MS/MS experiments, regardless of sequence. Our results confirm that substantial variation can be found at all levels, even when high-quality identifications are used and the experimental conditions as well as the peptide sequences are kept constant. Finally, we discuss the observed variability in light of ongoing efforts to create spectral libraries and predictive software for target selection in targeted proteomics.

Sequencing Lys-N proteolytic peptides by ESI and MALDI tandem mass spectrometry

In this study, we explored the MS/MS behavior of various synthetic peptides that possess a lysine residue at the N-terminal position. These peptides were designed to mimic peptides produced upon proteolysis by the Lys-N enzyme, a metalloendopeptidase issued from a Japanese fungus Grifola frondosa that was recently investigated in proteomic studies as an alternative to trypsin digestion, as a specific cleavage at the amide X-Lys chain is obtained that provides N-terminal lysine peptide fragments. In contrast to tryptic peptides exhibiting a lysine or arginine residue solely at the C-terminal position, and are thus devoid of such basic amino acids within the sequence, these Lys-N proteolytic peptides can contain the highly basic arginine residue anywhere within the peptide chain. The fragmentation patterns of such sequences with the ESI-QqTOF and MALDI-TOF/TOF mass spectrometers commonly used in proteomic bottom-up experiments were investigated.

Bioinformatics tools for the structural elucidation of multi-subunit protein complexes by mass spectrometric analysis of protein-protein cross-links

Multi-subunit protein complexes are involved in many essential biochemical processes including signal transduction, protein synthesis, RNA synthesis, DNA replication and protein degradation. An accurate description of the relative structural arrangement of the constituent subunits in such complexes is crucial for an understanding of the molecular mechanism of the complex as a whole. Many complexes, however, lie in the mega-Dalton range, and are not amenable to X-ray crystallographic or nuclear magnetic resonance analysis. Techniques that are suited to structural studies of such large complexes, such as cryo-electron microscopy, do not provide the resolution required for a mechanistic insight. Mass spectrometry (MS) has increasingly been applied to identify the residues that are involved in chemical cross-links in compound protein assemblies, and have provided valuable insight into the molecular arrangement, orientation and contact surfaces of subunits within such large complexes. This approach is known as MS3D, and involves the MS analysis of cross-linked di-peptides following the enzymatic cleavage of a chemically cross-linked complex. A major challenge of this approach is the identification of the cross-linked di-peptides in a composite mixture of peptides, as well as the identification of the residues involved in the cross-link. These analyses require bioinformatics tools with capabilities beyond that of general, MS-based proteomic analysis software. Many MS3D software tools have appeared, often designed for very specific experimental methods. Here, we provide a review of all major MS3D bioinformatics programmes, reviewing their applicability to different workflows, specific experimental requirements and the computational approach taken by each.

Effects of estrogen on beta-amyloid-induced cholinergic cell death in the nucleus basalis magnocellularis

Alzheimer disease is characterized by accumulation of β-amyloid (Aβ) and cognitive dysfunctions linked to early loss of cholinergic neurons. As estrogen-based hormone replacement therapy has beneficial effects on cognition of demented patients, and it may prevent memory impairments, we investigated the effect of estrogen-pretreatment on Aβ-induced cholinergic neurodegeneration in the nucleus basalis magnocellularis (NBM). We tested which Aβ species induces the more pronounced cholinotoxic effect in vivo. We injected different Aβ assemblies in the NBM of mice, and measured cholinergic cell and cortical fiber loss. Spherical Aβ oligomers had the most toxic effect. Pretreatment of ovariectomized mice with estrogen before Aβ injection decreased cholinergic neuron loss and partly prevented fiber degeneration. By using proteomics, we searched for proteins involved in estrogen-mediated protection and in Aβ toxicity 24 h following injection. The change in expression of, e.g., DJ-1, NADH ubiquinone oxidoreductase, ATP synthase, phosphatidylethanolamine-binding protein 1, protein phosphatase 2A and dimethylarginine dimethylaminohydrolase 1 support our hypothesis that Aβ induces mitochondrial dysfunction, decreases MAPK signaling, and increases NOS activation in NBM. On the other hand, altered expression of, e.g., MAP kinase kinase 1 and 2, protein phosphatase 1 and 2A by Aβ might increase MAPK suppression and NOS signaling in the cortical target area. Estrogen pretreatment reversed most of the changes in the proteome in both areas. Our experiments suggest that regulation of the MAPK pathway, mitochondrial pH and NO production may all contribute to Aβ toxicity, and their regulation can be prevented partly by estrogen pretreatment.

Low-molecular weight plasma proteome analysis using centrifugal ultrafiltration

The low-molecular weight fraction (LMF) of the human plasma proteome is an invaluable source of biological information, especially in the context of identifying plasma-based biomarkers of disease. This protocol outlines a standardized procedure for the rapid/reproducible LMF profiling of human plasma samples using centrifugal ultrafiltration fractionation, followed by 1D-SDS-PAGE separation and nano-LC-MS/MS. Ultrafiltration is a convective process that uses anisotropic semipermeable membranes to separate macromolecular species on the basis of size. We have optimized centrifugal ultrafiltration for plasma fractionation with respect to buffer and solvent composition, centrifugal force, duration and temperature to facilitate >95% recovery, and enrichment of low-M (r) components from human plasma. Using this protocol, >260 unique peptides can be identified from a single plasma profiling experiment using 100 μL of plasma (Greening and Simpson, J Proteomics 73:637-648, 2010). The efficacy of this method is demonstrated by the identification, for the first time, of several plasma proteins (e.g., protein KIAA0649 (Q9Y4D3), rheumatoid factor D5, serine protease inhibitor A3, and transmembrane adapter protein PAG) previously not reported in extant high-confidence Human Proteome Organization Plasma Proteome Project datasets.

Computational analysis of unassigned high-quality MS/MS spectra in proteomic data sets

In a typical shotgun proteomics experiment, a significant number of high-quality MS/MS spectra remain "unassigned." The main focus of this work is to improve our understanding of various sources of unassigned high-quality spectra. To achieve this, we designed an iterative computational approach for more efficient interrogation of MS/MS data. The method involves multiple stages of database searching with different search parameters, spectral library searching, blind searching for modified peptides, and genomic database searching. The method is applied to a large publicly available shotgun proteomic data set.

The relative charge ratio between C and N atoms in amide bond acts as a key factor to determine peptide fragment efficiency in different charge states

The influence of charge state on the peptide dissociation behavior in tandem mass spectrometry (MS/MS) is worthy of discussion. Comparative studies of singly- and doubly-protonated peptide molecules are performed to explore the effect and mechanism of charge state on peptide fragmentation. In view of the charge-directed cleavage of protonated peptides described in the mobile proton model, radiolytic oxidation was applied to change the charge distribution of peptides but retain the sequence. Experimental studies of collision energy-dependent fragmentation efficiencies coupled with quantum chemical calculations indicated that the cleavage of ARRA and its side-chain oxidation products with oxygen atoms added followed a trend that doubly-protonated peptides fragment more easily than singly-protonated forms, while the oxidation product with the guanidine group deleted showed the opposite trend. By analyzing the charge distribution around the amide bonds, we found that the relative charge ratios between C and N atoms (Q(C)/Q(N)) in the amide bonds provided a reasonable explanation for peptide fragmentation efficiencies. An increase of the Q(C)/Q(N) value of the amide bond means that a peptide fragments more easily, and vice versa. The results described in this paper provide an experimental and calculation strategy for predicting peptide fragmentation efficiency.

Harvest: an open-source tool for the validation and improvement of peptide identification metrics and fragmentation exploration

Background

Protein identification using mass spectrometry is an important tool in many areas of the life sciences, and in proteomics research in particular. Increasing the number of proteins correctly identified is dependent on the ability to include new knowledge about the mass spectrometry fragmentation process, into computational algorithms designed to separate true matches of peptides to unidentified mass spectra from spurious matches. This discrimination is achieved by computing a function of the various features of the potential match between the observed and theoretical spectra to give a numerical approximation of their similarity. It is these underlying "metrics" that determine the ability of a protein identification package to maximise correct identifications while limiting false discovery rates. There is currently no software available specifically for the simple implementation and analysis of arbitrary novel metrics for peptide matching and for the exploration of fragmentation patterns for a given dataset.

Results

We present Harvest: an open source software tool for analysing fragmentation patterns and assessing the power of a new piece of information about the MS/MS fragmentation process to more clearly differentiate between correct and random peptide assignments. We demonstrate this functionality using data metrics derived from the properties of individual datasets in a peptide identification context. Using Harvest, we demonstrate how the development of such metrics may improve correct peptide assignment confidence in the context of a high-throughput proteomics experiment and characterise properties of peptide fragmentation.

Conclusions

Harvest provides a simple framework in C++ for analysing and prototyping metrics for peptide matching, the core of the protein identification problem. It is not a protein identification package and answers a different research question to packages such as Sequest, Mascot, X!Tandem, and other protein identification packages. It does not aim to maximise the number of assigned peptides from a set of unknown spectra, but instead provides a method by which researchers can explore fragmentation properties and assess the power of novel metrics for peptide matching in the context of a given experiment. Metrics developed using Harvest may then become candidates for later integration into protein identification packages.

Mobile protons versus mobile radicals: gas-phase unimolecular chemistry of radical cations of cysteine-containing peptides

A combination of electrospray ionization (ESI), multistage, and high-resolution mass spectrometry experiments are used to examine the gas-phase fragmentation reactions of radical cations of cysteine containing di- and tripeptides. Two different chemical methods were used to form initial populations of radical cations in which the radical sites were located at different positions: (1) sulfur-centered cysteinyl radicals via bond homolysis of protonated S-nitrosocysteine containing peptides; and (2) alpha-carbon backbone-centered radicals via Siu's sequence of reactions (J. Am. Chem. Soc.2008, 130, 7862). Comparison of the fragmentation reactions of these regiospecifically generated radicals suggests that hydrogen atom transfer (HAT) between the alpha C-H of adjacent residues and the cysteinyl radical can occur. In addition, using accurate mass measurements, deuterium labeling, and comparison with an authentic sample, a novel loss of part of the N-terminal cysteine residue was shown to give rise to the protonated, truncated N-formyl peptide (an even-electron x(n) ion). DFT calculations were performed on the radical cation [GCG]*(+) to examine: the relative stabilities of isomers with different radical and protonation sites; the barriers associated with radical migration between four possible radical sites, [G*CG](+), [GC*G](+), [GCG*](+), and [GC(S*)G](+); and for dissociation from these sites to yield b(2)-type ions.

'Fixed charge' chemical derivatization and data dependant multistage tandem mass spectrometry for mapping protein surface residue accessibility

Protein surface accessible residues play an important role in protein folding, protein-protein interactions and protein-ligand binding. However, a common problem associated with the use of selective chemical labeling methods for mapping protein solvent accessible residues is that when a complicated peptide mixture resulting from a large protein or protein complex is analyzed, the modified peptides may be difficult to identify and characterize amongst the largely unmodified peptide population (i.e., the 'needle in a haystack' problem). To address this challenge, we describe here the development of a strategy involving the synthesis and application of a novel 'fixed charge' sulfonium ion containing lysine-specific protein modification reagent, S,S'-dimethylthiobutanoylhydroxysuccinimide ester (DMBNHS), coupled with capillary HPLC-ESI-MS, automated CID-MS/MS, and data-dependant neutral loss mode MS(3) in an ion trap mass spectrometer, to map the surface accessible lysine residues in a small model protein, cellular retinoic acid binding protein II (CRABP II). After reaction with different reagent:protein ratios and digestion with Glu-C, modified peptides are selectively identified and the number of modifications within each peptide are determined by CID-MS/MS, via the exclusive neutral loss(es) of dimethylsulfide, independently of the amino acid composition and precursor ion charge state (i.e., proton mobility) of the peptide. The observation of these characteristic neutral losses are then used to automatically 'trigger' the acquisition of an MS(3) spectrum to allow the peptide sequence and the site(s) of modification to be characterized. Using this approach, the experimentally determined relative solvent accessibilities of the lysine residues were found to show good agreement with the known solution structure of CRABP II.

Integrating Lys-N proteolysis and N-terminal guanidination for improved fragmentation and relative quantification of singly-charged ions

The study of isolated protein complexes has greatly benefited from recent advances in mass spectrometry instrumentation and quantitative, isotope labeling techniques. The comprehensive characterization of protein complex components and quantification of their relative abundance relies heavily upon maximizing protein and peptide sequence information obtained from MS and tandem MS studies. Recent work has shown that using a metalloendopeptidase, Lys-N, for proteomic analysis of biological protein mixtures produces complementary protein sequence information compared with trypsin digestion alone. Here, we have investigated the suitability of Lys-N proteolysis for use with MALDI mass spectrometry to characterize the yeast Arp2 complex and E. coli PAP I protein interactions. Although Lys-N digestion resulted in an average decrease in protein sequence coverage of approximately 30% compared with trypsin digestion, CID analysis of singly-charged Lys-N peptides yielded a more extensive b-ions series compared with complementary tryptic peptides. Taking advantage of this improved fragmentation pattern, we utilized differential (15)N/(14)N guanidination of Lys-N peptides and MALDI-MS/MS analysis to relatively quantify the changes in PAP I associations due to deletion of sprE, previously shown to regulate PAP I-dependent polyadenylation. Overall, this Lys-N/guanidination integrative approach is applicable for functional proteomic studies utilizing MALDI mass spectrometry analysis, as it provides an effective and economical mean for relative quantification of proteins in conjunction with increased sensitivity of detection and fragmentation efficiency.

Value of using multiple proteases for large-scale mass spectrometry-based proteomics

Large-scale protein sequencing methods rely on enzymatic digestion of complex protein mixtures to generate a collection of peptides for mass spectrometric analysis. Here we examine the use of multiple proteases (trypsin, LysC, ArgC, AspN, and GluC) to improve both protein identification and characterization in the model organism Saccharomyces cerevisiae. Using a data-dependent, decision tree-based algorithm to tailor MS(2) fragmentation method to peptide precursor, we identified 92 095 unique peptides (609 665 total) mapping to 3908 proteins at a 1% false discovery rate (FDR). These results were a significant improvement upon data from a single protease digest (trypsin) - 27 822 unique peptides corresponding to 3313 proteins. The additional 595 protein identifications were mainly from those at low abundances (i.e., < 1000 copies/cell); sequence coverage for these proteins was likewise improved nearly 3-fold. We demonstrate that large portions of the proteome are simply inaccessible following digestion with a single protease and that multiple proteases, rather than technical replicates, provide a direct route to increase both protein identifications and proteome sequence coverage.

Structural characterization of intact proteins is enhanced by prevalent fragmentation pathways rarely observed for peptides

While collisionally activated dissociation (CAD) pathways for peptides are well characterized, those of intact proteins are not. We systematically assigned CAD product ions of ubiquitin, myoglobin, and bovine serum albumin generated using high-yield, in-source fragmentation. Assignment of >98% of hundreds of product ions implies that the fragmentation pathways described are representative of the major pathways. Protein dissociation mechanisms were found to be modulated by both source declustering potential and precursor ion charge state. Like peptides, higher charge states of proteins fragmented at lower energies next to Pro, via mobile protons, while lower charge states fragmented at higher energies after Asp and Glu, via localized protons. Unlike peptides, however, predominant fragmentation channels of proteins occurred at intermediate charge states via non-canonical mechanisms and produced extensive internal fragmentation. The non-canonical mechanisms include prominent cleavages C-terminal to Pro and Asn, and N-terminal to Ile, Leu, and Ser; these cleavages, along with internal fragments, led to a 45% increase in sequence coverage, improving the specificity of top-down protein identification. Three applications take advantage of the different mechanisms of protein fragmentation. First, modulation of declustering potential selectively fragments different charge states, allowing the source region to be used as the first stage of a low-resolution tandem mass spectrometer, facilitating pseudo-MS(3) of product ions with known parent charge states. Second, development and integration of automated modulation of ion funnel declustering potential allows users access to a particular fragmentation mechanism, yielding facile cleavage on a liquid chromatography timescale. Third, augmentation of a top-down search engine improved protein characterization.

Statistical analysis of Peptide electron transfer dissociation fragmentation mass spectrometry

It is well established that protein sequence determination may be achieved by mass spectrometric analysis of protonated tryptic peptides subjected to collisional activation. When separated by nanoflow HPLC, a high percentage of peptides from complex mixtures of proteins can usually be identified. Recently, alternative, radical-driven fragmentation approaches of electron capture dissociation and the more common electron transfer dissociation (ETD) have been introduced and made widely available. In order to utilize these techniques in large scale proteomics studies, it is important to characterize the performance of these fragmentation processes on peptides formed by a range of enzymatic cleavages. In this study, we present a statistical analysis of the ion types that are observed from peptides produced by different enzymes and highlight the different characteristics of ETD spectra of doubly charged precursors in comparison to precursors of higher charge states.

Proteomics data repositories

The field of proteomics, particularly the application of MS analysis to protein samples, is well established and growing rapidly. Proteomic studies generate large volumes of raw experimental data and inferred biological results. To facilitate the dissemination of these data, centralized data repositories have been developed that make the data and results accessible to proteomic researchers and biologists alike. This review of proteomics data repositories focuses exclusively on freely available, centralized data resources that disseminate or store experimental MS data and results. The resources chosen reflect a current "snapshot" of the state of resources available with an emphasis placed on resources that may be of particular interest to yeast researchers. Resources are described in terms of their intended purpose and the features and functionality provided to users.

Understanding and exploiting Peptide fragment ion intensities using experimental and informatic approaches

Tandem mass spectrometry is a widely used tool in proteomics. This section will address the properties that describe how protonated peptides fragment when activated by collisions in a mass spectrometer and how that information can be used to identify proteins. A review of the mobile proton model is presented, along with a summary of commonly observed peptide cleavage enhancements, including the proline effect. The methods used to elucidate peptide dissociation chemistry by using both small groups of model peptides and large datasets are also discussed. Finally, the role of peak intensity in commercially available and developmental peptide identification algorithms is examined.

Correlation between y-type ions observed in ion trap and triple quadrupole mass spectrometers

Multiple reaction monitoring mass spectrometry (MRM-MS) is a technique for high-sensitivity targeted analysis. In proteomics, MRM-MS can be used to monitor and quantify a peptide based on the production of expected fragment peaks from the selected peptide precursor ion. The choice of which fragment ions to monitor in order to achieve maximum sensitivity in MRM-MS can potentially be guided by existing MS/MS spectra. However, because the majority of discovery experiments are performed on ion trap platforms, there is concern in the field regarding the generalizability of these spectra to MRM-MS on a triple quadrupole instrument. In light of this concern, many operators perform an optimization step to determine the most intense fragments for a target peptide on a triple quadrupole mass spectrometer. We have addressed this issue by targeting, on a triple quadrupole, the top six y-ion peaks from ion trap-derived consensus library spectra for 258 doubly charged peptides from three different sample sets and quantifying the observed elution curves. This analysis revealed a strong correlation between the y-ion peak rank order and relative intensity across platforms. This suggests that y-type ions obtained from ion trap-based library spectra are well-suited for generating MRM-MS assays for triple quadrupoles and that optimization is not required for each target peptide.

Peptide photodissociation with 157 nm light in a commercial tandem time-of-flight mass spectrometer

Photodissociation with 157 nm light was implemented in an ABI model 4700 matrix-assisted laser desorption ionization (MALDI) tandem time-of-flight (TOF) mass spectrometer for peptide analysis. With a homemade computer program to control the light timing based on the m/z of each precursor ion, the photodissociation setup was seamlessly automated with the mass spectrometer. Peptide photodissociation in this apparatus yielded fragments similar to those observed in previous experiments with a home-built tandem-TOF mass spectrometer. Peptides having arginine at their C-termini yielded high-energy x-, v-, and w- type fragments, while peptides with N-terminal arginine produced many a- and d- type ions. Abundant immonium ions were also generated. High-quality photodissociation spectra were obtained with as little as 5 fmol of peptides. In the analysis of various tryptic peptides, photodissociation provided much more sequence information than the conventional TOF-TOF collision induced dissociation (CID). Because of the high fragmentation efficiency, sensitivity was not sacrificed to achieve this.

Analysis of human C1q by combined bottom-up and top-down mass spectrometry: detailed mapping of post-translational modifications and insights into the C1r/C1s binding sites

C1q is a subunit of the C1 complex, a key player in innate immunity that triggers activation of the classical complement pathway. Featuring a unique structural organization and comprising a collagen-like domain with a high level of post-translational modifications, C1q represents a challenging protein assembly for structural biology. We report for the first time a comprehensive proteomics study of C1q combining bottom-up and top-down analyses. C1q was submitted to proteolytic digestion by a combination of collagenase and trypsin for bottom-up analyses. In addition to classical LC-MS/MS analyses, which provided reliable identification of hydroxylated proline and lysine residues, sugar loss-triggered MS(3) scans were acquired on an LTQ-Orbitrap (Linear Quadrupole Ion Trap-Orbitrap) instrument to strengthen the localization of glucosyl-galactosyl disaccharide moieties on hydroxylysine residues. Top-down analyses performed on the same instrument allowed high accuracy and high resolution mass measurements of the intact full-length C1q polypeptide chains and the iterative fragmentation of the proteins in the MS(n) mode. This study illustrates the usefulness of combining the two complementary analytical approaches to obtain a detailed characterization of the post-translational modification pattern of the collagen-like domain of C1q and highlights the structural heterogeneity of individual molecules. Most importantly, three lysine residues of the collagen-like domain, namely Lys(59) (A chain), Lys(61) (B chain), and Lys(58) (C chain), were unambiguously shown to be completely unmodified. These lysine residues are located about halfway along the collagen-like fibers. They are thus fully available and in an appropriate position to interact with the C1r and C1s protease partners of C1q and are therefore likely to play an essential role in C1 assembly.

Fragmentation of doubly-protonated Pro-His-Xaa tripeptides: formation of b(2)(2+) ions

When ionized by electrospray from acidic solutions, the tripeptides Pro-His-Xaa (Xaa = Gly, Ala, Leu) form abundant doubly-protonated ions, [M + 2H]2+. Collision-induced dissociation (CID) of these doubly-protonated species results, in part, in formation of b(2)(2+) ions, which fragment further by loss of CO to form a(2)(2+) ions; the latter fragment by loss of CO to form the Pro and His iminium [immonium is commonly used in peptide MS work] ions. Although larger doubly-charged b ions are known, this represents the first detailed study of b(2)(2+) ions in CID of small doubly protonated peptides. The most abundant CID products of the studied doubly-protonated peptides arise mainly in charge separation involving two primary fragmentation channels, formation of the b2/y1 pair and formation of the a1/y2 pair. Combined molecular dynamics and density functional theory calculations are used to gain insight into the structures and fragmentation pathways of doubly-protonated Pro-His-Gly including the energetics of potential protonation sites, backbone cleavages, post-cleavage charge-separation reactions and the isomeric structures of b(2)(2+) ions. Three possible structures are considered for the b(2)(2+) ions: the oxazolone, diketopiperazine, and fused ring isomers. The last is formed by cleavage of the His-Gly amide bond on a pathway that is initiated by nucleophilic attack of one of the His side-chain imidazole nitrogens. Our calculations indicate the b(2)(2+) ion population is dominated by the oxazolone and/or fused ring isomers.