Incorporating sequence information into the scoring function: a hidden Markov model for improved peptide identification

MOTIVATION: The identification of peptides by tandem mass spectrometry (MS/MS) is a central method of proteomics research, but due to the complexity of MS/MS data and the large databases searched, the accuracy of peptide identification algorithms remains limited. To improve the accuracy of identification we applied a machine-learning approach using a hidden Markov model (HMM) to capture the complex and often subtle links between a peptide sequence and its MS/MS spectrum. Model: Our model, HMM_Score, represents ion types as HMM states and calculates the maximum joint probability for a peptide/spectrum pair using emission probabilities from three factors: the amino acids adjacent to each fragmentation site, the mass dependence of ion types and the intensity dependence of ion types. The Viterbi algorithm is used to calculate the most probable assignment between ion types in a spectrum and a peptide sequence, then a correction factor is added to account for the propensity of the model to favor longer peptides. An expectation value is calculated based on the model score to assess the significance of each peptide/spectrum match. RESULTS: We trained and tested HMM_Score on three data sets generated by two different mass spectrometer types. For a reference data set recently reported in the literature and validated using seven identification algorithms, HMM_Score produced 43% more positive identification results at a 1% false positive rate than the best of two other commonly used algorithms, Mascot and X!Tandem. HMM_Score is a highly accurate platform for peptide identification that works well for a variety of mass spectrometer and biological sample types. AVAILABILITY: The program is freely available on ProteomeCommons via an OpenSource license. See http://bioinfo.unc.edu/downloads/ for the download link.     

A data-mining scheme for identifying peptide structural motifs responsible for different MS/MS fragmentation intensity patterns

Although tandem mass spectrometry (MS/MS) has become an integral part of proteomics, intensity patterns in MS/MS spectra are rarely weighted heavily in most widely used algorithms because they are not yet fully understood. Here a knowledge mining approach is demonstrated to discover fragmentation intensity patterns and elucidate the chemical factors behind such patterns. Fragmentation intensity information from 28 330 ion trap peptide MS/MS spectra of different charge states and sequences went through unsupervised clustering using a penalized K-means algorithm. Without any prior chemistry assumptions, four clusters with distinctive fragmentation patterns were obtained. A decision tree was generated to investigate peptide sequence motif and charge state status that caused these fragmentation patterns. This data-mining scheme is generally applicable for any large data sets. It bypasses the common prior knowledge constraints and reports on the overall peptide fragmentation behavior. It improves the understanding of gas-phase peptide dissociation and provides a foundation for new or improved protein identification algorithms.     

Charge variant native mass spectrometry benefits mass precision and dynamic range of monoclonal antibody intact mass analysis

ABSTRACT

The preponderance and diversity of charge variants in therapeutic monoclonal antibodies has implications for antibody efficacy and degradation. Understanding the extent and impact of minor antibody variants is of great interest, and it is also a critical regulatory requirement. Traditionally, a combination of approaches is used to characterize antibody charge heterogeneity, including ion exchange chromatography and independent mass spectrometric variant site mapping after proteolytic digestion. Here, we describe charge variant native mass spectrometry (CVMS), an integrated native ion exchange mass spectrometry-based charge variant analytical approach that delivers detailed molecular information in a single, semi-automated analysis. We utilized pure volatile salt mobile phases over a pH gradient that effectively separated variants based on minimal differences in isoelectric point. Characterization of variants such as deamidation, which are traditionally unattainable by intact mass due to their minimal molecular weight differences, were measured unambiguously by mass and retention time to allow confident MS1 identification. We demonstrate that efficient chromatographic separation allows introduction of the purified forms of the charge variant isoforms into the Orbitrap mass spectrometer. Our CVMS method allows confident assignment of intact monoclonal antibody isoforms of similar mass and relative abundance measurements across three orders of magnitude dynamic range.     

Increasing information from shotgun proteomic data by accounting for misassigned precursor ion masses

