Comprehensive Peptide Ion Structure Studies Using Ion Mobility Techniques: Part 1. An Advanced Protocol for Molecular Dynamics Simulations and Collision Cross-Section Calculation

Ionmob: a Python package for prediction of peptide collisional cross-section values

MOTIVATION

Including ion mobility separation (IMS) into mass spectrometry proteomics experiments is useful to improve coverage and throughput. Many IMS devices enable linking experimentally derived mobility of an ion to its collisional cross-section (CCS), a highly reproducible physicochemical property dependent on the ion's mass, charge and conformation in the gas phase. Thus, known peptide ion mobilities can be used to tailor acquisition methods or to refine database search results. The large space of potential peptide sequences, driven also by posttranslational modifications of amino acids, motivates an in silico predictor for peptide CCS. Recent studies explored the general performance of varying machine-learning techniques, however, the workflow engineering part was of secondary importance. For the sake of applicability, such a tool should be generic, data driven, and offer the possibility to be easily adapted to individual workflows for experimental design and data processing.

RESULTS

We created ionmob, a Python-based framework for data preparation, training, and prediction of collisional cross-section values of peptides. It is easily customizable and includes a set of pretrained, ready-to-use models and preprocessing routines for training and inference. Using a set of ≈21 000 unique phosphorylated peptides and ≈17 000 MHC ligand sequences and charge state pairs, we expand upon the space of peptides that can be integrated into CCS prediction. Lastly, we investigate the applicability of in silico predicted CCS to increase confidence in identified peptides by applying methods of re-scoring and demonstrate that predicted CCS values complement existing predictors for that task.

AVAILABILITY AND IMPLEMENTATION

The Python package is available at github: https://github.com/theGreatHerrLebert/ionmob.

Formation of Gaseous Peptide Ions from Electrospray Droplets: Competition between the Ion Evaporation Mechanism and Charged Residue Mechanism

The transfer of peptide ions from solution into the gas phase by electrospray ionization (ESI) is an integral component of mass spectrometry (MS)-based proteomics. The mechanisms whereby gaseous peptide ions are released from charged ESI nanodroplets remain unclear. This is in contrast to intact protein ESI, which has been the focus of detailed investigations using molecular dynamics (MD) simulations and other methods. Under acidic liquid chromatography/MS conditions, many peptides carry a solution charge of 3+ or 2+. Because of this pre-existing charge and their relatively small size, prevailing views suggest that peptides follow the ion evaporation mechanism (IEM). The IEM entails analyte ejection from ESI droplets, driven by electrostatic repulsion between the analyte and droplet. Surprisingly, recent peptide MD investigations reported a different behavior, that is, the release of peptide ions via droplet evaporation to dryness which represents the hallmark of the charged residue mechanism (CRM). Here, we resolved this conundrum by performing MD simulations on a common model peptide (bradykinin) in Rayleigh-charged aqueous droplets. The primary focus was on pH 2 conditions (bradykinin solution charge = 3+), but we also verified that our MD strategy captured pH-dependent charge state shifts seen in ESI-MS experiments. In agreement with earlier simulations, we found that droplets with initial radii of 1.5-3 nm predominantly release peptide ions via the CRM. In contrast, somewhat larger radii (4-5 nm) favor IEM behavior. It appears that these are the first MD data to unequivocally demonstrate the viability of peptide IEM events. Electrostatic arguments can account for the observed droplet size dependence. In summary, both CRM and IEM can be operative in peptide ESI-MS. The prevalence of one over the other mechanism depends on the droplet size distribution in the ESI plume.

Weight average approaches for predicting dynamical properties of biomolecules

Recent advances in atomistic molecular dynamics (MD) simulations of biomolecules allow us to explore their conformational spaces widely, observing large-scale conformational fluctuations or transitions between distinct structures. To reproduce or refine experimental data using MD simulations, structure ensembles, which are characterized by multiple structures and their statistical weights on the rugged free-energy landscapes, are often used. Here, we summarize weight average approaches for various experimental measurements. Weight average approaches are now applied to hybrid quantum mechanics/molecular mechanics MD simulations to predict fast vibrational motions in a protein with a high accuracy for better understanding of molecular functions from atomic structures.

