Threshold Collision-induced Dissociation Research Articles

Absolute Trends and Accurate and Precise Gas-Phase Binding Energies of 1-Alkyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid Clusters from Combined Independent and Competitive TCID Measurements.

Ionic liquid (IL) development efforts have focused on achieving desired properties via tuning of the IL through variation of the cations and anions. However, works geared toward a microscopic understanding of the nature and strength of the intrinsic cation-anion interactions of ILs have been rather limited such that the intrinsic strength of the cation-anion interactions in ILs is largely unknown. In previous work, we employed threshold collision-induced dissociation approaches supported and enhanced by electronic structure calculations to characterize the nature of the cation-anion interactions in and determine the bond dissociation energies (BDEs) of a series of four 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate anions, [2Cnmim:BF4]+. The cation was varied over the series: 1-ethyl-3-methylimidazolium, [C2mim]+, 1-butyl-3-methylimidazolium, [C4mim]+, 1-hexyl-3-methylimidazolium, [C6mim]+, and 1-octyl-3-methylimidazolium, [C8mim]+, to determine the structural and energetic effects of the size of the 1-alkyl substituent on the binding. The variation in the strength of binding determined for these [2Cnmim:BF4]+ clusters was found to be similar in magnitude to the average experimental uncertainty in these determinations. To definitively establish an absolute order of binding among these [2Cnmim:BF4]+ clusters, we extend this work here to include competitive binding studies of three mixed 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate anions, [Cn-2mim:BF4:Cnmim]+ for n = 4, 6, and 8. Importantly, the results of the present work simultaneously provide the absolute BDEs of these mixed [Cn-2mim:BF4:Cnmim]+ clusters and the absolute relative order of the intrinsic binding interactions as a function of the cation with significantly improved precision. Further, by combining the thermochemical results of the previous and present studies, the BDEs of the [2Cnmim:BF4]+ clusters are more accurately and precisely determined.

Gas-Phase Binding Energies and Dissociation Dynamics of 1-Alkyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid Clusters.

Ionic liquids (ILs) have become increasingly popular due to their useful and unique properties, yet there are still many unanswered questions regarding their fundamental interactions. In particular, details regarding the nature and strength of the intrinsic cation-anion interactions and how they influence the macroscopic properties of ILs are still largely unknown. Elucidating the molecular-level details of these interactions is essential to the development of better models for describing ILs and enabling the purposeful design of ILs with properties tailored for specific applications. Current uses of ILs are widespread and diverse and include applications for energy storage, electrochemistry, designer/green solvents, separations, and space propulsion. To advance the understanding of the energetics, conformations, and dynamics of gas-phase IL clustering relevant to space propulsion, threshold collision-induced dissociation approaches are used to measure the bond dissociation energies (BDEs) of the 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate, [2Cnmim:BF4]+. The cation, [Cnmim]+, is varied across the series, 1-ethyl-3-methylimidazolium [C2mim]+, 1-butyl-3-methylimidazolium [C4mim]+, 1-hexyl-3-methylimidazolium [C6mim]+, and 1-octyl-3-methylimidazolium [C8mim]+, to examine the structural and energetic effects of the size of the 1-alkyl substituent on binding. Complementary electronic structure calculations are performed to determine the structures and energetics of the [Cnmim]+ and [BF4]- ions and their binding preferences in the (Cnmim:BF4) ion pairs and [2Cnmim:BF4]+ clusters. Several levels of theory, B3LYP, B3LYP-GD3BJ, and M06-2X, using the 6-311+G(d,p) basis set for geometry optimizations and frequency analyses and the 6-311+G(2d,2p) basis set for energetics, are benchmarked to examine their abilities to properly describe the nature of the binding interactions and to reproduce the measured BDEs. The modest structural variation among these [Cnmim]+ cations produces only minor structural changes and variation in the measured BDEs of the [2Cnmim:BF4]+ clusters. Present findings indicate that the dominant cation-anion interactions involve the 3-methylimidazolium moieties and that these clusters are sufficiently small that differences in packing effects associated with the variable length of the 1-alkyl substituents are not yet significant.

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 Hydrated Thorium(IV) Trihydroxide Cation: Experimental and Theoretical Investigation of the Binding Energies for Th(OH)3+(H2O)n Complexes (n = 1-4).

