N3 and O2 protonated tautomeric conformations of 2'-deoxycytidine and cytidine coexist in the gas phase

OH Radical-Induced Oxidation in Nucleosides and Nucleotides Unraveled by Tandem Mass Spectrometry and Infrared Multiple Photon Dissociation Spectroscopy

OH⋅-induced oxidation products of DNA nucleosides and nucleotides have been structurally characterized by collision-induced dissociation tandem mass spectrometry (CID-MS2 ) and Infrared Multiple Photon Dissociation (IRMPD) spectroscopy. CID-MS2 results have shown that the addition of one oxygen atom occurs on the nucleobase moiety. The gas-phase geometries of +16 mass increment products of 2'-deoxyadenosine (dA(O)H+ ), 2'-deoxyadenosine 5'-monophosphate (dAMP(O)H+ ), 2'-deoxycytidine (dC(O)H+ ), and 2'-deoxycytidine 5'-monophosphate (dCMP(O)H+ ) are extensively investigated by IRMPD spectroscopy and quantum-chemical calculations. We show that a carbonyl group is formed at the C8 position after oxidation of 2'-deoxyadenosine and its monophosphate derivative. For 2'-deoxycytidine and its monophosphate derivative, the oxygen atom is added to the C5 position to form a C-OH group. IRMPD spectroscopy has been employed for the first time to provide direct structural information on oxidative lesions in DNA model systems.

Influence of 2'-Modifications (O-Methylation, Fluorination, and Stereochemical Inversion) on the Base Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of i-Motif Structures

Naturally occurring and chemically engineered modifications are among the most powerful strategies explored for fine-tuning the conformational characteristics and intrinsic stability of nucleic acids topologies. Modifications at the 2'-position of the ribose or 2'-deoxyribose moieties differentiate nucleic acid structures and have a significant impact on their electronic properties and base-pairing interactions. 2'-O-Methylation, a common post-transcriptional modification of tRNA, is directly involved in modulating specific anticodon-codon base-pairing interactions. 2'-Fluorinated and arabino nucleosides possess novel and beneficial medicinal properties and find use as therapeutics for treating viral diseases and cancer. However, the potential to deploy 2'-modified cytidine chemistries for tuning i-motif stability is largely unknown. To address this knowledge gap, the effects of 2'-modifications including O-methylation, fluorination, and stereochemical inversion on the base-pairing interactions of protonated cytidine nucleoside analogue base pairs, the core stabilizing interactions of i-motif structures, are examined using complementary threshold collision-induced dissociation techniques and computational methods. The 2'-modified cytidine nucleoside analogues investigated here include 2'-O-methylcytidine, 2'-fluoro-2'-deoxycytidine, arabinofuranosylcytosine, 2'-fluoro-arabinofuranosylcytosine, and 2',2'-difluoro-2'-deoxycytidine. All five 2'-modifications examined here are found to enhance the base-pairing interactions relative to the canonical DNA and RNA cytidine nucleosides with the greatest enhancements arising from 2'-O-methylation and 2',2'-difluorination, suggesting that these modifications should well be tolerated in the narrow grooves of i-motif conformations.

Structural Determination of Lysine-Linked Cisplatin Complexes via IRMPD Action Spectroscopy: NNs and NO- Binding Modes of Lysine to Platinum Coexist

Despite its success as an anticancer drug, cisplatin suffers from resistance and produces side effects. To overcome these limitations, amino-acid-linked cisplatin analogues have been investigated. Lysine-linked cisplatin, Lysplatin, (Lys)PtCl₂, exhibited outstanding reactivity toward DNA and RNA that differs from that of cisplatin. To gain insight into its differing reactivity, the structure of Lysplatin is examined here using infrared multiple photon dissociation (IRMPD) action spectroscopy. To probe the influence of the local chemical environment on structure, the deprotonated and sodium-cationized Lysplatin complexes are examined. Electronic structure calculations are performed to explore possible modes of binding of Lys to Pt, their relative stabilities, and to predict their infrared spectra. Comparisons of the measured IRMPD and predicted IR spectra elucidate the structures contributing to the experimental spectra. Coexistence of two modes of binding of Lys to Pt is found where Lys binds via the backbone and side-chain amino nitrogen atoms, NN_s, or to the backbone amino and carboxylate oxygen atoms, NO^-. Glycine-linked cisplatin and arginine-linked cisplatin complexes have previously been found to bind only via the NO^- binding mode. Present results suggest that the NN_s binding conformers may be key to the outstanding reactivity of Lysplatin toward DNA and RNA.

Influence of 5-Halogenation on the Base-Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of Synthetic i-Motif Structures for DNA Nanotechnology Applications

DNA nanotechnology has been employed to develop devices based on i-motif structures. The protonated cytosine-cytosine base pairs that stabilize i-motif conformations are favored under slightly acidic conditions. This unique property has enabled development of the first DNA molecular motor driven by pH changes. The ability to alter the stability and pH transition range of such DNA molecular motors is desirable. Understanding how i-motif structures are influenced by modifications, and which modifications enhance stability and/or affect the pH characteristics, are therefore of great interest. Here, the influence of 5-halogenation of the cytosine nucleobases on the base pairing of protonated cytidine nucleoside analogue base pairs is examined using complementary threshold collision-induced dissociation techniques and computational methods. The nucleoside analogues examined here include the 5-halogenated forms of the canonical DNA and RNA cytidine nucleosides. Comparisons among these systems and to the analogous canonical base pairs previously examined enable the influence of 5-halogenation and the 2'-hydroxy substituent on the base pairing to be elucidated. 5-Halogenation of the cytosine nucleobases is found to enhance the strength of base pairing of DNA base pairs and generally weakens the base pairing for RNA base pairs. Trends in the strength of base pairing indicate that both inductive and polarizability effects influence the strength of base pairing. Overall, the present results suggest that 5-halogenation, and in particular, 5-fluorination and 5-iodination, provide effective means of stabilizing DNA i-motif conformations for applications in nanotechnology, whereas only 5-iodination is effective for stabilizing RNA i-motif conformations but the enhancement in stability is less significant.

Structural determination of arginine-linked cisplatin complexes via IRMPD action spectroscopy: arginine binds to platinum via NO- binding mode

