Experimental and Theoretical Investigation into the Correlation between Mass and Ion Mobility for Choline and Other Ammonium Cations in N2

Structural studies of supramolecular complexes and assemblies by ion mobility mass spectrometry

Recent advances in instrumentation and development of computational strategies for ion mobility mass spectrometry (IM-MS) studies have contributed to an extensive growth in the application of this analytical technique to comprehensive structural description of supramolecular systems. Apart from the benefits of IM-MS for interrogation of intrinsic properties of noncovalent aggregates in the experimental gas-phase environment, its merits for the description of native structural aspects, under the premises of having maintained the noncovalent interactions innate upon the ionization process, have attracted even more attention and gained increasing interest in the scientific community. Thus, various types of supramolecular complexes and assemblies relevant for biological, medical, material, and environmental sciences have been characterized so far by IM-MS supported by computational chemistry. This review covers the state-of-the-art in this field and discusses experimental methods and accompanying computational approaches for assessing the reliable three-dimensional structural elucidation of supramolecular complexes and assemblies by IM-MS.

UV-A-Induced Photoisomerization and Photodimerization of Curcumin: An Ion Mobility Mass Spectrometry Study

Curcumin, the bioactive compound present in spice plant turmeric, has been shown to exhibit selective phototoxic activities toward mammalian cancer cells, and it is being used extensively as a photosensitizer (PS) in photodynamic therapies (PDT). However, so far, the fate of curcumin toward photochemical transformations is not well understood. Here we report our findings of a number of novel photochemical reaction channels of curcumin in water-methanol mixture, like photoisomerization, photodimerization, and photooxidation (H2-loss). The reaction was performed by irradiating the curcumin solution with ultraviolet (UV) light of wavelength 350 nm, which is abundant in the earth's troposphere. Product identification and structure elucidation are done by employing an integrated method of drift tube ion mobility mass spectrometry (DTIMS) in combination with high-performance liquid chromatography (HPLC) and collision-induced dissociation (CID) of the mass-selected molecular ions. Two photoisomers of curcumin produced as a result of trans-cis configurational changes about C═C double bonds in the excited state have been identified, and it has been shown that they could serve as the precursors for formation of isomeric dimers via [2 + 2] cycloaddition and H2-loss products. Comparisons of the experimentally measured collision cross-section (CCS) values of the reactant and product ions obtained by the DTIMS method with those predicted by the electronic structure theory are found to be very effective for the discrimination of the produced photoisomers. The observed photochemical reaction channels are potentially significant toward uses of curcumin as a photosensitizer in photodynamic therapy.

Do d(GCGAAGC) Cations Retain the Hairpin Structure in the Gas Phase? A Cyclic Ion Mobility Mass Spectrometry and Density Functional Theory Computational Study

d(GCGAAGC) is the smallest oligonucleotide with a well-defined hairpin structure in solution. We report a study of multiply protonated d(GCGAAGC) and its sequence-scrambled isomers, d(CGAAGCG), d(GCGAACG), and d(CGGAAGC), that were produced by electrospray ionization with the goal of investigating their gas-phase structures and dissociations. Cyclic ion mobility measurements revealed that dications of d(GCGAAGC) as well as the scrambled-sequence ions were mixtures of protomers and/or conformers that had collision cross sections (CCS) within a 439-481 Å2 range. Multiple ion conformers were obtained by electrospray under native conditions as well as from aqueous methanol. Arrival time distribution profiles were characteristic of individual isomeric heptanucleotides. Extensive Born-Oppenheimer molecular dynamics (BOMD) and density functional theory (DFT) calculations of d(GCGAAGC)2+ isomers indicated that hairpin structures were high-energy isomers of more compact distorted conformers. Protonation caused a break up of the C2···G6 pair that was associated with the formation of strong hydrogen bonds in zwitterionic phosphate anion-nucleobase cation motifs that predominated in low energy ions. Multiple components were also obtained for d(GCGAAGC)3+ trications under native and denaturing electrospray conditions. The calculated trication structures showed disruption of the G···C pairs in low energy zwitterions. A hairpin trication was calculated to be a high energy isomer. d(GCGAAGC)4+ tetracations were produced and separated by c-IMS as two major isomers. All low energy d(GCGAAGC)4+ ions obtained by DFT geometry optimizations were zwitterions in which all five purine bases were protonated, and the ion charge was balanced by a phosphate anion. Tetracations of the scrambled sequences were each formed as one dominant isomer. The CCS calculated with the MobCal-MPI method were found to closely match experimental values. Collision-induced dissociation (CID) spectra of multiply charged heptanucleotides showed nucleobase loss and backbone cleavages occurring chiefly at the terminal nucleosides. Electron-transfer-CID tandem mass spectra were used to investigate dissociations of different charge and spin states of charge-reduced heptanucleotide cation radicals.

MobCal-MPI 2.0: an accurate and parallelized package for calculating field-dependent collision cross sections and ion mobilities

Ion mobility spectrometry (IMS), which can be employed as either a stand-alone instrument or coupled to mass spectrometry, has become an important tool for analytical chemistry. Because of the direct relation between an ion's mobility and its structure, which is intrinsically related to its collision cross section (CCS), IMS techniques can be used in tandem with computational tools to elucidate ion geometric structure. Here, we present MobCal-MPI 2.0, a software package that demonstrates excellent accuracy (RMSE 2.16%) and efficiency in calculating low-field CCSs via the trajectory method (≤30 minutes on 8 cores for ions with ≤70 atoms). MobCal-MPI 2.0 expands on its predecessor by enabling the calculation of high-field mobilities through the implementation of the 2^nd order approximation to two-temperature theory (2TT). By further introducing an empirical correction to account for deviations between 2TT and experiment, MobCal-MPI 2.0 can compute accurate high-field mobilities that exhibit a mean deviation of <4% from experimentally measured values. Moreover, the velocities used to sample ion-neutral collisions were updated from a weighted to a linear grid, enabling the near-instantaneous evaluation of mobility/CCS at any effective temperature from a single set of N₂ scattering trajectories. Several enhancements made to the code are also discussed, including updates to the statistical analysis of collision event sampling and benchmarking of overall performance.