Threshold collision-induced dissociation of Th(OH)3+(H2O)n (n = 1-4) with xenon was performed using a guided ion beam tandem mass spectrometer. The primary dissociation pathway for all complexes is a loss of a single water molecule followed by the sequential loss of additional water molecules at higher collision energies. The data were analyzed using a statistical model after accounting for lifetime effects and reactant internal and kinetic energy distributions to obtain 0 K bond dissociation energies (BDEs). These were also converted using rigid rotor/harmonic oscillator approximations to yield thermodynamic values at room temperature. The 0 K BDEs of H2O ligands to Th(OH)3+ (IV) are experimentally determined for the first time as 106 ± 6, 89 ± 6, 76 ± 4, and 51 ± 4 kJ/mol for the first, second, third, and fourth water ligand added. These values agree reasonably well with values calculated at the B3LYP, B3PW91, and PBE0 levels of theory with aug-cc-pVQZ basis sets, whereas B3LYP-GD3BJ, MP2, and CCSD(T) single point energies with (without) counterpoise corrections systematically overestimate the bond energies by about 15 (20), 19 (25), and (18) kJ/mol, respectively.

Ion spectroscopy and guided ion beam studies of protonated asparaginyl-threonine decomposition: Influence of a hydroxyl containing C-Terminal residue on deamidation processes

The spontaneous deamidation of asparagine residues plays a significant role in various biological functions and degenerative, aging diseases. Here, we present a full description of the deamidation reaction (as well as other key dissociation processes) of protonated asparaginyl-threonine, [AsnThr + H]+, via complementary infrared multiple photon dissociation (IRMPD) ion spectroscopy, threshold collision-induced dissociation (TCID), and theoretical studies. IRMPD spectra allow for the clear identification of precursor and product ion structural conformations when compared to theoretically calculated spectra for likely structures. Analysis of kinetic energy dependent cross sections measured via TCID with xenon using a guided ion beam tandem mass spectrometer allows for characterization of the energies involved in the decomposition processes of interest. Threshold energies are compared to relative single point energies of major reaction species calculated at B3LYP, B3LYP-GD3BJ, B3P86, and MP2(full) levels of theory, thereby determining important rate-limiting steps involved in [AsnThr + H]+ decomposition. Our studies confirm the formation of a succinimide intermediate via deamidation of [AsnThr + H]+, an observation consistent with condensed-phase deamidation analyses. Interestingly, our spectroscopic results suggest deamidation does not produce furanone isomers even though theoretical results indicate this pathway (exhibited in gas-phase analyses of similar dipeptide systems) should be competitive. Dehydration of [AsnThr + H]+ is also observed, where theory suggests that oxazolone and oxazoline formation are competitive at threshold energies, but IRMPD analyses conclusively confirm the formation of the oxazoline structure. The comprehensive results presented (in addition to complementary studies discussed herein) allow for a valuable analysis of C-terminal residue side-chain effects on the deamidation process.

Thermodynamics and Reaction Mechanisms for Decomposition of a Simple Protonated Tripeptide, H+GAG: a Guided Ion Beam and Computational Study.

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated glycylalanylglycine (H+GAG) with Xe using a guided ion beam tandem mass spectrometer. Kinetic energy dependent cross sections for nine ionic products were observed and analyzed to provide 0K barriers for the six primary products: [b2]+, [y1 + 2H]+, [b3]+, CO loss, [y2 + 2H]+, and [a1]+; and three secondary products: [a2]+, [a3]+, and CH3CHNH2+, after accounting for multiple ion-molecule collisions, internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-GD3BJ/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the six primary and one secondary products (where the other two secondary products have mechanisms previously established). Geometry optimizations and single-point energy calculations were performed at several levels of theory. These theoretical energies are compared with experimental threshold energies and are found to give reasonably good agreement, with B3LYP-GD3BJ and M06-2X levels of theory performing better than other levels. The results obtained here are also compared with previous results for decomposition of H+GGG. The primary difference observed is a lowering of the threshold for the [b2]+ product ion and a concomitant suppression of the directly competing [y1 + 2H]+ product, the result of specific methylation of the [b2]+ product ion.