Cisplatin, (NH₃)₂PtCl₂, has been known as a successful metal-based anticancer drug for more than half a century. Its analogue, Argplatin, arginine-linked cisplatin, (Arg)PtCl₂, is being investigated because it exhibits reactivity towards DNA and RNA that differs from that of cisplatin. In order to understand the basis for its altered reactivity, the deprotonated and sodium cationized forms of Argplatin, [(Arg-H)PtCl₂]^- and [(Arg)PtCl₂ + Na]⁺, are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy in the IR fingerprint and hydrogen-stretching regions. Complementary electronic structure calculations are performed using density functional theory approaches to characterize the stable structures of these complexes and to predict their infrared spectra. Comparison of the theoretical IR spectra predicted for various stable conformations of these Argplatin complexes to their measured IRMPD spectra enables determination of the binding mode(s) of Arg to the Pt metal center to be identified. Arginine is found to bind to Pt in a bidentate fashion to the backbone amino nitrogen and carboxylate oxygen atoms in both the [(Arg-H)PtCl₂]^- and [(Arg)PtCl₂ + Na]⁺ complexes, the NO^- binding mode. The neutral side chain of Arg also interacts with the Pt center to achieve additional stabilization in the [(Arg-H)PtCl₂]^- complex. In contrast, Na⁺ binds to both chlorido ligands in the [(Arg)PtCl₂ + Na]⁺ complex and the protonated side chain of Arg is stabilized via hydrogen-bonding interactions with the carboxylate moiety. These findings are consistent with condensed-phase results, indicating that the NO^- binding mode of arginine to Pt is preserved in the electrospray ionization process even under variable pH and ionic strength.

Charge-Tagged Nucleosides in the Gas Phase: UV-Vis Action Spectroscopy and Structures of Cytidine Cations, Dications, and Cation Radicals

Cytidine ribonucleosides were furnished at O5' with fixed-charge 6-trimethylammoniumhexan-1-aminecarbonyl tags and studied by UV-vis photodissociation action spectroscopy in the gas phase to probe isolated nucleobase chromophores in their neutral, protonated, and hydrogen-adduct radical forms. The action spectrum of the doubly charged cytidine conjugate showed bands at 310 and 270 nm that were assigned to the N3- and O2-protonated cytosine tautomers formed by electrospray, respectively. In contrast, cytidine conjugates coordinated to dibenzo-18-crown-6-ether (DBCE) in a noncovalent complex were found to strongly favor protonation at N3, forming a single-ion tautomer. This allowed us to form cytidine N3-H radicals by electron transfer dissociation of the complex and study their action spectra. Cytidine radicals showed only very weak absorption in the visible region of the spectrum for dipole-disallowed transitions to the low (A and B) excited states. The main bands were observed at 360, 300, and 250 nm that were assigned with the help of theoretical vibronic spectra obtained by time-dependent density functional theory calculations of multiple (>300) radical vibrational configurations. Collision-induced dissociations of cytidine radicals proceeded by major cleavage of the N1-C1' glycosidic bond leading to loss of cytosine and competitive loss of N3-hydrogen atom. These dissociations were characterized by calculations of transition-state structures and energies using combined Born-Oppenheimer molecular dynamics and DFT calculations. Overall, cytidine radicals were found to be kinetically and thermodynamically more stable than previously reported analogous adenosine and guanosine radicals.

Influence of 5-Methylation and the 2'- and 3'-Hydroxy Substituents on the Base Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of i-Motif Structures

Repetitive nucleic acid sequences, which occur in abundance throughout the mammalian genome, are of enormous research interest due to their potential to adopt fascinating and unusual molecular structures such as the i-motif. In remarkable contrast to the DNA double helix, i-motif conformations are stabilized by protonated cytosine base pairs, (Cyt)H⁺(Cyt), that are centrally located in the core of the i-motif and intercalated vertically in an antiparallel fashion. An in-depth understanding of how modifications influence the stability of i-motif conformations is a prerequisite to understanding their biological functions and the development of effective means of tuning their stability for specific medical and technological applications. Here, the influence of the 2'- and 3'-hydroxy substituents of the sugar moieties and 5-methylation of the cytosine nucleobases on the base-pairing interactions of protonated cytidine nucleoside analogue base pairs, (xCyd)H⁺(xCyd), are examined by complementary threshold collision-induced dissociation techniques and computational methods. The xCyd nucleosides examined include the canonical DNA and RNA cytidine nucleosides, 2'-deoxycytidine (dCyd) and cytidine (Cyd), as well as several modified cytidine nucleoside analogues, 2',3'-dideoxycytidine (ddCyd), 5-methyl-2'-deoxycytidine (m⁵dCyd), and 5-methylcytidine (m⁵Cyd). Comparisons among these model base pairs indicate that the 2'- and 3'-hydroxy substituents of the sugar moieties have very little influence on the strength of the base-pairing interactions, whereas 5-methylation of the cytosine nucleobases is found to enhance the strength of the base-pairing interactions. The increase in stability resulting from 5-methylation is only modest but is more than twice as large for the DNA than RNA protonated cytidine base pair. Overall, present results suggest that canonical DNA i-motif conformations should be more stable than analogous RNA i-motif conformations and that 5-methylation of cytosine residues, a significant epigenetic marker, provides greater stabilization to DNA than RNA i-motif conformations.

Infrared Multiple Photon Dissociation Spectroscopy of Protonated Cyameluric Acid

The present study reports the first structural characterization of protonated cyameluric acid ([CA + H]⁺) in the gas phase, which paves the way for prospective bottom-up research on the condensed-phase chemistry of CA in the protonated form. A number of [CA + H]⁺ keto-enol isomers can a priori be produced as a result of protonation at available N and O positions of precursor neutral CA tautomers, yet ab initio computations predict different reduced [CA + H]⁺ isomer populations dominating the solution and gas phases that are involved in the ion generation process (i.e., electrospray ionization). Infrared multiple photon dissociation spectra were recorded in the 990-1900 and 3300-3650 cm^-1 regions and compared with theoretical [B3LYP/6-311++G(d,p)] IR absorption spectra of several [CA + H]⁺ isomers, providing a satisfactory agreement for the most stable monohydroxy form in the gas phase, [1358a]⁺, yet the contribution of its nearly isoenergetic OH rotamer, [1358b]⁺, cannot be neglected. This is indicative of the occurrence of [CA + H]⁺ isomer interconversion reactions, assisted by protic solvent molecules, during their transfer into the gas phase. The results suggest that available O positions on neutral CA are energetically favored protonation sites in the gas phase.

Structural and Energetic Effects of O2'-Ribose Methylation of Protonated Pyrimidine Nucleosides

