On the utility of ultrafast MS1-only proteomics in drug target discovery studies based on thermal proteome profiling method.

Advances in high-throughput high-resolution mass spectrometry and the development of thermal proteome profiling approach (TPP) have made it possible to accelerate a drug target search. Since its introduction in 2014, TPP quickly became a method of choice in chemical proteomics for identifying drug-to-protein interactions on a proteome-wide scale and mapping the pathways of these interactions, thus further elucidating the unknown mechanisms of action of a drug under study. However, the current TPP implementations based on tandem mass spectrometry (MS/MS), associated with employing lengthy peptide separation protocols and expensive labeling techniques for sample multiplexing, limit the scaling of this approach for the ever growing variety of drug-to-proteomes. A variety of ultrafast proteomics methods have been developed in the last couple of years. Among them, DirectMS1 provides MS/MS-free quantitative proteome-wide analysis in 5-min time scale, thus opening the way for sample-hungry applications, such as TPP. In this work, we demonstrate the first implementation of the TPP approach using the ultrafast proteome-wide analysis based on DirectMS1. Using a drug topotecan, which is a known topoisomerase I (TOP1) inhibitor, the feasibility of the method for identifying drug targets at the whole proteome level was demonstrated for an ovarian cancer cell line.

Proteomics-based scoring of cellular response to stimuli for improved characterization of signaling pathway activity.

Omics technologies focus on uncovering the complex nature of molecular mechanisms in cells and organisms, including biomarkers and drug targets discovery. Aiming at these tasks, we see that information extracted from omics data is still underused. In particular, characteristics of differentially regulated molecules can be combined in a single score to quantify the signaling pathway activity. Such a metric can be useful for comprehensive data interpretation to follow: (1) developing molecular responses in time; (2) potency of a drug on a certain cell culture; (3) ranking the signaling pathway activity in stimulated cells; and (4) integration of the omics data and assay-based measurements of cell viability, cytotoxicity, and proliferation. With recent advances in ultrafast mass spectrometry for quantitative proteomics allowing data collection in a few minutes, proteomics score for cellular response to stimuli can become a fast, accurate, and informative complement to bioassays. Here, we utilized an interquartile-based selection of differentially regulated features and a variety of schemes for quantifying cellular responses to come up with the quantitative metric for total cellular response and pathway activity. Validation was performed using antiproliferative and virus assays and label-free proteomics data collected for cancer cells subjected to drug stimulation.

IQMMA: Efficient MS1 Intensity Extraction Pipeline Using Multiple Feature Detection Algorithms for DDA Proteomics.

One of the key steps in data dependent acquisition (DDA) proteomics is detection of peptide isotopic clusters, also called "features", in MS1 spectra and matching them to MS/MS-based peptide identifications. A number of peptide feature detection tools became available in recent years, each relying on its own matching algorithm. Here, we provide an integrated solution, the intensity-based Quantitative Mix and Match Approach (IQMMA), which integrates a number of untargeted peptide feature detection algorithms and returns the most probable intensity values for the MS/MS-based identifications. IQMMA was tested using available proteomic data acquired for both well-characterized (ground truth) and real-world biological samples, including a mix of Yeast and E. coli digests spiked at different concentrations into the Human K562 digest used as a background, and a set of glioblastoma cell lines. Three open-source feature detection algorithms were integrated: Dinosaur, biosaur2, and OpenMS FeatureFinder. None of them was found optimal when applied individually to all the data sets employed in this work; however, their combined use in IQMMA improved efficiency of subsequent protein quantitation. The software implementing IQMMA is freely available at https://github.com/PostoenkoVI/IQMMA under Apache 2.0 license.

Massive Proteogenomic Reanalysis of Publicly Available Proteomic Datasets of Human Tissues in Search for Protein Recoding via Adenosine-to-Inosine RNA Editing.