Experimental and Computational Study of the Group 1 Metal Cation Chelates with Lysine: Bond Dissociation Energies, Structures, and Structural Trends.

The kinetic energy dependence of the collision-induced dissociation (CID) of Group 1 metal cations (M+ = Li+, Na+, K+, Rb+, and Cs+) chelated to the amino acid lysine (Lys) was measured by threshold CID using a guided ion beam tandem mass spectrometer. The simple loss of neutral lysine is the only dissociation channel observed with the heavier alkali metal cations, whereas CID of Li+(Lys) yields other competing channels including loss of NH3 (the dominant channel at low energy) and eight other reactions. Analysis of the kinetic energy-dependent cross sections yields experimental M+(Lys) bond dissociation energies (BDEs) of 376 ± 30, 219 ± 13, 160 ± 10, 141 ± 6, and 128 ± 4 kJ/mol for Li+, Na+, K+, Rb+, and Cs+, respectively. Computational searches yielded 18 distinct, low-energy structural families related to sites of M+ binding in M+(Lys) complexes and 10 distinct, low-energy structural families for neutral lysine. Among the four levels of theory and three basis sets used, four different ground conformers of M+(Lys) and four different ground conformers of lysine were found, including a ground conformer of K+(Lys) and Cs+(Lys), [Nε,CO(OH)], and its higher energy zwitterionic analogue, [Nε,CO2-], that better explains recent infrared multiple photon dissociation action spectroscopy results. Computational results for predicted ground structures of M+(Lys) complexes yielded computed BDEs in reasonable agreement with experiment.

Conformational Effects on Gas-Phase Acidities of Isomeric C3 and C5 Alkanols.

The competitive threshold collision-induced dissociation technique is used to examine conformational effects on the relative gas-phase acidities of selected alcohols. By use of HF and H2O as reference acids in a local thermochemical network to obtain absolute acidities, the measured 0 K gas-phase acidities for the propanol and pentanol isomers are Δacid H0(CH3CH2CH2O-H) = 1563.9 ± 2.9 kJ/mol, Δacid H0((CH3)2CHO-H) = 1568.2 ± 2.7 kJ/mol, Δacid H0(CH3(CH2)4O-H) = 1556.4 ± 2.9 kJ/mol, and Δacid H0((CH3)2CHCH2CH2O-H) = 1556.5 ± 3.0 kJ/mol. Conformational stabilization during deprotonation results in the observed acidity differences between isomers, which can be compared with the "intrinsic" acidity strength defined as deprotonation of the extended all- anti staggered conformations without relaxation. The intrinsic acidities for the propanol and pentanol isomers are 1567 and 1562 kJ/mol, respectively. The difference in intrinsic and observed acidity is largely due to the result of a twisted geometry of the alkoxide ion, stabilized by electrostatic interaction between the electronegative terminal O atom and a H atom on the γ-carbon. These interactions are primarily due to internal rotation about the Cα-Cβ bonds for n-propoxide and the primary pentoxides.

Protonated Asparaginyl-Alanine Decomposition: a TCID, SORI-CID, and Computational Analysis.

