Direct Dynamics Simulations of Collision- and Surface-Induced Dissociation of N-Protonated Glycine. Shattering Fragmentation

Modeling collisional kinetic energy damping, heating, and cooling of ions in mass spectrometers: a tutorial perspective

Many powerful methods in mass spectrometry rely on activation of ions by high-energy collisions with gas particles. For example, multiple Collision Induced Dissociation (CID) has been used for many years to determine structural information for ions ranging from small organics to large, native-like protein complexes. More recently, Collision Induced Unfolding (CIU) has proved to be a very powerful method for understanding high-order protein structure and detecting differences between similar proteins. Quantifying the thermochemistry underlying dissociation/unfolding in these experiments can be quite challenging without reliable models of ion heating and cooling. Established physical models of CID are valuable in predicting ion heating but do not explicitly include mechanisms for cooling, which may play a large part in CID/CIU in modern instruments. Ab initio and Molecular Dynamics methods are extremely computationally expensive for modeling CID/CIU of large analytes such as biomolecular ions. In this tutorial perspective, limiting behaviors of ion kinetic energy damping, heating, and cooling set by "extreme" cases are explored, and an Improved Impulsive Collision Theory and associated software ("Ion Simulations of the Physics of Activation", IonSPA) are introduced that can model all of these for partially inelastic collisions. Finally, examples of modeled collisional activation of native-like protein ions under realistic experimental conditions are discussed, with an outlook toward the use of IonSPA in accessing the thermochemical information hidden in CID breakdown curves and CIU fingerprints.

Surface-induced Dissociation Mass Spectrometry as a Structural Biology Tool

Native mass spectrometry (nMS) is evolving into a workhorse for structural biology. The plethora of online and offline preparation, separation, and purification methods as well as numerous ionization techniques combined with powerful new hybrid ion mobility and mass spectrometry systems has illustrated the great potential of nMS for structural biology. Fundamental to the progression of nMS has been the development of novel activation methods for dissociating proteins and protein complexes to deduce primary, secondary, tertiary, and quaternary structure through the combined use of multiple MS/MS technologies. This review highlights the key features and advantages of surface collisions (surface-induced dissociation, SID) for probing the connectivity of subunits within protein and nucleoprotein complexes and, in particular, for solving protein structure in conjunction with complementary techniques such as cryo-EM and computational modeling. Several case studies highlight the significant role SID, and more generally nMS, will play in structural elucidation of biological assemblies in the future as the technology becomes more widely adopted. Cases are presented where SID agrees with solved crystal or cryoEM structures or provides connectivity maps that are otherwise inaccessible by "gold standard" structural biology techniques.

How Symmetry Influences the Dissociation of Protonated Cyclic Peptides

Protonated cyclic dipeptides undergo collision-induced dissociation, and this reaction mechanism strongly depends on the symmetry and the nature of the residues. We review the main dissociation mechanism for a series of cyclic dipeptides, obtained through chemical dynamics simulations. The systems range from the symmetrical cyclo-(glycyl-glycyl), with two possible symmetrical protonation sites located on the peptide ring, to cyclo-(tyrosyl-prolyl), where the symmetry of protonation sites on the peptide ring is broken by the dissimilar nature of the different residues. Finally, cyclo-(phenylalanyl-histidyl) shows a completely asymmetric situation, with the proton located on one of the dipeptide side chains, which explains the peculiar fragmentation mechanism induced by shuttling the proton, whose efficiency is strongly dependent on the relative chirality of the residues.

PyFragMS-A Web Tool for the Investigation of the Collision-Induced Fragmentation Pathways