The proteogenomic search pipeline developed in this work has been applied for reanalysis of 40 publicly available shotgun proteomic datasets from various human tissues comprising more than 8000 individual LC-MS/MS runs, of which 5442 .raw data files were processed in total. This reanalysis was focused on searching for ADAR-mediated RNA editing events, their clustering across samples of different origins, and classification. In total, 33 recoded protein sites were identified in 21 datasets. Of those, 18 sites were detected in at least two datasets, representing the core human protein editome. In agreement with prior artworks, neural and cancer tissues were found to be enriched with recoded proteins. Quantitative analysis indicated that recoding the rate of specific sites did not directly depend on the levels of ADAR enzymes or targeted proteins themselves, rather it was governed by differential and yet undescribed regulation of interaction of enzymes with mRNA. Nine recoding sites conservative between humans and rodents were validated by targeted proteomics using stable isotope standards in the murine brain cortex and cerebellum, and an additional one was validated in human cerebrospinal fluid. In addition to previous data of the same type from cancer proteomes, we provide a comprehensive catalog of recoding events caused by ADAR RNA editing in the human proteome.

A Comprehensive Evaluation of Consensus Spectrum Generation Methods in Proteomics.

Spectrum clustering is a powerful strategy to minimize redundant mass spectra by grouping them based on similarity, with the aim of forming groups of mass spectra from the same repeatedly measured analytes. Each such group of near-identical spectra can be represented by its so-called consensus spectrum for downstream processing. Although several algorithms for spectrum clustering have been adequately benchmarked and tested, the influence of the consensus spectrum generation step is rarely evaluated. Here, we present an implementation and benchmark of common consensus spectrum algorithms, including spectrum averaging, spectrum binning, the most similar spectrum, and the best-identified spectrum. We have analyzed diverse public data sets using two different clustering algorithms (spectra-cluster and MaRaCluster) to evaluate how the consensus spectrum generation procedure influences downstream peptide identification. The BEST and BIN methods were found the most reliable methods for consensus spectrum generation, including for data sets with post-translational modifications (PTM) such as phosphorylation. All source code and data of the present study are freely available on GitHub at https://github.com/statisticalbiotechnology/representative-spectra-benchmark.

DirectMS1Quant: Ultrafast Quantitative Proteomics with MS/MS-Free Mass Spectrometry.

Recently, we presented the DirectMS1 method of ultrafast proteome-wide analysis based on minute-long LC gradients and MS1-only mass spectra acquisition. Currently, the method provides the depth of human cell proteome coverage of 2500 proteins at a 1% false discovery rate (FDR) when using 5 min LC gradients and 7.3 min runtime in total. While the standard MS/MS approaches provide 4000-5000 protein identifications within a couple of hours of instrumentation time, we advocate here that the higher number of identified proteins does not always translate into better quantitation quality of the proteome analysis. To further elaborate on this issue, we performed a one-on-one comparison of quantitation results obtained using DirectMS1 with three popular MS/MS-based quantitation methods: label-free (LFQ) and tandem mass tag quantitation (TMT), both based on data-dependent acquisition (DDA) and data-independent acquisition (DIA). For comparison, we performed a series of proteome-wide analyses of well-characterized (ground truth) and biologically relevant samples, including a mix of UPS1 proteins spiked at different concentrations into an Echerichia coli digest used as a background and a set of glioblastoma cell lines. MS1-only data was analyzed using a novel quantitation workflow called DirectMS1Quant developed in this work. The results obtained in this study demonstrated comparable quantitation efficiency of 5 min DirectMS1 with both TMT and DIA methods, yet the latter two utilized a 10-20-fold longer instrumentation time.

On the Feasibility of Using an Ultra-Fast DirectMS1 Method of Proteome-Wide Analysis for Searching Drug Targets in Chemical Proteomics.

Protein quantitation in tissue cells or physiological fluids based on liquid chromatography/mass spectrometry is one of the key sources of information on the mechanisms of cell functioning during chemotherapeutic treatment. Information on significant changes in protein expression upon treatment can be obtained by chemical proteomics and requires analysis of the cellular proteomes, as well as development of experimental and bioinformatic methods for identification of the drug targets. Low throughput of whole proteome analysis based on liquid chromatography and tandem mass spectrometry is one of the main factors limiting the scale of these studies. The method of direct mass spectrometric identification of proteins, DirectMS1, is one of the approaches developed in recent years allowing ultrafast proteome-wide analyses employing minute-scale gradients for separation of proteolytic mixtures. Aim of this work was evaluation of both possibilities and limitations of the method for identification of drug targets at the level of whole proteome and for revealing cellular processes activated by the treatment. Particularly, the available literature data on chemical proteomics obtained earlier for a large set of onco-pharmaceuticals using multiplex quantitative proteome profiling were analyzed. The results obtained were further compared with the proteome-wide data acquired by the DirectMS1 method using ultrashort separation gradients to evaluate efficiency of the method in identifying known drug targets. Using ovarian cancer cell line A2780 as an example, a whole-proteome comparison of two cell lysis techniques was performed, including the freeze-thaw lysis commonly employed in chemical proteomics and the one based on ultrasonication for cell disruption, which is the widely accepted as a standard in proteomic studies. Also, the proteome-wide profiling was performed using ultrafast DirectMS1 method for A2780 cell line treated with lonidamine, followed by gene ontology analyses to evaluate capabilities of the method in revealing regulation of proteins in the cellular processes associated with drug treatment.

