Peptide-spectrum Matches Research Articles

A multi-species benchmark for training and validating mass spectrometry proteomics machine learning models

Training machine learning models for tasks such as de novo sequencing or spectral clustering requires large collections of confidently identified spectra. Here we describe a dataset of 2.8 million high-confidence peptide-spectrum matches derived from nine different species. The dataset is based on a previously described benchmark but has been re-processed to ensure consistent data quality and enforce separation of training and test peptides.

NovoBoard: A Comprehensive Framework for Evaluating the False Discovery Rate and Accuracy of De Novo Peptide Sequencing

De novo peptide sequencing is one of the most fundamental research areas in mass spectrometry-based proteomics. Many methods have often been evaluated using a couple of simple metrics that do not fully reflect their overall performance. Moreover, there has not been an established method to estimate the false discovery rate (FDR) of de novo peptide-spectrum matches. Here we propose NovoBoard, a comprehensive framework to evaluate the performance of de novo peptide-sequencing methods. The framework consists of diverse benchmark datasets (including tryptic, nontryptic, immunopeptidomics, and different species) and a standard set of accuracy metrics to evaluate the fragment ions, amino acids, and peptides of the de novo results. More importantly, a new approach is designed to evaluate de novo peptide-sequencing methods on target-decoy spectra and to estimate and validate their FDRs. Our FDR estimation provides valuable information to assess the reliability of new peptides identified by de novo sequencing tools, especially when no ground-truth information is available to evaluate their accuracy. The FDR estimation can also be used to evaluate the capability of de novo peptide sequencing tools to distinguish between de novo peptide-spectrum matches and random matches. Our results thoroughly reveal the strengths and weaknesses of different de novo peptide-sequencing methods and how their performances depend on specific applications and the types of data.

Dear-PSM: A deep learning-based peptide search engine enables full database search for proteomics.

Peptide spectrum matching is the process of linking mass spectrometry data with peptide sequences. An experimental spectrum can match thousands of candidate peptides with variable modifications leading to an exponential increase in candidates. Completing the search within a limited time is a key challenge. Traditional searches expedite the process by restricting peptide mass errors and variable modifications, but this limits interpretive capability. To address this challenge, we propose Dear-PSM, a peptide search engine that supports full database searching. Dear-PSM does not restrict peptide mass errors, matching each spectrum to all peptides in the database and increasing the number of variable modifications per peptide from the conventional 3-20. Leveraging inverted index technology, Dear-PSM creates a high-performance index table of experimental spectra and utilizes deep learning algorithms for peptide validation. Through these techniques, Dear-PSM achieves a speed breakthrough 7 times faster than mainstream search engines on a regular desktop computer, with a remarkable 240-fold reduction in memory consumption. Benchmark test results demonstrate that Dear-PSM, in full database search mode, can reproduce over 90% of the results obtained by mainstream search engines when handling complex mass spectrometry data collected from different species using various instruments. Furthermore, it uncovers a substantial number of new peptides and proteins. Dear-PSM has been publicly released on the GitHub repository https://github.com/jianweishuai/Dear-PSM.

NeoMS: Mass Spectrometry-based Method for Uncovering Mutated MHC-I Neoantigens.

Major Histocompatibility Complex (MHC) molecules play a critical role in the immune system by presenting peptides on the cell surface for recognition by T-cells. Tumor cells often produce MHC peptides with amino acid mutations, known as neoantigens, which evade T-cell recognition, leading to rapid tumor growth. In immunotherapies such as TCR-T and CAR-T, identifying these mutated MHC peptide sequences is crucial. Current mass spectrometry-based peptide identification methods primarily rely on database searching, which fails to detect mutated peptides not present in human databases. In this paper, we propose a novel workflow called NeoMS, designed to efficiently identify both non-mutated and mutated MHC-I peptides from mass spectrometry data. NeoMS utilizes a tagging algorithm to generate an expanded sequence database that includes potential mutated proteins for each sample. Furthermore, it employs a machine learning-based scoring function for each peptide-spectrum match (PSM) to maximize search sensitivity. Finally, a rigorous target-decoy approach is implemented to control the false discovery rates (FDR) of the peptides with and without mutations separately. Experimental results for regular peptides demonstrate that NeoMS outperforms four benchmark methods. For mutated peptides, NeoMS successfully identifies hundreds of high-quality mutated peptides in a melanoma-associated sample, with their validity confirmed by further studies.