Dissociation induced by the accumulation of internal energy via collisions of ions with neutral molecules is one of the most important fragmentation techniques in mass spectrometry (MS), and the identification of small singly charged molecules is based mainly on the consideration of the fragmentation spectrum. Many research studies have been dedicated to the creation of databases of experimentally measured tandem mass spectrometry (MS/MS) spectra (such as MzCloud, Metlin, etc.) and developing software for predicting MS/MS fragments in silico from the molecular structure (such as MetFrag, CFM-ID, CSI:FingerID, etc.). However, the fragmentation mechanisms and pathways are still not fully understood. One of the limiting obstacles is that protomers (positive ions protonated at different sites) produce different fragmentation spectra, and these spectra overlap in the case of the presence of different protomers. Here, we are proposing to use a combination of two powerful approaches: computing fragmentation trees that carry information of all consecutive fragmentations and consideration of the MS/MS data of isotopically labeled compounds. We have created PyFragMS-a web tool consisting of a database of annotated MS/MS spectra of isotopically labeled molecules (after H/D and/or ¹⁶O/¹⁸O exchange) and a collection of instruments for computing fragmentation trees for an arbitrary molecule. Using PyFragMS, we investigated how the site of protonation influences the fragmentation pathway for small molecules. Also, PyFragMS offers capabilities for performing database search when MS/MS data of the isotopically labeled compounds are taken into account.

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.

Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories

In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.

From QCEIMS to QCxMS: A Tool to Routinely Calculate CID Mass Spectra Using Molecular Dynamics

Mass spectrometry (MS) is a powerful tool in chemical research and substance identification. For the computational modeling of electron ionization MS, we have developed the quantum-chemical electron ionization mass spectra (QCEIMS) program. Here, we present an extension of QCEIMS to calculate collision-induced dissociation (CID) spectra. The more general applicability is accounted for by the new name QCxMS, where "x" refers to EI or CID. To this end, fragmentation and rearrangement reactions are computed "on-the-fly" in Born-Oppenheimer molecular dynamics (MD) simulations with the semiempirical GFN2-xTB Hamiltonian, which provides an efficient quantum mechanical description of all elements up to Z = 86 (Rn). Through the explicit modeling of multicollision processes between precursor ions and neutral gas atoms as well as temperature-induced decomposition reactions, QCxMS provides detailed insight into the collision kinetics and fragmentation pathways. In combination with the CREST program to determine the preferential protonation sites, QCxMS becomes the first standalone MD-based program that can predict mass spectra based solely on molecular structures as input. We demonstrate this for six organic molecules with masses ranging from 159 to 296 Da, for which QCxMS yields CID spectra in reasonable agreement with experiments.

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.

Tandem surface-induced dissociation of protein complexes on an ultrahigh resolution platform

We describe instrumentation for conducting tandem surface-induced dissociation (tSID) of native protein complexes on an ultrahigh resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. The two stages of SID are accomplished with split lenses replacing the entrance lenses of the quadrupole mass filter (stage 1, referred to herein as SID-Q) and the collision cell (stage 2, Q-SID). After SID-Q, the scattered projectile ions and subcomplexes formed in transit traverse the 20 mm pre-filter prior to the mass-selecting quadrupole, providing preliminary insights into the SID fragmentation kinetics of noncovalent protein complexes. The isolated SID fragments (subcomplexes) are then fragmented by SID in the collision cell entrance lens (Q-SID), generating subcomplexes of subcomplexes. We show that the ultrahigh resolution of the FT-ICR can be used for deconvolving species overlapping in m/z, which are particularly prominent in tandem SID spectra due to the combination of symmetric charge partitioning and narrow product ion charge state distributions. Various protein complex topologies are explored, including homotetramers, homopentamers, a homohexamer, and a heterohexamer.

Direct Dynamics Simulations of the Thermal Fragmentation of a Protonated Peptide Containing Arginine

Arginine has significant effects on fragmentation patterns of the protonated peptide due to its high basicity guanidine tail. In this article, thermal dissociation of the singly protonated glycine-arginine dipeptide (GR-H⁺) was investigated by performing direct dynamics simulations at different vibrational temperatures of 2000-3500 K. Fourteen principal fragmentation mechanisms containing side-chain and backbone fragmentation were found and discussed in detail. The mechanism involving partial or complete loss of a guanidino group dominates side-chain fragmentation, while backbone fragmentation mainly involves the three cleavage sites of a1-x1⁺, a2⁺-x0, and b1-y1⁺. Fragmentation patterns for primary dissociation have been compared with experimental results, and the peak that was not identified by the experiment has been assigned by our simulation. Kinetic parameters for GR-H⁺ unimolecular dissociation may be determined by direct dynamics simulations, which are helpful in exploring the complex biomolecules.

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.

