Statistical Rate Theory and Kinetic Energy-Resolved Ion Chemistry: Theory and Applications

Impact of Protonation Sites on Collision-Induced Dissociation-MS/MS Using CIDMD Quantum Chemistry Modeling

Protonation is the most frequent adduct found in positive electrospray ionization collision-induced mass spectra (CID-MS/MS). In a parallel report Lee, J. J. Chem. Inf. Model. 2024, 10.1021/acs.jcim.4c00760, we developed a quantum chemistry framework to predict mass spectra by collision-induced dissociation molecular dynamics (CIDMD). As different protonation sites affect fragmentation pathways of a given molecule, the accuracy of predicting tandem mass spectra by CIDMD ultimately depends on the choice of its protomers. To investigate the impact of molecular protonation sites on MS/MS spectra, we compared CIDMD-predicted spectra to all available experimental MS/MS spectra by similarity matching. We probed 10 molecules with a total of 43 protomers, the largest study to date, including organic acids (sorbic acid, citramalic acid, itaconic acid, mesaconic acid, citraconic acid, and taurine) as well as aromatic amines including uracil, aniline, bufotenine, and psilocin. We demonstrated how different protomers can converge different fragmentation pathways to the same fragment ions but also may explain the presence of different fragment ions in experimental MS/MS spectra. For the first time, we used in silico MS/MS predictions to test the impact of solvents on proton affinities, comparing the gas phase and a mixture of acetonitrile/water (1:1). We also extended applications of in silico MS/MS predictions to investigate the impact of protonation sites on the energy barriers of isomerization between protomers via proton transfer. Despite our initial hypothesis that the thermodynamically most stable protomer should give the best match to the experiment, we found only weak inverse relationships between the calculated proton affinities and corresponding entropy similarities of experimental and CIDMD-predicted MS/MS spectra. CIDMD-predicted mechanistic details of fragmentation reaction pathways revealed a clear preference for specific protomer forms for several molecules. Overall, however, proton affinity was not a good predictor corresponding to the predicted CIDMD spectra. For example, for uracil, only one protomer predicted all experimental MS/MS fragment ions, but this protomer had neither the highest proton affinity nor the best MS/MS match score. Instead of proton affinity, the transfer of protons during the electrospray process from the initial protonation site (i.e., mobile proton model) better explains the differences between the thermodynamic rationale and experimental data. Protomers that undergo fragmentation with lower energy barriers have greater contributions to experimental MS/MS spectra than their thermodynamic Boltzmann populations would suggest. Hence, in silico predictions still need to calculate MS/MS spectra for multiple protomers, as the extent of distributions cannot be readily predicted.

Perspective: intrinsic interactions of metal ions with biological molecules as studied by threshold collision-induced dissociation and infrared multiple photon dissociation

In this perspective, gas-phase studies of group 1 monocations and group 12 dications with amino acids and small peptides are highlighted. Although the focus is on two experimental techniques, threshold collision-induced dissociation and infrared multiple photon dissociation action spectroscopy, these methods as well as complementary approaches are summarized. The synergistic interplay with theory, made particularly powerful by the small sizes of the systems explored and the absence of solvent and support, is also elucidated. Importantly, these gas-phase methods permit quantitative insight into the structures and thermodynamics of metal cations interacting with biological molecules. Periodic trends in how these interactions vary as the metal cations get heavier are discussed as are quantitative trends with changes in the amino acid side chain and effects of hydration. Such trends allow these results to transcend the limitations associated with the biomimetic model systems.

Magic cluster sizes of cationic and anionic sodium chloride clusters explained by statistical modeling of the complete phase space

