Development and validation of a spectral library searching method for peptide identification from MS/MS

A notable inefficiency of shotgun proteomics experiments is the repeated rediscovery of the same identifiable peptides by sequence database searching methods, which often are time-consuming and error-prone. A more precise and efficient method, in which previously observed and identified peptide MS/MS spectra are catalogued and condensed into searchable spectral libraries to allow new identifications by spectral matching, is seen as a promising alternative. To that end, an open-source, functionally complete, high-throughput and readily extensible MS/MS spectral searching tool, SpectraST, was developed. A high-quality spectral library was constructed by combining the high-confidence identifications of millions of spectra taken from various data repositories and searched using four sequence search engines. The resulting library consists of over 30,000 spectra for Saccharomyces cerevisiae. Using this library, SpectraST vastly outperforms the sequence search engine SEQUEST in terms of speed and the ability to discriminate good and bad hits. A unique advantage of SpectraST is its full integration into the popular Trans Proteomic Pipeline suite of software, which facilitates user adoption and provides important functionalities such as peptide and protein probability assignment, quantification, and data visualization. This method of spectral library searching is especially suited for targeted proteomics applications, offering superior performance to traditional sequence searching.     

Gaining knowledge from previously unexplained spectra-application of the PTM-Explorer software to detect PTM in HUPO BPP MS/MS data

A novel software tool named PTM-Explorer has been applied to LC-MS/MS datasets acquired within the Human Proteome Organisation (HUPO) Brain Proteome Project (BPP). PTM-Explorer enables automatic identification of peptide MS/MS spectra that were not explained in typical sequence database searches. The main focus was detection of PTMs, but PTM-Explorer detects also unspecific peptide cleavage, mass measurement errors, experimental modifications, amino acid substitutions, transpeptidation products and unknown mass shifts. To avoid a combinatorial problem the search is restricted to a set of selected protein sequences, which stem from previous protein identifications using a common sequence database search. Prior to application to the HUPO BPP data, PTM-Explorer was evaluated on excellently manually characterized and evaluated LC-MS/MS data sets from Alpha-A-Crystallin gel spots obtained from mouse eye lens. Besides various PTMs including phosphorylation, a wealth of experimental modifications and unspecific cleavage products were successfully detected, completing the primary structure information of the measured proteins. Our results indicate that a large amount of MS/MS spectra that currently remain unidentified in standard database searches contain valuable information that can only be elucidated using suitable software tools.     

PAS‐cal: A repetitive peptide sequence calibration standard for MALDI mass spectrometry

Mass spectrometers equipped with matrix‐assisted laser desorption/ionization (MALDI‐MS) require frequent multipoint calibration to obtain good mass accuracy over a wide mass range and across large numbers of samples. In this study, we introduce a new synthetic peptide mass calibration standard termed PAS‐cal tailored for MALDI‐MS based bottom‐up proteomics. This standard consists of 30 peptides between 8 and 37 amino acids long and each constructed to contain repetitive sequences of Pro, Ala and Ser as well as one C‐terminal arginine residue. MALDI spectra thus cover a mass range between 750 and 3200 m/z in MS mode and between 100 and 3200 m/z in MS/MS mode. Our results show that multipoint calibration of MS spectra using PAS‐cal peptides compares well to current commercial reagents for protein identification by PMF. Calibration of tandem mass spectra from LC‐MALDI experiments using the longest peptide, PAS‐cal37, resulted in smaller fragment ion mass errors, more matching fragment ions and more protein and peptide identifications compared to commercial standards, making the PAS‐cal standard generically useful for bottom‐up proteomics.     

GOAT – A simple LC‐MS/MS gradient optimization tool