The 2'-substituents distinguish DNA from RNA nucleosides. 2'-O-methylation occurs naturally in RNA and plays important roles in biological processes. Such 2'-modifications may alter the hydrogen-bonding interactions of the nucleoside and thus may affect the conformations of the nucleoside in an RNA chain. Structures of the protonated 2'-O-methylated pyrimidine nucleosides were examined by infrared multiple photon dissociation (IRMPD) action spectroscopy, assisted by electronic structure calculations. The glycosidic bond stabilities of the protonated 2'-O-methylated pyrimidine nucleosides, [Nuom+H]⁺, were also examined and compared to their DNA and RNA nucleoside analogues via energy-resolved collision-induced dissociation (ER-CID). The preferred sites of protonation of the 2'-O-methylated pyrimidine nucleosides parallel their canonical DNA and RNA nucleoside analogues, [dNuo+H]⁺ and [Nuo+H]⁺, yet their nucleobase orientation and sugar puckering differ. The glycosidic bond stabilities of the protonated pyrimidine nucleosides follow the order: [dNuo+H]⁺ < [Nuo+H]⁺ < [Nuom+H]⁺. The slightly altered structures help explain the stabilization induced by 2'-O-methylation of the pyrimidine nucleosides.

Impact of Sodium Cationization on Gas-Phase Conformations of DNA and RNA Cytidine Mononucleotides

Gas-phase conformations of the sodium-cationized forms of the 2'-deoxycytidine and cytidine mononucleotides, [pdCyd+Na]⁺ and [pCyd+Na]⁺, are examined by infrared multiple photon dissociation action spectroscopy. Complimentary electronic structure calculations at the B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) level of theory provide candidate conformations and their respective predicted IR spectra for comparison across the IR fingerprint and hydrogen-stretching regions. Comparisons of the predicted IR spectra and the measured infrared multiple photon dissociation action spectra provide insight into the impact of sodium cationization on intrinsic mononucleotide structure. Further, comparison of present results with those reported for the sodium-cationized cytidine nucleoside analogues elucidates the impact of the phosphate moiety on gas-phase structure. Across the neutral, protonated, and sodium-cationized cytidine mononucleotides, a preference for stabilization of the phosphate moiety and nucleobase orientation is observed, although the details of this stabilization differ with the state of cationization. Several low-energy conformations of [pdCyd+Na]⁺ and [pCyd+Na]⁺ involving several different orientations of the phosphate moiety and sugar puckering modes are observed experimentally.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

Structures and Relative Glycosidic Bond Stabilities of Protonated 2'-Fluoro-Substituted Purine Nucleosides

The 2'-substituent is the primary distinguishing feature between DNA and RNA nucleosides. Modifications to this critical position, both naturally occurring and synthetic, can produce biologically valuable nucleoside analogues. The unique properties of fluorine make it particularly interesting and medically useful as a synthetic nucleoside modification. In this work, the effects of 2'-fluoro modification on the protonated gas-phase purine nucleosides are examined using complementary tandem mass spectrometry and computational methods. Direct comparisons are made with previous studies on related nucleosides. Infrared multiple photon dissociation action spectroscopy performed in both the fingerprint and hydrogen-stretching regions allows for the determination of the experimentally populated conformations. The populated conformers of protonated 2'-fluoro-2'-deoxyadenosine, [Adofl+H]⁺, and 2'-fluoro-2'-deoxyguanosine, [Guofl+H]⁺, are highly parallel to their respective canonical DNA and RNA counterparts. Both N3 and N1 protonation sites are accessed by [Adofl+H]⁺, stabilizing syn and anti nucleobase orientations, respectively. N7 protonation and anti nucleobase orientation dominates in [Guofl+H]⁺. Spectroscopically observable intramolecular hydrogen-bonding interactions with fluorine allow more definitive sugar puckering determinations than possible for the canonical systems. [Adofl+H]⁺ adopts C2'-endo sugar puckering, whereas [Guofl+H]⁺ adopts both C2'-endo and C3'-endo sugar puckering. Energy-resolved collision-induced dissociation experiments with survival yield analyses provide relative glycosidic bond stabilities. The N-glycosidic bond stabilities of the protonated 2'-fluoro-substituted purine nucleosides are found to exceed those of their canonical analogues. Further, the N-glycosidic bond stability is found to increase with increasing electronegativity of the 2'-substituent, i.e., H < OH < F. The N-glycosidic bond stability is also greater for the adenine nucleoside analogues than the guanine nucleoside analogues.

Impact of the 2'- and 3'-Sugar Hydroxyl Moieties on Gas-Phase Nucleoside Structure

Modified nucleosides have been an important target for pharmaceutical development for the treatment of cancer, herpes simplex virus, and the human immunodeficiency virus (HIV). Amongst these nucleoside analogues, those based on 2',3'-dideoxyribose sugars are quite common, particularly in anti-HIV applications. The gas-phase structures of several protonated 2',3'-dideoxyribose nucleosides are examined in this work and compared with those of the analogous protonated DNA, RNA, and arabinose nucleosides to elucidate the influence of the 2'- and combined 2',3'-hydroxyl groups on intrinsic structure. Infrared multiple photon dissociation (IRMPD) action spectra are collected for the protonated 2',3'-dideoxy forms of adenosine, guanosine, cytidine, thymidine and uridine, [ddAdo+H]⁺, [ddGuo+H]⁺, [ddCyd+H]⁺, [ddThd+H]⁺, and [ddUrd+H]⁺, in the IR fingerprint and hydrogen-stretching regions. Molecular mechanics conformational searching followed by electronic structure calculations generates low-energy conformers of the protonated 2',3'-dideoxynucleosides and corresponding predicted linear IR spectra to facilitate interpretation of the measured IRMPD action spectra. These experimental IRMPD spectra and theoretical calculations indicate that the absence of the 2'- and 3'-hydroxyls largely preserves the protonation preferences of the canonical forms. The spectra and calculated structures indicate a slight preference for C3'-endo sugar puckering. The presence of the 3'- and further 2'-hydroxyl increases the available intramolecular hydrogen-bonding opportunities and shifts the sugar puckering modes for all nucleosides but the guanosine analogues to a slight preference for C2'-endo over C3'-endo. Graphical Abstract.

Structural characterization of nucleotide 5'-triphosphates by infrared ion spectroscopy and theoretical studies

