Peptide collision cross sections of 22 post-translational modifications

Recent advances have rekindled the interest in ion mobility as an additional dimension of separation in mass spectrometry (MS)-based proteomics. Ion mobility separates ions according to their size and shape in the gas phase. Here, we set out to investigate the effect of 22 different post-translational modifications (PTMs) on the collision cross section (CCS) of peptides. In total, we analyzed ~4300 pairs of matching modified and unmodified peptide ion species by trapped ion mobility spectrometry (TIMS). Linear alignment based on spike-in reference peptides resulted in highly reproducible CCS values with a median coefficient of variation of 0.26%. On a global level, we observed a redistribution in the m/z vs. ion mobility space for modified peptides upon changes in their charge state. Pairwise comparison between modified and unmodified peptides of the same charge state revealed median shifts in CCS between −1.4% (arginine citrullination) and +4.5% (O-GlcNAcylation). In general, increasing modified peptide masses were correlated with higher CCS values, in particular within homologous PTM series. However, investigating the ion populations in more detail, we found that the change in CCS can vary substantially for a given PTM and is partially correlated with the gas phase structure of its unmodified counterpart. In conclusion, our study shows PTM- and sequence-specific effects on the cross section of peptides, which could be further leveraged for proteome-wide PTM analysis. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1007/s00216-023-04957-4.


