Exploring the dynamic landscape of immunopeptidomics: Unravelling posttranslational modifications and navigating bioinformatics terrain.

Immunopeptidomics is becoming an increasingly important field of study. The capability to identify immunopeptides with pivotal roles in the human immune system is essential to shift the current curative medicine towards personalized medicine. Throughout the years, the field has matured, giving insight into the current pitfalls. Nowadays, it is commonly accepted that generalizing shotgun proteomics workflows is malpractice because immunopeptidomics faces numerous challenges. While many of these difficulties have been addressed, the road towards the ideal workflow remains complicated. Although the presence of Posttranslational modifications (PTMs) in the immunopeptidome has been demonstrated, their identification remains highly challenging despite their significance for immunotherapies. The large number of unpredictable modifications in the immunopeptidome plays a pivotal role in the functionality and these challenges. This review provides a comprehensive overview of the current advancements in immunopeptidomics. We delve into the challenges associated with identifying PTMs within the immunopeptidome, aiming to address the current state of the field.

A compositional data model to predict the isotope distribution for average peptides using a compositional spline model.

We propose an updated approach for approximating the isotope distribution of average peptides given their monoisotopic mass. Our methodology involves in-silico cleavage of the entire UNIPROT database of human-reviewed proteins using Trypsin, generating a theoretical peptide dataset. The isotope distribution is computed using BRAIN. We apply a compositional data modelling strategy that utilizes an additive log-ratio transformation for the isotope probabilities followed by a penalized spline regression. Furthermore, due to the impact of the number of sulphur atoms on the course of the isotope distribution, we develop separate models for peptides containing zero up to five sulphur atoms. Additionally, we propose three methods to estimate the number of sulphur atoms based on an observed isotope distribution. The performance of the spline models and the sulphur prediction approaches is evaluated using a mean squared error and a modified Pearson's χ2 goodness-of-fit measure on an experimental UPS2 data set. Our analysis reveals that the variability in spectral accuracy, that is, the variability between MS1 scans, contributes more to the errors than the approximation of the theoretical isotope distribution by our proposed average peptide model. Moreover, we find that the accuracy of predicting the number of sulphur atoms based on the observed isotope distribution is limited by measurement accuracy.

CPred: Charge State Prediction for Modified and Unmodified Peptides in Electrospray Ionization.

The mass-to-charge ratio serves as a critical parameter in peptide identification via mass spectrometry, enabling the precise determination of peptide masses and facilitating their differentiation based on unique charge characteristics, especially when peptides are ionized by tools like electrospray ionization, which produces multiply charged ions. We developed a neural network called CPred, which can accurately predict the charge state distribution from +1 to +7 for the modified and unmodified peptides. CPred was trained on the large-scale synthetic training data, consisting of tryptic and non-tryptic peptides, and various fragmentation methods. The model was further evaluated on independent, external test data sets. Results were evaluated through the Pearson correlation coefficient and showed high correlations of up to 0.9997117 between the predicted and acquired charge state distributions. The effect of specifying modifications in the neural network and feature importance was further investigated, revealing the value of modifications and vital peptide properties in holding on to protons. CPreds' accurate predictions of the charge state distribution can play an essential role in boosting confidence in peptide identifications during rescoring as a novel feature.

MIND4OLIGOS: Determining the Monoisotopic Mass of Oligonucleotides Observed in High-Resolution Mass Spectrometry.

Oligonucleotide therapeutics have emerged as an important class of drugs offering targeted therapeutic strategies that complement traditional modalities, such as monoclonal antibodies and small molecules. Their unique ability to precisely modulate gene expression makes them vital for addressing previously undruggable targets. A critical aspect of developing these therapies is characterizing their molecular composition accurately. This includes determining the monoisotopic mass of oligonucleotides, which is essential for identifying impurities, degradants, and modifications that can affect the drug efficacy and safety. Mass spectrometry (MS) plays a pivotal role in this process, yet the accurate interpretation of complex mass spectra remains challenging, especially for large molecules, where the monoisotopic peak is often undetectable. To address this issue, we have adapted the MIND algorithm, originally developed for top-down proteomics, for use with oligonucleotide data. This adaptation allows for the prediction of monoisotopic mass from the more readily detectable, most-abundant peak mass, enhancing the ability to annotate complex spectra of oligonucleotides. Our comprehensive validation of this modified algorithm on both in silico and real-world oligonucleotide data sets has demonstrated its effectiveness and reliability. To facilitate wider adoption of this advanced analytical technique, we have encapsulated the enhanced MIND algorithm in a user-friendly Shiny application. This online platform simplifies the process of annotating complex oligonucleotide spectra, making advanced mass spectrometry analysis accessible to researchers and drug developers. The application is available at https://valkenborg-lab.shinyapps.io/mind4oligos/.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania.

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