As one of the main components of sea salt aerosols, sodium chloride is involved in numerous atmospheric processes. Gas-phase clusters are ideal models to study fundamental physical and chemical properties of sodium chloride, which are significantly affected by the cluster size. Of particular interest are magic cluster sizes, which exhibit high intensities in mass spectra. In order to understand the origin of these magic cluster sizes, quantum chemical calculations at the CCSD(T)//DFT level are performed, yielding structures and binding energies of neutral (NaCl)x, anionic (NaCl)xCl- and cationic (NaCl)xNa+ clusters up to x = 8. Our calculations show that the clusters can easily isomerize, enabling dissociation into the lowest-energy isomers of the fragments. Energetics can explain the special stability of (NaCl)4Cl-, but (NaCl)4Na+ actually offers low-lying dissociation channels, despite being a magic cluster size. Collision-induced dissociation experiments reveal that the loss of neutral clusters (NaCl)x, x = 2, 4, is in most cases more favorable than the loss of NaCl or the atomic ion, i.e. sodium chloride clusters actually fragment via the cleavage of the entire cluster, not by evaporating small cluster building blocks. This is rationalized by the calculated high stability of even-numbered neutral clusters (NaCl)x, especially x = 2, 4. Analysis of the density of states and rate constants calculated with a modified Rice-Ramsperger-Kassel-Marcus (RRKM) equation called AWATAR - considering all energetically accessible isomers of reactants and fragments - shows that entropic effects are responsible for the magic cluster character of (NaCl)4Na+. In particular, low-lying vibrational modes provide a high density of states of the near-planar cluster. Together with the small contribution of an atomic ion to the sum of states in a loose transition state for dissociation, this leads to a very small unimolecular rate constant for dissociation into (NaCl)4 and Na+, which is the lowest energy fragmentation pathway. Thus, entropic effects may override energetics for certain magic cluster sizes.

Very Low-Pressure CID Experiments: High Energy Transfer and Fragmentation Pattern at the Single Collision Regime

We have performed CID experiments on a triple quadrupole instrument, lowering the collision gas pressure by 50 times compared to its conventional value. The results show that at very low-collision gas pressure, single collisions dominate the spectra. Indirectly, these results suggest that under conventional conditions, 20-50 collisions may be typical in CID experiments. The results show a marked difference between low- and high-pressure CID spectra, the latter being characterized in terms of 'slow heating' and predominance of consecutive reactions. The results indicate that under single collision conditions, the collisional energy transfer efficiency is very high: nearly 100% of the center of mass kinetic energy is converted to internal energy.

Threshold collision induced dissociation of protonated water clusters

We report threshold collision induced dissociation experiments on protonated water clusters thermalized at low temperature for sizes n = 19-23. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model allows us to account for dissociation cascades and provides values for the dissociation energies of each cluster. These values, averaging around 0.47 eV, are in good agreement with theoretical predictions at various levels of theory. Furthermore, the dissociation energies show a trend for the n = 21 magic and n = 22 anti-magic numbers relative to their neighbours, which is also in agreement with theory. These results provide further evidence to resolve the disagreement between previously published experimental values. A careful quantitative treatment of cascade dissociation in this model introduces interdependence between the dissociation energies of neighboring sizes, which reduces the number of free fitting parameters and improves both reliability and uncertainties on absolute dissociation energies deduced from experiments.

Non-ergodic fragmentation upon collision-induced activation of cysteine-water cluster cations

Cysteine-water cluster cations Cys(H2O)3,6+ and Cys(H2O)3,6H+ are assembled in He droplets and probed by tandem mass spectrometry with collision-induced activation. Benchmark experimental data for this biologically important system are complemented with theory to elucidate the details of the collision-induced activation process. Experimental energy thresholds for successive release of water are compared to water dissociation energies from DFT calculations showing that clusters do not only fragment exclusively by sequential emission of single water molecules but also by the release of small water clusters. Release of clustered water is observed also in the ADMP (atom centered density matrix propagation) molecular dynamics model of small Cys(H2O)3+ and Cys(H2O)3H+ clusters. For large clusters Cys(H2O)6+ and Cys(H2O)6H+ the less computationally demanding statistical Microcanonical Metropolis Monte-Carlo method (M3C) is used to model the experimental fragmentation patterns. We are able to detail the energy redistribution in clusters upon collision activation. In the present case, about two thirds of the collision energy redistribute via an ergodic process, while the remaining one third is transferred into a non-ergodic channel leading to ejection of a single water molecule from the cluster. In contrast to molecular fragmentation, which can be well described by statistical models, modelling of collision-induced activation of weakly bound clusters requires inclusion of non-ergodic processes.

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

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

Water Attachment onto Size-Selected Cationic Pyrene Clusters