An algorithm for decoy-free false discovery rate estimation in XL-MS/MS proteomics.

Cross-linking tandem mass spectrometry (XL-MS/MS) is an established analytical platform used to determine distance constraints between residues within a protein or from physically interacting proteins, thus improving our understanding of protein structure and function. To aid biological discovery with XL-MS/MS, it is essential that pairs of chemically linked peptides be accurately identified, a process that requires: (i) database search, that creates a ranked list of candidate peptide pairs for each experimental spectrum and (ii) false discovery rate (FDR) estimation, that determines the probability of a false match in a group of top-ranked peptide pairs with scores above a given threshold. Currently, the only available FDR estimation mechanism in XL-MS/MS is the target-decoy approach (TDA). However, despite its simplicity, TDA has both theoretical and practical limitations that impact the estimation accuracy and increase run time over potential decoy-free approaches (DFAs). We introduce a novel decoy-free framework for FDR estimation in XL-MS/MS. Our approach relies on multi-sample mixtures of skew normal distributions, where the latent components correspond to the scores of correct peptide pairs (both peptides identified correctly), partially incorrect peptide pairs (one peptide identified correctly, the other incorrectly), and incorrect peptide pairs (both peptides identified incorrectly). To learn these components, we exploit the score distributions of first- and second-ranked peptide-spectrum matches for each experimental spectrum and subsequently estimate FDR using a novel expectation-maximization algorithm with constraints. We evaluate the method on ten datasets and provide evidence that the proposed DFA is theoretically sound and a viable alternative to TDA owing to its good performance in terms of accuracy, variance of estimation, and run time. https://github.com/shawn-peng/xlms.

Comparative analysis of plasma affinity depletion methods: Impact on protein composition and phosphopeptide abundance in human plasma.

Affinity-based protein depletion and TiO2 enrichment methods play a crucial role in detection of low-abundant proteins and phosphopeptides enrichment, respectively. Here, we assessed the effectiveness of HSA/IgG (HU2) and Human 7 (HU7) depletion methods and their impact on phosphopeptides coverage through comparative proteome analysis, utilizing in-solution digestion and nano-LC-Orbitrap mass spectrometry (MS). Our results demonstrated that both HU2 and HU7 affinity depletion significantly decreased high-abundant proteins by 1.5-7.8-fold (p<0.001). A total of 1491 proteins were identified, with 48 proteins showing significant expression in the depleted groups. Notably, cadherin-13, neutrophil defensin 1, APM1, and desmoplakin variant protein were exclusively detected in the HU2/HU7-depleted groups. Furthermore, study on effect of depletion on phosphopeptides revealed an increase in tandem MS spectral counts with notable decrease (∼50%) in peptide spectrum matching in depleted groups, which was attributed to significant reduction in protein counts. Our post translation modification workflow for phosphoproteomics detected 42 phosphorylated peptides, corresponding to 12 phosphoproteins with unique peptide match ≥2 (high false discovery rates confidence). Among them, 10 phosphorylated proteins are highly expressed in depleted groups. Overall, these findings offer valuable insights in selection of protein depletion methods for comprehensive plasma proteomics analysis.

High-Throughput Proteomics Enabled by a Fully Automated Dual-Trap and Dual-Column LC-MS.

This Technical Note describes a dual-column liquid chromatography system coupled to mass spectrometry (LC-MS) for high-throughput bottom-up proteomic analysis. This system made full use of two 2-position 10-port valves and a binary pump with an integrated loading pump of a commercial LC instrument to provide successive operation of two parallel subsystems. Each subsystem consisted of a set of trap columns and an analytical column. A T-junction union was used to split the mobile phase from the loading pump into two parts. This allowed one set of columns to be washed and equilibrated, followed by the injection of the next sample, while the previous sample was eluting and being analyzed on the other set of columns, thereby greatly increasing the analysis throughput. This approach showed high reproducibility for the analysis of HeLa tryptic digests with average relative standard deviation (RSD) values of 1.75%, 6.90%, and 5.19% for the identification number of proteins, peptides, and peptide-spectrum matches (PSMs), respectively, across 10 consecutive runs. The capacity for peptide and protein identification, as well as proteome depth, of the dual-column LC system was comparable to a conventional single-column system. Due to its simple equipment requirements and set up process, this method should be highly accessible for other laboratories.