Theoretical investigation of the dissociation chemistry of formyl halides in the gas phase

Halogen substituted analogues of formaldehyde, HXCO (X = F, Cl, Br, and I), play a crucial role in the degradation of stratospheric ozone.

Collisional dynamics simulations revealing fragmentation properties of Zn(ii)-bound poly-peptide

Collisional simulations show how peptide fragmentation is modified by the presence of Zn(ii).

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.

Experimental and theoretical investigation of overall energy deposition in surface-induced unfolding of protein ions

Recent advances in native mass spectrometry have enabled its use to probe the structure of and interactions within biomolecular complexes. Surface-induced dissociation, in which inter- and intramolecular interactions are disrupted following an energetic ion-surface collision, is a method that can directly interrogate the topology of protein complexes. However, a quantitative relationship between the ion kinetic energy at the moment of surface collision and the internal energy deposited into the ion has not yet been established for proteins. The factors affecting energy deposition in surface-induced unfolding (SIU) of protein monomers were investigated and a calibration relating laboratory-frame kinetic energy to internal energy developed. Protein monomers were unfolded by SIU and by collision-induced unfolding (CIU). CIU and SIU cause proteins to undergo the same unfolding transitions at different values of laboratory-frame kinetic energy. There is a strong correlation between the SIU and CIU energies, demonstrating that SIU, like CIU, can largely be understood as a thermal process. The change in internal energy in CIU was modeled using a Monte Carlo approach and theory. Computed values of the overall efficiency were found to be approximately 25% and used to rescale the CIU energy axis and relate nominal SIU energies to internal energy. The energy deposition efficiency in SIU increases with mass and kinetic energy from a low of ∼20% to a high of ∼68%, indicating that the effective mass of the surface increases along with the mass of the ion. The effect of ion structure on energy deposition was probed using multiple stages of ion activation. Energy deposition in SIU strongly depends on structure, decreasing as the protein is elongated, due to decreased effective protein-surface collisional cross section and increased transfer to rotational modes.

Resolution and Assignment of Differential Ion Mobility Spectra of Sarcosine and Isomers

Due to their central role in biochemical processes, fast separation and identification of amino acids (AA) is of importance in many areas of the biomedical field including the diagnosis and monitoring of inborn errors of metabolism and biomarker discovery. Due to the large number of AA together with their isomers and isobars, common methods of AA analysis are tedious and time-consuming because they include a chromatographic separation step requiring pre- or post-column derivatization. Here, we propose a rapid method of separation and identification of sarcosine, a biomarker candidate of prostate cancer, from isomers using differential ion mobility spectrometry (DIMS) interfaced with a tandem mass spectrometer (MS/MS) instrument. Baseline separation of protonated sarcosine from α- and β-alanine isomers can be easily achieved. Identification of DIMS peak is performed using an isomer-specific activation mode where DIMS- and mass-selected ions are irradiated at selected wavenumbers allowing for the specific fragmentation via an infrared multiple photon dissociation (IRMPD) process. Two orthogonal methods to MS/MS are thus added, where the MS/MS(IRMPD) is nothing but an isomer-specific multiple reaction monitoring (MRM) method. The identification relies on the comparison of DIMS-MS/MS(IRMPD) chromatograms recorded at different wavenumbers. Based on the comparison of IR spectra of the three isomers, it is shown that specific depletion of the two protonated α- and β-alanine can be achieved, thus allowing for clear identification of the sarcosine peak. It is also demonstrated that DIMS-MS/MS(IRMPD) spectra in the carboxylic C=O stretching region allow for the resolution of overlapping DIMS peaks. Graphical Abstract ᅟ.

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.

Threshold for shattering fragmentation in collision-induced dissociation of the doubly protonated tripeptide TIK(H+)2

Determination of shattering threshold for unimolecular dissociation of a model tripeptide.

Quantum chemical calculation of electron ionization mass spectra for general organic and inorganic molecules