The isotope distribution: A rose with thorns.

The isotope distribution, which reflects the number and probabilities of occurrence of different isotopologues of a molecule, can be theoretically calculated. With the current generation of (ultra)-high-resolution mass spectrometers, the isotope distribution of molecules can be measured with high sensitivity, resolution, and mass accuracy. However, the observed isotope distribution can differ substantially from the expected isotope distribution. Although differences between the observed and expected isotope distribution can complicate the analysis and interpretation of mass spectral data, they can be helpful in a number of specific applications. These applications include, yet are not limited to, the identification of peptides in proteomics, elucidation of the elemental composition of small organic molecules and metabolites, as well as wading through peaks in mass spectra of complex bioorganic mixtures such as petroleum and humus. In this review, we give a nonexhaustive overview of factors that have an impact on the observed isotope distribution, such as elemental isotope deviations, ion sampling, ion interactions, electronic noise and dephasing, centroiding, and apodization. These factors occur at different stages of obtaining the isotope distribution: during the collection of the sample, during the ionization and intake of a molecule in a mass spectrometer, during the mass separation and detection of ionized molecules, and during signal processing.

Machine learning approach for the prediction of the number of sulphur atoms in peptides using the theoretical aggregated isotope distribution.

RATIONALE: The observed isotope distribution is an important attribute for the identification of peptides and proteins in mass spectrometry-based proteomics. Sulphur atoms have a very distinctive elemental isotope definition, and therefore, the presence of sulphur atoms has a substantial effect on the isotope distribution of biomolecules. Hence, knowledge of the number of sulphur atoms can improve the identification of peptides and proteins. METHODS: In this paper, we conducted a theoretical investigation on the isotope properties of sulphur-containing peptides. We proposed a gradient boosting approach to predict the number of sulphur atoms based on the aggregated isotope distribution. We compared prediction accuracy and assessed the predictive power of the features using the mass and isotope abundance information from the first three, five and eight aggregated isotope peaks. RESULTS: Mass features alone are not sufficient to accurately predict the number of sulphur atoms. However, we reach near-perfect prediction when we include isotope abundance features. The abundance ratios of the eighth and the seventh, the fifth and the fourth, and the third and the second aggregated isotope peaks are the most important abundance features. The mass difference between the eighth, the fifth or the third aggregated isotope peaks and the monoisotopic peak are the most predictive mass features. CONCLUSIONS: Based on the validation analysis it can be concluded that the prediction of the number of sulphur atoms based on the isotope profile fails, because the isotope ratios are not measured accurately. These results indicate that it is valuable for future instrument developments to focus more on improving spectral accuracy to measure peak intensities of higher-order isotope peaks more accurately.

CONSTANd: An Efficient Normalization Method for Relative Quantification in Small- and Large-Scale Omics Experiments in R BioConductor and Python.

For differential expression studies in all omics disciplines, data normalization is a crucial step that is often subject to a balance between speed and effectiveness. To keep up with the data produced by high-throughput instruments, researchers require fast and easy-to-use yet effective methods that fit into automated analysis pipelines. The CONSTANd normalization method meets these criteria, so we have made its source code available for R/BioConductor and Python. We briefly review the method and demonstrate how it can be used in different omics contexts for experiments of any scale. Widespread adoption across omics disciplines would ease data integration in multiomics experiments.

Predicting the number of sulfur atoms in peptides and small proteins based on the observed aggregated isotope distribution.

RATIONALE: Identification of peptides and proteins is a challenging task in mass spectrometry-based proteomics. Knowledge of the number of sulfur atoms can improve the identification of peptides and proteins. METHODS: In this article, we propose a method for the prediction of S-atoms based on the aggregated isotope distribution. The Mahalanobis distance is used as dissimilarity measure to compare mass- and intensity-based features from the observed and theoretical isotope distributions. RESULTS: The relative abundance of the second and the third aggregated isotopic variants (as compared to the monoisotopic one) and the mass difference between the second and third aggregated isotopic variants are the most important features to predict the number of S-atoms. CONCLUSIONS: The mass and intensity accuracies of the observed aggregated isotopic variants are insufficient to accurately predict the number of atoms. However, using a limited set of predictions for a peptide, rather than predicting a single number of S-atoms, has a reasonably high prediction accuracy.