Modern nano‐HPLC systems are capable of extremely precise control of solvent gradients, allowing high‐resolution separation of peptides. Most proteomics laboratories use a simple linear analytical gradient for nano‐LC‐MS/MS experiments, though recent evidence indicates that optimized non‐linear gradients result in increased peptide and protein identifications from cell lysates. In concurrent work, we examined non‐linear gradients for the analysis of samples fractionated at the peptide level, where the distribution of peptide retention times often varies by fraction. We hypothesized that greater coverage of these samples could be achieved using per‐fraction optimized gradients. We demonstrate that the optimized gradients improve the distribution of peptides throughout the analysis. Using previous generation MS instrumentation, a considerable gain in peptide and protein identifications can be realized. With current MS platforms that have faster electronics and achieve shorter duty cycle, the improvement in identifications is smaller. Our gradient optimization method has been implemented in a simple graphical tool (GOAT) that is MS‐vendor independent, does not require peptide ID input, and is freely available for non‐commercial use at http://proteomics.swmed.edu/goat/     

pepgrep: A Tool for Peptide MS/MS Pattern Matching

Typically, detection of protein sequences in collision-induced dissociation (CID) tandem MS (MS2) dataset is performed by mapping identified peptide ions back to protein sequence by using the protein database search (PDS) engine. Finding a particular peptide sequence of interest in CID MS2 records very often requires manual evaluation of the spectrum, regardless of whether the peptide-associated MS2 scan is identified by PDS algorithm or not. We have developed a compact cross-platform database-free command-line utility, pepgrep, which helps to find an MS2 fingerprint for a selected peptide sequence by pattern-matching of modelled MS2 data using Peptide-to-MS2 scoring algorithm. pepgrep can incorporate dozens of mass offsets corresponding to a variety of post-translational modifications (PTMs) into the algorithm. Decoy peptide sequences are used with the tested peptide sequence to reduce false-positive results. The engine is capable of screening an MS2 data file at a high rate when using a cluster computing environment. The matched MS2 spectrum can be displayed by using built-in graphical application programming interface (API) or optionally recorded to file. Using this algorithm, we were able to find extra peptide sequences in studied CID spectra that were missed by PDS identification. Also we found pepgrep especially useful for examining a CID of small fractions of peptides resulting from, for example, affinity purification techniques. The peptide sequences in such samples are less likely to be positively identified by using routine protein-centric algorithm implemented in PDS. The software is freely available at http://bsproteomics.essex.ac.uk:8080/data/download/pepgrep-1.4.tgz.     

LFQuant: A label‐free fast quantitative analysis tool for high‐resolution LC‐MS/MS proteomics data

Database searching based methods for label‐free quantification aim to reconstruct the peptide extracted ion chromatogram based on the identification information, which can limit the search space and thus make the data processing much faster. The random effect of the MS/MS sampling can be remedied by cross‐assignment among different runs. Here, we present a new label‐free fast quantitative analysis tool, LFQuant, for high‐resolution LC‐MS/MS proteomics data based on database searching. It is designed to accept raw data in two common formats (mzXML and Thermo RAW), and database search results from mainstream tools (MASCOT, SEQUEST, and X!Tandem), as input data. LFQuant can handle large‐scale label‐free data with fractionation such as SDS‐PAGE and 2D LC. It is easy to use and provides handy user interfaces for data loading, parameter setting, quantitative analysis, and quantitative data visualization. LFQuant was compared with two common quantification software packages, MaxQuant and IDEAL‐Q, on the replication data set and the UPS1 standard data set. The results show that LFQuant performs better than them in terms of both precision and accuracy, and consumes significantly less processing time. LFQuant is freely available under the GNU General Public License v3.0 at http://sourceforge.net/projects/lfquant/ .     

STEPS: A grid search methodology for optimized peptide identification filtering of MS/MS database search results