The molecular family of nucleotide triphosphates (NTPs), with adenosine 5'-triphosphate (ATP) as its best-known member, is of high biochemical importance as their phosphodiester bonds form Nature's main means to store and transport energy. Here, gas-phase IR spectroscopic studies and supporting theoretical studies have been performed on adenosine 5'-triphosphate, cytosine 5'-triphosphate and guanosine 5'-triphosphate to elucidate the intrinsic structural properties of NTPs, focusing on the influence of the nucleobase and the extent of deprotonation. Mass spectrometric studies involving collision induced dissociation showed similar fragmentation channels for the three studied NTPs within a selected charge state. The doubly charged anions exhibit fragmentation similar to the energy-releasing hydrolysis reaction in nature, while the singly charged anions show different dominant fragmentation channels, suggesting that the charge state plays a significant role in the favorability of the hydrolysis reaction. A combination of infrared ion spectroscopy and quantum-chemical computations indicates that the singly charged anions of all NTPs are preferentially deprotonated at their β-phosphates, while the doubly-charged anions are dominantly αβ-deprotonated. The assigned three-dimensional structure differs for ATP and CTP on the one hand and GTP on the other, in the sense that ATP and CTP show no interaction between nucleobase and phosphate tail, while in GTP they are hydrogen bonded. This can be rationalized by considering the structure and geometry of the NTPs where the final three dimensional structure depends on a subtle balance between hydrogen bond strength, flexibility and steric hindrance.

Modified Quadrupole Ion Trap Mass Spectrometer for Infrared Ion Spectroscopy: Application to Protonated Thiated Uridines

Modifications to a Paul-type quadrupole ion trap mass spectrometer providing optical access to the trapped ion cloud as well as hardware and software for coupling to a table-top IR optical parametric oscillator laser (OPO) are detailed. Critical experimental parameters for infrared multiple photon dissociation (IRMPD) on this instrument are characterized. IRMPD action spectra, collected in the hydrogen-stretching region with this instrument, complemented by spectra in the IR fingerprint region acquired at the FELIX facility, are employed to characterize the structures of the protonated forms of 2-thiouridine, [s²Urd+H]⁺, and 4-thiouridine, [s⁴Urd+H]⁺. The measured spectra are compared with predicted linear IR spectra calculated at the B3LYP/6-311+G(d,p) level of theory to determine the conformers populated in the experiments. This comparison indicates that thiation at the 2- or 4-positions shifts the protonation preference between the 2,4-H tautomer and 4-protonation in opposite directions versus canonical uridine, which displays a roughly equal preference for the 2,4-H tautomer and O4 protonation. As found for canonical uridine, protonation leads to a mixture of conformers exhibiting C2'-endo and C3'-endo sugar puckering with an anti nucleobase orientation being populated for both 2- and 4-thiated uridine. Graphical Abstract ᅟ.

Structural and Energetic Effects of O2'-Ribose Methylation of Protonated Purine Nucleosides

The chemical difference between DNA and RNA nucleosides is their 2'-hydrogen versus 2'-hydroxyl substituents. Modification of the ribosyl moiety at the 2'-position and 2'-O-methylation in particular, is common among natural post-transcriptional modifications of RNA. 2'-Modification may alter the electronic properties and hydrogen-bonding characteristics of the nucleoside and thus may lead to enhanced stabilization or malfunction. The structures and relative glycosidic bond stabilities of the protonated forms of the 2'-O-methylated purine nucleosides, 2'-O-methyladenosine (Adom) and 2'-O-methylguanosine (Guom), were examined using two complementary tandem mass spectrometry approaches, infrared multiple photon dissociation action spectroscopy and energy-resolved collision-induced dissociation. Theoretical calculations were also performed to predict the structures and relative stabilities of stable low-energy conformations of the protonated forms of the 2'-O-methylated purine nucleosides and their infrared spectra in the gas phase. Low-energy conformations highly parallel to those found for the protonated forms of the canonical DNA and RNA purine nucleosides are also found for the protonated 2'-O-methylated purine nucleosides. Importantly, the preferred site of protonation, nucleobase orientation, and sugar puckering are preserved among the DNA, RNA, and 2'-O-methylated variants of the protonated purine nucleosides. The 2'-substituent does however influence hydrogen-bond stabilization as the 2'-O-methyl and 2'-hydroxyl substituents enable a hydrogen-bonding interaction between the 2'- and 3'-substituents, whereas a 2'-hydrogen atom does not. Further, 2'-O-methylation reduces the number of stable low-energy hydrogen-bonded conformations possible and importantly inverts the preferred polarity of this interaction versus that of the RNA analogues. Trends in the CID_50% values extracted from survival yield analyses of the 2'-O-methylated and canonical DNA and RNA forms of the protonated purine nucleosides are employed to elucidate their relative glycosidic bond stabilities. The glycosidic bond stability of Adom is found to exceed that of its DNA and RNA analogues. The glycosidic bond stability of Guom is also found to exceed that of its DNA analogue; however, this modification weakens this bond relative to its RNA counterpart. The glycosidic bond stability of the protonated purine nucleosides appears to be correlated with the hydrogen-bond stabilization of the sugar moiety.

Protonation of Nucleobases in Single- and Double-Stranded DNA

Single-stranded model oligodeoxyribonucleotides, each containing a single protonatable base-cytosine, adenine, guanine, or 5-methylcytosine-centrally located in a background of non-protonatable thymine residues, were acid-titrated in aqueous solution, with UV monitoring. The basicity of the central base was shown to depend on the type of the central base and its nearest neighbours and to rise with increasing oligonucleotide length and decreasing ionic strength of the solution. More complex model oligonucleotides, each containing a centrally located 5-methylcytosine base, were comparatively evaluated in single-stranded and double-stranded form, by UV spectroscopy and high-field NMR. The N³ protonation of the 5-methylcytosine moiety in the double-stranded case occurred at much lower pH, at which the duplex was already experiencing general dissociation, than in the single-stranded case. The central guanine:5-methylcytosine base pair remained intact up to this point, possibly due to an unusual alternative protonation on O² of the 5-methylcytosine moiety, already taking place at neutral or weakly basic pH, as indicated by UV spectroscopy, thus suggesting that 5-methylcytosine sites in double-stranded DNA might be protonated to a significant extent under physiological conditions.

Formic acid catalyzed isomerization of protonated cytosine: a lower barrier reaction for tautomer production of potential biological importance