Deamidation of asparagine residues, one of the fastest known post-translational modifications in proteins, plays a significant role in various biological functions and degenerative, aging diseases. Here, we present a full description of deamidation (as well as other key dissociation processes) from protonated asparaginyl-alanine, H+(AsnAla), by studying its kinetic energy-dependent threshold collision-induced dissociation (TCID) with Xe using a guided ion beam tandem mass spectrometer. Relative thresholds compare favorably with those acquired by sustained off-resonance irradiation-CID of H+(AsnAla) with Ar in a Fourier transform ion cyclotron resonance mass spectrometer. Absolute threshold energies from the TCID studies are compared to relative single point energies of major reaction species calculated at the B3LYP, B3LYP-GD3BJ, B3P86, MP2(full), and M06-2X levels of theory. Relative energies of key TSs and products allow for the characterization of the important rate-limiting steps involved in H+(AsnAla) decomposition. The influence of water solvation on key TSs is also explored computationally, where bridging the gap between gas-phase and solvated studies is an important aspect of the biological relevance of this analysis. The comprehensive results presented (in addition to complementary studies discussed herein) allow for an insightful comparison to previous deamidation studies such that effects of the C-terminal residue side chain can be elucidated. Graphical abstract ᅟ.

Threshold collision-induced dissociation and theoretical study of protonated azobenzene.

Protonated azobenzene (AB), H+(C6H5N2C6H5), has been studied using threshold collision-induced dissociation in a guided ion beam tandem mass spectrometer. Product channels observed are C6H5N2+ + C6H6 and C6H5+ + N2 + C6H6. The experimental kinetic energy-dependent cross sections were analyzed using a statistical model that accounts for internal and kinetic energy distributions of the reactants, multiple collisions, and kinetic shifts. From this analysis, the activation energy barrier height of 2.02 ± 0.11 eV for benzene loss is measured. To identify the transition states (TSs) and intermediates (IMs) for these dissociations, relaxed potential energy surface (PES) scans were performed at the B3LYP/aug-cc-pVTZ level of theory. The PES indicates that there is a substantial activation energy along the dissociation reaction coordinate that is the rate-limiting step for benzene loss and at some levels of theory, for subsequent N2 loss as well. Relative energies of the reactant, TSs, IMs, and products were calculated at B3LYP, wB97XD, M06, PBEPBE, and MP2(full) levels of theory using both 6-311++G(2d,2p) and aug-cc-pVTZ basis sets. Comparison of the experimental results with theoretical values from various computational methods indicates how well these theoretical methods can predict thermochemical properties. In addition to these density functional theory and MP2 methods, several high accuracy multi-level calculations such as CBS-QB3, G3, G3MP2, G3B3MP2, G4, and G4MP2 were performed to determine the thermochemical properties of AB including the proton affinity and gas-phase basicity, and to compare the performance of different theoretical methods.

Non-adiabatic behavior in the homolytic and heterolytic bond dissociation of protonated hydrazine: A guided ion beam and theoretical investigation.

Threshold collision-induced dissociation using a guided ion beam tandem mass spectrometer was performed on protonated hydrazine and its perdeuterated variant. The dominant dissociation pathways observed were endothermic homolytic and heterolytic cleavages of the N-N bond. The data were analyzed using a statistical model after accounting for internal and kinetic energy distributions, multiple collisions, and kinetic shifts to obtain 0 K bond dissociation energies. Comparison with literature thermochemistry demonstrates that both channels behave non-adiabatically. Heterolytic bond cleavage yields NH2+ + NH3 products, but the NH2+ fragment is in the spin-restricted excited 1A1 state and not in the spin-forbidden ground 3B1 state, whereas homolytic bond cleavage leads to dissociation to the NH3+ + NH2 product asymptote with NH2 in its excited 2A1 state rather than the energetically favored 2B1 state. The rationale for the non-adiabatic behavior observed in the homolytic bond cleavage is revealed by detailed theoretical calculations of the relevant potential energy surfaces and the relevant occupied valence molecular orbitals. These calculations suggest that the non-adiabatic behavior results from conservation of the σ and π character of the binding and lone pair electrons on the nitrogen atoms.

Binding energies of hydrated cobalt hydroxide ion complexes: A guided ion beam and theoretical investigation.

