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

Calculation of Mass Spectra with the QCxMS Method for Negatively and Multiply Charged Molecules

Analysis and validation of a mass spectrometry (MS) experiment are usually performed by comparison to reference spectra. However, if references are missing, measured spectra cannot be properly matched. To close this gap, the Quantum Chemical Mass Spectrometry (QCxMS) program has been developed. It enables fully automatic calculations of electron ionization (EI) and positive ion collision-induced dissociation (CID) mass spectra of singly charged molecular ions. In this work, the extension to negative and multiple ion charge for the CID run mode is presented. QCxMS is now capable of calculating structures carrying any charge, without the need for pretabulated fragmentation pathways or machine learning of database spectra. Mass spectra of four single negatively charged and two multiple positively charged organic ions with molecular sizes from 12 to 92 atoms were computed and compared to reference spectra. The underlying Born-Oppenheimer molecular dynamics (MD) calculations were conducted using the semiempirical quantum mechanical GFN2-xTB method, while for some small molecules, ab initio DFT-based MD simulations were performed. Detailed insights into the fragmentation pathways were gained, and the effects of the computed charge assignments on the resulting spectrum are discussed. Especially for the negative ion mode, the influence of the deprotonation site to create the anion was found to be substantial. Doubly charged fragments could successfully be calculated fully automatically for the first time, while higher charged structures introduced severe assignment problems. Overall, this extension of the QCxMS program further enhances its applicability and underlines its value as a sophisticated toolkit for CID-based tandem MS structure elucidation.

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.

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.

Unimolecular Fragmentation Properties of Thermometer Ions from Chemical Dynamics Simulations

Thermometer ions are widely used to calibrate the internal energy of the ions produced by electrospray ionization in mass spectrometry. Typically, benzylpyridinium ions with different substituents are used. More recently, benzhydrylpyridinium ions were proposed for their lower bond dissociation energies. Direct dynamics simulations using M06-2X/6-31G(d), DFTB, and PM6-D3 are performed to characterize the activation energies of two representative systems: para-methylbenzylpyridinium ion (p-Me-BnPy⁺) and methyl,methylbenzhydrylpyridinium ion (Me,Me-BhPy⁺). Simulation results are used to calculate rate constants for the two systems. These rate constants and their uncertainties are used to find the Arrhenius activation energies and RRK fitted threshold energies which give reasonable agreement with calculated bond dissociation energies at the same level of theory. There is only one fragmentation mechanism observed for both systems, which involves C-N bond dissociation via a loose transition state, to generate either benzylium or benzhydrylium ion and a neutral pyridine molecule. For p-Me-BnPy⁺ using DFTB and PM6-D3 the formation of tropylium ion, from rearrangement of benzylium ion, was observed but only at higher excitation energies and for longer simulation times. These observations suggest that there is no competition between reaction pathways that could affect the reliability of internal energy calibrations. Finally, we suggest using DFTB with a modified-Arrhenius model in future studies.

Observation of Enhanced Dissociative Photochemistry in the Non-Native Nucleobase 2-Thiouracil

We present the first study to measure the dissociative photochemistry of 2-thiouracil (2-TU), an important nucleobase analogue with applications in molecular biology and pharmacology. Laser photodissociation spectroscopy is applied to the deprotonated and protonated forms of 2-TU, which are produced in the gas-phase using electrospray ionization mass spectrometry. Our results show that the deprotonated form of 2-thiouracil ([2-TU-H]⁻) decays predominantly by electron ejection and hence concomitant production of the [2-TU-H]· free-radical species, following photoexcitation across the UVA-UVC region. Thiocyanate (SCN⁻) and a m/z 93 fragment ion are also observed as photodecay products of [2-TU-H]⁻ but at very low intensities. Photoexcitation of protonated 2-thiouracil ([2-TU·H]⁺) across the same UVA-UVC spectral region produces the m/z 96 cationic fragment as the major photofragment. This ion corresponds to ejection of an HS· radical from the precursor ion and is determined to be a product of direct excited state decay. Fragment ions associated with decay of the hot ground state (i.e., the ions we would expect to observe if 2-thiouracil was behaving like UV-dissipating uracil) are observed as much more minor products. This behaviour is consistent with enhanced intersystem crossing to triplet excited states compared to internal conversion back to the ground state. These are the first experiments to probe the effect of protonation/deprotonation on thionucleobase photochemistry, and hence explore the effect of pH at a molecular level on their photophysical properties.

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.

Collision-Induced Dissociation Studies of Protonated Ions of Alkylated Thymidine and 2'-Deoxyguanosine