We introduce a fully stand-alone version of the Quantum Chemistry Electron Ionization Mass Spectra (QCEIMS) program [S. Grimme, Angew. Chem. Int. Ed., 2013, 52, 6306] allowing efficient simulations for molecules composed of elements with atomic numbers up to Z = 86. The recently developed extended tight-binding semi-empirical method GFN-xTB has been combined with QCEIMS, thereby eliminating dependencies on third-party electronic structure software. Furthermore, for reasonable calculations of ionization potentials, as required by the method, a second tight-binding variant, IPEA-xTB, is introduced here. This novel combination of methods allows the automatic, fast and reasonably accurate computation of electron ionization mass spectra for structurally different molecules across the periodic table. In order to validate and inspect the transferability of the method, we perform large-scale simulations for some representative organic, organometallic, and main-group inorganic systems. Theoretical spectra for 23 molecules are compared directly to experimental data taken from standard databases. For the first time, realistic quantum chemistry based EI-MS for organometallic systems like ferrocene or copper(ii)acetylacetonate are presented. Compared to previously used semiempirical methods, GFN-xTB is faster, more robust, and yields overall higher quality spectra. The partially analysed theoretical reaction and fragmentation mechanisms are chemically reasonable and reveal in unprecedented detail the extreme complexity of high energy gas phase ion chemistry including complicated rearrangement reactions prior to dissociation.

Classical Dynamics Simulations of Dissociation of Protonated Tryptophan in the Gas Phase

Gas phase decomposition of protonated amino acids are of great interest due to their role in understanding protein and peptide chemistry. Several experimental and theoretical studies have been reported in the literature on this subject. In the present work, decomposition of the aromatic amino acid protonated tryptophan was studied by on-the-fly classical chemical dynamics simulations using density functional theory. Mass spectrometry and electronic structure theory studies have shown multiple dissociation pathways for this biologically relevant molecule. Unlike aliphatic amino acids, protonated tryptophan dissociates via NH₃ elimination rather than the usual iminium ion formation by combined removal of H₂O and CO molecules. Also, a major fragmentation pathway in the present work involves C_α-C_β bond fission. Results of the chemical dynamics simulations reported here are in overall agreement with experiments, and detailed atomic level mechanisms are presented.

On-the-Fly ab Initio Semiclassical Calculation of Glycine Vibrational Spectrum

We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50 000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.

A direct dynamics study of the deprotonated guanine·cytosine base pair: intra-base pair proton transfer, thermal dissociation vs. collision-induced dissociation, and comparison with experiment

Dynamics of intra-base pair hydrogen bonds upon random thermal excitation vs. non-random collisional activation.

Dynamics of Protonated Peptide Ion Collisions with Organic Surfaces: Consonance of Simulation and Experiment

In this Perspective, mass spectrometry experiments and chemical dynamics simulations are described that have explored the atomistic dynamics of protonated peptide ions, peptide-H(+), colliding with organic surfaces. These studies have investigated the energy transfer and fragmentation dynamics for peptide-H(+) surface-induced dissociation (SID), peptide-H(+) physisorption on the surface, soft landing (SL), and peptide-H(+) reaction with the surface, reactive landing (RL). SID provides primary structures of biological ions and information regarding their fragmentation pathways and energetics. Two SID mechanisms are found for peptide-H(+) fragmentation. A traditional mechanism in which peptide-H(+) is vibrationally excited by its collision with the surface, rebounds off the surface and then dissociates in accord with the statistical, RRKM unimolecular rate theory. The other, shattering, is a nonstatistical mechanism in which peptide-H(+) fragments as it collides with the surface, dissociating via many pathways and forming many product ions. Shattering is important for collisions with diamond and perfluorinated self-assembled monolayer (F-SAM) surfaces, increasing in importance with the peptide-H(+) collision energy. Chemical dynamics simulations also provide important mechanistic insights on SL and RL of biological ions on surfaces. The simulations indicate that SL occurs via multiple mechanisms consisting of sequences of peptide-H(+) physisorption on and penetration in the surface. SL and RL have a broad range of important applications including preparation of protein or peptide microarrays, development of biocompatible substrates and biosensors, and preparation of novel synthetic materials, including nanomaterials. An important RL mechanism is intact deposition of peptide-H(+) on the surface.