We report measurements of the attachment rates of water molecules onto mass-selected cationic pyrene clusters for size from n = 4 to 13 pyrene units and for different collision energies. Comparison of the attachment rates with the collision rates measured in collision-induced dissociation experiments provides access to the values of the sticking coefficient. The strong dependence of the attachment rates on size and collision energy is rationalized through a model in which we use a Langevin-type collision rate and adjust on experimental data the statistical dissociation rate of the water molecule from the cluster after attachment. This allows us to extrapolate our results to the conditions of isolation and long time scales encountered in astrophysical environments.

Model Compounds for Intermediates and Transition States in Sonogashira and Negishi Coupling: d8-d10 Bonds in Large Heterobimetallic Complexes Are Weaker than Computational Chemistry Predicts

We report an experimental study, with accompanying DFT calculations, on a series of heterobimetallic complexes with Pd(II) and Cu(I), Ag(I), or Au(I). The isolable pincer complexes are models for the intermediates and transition state for the transmetalation step in Sonogashira and Negishi coupling reactions, among which, according to the DFT calculations, only the transition state has the two metal centers within bonding distance. Furthermore, we report a substituted version of an analogous heterobimetallic complex in which a competing dissociation sets an upper-bound on the strength of the d⁸-d¹⁰ metal-metal bond. Analysis of the structures in the solid state and in solution, and the competitive dissociation experiment in the gas phase, indicate that the dispersion-corrected DFT method used in the study appears to overestimate the strength of the metal-metal interaction, thus distorting the shape of the computed potential energy surface systematically for transmetalation.

Reactions of Atomic Thorium and Uranium Cations with SF6 Studied by Guided Ion Beam Tandem Mass Spectrometry

The fundamental chemistry of the thorium and uranium fluorides continues to be an area of interest because of the use of thorium and uranium fluoride compounds in nuclear fuel systems. Here, we study the reaction of thorium cations with sulfur hexafluoride for the first time and revisit the reaction of uranium cations with sulfur hexafluoride. By using guided ion beam tandem mass spectrometry, we explore the reaction pathways that become accessible well above thermal energies (E ∼ 0.04 eV). Overall, we find that both Th⁺ and U⁺ react very efficiently with SF₆, approaching the collision limit at both thermal and elevated energies. The primary products observed at low energies include Th_1-3⁺, UF_1-4⁺, and SF_1-4⁺, all of which are formed in barrierless, exothermic processes. SF₅⁺ was also observed, although the pressure dependence of this channel reveals that SF₅⁺ forms exothermically through secondary reactions, which the energy dependences suggest result from reactions between ThF₂⁺ and UF₃⁺ with SF₆. At higher energies, both AnF₃⁺ products are observed to decay to AnF⁺ + F₂, and both SF₄⁺ and SF₂⁺ exhibit cross sections with endothermic features. For both systems, the rise in SF₄⁺ can be attributed to a secondary collision between AnF⁺ with SF₆ on the basis of the pressure dependence of the SF₄⁺ channel at higher energies, and the rise in SF₂⁺ appears to result from the decomposition of SF₃⁺ to SF₂⁺ + F.

Origin of the different reactivity of the high-valent coinage-metal complexes [RCuiii Me3 ]- and [RAgiii Me3 ]- (R=allyl)

High-valent tetraalkylcuprates(iii) and -argentates(iii) are key intermediates of copper- and silver-mediated C-C coupling reactions. Here, we investigate the previously reported contrasting reactivity of [RMiii Me3 ]- complexes (M=Cu, Ag and R=allyl) with energy-dependent collision-induced dissociation experiments, advanced quantum-chemical calculations and kinetic computations. The gas-phase fragmentation experiments confirmed the preferred formation of the [RCuMe]- anion upon collisional activation of the cuprate(iii) species, consistent with a homo-coupling reaction, whereas the silver analogue primarily yielded [AgMe2 ]- , consistent with a cross-coupling reaction. For both complexes, density functional theory calculations identified one mechanism for homo coupling and four different ones for cross coupling. Of these pathways, an unprecedented concerted outer-sphere cross coupling is of particular interest, because it can explain the formation of [AgMe2 ]- from the argentate(iii) species. Remarkably, the different C-C coupling propensities of the two [RMiii Me3 ]- complexes become only apparent when properly accounting for the multi-configurational character of the wave function for the key transition state of [RAgMe3 ]- . Backed by the obtained detailed mechanistic insight for the gas-phase reactions, we propose that the previously observed cross-coupling reaction of the silver complex in solution proceeds via the outer-sphere mechanism.