Exploring Isomerism in Isolated Cyclodextrin Oligomers through Trapped Ion Mobility Mass Spectrometry

Cyclodextrin (CD) macrocycles are used to create a wide range of supramolecular architectures which are also of interest in applications such as selective gas adsorption, drug delivery, and catalysis. However, predicting their assemblies and identifying the possible isomers in CD oligomers have always remained challenging due to their dynamic nature. Herein, we interacted CDs (α, β, and γ) with a divalent metal ion, Cu²⁺, to create a series of Cu²⁺-linked CD oligomers, from dimers to pentamers. We characterized these oligomers using electrospray ionization mass spectrometry and probed isomerism in each of these isolated oligomers using high resolution trapped ion mobility spectrometry. Using this technique, we separated multiple isomers for each of the Cu²⁺-interlinked CD oligomers and estimated their relative population, which was not accessible previously using other characterization techniques. We further carried out structural analysis of the observed isomers by comparing the experimental collision cross sections (CCSs) to that of modeled structures. We infer that the isomeric heterogeneity reflects size-specific packing patterns of individual CDs (e.g., close-packed/linear). In some cases, we also reveal the existence of kinetically trapped structures in the gas phase and study their transformation to thermodynamically controlled forms by examining the influence of activation of the ions on isomer interconversion.

Protonated Forms of Naringenin and Naringenin Chalcone: Proteiform Bioactive Species Elucidated by IRMPD Spectroscopy, IMS, CID-MS, and Computational Approaches

Naringenin (Nar) and its structural isomer, naringenin chalcone (ChNar), are two natural phytophenols with beneficial health effects belonging to the flavonoids family. A direct discrimination and structural characterization of the protonated forms of Nar and ChNar, delivered into the gas phase by electrospray ionization (ESI), was performed by mass spectrometry-based methods. In this study, we exploit a combination of electrospray ionization coupled to (high-resolution) mass spectrometry (HR-MS), collision-induced dissociation (CID) measurements, IR multiple-photon dissociation (IRMPD) action spectroscopy, density functional theory (DFT) calculations, and ion mobility-mass spectrometry (IMS). While IMS and variable collision-energy CID experiments hardly differentiate the two isomers, IRMPD spectroscopy appears to be an efficient method to distinguish naringenin from its related chalcone. In particular, the spectral range between 1400 and 1700 cm^-1 is highly specific in discriminating between the two protonated isomers. Selected vibrational signatures in the IRMPD spectra have allowed us to identify the nature of the metabolite present in methanolic extracts of commercial tomatoes and grapefruits. Furthermore, comparisons between experimental IRMPD and calculated IR spectra have clarified the geometries adopted by the two protonated isomers, allowing a conformational analysis of the probed species.

Energy-Resolved Ion Mobility Spectrometry: Composite Breakdown Curves for Distinguishing Isomeric Product Ions

Identification of lipopeptides (LpAA) synthesized from bacteria involves the study of structural characterization. Twenty LpAA have been characterized using commercial tandem high-resolution mass spectrometers in negative electrospray, employing nonresonant excitation in "RF only" collision cells and generally behave identically. However, [LpAA-H]^- (AA = Asp or Glu) shows surprising fragmentation pathways, yielding a complementary fatty acid carboxylate and dehydrated amino acid fragment anions. In this study, the dissociation mechanisms of [C12Glu-H]^- were determinate using energy-resolved mass spectrometry (ERMS). Product ion breakdown profiles are, generally, unimodal with full width at half-maximum (fwhm) increasing as product ion m/z ratios decrease, except for the two product ions of interest (fatty acid carboxylate and dehydrated glutamate) characterized by broad and composite profiles. Such behavior was already shown for other ions using a custom-built guided ion beam mass spectrometer. In this study, we investigate the meaning of these particular profiles from an ERMS breakdown, using fragmentation mechanisms depending on the collision energy. ERMS on line with ion mobility spectrometry (IMS), here called ER-IMS, provides a way to probe such questions. Broad or composite profiles imply that the corresponding product ions may be generated by two (or more) pathways, resulting in common or isomeric product ion structures. ER-IMS analysis indicates that the fatty acid carboxylate product ion is produced with a common structure through different pathways, while dehydrated glutamate has two isomeric forms depending on the mechanism involved. Drift time values correlate with the calculated collision cross section that confirms the product ion structures and fragmentation mechanisms.

The hitchhiker's guide to dynamic ion-solvent clustering: applications in differential ion mobility spectrometry

This article highlights the fundamentals of ion-solvent clustering processes that are pertinent to understanding an ion's behaviour during differential mobility spectrometry (DMS) experiments. We contrast DMS with static-field ion mobility, where separation is affected by mobility differences under the high-field and low-field conditions of an asymmetric oscillating electric field. Although commonly used in mass spectrometric (MS) workflows to enhance signal-to-noise ratios and remove isobaric contaminants, the chemistry and physics that underpins the phenomenon of differential mobility has yet to be fully fleshed out. Moreover, we are just now making progress towards understanding how the DMS separation waveform creates a dynamic clustering environment when the carrier gas is seeded with the vapour of a volatile solvent molecule (e.g., methanol). Interestingly, one can correlate the dynamic clustering behaviour observed in DMS experiments with gas-phase and solution-phase molecular properties such as hydrophobicity, acidity, and solubility. However, to create a generalized, global model for property determination using DMS data one must employ machine learning. In this article, we provide a first-principles description of differential ion mobility in a dynamic clustering environment. We then discuss the correlation between dynamic clustering propensity and analyte physicochemical properties and demonstrate that analytes exhibiting similar ion-solvent interactions (e.g., charge-dipole) follow well-defined trends with respect to DMS clustering behaviour. Finally, we describe how supervised machine learning can be used to create predictive models of molecular properties using DMS data. We additionally highlight open questions in the field and provide our perspective on future directions that can be explored.