Tautomerism in nucleotide bases is one of the possible mechanisms of DNA mutation. In spite of numerous studies on the structure and energy of protonated cytosine tautomers, little information is available on the process of their intra- and intermolecular tautomerizations. The catalytic ability of H₂O, HCOOH, and the HCOOHH₂O group to facilitate the tautomerism of the Cyt2t⁺ to CytN3⁺ isomer has been studied. It is shown that the activation free energies of tautomerism in the gas phase are 161.17, 58.96, 26.06, and 15.69 kJ mol^-1, respectively, when the reaction is carried out in the absence and presence of H₂O, HCOOH, or the HCOOHH₂O group. The formation of a doubly hydrogen bonded transition state is central to lowering the activation free energy and facilitating the intramolecular hydrogen atom transfer that is required for isomerization. In the aqueous phase, although the solvent effects of water significantly decrease the activation free energy of intramolecular tautomerization, the isomerization of the Cyt2t⁺ to CytN3⁺ isomer remains unfavorable, and the HCOOH and HCOOHH₂O group mediated mechanisms are still more favorable. Meanwhile, conventional transition state theory (CTST) followed by Wigner tunneling correction is then applied to estimate the rate constants. The rate constant with Wigner tunneling correction for direct tautomerization is obviously smaller than that of HCOOH-mediated tautomerization, which is the most plausible mechanism. Finally, another important finding is that the product complex (CytN3⁺HCOOH) is in the rapid tautomeric equilibrium with the reaction complex (Cyt2t⁺HCOOH) (τ_99.9% = 3.84 × 10^-12 s), which is implemented by the mechanism of the concerted synchronous double proton transfer. Its lifetime of the formed CytN3⁺HCOOH complex (τ = 8.33 × 10^-9 s) is almost one order of magnitude larger than the time required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication (several ns), which is further dissociated into the CytN3⁺ and HCOOH monomers. The results of the present study demonstrate the feasibility of acid catalysis for DNA base isomerization reactions that would otherwise be forbidden.

The intrinsic basicity of the phosphate backbone exceeds that of uracil and thymine residues: protonation of the phosphate moiety is preferred over the nucleobase for pdThd and pUrd

The gas-phase conformations of the protonated forms of thymidine-5'-monophosphate and uridine-5'-monophosphate, [pdThd+H]⁺ and [pUrd+H]⁺, are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and electronic structure calculations. The IRMPD action spectra of [pdThd+H]⁺ and [pUrd+H]⁺ are measured over the IR fingerprint and hydrogen-stretching regions using the FELIX free electron laser and an OPO/OPA laser system. Low-energy conformations of [pdThd+H]⁺ and [pUrd+H]⁺ and their relative stabilities are computed at the MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparisons of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers indicate that the dominant conformers of [pdThd+H]⁺ and [pUrd+H]⁺ populated in the experiments are protonated at the phosphate oxo oxygen atom, with a syn nucleobase orientation that is stabilized by strong P[double bond, length as m-dash]OH⁺O2 and P-OHO4' hydrogen-bonding interactions, and C2'-endo sugar puckering. Minor abundance of conformers protonated at the O2 carbonyl of the nucleobase residue may also contribute for [pdThd+H]⁺, but do not appear to be important for [pUrd+H]⁺. Comparisons to previous IRMPD spectroscopy investigations of the protonated forms of thymidine and uridine, [dThd+H]⁺ and [Urd+H]⁺, and the deprotonated forms of pdThd and pUrd, [pdThd-H]^- and [pUrd-H]^-, provide insight into the effects of the phosphate moiety and protonation on the conformational features of the nucleobase and sugar moieties. Most interestingly, the thymine and uracil nucleobases remain in their canonical forms for [pdThd+H]⁺ and [pUrd+H]⁺, unlike [dThd+H]⁺ and [Urd+H]⁺, where protonation occurs on the nucleobases and induces tautomerization of the thymine and uracil residues.

Protonation induces base rotation of purine nucleotides pdGuo and pGuo

Infrared multiple photon dissociation (IRMPD) action spectra of the protonated forms of 2'-deoxyguanosine-5'-monophosphate and guanosine-5'-monophosphate, [pdGuo+H](+) and [pGuo+H](+), are measured over the IR fingerprint and hydrogen-stretching regions using the FELIX free electron laser and an OPO/OPA laser system. Electronic structure calculations are performed to generate low-energy conformations of [pdGuo+H](+) and [pGuo+H](+) and determine their relative stabilities at the B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparative analyses of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers are performed to determine the most favorable site of protonation and the conformers present in the experiments. These comparisons and the computed energetics find that N7 protonation is considerably preferred over O6 and N3, and the N7 protonated ground-state conformers of [pdGuo+H](+) and [pGuo+H](+) are populated in the experiments. The 2'-hydroxyl substituent does not significantly impact the stable low-energy conformers of [pdGuo+H](+)vs. those of [pGuo+H](+). The effect of the 2'-hydroxyl substituent is primarily reflected in the relative intensities of the measured IRMPD bands, as the IRMPD profiles of [pdGuo+H](+) and [pGuo+H](+) are quite similar. Comparisons to previous IRMPD spectroscopy investigations of the protonated forms of the guanine nucleosides, [dGuo+H](+) and [Guo+H](+), and deprotonated forms of the guanine nucleotides, [pdGuo-H](-) and [pGuo-H](-), provide insight into the effects of the phosphate moiety and protonation on the conformational features of the nucleobase and sugar moieties. Protonation is found to induce base rotation of the guanine residue to an anti orientation vs. the syn orientation found for the deprotonated forms of the guanine nucleotides.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

N3 and O2 Protonated Conformers of the Cytosine Mononucleotides Coexist in the Gas Phase