APIR: Aggregating Universal Proteomics Database Search Algorithms for Peptide Identification with FDR Control.

Advances in mass spectrometry (MS) have enabled high-throughput analysis of proteomes in biological systems. The state-of-the-art MS data analysis relies on database search algorithms to quantify proteins by identifying peptide-spectrum matches (PSMs), which convert mass spectra to peptide sequences. Different database search algorithms use distinct search strategies and thus may identify unique PSMs. However, no existing approaches can aggregate all user-specified database search algorithms with a guaranteed increase in the number of identified peptides and a control on the false discovery rate (FDR). To fill in this gap, we proposed a statistical framework, Aggregation of Peptide Identification Results (APIR), that is universally compatible with all database search algorithms. Notably, under an FDR threshold, APIR is guaranteed to identify at least as many, if not more, peptides as individual database search algorithms do. Evaluation of APIR on a complex proteomics standard dataset showed that APIR outpowers individual database search algorithms and empirically controls the FDR. Real data studies showed that APIR can identify disease-related proteins and post-translational modifications missed by some individual database search algorithms. The APIR framework is easily extendable to aggregating discoveries made by multiple algorithms in other high-throughput biomedical data analysis, e.g., differential gene expression analysis on RNA sequencing data. The APIR R package is available at https://github.com/yiling0210/APIR.

Rescoring Peptide Spectrum Matches: Boosting Proteomics Performance by Integrating Peptide Property Predictors Into Peptide Identification

Rescoring of peptide spectrum matches originating from database search engines enabled by peptide property predictors is exceeding the performance of peptide identification from traditional database search engines. In contrast to the peptide spectrum match scores calculated by traditional database search engines, rescoring peptide spectrum matches generates scores based on comparing observed and predicted peptide properties, such as fragment ion intensities and retention times. These newly generated scores enable a more efficient discrimination between correct and incorrect peptide spectrum matches. This approach was shown to lead to substantial improvements in the number of confidently identified peptides, facilitating the analysis of challenging datasets in various fields such as immunopeptidomics, metaproteomics, proteogenomics, and single-cell proteomics. In this review, we summarize the key elements leading up to the recent introduction of multiple data-driven rescoring pipelines. We provide an overview of relevant post-processing rescoring tools, introduce prominent data-driven rescoring pipelines for various applications, and highlight limitations, opportunities, and future perspectives of this approach and its impact on mass spectrometry-based proteomics.

Fragment ion intensity prediction improves the identification rate of non-tryptic peptides in timsTOF

Immunopeptidomics is crucial for immunotherapy and vaccine development. Because the generation of immunopeptides from their parent proteins does not adhere to clear-cut rules, rather than being able to use known digestion patterns, every possible protein subsequence within human leukocyte antigen (HLA) class-specific length restrictions needs to be considered during sequence database searching. This leads to an inflation of the search space and results in lower spectrum annotation rates. Peptide-spectrum match (PSM) rescoring is a powerful enhancement of standard searching that boosts the spectrum annotation performance. We analyze 302,105 unique synthesized non-tryptic peptides from the ProteomeTools project on a timsTOF-Pro to generate a ground-truth dataset containing 93,227 MS/MS spectra of 74,847 unique peptides, that is used to fine-tune the deep learning-based fragment ion intensity prediction model Prosit. We demonstrate up to 3-fold improvement in the identification of immunopeptides, as well as increased detection of immunopeptides from low input samples.

Comprehensive shotgun proteomic characterization and virulence factors of seafood spoilage bacteria

This article summarizes the characterization, by shotgun proteomics, of 11 bacterial strains identified as responsible for seafood spoilage. A total of 4455 peptide spectrum matches, corresponding to 4299 peptides and 3817 proteins were identified. Analyses of data determined the functional pathways they are involved in. The proteins identified were integrated into a protein-protein network that involves 371 nodes and 3016 edges. Those proteins are implicated in energy pathways, peptidoglycan biosynthesis, spermidine/putrescine metabolism. An additional 773 peptides were characterized as virulence factors, that participates in bacterial pathogenesis; while 14 peptides were defined as biomarkers, as they can be used to differentiate the bacterial species present. This report represents the most extensive proteomic repository available in the field of seafood spoilage bacteria; the data substantially advances the understanding of seafood decay, as well as provides fundamental bases for the recognition of the bacteria existent in seafood that cause spoilage during food processing/storage.