Sequence-Specific Model for Predicting Peptide Collision Cross Section Values in Proteomic Ion Mobility Spectrometry

The contribution of peptide amino acid sequence to collision cross section values (CCS) has been investigated using a dataset of ∼134 000 peptides of four different charge states (1+ to 4+). The migration data were acquired using a two-dimensional liquid chromatography (LC)/trapped ion mobility spectrometry/quadrupole/time-of-flight mass spectrometry (MS) analysis of HeLa cell digests created using seven different proteases and was converted to CCS values. Following the previously reported modeling approaches using intrinsic size parameters (ISP), we extended this methodology to encode the position of individual residues within a peptide sequence. A generalized prediction model was built by dividing the dataset into eight groups (four charges for both tryptic/nontryptic peptides). Position-dependent ISPs were independently optimized for the eight subsets of peptides, resulting in prediction accuracy of ∼0.981 for the entire population of peptides. We find that ion mobility is strongly affected by the peptide's ability to solvate the positively charged sites. Internal positioning of polar residues and proline leads to decreased CCS values as they improve charge solvation; conversely, this ability decreases with increasing peptide charge due to electrostatic repulsion. Furthermore, higher helical propensity and peptide hydrophobicity result in a preferential formation of extended structures with higher than predicted CCS values. Finally, acidic/basic residues exhibit position-dependent ISP behavior consistent with electrostatic interaction with the peptide macrodipole, which affects the peptide helicity. The MS raw data files have been deposited with the ProteomeXchange Consortium via the jPOST partner repository (http://jpostdb.org) with the dataset identifiers PXD021440/JPST000959, PXD022800/JPST001017, and PXD026087/ JPST001176.

Sub-ambient pressure IR-MALDI ion mobility spectrometer for the determination of low and high field mobilities

A new ion mobility (IM) spectrometer, enabling mobility measurements in the pressure range between 5 and 500 mbar and in the reduced field strength range E/N of 5-90 Td, was developed and characterized. Reduced mobility (K₀) values were studied under low E/N (constant value) as well as high E/N (deviation from low field K₀) for a series of molecular ions in nitrogen. Infrared matrix-assisted laser desorption ionization (IR-MALDI) was used in two configurations: a source working at atmospheric pressure (AP) and, for the first time, an IR-MALDI source working with a liquid (aqueous) matrix at sub-ambient/reduced pressure (RP). The influence of RP on IR-MALDI was examined and new insights into the dispersion process were gained. This enabled the optimization of the IM spectrometer for best analytical performance. While ion desolvation is less efficient at RP, the transport of ions is more efficient, leading to intensity enhancement and an increased number of oligomer ions. When deciding between AP and RP IR-MALDI, a trade-off between intensity and resolving power has to be considered. Here, the low field mobility of peptide ions was first measured and compared with reference values from ESI-IM spectrometry (at AP) as well as collision cross sections obtained from molecular dynamics simulations. The second application was the determination of the reduced mobility of various substituted ammonium ions as a function of E/N in nitrogen. The mobility is constant up to a threshold at high E/N. Beyond this threshold, mobility increases were observed. This behavior can be explained by the loss of hydrated water molecules.

Investigating the interactions of the first 17 amino acid residues of Huntingtin with lipid vesicles using mass spectrometry and molecular dynamics