The gas-phase conformations of the protonated forms of the DNA and RNA cytosine mononucleotides, [pdCyd+H]⁺ and [pCyd+H]⁺, are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy over the IR fingerprint and hydrogen-stretching regions complemented by electronic structure calculations. The low-energy conformations of [pdCyd+H]⁺ and [pCyd+H]⁺ and their relative stabilities are computed at the B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparisons of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers allow the conformers present in the experiments to be determined. Similar to that found in previous IRMPD action spectroscopy studies of the protonated forms of the cytosine nucleosides, [dCyd+H]⁺ and [Cyd+H]⁺, both N3 and O2 protonated cytosine mononucleotides exhibiting an anti orientation of cytosine are found to coexist in the experimental population. The 2'-hydroxyl substituent does not significantly influence the most stable conformations of [pCyd+H]⁺ versus those of [pdCyd+H]⁺, as the IRMPD spectral profiles of [pdCyd+H]⁺ and [pCyd+H]⁺ are similar. However, the presence of the 2'-hydroxyl substituent does influence the relative intensities of the measured IRMPD bands. Comparisons to IRMPD spectroscopy studies of the deprotonated forms of the cytosine mononucleotides, [pdCyd-H]^- and [pCyd-H]^-, provide insight into the effects of protonation versus deprotonation on the conformational features of the nucleobase and sugar moieties. Likewise, comparisons to results of IRMPD spectroscopy studies of the protonated cytosine nucleosides provide insight into the influence of the phosphate moiety on structure. Comparison with previous ion mobility results shows the superiority of IRMPD spectroscopy for distinguishing various protonation sites. Graphical Abstract ᅟ.

Effects of sodium cationization versus protonation on the conformations and N-glycosidic bond stabilities of sodium cationized Urd and dUrd: solution conformation of [Urd+Na]+ is preserved upon ESI

Uridine (Urd) is one of the naturally occurring pyrimidine nucleosides of RNA. 2'-Deoxyuridine (dUrd) is a naturally occurring modified form of Urd, but is not one of the canonical DNA nucleosides. In order to understand the effects of sodium cationization on the conformations and energetics of Urd and dUrd, infrared multiple photon dissociation (IRMPD) action spectroscopy experiments and density functional theory (DFT) calculations are performed. By comparing the calculated IR spectra of [Urd+Na]⁺ and [dUrd+Na]⁺ with the measured IRMPD spectra, the stable low-energy conformers populated in the experiments are determined. Anti oriented bidentate O2 and O2' binding conformers of [Urd+Na]⁺ are the dominant conformers populated in the experiments, whereas syn oriented tridentate O2, O4', and O5' binding conformers of [dUrd+Na]⁺ are dominantly populated in the experiments. The 2'-hydroxyl substituent of Urd stabilizes the anti oriented O2 binding conformers of [Urd+Na]⁺. Significant differences between the measured IRMPD and calculated IR spectra for complexes of [Urd+Na]⁺ and [dUrd+Na]⁺ involving minor tautomeric forms of the nucleobase make it obvious that none are populated in the experiments. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized Urd and dUrd follow the order: [dUrd+H]⁺ < [Urd+H]⁺ < [dUrd+Na]⁺ < [Urd+Na]⁺. The 2'-deoxy modification is found to weaken the glycosidic bond of dUrd versus that of Urd for the sodium cationized uridine nucleosides.

Gas-Phase Conformations and N-Glycosidic Bond Stabilities of Sodium Cationized 2'-Deoxyguanosine and Guanosine: Sodium Cations Preferentially Bind to the Guanine Residue

2'-Deoxyguanosine (dGuo) and guanosine (Guo) are fundamental building blocks of DNA and RNA nucleic acids. In order to understand the effects of sodium cationization on the gas-phase conformations and stabilities of dGuo and Guo, infrared multiple photon dissociation (IRMPD) action spectroscopy experiments and complementary electronic structure calculations are performed. The measured IRMPD spectra of [dGuo+Na]⁺ and [Guo+Na]⁺ are compared to calculated IR spectra predicted for the stable low-energy structures computed for these species to determine the most favorable sodium cation binding sites, identify the structures populated in the experiments, and elucidate the influence of the 2'-hydroxyl substituent on the structures and IRMPD spectral features. These results are compared with those from a previous IRMPD study of the protonated guanine nucleosides to elucidate the differences between sodium cationization and protonation on structure. Energy-resolved collision-induced dissociation (ER-CID) experiments and survival yield analyses of protonated and sodium cationized dGuo and Guo are performed to compare the effects of these cations toward activating the N-glycosidic bonds of these nucleosides. For both [dGuo+Na]⁺ and [Guo+Na]⁺, the gas-phase structures populated in the experiments are found to involve bidentate binding of the sodium cation to the O6 and N7 atoms of guanine, forming a 5-membered chelation ring, with guanine found in both anti and syn orientations and C2'-endo (²T₃ or ³T₂) puckering of the sugar. The ER-CID results, IRMPD yields and the computed C1'-N9 bond lengths indicate that sodium cationization activates the N-glycosidic bond less effectively than protonation for both dGuo and Guo. The 2'-hydroxyl substituent of Guo is found to impact the preferred structures very little except that it enables a 2'OH···3'OH hydrogen bond to be formed, and stabilizes the N-glycosidic bond relative to that of dGuo in both the sodium cationized and protonated complexes.

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.

Enhanced Basicity of Push-Pull Nitrogen Bases in the Gas Phase

Nitrogen bases containing one or more pushing amino-group(s) directly linked to a pulling cyano, imino, or phosphoimino group, as well as those in which the pushing and pulling moieties are separated by a conjugated spacer (C═X)_n, where X is CH or N, display an exceptionally strong basicity. The n-π conjugation between the pushing and pulling groups in such systems lowers the basicity of the pushing amino-group(s) and increases the basicity of the pulling cyano, imino, or phosphoimino group. In the gas phase, most of the so-called push-pull nitrogen bases exhibit a very high basicity. This paper presents an analysis of the exceptional gas-phase basicity, mostly in terms of experimental data, in relation with structure and conjugation of various subfamilies of push-pull nitrogen bases: nitriles, azoles, azines, amidines, guanidines, vinamidines, biguanides, and phosphazenes. The strong basicity of biomolecules containing a push-pull nitrogen substructure, such as bioamines, amino acids, and peptides containing push-pull side chains, nucleobases, and their nucleosides and nucleotides, is also analyzed. Progress and perspectives of experimental determinations of GBs and PAs of highly basic compounds, termed as "superbases", are presented and benchmarked on the basis of theoretical calculations on existing or hypothetical molecules.

Influence of Sodium Cationization versus Protonation on the Gas-Phase Conformations and Glycosidic Bond Stabilities of 2'-Deoxyadenosine and Adenosine