Classification with a Deferral Option and Low-Trust Filtering for Automated Seizure Detection.

Wearable technology will become available and allow prolonged electroencephalography (EEG) monitoring in the home environment of patients with epilepsy. Neurologists analyse the EEG visually and annotate all seizures, which patients often under-report. Visual analysis of a 24-h EEG recording typically takes one to two hours. Reliable automated seizure detection algorithms will be crucial to reduce this analysis. We investigated such algorithms on a dataset of behind-the-ear EEG measurements. Our first aim was to develop a methodology where part of the data is deferred to a human expert, who performs perfectly, with the goal of obtaining an (almost) perfect detection sensitivity (DS). Prediction confidences are determined by temperature scaling of the classification model outputs and trust scores. A DS of approximately 90% (99%) can be achieved when deferring around 10% (40%) of the data. Perfect DS can be achieved when deferring 50% of the data. Our second contribution demonstrates that a common modelling strategy, where predictions from several short EEG segments are combined to obtain a final prediction, can be improved by filtering out untrustworthy segments with low trust scores. The false detection rate shows a relative decrease between 21% and 43%, and the DS shows a small increase or decrease.

IsoSpec2: Ultrafast Fine Structure Calculator.

High-resolution mass spectrometry becomes increasingly available with its ability to resolve the fine isotopic structure of measured analytes. It allows for high-sensitivity spectral deconvolution, leading to less false-positive identifications. Analytes can be identified by comparing their theoretical isotopic signal with the observed peaks. Necessary calculations are, however, computationally demanding and lead to long processing times. For wheat (trictum oestivum) alone, Uniprot holds more than 142 000 candidate protein sequences. This is doubled upon sequence reversal for identification FDR estimation and further multiplied by performing in silico digestion into peptides. The same peptide might originate from more than one protein, which reduces the overall number of sequences to be calculated. However, it is still huge. IsoSpec2 can perform these calculations fast. Compared to IsoSpec1, the algorithm is simpler, orders of magnitude faster, and offers more flexibility for the developers of algorithms for raw data analysis. It is freely available under a 2-clause BSD license, with bindings for the C++, C, R, and Python programming languages.

The challenges of peptidomics in complementing proteomics in a clinical context.

Naturally occurring peptides, including growth factors, hormones, and neurotransmitters, represent an important class of biomolecules and have crucial roles in human physiology. The study of these peptides in clinical samples is therefore as relevant as ever. Compared to more routine proteomics applications in clinical research, peptidomics research questions are more challenging and have special requirements with regard to sample handling, experimental design, and bioinformatics. In this review, we describe the issues that confront peptidomics in a clinical context. After these hurdles are (partially) overcome, peptidomics will be ready for a successful translation into medical practice.

Starch biosynthesis contributes to the maintenance of photosynthesis and leaf growth under drought stress in maize.

To understand the growth response to drought, we performed a proteomics study in the leaf growth zone of maize (Zea mays L.) seedlings and functionally characterized the role of starch biosynthesis in the regulation of growth, photosynthesis and antioxidant capacity, using the shrunken-2 mutant (sh2), defective in ADP-glucose pyrophosphorylase. Drought altered the abundance of 284 proteins overrepresented for photosynthesis, amino acid, sugar and starch metabolism, and redox-regulation. Changes in protein levels correlated with enzyme activities (increased ATP synthase, cysteine synthase, starch synthase, RuBisCo, peroxiredoxin, glutaredoxin, thioredoxin and decreased triosephosphate isomerase, ferredoxin, cellulose synthase activities, respectively) and metabolite concentrations (increased ATP, cysteine, glycine, serine, starch, proline and decreased cellulose levels). The sh2 mutant showed a reduced increase of starch levels under drought conditions, leading to soluble sugar starvation at the end of the night and correlating with an inhibition of leaf growth rates. Increased RuBisCo activity and pigment concentrations observed in WT, in response to drought, were lacking in the mutant, which suffered more oxidative damage and recovered more slowly after re-watering. These results demonstrate that starch biosynthesis contributes to maintaining leaf growth under drought stress and facilitates enhanced carbon acquisition upon recovery.

PRiSM: A prototype for exhaustive, restriction-free database searching for mass spectrometry-based identification.