Tautomers and Rotamers of Curcumin: A Combined UV Spectroscopy, High-Performance Liquid Chromatography, Ion Mobility Mass Spectrometry, and Electronic Structure Theory Study

The structures of tautomers and rotameric forms of curcumin, the bioactive compound present in spice plant turmeric, have been investigated using ion mobility mass spectrometry (IMMS) in conjunction with high-performance liquid chromatography (HPLC) and UV-visible spectroscopy. Two tautomeric forms of this β-diketone compound, keto-enol and diketo, have been chromatographically separated, and the electronic absorption spectra for these two tautomeric forms in methanol solution have been recorded separately for the first time. The molecular identity of the HPLC-separated solution fractions is established unambiguously by recording the mass and fragmentation spectra simultaneously. The ion mobility spectrum for the deprotonated curcumin anion, [Cur-H]^-, corresponding to the diketo tautomer, displays only one peak (P), and the collision cross-section (CCS) value is measured to be 185.9 Ų. However, the ion mobility spectrum corresponding to the HPLC-separated keto-enol tautomer shows three distinctly separated peaks, P, Q, and R, with CCS values of 185.9, 194.8, and 203.4 Ų, respectively, whereby peak R appears to be the most intense one, followed by peaks P and Q. The theoretically calculated CCS values of different isomers of [Cur-H]^-, optimized by electronic structure theory methods, display satisfactory correlation with the experimentally observed values, corroborating our assignments. The spectral attributes also indicate the occurrence of structural rearrangements in the electrospray ionization process. With the aid of the electronic structure calculation, low-energy pathways for the occurrence of the structural isomerization to surpass the energy barrier are suggested, which are consistent with the assignments of the peaks observed in the IM spectra.

Molecular Basis for the Remarkably Different Gas-Phase Behavior of Deprotonated Thyroid Hormones Triiodothyronine (T3) and Reverse Triiodothyronine (rT3): A Clue for Their Discrimination?

Thyroid hormones are biologically active small molecules responsible for growth and development regulation, basal metabolic rate, and lipid and carbohydrate metabolism. Liquid chromatography mass spectrometry (LC–MS) can be used to quantify thyroid hormones blood level with high speed and selectivity, aiming to improve the diagnosis and treatment of the severe pathological conditions in which they are implicated, i.e., hypo- and hyperthyroidism. In this work, the gas-phase behavior of the isomeric thyroid hormones triiodothyronine (T3) and reverse triiodothyronine (rT3) in their deprotonated form was studied at a molecular level using MS-based techniques. Previously reported collision-induced dissociation experiments yielded distinct spectra despite the high structural similarity of the two compounds, suggesting different charge sites to be responsible. Infrared multiple photon dissociation spectroscopy on [T3-H]⁻ and [rT3-H]⁻ was performed, and the results were interpreted using DFT and MP2 calculations, assessing the prevalence of T3 in the carboxylate form and rT3 as a phenolate isomer. The different deprotonation sites of the two isomers were also found to drive their ion-mobility behavior. In fact, [T3-H]⁻ and [rT3-H]⁻ were successfully separated. Drift times were correlated with collisional cross section values of 209 and 215 Ų for [T3-H]⁻ and [rT3-H]⁻, respectively. Calculations suggested the charge site to be the main parameter involved in the different mobilities of the two anions. Finally, bare [T3-H]⁻ and [rT3-H]⁻ were made to react with neutral acetylacetone and trifluoroacetic acid, confirming rT3 to be more acidic than T3 in agreement with the calculated gas-phase acidities of T3 and rT3 equal to 1345 and 1326 kJ mol^–1, respectively.

Current status and future prospects for ion-mobility mass spectrometry in the biopharmaceutical industry

Detailed characterization of protein reagents and biopharmaceuticals is key in defining successful drug discovery campaigns, aimed at bringing molecules through different discovery stages up to development and commercialization. There are many challenges in this process, with complex and detailed analyses playing paramount roles in modern industry. Mass spectrometry (MS) has become an essential tool for characterization of proteins ever since the onset of soft ionization techniques and has taken the lead in quality assessment of biopharmaceutical molecules, and protein reagents, used in the drug discovery pipeline. MS use spans from identification of correct sequences, to intact molecule analyses, protein complexes and more recently epitope and paratope identification. MS toolkits could be incredibly diverse and with ever evolving instrumentation, increasingly novel MS-based techniques are becoming indispensable tools in the biopharmaceutical industry. Here we discuss application of Ion Mobility MS (IMMS) in an industrial setting, and what the current applications and outlook are for making IMMS more mainstream.

Modular Ion Mobility Calibrants for Organometallic Anions Based on Tetraorganylborate Salts

Organometallics are widely used in catalysis and synthesis. Their analysis relies heavily on mass spectrometric methods, among which traveling-wave ion mobility spectrometry (TWIMS) has gained increasing importance. Collision cross sections (CCS) obtainable by TWIMS significantly aid the structural characterization of ions in the gas phase, but for organometallics, their accuracy has been limited by the lack of appropriate calibrants. Here, we propose tetraorganylborates and their alkali-metal bound oligomers [M_n-1(BR₄)_n]^- (M = Li, Na, K, Rb, Cs; R = aryl, Et; n = 1-6) as calibrants for electrospray ionization (ESI) TWIMS. These species chemically resemble typical organometallics and readily form upon negative-ion mode ESI of solutions of alkali-metal tetraorganylborates. By combining different tetraorganylborate salts, we have generated a large number of anions in a modular manner and determined their CCS values by drift-tube ion mobility spectrometry (DTIMS) (^DTCCS_He = 81-585, ^DTCCS_N2 = 130-704 Ų). In proof-of-concept experiments, we then applied these ^DTCCS values to the calibration of a TWIMS instrument and analyzed phenylcuprate and argentate anions, [Li_n-1M_nPh_2n]^- and [M_nPh_n+1]^- (M = Cu, Ag), as prototypical reactive organometallics. The ^TWCCS_N2 values derived from TWIMS measurements are in excellent agreement with those determined by DTIMS (<2% relative difference), demonstrating the effectiveness of the proposed calibration scheme. Moreover, we used theoretical methods to predict the structures and CCS values of the anions considered. These predictions are in good agreement with the experimental results and give further insight into the trends governing the assembly of tetraorganylborate, cuprate, and argentate oligomers.