The influence of noncovalent interactions with a sodium cation on the gas-phase structures and N-glycosidic bond stabilities of 2'-deoxyadenosine (dAdo) and adenosine (Ado), [dAdo+Na](+) and [Ado+Na](+), are probed via infrared multiple photon dissociation (IRMPD) action spectroscopy and energy-resolved collision-induced dissociation (ER-CID) experiments. ER-CID experiments are also performed on the protonated forms of these nucleosides, [dAdo+H](+) and [Ado+H](+), for comparison purposes. Complementary electronic structure calculations are performed to determine the structures and relative stabilities of the stable low-energy conformations of the sodium cationized nucleoside complexes and to predict their IR spectra. Comparison between the measured IRMPD action spectra and calculated IR spectra enables the conformations of the sodium cationized nucleosides present in the experiments to be elucidated. The influence of sodium cationization versus protonation on the structures and IR spectra is elucidated by comparison to IRMPD and theoretical results previously reported for the protonated forms of these nucleosides. The influence of sodium cationization versus protonation on the glycosidic bond stability of the adenine nucleosides is determined by comparison of the ER-CID behavior of these systems. All structures present in the experiments are found to involve tridentate binding of Na(+) to the N3, O4', and O5' atoms forming favorable 5- and 6-membered chelation rings, which requires that adenine rotate to a syn configuration. This mode of sodium cation binding results in moderate flexibility of the sugar moiety such that the sugar puckering of the conformations present varies between C2'-endo and O4'-endo. Sodium cationization is found to be less effective toward activating the N-glycosidic bond than protonation for both dAdo and Ado. Both the IRMPD yields and ER-CID behavior indicate that the 2'-hydroxyl substituent of Ado stabilizes the N-glycosidic bond relative to that of dAdo.

Tautomerization lowers the activation barriers for N-glycosidic bond cleavage of protonated uridine and 2'-deoxyuridine

The gas-phase conformations of protonated uridine, [Urd+H](+), and its 2'-deoxy form, protonated 2'-deoxyuridine, [dUrd+H](+), have been examined in detail previously by infrared multiple photon dissociation action spectroscopy techniques. Both 2,4-dihydroxy tautomers and O4 protonated conformers of [Urd+H](+) and [dUrd+H](+) were found to coexist in the experiments with the 2,4-dihydroxy tautomers dominating the population. In the present study, the kinetic energy dependence of the collision-induced dissociation behavior of [Urd+H](+) and [dUrd+H](+) are examined using a guided ion beam tandem mass spectrometer to probe the mechanisms and energetics for activated dissociation of these protonated nucleosides. The primary dissociation pathways observed involve N-glycosidic bond cleavage leading to competitive elimination of protonated or neutral uracil. The potential energy surfaces (PESs) for these N-glycosidic bond cleavage pathways are mapped out via electronic structure calculations for the mixture of 2,4-dihydroxy tautomers and O4 protonated conformers of [Urd+H](+) and [dUrd+H](+) populated in the experiments. The calculated activation energies (AEs) and heats of reaction (ΔHrxns) for N-glycosidic bond cleavage at both the B3LYP and MP2(full) levels of theory are compared to the measured values. The agreement between experiment and theory indicates that B3LYP provides better estimates of the energetics of the species along the PESs for N-glycosidic bond cleavage than MP2, and that the 2,4-dihydroxy tautomers, which are stabilized by strong hydrogen-bonding interactions, predominantly influence the observed threshold dissociation behavior of [Urd+H](+) and [dUrd+H](+).

Diverse mixtures of 2,4-dihydroxy tautomers and O4 protonated conformers of uridine and 2'-deoxyuridine coexist in the gas phase

The gas-phase conformations of protonated uridine, [Urd+H](+), and its modified form, protonated 2'-deoxyuridine, [dUrd+H](+), generated by electrospray ionization are investigated using infrared multiple photon dissociation (IRMPD) action spectroscopy techniques. IRMPD action spectra of [Urd+H](+) and [dUrd+H](+) are measured over the IR fingerprint and hydrogen-stretching regions. [Urd+H](+) and [dUrd+H](+) exhibit very similar IRMPD spectral profiles. However, the IRMPD yields of [Urd+H](+) exceed those of [dUrd+H](+) in both the IR fingerprint and hydrogen-stretching regions. The measured spectra are compared to the linear IR spectra predicted for the stable low-energy structures of these species computed at the B3LYP/6-311+G(d,p) level of theory to determine the tautomeric conformations populated by electrospray ionization. Both B3LYP and MP2 methods find O4 and O2 protonated canonical as well as 2,4-dihydroxy tautomers among the stable low-energy structures of [Urd+H](+) and [dUrd+H](+). Comparison between the measured IRMPD and calculated linear IR spectra suggests that these species exist in their ring-closed forms and that both 2,4-dihydroxy tautomers as well as O4 protonated canonical conformers coexist in the population generated by electrospray ionization for both [Urd+H](+) and [dUrd+H](+). The 2'-deoxy modification of [dUrd+H](+) reduces the variety of 2,4-dihydroxy tautomers populated in the experiments vs. those of [Urd+H](+).

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.

N3 Protonation Induces Base Rotation of 2'-Deoxyadenosine-5'-monophosphate and Adenosine-5'-monophosphate

Infrared multiple photon dissociation (IRMPD) action spectroscopy experiments combined with theoretical calculations are performed to investigate the stable gas-phase conformations of the protonated adenine mononucleotides, [pdAdo+H](+) and [pAdo+H](+). Conformations that are present in the experiments are elucidated via comparative analyses of the experimental IRMPD spectra and the B3LYP/6-311+G(d,p) IR spectra predicted for the conformers optimized at this level of theory. N3 protonation is preferred as it induces base rotation, which allows a strong hydrogen bond to be formed between the excess proton of adenine and the phosphate moiety. In contrast, both N1 and N7 protonation are predicted to be >35 kJ/mol less favorable than N3 protonation. Only N3 protonated conformers are present in the experiments in measurable abundance. Both the low-energy conformers computed and the experimental IRMPD spectra of [pdAdo+H](+) and [pAdo+H](+) indicate that the 2'-hydroxyl moiety does not significantly impact the structure of the most stable conformer or the IRMPD spectral profile of [pAdo+H](+) vs that of [pdAdo+H](+). However, the 2'-hydroxyl leads to a 3-fold enhancement in the IRMPD yield of [pAdo+H](+) in the fingerprint region. Comparison of present results to those reported in a previous IRMPD study of the analogous protonated adenine nucleosides allows the effects of the phosphate moiety on the gas-phase conformations to be elucidated.

2,4-Dihydroxy and O2 Protonated Tautomers of dThd and Thd Coexist in the Gas Phase: Methylation Alters Protonation Preferences versus dUrd and Urd