MS2Rescore 3.0 Is a Modular, Flexible, and User-Friendly Platform to Boost Peptide Identifications, as Showcased with MS Amanda 3.0.

Rescoring of peptide-spectrum matches (PSMs) has emerged as a standard procedure for the analysis of tandem mass spectrometry data. This emphasizes the need for software maintenance and continuous improvement for such algorithms. We introduce MS2Rescore 3.0, a versatile, modular, and user-friendly platform designed to increase peptide identifications. Researchers can install MS2Rescore across various platforms with minimal effort and benefit from a graphical user interface, a modular Python API, and extensive documentation. To showcase this new version, we connected MS2Rescore 3.0 with MS Amanda 3.0, a new release of the well-established search engine, addressing previous limitations on automatic rescoring. Among new features, MS Amanda now contains additional output columns that can be used for rescoring. The full potential of rescoring is best revealed when applied on challenging data sets. We therefore evaluated the performance of these two tools on publicly available single-cell data sets, where the number of PSMs was substantially increased, thereby demonstrating that MS2Rescore offers a powerful solution to boost peptide identifications. MS2Rescore's modular design and user-friendly interface make data-driven rescoring easily accessible, even for inexperienced users. We therefore expect the MS2Rescore to be a valuable tool for the wider proteomics community. MS2Rescore is available at https://github.com/compomics/ms2rescore.

Unraveling the Intraday Variations in the Tear Fluid Proteome.

Tear fluid is a complex and dynamic biological fluid that plays essential roles in maintaining ocular homeostasis and protecting against the external environment. Owing to the small sample volume, studying the tear proteome is challenging. However, advances in high-resolution mass spectrometry have expanded tear proteome profiling, revealing >500 unique proteins. Tears are emerging as a noninvasive source of biomarkers for both ocular and systemic diseases; nevertheless, intraday variability of proteins in tear fluid remains questionable. This study investigates intraday variations in the tear fluid proteome to identify stable proteins that could act as candidate biomarkers. Tear samples from 15 individuals at four time points (10 am, 12 pm, 2 pm, and 4 pm) were analyzed using mass spectrometry to evaluate protein variation during these intervals. Technical variation was assessed by analyzing pooled samples and was subtracted from the total variation to isolate biological variability. Owing to high technical variation, low-abundant proteins were filtered, and only 115 proteins met the criteria for further analysis. These criteria include being detected at all four time points in at least eight subjects, having a mean peptide-spectrum match count greater than 5, and having a technical variation less than 0.10. Lactotransferrin, lipocalin-1, and several immunoglobulins were among the 51 stable proteins (mean biological coefficient of variation < 0.10). Additionally, 43 proteins displayed significant slopes across the 4 time points, with 17 increasing and 26 decreasing over time. These findings contribute to the understanding of tear fluid dynamics and further expand our knowledge of the tear proteome.

PXg: Comprehensive Identification of Noncanonical MHC-I–Associated Peptides From De Novo Peptide Sequencing Using RNA-Seq Reads

Discovering noncanonical peptides has been a common application of proteogenomics. Recent studies suggest that certain noncanonical peptides, known as noncanonical major histocompatibility complex-I (MHC-I)-associated peptides (ncMAPs), that bind to MHC-I may make good immunotherapeutic targets. De novo peptide sequencing is a great way to find ncMAPs since it can detect peptide sequences from their tandem mass spectra without using any sequence databases. However, this strategy has not been widely applied for ncMAP identification because there is not a good way to estimate its false-positive rates. In order to completely and accurately identify immunopeptides using de novo peptide sequencing, we describe a unique pipeline called proteomics X genomics. In contrast to current pipelines, it makes use of genomic data, RNA-Seq abundance and sequencing quality, in addition to proteomic features to increase the sensitivity and specificity of peptide identification. We show that the peptide-spectrum match quality and genetic traits have a clear relationship, showing that they can be utilized to evaluate peptide-spectrum matches. From 10 samples, we found 24,449 canonical MHC-I–associated peptides and 956 ncMAPs by using a target-decoy competition. Three hundred eighty-seven ncMAPs and 1611 canonical MHC-I–associated peptides were new identifications that had not yet been published. We discovered 11 ncMAPs produced from a squirrel monkey retrovirus in human cell lines in addition to the two ncMAPs originating from a complementarity determining region 3 in an antibody thanks to the unrestricted search space assumed by de novo sequencing. These entirely new identifications show that proteomics X genomics can make the most of de novo peptide sequencing’s advantages and its potential use in the search for new immunotherapeutic targets.