Mechanistic details of energy transfer and soft landing in ala2-H(+) collisions with a F-SAM surface

Previous chemical dynamics simulations (Phys. Chem. Chem. Phys., 2014, 16, 23769-23778) were analyzed to delineate atomistic details for collision of N-protonated dialanine (ala2-H(+)) with a C8 perfluorinated self-assembled monolayer (F-SAM) surface. Initial collision energies Ei of 5-70 eV and incident angles θi of 0° and 45°, with the surface normal, were considered. Four trajectory types were identified: (1) direct scattering; (2) temporary sticking/physisorption on top of the surface; (3) temporary penetration of the surface with additional physisorption on the surface; and (4) trapping on/in the surface, by physisorption or surface penetration, when the trajectory is terminated. Direct scattering increases from 12 to 100% as Ei is increased from 5 to 70 eV. For the direct scattering at 70 eV, at least one ala2-H(+) heavy atom penetrated the surface for all of the trajectories. For ∼33% of the trajectories all eleven of the ala2-H(+) heavy atoms penetrated the F-SAM at the time of deepest penetration. The importance of trapping decreased with increase in Ei, decreasing from 84 to 0% with Ei increase from 5 to 70 eV at θi = 0°. Somewhat surprisingly, the collisional energy transfers to the F-SAM surface and ala2-H(+) are overall insensitive to the trajectory type. The energy transfer to ala2-H(+) is primarily to vibration, with the transfer to rotation ∼10% or less. Adsorption and then trapping of ala2-H(+) is primarily a multi-step process, and the following five trapping mechanisms were identified: (i) physisorption-penetration-physisorption (phys-pen-phys); (ii) penetration-physisorption-penetration (pen-phys-pen); (iii) penetration-physisorption (pen-phys); (iv) physisorption-penetration (phys-pen); and (v) only physisorption (phys). For Ei = 5 eV, the pen-phys-pen, pen-phys, phys-pen, and phys trapping mechanisms have similar probabilities. For 13.5 eV, the phys-pen mechanism, important at 5 eV, is unimportant. The radius of gyration of ala2-H(+) was calculated once it is trapped on/in the F-SAM surface and trapping decreases the ion's compactness, in part by breaking hydrogen bonds. The ala2-H(+) + F-SAM simulations are compared with the penetration and trapping dynamics found in previous simulations of projectile + organic surface collisions.

Unimolecular dissociation of peptides: statistical vs. non-statistical fragmentation mechanisms and time scales

In the present work we have investigated mechanisms of gas phase unimolecular dissociation of a relatively simple dipeptide, the di-proline anion, by means of chemical dynamics simulations, using the PM3 semi-empirical Hamiltonian. In particular, we have considered two activation processes that are representative limits of what occurs in collision induced dissociation experiments: (i) thermal activation, corresponding to several low energy collisions, in which the system is prepared with a microcanonical distribution of energy; (ii) collisional activation where a single shock of hundreds of kcal mol⁻¹ (300 kcal mol⁻¹ in the present case) can transfer sufficient energy to allow dissociation. From these two activation processes we obtained different product abundances, and for one particular fragmentation pathway a clear mechanistic difference for the two activation processes. This mechanism corresponds to the leaving of an OH⁻ group and subsequent formation of water by taking a proton from the remaining molecule. This last reaction is always observed in thermal activation while in collisional activation it is less favoured and the formation of OH⁻ as a final product is observed. More importantly, we show that while in thermal activation unimolecular dissociation follows exponential decay, in collision activation the initial population decays with non-exponential behaviour. Finally, from the thermal activation simulations it was possible to obtain rate constants as a function of temperature that show Arrhenius behaviour. Thus activation energies have also been extracted from these simulations.

Chemical dynamics simulations of energy transfer, surface-induced dissociation, soft-landing, and reactive-landing in collisions of protonated peptide ions with organic surfaces

Different simulation approaches like MM, QM + MM, and QM/MM, were used to study surface-induced dissociation, soft-landing, and reactive-landing for the peptide-H⁺ + surface collisions.

Elucidating collision induced dissociation products and reaction mechanisms of protonated uracil by coupling chemical dynamics simulations with tandem mass spectrometry experiments