Introduction
Post-translational modifications (PTMs) are key regulators of protein activity and function in health and disease.Not at least because most PTMs involve a specific shift in the molecular weight of modified amino acids, mass spectrometry (MS)-based proteomics has evolved as the method of choice for the investigation of PTMs on a proteome-wide scale [1][2][3].MS-based proteomics analyzes complex mixtures of (modified) peptides derived from tryptic digests by liquid chromatography coupled to high-resolution MS.Since the modification mass shift also transfers to fragment ions, this allows identifying modified peptide sequences and also often localizing the PTM to a specific amino acid.
Recently, ion mobility spectrometry (IMS) has become a popular extension of the proteomics toolbox that adds one more dimension of separation [4][5][6][7][8][9][10].IMS distinguishes ions in the gas phase by their size and shape, which can be inferred from time-and field-dispersive ion mobility measurements in the form of a collision cross section (CCS) or, more precisely, the momentum transfer collision integral [11][12][13].In the CCS vs m/z dimension, electrosprayed tryptic peptides typically split into distinct populations according to their charge state [14].A more fine-structured heterogeneity, most prominently for triply charged species, has been attributed to either more extended or more compact structures that are determined by the linear peptide sequence [15,16].As a result, even isobaric and isomeric peptides can have distinct cross sections [17,18].
Trapped ion mobility spectrometry (TIMS) is a relatively new type of IMS that inverses the concept of classical drift tube IMS by holding ions in an electric field against an opposing gas flow [19][20][21].Lowering the electric field strength releases ions sequentially from the TIMS device to the downstream mass analyzer as a function of increasing ion mobility (or decreasing CCS).The PASEF acquisition mode synchronizes the separation with a quadrupole timeof-flight mass analyzer, which greatly enhances the speed and sensitivity of peptide sequencing [5,22,23].We have recently shown that an intriguing feature of this setup is that it enables peptide CCS measurements at very large scale and with high precision, sufficient to train a deep learning model to predict peptide cross sections based solely given the linear amino acid sequence and charge state as an input [16].
The advantages of IMS in terms of speed, sensitivity, and specificity should equally apply to or even be enhanced in the analysis of PTMs, given the additional complexity arising from isobaric modifications and positional isomers [24][25][26][27].This motivated researchers since the beginning of IMS to study the effect of specific modifications on the gas phase characteristics of model peptides.By far the most studied example is phosphorylation of serine and threonine residues, which, interestingly, often leads to more compact gas phase structures compared to unmodified peptides -despite the increase in mass [26,[28][29][30].Another field of interest is N-glycosylated peptides, for which the glycan moiety results in a distinct separation from unmodified peptides in the ion mobility dimension, and even differentiation between isomeric localization variants has been demonstrated [27,31,32].To model the effect of PTMs on peptide CCS values, Clemmer and co-workers generated a dataset of cysteinepalmitoylated as well as cysteine-carboxyamidomethylated peptides and derived modification-specific intrinsic size parameters [33].Kaszycki and Shvartsburg extended this approach to predict the intrinsic size parameters of 100 different PTMs [34].More recent large-scale studies further highlighted additional sequence-dependent effects on peptide cross sections resulting from intramolecular interactions [16,26,35].
The further exploration of the potential of IMS for PTM analysis and the development of accurate CCS prediction models for a wide range of PTMs, similar to those for unmodified peptides [16,35,36], is currently limited by the availability of comprehensive experimental data.Here, we investigate the effect of 22 different naturally occurring PTMs on the collision cross section of peptides, including 21 sets of modified peptides synthesized as part of the Pro-teomeTools project [37,38] and an additional set of O-Glc-NAcylated peptides.

Synthetic peptides
A library of ~4500 lyophilized synthetic peptides in 96-well format representing 21 naturally occurring post-translational modifications on K, R, P, and Y residues and their unmodified counterparts were obtained from the ProteomeTools project [38].Additionally, we purchased synthetic O-Glc-NAcylated peptides and corresponding unmodified peptides from JPT Peptide Technologies GmbH (SpikeMix PTM-Kit 57 and 55).The pre-pooled peptides were reconstituted in 2% acetonitrile/0.1% formic acid to a final concentration of ~100 fmol/µL.As a reference, we spiked each sample with a retention time standard (Biognosys iRT) in a ratio of 1:40 (vol/vol) [39].

Liquid chromatography and mass spectrometry
All solvents were HPLC grade and purchased from Sigma-Aldrich.Nanoflow reversed-phase liquid chromatography was performed on a nanoElute system (Bruker Daltonics).Peptides were separated with a 120-min gradient at a flow rate of 0.3 µL/min at 60 °C on a homemade 50 cm × 75 µm column with a pulled emitter tip, packed with 1.9-μm ReproSil-Pur C18 -AQ beads.Mobile phases A and B were water/0.1% formic acid and acetonitrile/0.1% formic acid.The LC system was connected online to a TIMS-quadrupole time-of-flight mass spectrometer (Bruker timsTOF Pro or timsTOF HT) via a CaptiveSpray nano-electrospray source [5].TIMS analysis was performed in a range from 1/K 0 = 1.5 to 0.6 Vs cm −2 , while accumulating and analyzing in parallel for 100 ms each.The ion mobility gas was nitrogen from ambient air without temperature control, and the pressure at the TIMS entrance was kept at ~2.4 mbar.Detailed experimental parameters are provided in Supplementary Tables 1 and 2. To calibrate the ion mobility dimension, we added low-concentration Agilent ESI LC/MS to the inlet filter of the CaptiveSpray ion source and fitted the TIMS elution voltages to 1/K 0 values using a linear model for at least three ions (m/z, 1/K 0 : 622.0289, 0.9848 Vs cm −2 ; 922.0097, 1.1895 Vs cm −2 ; and 1221.9906,1.3820 Vs cm −2 ).All data were acquired in dda-PASEF mode and suitable precursor ions for fragmentation were selected by their relative position in the 1/K 0 vs. m/z plane.The quadrupole isolation window was set to 2 Th for m/z < 700 and 3 Th for m/z > 700.To prevent repeated fragmentation of precursors for which the dda-PASEF target intensity value of 20,000 a.u.has been already achieved, we defined an exclusion time of 0.4 min.

Data processing
The MS raw data were processed with MaxQuant (version 2.1.4.0) using the default parameters for the analysis of timsTOF data [40,41].MS/MS spectra were searched against concatenated peptide sequences as provided by the ProteomeTools consortium and JPT Peptide Technologies, both supplemented with the iRT peptide sequences.The digestion mode was set to "specific" according to the cleavage rules for Trypsin.We analyzed each PTM pool separately, defining methionine oxidation and each peptide pool`s respective modification as variable modifications (Supplementary Table 3) and cysteine carbamidomethylation as a fixed modification.The false discovery rate on the peptide spectrum match and protein level was controlled by a target-decoy approach <1%.For modified peptides, we required a minimum Andromeda score [42] of 40 and a minimum delta score of 6 [43].
To account for time-dependent drifts in measured ion mobility (1/K 0 ) and retention time (RT) values, we linearly aligned all experiments by a run-wise alignment value (see also Results).To determine appropriate reference values for this, we first measured three replicate injections of the 11-iRT peptide mixture, calibrating the TIMS dimension before each injection as described above, and calculated the mean 1/K 0 (or RT) value for each iRT peptide (Supplementary Table 4).As we spiked iRT peptides in all our measurements, we could then extract the experimental 1/K 0 for each identified iRT peptide in each experiment.From this, and separately for each run, we calculated the median deviation by subtracting each measured iRT peptide 1/K 0 (and RT) value from its respective reference value.This resulted in one value per run, which we used to align all runs by subtracting it from all measured 1/K 0 values, run by run.The aligned 1/K 0 values were converted to CCS values using the Mason-Schamp equation.For further analysis, we kept only peptide sequences identified with a single modification site not localized on the C-terminus, excluding oxidized peptides.In cases, in which MaxQuant listed multiple "evidences" for a modified peptide sequence, we retained only the most abundant feature (highest intensity value) for charge state 2 and 3 for our analysis.CCS and RT values represent mean values calculated from one to three technical replicates.Differences in the collision cross section of modified peptides and their unmodified counterparts were calculated relative to the cross section of the unmodified peptide (ΔCCS = (CCS modified − CCS unmodified )/CCS unmodified ).

Modeling of peptide collision cross sections
The potential energy surface of the selected doubly protonated peptides was explored by molecular dynamics (MD) simulations using the OPLS all-atom force field [44,45] in conjunction with the GROMACS suite of programs.For simulation of peptides with succinyl-modified lysine residues, we used the parameters tabulated for aliphatic molecular systems [46].Extra protons were placed on the most basic sites of the peptides, i.e., the lysine residues that are known to sequester protons.We further used a deprotonated carboxy-terminus and a protonated amino-terminus.For the succinylated sequences, which lack one of the basic lysine residues, we neutralized the carboxy-terminus.Note that the semi-empirical calculations performed subsequent to the MD simulations would allow proton transfer processes to take place if they were to lead to energetically more stable structures.During the MD simulations, we used simulated annealing techniques to produce candidate structures for further refinement.These simulations were carried out for a duration of 100 ns at simulation temperatures of up to 600 K to effectively overcome barriers on the potential energy surface.Snapshots were saved every 100 ps during these simulations and analyzed by a conformer family search program which assigns structures into families within which the most important characteristic torsion angles are similar.Note that this approach was extensively used to discuss gasphase fragmentation pathways of protonated peptides [46][47][48].A full geometry optimization was subsequently carried out for the most stable structure of each conformer family at the PM6 level of theory [49] using the MOPAC suite of programs [50].The optimized structure was then used as input for a cross section calculation using our projection superposition approximation (PSA) method [51] in nitrogen gas [52].We restricted the putative assignment of gas phase structures to experimental mobility spectra to the respective lowest-energy conformers.

Data availability
The mass spectrometry proteomics data underlying this study have been deposited to the ProteomeXchange Consortium [53,54] via the PRIDE partner repository with the dataset identifier PXD042416.A summary result file is provided as Supplementary File 2. Code to reproduce the data analysis and visualization in this study can be accessed via https:// github.com/ Meier Lab/ 2022_ 22PTM.

Constructing a peptide CCS dataset for 22 PTMs
To generate a high-quality CCS dataset of a wide range of naturally occurring peptide modifications, we analyzed pooled libraries of synthetic peptides by nanoflow reversedphase liquid chromatography and TIMS-quadrupole timeof-flight mass spectrometry (Fig. 1a).In dda-PASEF mode, the mass spectrometer selects suitable precursor ions from a survey TIMS-MS scan and targets them for fragmentation in the subsequent PASEF-MS/MS scans.As quadrupole and collision cell are positioned downstream of the TIMS analyzer, precursor and fragment spectra are linked through their position in the ion mobility spectrum.We processed this data in the MaxQuant software to assemble three-dimensional features in ion mobility, m/z and retention time dimension and match the associated MS/MS spectra to (modified) peptide sequences [41].The inverse reduced ion mobility value (1/K 0 ) determined from the mobility spectrum of each feature can then be converted into ion-nitrogen TIMS CCS N2 values using the Mason-Schamp equation [55].
Figure 1b shows an overview of all PTMs in our dataset.The ProteomeTools library contributed matching modified and unmodified tryptic peptides derived from human protein sequences, which were selected for their favorable LC-MS properties and synthesizability as described in more detail in the original publications [37,38].The library contains four different sets of base peptide sequences that carry N-terminal or internal modifications on one of four amino acids (lysine, arginine, tyrosine, and proline).All peptide sequences with the same modification are combined, resulting in 21 pools of modified peptides and four matching pools of unmodified peptides, which we measured in randomized order and in triplicate.In addition, we analyzed a pool of O-Glc-NAcylated peptides (serine/threonine) and their unmodified counterparts in the same way.Overall, lysine modifications represent the largest group in our dataset, including three homologous series for acylation with aliphatic residues (formylation to butyrylation), carboxylic acid residues (malonylation to glutarylation) and methylation.Further, lysine biotinylation is included as well as the GlyGly remnant of ubiquitination and hydroxylated proline.Tyrosine modifications in our dataset are phosphorylation and nitration.The In total, we compiled 84 raw files, which resulted in ~58,000 peptide spectrum matches mapping to ~5000 unique combinations of peptide sequence, charge state, and modification.Of these, 74% and 26% were detected as doubly and triply charged species, and only 1% in charge state 4. The median Andromeda score was 105 with a very high localization probability close to 1 for all modifications except for O-GlcNAc, which is labile in collision induced dissociation experiments (Supplementary Fig. 1).Plotting the m/z vs. CCS distribution of modified and unmodified peptides shows the expected clustering by charge state distributed over an m/z range of about 300-1200 and a CCS range of about 300-700 Å 2 (Supplementary Fig. 2).For further analysis, we extracted the CCS value of the most abundant feature for each modified peptide sequence, while keeping doubly and triply charged ions separate because of their distinct ion mobility.This yielded approximately 4300 matched pairs of modified and unmodified peptides and about 115 to 292 pairs per modification.

Precision of TIMS CCS measurements
Experimental ion mobility values depend on the analyte itself as well as the electric field and experimental parameters such as temperature and the nature of the mobility gas [13,56].To make them comparable between experiments, TIMS is usually calibrated by a linear regression of known 1/K 0 values to the elution voltage, resulting in good agreement with conventional drift tube experiments [21].
In our previous study, we demonstrated that TIMS CCS N2 values from multiple experiments can additionally be linearly aligned based on overlapping peptide identifications, leading to a remarkable reproducibility over long periods of time and across instruments [16].To allow a similar alignment in our dataset of non-overlapping peptide pools and avoid external calibration before each injection, we here spiked a standard of eleven synthetic iRT peptides into each sample.The peptides eluted evenly distributed over the 2-h chromatographic gradient and were detected as doubly protonated species with CCS values ranging from 320 to 430 Å 2 .Our analysis revealed a time-dependent shift of their measured CCS values over the course of the experiment (84 LC-MS injections, acquired on two independent instruments), which could be attributed to changes in experimental conditions such as ambient temperature and pressure that are typically not controlled in TIMS experiments [11] (Fig. 2a).A linear alignment to the average iRT 1/K 0 values from three reference experiments successfully corrected these drifts (Methods, Fig. 2a lower panel), while using the median deviation as a correction factor makes it robust against outliers (e.g., "DGLDAA…" in Fig. 2a).Across all experiments, the median coefficient of variation of the iRT peptide CCS values was 1.41% before and 0.30% after the alignment.
Next, we determined the precision of the CCS measurements for all other peptides in our dataset (Fig. 2b).Because of the high sequencing rate of dda-PASEF and the relatively low complexity of the synthetic peptide pools, 76% of all peptides were identified in all three replicates.For these, the coefficient of variation was significantly improved (p < 0.001, Kolmogorov-Smirnov test) from 1.23% before to only 0.26% after linear alignment.This indicates an excellent reproducibility of our TIMS measurements and is in line with previous reports on this instrument platform [5,16].

Global view on CCS values of modified peptides
The fact that charge is a major determinant of ion mobility prompted us to investigate the occurrence of different charge states in our dataset.Figure 3a provides an overview of predominant charge states as well as the relative abundance fraction of charge states for all peptides.Tryptic peptides generally take up two to four protons in the electrospray process, depending on the length of the amino acid sequence and the number of basic residues.PTMs can alter the pK a and gas phase basicity of the modified amino acid and hence the charge state distribution.This effect was most striking for lysine modifications, as doubly and triply charged species appeared in roughly equal abundance for the unmodified peptides, whereas the various acylations shifted the charge distribution almost completely to charge 2. By contrast, lysine methylations as well as the GlyGly residue retain basic properties at the lysine site and thus showed only little effect on the relative abundance distribution, but rather tipped the predominant charge state to 3. We observed similar trends for arginine methylation and, as expected, citrullination reduced the charge state.Hydroxylation of proline, O-GlcNAcylation of serine and threonine as well as nitration and phosphorylation of tyrosine did not alter the charge state distribution as compared with their unmodified counterparts.Overall, these results are in agreement with the preceding analysis of the 21 PTM library on a different instrument platform [38].
To illustrate the global effect of charge state alterations on the ion mobility of modified peptides, we plotted them in their predominant charge state in the m/z vs. CCS space (Fig. 3b).The unmodified peptides of the lysine pool (left panel) are not strictly tryptic peptides because of their internal lysine.Nevertheless, they followed the well-characterized distribution of ion mobility and charge state occupancy.Malonylation, as an example for the group of acylations, thinned the population of charge 3 species and led to a more dense population of charge 2 peptides.Conversely, GlyGlycation caused a distinct shift in the opposite direction and predominantly populated the area of triply charged species.Thus, simple charge state alterations can already contribute to the separation of peptides and modifications in the ion mobility dimension.

Pairwise comparison of modified and unmodified peptides
A distinct feature of our dataset is the large number of matching pairs of modified and unmodified peptide sequences.Thus, having determined the position of modified peptide populations in the CCS vs. m/z space, we next performed a pairwise analysis of the respective counterparts.Intuitively, one could expect increasing CCS values throughout, as all modifications in our dataset introduce an additional functional group to the peptide (with the only exception of citrullination).However, we observed relative differences (∆CCS = (CCS modified − CCS unmodified )/CCS unmodified ) ranging all the way from about −10 to +10%, depending on the type of modification, the modified amino acid, and the charge state (Fig. 4a).As a proxy for the experimental precision in the dataset, we indicate the +/− three-fold coefficient of variation (0.78%) interval in the boxplot.The median ∆CCS shifts of specific modifications ranged from almost no difference (0% for lysine methylation) to slightly negative values To delve deeper into factors that determine CCS values of modified peptides, it is insightful to first consider the two different homologous series in the subset of lysine modifications: acylations with mono-(formyl to butyryl) and dicarboxylic acids (malonyl to glutaryl).Both series show a nearly linear increase in CCS with longer acyl chains, consistent for both charges 2 and 3 (Fig. 4a).This is in line with the intuition of increasing size with increasing chain length.To our surprise, and despite the fact that both series differ by one carboxylic acid, the median shifts in CCS of, for example, acetyl and malonyl were almost identical (−0.2%, −0.5% respectively).Moreover, they were similar to those of lysine methylation (0%), an even smaller moiety.A possible explanation for this observation are the different chemical properties of the modifications.While methylation tends to increase the basicity at the modified lysine, the acetyl group is electron withdrawing and malonylation adds an acidic functionality.Peptide conformations in the gas phase can be partially explained by Coulomb interactions and intramolecular charge solvation [14].Analyzing a larger dataset of peptide CCS values, Chang et al. concluded that internal acidic residues facilitate more compact conformations and, conversely, internal basic residues result in more extended conformations with larger cross sections [35].Our data is in line with this model and suggests that the carboxyl residue of malonyl compensates for its increased size.
Next, we plotted the median ∆CCS values as a function of the modification's molecular weight (Fig. 4b, Suppl.Fig. 3a).This analysis reproduced the large effect size for O-GlcNAcylation and biotinylation and showed an overall correlation of r 2 = 0.74 between ∆CCS and ∆M.However, when excluding the latter from the analysis, r 2 dropped to only 0.38, indicating that mass alone is a poor predictor of ∆CCS values for this group of chemically diverse modifications.In line with the results above, the homologous lysine modification series appear on parallel lines in this plot (dashed gray lines).In addition, we observed modifications with different molecular weight but similar ∆CCS such as the methyl-acetyl-malonyl example above (7,21,and 43 Da).Conversely, other modifications such as malonyl and hydroxyisobutyryl have a similar molecular weight but different ∆CCS (−0.5% vs. +1.4%).As the example of lysineand arginine-dimethylation shows, this effect is not limited to the chemistry of the modification itself, but can also depend on the modified amino acid.Performing a similar correlation analysis of ∆CCS and changes in the chromatographic retention time (Fig. 4c, Suppl.Fig. 3b) revealed some degree of orthogonality between LC and ion mobility.O-GlcNAcylation, phosphorylation, and GlyGlycation, for example, only slightly altered the retention time, while resulting in relatively large CCS shifts.
Although our analysis highlighted conceivable trends for each modification in our data, we also noted a relatively high variance within the modification groups (Fig. 4a).In most cases, the interquartile range of the pairwise analysis was several-fold larger than the precision of our experiments.Furthermore, some modifications showed less variance compared to others and the variance in charge 3 species was generally higher.This hints towards sequence-dependent effects within the peptide pools that modulate the effect of modifications on peptide cross sections.

Resolving sequence-dependent CCS determinants
To dissect the high variance within the modification groups, we next resolved our data by peptide sequences.This is possible because all peptides from one amino acid pool have the same unmodified base sequence.Figure 5a shows absolute ΔCCS values of all doubly charged peptides in the subset of lysine modifications color-coded by their respective base peptide sequence.Strikingly, the rank order of ΔCCS values remained largely unchanged throughout the homologous series of acyl modifications.In other words, while a particular modification could alter a peptide's cross section from +10 Å 2 to −25 Å 2 (formylation) depending on its amino acid sequence, elongating the modification, e.g., from formyl to butyryl, resulted in consistent increments across all peptide sequences.This observation also applied to triply charged peptide ions (Suppl.Fig. 4) and only few peptides deviated from this trend.Even for the chemically rather distant biotinylation, we observed a similar behavior, although more peptides swapped positions, in particular those with larger ΔCCS values.In contrast, when extending the line plot to lysine methylations and GlyGly, the trend was interrupted.However, the latter modifications are also distinct from the former acyl-type modifications in terms of charge state distribution (Fig. 3).This suggests that the observed grouping is driven by intramolecular charge localization and solvation.Similarly, we found that for individual sequences, the absolute ΔCCS as well as the shift relative to other sequences can be discordant for charge states (Fig. 5b).These results raised the question of whether there are commonalities in peptide sequences that undergo similar changes in their gas phase structure upon modification.To this end, we plotted the m/z vs. CCS distribution of the unmodified peptides and overlaid the corresponding ΔCCS values for different modifications.As a visual aid, we divided them into extended (larger CCS values) and compact structures (smaller CCS values) by fitting a linear model to each charge state.Figure 5c shows lysine succinylation and trimethylation as examples.For succinylation, close inspection of the doubly and triply charged ion populations revealed a tendency towards negative ΔCCS values for extended structures and vice versa.We also observed a compaction of extended structures for triply charged trimethylated peptides, indicating that the internally localized charge destabilizes extended conformations.This sequence dependency of trimethylation was less pronounced for doubly charged species, in line with their narrower ΔCCS distribution in Fig. 4a.
To investigate this further, we selected two peptide sequences from our lysine succinylation data and modeled the gas phase structures of their doubly charged ions using a computational approach based on molecular dynamics (Methods, Suppl.Fig. 5).Our modeling recapitulated the expected charge reduction at the modified lysine residue, which for the succinylated peptides resulted in proton localization at the peptide N-terminus and the C-terminal lysine residue.By contrast, one acidic proton was assumed to be sequestered at the internal lysine residues for the unmodified peptides.For "VGID…K(succinyl)LK," the proton localization to the terminal residues resulted in re-folding to a more extended structure as compared to its unmodified counterpart (+32 Å 2 , +6.7%), while the lysine succinylation prevented interactions between the terminal residues.This was in good agreement with our experiment (+32 Å 2 ).Conversely, for the modeled gas phase structure of the doubly charged "GTI… K(succinyl)…AK," proton relocalization and succinylation of the internal lysine did not prevent the adoption of folded peptide conformations by interaction of the terminal residues.Consequently, only minor changes of the gas phase collisional cross sections were computed (+8 Å 2 , +1.9%), while we even found a negative ΔCCS value (−25 Å 2 ) experimentally.
These results led us to hypothesize that the gas phase conformation of the unmodified peptide ion is indicative of the relative effect of specific modifications.To test this on all our data, we plotted the ΔCCS value for each modification as a function of the corresponding residuals of the linear fit (indicating whether the unmodified peptide adapts a more compact or extended structure) (Suppl.Figs. 6 and 7).The slope of the resulting linear trend lines can be interpreted as "ΔCCS gradients" within the CCS vs. m/z ion populations (Fig. 5d).Indeed, all lysine acylations and citrullination followed the trend described above for succinylation, while GlyGly and methylated peptides resembled trimethylation.In particular for doubly charged peptides with lysine acylations, the ion populations in CCS vs. m/z appeared narrower with respect to their trendlines (Supplementary Fig. 8).Taken together, our data suggests that the ΔCCS value associated with a modification is indeed correlated with the "starting conformation" of the unmodified peptide and hence at least partially dependent on the amino acid sequence, rather than a fixed increment determined by its chemical composition.

Conclusions
Recent advances in the application of ion mobility spectrometry to MS-based proteomics also promise new opportunities for the proteome-wide characterization of posttranslational modifications.In particular, the combination of TIMS and PASEF has enabled the precise measurement of CCS values on the scale of hundreds of thousands to more than a million data points and contributed to a better understanding of sequence-and position-dependent determinants of peptide cross sections [16,35,57].To extend this work beyond unmodified peptides, here, we used synthetic peptide libraries with known ground truth to characterize the effect of 22 different PTMs on peptide CCS values.
Our study provides data on 115 to 292 matched modified and unmodified peptides per modification with a precision <1% after linear alignment, which is on par with previous studies on the same instrument platform [5,16].Limitations of this approach include its reliance on the accuracy of the software-based feature detection as well as the fact that the accuracy of TIMS CCS measurements, particularly for complex samples, can be affected by exceeding the local charge capacity in the TIMS cartridge.However, the precision of the dataset enabled us to investigate different layers of modification-specific effects on peptide CCS values.On a global level, we observed major shifts in the m/z vs. ion mobility distribution for modified peptides, which we attributed to changes in their predominant charge state.In proteomics practice, such effects can be important, for example, to optimize the precursor selection scheme in dia-PASEF experiments [24,25] or to bias data-dependent acquisition towards modified peptides [27,58].
Turning to pairwise comparisons of modified peptides and their unmodified counterparts, we observed median ∆CCS values in a range of −1.4 to 4.8%.Surprisingly, despite the correlation between ion mass and mobility, the modification mass alone proved to be a poor predictor of ∆CCS values for most modifications in our dataset.In parts, we could rationalize these observations by the counteracting effects of increased modification size on the one hand and intramolecular charge solvation on the other hand.In addition, our data revealed substantial sequence-dependent effects on the cross section of modified peptides.This is in line with another recent study focused on phosphorylated peptides [26].All in all, our study adds to the increasing body of work indicating that peptide cross sections are determined by the amino acid composition [59,60] as well as their linear sequence [16,35,61].
Accurate predictions of peptide properties such as retention time, MS/MS spectra, and CCS values are increasingly used in MS-based proteomics [62,63].In this context, synthetic peptides can provide important training data as they have a known ground truth [37].This applies in particular to modified peptides, which are not always readily accessible from biological sources via efficient and affordable enrichment protocols.We envision that our high-quality dataset fills this gap and helps to extend CCS prediction algorithms to various post-translational modifications, for example, via transfer learning [36].

Fig. 1
Fig. 1 Overview of 22 post-translational modifications (PTMs) analyzed by trapped ion mobility spectrometry.A Schematic workflow for the analysis of synthetic peptide pools by liquid chromatography and trapped ion mobility-quadrupole time-of-flight mass spectrom-

Fig. 2
Fig. 2 Reproducibility and cross-run alignment of peptide collision cross sections.A CCS values of eleven reference peptides spiked into 84 consecutive LC-MS experiments.Data are shown in chronological order before (top) and after (below) performing a linear alignment to correct for drifts in the ion mobility measurement (Methods).B Coefficients of variation of CCS values of (modified) peptides measured in three replicate injections in the full dataset (top: raw data, bottom: aligned data).N = 4213, bin width = 0.05, 1% outliers out of bounds not shown

Fig. 3
Fig. 3 Analysis of charge state and collision cross section (CCS) distribution for different PTMs.A Number of identified peptides per modification in their predominant charge state (top).Relative abundance of different charge states for all identified peptides (bottom).For R and K, only unmodified peptide sequences with at least one internal R or K are shown.B Peptide m/z vs. CCS distribution of selected modifications from our dataset of 22 PTMs.The scatter plots show peptides in their predominant charge state

Fig. 4
Fig. 4 Pairwise comparison of modified and unmodified peptide collision cross section (CCS).A Relative difference of CCS values of matching pairs, separated by charge state (ΔCCS).The number of pairs is plotted on the top axis.The median relative deviation (MED) indicates the combined value for charges 2 and 3, which was tested for significance using a one-sample Wilcoxon test with Benjamini-Hochberg correction (*p<0.05;**p<0.01;***p<0.001;****p<1e-04).Dashed lines indicate the 3-fold median coefficient of variation in the full dataset as a proxy for experimental precision.Boxplot elements: Interquartile range within boxes; median indicated by horizontal line; whiskers spanning 1.5-fold the interquartile range.B Median ΔCCS values of all investigated modifications as a function of the molecular weight change (ΔM).Dashed lines visualize homologous series.See also Supplementary Fig. 3 for absolute ΔCCS and relative ΔM values.C Same as B, but as a function of the observed shift in retention time (ΔRT)

Fig. 5
Fig. 5 Sequence-dependent fine structure of PTM-induced CCS shifts.A ΔCCS values of individual peptide sequences for lysine modifications at charge 2. Only peptides which were detected for all shown modifications are included (N = 63).Peptides with the same unmodified sequence are connected and color-coded according to their rank order for formylation, from highest (pink) to lowest ΔCCS (orange).B Heat maps of ΔCCS values comparing charges 2 and 3 for both succinylation (N = 89) trimethylation (N = 102).C Peptide m/z vs. CCS distribution of unmodified reference peptides for succinylation (N = 244) and trimethylation (N = 246), overlaid with each peptides PTM-induced CCS shift.Linear trend lines are fitted to both charge states.D ΔCCS gradient within ion clouds for each investigated PTM, separated by charge.Each peptide's residual to a linear regression, as depicted in panel C, was plotted against its ΔCCS value.The slope of a subsequent linear regression represents a ΔCCS gradient within the respective ion clouds (panel C)

Andreas
Will is currently a doctoral researcher at Jena University Hospital jointly supervised by Prof. Florian Meier and Prof. Regine Heller.His research focuses on the application of trapped ion mobility-mass spectrometry to detect post-translational modifications, with a particular interest in the role of O-GlcNAcylation during aging and cellular senescence.Denys Oliinyk is currently a PhD student at Jena University Hospital in Prof. Florian Meier's group.His expertise is in the development of mass spectrometry-based methods to analyze protein phosphorylation networks in human health and disease at the systems scale.Christian Bleiholder is Associate Professor of Chemistry and Biochemistry at Florida State University with a D.Sc. in chemistry (University of Heidelberg & Ger man Cancer Research Center).His current research focuses on the development and application of ion mobility spectrometry/mass spectrometry methods.His awards include a fellowship from the Alexandervon-Humboldt Foundation and a CAREER award from the National Science Foundation.Florian Meier is Junior Professor and independent research group leader at Jena University Hospital with a PhD in biochemistry from Ludwig Maximilian University Munich.His interest is in the development of (ion mobility) mass spectrometry-based omics technologies and their applications in biomedical research.His awards include the 2019 Wolfgang Paul Study Award of the German Mass Spectrometry Society.