For bottom‐up proteomics, there are wide variety of database‐searching algorithms in use for matching peptide sequences to tandem MS spectra. Likewise, there are numerous strategies being employed to produce a confident list of peptide identifications from the different search algorithm outputs. Here we introduce a grid‐search approach for determining optimal database filtering criteria in shotgun proteomics data analyses that is easily adaptable to any search. Systematic Trial and Error Parameter Selection‐–referred to as STEPS‐–utilizes user‐defined parameter ranges to test a wide array of parameter combinations to arrive at an optimal “parameter set” for data filtering, thus maximizing confident identifications. The benefits of this approach in terms of numbers of true‐positive identifications are demonstrated using datasets derived from immunoaffinity‐depleted blood serum and a bacterial cell lysate, two common proteomics sample types.     

MassSieve: Panning MS/MS peptide data for proteins

We present MassSieve, a Java‐based platform for visualization and parsimony analysis of single and comparative LC‐MS/MS database search engine results. The success of mass spectrometric peptide sequence assignment algorithms has led to the need for a tool to merge and evaluate the increasing data set sizes that result from LC‐MS/MS‐based shotgun proteomic experiments. MassSieve supports reports from multiple search engines with differing search characteristics, which can increase peptide sequence coverage and/or identify conflicting or ambiguous spectral assignments.     

PROTEOMER: A workflow‐optimized laboratory information management system for 2‐D electrophoresis‐centered proteomics

In recent years proteomics became increasingly important to functional genomics. Although a large amount of data is generated by high throughput large‐scale techniques, a connection of these mostly heterogeneous data from different analytical platforms and of different experiments is limited. Data mining procedures and algorithms are often insufficient to extract meaningful results from large datasets and therefore limit the exploitation of the generated biological information. In our proteomic core facility, which almost exclusively focuses on 2‐DE/MS‐based proteomics, we developed a proteomic database custom tailored to our needs aiming at connecting MS protein identification information to 2‐DE derived protein expression profiles. The tools developed should not only enable an automatic evaluation of single experiments, but also link multiple 2‐DE experiments with MS‐data on different levels and thereby helping to create a comprehensive network of our proteomics data. Therefore the key feature of our “PROTEOMER” database is its high cross‐referencing capacity, enabling integration of a wide range of experimental data. To illustrate the workflow and utility of the system, two practical examples are provided to demonstrate that proper data cross‐referencing can transform information into biological knowledge.     

Characterization of a high field Orbitrap mass spectrometer for proteome analysis

The field of proteomics continues to be driven by improvements in analytical technology, notably in peptide separation, quantitative MS, and informatics. In this study, we have characterized a hybrid linear ion trap high field Orbitrap mass spectrometer (Orbitrap Elite) for proteomic applications. The very high resolution available on this instrument allows 95% of all peptide masses to be measured with sub‐ppm accuracy that in turn improves protein identification by database searching. We further confirm again that mass accuracy in tandem mass spectra is a valuable parameter for improving the success of protein identification. The new CID rapid scan type of the Orbitrap Elite achieves similar performance as higher energy collision induced dissociation fragmentation and both allow the identification of hundreds of proteins from as little as 0.1 ng of protein digest on column. The new instrument outperforms its predecessor the Orbitrap Velos by a considerable margin on each metric assessed that makes it a valuable and versatile tool for MS‐based proteomics.     

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.     

Data‐independent acquisition‐based SWATH‐MS for quantitative proteomics: a tutorial

Many research questions in fields such as personalized medicine, drug screens or systems biology depend on obtaining consistent and quantitatively accurate proteomics data from many samples. SWATH‐MS is a specific variant of data‐independent acquisition (DIA) methods and is emerging as a technology that combines deep proteome coverage capabilities with quantitative consistency and accuracy. In a SWATH‐MS measurement, all ionized peptides of a given sample that fall within a specified mass range are fragmented in a systematic and unbiased fashion using rather large precursor isolation windows. To analyse SWATH‐MS data, a strategy based on peptide‐centric scoring has been established, which typically requires prior knowledge about the chromatographic and mass spectrometric behaviour of peptides of interest in the form of spectral libraries and peptide query parameters. This tutorial provides guidelines on how to set up and plan a SWATH‐MS experiment, how to perform the mass spectrometric measurement and how to analyse SWATH‐MS data using peptide‐centric scoring. Furthermore, concepts on how to improve SWATH‐MS data acquisition, potential trade‐offs of parameter settings and alternative data analysis strategies are discussed.     