RATIONALE: The current methods for identifying peptides in mass spectral product ion data still struggle to do so for the majority of spectra. Based on the experimental setup and other assumptions, such methods restrict the search space to speed up computations, but at the cost of creating blind spots. The proteomics community would greatly benefit from a method that is capable of covering the entire search space without using any restrictions, thus establishing a baseline for identification. METHODS: We conceived the "mass pattern paradigm" (MPP) that enables the creation of such an identification method, and we implemented it into a prototype database search engine "PRiSM" (PRotein-Spectrum Matching). We then assessed its operational characteristics by applying it to publicly available high-precision mass spectra of low and high identification difficulty. We used those characteristics to gain theoretical insights into trade-offs between sensitivity and speed when trying to establish a baseline for identification. RESULTS: Of 100 low difficulty spectra, PRiSM and SEQUEST agree on 84 identifications (of which 75 are statistically significant). Of 15 of 100 spectra not identified in a previous study (using SEQUEST), 13 are considered reliable after visual inspection and represent 3 proteins (out of 9 in total) not detected previously. CONCLUSIONS: Despite leaving noise intact, the simple PRiSM prototype can make statistically reliable identifications, while controlling the false discovery rate by fitting a null distribution. It also identifies some spectra previously unidentifiable in an "extremely open" SEQUEST search, paving the way to establishing a baseline for identification in proteomics.

Masserstein: Linear regression of mass spectra by optimal transport.

RATIONALE: The linear regression of mass spectra is a computational problem defined as fitting a linear combination of reference spectra to an experimental one. It is typically used to estimate the relative quantities of selected ions. In this work, we study this problem in an abstract setting to develop new approaches applicable to a diverse range of experiments. METHODS: To overcome the sensitivity of the ordinary least-squares regression to measurement inaccuracies, we base our methods on a non-conventional spectral dissimilarity measure, known as the Wasserstein or the Earth Mover's distance. This distance is based on the notion of the cost of transporting signal between mass spectra, which renders it naturally robust to measurement inaccuracies in the mass domain. RESULTS: Using a data set of 200 mass spectra, we show that our approach is capable of estimating ion proportions accurately without extensive preprocessing of spectra required by other methods. The conclusions are further substantiated using data sets simulated in a way that mimics most of the measurement inaccuracies occurring in real experiments. CONCLUSIONS: We have developed a linear regression algorithm based on the notion of the cost of transporting signal between spectra. Our implementation is available in a Python 3 package called masserstein, which is freely available at https://github.com/mciach/masserstein.

Machine learning analysis of motor evoked potential time series to predict disability progression in multiple sclerosis.

BACKGROUND: Evoked potentials (EPs) are a measure of the conductivity of the central nervous system. They are used to monitor disease progression of multiple sclerosis patients. Previous studies only extracted a few variables from the EPs, which are often further condensed into a single variable: the EP score. We perform a machine learning analysis of motor EP that uses the whole time series, instead of a few variables, to predict disability progression after two years. Obtaining realistic performance estimates of this task has been difficult because of small data set sizes. We recently extracted a dataset of EPs from the Rehabiliation & MS Center in Overpelt, Belgium. Our data set is large enough to obtain, for the first time, a performance estimate on an independent test set containing different patients. METHODS: We extracted a large number of time series features from the motor EPs with the highly comparative time series analysis software package. Mutual information with the target and the Boruta method are used to find features which contain information not included in the features studied in the literature. We use random forests (RF) and logistic regression (LR) classifiers to predict disability progression after two years. Statistical significance of the performance increase when adding extra features is checked. RESULTS: Including extra time series features in motor EPs leads to a statistically significant improvement compared to using only the known features, although the effect is limited in magnitude (ΔAUC = 0.02 for RF and ΔAUC = 0.05 for LR). RF with extra time series features obtains the best performance (AUC = 0.75±0.07 (mean and standard deviation)), which is good considering the limited number of biomarkers in the model. RF (a nonlinear classifier) outperforms LR (a linear classifier). CONCLUSIONS: Using machine learning methods on EPs shows promising predictive performance. Using additional EP time series features beyond those already in use leads to a modest increase in performance. Larger datasets, preferably multi-center, are needed for further research. Given a large enough dataset, these models may be used to support clinicians in their decision making process regarding future treatment.

QCQuan: A Web Tool for the Automated Assessment of Protein Expression and Data Quality of Labeled Mass Spectrometry Experiments.