The sequential bond energies of CoOH+(H2O)x complexes, where x = 1-4, are measured by threshold collision-induced dissociation using a guided ion beam tandem mass spectrometer. The primary dissociation pathway for all reactants consists of loss of a single water molecule. This is followed by the sequential loss of additional water molecules at higher collision energies for the x = 2-4 complexes, whereas the x = 1 reactant loses the OH ligand competitively with the H2O ligand. The kinetic energy dependent cross sections for dissociation of CoOH+(H2O)x complexes are modeled to obtain 0 and 298 K binding energies. Our experimental results agree well with theoretically determined bond dissociation energies (BDEs) at the B3LYP, B3LYP-GD3BJ, B3P86, and MP2(full) levels of theory with a 6-311+G(2d,2p) basis set using geometries and vibrational frequencies determined at the B3LYP/6-311+G(d,p) level. Thermochemical information for the loss of OH from CoOH+(H2O)x where x = 0-4 is also derived by combining the present experimental HO-Co+(H2O) and water loss BDEs from CoOH+(H2O)x with those for Co+(H2O)y from the literature. These BDEs are also compared to theory with mixed results.

Second–shell thermochemistry for hydration of strontium dications as determined by threshold collision-induced dissociation and computations

This work reports a continuation of the guided ion beam tandem mass spectrometer (GIBMS) study of the threshold collision-induced dissociation of Sr2+(H2O)x. The present study measures 0K hydration energies for Sr2+(H2O)x, where x=7–9, which can be added to those previously measured for x=1–6. The full range of these complexes is also examined theoretically employing B3LYP/SDD/def2-TZVP and B3LYP/DHF/def2-TZVPP geometry optimizations with additional B3LYP, B3P86, M06, and MP2(full) single point energy calculations performed with the larger SDD/def2-TZVPP and DHF/def2-TZVPP basis sets, respectively. The transition between first and second hydration shells is explored from both experimental and theoretical findings. Neither experimental nor theoretical results are completely consistent with a single coordination number for the first hydration shell at 298K, but rather suggest mixed coordination numbers ranging from 6 to 8 water molecules.

How Hot are Your Ions Really? A Threshold Collision-Induced Dissociation Study of Substituted Benzylpyridinium "Thermometer" Ions.

The first absolute experimental bond dissociation energies (BDEs) for the main heterolytic bond cleavages of four benzylpyridinium "thermometer" ions are measured using threshold collision-induced dissociation in a guided ion beam tandem mass spectrometer. In this experiment, substituted benzylpyridinium ions are introduced into the apparatus using an electrospray ionization source, thermalized, and collided with Xe at varied kinetic energies to determine absolute cross-sections for these reactions. Various effects are accounted for, including kinetic shifts, multiple collisions, and internal and kinetic energy distributions. These experimentally measured 0K BDEs are compared with computationally predicted values at the B3LYP-GD3BJ, M06-GD3, and MP2(full) levels of theory with a 6-311+G(2d,2p) basis set using vibrational frequencies and geometries determined at the B3LYP/6-311+G(d,p) level. Additional dissociation pathways are observed for nitrobenzylpyridinium experimentally and investigated using these same levels of theory. Experimental BDEs are also compared against values in the literature at the AM1, HF, B3LYP, B3P86, and CCSD(T) levels of theory. Of the calculated values obtained in this work, the MP2(full) level of theory with counterpoise corrections best reproduces the experimental results, as do the similar literature CCSD(T) values. Lastly, the survival yield method is used to determine the characteristic temperature (Tchar) of the electrospray source prior to the thermalization region and to confirm efficient thermalization. Graphical Abstract ᅟ.

Thermochemical Investigations of Hydrated Nickel Dication Complexes by Threshold Collision-Induced Dissociation and Theory.