FragmentationAnalyzer: An open‐source tool to analyze MS/MS fragmentation data

A thorough understanding of the fragmentation processes in MS/MS can be a powerful tool in assessing the resulting peptide and protein identifications. We here present the freely available, open‐source FragmentationAnalyzer tool ( http://fragmentation‐analyzer.googlecode.com ) that makes it straightforward to analyze large MS/MS data sets for specific types of identified peptides, using a common set of peptide properties. This enables the detection of fragmentation pattern nuances related to specific instruments or due to the presence of post‐translational modifications.     

A reproducibility-based evaluation procedure for quantifying the differences between MS/MS peak intensity normalization methods

The identification of peptides and proteins from fragmentation mass spectra is a very common approach in the field of proteomics. Contemporary high-throughput peptide identification pipelines can quickly produce large quantities of MS/MS data that contain valuable knowledge about the actual physicochemical processes involved in the peptide fragmentation process, which can be extracted through extensive data mining studies. As these studies attempt to exploit the intensity information contained in the MS/MS spectra, a critical step required for a meaningful comparison of this information between MS/MS spectra is peak intensity normalization. We here describe a procedure for quantifying the efficiency of different published normalization methods in terms of the quartile coefficient of dispersion (qcod) statistic. The quartile coefficient of dispersion is applied to measure the dispersion of the peak intensities between redundant MS/MS spectra, allowing the quantification of the differences in computed peak intensity reproducibility between the different normalization methods. We demonstrate that our results are independent of the data set used in the evaluation procedure, allowing us to provide generic guidance on the choice of normalization method to apply in a certain MS/MS pipeline application.     

Validation of an LC‐MS based approach for profiling histones in chronic lymphocytic leukemia

The in vitro evaluation of histones and their PTMs has drawn substantial interest in the development of epigenetic therapies. The differential expression of histone isoforms may serve as a potential marker in the classification of diseases affected by chromatin abnormalities. In this study, protein profiling by LC and MS was used to explore differences in histone composition in primary chronic lymphocytic leukemia (CLL) cells. Extensive method validations were performed to determine the experimental variances that would impact histone relative abundance. The resulting data demonstrated that the proposed methodology was suitable for the analysis of histone profiles. In 4 normal individuals and 40 CLL patients, a significant decrease in the relative abundance of histone H2A variants (H2AFL and H2AFA/M*) was observed in primary CLL cells as compared to normal B cells. Protein identities were determined using high mass accuracy MS and shotgun proteomics.     

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.     

SwissPIT: An workflow‐based platform for analyzing tandem‐MS spectra using the Grid

The identification and characterization of peptides from MS/MS data represents a critical aspect of proteomics. It has been the subject of extensive research in bioinformatics resulting in the generation of a fair number of identification software tools. Most often, only one program with a specific and unvarying set of parameters is selected for identifying proteins. Hence, a significant proportion of the experimental spectra do not match the peptide sequences in the screened database due to inappropriate parameters or scoring schemes. The Swiss protein identification toolbox (swissPIT) project provides the scientific community with an expandable multitool platform for automated in‐depth analysis of MS data also able to handle data from high‐throughput experiments. With swissPIT many problems have been solved: The missing standards for input and output formats (A), creation of analysis workflows (B), unified result visualization (C), and simplicity of the user interface (D). Currently, swissPIT supports four different programs implementing two different search strategies to identify MS/MS spectra. Conceived to handle the calculation‐intensive needs of each of the programs, swissPIT uses the distributed resources of a Swiss‐wide computer Grid (http://www.swing‐grid.ch).