Probing d- and l-Adrenaline Binding to β2-Adrenoreceptor Peptide Motifs by Gas-Phase Photodissociation Cross-Linking and Ion Mobility Mass Spectrometry

Diazirine-tagged d- and l-adrenaline derivatives formed abundant noncovalent gas-phase ion complexes with peptides N-Ac-SSIVSFY-NH₂ (peptide S) and N-Ac-VYILLNWIGY-NH₂ (peptide V) upon electrospray ionization. These peptide sequences represent the binding motifs in the β₂-adrenoreceptor. The structures of the gas-phase complexes were investigated by selective laser photodissociation of the diazirine chromophore at 354 nm, which resulted in a loss of N₂ and formation of a transient carbene intermediate in the adrenaline ligand without causing its expulsion. The photolyzed complexes were analyzed by collision-induced dissociation (CID-MS³ and CID-MS⁴) in an attempt to detect cross-links and establish the binding sites. However, no cross-linking was detected in the complexes regardless of the peptide and d- or l-configuration in adrenaline. Cyclic ion mobility measurements were used to obtain collision cross sections (CCS) in N₂ for the peptide S complexes. These showed identical values, 334 ± 0.9 Ų, for complexes of the l- and d-adrenaline derivatives, respectively. Identical CCS were also obtained for peptide S complexes with natural l- and d-adrenaline, 317 ± 1.2 Ų, respectively. Born-Oppenheimer molecular dynamics (BOMD) in combination with full geometry optimization by density functional theory calculations provided structures for the complexes that were used to calculate theoretical CCS with the ion trajectory method. A close match (337 Ų) was found for a single low Gibbs energy structure that displayed a binding pocket with Ser 2 and Ser 5 residues forming hydrogen bonds to the adrenaline catechol hydroxyls. Analysis of the BOMD trajectories revealed a small number of contacts between the incipient carbene carbon atom in the ligand and X-H bonds in the peptide, which was consistent with the lack of cross-linking. Temperature dependence of the internal dynamics of peptide S-adrenaline complexes as well as the specifics of the adrenaline carbene reactions are discussed. In particular, peptide amide hydrogen transfer to the carbene carbon atom was calculated to require crossing a potential energy barrier, which may hamper cross-linking in competition with carbene internal rearrangements.

Reevaluating the Role of Polarizability in Ion Mobility Spectrometry

With the expanding commercial availability of gas-phase separation systems that incorporate gas-phase mobility, there is a concurrent rise in efforts to cast the gas-phase mobility coefficient in terms of an ion-neutral collision cross-section (CCS). The motivating factors for this trend are varied, but many aim to complement experimental results with computationally generated CCS values from in silico structural approximations. Unfortunately, the current paradigm for relating experimental mobility results to computationally derived structures relies upon empirical approaches, including a myriad of variables that do not realistically bound the comparison. In this Critical Insight, we advocate for the development of a self-consistent experimental and computational framework that uses laboratory results to constrain the scope of the modeling effort. This paper aims to prompt discussion, challenge assumptions, and promote the development of more efficient, accurate computational techniques within the gas-phase ion measurement community. Specifically, we postulate whether experimental deviations from Langevin's polarization limit (K_pol) are suitable to estimate the relative contributions of hard-sphere collisions and long-range interactions within CCS values. Not surprisingly, different molecule classes exhibit different trends in the K/K_pol ratio when normalized for reduced mass, and the most common IMS calibrants (e.g., tune mix, polyalanine, tetraalkylammonium salts) follow different polarizability trends than many of the analytes probed in the literature. Succinctly, if gas-phase ion structure is largely invariant based upon the colliding neutral and newly developed experimental efforts can quantitatively capture ion polarizability, then modeling efforts describing a target analyte must be self-consistent as the collision neutral is changed in silico.

Evaluation of ion mobility in capillary electrophoresis coupled to mass spectrometry for the identification in metabolomics

Currently, feature annotation remains one of the main challenges in untargeted metabolomics. In this context, the information provided by high-resolution mass spectrometry (HRMS) in addition to accurate mass can improve the quality of metabolite annotation, and MS/MS fragmentation patterns are widely used. Accurate mass and a separation index, such as retention time or effective mobility (μ_eff ), in chromatographic and electrophoretic approaches, respectively, must be used for unequivocal metabolite identification. The possibility of measuring collision cross-section (CCS) values by using ion mobility (IM) is becoming increasingly popular in metabolomic studies thanks to the new generation of IM mass spectrometers. Based on their similar separation mechanisms involving electric field and the size of the compounds, the complementarity of ^DT CCS_N2 and μ_eff needs to be evaluated. In this study, a comparison of ^DT CCS_N2 and μ_eff was achieved in the context of feature identification ability in untargeted metabolomics by capillary zone electrophoresis (CZE) coupled with HRMS. This study confirms the high correlation of ^DT CCS_N2 with the mass of the studied metabolites as well as the orthogonality between accurate mass and μ_eff , making this combination particularly interesting for the identification of several endogenous metabolites. The use of IM-MS remains of great interest for facilitating the annotation of neutral metabolites present in the electroosmotic flow (EOF) that are poorly or not separated by CZE.

Isomers and Rotamers of DCM in Methanol and in Gas Phase Probed by Ion Mobility Mass Spectrometry in Combination with High Performance Liquid Chromatography