The experimental bond energies of Ni2+(H2O)x complexes, where x = 4-11, are determined by threshold collision-induced dissociation using a guided ion beam tandem mass spectrometer with an electrospray ionization source. The electrospray ionization source produces a distribution of Ni2+(H2O)x complexes, where an in-source fragmentation technique is employed to access the x = 4-6 complexes and control the population of excited isomers. The kinetic energy-dependent cross sections are modeled to yield 0 K bond energies for sequential loss of neutral water molecules, which are converted to 298 K binding energies. Analysis of the primary and secondary water losses from the Ni2+(H2O)x reactant ion complexes, x = 4-11, provide accurate thermochemistry for the hydration energies of Ni2+ and yield the first experimental values for x = 4 and 5 binding energies. Speculative thermochemistry for excited isomers of the x = 4-6 complexes is also obtained. Quantum chemical calculations explore the relative energies of possible geometries. Theoretical bond energies for ground structures are used for direct comparison with experimental values. Our experimental results agree well with previously calculated and experimentally obtained binding enthalpies as well as with the more extensive quantum chemical calculations performed here.

Thermodynamics and Reaction Mechanisms of Decomposition of the Simplest Protonated Tripeptide, Triglycine: A Guided Ion Beam and Computational Study.

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated triglycine (H+GGG) with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Kinetic energy-dependent cross-sections for 10 ionic products are observed and analyzed to provide 0 K barriers for six primary products: [b2]+, [y1 + 2H]+, [b3]+, CO loss, [y2 + 2H]+, and [a1]+; three secondary products: [a2]+, [a3]+, and [y2 + 2H - CO]+; and two tertiary products: high energy [y1 + 2H]+ and [a2 - CO]+ after accounting for multiple ion-molecule collisions, internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-D3/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the six primary and one secondary products. Geometry optimizations and single point energy calculations were performed at several levels of theory. These theoretical energies are compared with experimental energies and are found to give reasonably good agreement, in particular for the M06-2X level of theory. This good agreement between experiment and theory validates the reaction mechanisms explored computationally here and elsewhere and allows identification of the product structures formed at threshold energies. The present work presents the first measurement of absolute experimental threshold energies of important sequence ions and non-sequence ions: [y1 + 2H]+, [b3]+, CO loss, [a1]+, and [a3]+, and refines those for [b2]+ and [y2 + 2H]+ previously measured. Graphical Abstract ᅟ.

Threshold collision-induced dissociation of protonated hydrazine and dimethylhydrazine clustered with water.

Threshold collision-induced dissociation using a guided ion beam tandem mass spectrometer is performed on (N2H4)H+(H2O)n, where n = 1 and 2, and on the protonated unsymmetrical 1,1-dimethylhydrazine one-water complex. The primary dissociation pathway for all clusters is a loss of a single water molecule, which for n = 2 is followed by the sequential loss of an additional water molecule at higher collision energies. The data are analyzed using a statistical model after accounting for internal and kinetic energy distributions, multiple collisions, and kinetic shifts to obtain 0 K bond dissociation energies (BDEs). These are also converted using a rigid rotor/harmonic oscillator approximation to yield thermodynamic values at room temperature. Experimental BDEs compare favorably to theoretical BDEs determined at the B3LYP, M06, mPW1PW91, PBE0, MP2(full), and CCSD(T) levels of theory with a 6-311+G(2d,2p) basis set both with and without empirical dispersion. These calculations also allow visualization of the structures of these complexes, which are simple hydrogen-bonded donor-acceptors.

Threshold Collision-Induced Dissociation of Proton-Bound Hydrazine and Dimethylhydrazine Clusters.