Sodium Binding Interactions with Aliphatic Amino Acids: A Guided Ion Beam and Computational Study

Metal binding affinities play a vital role in medicinal, biological, and industrial applications. In particular, metal cation-amino acid (AA) interactions contribute to protein stability such that analyzing analogous prototypical interactions is important. Here, we present a full description of the interactions of sodium cations (Na⁺) and six aliphatic amino acids (AA), where AA = glycine (Gly), alanine (Ala), homoalanine (hAla), valine (Val), leucine (Leu), and isoleucine (Ile). Experimentally, these interactions are evaluated using threshold collision-induced dissociation carried out in a guided ion beam tandem mass spectrometer, allowing for the determination of the kinetic-energy-dependent behavior of Na⁺-AA dissociation. Analysis of these dissociation cross sections, after accounting for multiple ion-molecule collisions, internal energy of reactant ions, and unimolecular decay rates, allows the determination of absolute Na⁺-AA bond dissociation energies (BDEs) in kJ/mol of Gly (164.0), Ala (166.9), hAla (167.9), Val (172.7), Leu (173.7), and Ile (174.6). These are favorably compared to quantum chemical calculations conducted at the B3LYP, B3P86, MP2(full), B3LYP-GD3BJ, and M06-2X levels of theory. Our combination of structural and energetic analyses provides information regarding the specific factors responsible for Na⁺ interactions with amino acids. Specifically, we find that the BDEs increase linearly with increasing polarizability of the amino acid.

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.

Modeling Gas-Phase Unimolecular Dissociation for Bond Dissociation Energies: Comparison of Statistical Rate Models within RRKM Theory

The Rice-Ramsperger-Kassel-Marcus (RRKM) theory provides a simple yet powerful rate theory for calculating microcanonical rate constants. In particular, it has found widespread use in combination with gas-phase kinetic experiments of unimolecular dissociations to extract experimental bond dissociation energies (BDEs). We have previously found several discrepancies between the computed BDE values and the respective experimental ones, obtained with our empirical rate model, named L-CID. To investigate the reliability of our rate model, we conducted a theoretical analysis and comparison of the performance of conventional rate models and L-CID within the RRKM framework. Using the previously published microcanonical rate data as well as reaction cross-section data, we show that the BDE values obtained with the L-CID model agree with the ones from the other rate models within the expected uncertainty bounds. Based on this agreement, we discuss the possible rationalization of the good performance of the L-CID model.

Bond Dissociation Energies in the Gas Phase for Large Molecular Ions by Threshold Collision-Induced Dissociation Experiments: Stretching the Limits

Accurate bond dissociation energies for large molecules are difficult to obtain by either experimental or computational methods. The former methods are hampered by a range of physical and practical limitations in gas-phase measurement techniques, while the latter require incorporation of multiple approximations whose impact on accuracy may not always be clear. When internal benchmarks are not available, one hopes that experiment and theory can mutually support each other. A recent report found, however, a large discrepancy between gas-phase bond dissociation energies, measured mass spectrometrically, and the corresponding quantities computed using density functional theory (DFT)-D3 and DLPNO-CCSD(T) methods. With the widespread application of these computational methods to large molecular systems, the discrepancy needs to be resolved. We report a series of experimental studies that validate the mass spectrometric methods from small to large ions and find that bond dissociation energies extracted from threshold collision-induced dissociation experiments on large ions do indeed behave correctly. The implications for the computational studies are discussed.

Threshold collision induced dissociation of pyrene cluster cations

We report threshold collision induced dissociation experiments on cationic pyrene clusters, for sizes n = 2-6. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model can account for the dissociation cascades and provides values for the dissociation energies. These values, of the order of 0.7 eV-1 eV, are in excellent agreement with those previously derived from thermal evaporation. They confirm the charge resonance stability enhancement predicted by theoretical calculations. In addition, remarkable agreement is obtained with theoretical predictions for the two smaller sizes n = 2 and 3. For the larger sizes, the agreement remains good, although the theoretical values obtained for the most stable structures are systematically higher by 0.2 eV. This offset could be attributed to approximations in the calculations. Still, there is an indication in the results of an incomplete description of the role of isomerization and/or direct dissociation upon collisions. Finally, by-product clusters containing dehydrogenated species are found to dissociate at energies comparable to the non-dehydrogenated ones, which shows no evidence for covalent bonds within the clusters.