The gas-phase structures of protonated thymidine, [dThd + H](+), and its modified form, protonated 5-methyluridine, [Thd + H](+), are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy combined with electronic structure calculations. IRMPD action spectra are measured over the ranges extending from ~600 to 1900 cm(-1) and ~2800 to 3800 cm(-1) using the FELIX free electron laser and an optical parametric oscillator/amplifier (OPO/OPA) laser system, respectively. Comparisons between the B3LYP/6-311+G(d,p) linear IR spectra calculated for the stable low-energy conformers and the measured IRMPD spectra are used to determine the most favorable tautomeric conformations of [dThd + H](+) and [Thd + H](+) and to identify those populated in the experiments. Both B3LYP and MP2 levels of theory predict a minor 2,4-dihydroxy tautomer as the ground-state conformer of [dThd + H](+) and [Thd + H](+) indicating that the 2'-hydroxyl substituent of Thd does not exert a significant impact on the structural features. [dThd + H](+) and [Thd + H](+) share parallel IRMPD spectral profiles and yields in both the FELIX and OPO regions. Comparisons between the measured IRMPD and calculated IR spectra suggest that minor 2,4-dihydroxy tautomers and O2 protonated conformers of [dThd + H](+) and [Thd + H](+) are populated in the experiments. Comparison of this work to our previous IRMPD spectroscopy study of protonated 2'-deoxyuridine and uridine suggests that the 5-methyl substituent alters the preferences of O2 versus O4 protonation.

Interaction of Cu(+) with cytosine and formation of i-motif-like C-M(+)-C complexes: alkali versus coinage metals

The Watson-Crick structure of DNA is among the most well-known molecular structures of our time. However, alternative base-pairing motifs are also known to occur, often depending on base sequence, pH, or the presence of cations. Pairing of cytosine (C) bases induced by the sharing of a single proton (C-H(+)-C) may give rise to the so-called i-motif, which occurs primarily in expanded trinucleotide repeats and the telomeric region of DNA, particularly at low pH. At physiological pH, silver cations were recently found to stabilize C dimers in a C-Ag(+)-C structure analogous to the hemiprotonated C-dimer. Here we use infrared ion spectroscopy in combination with density functional theory calculations at the B3LYP/6-311G+(2df,2p) level to show that copper in the 1+ oxidation state induces an analogous formation of C-Cu(+)-C structures. In contrast to protons and these transition metal ions, alkali metal ions induce a different dimer structure, where each ligand coordinates the alkali metal ion in a bidentate fashion in which the N3 and O2 atoms of both cytosine ligands coordinate to the metal ion, sacrificing hydrogen-bonding interactions between the ligands for improved chelation of the metal cation.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated adenine nucleosides: N3 protonation induces base rotation and enhances N-glycosidic bond stability

N3 protonation induces base rotation and stabilizes the syn orientation of the adenine nucleobase of [dAdo+H]⁺ and [Ado+H]⁺via formation of a strong intramolecular N3H⁺⋯O5′ hydrogen-bonding interaction, which in turn influences the mechanisms and energetics for N-glycosidic bond cleavage.

Structure of protonated thymidine characterized by infrared multiple photon dissociation and quantum calculations

RATIONALE

Many fundamental studies are motivated by the probable relationship between the presence of rare enol tautomers of nucleobases and point mutation developing during nucleic acid replication. The evaluation of the tautomeric behaviour of nucleobases is therefore of fundamental importance. This can be probed in the gas phase by combining action spectroscopy and mass spectrometry.

METHODS

Experimental Infrared Multiple Photon Dissociation spectra in the fingerprint region of electrospray-generated and subsequently selected ions were recorded at the CLIO free electron laser (FEL) facility, by coupling the FEL to a quadrupole ion trap, and compared to calculated harmonic vibrational infrared spectra of the different low-lying isomers computed at the B3LYP/6-31++G(d,p) level. Relative energies were refined using the extended basis set 6-311++G(3df,2p).

RESULTS

The Density Functional Theory (DFT) study shows that, as for protonated thymine, the global energy minimum of protonated thymidine corresponds to an enol tautomer, whose infrared absorption spectrum is found to be in very good agreement with the experimental IRMPD spectrum. A very weak IRMPD signal observed at ~1780 cm(-1) is very likely the signature of an oxo tautomer. Consequently, as for thymine, protonated thymidine generated by electrospray corresponds to a mixture of at least two tautomeric forms.

CONCLUSIONS

Tautomerization can be characterized by IRMPD spectroscopy. Interestingly, the dominant enolic tautomeric form(s) presently observed cannot be directly generated from the most stable neutral tautomer of the thymine residue.

Base-Pairing Energies of Protonated Nucleoside Base Pairs of dCyd and m(5)dCyd: Implications for the Stability of DNA i-Motif Conformations

Hypermethylation of cytosine in expanded (CCG)n•(CGG)n trinucleotide repeats results in Fragile X syndrome, the most common cause of inherited mental retardation. The (CCG)n•(CGG)n repeats adopt i-motif conformations that are preferentially stabilized by base-pairing interactions of protonated base pairs of cytosine. Here we investigate the effects of 5-methylation and the sugar moiety on the base-pairing energies (BPEs) of protonated cytosine base pairs by examining protonated nucleoside base pairs of 2'-deoxycytidine (dCyd) and 5-methyl-2'-deoxycytidine (m(5)dCyd) using threshold collision-induced dissociation techniques. 5-Methylation of a single or both cytosine residues leads to very small change in the BPE. However, the accumulated effect may be dramatic in diseased state trinucleotide repeats where many methylated base pairs may be present. The BPEs of the protonated nucleoside base pairs examined here significantly exceed those of Watson-Crick dGuo•dCyd and neutral dCyd•dCyd base pairs, such that these base-pairing interactions provide the major forces responsible for stabilization of DNA i-motif conformations. Compared with isolated protonated nucleobase pairs of cytosine and 1-methylcytosine, the 2'-deoxyribose sugar produces an effect similar to the 1-methyl substituent, and leads to a slight decrease in the BPE. These results suggest that the base-pairing interactions may be slightly weaker in nucleic acids, but that the extended backbone is likely to exert a relatively small effect on the total BPE. The proton affinity (PA) of m(5)dCyd is also determined by competitive analysis of the primary dissociation pathways that occur in parallel for the protonated (m(5)dCyd)H(+)(dCyd) nucleoside base pair and the absolute PA of dCyd previously reported.