An integrated method of ion mobility mass spectrometry and high-performance liquid chromatography (HPLC) has been used to investigate the isomeric distribution of a popular fluorescent dye DCM (4-(dicyanomethylene)-2-methyl-6-(4-dimethylaminostyryl)-4H-pyran) in methanol solution. Chromatographic separation of DCM isomers in methanol has been performed by probing the molecular mass (DCMH⁺), and two distinctly separated peaks are observed at retention times 3.73 (peak-I) and 3.87 (peak-II) min, where the latter one appears nearly twice as intense as the former. However, peak-I appears much weaker compared to peak-II if the chromatogram is recorded by optical probing at the absorption maximum of this dye (467 nm). The ion mobility (IM) spectra of DCMH⁺ ions corresponding to each of the LC-separated factions show three common peaks A, B, and C, with collision cross-section (CCS) values of 174, 185, and 197 Ų, respectively, but their relative intensities in the two IM spectra appear in opposite sequences. The three IM peaks have been assigned by considering the theoretically calculated CCS values of 13 possible isomers of DCMH⁺ ions. The IM spectral features also reveal that isomeric interconversions occur during the ESI process. Electronic structure calculations have been used to optimize the geometries of the four isomers of solvated DCM and the corresponding protomeric structures of DCMH⁺. The isomerization pathways and associated energy barriers have also been calculated. The gas-phase protomers are found to follow a completely different sequence of stability as compared to the neutral isomers. The analysis reveals that peak-I corresponds to one of the cis isomers, whereas peak-II arises due to cumulative contributions of the other three isomers. The absorption spectrum of DCM in methanol is simulated from the computed spectral profiles of the isomers which indicates a distribution of trans1, trans2, cis1, and cis2 isomers as 33.5, 61.5, 2.0, and 3.0%, respectively. The fragmentation behavior of DCMH⁺ ions in a collision-induced dissociation experiment has been found to be isomer dependent.

Breaking Down Structural Diversity for Comprehensive Prediction of Ion-Neutral Collision Cross Sections

Identification of unknowns is a bottleneck for large-scale untargeted analyses like metabolomics or drug metabolite identification. Ion mobility-mass spectrometry (IM-MS) provides rapid two-dimensional separation of ions based on their mobility through a neutral buffer gas. The mobility of an ion is related to its collision cross section (CCS) with the buffer gas, a physical property that is determined by the size and shape of the ion. This structural dependency makes CCS a promising characteristic for compound identification, but this utility is limited by the availability of high-quality reference CCS values. CCS prediction using machine learning (ML) has recently shown promise in the field, but accurate and broadly applicable models are still lacking. Here we present a novel ML approach that employs a comprehensive collection of CCS values covering a wide range of chemical space. Using this diverse database, we identified the structural characteristics, represented by molecular quantum numbers (MQNs), that contribute to variance in CCS and assessed the performance of a variety of ML algorithms in predicting CCS. We found that by breaking down the chemical structural diversity using unsupervised clustering based on the MQNs, specific and accurate prediction models for each cluster can be trained, which showed superior performance than a single model trained with all data. Using this approach, we have robustly trained and characterized a CCS prediction model with high accuracy on diverse chemical structures. An all-in-one web interface (https://CCSbase.net) was built for querying the CCS database and accessing the predictive model to support unknown compound identifications.

A parallelized molecular collision cross section package with optimized accuracy and efficiency

A new parallelized calculation package predicts collision cross sections with high accuracy and efficiency.

Probing the structure of giant fullerenes by high resolution trapped ion mobility spectrometry

We present high-resolution trapped ion mobility spectrometry (TIMS) measurements for fullerene ions in molecular nitrogen.

Differential-Mobility Spectrometry of 1-Deoxysphingosine Isomers: New Insights into the Gas Phase Structures of Ionized Lipids

Separation and structural identification of lipids remain a major challenge for contemporary lipidomics. Regioisomeric lipids differing only in position(s) of unsaturation are often not differentiated by conventional liquid chromatography-mass spectrometry approaches leading to the incomplete, or sometimes incorrect, assignation of molecular structure. Here we describe an investigation of the gas phase separations by differential-mobility spectrometry (DMS) of a series of synthetic analogues of the recently described 1-deoxysphingosine. The dependence of the DMS behavior on the position of the carbon-carbon double bond within the ionized lipid is systematically explored and compared to trends from complementary investigations, including collision cross-sections measured by drift tube ion mobility, reaction efficiency with ozone, and molecular dynamics simulations. Consistent trends across these modes of interrogation point to the importance of direct, through-space interactions between the charge site and the carbon-carbon double bond. Differences in the geometry and energetics of this intramolecular interaction underpin DMS separations and influence reactivity trends between regioisomers. Importantly, the disruption and reformation of these intramolecular solvation interactions during DMS are proposed to be the causative factor in the observed separations of ionized lipids which are shown to have otherwise identical collision cross-sections. These findings provide key insights into the strengths and limitations of current ion-mobility technologies for lipid isomer separations and can thus guide a more systematic approach to improved analytical separations in lipidomics.

Identification of ortho-Substituted Benzoic Acid/Ester Derivatives via the Gas-Phase Neighboring Group Participation Effect in (+)-ESI High Resolution Mass Spectrometry

Benzoic acid/ester/amide derivatives are common moieties in pharmaceutical compounds and present a challenge in positional isomer identification by traditional tandem mass spectrometric analysis. A method is presented for exploiting the gas-phase neighboring group participation (NGP) effect to differentiate ortho-substituted benzoic acid/ester derivatives with high resolution mass spectrometry (HRMS¹). Significant water/alcohol loss (>30% abundance in MS¹ spectra) was observed for ortho-substituted nucleophilic groups; these fragment peaks are not observable for the corresponding para and meta-substituted analogs. Experiments were also extended to the analysis of two intermediates in the synthesis of suvorexant (Belsomra) with additional analysis conducted with nuclear magnetic resonance (NMR), density functional theory (DFT), and ion mobility spectrometry-mass spectrometry (IMS-MS) studies. Significant water/alcohol loss was also observed for 1-substituted 1, 2, 3-triazoles but not for the isomeric 2-substituted 1, 2, 3-triazole analogs. IMS-MS, NMR, and DFT studies were conducted to show that the preferred orientation of the 2-substituted triazole rotamer was away from the electrophilic center of the reaction, whereas the 1-subtituted triazole was oriented in close proximity to the center. Abundance of NGP product was determined to be a product of three factors: (1) proton affinity of the nucleophilic group; (2) steric impact of the nucleophile; and (3) proximity of the nucleophile to carboxylic acid/ester functional groups. Graphical Abstract ᅟ.

Structural characterization of small molecular ions by ion mobility mass spectrometry in nitrogen drift gas: improving the accuracy of trajectory method calculations

An accurate theoretical collision cross section calculation method in nitrogen was developed for reliable structural ion mobility mass spectrometry.

High resolution-ion mobility mass spectrometry as an additional powerful tool for structural characterization of mycotoxin metabolites

This work was designed as a proof of concept, to demonstrate the successful use of the comparison between theoretical and experimental collision cross section (CCS) values to support the identification of isomeric forms. To this purpose, thirteen mycotoxins were considered and analyzed using drift time ion mobility mass spectrometry. A good linear correlation (r² = 0.962) between theoretical and experimental CCS was found. The average ΔCCS was 3.2%, fully consistent with the acceptability threshold value commonly set at 5%. The agreement between theoretical and experimental CCS obtained for mycotoxin glucuronides suggested the potential of the CCS matching in supporting the annotation procedure.

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).

