Characterization of Protonated Model Disaccharides from Tandem Mass Spectrometry and Chemical Dynamics Simulations

Quantum Chemistry-based Molecular Dynamics Simulations as a Tool for the Assignment of ESI-MS/MS Spectra of Drug Molecules

In organic mass spectrometry, fragment ions provide important information on the analyte as a central part of its structure elucidation. With increasing molecular size and possible protonation sites, the potential energy surface (PES) of the analyte can become very complex, which results in a large number of possible fragmentation patterns. Quantum chemical (QC) calculations can help here, enabling the fast calculation of the PES and thus enhancing the mass spectrometry-based structure elucidation processes. In this work, the previously unknown fragmentation pathways of the two drug molecules Nateglinide (45 atoms) and Zopiclone (51 atoms) were investigated using a combination of generic formalisms and calculations conducted with the Quantum Chemical Mass Spectrometry (QCxMS) program. The computations of the de novo fragment spectra were conducted with the semi-empirical GFNn-xTB (n=1, 2) methods and compared against Orbitrap measured electrospray ionization (ESI) spectra in positive ion mode. It was found that the unbiased QC calculations are particularly suitable to predict non-evident fragment ion structures, sometimes contrasting the accepted generic formulation of fragment ion structures from electron migration rules, where the "true" ion fragment structures are approximated. For the first time, all fragment and intermediate structures of these large-sized molecules could be elucidated completely and routinely using this merger of methods, finding new undocumented mechanisms, that are not considered in common rules published so far. Given the importance of ESI for medicinal chemistry, pharmacokinetics, and metabolomics, this approach can significantly enhance the mass spectrometry-based structure elucidation processes and contribute to the understanding of previously unknown fragmentation pathways.

Collision-induced dissociation of protonated uracil water clusters probed by molecular dynamics simulations

Collision-induced dissociation experiments of hydrated molecular species can provide a wealth of important information. However, they often need a theoretical support to extract chemical information. In the present article, in order to provide a detailed description of recent experimental measurements [Braud et al., J. Chem. Phys., 2019, 150, 014303], collision simulations between low-energy protonated uracil water clusters (H₂O)_1-7,11,12UH⁺ and an Ar atom were performed using a quantum mechanics/molecular mechanics formalism based on the self-consistent-charge density-functional based tight-binding method. The theoretical proportion of formed neutral vs. protonated uracil containing clusters, total fragmentation cross sections as well as the mass spectra of charged fragments are consistent with the experimental data which highlights the accuracy of the present simulations. They allow to probe which fragments are formed on the short time scale and rationalize the location of the excess proton on these fragments. We demonstrate that this latter property is highly influenced by the nature of the aggregate undergoing the collision. Analyses of the time evolution of the fragments populations and of their relative abundances demonstrate that, up to 7 water molecules, a direct dissociation mechanism occurs after collision whereas for 11 and 12 water molecules a statistical mechanism is more likely to participate. Although scarce in the literature, the present simulations appear as a useful tool to complement collision-induced dissociation experiments of hydrated molecular species.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Size Dependent Fragmentation Chemistry of Short Doubly Protonated Tryptic Peptides

Tandem mass spectrometry of electrospray ionized multiply charged peptide ions is commonly used to identify the sequence of peptide(s) and infer the identity of source protein(s). Doubly protonated peptide ions are consistently the most efficiently sequenced ions following collision-induced dissociation of peptides generated by tryptic digestion. While the broad characteristics of longer (N ≥ 8 residue) doubly protonated peptides have been investigated, there is comparatively little data on shorter systems where charge repulsion should exhibit the greatest influence on the dissociation chemistry. To address this gap and further understand the chemistry underlying collisional-dissociation of doubly charged tryptic peptides, two series of analytes ([G_xR+2H]²⁺ and [A_xR+2H]²⁺, x = 2-5) were investigated experimentally and with theory. We find distinct differences in the preference of bond cleavage sites for these peptides as a function of size and to a lesser extent composition. Density functional calculations at two levels of theory predict that the threshold relative energies required for bond cleavages at the same site for peptides of different size are quite similar (for example, b₂-y_N-2). In isolation, this finding is inconsistent with experiment. However, the predicted extent of entropy change of these reactions is size dependent. Subsequent RRKM rate constant calculations provide a far clearer picture of the kinetics of the competing bond cleavage reactions enabling rationalization of experimental findings. The M06-2X data were substantially more consistent with experiment than were the B3LYP data.