Boosting MS1-only Proteomics with Machine Learning Allows 2000 Protein Identifications in Single-Shot Human Proteome Analysis Using 5 min HPLC Gradient.

Proteome-wide analyses rely on tandem mass spectrometry and the extensive separation of proteolytic mixtures. This imposes considerable instrumental time consumption, which is one of the main obstacles in the broader acceptance of proteomics in biomedical and clinical research. Recently, we presented a fast proteomic method termed DirectMS1 based on ultrashort LC gradients as well as MS1-only mass spectra acquisition and data processing. The method allows significant reduction of the proteome-wide analysis time to a few minutes at the depth of quantitative proteome coverage of 1000 proteins at 1% false discovery rate (FDR). In this work, to further increase the capabilities of the DirectMS1 method, we explored the opportunities presented by the recent progress in the machine-learning area and applied the LightGBM decision tree boosting algorithm to the scoring of peptide feature matches when processing MS1 spectra. Furthermore, we integrated the peptide feature identification algorithm of DirectMS1 with the recently introduced peptide retention time prediction utility, DeepLC. Additional approaches to improve the performance of the DirectMS1 method are discussed and demonstrated, such as using FAIMS for gas-phase ion separation. As a result of all improvements to DirectMS1, we succeeded in identifying more than 2000 proteins at 1% FDR from the HeLa cell line in a 5 min gradient LC-FAIMS/MS1 analysis. The data sets generated and analyzed during the current study have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD023977.

Biosaur: An open-source Python software for liquid chromatography-mass spectrometry peptide feature detection with ion mobility support.

RATIONALE: One of the important steps in initial data processing of peptide mass spectra is the detection of peptide features in full-range mass spectra. Ion mobility offers advantages over previous methods performing this detection by providing an additional structure-specific separation dimension. However, there is a lack of open-source software that utilizes these advantages and detects peptide features in mass spectra acquired along with ion mobility data using new instruments such as timsTOF and/or FAIMS-Orbitrap. METHODS: Recently, a utility called Dinosaur was presented, which provides an efficient way for feature detection in peptide ion mass spectra. In this work we extended its functionality by developing Biosaur software to fully employ the additional information provided by ion mobility data. Biosaur was developed using the Python 3.8 programming language. RESULTS: Biosaur supports the processing of data acquired using mass spectrometers with ion mobility capabilities, specifically timsTOF and FAIMS. In addition, it processes mass spectra obtained in negative ion mode and reports cosine correlation table for peptide features which is useful for differentiation between in-source fragments and semi-tryptic peptides. CONCLUSIONS: Biosaur is a utility for detecting peptide features in liquid chromatography-mass spectra with ion mobility and negative ion supports. The software is distributed with an open-source APACHE 2.0 license and is freely available on Github: https://github.com/abdrakhimov1/Biosaur.

HUMOS: How to Understand My Orbitrap Spectrum?-An Interactive Web-Based Tool to Teach the Basics of Mass-Spectrometry-Based Proteomics.