Mass spectrometry and computational study of collision-induced dissociation of 9-methylguanine–1-methylcytosine base-pair radical cation: intra-base-pair proton transfer and hydrogen transfer, non-statistical dissociation, and reaction with a water ligand

A combined experimental and theoretical study is presented on the collision-induced dissociation of 9-methylguanine–1-methylcytosine base-pair radical cation ([9MG·1MC]˙⁺) and its monohydrate ([9MG·1MC]˙⁺·H₂O) with Xe and Ar gases.

Collision Induced Dissociation of Benzylpyridinium-Substituted Porphyrins: Towards a Thermometer Scale for Multiply Charged Ions?

We have determined breakdown curves for a range of multiply charged benzylpyridinium-substituted porphyrin cations by collision induced dissociation measurements (CID) as mediated by resonant pulsed radio-frequency (rf) excitation in a helium-filled linear ion trap. Measurements were compared with the predictions of DFT calculations. We find a linear correlation between experimental fragmentation thresholds (in instrumental units of "normalized collision energy") and theoretical dissociation energies, suggesting that these species can be used as calibrants to gauge the fragmentation energetics of closely related systems. We have confirmed this by also studying the fragmentation thresholds of metalloporphyrin-based ions - including multiply negatively charged metalloporphyrin oligomers. Unfortunately, the slope of the linear correlation obtained for benzylpyridinium-substituted porphyrin multications differs significantly from that obtained by us for a set of smaller, singly charged substituted benzylpyridines put forward as "thermometer" ions in previous work. Multiplying the threshold energies in an ad hoc fashion by the ion charge basically reconciles both calibration curves. We conclude that one should use caution when applying small, singly charged benzylpyridines as calibrants to gauge the CID of large, multiply charged ions in ion-trap mass spectrometers. Graphical Abstract.

N-Protonated Isomers and Coulombic Barriers to Dissociation of Doubly Protonated Ala8Arg

Collision-induced dissociation (or tandem mass spectrometry, MS/MS) of a protonated peptide results in a spectrum of fragment ions that is useful for inferring amino acid sequence. This is now commonplace and a foundation of proteomics. The underlying chemical and physical processes are believed to be those familiar from physical organic chemistry and chemical kinetics. However, first-principles predictions remain intractable because of the conflicting necessities for high accuracy (to achieve qualitatively correct kinetics) and computational speed (to compensate for the high cost of reliable calculations on such large molecules). To make progress, shortcuts are needed. Inspired by the popular mobile proton model, we have previously proposed a simplified theoretical model in which the gas-phase fragmentation pattern of protonated peptides reflects the relative stabilities of N-protonated isomers, thus avoiding the need for transition-state information. For singly protonated Ala _n (n = 3-11), the resulting predictions were in qualitative agreement with the results from low-energy MS/MS experiments. Here, the comparison is extended to a model tryptic peptide, doubly protonated Ala₈Arg. This is of interest because doubly protonated tryptic peptides are the most important in proteomics. In comparison with experimental results, our model seriously overpredicts the degree of backbone fragmentation at N₉. We offer an improved model that corrects this deficiency. The principal change is to include Coulombic barriers, which hinder the separation of the product cations from each other. Coulombic barriers may be equally important in MS/MS of all multiply charged peptide ions. Graphical Abstract ᅟ.

Potential energy surface for the reaction Sm+ + CO2 → SmO+ + CO: guided ion beam and theoretical studies

The potential energy surface (PES) for the oxidation of samarium cations by carbon dioxide is explored both experimentally and theoretically. Using guided ion beam tandem mass spectrometry, several reactions are examined as a function of kinetic energy. These include the title reaction as well as its reverse along with the collision-induced dissociation of Sm⁺(CO₂) and OSm⁺(CO) with Xe. Analysis of the kinetic energy dependent cross sections yields barriers for the forward and reverse oxidation reaction of 1.77 ± 0.11 and 2.04 ± 0.13 eV, respectively, and Sm⁺-OCO and OSm⁺-CO bond dissociation energies (BDEs) of 0.42 ± 0.03 and 0.97 ± 0.07 eV, respectively. BDEs for Sm⁺(CO₂)_x for x = 2 and 3 are also determined as 0.40 ± 0.13 and 0.48 ± 0.12 eV, respectively. The PESs for the title reaction along the sextet and octet spin surfaces are also examined theoretically at the MP2 and CCSD(T) levels using both effective core potential and all-electron basis sets. Reasonable agreement between theory and experiment is obtained for the experimentally characterized intermediates, although all-electron basis sets and spin-orbit effects are needed for quantitative agreement. The observed barrier for oxidation is shown to likely correspond to the energy of the crossing between surfaces corresponding to the ground state electronic configuration of Sm⁺ (⁸F,4f⁶6s¹) and an excited surface having two electrons in the valence space (excluding 4f), which are needed to form the strong SmO⁺ bond.

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 Fragmentation of Aluminum Oxide Nitrate Anions Driven by Reactive Oxygen Radical Ligands