Linkage Memory in Underivatized Protonated Carbohydrates

Carbohydrates are among the most complex class of biomolecules, and even subtle variations in their structures are attributed to diverse biological functions. Mass spectrometry has been essential for large scale glycomics and glycoproteomics studies, but the gas-phase structures and sometimes anomalous fragmentation properties of carbohydrates present long-standing challenges. Here we investigate the gas-phase properties of a panel of isomeric protonated disaccharides differing in their linkage configurations. Multiple conformations were evident for most of the structures based on their fragment ion abundances by tandem mass spectrometry, their ion mobilities in several gases, and their deuterium uptake kinetics by gas-phase hydrogen-deuterium exchange. Most notably, we find that the properties of the Y-ion fragments are characteristically influenced by the precursor carbohydrate's linkage configuration. This study reveals how protonated carbohydrate fragment ions can retain "linkage memory" that provides structural insight into their intact precursor.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Fragmentation Spectra Prediction and DNA Adducts Structural Determination

In this work, chemical dynamics simulations were optimized and used to predict fragmentation mass spectra for DNA adduct structural determination. O⁶-methylguanine (O⁶-Me-G) was used as a simple model adduct to calculate theoretical spectra for comparison with measured high-resolution fragmentation data. An automatic protocol was established to consider the different tautomers accessible at a given energy and obtain final theoretical spectra by insertion of an initial tautomer. In the work reported here, the most stable tautomer was chosen as the initial structure, but in general, any structure could be considered. Allowing for the formation of the various possible tautomers during simulation calculations was found to be important to getting a more complete fragmentation spectrum. The calculated theoretical results reproduce the experimental peaks such that it was possible to determine reaction pathways and product structures. The calculated tautomerization network was crucial to correctly identifying all the observed ion peaks, showing that a mobile proton model holds not only for peptide fragmentation but also for nucleobases. Finally, first principles results were compared to simple machine learning fragmentation models.

ESI-QTof-MS characterization of hirsutinolide and glaucolide sesquiterpene lactones: Fragmentation mechanisms and differentiation based on Na+ /H+ adducts interactions in complex mixture

Sesquiterpene lactones (SL) have been reported with various biological effects. Among the described SL skeletons, hirsutinolide and glaucolide have not been extensively studied by mass spectrometry (MS), especially how to distinguish them in organic matrices. Thus, this paper reports (1) a strategy of their differentiation based on MS behavior during the ionization and (2) a proposal of the fragmentation pattern for both SL-subtypes. ESI(+)-HRMS data of four isolated SL (hirsutinolides 1 and 3; glaucolides 2 and 4) were recorded by direct and UPLC water-sample combined injections. These analyses revealed that hirsutinolides and glaucolides formed [M+Na]⁺ ion during the operation of the direct MS injection, and ([M+Na]⁺ and [M+H-H₂ O]⁺ ) and [M+H]⁺ ions were respectively observed for hirsutinolides and glaucolides during the operation of combined UPLC water and sample MS injection. Computational simulations showed that the complex hirsutinolide (1)-Na⁺ formed with a lower preparation energy compared with the complex glaucolide (2)-Na⁺ . However, despite their different behavior during the ionization process, ESI(+)-HRMS/MS analyses of 1-4 gave similar fragmentation patterns at m/z 277, 259, 241, and 231 that can be used as diagnostic ions for both skeletons. Moreover, the differentiation strategy based on the nature of the complex SL-adducts and their MS/MS fragmentation pattern were successfully applied for the chemical characterization of the extract from Vernonanthura tweedieana using UPLC-ESI-HRMS/MS. Among the characterized metabolites, SL with hirsutinolide and glaucolide skeletons showed the aforementioned diagnostic fragments and an ionization behavior that was similar to those observed during the water-sample combined injection.

Advancing Solutions to the Carbohydrate Sequencing Challenge