In this study we have coupled mixed quantum-classical (quantum mechanics/molecular mechanics) direct chemical dynamics simulations with electrospray ionization/tandem mass spectrometry experiments in order to achieve a deeper understanding of the fragmentation mechanisms occurring during the collision induced dissociation of gaseous protonated uracil. Using this approach, we were able to successfully characterize the fragmentation pathways corresponding to ammonia loss (m/z 96), water loss (m/z 95) and cyanic or isocyanic acid loss (m/z 70). Furthermore, we also performed experiments with isotopic labeling completing the fragmentation picture. Remarkably, fragmentation mechanisms obtained from chemical dynamics simulations are consistent with those deduced from isotopic labeling.

A RRKM study and a DFT assessment on gas-phase fragmentation of formamide-M(2+) (M = Ca, Sr)

A kinetic study of the unimolecular reactivity of formamide-M(2+) (M = Ca, Sr) systems was carried out by means of RRKM statistical theory using high-level DFT. The results predict M(2+), [M(NH2)](+) and [HCO](+) as the main products, together with an intermediate that could eventually evolve to produce [M(NH3)](2+) and CO, for high values of internal energy. In this framework, we also evaluated the influence of the external rotational energy on the reaction rate constants. In order to find a method to perform reliable electronic structure calculations for formamide-M(2+) (M = Ca, Sr) at a relatively low computational cost, an assessment of different methods was performed. In the first assessment twenty-one functionals, belonging to different DFT categories, and an MP2 wave function method using a small basis set were evaluated. CCSD(T)/cc-pWCVTZ single point calculations were used as reference. A second assessment has been performed on geometries and energies. We found BLYP/6-31G(d) and G96LYP/6-31+G(d,p) as the best performing methods, for formamide-Ca(2+) and formamide-Sr(2+), respectively. Furthermore, a detailed assessment was done on RRKM reactivity and G96LYP/6-31G(d) provided results in agreement with higher level calculations. The combination of geometrical, energetics and kinetics (RRKM) criteria to evaluate DFT functionals is rather unusual and provides an original assessment procedure. Overall, we suggest using G96LYP as the best performing functional with a small basis set for both systems.

Potential energy surfaces for the HBr(+) + CO2 → Br + HOCO(+) reaction in the HBr(+) (2)Π3/2 and (2)Π1/2 spin-orbit states

Quantum mechanical (QM) + molecular mechanics (MM) models are developed to represent potential energy surfaces (PESs) for the HBr(+) + CO2 → Br + HOCO(+) reaction with HBr(+) in the (2)Π3/2 and (2)Π1/2 spin-orbit states. The QM component is the spin-free PES and spin-orbit coupling for each state is represented by a MM-like analytic potential fit to spin-orbit electronic structure calculations. Coupled-cluster single double and perturbative triple excitation (CCSD(T)) calculations are performed to obtain "benchmark" reaction energies without spin-orbit coupling. With zero-point energies removed, the "experimental" reaction energy is 44 ± 5 meV for HBr(+)((2)Π3/2) + CO2 → Br((2)P3/2) + HOCO(+), while the CCSD(T) value with spin-orbit effects included is 87 meV. Electronic structure calculations were performed to determine properties of the BrHOCO(+) reaction intermediate and [HBr⋯OCO](+) van der Waals intermediate. The results of different electronic structure methods were compared with those obtained with CCSD(T), and UMP2/cc-pVTZ/PP was found to be a practical and accurate QM method to use in QM/MM direct dynamics simulations. The spin-orbit coupling calculations show that the spin-free QM PES gives a quite good representation of the shape of the PES originated by (2)Π3/2HBr(+). This is also the case for the reactant region of the PES for (2)Π1/2 HBr(+), but spin-orbit coupling effects are important for the exit-channel region of this PES. A MM model was developed to represent these effects, which were combined with the spin-free QM PES.

Dynamics of energy transfer and soft-landing in collisions of protonated dialanine with perfluorinated self-assembled monolayer surfaces