The Orbitrap mass analyzer can provide high mass accuracy and throughput, which has significantly improved proteome research and made this type of instrumentation one of the most frequently applied in proteomics. The efficient use of Orbitrap mass spectrometers requires training. Students in the field of proteomics can benefit from a deeper understanding of the Orbitrap technology to comprehend mass spectral interpretation, troubleshooting, and judgment of experimental settings. Unfortunately, the cost of high-end mass spectrometers limits the implementation of this type of equipment in educational laboratories. Guided by these concerns, we developed an eLearning web application called HUMOS aimed to help teach Orbitrap mass spectrometry. Although a typical proteomics experiment includes the use of several different technologies, such as liquid chromatography, mass spectrometry, and bioinformatics, the learning objectives of HUMOS are focused on mass spectrometry. HUMOS models a mass spectrum of a peptide mixture, allowing us to investigate the influence of mass spectral acquisition parameters. By changing the parameters and observing the differences, students can learn more about the mass spectral resolution, duty cycle, throughput of the analysis, ion accumulation, and spectral dynamic range and get familiar with advanced spectral acquisition methods, such as BoxCar. HUMOS is an open-source software published under the Apache license; the live installation is available at http://humos.bmb.sdu.dk.

DirectMS1: MS/MS-Free Identification of 1000 Proteins of Cellular Proteomes in 5 Minutes.

Proteome characterization relies heavily on tandem mass spectrometry (MS/MS) and is thus associated with instrumentation complexity, lengthy analysis time, and limited duty cycle. It was always tempting to implement approaches that do not require MS/MS, yet they were constantly failing to achieve a meaningful depth of quantitative proteome coverage within short experimental times, which is particularly important for clinical or biomarker-discovery applications. Here, we report on the first successful attempt to develop a truly MS/MS-free method, DirectMS1, for bottom-up proteomics. The method is compared with the standard MS/MS-based data-dependent acquisition approach for proteome-wide analysis using 5 min LC gradients. Specifically, we demonstrate identification of 1 000 protein groups for a standard HeLa cell line digest. The amount of loaded sample was varied in a range from 1 to 500 ng, and the method demonstrated 10-fold higher sensitivity. Combined with the recently introduced Diffacto approach for relative protein quantification, DirectMS1 outperforms most popular MS/MS-based label-free quantitation approaches because of significantly higher protein sequence coverage.

Rpn4 and proteasome-mediated yeast resistance to ethanol includes regulation of autophagy.

Distilled spirits production using Saccharomyces cerevisiae requires understanding of the mechanisms of yeast cell response to alcohol stress. Reportedly, specific mutations in genes of the ubiquitin-proteasome system, e.g., RPN4, may result in strains exhibiting hyper-resistance to different alcohols. To study the Rpn4-dependent yeast response to short-term ethanol exposure, we performed a comparative analysis of the wild-type (WT) strain, strain with RPN4 gene deletion (rpn4-Δ), and a mutant strain with decreased proteasome activity and consequent Rpn4 accumulation due to PRE1 deregulation (YPL). The stress resistance tests demonstrated an increased sensitivity of mutant strains to ethanol compared with WT. Comparative proteomics analysis revealed significant differences in molecular responses to ethanol between these strains. GO analysis of proteins upregulated in WT showed enrichments represented by oxidative and heat responses, protein folding/unfolding, and protein degradation. Enrichment of at least one of these responses was not observed in the mutant strains. Moreover, activity of autophagy was not increased in the RPN4 deletion strain upon ethanol stress which agrees with changes in mRNA levels of ATG7 and PRB1 genes of the autophagy system. Activity of the autophagic system was clearly induced and accompanied with PRB1 overexpression in the YPL strain upon ethanol stress. We demonstrated that Rpn4 stabilization contributes to the PRB1 upregulation. CRISPR-Cas9-mediated repression of PACE-core Rpn4 binding sites in the PRB1 promoter inhibits PRB1 induction in the YPL strain upon ethanol treatment and results in YPL hypersensitivity to ethanol. Our data suggest that Rpn4 affects the autophagic system activity upon ethanol stress through the PRB1 regulation. These findings can be a basis for creating genetically modified yeast strains resistant to high levels of alcohol, being further used for fermentation in ethanol production.

Scavager: A Versatile Postsearch Validation Algorithm for Shotgun Proteomics Based on Gradient Boosting.

Shotgun proteomics workflows for database protein identification typically include a combination of search engines and postsearch validation software based mostly on machine learning algorithms. Here, a new postsearch validation tool called Scavager employing CatBoost, an open-source gradient boosting library, which shows improved efficiency compared with the other popular algorithms, such as Percolator, PeptideProphet, and Q-ranker, is presented. The comparison is done using multiple data sets and search engines, including MSGF+, MSFragger, X!Tandem, Comet, and recently introduced IdentiPy. Implemented in Python programming language, Scavager is open-source and freely available at https://bitbucket.org/markmipt/scavager.