Although mass spectrometers are capable of providing high mass accuracy data, assignment of true monoisotopic precursor ion mass is complicated during data-dependent ion selection for LC-MS/MS analysis of complex mixtures. The complication arises when chromatographic peak widths for a given analyte exceed the time required to acquire a precursor ion mass spectrum. The result is that many measured monoisotopic masses are misassigned due to calculation from a single mass spectrum with poor ion statistics based on only a fraction of the total available ions for a given analyte. Such data in turn produces errors in automated database searches, where precursor m/z value is one search parameter. We propose here a postacquisition approach to correct misassigned monoisotopic m/z values that involves peak detection over the entire elution profile and correction of the precursor ion monoisotopic mass. As a result of using this approach to reprocess shotgun proteomic data we increased peptide sequence assignments by 10% while reducing the estimated false positive ratio from 1 to 0.2%. We also show that 4% of the salvaged identifications may be accounted for by correction of mixed tandem mass spectra resulting from fragmentation of multiple peptides simultaneously, a situation which we refer to as accidental CID.     

Adaptive discriminant function analysis and reranking of MS/MS database search results for improved peptide identification in shotgun proteomics

Robust statistical validation of peptide identifications obtained by tandem mass spectrometry and sequence database searching is an important task in shotgun proteomics. PeptideProphet is a commonly used computational tool that computes confidence measures for peptide identifications. In this paper, we investigate several limitations of the PeptideProphet modeling approach, including the use of fixed coefficients in computing the discriminant search score and selection of the top scoring peptide assignment per spectrum only. To address these limitations, we describe an adaptive method in which a new discriminant function is learned from the data in an iterative fashion. We extend the modeling framework to go beyond the top scoring peptide assignment per spectrum. We also investigate the effect of clustering the spectra according to their spectrum quality score followed by cluster-specific mixture modeling. The analysis is carried out using data acquired from a mixture of purified proteins on four different types of mass spectrometers, as well as using a complex human serum data set. A special emphasis is placed on the analysis of data generated on high mass accuracy instruments.     

Improved de novo peptide sequencing using LC retention time information

Background

Liquid chromatography combined with tandem mass spectrometry is an important tool in proteomics for peptide identification. Liquid chromatography temporally separates the peptides in a sample. The peptides that elute one after another are analyzed via tandem mass spectrometry by measuring the mass-to-charge ratio of a peptide and its fragments. De novo peptide sequencing is the problem of reconstructing the amino acid sequences of a peptide from this measurement data. Past de novo sequencing algorithms solely consider the mass spectrum of the fragments for reconstructing a sequence.

Results

We propose to additionally exploit the information obtained from liquid chromatography. We study the problem of computing a sequence that is not only in accordance with the experimental mass spectrum, but also with the chromatographic retention time. We consider three models for predicting the retention time and develop algorithms for de novo sequencing for each model.

Conclusions

Based on an evaluation for two prediction models on experimental data from synthesized peptides we conclude that the identification rates are improved by exploiting the chromatographic information. In our evaluation, we compare our algorithms using the retention time information with algorithms using the same scoring model, but not the retention time.

     

QuickMod: A tool for open modification spectrum library searches

MS2 library spectra are rich in reproducible information about peptide fragmentation patterns compared to theoretical spectra modeled by a sequence search tool. So far, spectrum library searches are mostly applied to detect peptides as they are present in the library. However, they also allow finding modified variants of the library peptides if the search is done with a large precursor mass window and an adapted Spectrum-Spectrum Match (SSM) scoring algorithm. We perform a thorough evaluation on the use of library spectra as opposed to theoretical peptide spectra for the identification of PTMs, analyzing spectra of a well-annotated modification-rich test data set compiled from public data repositories. These initial studies motivate the development of our modification tolerant spectrum library search tool QuickMod, designed to identify modified variants of the peptides listed in the spectrum library without any prior input from the user estimating the modifications present in the sample. We built the search algorithm of QuickMod after carefully testing different SSM similarity scores. The final spectrum scoring scheme uses a support vector machine (SVM) on a selection of scoring features to classify correct and incorrect SSM. After identification of a list of modified peptides at a given False Discovery Rate (FDR), the modifications need to be positioned on the peptide sequence. We present a rapid modification site assignment algorithm and evaluate its positioning accuracy. Finally, we demonstrate that QuickMod performs favorably in terms of speed and identification rate when compared to other software solutions for PTM analysis.     

Low energy peptide fragmentations in an ESI-Q-Tof type mass spectrometer