Threshold collision-induced dissociation (TCID) using a guided ion beam tandem mass spectrometer is performed on (N2H4)nH+ where n = 2-4 and on the proton-bound unsymmetrical 1,1-dimethylhydrazine (UDMH) dimer complex. The primary dissociation pathway for all reactants consists of loss of a single hydrazine (or UDMH) molecule followed by the sequential loss of additional hydrazine molecules at higher collision energies for n = 3 and 4. The data were analyzed using a statistical model after accounting for internal and kinetic energy distributions, multiple collisions, and kinetic shifts to obtain 0 K bond dissociation energies (BDEs). These are also converted to values at room temperature by using a rigid rotor/harmonic oscillator approximation and theoretical molecular constants. Experimental BDEs are compared to theoretical BDEs determined at the B3LYP, M06, mPW1PW91, PBE0, MP2(full), and CCSD(T) levels of theory with and without empirical dispersion with a 6-311+G(2d,2p) basis set. The structures of all clusters are explored and exhibit extensive hydrogen bonding.

O2 Protonation Controls Threshold Behavior for N-Glycosidic Bond Cleavage of Protonated Cytosine Nucleosides.

IRMPD action spectroscopy studies of protonated 2'-deoxycytidine and cytidine, [dCyd+H](+) and [Cyd+H](+), have established that both N3 and O2 protonated conformers coexist in the gas phase. Threshold collision-induced dissociation (CID) of [dCyd+H](+) and [Cyd+H](+) is investigated here using guided ion beam tandem mass spectrometry techniques to elucidate the mechanisms and energetics for N-glycosidic bond cleavage. N-Glycosidic bond cleavage is observed as the major dissociation pathways resulting in competitive elimination of either protonated or neutral cytosine for both protonated cytosine nucleosides. Electronic structure calculations are performed to map the potential energy surfaces (PESs) for both N-glycosidic bond cleavage pathways observed. The molecular parameters derived from theoretical calculations are employed for thermochemical analysis of the energy-dependent CID data to determine the minimum energies required to cleave the N-glycosidic bond along each pathway. B3LYP and MP2(full) computed activation energies for N-glycosidic bond cleavage associated with elimination of protonated and neutral cytosine, respectively, are compared to measured values to evaluate the efficacy of these theoretical methods in describing the dissociation mechanisms and PESs for N-glycosidic bond cleavage. The 2'-hydroxyl of [Cyd+H](+) is found to enhance the stability of the N-glycosidic bond vs that of [dCyd+H](+). O2 protonation is found to control the threshold energies for N-glycosidic bond cleavage as loss of neutral cytosine from the O2 protonated conformers is found to require ∼25 kJ/mol less energy than the N3 protonated analogues, and the activation energies and reaction enthalpies computed using B3LYP exhibit excellent agreement with the measured thresholds for the O2 protonated conformers.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated 2'-deoxyguanosine and guanosine.

Experimental and theoretical investigations suggest that hydrolysis of N-glycosidic bonds generally involves a concerted SN2 or a stepwise SN1 mechanism. While theoretical investigations have provided estimates for the intrinsic activation energies associated with N-glycosidic bond cleavage reactions, experimental measurements to validate the theoretical studies remain elusive. Here we report experimental investigations for N-glycosidic bond cleavage of the protonated guanine nucleosides, [dGuo+H](+) and [Guo+H](+), using threshold collision-induced dissociation (TCID) techniques. Two major dissociation pathways involving N-glycosidic bond cleavage, resulting in production of protonated guanine or the elimination of neutral guanine are observed in competition for both [dGuo+H](+) and [Guo+H](+). The detailed mechanistic pathways for the N-glycosidic bond cleavage reactions observed are mapped via electronic structure calculations. Excellent agreement between the measured and B3LYP calculated activation energies and reaction enthalpies for N-glycosidic bond cleavage of [dGuo+H](+) and [Guo+H](+) in the gas phase is found indicating that these dissociation pathways involve stepwise E1 mechanisms in analogy to the SN1 mechanisms that occur in the condensed phase. In contrast, MP2 is found to significantly overestimate the activation energies and slightly overestimate the reaction enthalpies. The 2'-hydroxyl substituent is found to stabilize the N-glycosidic bond such that [Guo+H](+) requires ∼25 kJ mol(-1) more than [dGuo+H](+) to activate the glycosidic bond.