Effects of Charge State, Charge Distribution, and Structure on the Ion Mobility of Protein Ions in Helium Gas: Results from Trajectory Method Calculations

Collision cross section (Ω) values of gas-phase ions of proteins and protein complexes are used to probe the structures of the corresponding species in solution. Ions of many proteins exhibit increasing Ω-values with increasing charge state but most Ω-values calculated for protein ions have used simple collision models that do not explicitly account for charge. Here we use a combination of ion mobility mass spectrometry experiments with helium gas and trajectory method calculations to characterize the extents to which increases in experimental Ω-values with increasing charge state may be attributed to increased momentum transfer concomitant with enhanced long-range interactions between the protein ion and helium atoms. Ubiquitin and C-to-N terminally linked diubiquitin ions generated from different solution conditions exhibit more than a 2-fold increase in Ω with increasing charge state. For native and energy-relaxed models of the proteins and most methods for distributing charge, Ω-values calculated using the trajectory method increase by less than 1% over the range of charge states observed from typical solution conditions used for native mass spectrometry. However, the calculated Ω-values increase by 10% to 15% over the full range of charge states observed from all solution conditions. Therefore, contributions from enhanced ion-induced dipole interactions with increasing charge state are significant but without additional structural changes can account for only a fraction of the increase in Ω observed experimentally. On the basis of these results, we suggest guidelines for calculating Ω-values in the context of applications in biophysics and structural biology.

Benchmark Comparison for a Multi-Processing Ion Mobility Calculator in the Free Molecular Regime

A benchmark comparison between two ion mobility and collision cross-section (CCS) calculators, MOBCAL and IMoS, is presented here as a standard to test the efficiency and performance of both programs. Utilizing 47 organic ions, results are in excellent agreement between IMoS and MOBCAL in He and N_2, when both programs use identical input parameters. Due to a more efficiently written algorithm and to its parallelization, IMoS is able to calculate the same CCS (within 1%) with a speed around two orders of magnitude faster than its MOBCAL counterpart when seven cores are used. Due to the high computational cost of MOBCAL in N₂, reaching tens of thousands of seconds even for small ions, the comparison between IMoS and MOBCAL is stopped at 70 atoms. Large biomolecules (>10000 atoms) remain computationally expensive when IMoS is used in N₂ (even when employing 16 cores). Approximations such as diffuse trajectory methods (DHSS, TDHSS) with and without partial charges and projected area approximation corrections can be used to reduce the total computational time by several folds without hurting the accuracy of the solution. These latter methods can in principle be used with coarse-grained model structures and should yield acceptable CCS results. Graphical Abstract ᅟ.

Investigating Differences in Gas-Phase Conformations of 25-Hydroxyvitamin D3 Sodiated Epimers using Ion Mobility-Mass Spectrometry and Theoretical Modeling

Drift tube ion mobility coupled with mass spectrometry was used to investigate the gas-phase structure of 25-hydroxyvitamin D3 (25OHD3) and D2 (25OHD2) epimers, and to evaluate its potential in rapid separation of these compounds. Experimental results revealed two distinct drift species for the 25OHD3 sodiated monomer, whereas only one of these conformations was observed for its epimer (epi25OHD3). The unique species allowed 25OHD3 to be readily distinguished, and the same pattern was observed for 25OHD2 epimers. Theoretical modeling of 25OHD3 epimers identified energetically stable gas-phase structures, indicating that both compounds may adopt a compact "closed" conformation, but that 25OHD3 may also adopt a slightly less energetically favorable "open" conformation that is not accessible to its epimer. Calculated theoretical collision cross-sections for these structures agreed with experimental results to <2%. Experimentation indicated that additional energy in the ESI source (i.e., increased temperature, spray voltage) affected the ratio of 25OHD3 conformations, with the less energetically favorable "open" conformation increasing in relative intensity. Finally, LC-IM-MS results yielded linear quantitation of 25OHD3, in the presence of the epimer interference, at biologically relevant concentrations. This study demonstrates that ion mobility can be used in tandem with theoretical modeling to determine structural differences that contribute to drift separation. These separation capabilities provide potential for rapid (<60 ms) identification of 25OHD3 and 25OHD2 in mixtures with their epimers. Graphical Abstract ᅟ.

Molecular Structures and Momentum Transfer Cross Sections: The Influence of the Analyte Charge Distribution