Validation of Peptide Identification Results in Proteomics Using Amino Acid Counting.

The efficiency of proteome analysis depends strongly on the configuration parameters of the search engine. One of the murkiest and nontrivial among them is the list of amino acid modifications included for the search. Here, an approach called AA_stat is presented for uncovering the unexpected modifications of amino acid residues in the protein sequences, as well as possible artifacts of data acquisition or processing, in the results of proteome analyses. The approach is based on comparing the amino acid frequencies of different mass shifts observed using the open search method introduced recently. In this work, the proposed approach is applied to publicly available proteomic data is applied and its feasibility for discovering unaccounted modifications or possible pitfalls of the identification workflow is demonstrated. The algorithm is implemented in Python as an open-source command-line tool available at https://bitbucket.org/J_Bale/aa_stat/.

IdentiPy: An Extensible Search Engine for Protein Identification in Shotgun Proteomics.

We present an open-source, extensible search engine for shotgun proteomics. Implemented in Python programming language, IdentiPy shows competitive processing speed and sensitivity compared with the state-of-the-art search engines. It is equipped with a user-friendly web interface, IdentiPy Server, enabling the use of a single server installation accessed from multiple workstations. Using a simplified version of X!Tandem scoring algorithm and its novel "autotune" feature, IdentiPy outperforms the popular alternatives on high-resolution data sets. Autotune adjusts the search parameters for the particular data set, resulting in improved search efficiency and simplifying the user experience. IdentiPy with the autotune feature shows higher sensitivity compared with the evaluated search engines. IdentiPy Server has built-in postprocessing and protein inference procedures and provides graphic visualization of the statistical properties of the data set and the search results. It is open-source and can be freely extended to use third-party scoring functions or processing algorithms and allows customization of the search workflow for specialized applications.

MS/MS-Free Protein Identification in Complex Mixtures Using Multiple Enzymes with Complementary Specificity.

In this work, we present the results of evaluation of a workflow that employs a multienzyme digestion strategy for MS1-based protein identification in "shotgun" proteomic applications. In the proposed strategy, several cleavage reagents of different specificity were used for parallel digestion of the protein sample followed by MS1 and retention time (RT) based search. Proof of principle for the proposed strategy was performed using experimental data obtained for the annotated 48-protein standard. By using the developed approach, up to 90% of proteins from the standard were unambiguously identified. The approach was further applied to HeLa proteome data. For the sample of this complexity, the proposed MS1-only strategy determined correctly up to 34% of all proteins identified using standard MS/MS-based database search. It was also found that the results of MS1-only search were independent of the chromatographic gradient time in a wide range of gradients from 15-120 min. Potentially, rapid MS1-only proteome characterization can be an alternative or complementary to the MS/MS-based "shotgun" analyses in the studies, in which the experimental time is more important than the depth of the proteome coverage.

Comparative evaluation of label-free quantification methods for shotgun proteomics.

RATIONALE: Label-free quantification (LFQ) is a popular strategy for shotgun proteomics. A variety of LFQ algorithms have been developed recently. However, a comprehensive comparison of the most commonly used LFQ methods is still rare, in part due to a lack of clear metrics for their evaluation and an annotated and quantitatively well-characterized data set. METHODS: Five LFQ methods were compared: spectral counting based algorithms SIN , emPAI, and NSAF, and approaches relying on the extracted ion chromatogram (XIC) intensities, MaxLFQ and Quanti. We used three criteria for performance evaluation: coefficient of variation (CV) of protein abundances between replicates; analysis of variance (ANOVA); and the root-mean-square error of logarithmized calculated concentration ratios, referred to as standard quantification error (SQE). Comparison was performed using a quantitatively annotated publicly available data set. RESULTS: The best results in terms of inter-replicate reproducibility were observed for MaxLFQ and NSAF, although they exhibited larger standard quantification errors. Using NSAF, all quantitatively annotated proteins were correctly identified in the Bonferronni-corrected results of the ANOVA test. SIN was found to be the most accurate in terms of SQE. Finally, the current implementations of XIC-based LFQ methods did not outperform the methods based on spectral counting for the data set used in this study. CONCLUSIONS: Surprisingly, the performances of XIC-based approaches measured using three independent metrics were found to be comparable with more straightforward and simple MS/MS-based spectral counting approaches. The study revealed no clear leader among the latter. Copyright © 2017 John Wiley & Sons, Ltd.