Gas-phase metal nitrate anions are known to yield a variety of interesting metal oxides upon fragmentation. The aluminum nitrate anion complexes, Al(NO3)4(-) and AlO(NO3)3(-) were generated by electrospray ionization and studied with collision-induced dissociation and energy-resolved mass spectrometry. Four different decomposition processes were observed, the loss of NO3(-), NO3(•), NO2(•), and O2. The oxygen radical ligand in AlO(NO3)3(-) is highly reactive and drives the formation of AlO(NO3)2(-) upon loss of NO3(•), AlO2(NO3)2(-) upon NO2(•) loss, or Al(NO2)(NO3)2(-) upon abstraction of an oxygen atom from a neighboring nitrate ligand followed by loss of O2. The AlO2(NO3)2(-) fragment also undergoes elimination of O2. The mechanism for O2 elimination requires oxygen atom abstraction from a nitrate ligand in both AlO(NO3)3(-) and AlO2(NO3)2(-), revealing the hidden complexity in the fragmentation of these clusters.

Cationic Noncovalent Interactions: Energetics and Periodic Trends

In this review, noncovalent interactions of ions with neutral molecules are discussed. After defining the scope of the article, which excludes anionic and most protonated systems, methods associated with measuring thermodynamic information for such systems are briefly recounted. An extensive set of tables detailing available thermodynamic information for the noncovalent interactions of metal cations with a host of ligands is provided. Ligands include small molecules (H2, NH3, CO, CS, H2O, CH3CN, and others), organic ligands (O- and N-donors, crown ethers and related molecules, MALDI matrix molecules), π-ligands (alkenes, alkynes, benzene, and substituted benzenes), miscellaneous inorganic ligands, and biological systems (amino acids, peptides, sugars, nucleobases, nucleosides, and nucleotides). Hydration of metalated biological systems is also included along with selected proton-based systems: 18-crown-6 polyether with protonated peptides and base-pairing energies of nucleobases. In all cases, the literature thermochemistry is evaluated and, in many cases, reanchored or adjusted to 0 K bond dissociation energies. Trends in these values are discussed and related to a variety of simple molecular concepts.

Energy-resolved collision-induced dissociation of non-covalent ions: charge- and guest-dependence of decomplexation reaction efficiencies

The influence of the charge state on the CID reactions undergone by host–guest complexes is monitored by mass spectrometry and described by computational chemistry.

Guided ion beam and computational studies of the decomposition of a model thiourea protein cross-linker

The dissociation of protonated methyl-d3 thiourea-4-butyric acid methyl amide (1), a model of thiourea-based protein cross-linking compounds, is examined both experimentally and computationally. Using a guided ion beam tandem mass spectrometer (GIBMS), the threshold collision-induced dissociation (TCID) of [1 + H](+) with Xe is examined as a function of collision energy. Analysis of the kinetic energy-dependent CID cross sections provides the 0 K barriers for four primary and four secondary dissociation pathways, after accounting for competition between channels, sequential dissociations, unimolecular decay rates, internal energy of reactant ions, and multiple ion-neutral collisions. Computations are used to explore the pathways for the various processes and elucidation of their rate-limiting transition states. These results indicate that dissociation is initiated by migration of the excess proton from sulfur to one of three nitrogen atoms in 1, similar to the "mobile proton" model of peptide fragmentation. The computational energies for the rate-limiting transition states are generally in good agreement with the experimentally derived threshold energies, with MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) results being particularly favorable. This good comparison validates the mechanisms explored theoretically and allows identification of the structures of the various product ions and neutrals.