Structure elucidation by ion mobility spectrometry-mass spectrometry methods is based on the comparison of an experimentally measured momentum transfer cross-section to cross-sections calculated for model structures. Thus, it is imperative that the calculated cross-section must be accurate. However, it is not fully understood how important it is to accurately model the charge distribution of an analyte ion when calculating momentum transfer cross-sections. Here, we calculate and compare momentum transfer cross-sections for carbon clusters that differ in mass, charge state, and mode of charge distribution, and vary temperature and polarizability of the buffer gas. Our data indicate that the detailed distribution of the ion charge density is intimately linked to the contribution of glancing collisions to the momentum transfer cross-section. The data suggest that analyte ions with molecular mass ~3 kDa or momentum transfer cross-section 400-500 Ų would be significantly influenced by the charge distribution in nitrogen buffer gas. Our data further suggest that accurate structure elucidation on the basis of IMS-MS data measured in nitrogen buffer gas must account for the molecular charge distribution even for systems as large as C₉₆₀ (~12 kDa) when localized charges are present and/or measurements are conducted under cryogenic temperatures. Finally, our data underscore that accurate structure elucidation is unlikely if ion mobility data recorded in one buffer gas is converted into other buffer gases when electronic properties of the buffer gases differ. Graphical Abstract ᅟ.

Experimental and Theoretical Investigation of Sodiated Multimers of Steroid Epimers with Ion Mobility-Mass Spectrometry

Ion mobility-mass spectrometry (IM-MS) has recently seen increased use in the analysis of small molecules, especially in the field of metabolomics, for increased breadth of information and improved separation of isomers. In this study, steroid epimers androsterone and trans-androsterone were analyzed with IM-MS to investigate differences in their relative mobilities. Although sodiated monomers exhibited very similar collision cross-sections (CCS), baseline separation was observed for the sodiated dimer species (R_S = 1.81), with measured CCS of 242.6 and 256.3 Ų, respectively. Theoretical modeling was performed to determine the most energetically stable structures of solution-phase and gas-phase monomer and dimer structures. It was revealed that these epimers differ in their preferred dimer binding mode in solution phase: androsterone adopts a R=O - Na⁺ - OH-R' configuration, whereas trans-androsterone adopts a R=O - Na⁺ - O=R' configuration. This difference contributes to a significant structural variation, and subsequent CCS calculations based on these structures relaxed in the gas phase were in agreement with experimentally measured values (ΔCCS ~ 5%). Additionally, these calculations accurately predicted the relative difference in mobility between the epimers. This study illustrates the power of combining experimental and theoretical results to better elucidate gas-phase structures. Graphical Abstract ᅟ.

Collision cross sections and ion structures: development of a general calculation method via high-quality ion mobility measurements and theoretical modeling

Theoretical collision cross section calculations revisited for reliable ion structural studies.

Gas molecule scattering & ion mobility measurements for organic macro-ions in He versus N2 environments

A pending issue in linking ion mobility measurements to ion structures is that the collisional cross section (CCS, the measured structural parameter in ion mobility spectrometry) of an ion is strongly dependent upon the manner in which gas molecules effectively impinge on and are reemitted from ion surfaces (when modeling ions as fixed structures). To directly examine the gas molecule impingement and reemission processes and their influence, we measured the CCSs of positively charged ions of room temperature ionic liquids 1-ethyl-3-methylimidazolium dicyanamide (EMIM-N(CN)2) and 1-ethyl-3-methylimidazolium tetrafluoroborate (EMIM-BF4) in N2 using a differential mobility analyzer-mass spectrometer (DMA-MS) and in He using a drift tube mobility spectrometer-mass spectrometer (DT-MS). Cluster ions, generated via electrosprays, took the form (AB)N(A)z, spanning up to z = 20 and with masses greater than 100 kDa. As confirmed by molecular dynamics simulations, at the measurement temperature (∼300 K), such cluster ions took on globular conformations in the gas phase. Based upon their attained charge levels, in neither He nor N2 did the ion-induced dipole potential significantly influence gas molecule-ion collisions. Therefore, differences in the CCSs measured for ions in the two different gases could be primarily attributed to differences in gas molecule behavior upon collision with ions. Overwhelmingly, by comparison of predicted CCSs with selected input impingement-reemission laws to measurements, we find that in N2, gas molecules collide with ions diffusely--they are reemitted at random angles relative to the gas molecule incoming angle--and inelastically. Meanwhile, in He, gas molecules collide specularly and elastically and are emitted from ion surfaces at determined angles. The results can be rationalized on the basis of the momentum transferred per collision; in the case of He, individual gas molecule collisions minimally perturb the atoms within a cluster ion (internal motion), while in the case of N2, individual gas molecules have sufficiently large momentum to alter the internal motion in organic ions.

Combining density functional theory (DFT) and collision cross-section (CCS) calculations to analyze the gas-phase behaviour of small molecules and their protonation site isomers

Computational methods are employed to study the protomers in ESI-IM-MS.

Linking molecular models with ion mobility experiments. Illustration with a rigid nucleic acid structure

Ion mobility spectrometry experiments allow the mass spectrometrist to determine an ion's rotationally averaged collision cross section Ω(EXP). Molecular modelling is used to visualize what ion three-dimensional structure(s) is(are) compatible with the experiment. The collision cross sections of candidate molecular models have to be calculated, and the resulting Ω(CALC) are compared with the experimental data. Researchers who want to apply this strategy to a new type of molecule face many questions: (1) What experimental error is associated with Ω(EXP) determination, and how to estimate it (in particular when using a calibration for traveling wave ion guides)? (2) How to generate plausible 3D models in the gas phase? (3) Different collision cross section calculation models exist, which have been developed for other analytes than mine. Which one(s) can I apply to my systems? To apply ion mobility spectrometry to nucleic acid structural characterization, we explored each of these questions using a rigid structure which we know is preserved in the gas phase: the tetramolecular G-quadruplex [dTGGGGT]4, and we will present these detailed investigation in this tutorial.

A local collision probability approximation for predicting momentum transfer cross sections

The local collision probability approximation (LCPA) method is introduced to compute molecular momentum transfer cross sections for comparison to ion mobility experiments.

Evaluation of drift gas selection in complex sample analyses using a high performance drift tube ion mobility-QTOF mass spectrometer

Drift gas selection in an ion mobility-QTOF can affect the separation of targeted isomeric compounds and collision cross section values.