Mass spectrometry and tandem MS (MS/MS) have been widely employed for the identification and quantification of damaged nucleosides in DNA, including those induced by alkylating agents. Upon collisional activation, protonated ions of alkylated nucleosides frequently undergo facile neutral loss of a 2-deoxyribose in MS/MS, and further cleavage of the resulting protonated nucleobases in MS³ can sometimes be employed for differentiating regioisomeric alkylated DNA lesions. Herein, we investigated systematically the collision-induced dissociation (CID) of the protonated ions of O⁴-alkylthymidine (O⁴-alkyldT), O²-alkyldT, O⁶-alkyl-2'-deoxyguanosine (O⁶-alkyldG), and N²-alkyldG through MS³ analysis. The MS³ of O²- and O⁴-MedT exhibit different fragmentation patterns from each other and from other O²- and O⁴-alkyldT adducts carrying larger alkyl groups. Meanwhile, elimination of alkene via a six-membered ring transition state is the dominant fragmentation pathway for O²-alkyldT, O⁴-alkyldT, and O⁶-alkyldG adducts carrying larger alkyl groups, whereas O⁶-MedG mainly undergoes elimination of ammonia. The breakdown of N²-alkyldG is substantially influenced by the structure of the alkyl group, where the relative ease in eliminating ammonia and alkene is modulated by the chain length and branching of the alkyl groups. We also rationalize our observations with density functional theory (DFT) calculations.

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

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.

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.

Gas-Phase Oxidation of the Protonated Uracil-5-yl Radical Cation

This study targets the kinetics and product detection of the gas-phase oxidation reaction of the protonated 5-dehydrouracil (uracil-5-yl) distonic radical cation using ion-trap mass spectrometry. Protonated 5-dehydrouracil radical ions (5-dehydrouracilH⁺ radical ion, m/z 112) are produced within an ion trap by laser photolysis of protonated 5-iodouracil. Storage of the 5-dehydrouracilH⁺ radical ion in the presence of controlled concentration of O₂ reveals two main products. The major reaction product pathway is assigned as the formation of protonated 2-hydroxypyrimidine-4,5-dione (m/z 127) + ^•OH. A second product ion (m/z 99), putatively assigned as a five-member-ring ketone structure, is tentatively explained as arising from the decarbonylation (-CO) of protonated 2-hydroxypyrimidine-4,5-dione. Because protonation of the 5-dehydrouracil radical likely forms a dienol structure, the O₂ reaction at the 5 position is ortho to an -OH group. Following this addition of O₂, the peroxyl-radical intermediate isomerizes by H atom transfer from the -OH group. The ensuing hydroperoxide then decomposes to eliminate ^•OH radical. It is shown that this elimination of ^•OH radical (-17 Da) is evidence for the presence of an -OH group ortho to the initial phenyl radical site, in good accord with calculations. The subsequent CO loss mechanism, to form the aforementioned five-member-ring structure, is unclear, but some pathways are discussed. By following the kinetics of the reaction, the room temperature second-order rate coefficient of the 5-dehydrouracilH⁺ distonic radical cation with molecular oxygen is measured at 7.2 × 10^-11 cm³ molecule^-1 s^-1, Φ = 12% (with ±50% total accuracy). For aryl radical reactions with O₂, the presence of the ^•OH elimination product pathway, following the peroxyl-radical formation, is an indicator of an -OH group ortho to the radical site.

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.

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

Gas-phase interactions of organotin compounds with cysteine

The gas-phase interactions of cysteine with di-organotin and tri-organotin compounds have been studied by mass spectrometry experiments and quantum calculations. Positive-ion electrospray spectra show that the interaction of di- and tri-organotins with cysteine results in the formation of [(R)₂ Sn(Cys-H)]⁺ and [(R)₃ Sn(Cys)]⁺ ions, respectively. MS/MS spectra of [(R)₂ Sn(Cys-H)]⁺ complexes are characterized by numerous fragmentation processes, notably associated with elimination of NH₃ and (C,H₂ ,O₂ ). Several dissociation routes are characteristic of each given organic species. Upon collision, both the [(R)₃ Sn(Gly)]⁺ and [(R)₃ Sn(Cys)]⁺ complexes are associated with elimination of the intact amino acid, leading to the formation of [(R)₃ Sn]⁺ cation. But for the latter complex, two additional fragmentation processes are observed, associated with the elimination of NH₃ and C₃ H₄ O₂ S. Calculations indicate that the interaction between organotins and cysteine is predominantly electrostatic but also exhibits a considerable covalent character, which is slightly more pronounced in tri-organotin complexes. A preferred bidentate interaction of the type -η² -S-NH₂ , with sulfur and the amino group, is observed. As for the [(R)₃ Sn(Cys)]⁺ complexes, their stability is due to the combination of the hydrogen bond taking place between the amino group and the sulfur lone pair and the interaction between the carboxylic oxygen atom and the metal. Copyright © 2016 John Wiley & Sons, Ltd.

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.

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.