Yeast Ribonucleotide Reductase Is a Direct Target of the Proteasome and Provides Hyper Resistance to the Carcinogen 4-NQO.

Various external and internal factors damaging DNA constantly disrupt the stability of the genome. Cells use numerous dedicated DNA repair systems to detect damage and restore genomic integrity in a timely manner. Ribonucleotide reductase (RNR) is a key enzyme providing dNTPs for DNA repair. Molecular mechanisms of indirect regulation of yeast RNR activity are well understood, whereas little is known about its direct regulation. The study was aimed at elucidation of the proteasome-dependent mechanism of direct regulation of RNR subunits in Saccharomyces cerevisiae. Proteome analysis followed by Western blot, RT-PCR, and yeast plating analysis showed that upregulation of RNR by proteasome deregulation is associated with yeast hyper resistance to 4-nitroquinoline-1-oxide (4-NQO), a UV-mimetic DNA-damaging drug used in animal models to study oncogenesis. Inhibition of RNR or deletion of RNR regulatory proteins reverses the phenotype of yeast hyper resistance to 4-NQO. We have shown for the first time that the yeast Rnr1 subunit is a substrate of the proteasome, which suggests a common mechanism of RNR regulation in yeast and mammals.

AA_stat: Intelligent profiling of in vivo and in vitro modifications from open search results.

Characterization of post-translational modifications is among the most challenging tasks in tandem mass spectrometry-based proteomics which has yet to find an efficient solution. The ultra-tolerant (open) database search attempts to meet this challenge. However, interpretation of the mass shifts observed in open search still requires an effective and automated solution. We have previously introduced the AA_stat tool for analysis of amino acid frequencies at different mass shifts and generation of hypotheses on unaccounted in vitro modifications. Here, we report on the new version of AA_stat, which now complements amino acid frequency statistics with a number of new features: (1) MS/MS-based localization of mass shifts and localization scoring, including shifts which are the sum of modifications; (2) inferring fixed modifications to increase method sensitivity; (3) inferring monoisotopic peak assignment errors and variable modifications based on abundant mass shift localizations to increase the yield of closed search; (4) new mass calibration algorithm to account for partial systematic shifts; (5) interactive integration of all results and a rated list of possible mass shift interpretations. With these options, we improve interpretation of open search results and demonstrate the utility of AA_stat for profiling of abundant and rare amino acid modifications. AA_stat is implemented in Python as an open-source tool available at https://github.com/SimpleNumber/aa_stat. SIGNIFICANCE: Mass spectrometry-based PTM characterization has a long history, yet most of the methods rely on a priori knowledge of modifications of interest and do not provide a whole proteome modification landscape in a blind manner. The open database search is an efficient attempt to address this challenge by identifying peptides with mass shifts corresponding to possible modifications. Then, interpreting these mass shifts is required. Therefore, development of bioinformatics software for post-processing of the open search results, which is capable of detection and accurate annotation of new or unexpected modifications, from characterization of sample preparation efficiency and quality control to discovery of rare post-translational modifications, is of high importance.

Improving the Protein Inference from Bottom-Up Proteomic Data Using Identifications from MS1 Spectra.

Protein inference is one of the crucial steps in proteome characterization using a bottom-up approach. Multiple algorithms to solve the problem are focused on extensive analysis of shared peptides identified from fragmentation mass spectra (MS/MS). However, many protein homologues with a similar amino acid sequence typically have identical lists of identified peptides due to the problem of proteome undersampling in a bottom-up approach and, thus, cannot be distinguished by existing protein inference methods. Here, we propose the use of peptide feature information extracted from precursor mass spectra to assist in identification of proteins otherwise indistinguishable from MS/MS. The proposed method was integrated with a protein inference algorithm based on the parsimony principle and built-in in the postsearch utility Scavager. The results demonstrate increasing accuracy and efficiency of homologous protein identifications for the well characterized data sets including the one with known protein sequences from iPRG-2016 study.