Characterization of polylactides with different stereoregularity using electrospray ionization ion mobility mass spectrometry

We investigated the effect of stereoregularity on the gas-phase conformations of linear and cyclic polylactides (PLA) using electrospray ionization ion mobility mass spectrometry (ESI-IM-MS) combined with molecular dynamics simulations. IM-MS analysis of PLA ions shows intriguing difference between the collision cross section (ΩD) value of poly-L-lactide (PLLA) and poly-LD-lactide (PLDLA) ions with respect to their chain architecture and stereoregularity. In the singly sodiated linear PLA (l-PLA∙Na(+)) case, both l-PLLA and l-PLDLA up to 11mer have very similar ΩD values, but the ΩD values of l-PLLA are greater than that of l-PLDLA ions for larger ions. In the case of cyclic PLA (c-PLA), c-PLLA∙Na(+) is more compact than c-PLDLA∙Na(+) for short PLA ions. However, c-PLLA exhibits larger ΩD value than c-PLDLA for PLA ions longer than 13mer. The origin of difference in the ΩD values was investigated using theoretical investigation of PLAs in the gas phase. The gas-phase conformation of PLA ions is influenced by Na(+)-oxygen coordination and the weak intramolecular hydrogen bond interaction, which are more effectively formed in more flexible chains. Therefore, the less flexible PLLA has a larger ΩD value than PLDLA. However, for short c-PLA, concomitant maximization of both Na(+)-oxygen coordination and hydrogen bond interaction is difficult due to the constricted chain freedom, which makes the ΩD value of PLAs in this range show a different trend compared with other PLA ions. Our study facilitates the understanding of correlation between stereoregularity of PLAs and their structure, providing potential utility of IM-MS to characterize stereoisomers of polymers.

Ion-molecule clustering in differential mobility spectrometry: lessons learned from tetraalkylammonium cations and their isomers

Differential mobility spectrometry (DMS) can distinguish ions based upon the differences in their high- and low-field ion mobilities as they experience the asymmetric waveform applied to the DMS cell. These mobilities are known to be influenced by the ions' structure, m/z, and charge distribution (i.e., resonance structures) within the ions themselves, as well as by the gas-phase environment of the DMS cell. While these associations have been developed over time through empirical observations, the exact role of ion structures or their interactions with clustering molecules remains generally unknown. In this study, that relationship is explored by observing the DMS behaviors of a series of tetraalkylammonium ions as a function of their structures and the gas-phase environment of the DMS cell. To support the DMS experiments, the basin-hopping search strategy was employed to identify candidate cluster structures for density functional theory treatment. More than a million cluster structures distributed across 72 different ion-molecule cluster systems were sampled to determine global minimum structures and cluster binding energies. This joint computational and experimental approach suggests that cluster geometry, in particular ion-molecule intermolecular separation, plays a critical role in DMS.

Understanding gas phase modifier interactions in rapid analysis by differential mobility-tandem mass spectrometry

A systematic study involving the use and optimization of gas-phase modifiers in quantitative differential mobility-mass spectrometry (DMS-MS) analysis is presented using nucleoside-adduct biomarkers of DNA damage as an important reference point for analysis in complex matrices. Commonly used polar protic and polar aprotic modifiers have been screened for use against two deoxyguanosine adducts of DNA: N-(deoxyguanosin-8-yl)-4-aminobiphenyl (dG-C8-4-ABP) and N-(deoxyguanosin-8-y1)-2-amino-l-methyl-6-phenylimidazo[4,5-b]pyridine (dG-C8-PhIP). Particular attention was paid to compensation voltage (CoV) shifts, peak shapes, and product ion signal intensities while optimizing the DMS-MS conditions. The optimized parameters were then applied to rapid quantitation of the DNA adducts in calf thymus DNA. After a protein precipitation step, adduct levels corresponding to less than one modification in 10(6) normal DNA bases were detected using the DMS-MS platform. Based on DMS fundamentals and ab initio thermochemical results, we interpret the complexity of DMS modifier responses in terms of thermal activation and the development of solvent shells. At very high bulk gas temperature, modifier dipole moment may be the most important factor in cluster formation and cluster geometry, but at lower temperatures, multi-neutral clusters are important and less predictable. This work provides a useful protocol for targeted DNA adduct quantitation and a basis for future work on DMS modifier effects.

Baseline resolution of isomers by traveling wave ion mobility mass spectrometry: investigating the effects of polarizable drift gases and ionic charge distribution

We have studied the behavior of isomers and analogues by traveling wave ion mobility mass spectrometry (TWIM-MS) using drift-gases with varying masses and polarizabilities. Despite the reduced length of the cell (18 cm), a pair of constitutional isomers, N-butylaniline and para-butylaniline, with theoretical collision cross-section values in helium (ΩHe ) differing by as little as 1.2 Å(2) (1.5%) but possessing contrasting charge distribution, showed baseline peak-to-peak resolution (Rp-p ) for their protonated molecules, using carbon dioxide (CO2), nitrous oxide (N2O) and ethene (C2H4 ) as the TWIM drift-gas. Near baseline Rp-p was also obtained in CO2 for a group of protonated haloanilines (para-chloroaniline, para-bromoaniline and para-iodoaniline) which display contrasting masses and theoretical ΩHe , which differ by as much as 15.7 Å(2) (19.5%) but similar charge distributions. The deprotonated isomeric pair of trans-oleic acid and cis-oleic acid possessing nearly identical theoretical ΩHe and ΩN2 as well as similar charge distributions, remained unresolved. Interestingly, an inversion of drift-times were observed for the 1,3-dialkylimidazolium ions when comparing He, N2 and N2O. Using density functional theory as a means of examining the ions electronic structure, and He and N2-based trajectory method algorithm, we discuss the effect of the long-range charge induced dipole attractive and short-range Van der Waals forces involved in the TWIM separation in drift-gases of differing polarizabilities. We therefore propose that examining the electronic structure of the ions under investigation may potentially indicate whether the use of more polarizable drift-gases could improve separation and the overall success of TWIM-MS analysis.