Efficient peptide sequencing relies on both high quality MS/MS data acquisition and exhaustive knowledge of gas-phase dissociation mechanisms. We report our contribution to the elaboration of more comprehensive fragmentation models required for efficient automated MS/MS spectra interpretation. Following a statistical approach, various peptides (296 sequences of variable compositions and lengths) were prepared and subjected to low-energy collision-induced dissociations (CID) in an electrospray hybrid instrument (ESI-Q-q-Tof type mass spectrometer) that has retained relatively limited attention so far. Besides, our studies were focused on low molecular weight singly charged peptides that often failed to be identified by sequencing algorithms. Only half of the studied compounds showed charge directed dissociations in accordance with the mobile proton model producing fragment ions directly related to the primary sequence. For the peptides that did not exhibit the expected fragment ion series, alternative dissociation behaviors issued from complex rearrangements were evidenced.     

Effect of iTRAQ Labeling on the Relative Abundance of Peptide Fragment Ions Produced by MALDI-MS/MS

The identification of proteins in proteomics experiments is usually based on mass information derived from tandem mass spectrometry data. To improve the performance of the identification algorithms, additional information available in the fragment peak intensity patterns has been shown to be useful. In this study, we consider the effect of iTRAQ labeling on the fragment peak intensity patterns of singly charged peptides from MALDI tandem MS data. The presence of an iTRAQ-modified basic group on the N-terminus leads to a more pronounced set of b-ion peaks and distinct changes in the abundance of specific peptide types. We performed a simple intensity prediction by using a decision-tree machine learning approach and were able to show that the relative ion abundance in a spectrum can be correctly predicted and distinguished from closely related sequences. This information will be useful for the development of improved method-specific intensity-based protein identification algorithms.     

The Effect of Peptide Identification Search Algorithms on MS2-Based Label-Free Protein Quantification

Abstract Several approaches exist for the quantification of proteins in complex samples processed by liquid chromatography-mass spectrometry followed by fragmentation analysis (MS2). One of these approaches is label-free MS2-based quantification, which takes advantage of the information computed from MS2 spectrum observations to estimate the abundance of a protein in a sample. As a first step in this approach, fragmentation spectra are typically matched to the peptides that generated them by a search algorithm. Because different search algorithms identify overlapping but non-identical sets of peptides, here we investigate whether these differences in peptide identification have an impact on the quantification of the proteins in the sample. We therefore evaluated the effect of using different search algorithms by examining the reproducibility of protein quantification in technical repeat measurements of the same sample. From our results, it is clear that a search engine effect does exist for MS2-based label-free protein quantification methods. As a general conclusion, it is recommended to address the overall possibility of search engine-induced bias in the protein quantification results of label-free MS2-based methods by performing the analysis with two or more distinct search engines.     

Modeling and characterization of multi-charge mass spectra for peptide sequencing

Peptide sequencing using tandem mass spectrometry data is an important and challenging problem in proteomics. We address the problem of peptide sequencing for multi-charge spectra. Most peptide sequencing algorithms currently consider only charge one or two ions even for higher-charge spectra. We give a characterization of multi-charge spectra by generalizing existing models. Using our models, we analyzed spectra from Global Proteome Machine (GPM) [Craig R, Cortens JP, Beavis RC, J Proteome Res 3:1234-1242, 2004.] (with charges 1-5), Institute for Systems Biology (ISB) [Keller A, Purvine S, Nesvizhskii AI, Stolyar S, Goodlett DR, Kolker E, OMICS 6:207-212, 2002.] and Orbitrap (both with charges 1-3). Our analysis for the GPM dataset shows that higher charge peaks contribute significantly to prediction of the complete peptide. They also help to explain why existing algorithms do not perform well on multi-charge spectra. Based on these analyses, we claim that peptide sequencing algorithms can achieve higher sensitivity results if they also consider higher charge ions. We verify this claim by proposing a de novo sequencing algorithm called the greedy best strong tag (GBST) algorithm that is simple but considers higher charge ions based on our new model. Evaluation on multi-charge spectra shows that our simple GBST algorithm outperforms Lutefisk and PepNovo, especially for the GPM spectra of charge three or more.     

Implementation and application of a versatile clustering tool for tandem mass spectrometry data

High-throughput proteomics experiments typically generate large amounts of peptide fragmentation mass spectra during a single experiment. There is often a substantial amount of redundant fragmentation of the same precursors among these spectra, which is usually considered a nuisance. We here discuss the potential of clustering and merging redundant spectra to turn this redundancy into a useful property of the dataset. To this end, we have created the first general-purpose, freely available open-source software application for clustering and merging MS/MS spectra. The application also introduces a novel approach to calculating the similarity of fragmentation mass spectra that takes into account the increased precision of modern mass spectrometers, and we suggest a simple but effective improvement to single-linkage clustering. The application and the novel algorithms are applied to several real-life proteomic datasets and the results are discussed. An analysis of the influence of the different algorithms available and their parameters is given, as well as a number of important applications of the overall approach.     