Carbohydrates possess a variety of distinct features with stereochemistry playing a particularly important role in distinguishing their structure and function. Monosaccharide building blocks are defined by a high density of chiral centers. Additionally, the anomericity and regiochemistry of the glycosidic linkages carry important biological information. Any carbohydrate-sequencing method needs to be precise in determining all aspects of this stereodiversity. Recently, several advances have been made in developing fast and precise analytical techniques that have the potential to address the stereochemical complexity of carbohydrates. This perspective seeks to provide an overview of some of these emerging techniques, focusing on those that are based on NMR and MS-hybridized technologies including ion mobility spectrometry and IR spectroscopy.

A Trajectory-Based Method to Explore Reaction Mechanisms

The tsscds method, recently developed in our group, discovers chemical reaction mechanisms with minimal human intervention. It employs accelerated molecular dynamics, spectral graph theory, statistical rate theory and stochastic simulations to uncover chemical reaction paths and to solve the kinetics at the experimental conditions. In the present review, its application to solve mechanistic/kinetics problems in different research areas will be presented. Examples will be given of reactions involved in photodissociation dynamics, mass spectrometry, combustion chemistry and organometallic catalysis. Some planned improvements will also be described.

Dissection of Fragmentation Pathways in Protonated N-Acetylhexosamines

Structural characterization of carbohydrates by mass spectrometry necessitates a detailed understanding of their gas phase behavior, particularly for protonated carbohydrates that can undergo complex structural rearrangements during fragmentation. Here we utilize tandem mass spectrometry, isotopic labeling, gas-phase hydrogen/deuterium exchange, and ion mobility measurements to characterize structures of the various product ions of protonated N-acetylhexosamines. Following the facile loss of the reducing end hydroxyl group, we identify two primary fragmentation pathways. Detailed mapping of each step in the fragmentation pathway provides new insight into the mechanisms that drive collision-induced dissociation of protonated carbohydrates. Several of the smaller fragment ions are mixtures of structural isomers, and the relative distributions of these structures reveals information about the stereochemistry of the precursor molecule.

Direct Experimental Characterization of Glycosyl Cations by Infrared Ion Spectroscopy

Glycosyl cations are crucial intermediates formed during enzymatic and chemical glycosylation. The intrinsic high reactivity and short lifetime of these reaction intermediates make them very challenging to characterize using spectroscopic techniques. Herein, we report the use of collision induced dissociation tandem mass spectrometry to generate glycosyl cations in the gas phase followed by infrared ion spectroscopy using the FELIX infrared free electron laser. The experimentally observed IR spectra were compared to DFT calculated spectra enabling the detailed structural elucidation of elusive glycosyl oxocarbenium and dioxolenium ions.

Chemical dynamics simulations of CID of peptide ions: comparisons between TIK(H+)2 and TLK(H+)2 fragmentation dynamics, and with thermal simulations

Gas phase unimolecular fragmentation of the two model doubly protonated tripeptides threonine-isoleucine-lysine (TIK) and threonine-leucine-lysine (TLK) is studied using chemical dynamics simulations. Attention is focused on different aspects of collision induced dissociation (CID): fragmentation pathways, energy transfer, theoretical mass spectra, fragmentation mechanisms, and the possibility of distinguishing isoleucine (I) and leucine (L). Furthermore, discussion is given regarding the differences between single collision CID activation, which results from a localized impact between the ions and a colliding molecule N₂, and previous thermal activation simulation results; Z. Homayoon, S. Pratihar, E. Dratz, R. Snider, R. Spezia, G. L. Barnes, V. Macaluso, A. Martin-Somer and W. L. Hase, J. Phys. Chem. A, 2016, 120, 8211-8227. Upon thermal activation unimolecular fragmentation is statistical and in accord with RRKM unimolecular rate theory. Simulations show that in collisional activation some non-statistical fragmentation occurs, including shattering, which is not present when the ions dissociate statistically. Products formed by non-statistical shattering mechanisms may be related to characteristic mass spectrometry peaks which distinguish the two isomers I and L.

Collision-induced dissociation of sodiated glucose, galactose, and mannose, and the identification of anomeric configurations

Different dehydration barrier heights of cis and trans configurations between O1 and O2 provide a simple and fast anomeric configuration identification.