Introducing π-HelixNovo for practical large-scale de novo peptide sequencing.

De novo peptide sequencing is a promising approach for novel peptide discovery, highlighting the performance improvements for the state-of-the-art models. The quality of mass spectra often varies due to unexpected missing of certain ions, presenting a significant challenge in de novo peptide sequencing. Here, we use a novel concept of complementary spectra to enhance ion information of the experimental spectrum and demonstrate it through conceptual and practical analyses. Afterward, we design suitable encoders to encode the experimental spectrum and the corresponding complementary spectrum and propose a de novo sequencing model $\pi$-HelixNovo based on the Transformer architecture. We first demonstrated that $\pi$-HelixNovo outperforms other state-of-the-art models using a series of comparative experiments. Then, we utilized $\pi$-HelixNovo to de novo gut metaproteome peptides for the first time. The results show $\pi$-HelixNovo increases the identification coverage and accuracy of gut metaproteome and enhances the taxonomic resolution of gut metaproteome. We finally trained a powerful $\pi$-HelixNovo utilizing a larger training dataset, and as expected, $\pi$-HelixNovo achieves unprecedented performance, even for peptide-spectrum matches with never-before-seen peptide sequences. We also use the powerful $\pi$-HelixNovo to identify antibody peptides and multi-enzyme cleavage peptides, and $\pi$-HelixNovo is highly robust in these applications. Our results demonstrate the effectivity of the complementary spectrum and take a significant step forward in de novo peptide sequencing.

The proteomic landscape of sperm surface deciphers its maturational and functional aspects in buffalo.

Buffalo is a dominant dairy animal in many agriculture-based economies. However, the poor reproductive efficiency (low conception rate) of the buffalo bulls constrains the realization of its full production potential. This in turn leads to economic and welfare issues, especially for the marginal farmers in such economies. The mammalian sperm surface proteins have been implicated in the regulation of survival and function of the spermatozoa in the female reproductive tract (FRT). Nonetheless, the lack of specific studies on buffalo sperm surface makes it difficult for researchers to explore and investigate the role of these proteins in the regulation of mechanisms associated with sperm protection, survival, and function. This study aimed to generate a buffalo sperm surface-specific proteomic fingerprint (LC-MS/MS) and to predict the functional roles of the identified proteins. The three treatments used to remove sperm surface protein viz. Elevated salt, phosphoinositide phospholipase C (PI-PLC) and in vitro capacitation led to the identification of N = 1,695 proteins (≥1 high-quality peptide-spectrum matches (PSMs), p < 0.05, and FDR<0.01). Almost half of these proteins (N = 873) were found to be involved in crucial processes relevant in the context of male fertility, e.g., spermatogenesis, sperm maturation and protection in the FRT, and gamete interaction or fertilization, amongst others. The extensive sperm-surface proteomic repertoire discovered in this study is unparalleled vis-à-vis the depth of identification of reproduction-specific cell-surface proteins and can provide a potential framework for further studies on the functional aspects of buffalo spermatozoa.

Making MS Omics Data ML-Ready: SpeCollate Protocols.

The increasing complexity and volume of mass spectrometry (MS) data have presented new challenges and opportunities for proteomics data analysis and interpretation. In this chapter, we provide a comprehensive guide to transforming MS data for machine learning (ML) training, inference, and applications. The chapter is organized into three parts. The first part describes the data analysis needed for MS-based experiments and a general introduction to our deep learning model SpeCollate-which we will use throughout the chapter for illustration. The second part of the chapter explores the transformation of MS data for inference, providing a step-by-step guide for users to deduce peptides from their MS data. This section aims to bridge the gap between data acquisition and practical applications by detailing the necessary steps for data preparation and interpretation. In the final part, we present a demonstrative example of SpeCollate, a deep learning-based peptide database search engine that overcomes the problems of simplistic simulation of theoretical spectra and heuristic scoring functions for peptide-spectrum matches by generating joint embeddings for spectra and peptides. SpeCollate is a user-friendly tool with an intuitive command-line interface to perform the search, showcasing the effectiveness of the techniques and methodologies discussed in the earlier sections and highlighting the potential of machine learning in the context of mass spectrometry data analysis. By offering a comprehensive overview of data transformation, inference, and ML model applications for mass spectrometry, this chapter aims to empower researchers and practitioners in leveraging the power of machine learning to unlock novel insights and drive innovation in the field of mass spectrometry-based omics.

