Protonated Urea Collision-Induced Dissociation. Comparison of Experiments and Chemical Dynamics Simulations

Robustness of Threshold Collision-Induced Dissociation Simulations for Bond Dissociation Energies

The threshold collision-induced dissociation (T-CID) method is the workhorse for gas-phase bond dissociation energy (BDE) measurements. However, T-CID does not measure BDEs directly; instead, BDEs are obtained by fitting simulated data to the experimental data. We previously observed several large discrepancies between the computed and experimental BDEs. To analyze the reliability of the experimental values, we previously reported a study of the dissociation rate models in the simulation. Here, we report a study of the collision simulation part, specifically in the L-CID (ligand CID) program. We show that the BDE values are robust even to intentionally introduced mistakes in the simulations, varying in most cases by less than 3 kcal mol-1. The most significant exception is the collisional energy transfer (CET) simulation, which led to deviations larger than 10 kcal mol-1. However, we found that the BDEs obtained with explicitly simulated CET distributions deviated by only 3 kcal mol-1 from those simulated with the original model. Collectively, our results suggest that the T-CID-derived BDE values are robust and are likely to be accurate.

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.

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.

Current Challenges and Recent Developments in Mass Spectrometry-Based Metabolomics

High-resolution mass spectrometry (MS) has advanced the study of metabolism in living systems by allowing many metabolites to be measured in a single experiment. Although improvements in mass detector sensitivity have facilitated the detection of greater numbers of analytes, compound identification strategies, feature reduction software, and data sharing have not kept up with the influx of MS data. Here, we discuss the ongoing challenges with MS-based metabolomics, including de novo metabolite identification from mass spectra, differentiation of metabolites from environmental contamination, chromatographic separation of isomers, and incomplete MS databases. Because of their popularity and sensitive detection of small molecules, this review focuses on the challenges of liquid chromatography-mass spectrometry-based methods. We then highlight important instrumentational, experimental, and computational tools that have been created to address these challenges and how they have enabled the advancement of metabolomics research.

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.

Gas-Phase Synthesis of 3-Vinylcyclopropene via the Crossed Beam Reaction of the Methylidyne Radical (CH; X2 Π) with 1,3-Butadiene (CH2 CHCHCH2 ; X1 Ag )

The crossed molecular beam reactions of the methylidyne radical (CH; X² Π) with 1,3-butadiene (CH₂ CHCHCH₂ ; X¹ A_g ) along with their (partially) deuterated counterparts were performed at collision energies of 20.8 kJ mol^-1 under single collision conditions. Combining our laboratory data with ab initio calculations, we reveal that the methylidyne radical may add barrierlessly to the terminal carbon atom and/or carbon-carbon double bond of 1,3-butadiene, leading to doublet C₅ H₇ intermediates with life times longer than the rotation periods. These collision complexes undergo non-statistical unimolecular decomposition through hydrogen atom emission yielding the cyclic cis- and trans-3-vinyl-cyclopropene products with reaction exoergicities of 119±42 kJ mol^-1 . Since this reaction is barrierless, exoergic, and all transition states are located below the energy of the separated reactants, these cyclic C₅ H₆ products are predicted to be accessed even in low-temperature environments, such as in hydrocarbon-rich atmospheres of planets and cold molecular clouds such as TMC-1.

Insights into Cisplatin Binding to Uracil and Thiouracils from IRMPD Spectroscopy and Tandem Mass Spectrometry

The monofunctional primary complexes cis-[PtCl(NH₃)₂(L)]⁺, formed by the reaction of cisplatin, a major chemotherapeutic agent, with four nucleobases L, i.e., uracil (U), 2-thiouracil (2SU), 4-thiouracil (4SU), and 2,4-dithiouracil (24dSU), have been studied by a combination of infrared multiple photon dissociation (IRMPD) action spectroscopy in both the fingerprint (900-1900 cm^-1) and the N-H/O-H stretching (3000-3800 cm^-1) ranges, energy-resolved collision-induced dissociation (CID) mass spectrometry, and density functional calculations at the B3LYP/LACVP/6-311G** level. On the basis of the comparison across the experimental features and the linear IR spectra of conceivable structures, the cisplatin residue is found to promote a monodentate interaction preferentially with the O4(S4) atoms of the canonical forms of U, 4SU, and 24dSU and to the S2 atom of 2SU, yielding the most stable structures. Additional absorptions reveal the presence of minor, alternative tautomers in the sampled ion populations of 2SU and 24dSU, underlying the ability of cisplatin to increase the prospect of (therapeutically beneficial) nucleic acid strand disorder. Implication of these evidence may provide insights into drug mechanism and design.

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.

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

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.

Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces

In this Introduction, we show the basic problems of non-statistical and non-equilibrium phenomena related to the papers collected in this themed issue. Over the past few years, significant advances in both computing power and development of theories have allowed the study of larger systems, increasing the time length of simulations and improving the quality of potential energy surfaces. In particular, the possibility of using quantum chemistry to calculate energies and forces 'on the fly' has paved the way to directly study chemical reactions. This has provided a valuable tool to explore molecular mechanisms at given temperatures and energies and to see whether these reactive trajectories follow statistical laws and/or minimum energy pathways. This themed issue collects different aspects of the problem and gives an overview of recent works and developments in different contexts, from the gas phase to the condensed phase to excited states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Gas-phase reactivity of [Ca(formamide)]2+ complex: an example of different dynamical behaviours

In the present contribution, we have summarized our recent work on the comprehension of [Ca(formamide)]²⁺ complex gas-phase unimolecular dissociation. By using different theoretical approaches, we were able to revise the original (and typical for such kind of problems) understanding given in terms of stationary points on the potential energy surface, which did not provide a satisfactory explanation of the experimentally observed reactivity. In particular, we point out how non-statistical and non-intrinsic reaction coordinate mechanisms are of fundamental importance.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

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.

The chemistry of high valent tungsten chlorides with N-substituted ureas, including urea self-protonation reactions triggered by WCl6

N-Substituted ureas are oxidized by WCl₆in dichloromethane, resulting in protonated urea saltsviaactivation of C–H and/or N–H bonds.

On the gas phase fragmentation of protonated uracil: a statistical perspective

The major fragmentation pathways of protonated uracil calculated in this work using statistical methods agree with mass spectrometry experiments.

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.

Gas-phase basicities of polyfunctional molecules. Part 4: Carbonyl groups as basic sites

This article constitutes the fourth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435). This fourth part is devoted to carbonyl containing polyfunctional molecules. After a short reminder of the methods of determination of gas-phase basicity and the underlying physicochemical concepts, specific examples are examined under two major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) ketones, diketones, ketoalcohols, and ketoethers are considered. A second chapter describes the protonation energetic of gaseous acids and derivatives including diacids, diesters, diamides, anhydrides, imides, ureas, carbamates, amino acid derivatives, and peptides. Experimental data were re-evaluated according to the presently adopted basicity scale. Structural and energetic information given by G3 and G4 quantum chemistry computations on typical systems are presented.

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.

A crystallographically characterized salt of self-generated N-protonated tetraethylurea

The first crystallographically characterizedN-protonated urea (1,1,3,3-tetraethylurea, teu) was obtained by a WCl₆-directed electron transfer/C–H activation process.

Comparing femtosecond multiphoton dissociative ionization of tetrathiafulvene with imaging photoelectron photoion coincidence spectroscopy

In this paper we describe femtosecond photoionization and the imaging photoelectron photoion coincidence spectroscopy of tetrathiafulvene, TTF. Femtosecond photoionization of TTF results in the absorption of up to twelve 808 nm photons leading to ion internal energies up to 12.1 eV as deduced from the photoelectron spectrum. Within this internal energy a variety of dissociation channels are accessible. In order to disentangle the complex ionic dissociation, we utilized the imaging photoelectron photoion coincidence (iPEPICO) technique. Above the dissociation threshold, iPEPICO results show that the molecular ion (m/z = 204) dissociates into seven product ions, six of which compete in a 1.0 eV internal energy window and are formed with the same appearance energy. Ab initio calculations are reported on the possible fragment ion structures of five dissociation channels as well as trajectories showing the loss of C2H2 and C2H2S from high internal energy TTF cations. A three-channel dissociation model is used to fit the PEPICO data in which two dissociation channels are treated as simple dissociations (one with a reverse barrier), while the rest involve a shared barrier. The two lower energy dissociation channels, m/z = 146 and the channel leading to m/z = 178, 171, 159, 140, and 127, have E0 values of 2.77 ± 0.10 and 2.38 ± 0.10 eV, respectively, and are characterized by ΔS(‡)(600 K) values of -9 ± 6 and 1 ± 6 J K(-1) mol(-1), respectively. Competing with them at higher internal energy is the cleavage of the central bond to form the m/z = 102 fragment ion, with an E0 value of 3.65 ± 0.10 eV and ΔS(‡)(600 K) = 83 ± 10 J K(-1) mol(-1).

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.