An IMS-IMS threshold method for semi-quantitative determination of activation barriers: Interconversion of proline cis↔trans forms in triply protonated bradykinin

Collisional activation of selected conformations by multidimensional ion mobility spectrometry (IMS-IMS), combined with mass spectrometry (MS), is described as a method to determine semi-quantitative activation energies for interconversion of different structures of the nonapeptide bradykinin (BK, Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe-Arg). This analysis is based on a calibration involving collision-induced dissociation measurements of ions with known dissociation energies (i.e., "thermometer" ions) such as leucine enkephalin, BK, and amino acid-metal cation systems. The energetic barriers between six conformations of [BK+3H]3+ range from 0.23 ±0.01 to 0.55 ±0.03 eV. Prior results indicate that the major peaks in the IMS distributions correspond to specific combinations of cis and trans configurations of the three proline residues in the peptide sequence. The analysis allows us to directly assess pathways for specific transitions. The combination of structural assignments, experimentally determined barrier heights, onset of the quasi-equilibrium region, and dissociation threshold are used to derive a semi-quantitative potential energy surface for main features of [BK+3H]3+.

Probing the stability of insulin oligomers using electrospray ionization ion mobility mass spectrometry

The peptide hormone insulin is central to regulating carbohydrate and fat metabolism in the body by controlling blood sugar levels. Insulin's most active form is the monomer and the extent of insulin oligomerization is related to insulin's activity of controlling blood sugar levels. Electrospray ionization (ESI) of human insulin produced a series of oligomers from the monomer to the undecamer identified using quadrupole ion mobility time-of-flight mass spectrometry. Previous research suggested that only the monomer, dimer and hexamer are native forms of insulin in solution and the range of oligomers observed in the gas-phase are ESI artifacts. Here the properties of three distinct oligomer bands I, II and III, where both the charge state and number of insulin units of the oligomer increase incrementally, were investigated. When Zn(ii) was added to the insulin sample the same oligomers were observed but with 0-6 Zn(ii) ions bound to each of the oligomers. The oligomers of bands I, II and III were characterized by comparing their drift times, collision cross- sections, relative intensities, collision-induced dissociation (CID) patterns and relative breakdown energies. Insulin oligomers of band I dissociated primarily by releasing either the 2+ or 3+ monomer accompanied by an oligomer that conserved the mass, charge and Zn(ii) of the precursor. Insulin oligomers of bands II and III dissociated primarily by releasing the 2+ monomer accompanied by an oligomer which conserved the mass, charge and Zn(ii) of the precursor. Comparison of CID patterns and breakdown energies showed all the oligomers in band II required higher collision energies to dissociate than the oligomers in band I, and the oligomers of band III required higher energies to dissociate than oligomers of band II. These results show that the amount of excess charge on the oligomer in respect to the number of insulin monomers in the oligomer affects their stability.

Thermodynamics and mechanism of protonated cysteine decomposition: a guided ion beam and computational study

A quantitative molecular description of the decomposition of protonated cysteine, H(+)Cys, is provided by studying the kinetic energy dependence of threshold collision-induced dissociation (CID) with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Primary dissociation channels are deamidation (yielding both NH3 loss and NH4(+) formation) and (H2O + CO) loss reactions, followed by an additional six subsequent decompositions. Analysis of the kinetic energy-dependent CID cross sections provides the 0 K barriers for six different reactions after accounting for unimolecular decay rates, internal energy of reactant ions, multiple ion-molecule collisions, and competition among the decay channels. To identify the mechanisms associated with these reactions, quantum chemical calculations performed at the B3LYP/6-311 + G(d,p) level were used to locate the transition states (TSs) and intermediates for these processes. Single point energies of the reactants, products, and key optimized TSs and intermediates are calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311 + G(2d,2p) basis set. The computational characterization of the elementary steps of these reactions, including the structures of the final products, is validated by quantitative agreement with the experimental energetics. In agreement with previous work, deamidation is facilitated by anchimeric assistance of the thio group, which also leads to an interesting rearrangement of the intact amino acid identified computationally.

Gas-phase perspective on the thermodynamics and kinetics of heterogeneous catalysis

Gas-phase studies of small transition metal cluster cations provide thermochemistry of utility to surface science and heterogeneous catalysis.