The first 17 amino acid residues of Huntingtin protein (Nt17 of htt) are thought to play an important role in the protein's function; Nt17 is one of two membrane binding domains in htt. In this study the binding ability of Nt17 peptide with vesicles comprised of two subclasses of phospholipids is studied using electrospray ionization - mass spectrometry (ESI-MS) and molecular dynamics (MD) simulations. Overall, the peptide is shown to have a greater propensity to interact with vesicles of phosphatidylcholine (PC) rather than phosphatidylethanolamine (PE) lipids. Mass spectra show an increase in lipid-bound peptide adducts where the ordering of the number of such specie is 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) > 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) > 1-palmitoyl-2-oleoyl-sn-glycero-3 phosphoethanolamine (POPE). MD simulations suggest that the compactness of the bilayer plays a role in governing peptide interactions. The peptide shows greater disruption of the DOPC bilayer order at the surface and interacts with the hydrophobic tails of lipid molecules via hydrophobic residues. Conversely, the POPE vesicle remains ordered and lipids display transient interactions with the peptide through the formation of hydrogen bonds with hydrophilic residues. The POPC system displays intermediate behavior with regard to the degree of peptide-membrane interaction. Finally, the simulations suggest a helix stabilizing effect resulting from the interactions between hydrophobic residues and the lipid tails of the DOPC bilayer.

Comparison of Peptide Ion Conformers Arising from Non-Helical and Helical Peptides Using Ion Mobility Spectrometry and Gas-Phase Hydrogen/Deuterium Exchange

The dominant gas-phase conformer of [M+3H]³⁺ ions of the model peptide acetyl-PSSSSKSSSSKSSSSKSSSSK has been examined with ion mobility spectrometry (IMS), gas-phase hydrogen deuterium exchange (HDX), and mass spectrometry (MS) techniques. The [M+3H]³⁺ peptide ions are observed predominantly as a relatively compact conformer type. Upon subjecting these ions to electron transfer dissociation (ETD), the level of protection for each amino acid residue in the peptide sequence is assessed. The overall per-residue deuterium uptake is observed to be relatively more efficient for the neutral residues than for the model peptide acetyl-PAAAAKAAAAKAAAAKAAAAK. In comparison, the N-terminal and C-terminal regions of the serine peptide show greater relative protection compared with interior residues. Molecular dynamics (MD) simulations have been used to generate candidate structures for collision cross section and HDX reactivity matching. Hydrogen accessibility scoring (HAS) for select structural candidates from MD simulations has been used to suggest conformer types that could contribute to the observed HDX patterns. The results are discussed with respect to recent studies employing extensive MD simulations of gas-phase structure establishment of a peptide system. Graphical Abstract ᅟ.

Ion mobility mass spectrometry with molecular modelling to reveal bioactive isomer conformations and underlying relationship with isomerization

RATIONALE

In medicine and drug development, molecular modelling is an important tool. It is attractive to develop a platform connecting the theoretical structural modelling and the results from experimental measurement. In addition, the separation and structural analysis of bioactive constituent isomers are still challenging tasks.

METHODS

Drift tube ion mobility (IM) mass spectrometry (MS) provides the experimental collision cross section (CCS) which contains the structural information. The experimental CCS can be compared with the calculated CCS of the molecular modelling structures. This technique is especially useful for bioactive constituents in herbal medicine because active isomers with the same chemical formula are common in these samples. IM helps separate and identify these isomers and reveals details about their structures and conformations.

RESULTS

Two model bioactive constituents, caffeoylquinic acids (CQAs) and dicaffeoylquinic acids (di-CQAs), were selected to systematically investigate the influence of solution, ion source conditions and ion heating on the isomer CCS distributions. By comparing the calculated CCS with the experimental value, we identified the favorable conformations of CQAs. The most compact conformation of a CQA was less likely to isomerize than the more extended conformation. It was found that the isomerization tendency was in accord with the conformation favorability.

CONCLUSIONS

This study offers an effective approach to predict and demystify the conformation and isomerization of the active constituents in herbal medicines.

Comprehensive Peptide Ion Structure Studies Using Ion Mobility Techniques: Part 3. Relating Solution-Phase to Gas-Phase Structures