Quantifying the Impact of the Peptide Identification Framework on the Results of Fast Photochemical Oxidation of Protein Analysis.

Fast Photochemical Oxidation of Proteins (FPOP) is a promising technique for studying protein structure and dynamics. The quality of insight provided by FPOP depends on the reliability of the determination of the modification site. This study investigates the performance of two search engines, Mascot and PEAKS, for the data processing of FPOP analyses. Comparison of Mascot and PEAKS of the hemoglobin--haptoglobin Bruker timsTOF data set (PXD021621) revealed greater consistency in the Mascot identification of modified peptides, with around 26% of the IDs being mutual for all three replicates, compared to approximately 22% for PEAKS. The intersection between Mascot and PEAKS results revealed a limited number (31%) of shared modified peptides. Principal Component Analysis (PCA) using the peptide-spectrum match (PSM) score, site probability, and peptide intensity was applied to evaluate the results, and the analyses revealed distinct clusters of modified peptides. Mascot showed the ability to assess confident site determination, even with lower PSM scores. However, high PSM scores from PEAKS did not guarantee a reliable determination of the modification site. Fragmentation coverage of the modification position played a crucial role in Mascot assignments, while the AScore localizations from PEAKS often become ambiguous because the software employs MS/MS merging.

Test-Time Training for Deep MS/MS Spectrum Prediction Improves Peptide Identification.

In bottom-up proteomics, peptide-spectrum matching is critical for peptide and protein identification. Recently, deep learning models have been used to predict tandem mass spectra of peptides, enabling the calculation of similarity scores between the predicted and experimental spectra for peptide-spectrum matching. These models follow the supervised learning paradigm, which trains a general model using paired peptides and spectra from standard data sets and directly employs the model on experimental data. However, this approach can lead to inaccurate predictions due to differences between the training data and the experimental data, such as sample types, enzyme specificity, and instrument calibration. To tackle this problem, we developed a test-time training paradigm that adapts the pretrained model to generate experimental data-specific models, namely, PepT3. PepT3 yields a 10-40% increase in peptide identification depending on the variability in training and experimental data. Intriguingly, when applied to a patient-derived immunopeptidomic sample, PepT3 increases the identification of tumor-specific immunopeptide candidates by 60%. Two-thirds of the newly identified candidates are predicted to bind to the patient's human leukocyte antigen isoforms. To facilitate access of the model and all the results, we have archived all the intermediate files in Zenodo.org with identifier 8231084.

Deep Learning Prediction Boosts Phosphoproteomics-Based Discoveries Through Improved Phosphopeptide Identification

Shotgun phosphoproteomics enables high-throughput analysis of phosphopeptides in biological samples. One of the primary challenges associated with this technology is the relatively low rate of phosphopeptide identification during data analysis. This limitation hampers the full realization of the potential offered by shotgun phosphoproteomics. Here we present DeepRescore2, a computational workflow that leverages deep learning-based retention time and fragment ion intensity predictions to improve phosphopeptide identification and phosphosite localization. Using a state-of-the-art computational workflow as a benchmark, DeepRescore2 increases the number of correctly identified peptide-spectrum matches by 17% in a synthetic dataset and identifies 19% to 46% more phosphopeptides in biological datasets. In a liver cancer dataset, 30% of the significantly altered phosphosites between tumor and normal tissues and 60% of the prognosis-associated phosphosites identified from DeepRescore2-processed data could not be identified based on the state-of-the-art workflow. Notably, DeepRescore2-processed data uniquely identifies EGFR hyperactivation as a new target in poor-prognosis liver cancer, which is validated experimentally. Integration of deep learning prediction in DeepRescore2 improves phosphopeptide identification and facilitates biological discoveries.