Peptide identification using vectors of small fragment ions

Traditionally, peptide identification using fragmentation spectra relies on extracting the maximum amount of information from spectra. Using different combinations of small ion masses, we show that identifying a small number of fragment ions in a spectrum is sufficient for peptide identification. We consider y2-, y3-, b2-, and b3-ions and find the combination of b2-y2 to be sufficient for many peptides. Adding either the y3- or the b3-ion increases specificity and allows reliable peptide identification in the human proteome. Fragmentation spectra and peptides are represented as n-dimensional vectors, where n is given by the number of fragment ions considered, and the peptide mass. The identification score is given by the Euclidian distance between the spectra and the matching peptide in n-dimensional space. We show that this approach, using minimal information, allows for precise and fast peptide identification.     

Multi-spectra peptide sequencing and its applications to multistage mass spectrometry

Despite a recent surge of interest in database-independent peptide identifications, accurate de novo peptide sequencing remains an elusive goal. While the recently introduced spectral network approach resulted in accurate peptide sequencing in low-complexity samples, its success depends on the chance of presence of spectra from overlapping peptides. On the other hand, while multistage mass spectrometry (collecting multiple MS 3 spectra from each MS 2 spectrum) can be applied to all spectra in a complex sample, there are currently no software tools for de novo peptide sequencing by multistage mass spectrometry. We describe a rigorous probabilistic framework for analyzing spectra of overlapping peptides and show how to apply it for multistage mass spectrometry. Our software results in both accurate de novo peptide sequencing from multistage mass spectra (despite the inferior quality of MS 3 spectra) and improved interpretation of spectral networks. We further study the problem of de novo peptide sequencing with accurate parent mass (but inaccurate fragment masses), the protocol that may soon become the dominant mode of spectral acquisition. Most existing peptide sequencing algorithms (based on the spectrum graph approach) do not track the accurate parent mass and are thus not equipped for solving this problem. We describe a de novo peptide sequencing algorithm aimed at this experimental protocol and show that it improves the sequencing accuracy on both tandem and multistage mass spectrometry.     

Evaluation of algorithms for protein identification from sequence databases using mass spectrometry data

In this work, the commonly used algorithms for mass spectrometry based protein identification, Mascot, MS-Fit, ProFound and SEQUEST, were studied in respect to the selectivity and sensitivity of their searches. The influence of various search parameters were also investigated. Approximately 6600 searches were performed using different search engines with several search parameters to establish a statistical basis. The applied mass spectrometric data set was chosen from a current proteome study. The huge amount of data could only be handled with computational assistance. We present a software solution for fully automated triggering of several peptide mass fingerprinting (PMF) and peptide fragmentation fingerprinting (PFF) algorithms. The development of this high-throughput method made an intensive evaluation based on data acquired in a typical proteome project possible. Previous evaluations of PMF and PFF algorithms were mainly based on simulations.     

Learning score function parameters for improved spectrum identification in tandem mass spectrometry experiments

The identification of proteins from spectra derived from a tandem mass spectrometry experiment involves several challenges: matching each observed spectrum to a peptide sequence, ranking the resulting collection of peptide-spectrum matches, assigning statistical confidence estimates to the matches, and identifying the proteins. The present work addresses algorithms to rank peptide-spectrum matches. Many of these algorithms, such as PeptideProphet, IDPicker, or Q-ranker, follow a similar methodology that includes representing peptide-spectrum matches as feature vectors and using optimization techniques to rank them. We propose a richer and more flexible feature set representation that is based on the parametrization of the SEQUEST XCorr score and that can be used by all of these algorithms. This extended feature set allows a more effective ranking of the peptide-spectrum matches based on the target-decoy strategy, in comparison to a baseline feature set devoid of these XCorr-based features. Ranking using the extended feature set gives 10-40% improvement in the number of distinct peptide identifications relative to a range of q-value thresholds. While this work is inspired by the model of the theoretical spectrum and the similarity measure between spectra used specifically by SEQUEST, the method itself can be applied to the output of any database search. Further, our approach can be trivially extended beyond XCorr to any linear operator that can serve as similarity score between experimental spectra and peptide sequences.     

Protein labeling by iTRAQ: a new tool for quantitative mass spectrometry in proteome research