Molecular dynamics (MD) simulations have been utilized to study peptide ion conformer establishment during the electrospray process. An explicit water model is used for nanodroplets containing a model peptide and hydronium ions. Simulations are conducted at 300 K for two different peptide ion charge configurations and for droplets containing varying numbers of hydronium ions. For all conditions, modeling has been performed until production of the gas-phase ions and the resultant conformers have been compared to proposed gas-phase structures. The latter species were obtained from previous studies in which in silico candidate structures were filtered according to ion mobility and hydrogen-deuterium exchange (HDX) reactivity matches. Results from the present study present three key findings namely (1) the evidence from ion production modeling supports previous structure refinement studies based on mobility and HDX reactivity matching, (2) the modeling of the electrospray process is significantly improved by utilizing initial droplets existing below but close to the calculated Rayleigh limit, and (3) peptide ions in the nanodroplets sample significantly different conformers than those in the bulk solution due to altered physicochemical properties of the solvent. Graphical Abstract ᅟ.

Magnifying ion mobility spectrometry-mass spectrometry measurements for biomolecular structure studies

Ion mobility spectrometry-mass spectrometry (IMS-MS) provides information about the structures of gas-phase ions in the form of a collision cross section (CCS) with a neutral buffer gas. Indicating relative ion size, a CCS value alone is of limited utility. Although such information can be used to propose different conformer types, finer details of structure are not captured. The increased accessibility of IMS-MS measurements with commercial instrumentation in recent years has ballooned its usage in combination with separate measurements to provide enhanced data from which greater structural inferences can be drawn. This short review presents recent outstanding developments in scientific research that employs complementary measurements that when combined with IMS-MS data are used to characterize the structures of a wide range of compounds.

Adding a new separation dimension to MS and LC-MS: What is the utility of ion mobility spectrometry?

Ion mobility spectrometry is an analytical technique known for more than 100 years, which entails separating ions in the gas phase based on their size, shape, and charge. While ion mobility spectrometry alone can be useful for some applications (mostly security analysis for detecting certain classes of narcotics and explosives), it becomes even more powerful in combination with mass spectrometry and high-performance liquid chromatography. Indeed, the limited resolving power of ion mobility spectrometry alone can be tackled when combining this analytical strategy with mass spectrometry or liquid chromatography with mass spectrometry. Over the last few years, the hyphenation of ion mobility spectrometry to mass spectrometry or liquid chromatography with mass spectrometry has attracted more and more interest, with significant progresses in both technical advances and pioneering applications. This review describes the theoretical background, available technologies, and future capabilities of these techniques. It also highlights a wide range of applications, from small molecules (natural products, metabolites, glycans, lipids) to large biomolecules (proteins, protein complexes, biopharmaceuticals, oligonucleotides).

Comprehensive Gas-Phase Peptide Ion Structure Studies Using Ion Mobility Techniques: Part 2. Gas-Phase Hydrogen/Deuterium Exchange for Ion Population Estimation

Gas-phase hydrogen/deuterium exchange (HDX) using D2O reagent and collision cross-section (CCS) measurements are utilized to monitor the ion conformers of the model peptide acetyl-PAAAAKAAAAKAAAAKAAAAK. The measurements are carried out on a home-built ion mobility instrument coupled to a linear ion trap mass spectrometer containing electron transfer dissociation (ETD) capabilities. ETD is utilized to obtain per-residue deuterium uptake data for select ion conformers, and a new algorithm is presented for interpreting the HDX data. Using molecular dynamics (MD) production data and a hydrogen accessibility scoring (HAS)-number of effective collisions (NEC) model, hypothetical HDX behavior is attributed to various in-silico candidate (CCS match) structures. The HAS-NEC model is applied to all candidate structures, and non-negative linear regression is employed to determine structure contributions resulting in the best match to deuterium uptake. The accuracy of the HAS-NEC model is tested with the comparison of predicted and experimental isotopic envelopes for several of the observed c-ions. It is proposed that gas-phase HDX can be utilized effectively as a second criterion (after CCS matching) for filtering suitable MD candidate structures. In this study, the second step of structure elucidation, 13 nominal structures were selected (from a pool of 300 candidate structures) and each with a population contribution proposed for these ions. Graphical Abstract ᅟ.