Chemical dynamics simulations are reported which provide atomistic details of collisions of protonated dialanine, ala2-H(+), with a perfluorinated octanethiolate self-assembled monolayer (F-SAM) surface. The simulations are performed at collision energies Ei of 5.0, 13.5, 22.5, 30.00, and 70 eV, and incident angles 0° (normal) and 45° (grazing). Excellent agreement with experiment (J. Am. Chem. Soc., 2000, 122, 9703-9714) is found for both the average fraction and distribution of the collision energy transferred to the ala2-H(+) internal degrees of freedom. The dominant pathway for this energy transfer is to ala2-H(+) vibration, but for Ei = 5.0 eV ∼20% of the energy transfer is to ala2-H(+) rotation. Energy transfer to ala2-H(+) rotation decreases with increase in Ei and becomes negligible at high Ei. Three types of collisions are observed in the simulations: i.e. those for which ala2-H(+) (1) directly scatters off the F-SAM surface; (2) sticks/physisorbs on/in the surface, but desorbs within the 10 ps numerical integration of the simulations; and (3) remains trapped (i.e. soft-landed) on/in the surface when the simulations are terminated. Penetration of the F-SAM by ala2-H(+) is important for the latter two types of events. The trapped trajectories are expected to have relatively long residence times on the surface, since a previous molecular dynamics simulation (J. Phys. Chem. B, 2014, 118, 5577-5588) shows that thermally accommodated ala2-H(+) ions have an binding energy with the F-SAM surface of at least ∼15 kcal mol(-1).

Intermolecular potential for binding of protonated peptide ions with perfluorinated hydrocarbon surfaces

An analytic potential energy function was developed to model both short-range and long-range interactions between protonated peptide ions and perfluorinated hydrocarbon chains. The potential function is defined as a sum of two-body potentials of the Buckingham form. The parameters of the two-body potentials were obtained by fits to intermolecular potential energy curves (IPECs) calculated for CF4, which represents the F and C atoms of a perfluoroalkane chain, interacting with small molecules chosen as representatives of the main functional groups and atoms present in protonated peptide ions: specifically, CH4, NH3, NH4(+), and HCOOH. The IPECs were calculated at the MP2/aug-cc-pVTZ level of theory, with basis set superposition error (BSSE) corrections. Good fits were obtained for an energy range extending up to about 400 kcal/mol. It is shown that the pair potentials derived from the NH3/CF4 and HCOOH/CF4 fits reproduce acceptably well the intermolecular interactions in HCONH2/CF4, which indicates that the parameters obtained for the amine and carbonyl atoms may be transferable to the corresponding atoms of the amide group. The derived potential energy function may be used in chemical dynamics simulations of collisions of peptide-H(+) ions with perfluorinated hydrocarbon surfaces.

Collision induced dissociation of doubly-charged ions: Coulomb explosion vs. neutral loss in [Ca(urea)]2+ gas phase unimolecular reactivity via chemical dynamics simulations

In this paper we report different theoretical approaches to study the gas-phase unimolecular dissociation of the doubly-charged cation [Ca(urea)](2+), in order to rationalize recent experimental findings. Quantum mechanical plus molecular mechanical (QM/MM) direct chemical dynamics simulations were used to investigate collision induced dissociation (CID) and rotational-vibrational energy transfer for Ar + [Ca(urea)](2+) collisions. For the picosecond time-domain of the simulations, both neutral loss and Coulomb explosion reactions were found and the differences in their mechanisms elucidated. The loss of neutral urea subsequent to collision with Ar occurs via a shattering mechanism, while the formation of two singly-charged cations follows statistical (or almost statistical) dynamics. Vibrational-rotational energy transfer efficiencies obtained for trajectories that do not dissociate during the trajectory integration were used in conjunction with RRKM rate constants to approximate dissociation pathways assuming complete intramolecular vibrational energy redistribution (IVR) and statistical dynamics. This statistical limit predicts, as expected, that at long time the most stable species on the potential energy surface (PES) dominate. These results, coupled with experimental CID from which both neutral loss and Coulomb explosion products were obtained, show that the gas phase dissociation of this ion occurs by multiple mechanisms leading to different products and that reactivity on the complicated PES is dynamically complex.