A novel, MS-based approach for the relative quantification of proteins, relying on the derivatization of primary amino groups in intact proteins using isobaric tag for relative and absolute quantitation (iTRAQ) is presented. Due to the isobaric mass design of the iTRAQ reagents, differentially labeled proteins do not differ in mass; accordingly, their corresponding proteolytic peptides appear as single peaks in MS scans. Because quantitative information is provided by isotope-encoded reporter ions that can only be observed in MS/MS spectra, we analyzed the fragmentation behavior of ESI and MALDI ions of peptides generated from iTRAQ-labeled proteins using a TOF/TOF and/or a QTOF instrument. We observed efficient liberation of reporter ions for singly protonated peptides at low-energy collision conditions. In contrast, increased collision energies were required to liberate the iTRAQ label from lysine side chains of doubly charged peptides and, thus, to observe reporter ions suitable for relative quantification of proteins with high accuracy. We then developed a quantitative strategy that comprises labeling of intact proteins by iTRAQ followed by gel electrophoresis and peptide MS/MS analyses. As proof of principle, mixtures of five different proteins in various concentration ratios were quantified, demonstrating the general applicability of the approach presented here to quantitative MS-based proteomics.     

The utility of ETD mass spectrometry in proteomic analysis

Mass spectrometry has played an integral role in the identification of proteins and their post-translational modifications (PTM). However, analysis of some PTMs, such as phosphorylation, sulfonation, and glycosylation, is difficult with collision-activated dissociation (CAD) since the modification is labile and preferentially lost over peptide backbone fragmentation, resulting in little to no peptide sequence information. The presence of multiple basic residues also makes peptides exceptionally difficult to sequence by conventional CAD mass spectrometry. Here we review the utility of electron transfer dissociation (ETD) mass spectrometry for sequence analysis of post-translationally modified and/or highly basic peptides. Phosphorylated, sulfonated, glycosylated, nitrosylated, disulfide bonded, methylated, acetylated, and highly basic peptides have been analyzed by CAD and ETD mass spectrometry. CAD fragmentation typically produced spectra showing limited peptide backbone fragmentation. However, when these peptides were fragmented using ETD, peptide backbone fragmentation produced a complete or almost complete series of ions and thus extensive peptide sequence information. In addition, labile PTMs remained intact. These examples illustrate the utility of ETD as an advantageous tool in proteomic research by readily identifying peptides resistant to analysis by CAD. A further benefit is the ability to analyze larger, non-tryptic peptides, allowing for the detection of multiple PTMs within the context of one another.     

Improved classification of mass spectrometry database search results using newer machine learning approaches

Manual analysis of mass spectrometry data is a current bottleneck in high throughput proteomics. In particular, the need to manually validate the results of mass spectrometry database searching algorithms can be prohibitively time-consuming. Development of software tools that attempt to quantify the confidence in the assignment of a protein or peptide identity to a mass spectrum is an area of active interest. We sought to extend work in this area by investigating the potential of recent machine learning algorithms to improve the accuracy of these approaches and as a flexible framework for accommodating new data features. Specifically we demonstrated the ability of boosting and random forest approaches to improve the discrimination of true hits from false positive identifications in the results of mass spectrometry database search engines compared with thresholding and other machine learning approaches. We accommodated additional attributes obtainable from database search results, including a factor addressing proton mobility. Performance was evaluated using publically available electrospray data and a new collection of MALDI data generated from purified human reference proteins.     

Combinatorial approaches for mass spectra recalibration

Mass spectrometry has become one of the most popular analysis techniques in Proteomics and Systems Biology. With the creation of larger datasets, the automated recalibration of mass spectra becomes important to ensure that every peak in the sample spectrum is correctly assigned to some peptide and protein. Algorithms for recalibrating mass spectra have to be robust with respect to wrongly assigned peaks, as well as efficient due to the amount of mass spectrometry data. The recalibration of mass spectra leads us to the problem of finding an optimal matching between mass spectra under measurement errors. We have developed two deterministic methods that allow robust computation of such a matching: The first approach uses a computational geometry interpretation of the problem, and tries to find two parallel lines with constant distance that stab a maximal number of points in the plane. The second approach is based on finding a maximal common approximate subsequence, and improves existing algorithms by one order of magnitude exploiting the sequential nature of the matching problem. We compare our results to a computational geometry algorithm using a topological line-sweep.