In the context of omics disciplines and especially proteomics and biomarker discovery, the analysis of a clinical sample using label-based tandem mass spectrometry (MS) can be affected by sample preparation effects or by the measurement process itself, resulting in an incorrect outcome. Detection and correction of these mistakes using state-of-the-art methods based on mixed models can use large amounts of (computing) time. MS-based proteomics laboratories are high-throughput and need to avoid a bottleneck in their quantitative pipeline by quickly discriminating between high- and low-quality data. To this end we developed an easy-to-use web-tool called QCQuan (available at qcquan.net ) which is built around the CONSTANd normalization algorithm. It automatically provides the user with exploratory and quality control information as well as a differential expression analysis based on conservative, simple statistics. In this document we describe in detail the scientifically relevant steps that constitute the workflow and assess its qualitative and quantitative performance on three reference data sets. We find that QCQuan provides clear and accurate indications about the scientific value of both a high- and a low-quality data set. Moreover, it performed quantitatively better on a third data set than a comparable workflow assembled using established, reliable software.

masstodon: A Tool for Assigning Peaks and Modeling Electron Transfer Reactions in Top-Down Mass Spectrometry.

Top-down mass spectrometry methods are becoming continuously more popular in the effort to describe the proteome. They rely on the fragmentation of intact protein ions inside the mass spectrometer. Among the existing fragmentation methods, electron transfer dissociation is known for its precision and wide coverage of different cleavage sites. However, several side reactions can occur under electron transfer dissociation (ETD) conditions, including nondissociative electron transfer and proton transfer reaction. Evaluating their extent can provide more insight into reaction kinetics as well as instrument operation. Furthermore, preferential formation of certain reaction products can reveal important structural information. To the best of our knowledge, there are currently no tools capable of tracing and analyzing the products of these reactions in a systematic way. In this Article, we present in detail masstodon: a computer program for assigning peaks and interpreting mass spectra. Besides being a general purpose tool, masstodon also offers the possibility to trace the products of reactions occurring under ETD conditions and provides insights into the parameters driving them. It is available free of charge under the GNU AGPL V3 public license.

MIND: A Double-Linear Model To Accurately Determine Monoisotopic Precursor Mass in High-Resolution Top-Down Proteomics.

Top-down proteomics approaches are becoming ever more popular, due to the advantages offered by knowledge of the intact protein mass in correctly identifying the various proteoforms that potentially arise due to point mutation, alternative splicing, post-translational modifications, etc. Usually, the average mass is used in this context; however, it is known that this can fluctuate significantly due to both natural and technical causes. Ideally, one would prefer to use the monoisotopic precursor mass, but this falls below the detection limit for all but the smallest proteins. Methods that predict the monoisotopic mass based on the average mass are potentially affected by imprecisions associated with the average mass. To address this issue, we have developed a framework based on simple, linear models that allows prediction of the monoisotopic mass based on the exact mass of the most-abundant (aggregated) isotope peak, which is a robust measure of mass, insensitive to the aforementioned natural and technical causes. This linear model was tested experimentally, as well as in silico, and typically predicts monoisotopic masses with an accuracy of only a few parts per million. A confidence measure is associated with the predicted monoisotopic mass to handle the off-by-one-Da prediction error. Furthermore, we introduce a correction function to extract the "true" (i.e., theoretically) most-abundant isotope peak from a spectrum, even if the observed isotope distribution is distorted by noise or poor ion statistics. The method is available online as an R shiny app: https://valkenborg-lab.shinyapps.io/mind/.

Radical solutions: Principles and application of electron-based dissociation in mass spectrometry-based analysis of protein structure.

In recent years, electron capture (ECD) and electron transfer dissociation (ETD) have emerged as two of the most useful methods in mass spectrometry-based protein analysis, evidenced by a considerable and growing body of literature. In large part, the interest in these methods is due to their ability to induce backbone fragmentation with very little disruption of noncovalent interactions which allows inference of information regarding higher order structure from the observed fragmentation behavior. Here, we review the evolution of electron-based dissociation methods, and pay particular attention to their application in "native" mass spectrometry, their mechanism, determinants of fragmentation behavior, and recent developments in available instrumentation. Although we focus on the two most widely used methods-ECD and ETD-we also discuss the use of other ion/electron, ion/ion, and ion/neutral fragmentation methods, useful for interrogation of a range of classes of biomolecules in positive- and negative-ion mode, and speculate about how this exciting field might evolve in the coming years.