REMPI Spectroscopy of Jet-Cooled Guanine

Experimental and theoretical study of the cytosine tautomerism through excited states

CONTEXT

The irradiation of water solution of cytosine with UV light (λmax = 254 nm) shows oxo-hydroxy tautomerism with a rate constant of 6.297 × 10-3 min-1. The order of the reaction implies a tautomeric conversion. After removing the UV light source, we observed a dark reaction with a rate constant of 1.473 × 10-3 min-1 which leads to a restoration of the initial tautomer as before the irradiation. The mechanism of oxo-hydroxy tautomerism of cytosine in water solution was studied in the excited state. It was found that the transformations occur along the 1πσ* excited-state reaction paths which link the Franck-Condon geometries of the tautomers and the conical intersections S0/S1 connected with the H-detachment processes of the corresponding bonds. Furthermore, we established that the conical intersections S0/S1 are also mutually accessible along the 1πσ* excited-state reaction paths.

METHODS

The ground-state equilibrium geometries were optimized at the B3LYP/aug-cc-pVDZ level of theory in water environment according to PCM as well as at the CC2/aug-cc-pVDZ level in the gas phase. The TD B3LYP and CC2 methods were applied for the study of the excited states. The tautomerization mechanisms were studied with the use of the linear interpolation in internal coordinates approach using the optimized geometries of tautomers minima and conical intersections S0/S1 at the CASSCF(6,6)/6-31G* level. All calculations were performed with the GAUSSIAN 16 commercial software.

Ionization energies and ionization-induced structural changes in 2-phenylethylamine and its monohydrate

We report the resonance-enhanced two-photon ionization combined with various detection approaches and quantum chemical calculations of biologically relevant neurotransmitter prototypes, the most stable conformer of 2-phenylethylamine (PEA), and its monohydrate, PEA-H₂O, to reveal the possible interactions between the phenyl ring and amino group in the neutral and ionic species. Extracting the ionization energies (IEs) and appearance energy was achieved by measuring the photoionization and photodissociation efficiency curves of the PEA parent and photofragment ions, together with velocity and kinetic energy-broadened spatial map images of photoelectrons. We obtained coinciding upper bounds for the IEs for PEA and PEA-H₂O of 8.63 ± 0.03 and 8.62 ± 0.04 eV, within the range predicted by quantum calculations. The computed electrostatic potential maps show charge separation, corresponding to a negative charge on phenyl and a positive charge on the ethylamino side chain in the neutral PEA and its monohydrate; in the cations, the charge distributions naturally become positive. The significant changes in geometries upon ionization include switching of the amino group orientation from pyramidal to nearly planar in the monomer but not in the monohydrate, lengthening of the N-H⋯π hydrogen bond (HB) in both species, C_α-C_β bond in the side chain of the PEA⁺ monomer, and the intermolecular O-H⋯N HB in PEA-H₂O cations, leading to distinct exit channels.

New light on the use of Theobroma cacao by Late Classic Maya

In order to address the distribution of and access to cacao, 54 sherds from Late Classic Period Maya residential and civic contexts around El Pilar (Belize/Guatemala) were tested for the presence of cacao. Positive identification of cacao requires that the technique of laser mass spectrometry detect a significant amount of the key biomarker of theophylline to signify cacao. Results show that cacao was culturally significant and widespread and found in civic and residential units regardless of size and location.

Cacao seeds, Theobroma cacao, provide the basis for a ceremonially important Mesoamerican food. Past efforts to identify cacao in ceramics focused on highly decorative vessel forms associated with elite ceremonial contexts, creating assumptions as to how cacao was distributed and who could access it. This study examines 54 archaeological ceramic sherds from El Pilar (Belize/Guatemala) of Late Classic (600 to 900 CE) residential and civic contexts representing a cross-section of ancient Maya inhabitants. Identification of cacao in ancient sherds has depended on the general presence of theobromine; we used the discrete presence of theophylline, a unique key biomarker for cacao in the region. Analysis was done by grinding off all outside surfaces to reduce contamination, pulverizing the inner clay matrix, extracting absorbed molecules, and concentrating the extractions. In order to obtain especially high selectivity and low limits of detection, our study utilized the technique of resonance-enhanced multiphoton ionization coupled with laser-desorption jet-cooling mass spectrometry. This technique isolates molecules in the cold gas phase where they can be selectively ionized through a resonant two-photon process. Of the sherds analyzed, 30 samples (56%) were found to contain significant amounts of theophylline and thus test positive for cacao. Importantly, cacao is present in all contexts, common to all Maya residents near and far from centers.

Guanine Tautomerism in Ionic Complexes with Ag+ Investigated by IRMPD Spectroscopy and Mass Spectrometry

In this paper, we present the IRMPD spectra of three ionic complexes between guanine (G) and silver (Ag⁺): [GAg-H₂O]⁺, [GAgG]⁺ produced in the electrospray ionization source of the mass spectrometer, and [GAg]⁺ produced by collision induced dissociation of the [GAgG]⁺ complex. On the basis of the comparison of theoretically calculated IR spectra, we show that there are two isomers of each complex containing two different keto-amino (KA) tautomers of G (GKA(1,9) and GKA(1,7)). The observed isomers are the most stable structures in aqueous solution, and their experimentally estimated relative populations are in better agreement with the calculated relative populations in solution than in the gas phase, both at 298 K. We concluded that these observations suggest that GKA(1,9) and GKA(1,7) coexist in solution according to previous theoretical reports (Colominas, C.; et al. J. Am. Chem. Soc. 1996, 118, 6811). We were unable to find any evidence of the presence of the GEA(9), GKA(3,7), GKA(3,9), or GKA(7,9), whose relative stabilities in solution are strongly dependent on the theoretical method used to account for the solvent effect (Hanus, M.; et al. J. Am. Chem. Soc. 2003, 125, 7678).

Exploration of the theobromine-water dimer: comparison with DNA microhydration

Understanding the molecular basis of the appearance of life on Earth is an exciting research field. Many factors may have influenced the election of the molecules used by living beings and evolution may have modified those original compounds. In an attempt to understand the role played by intermolecular interactions in the election of CGAT as the alphabet of life, we present here a thorough experimental and computational study on the interaction of theobromine with water. Theobromine is a xanthine derivative, structurally related to the nucleobases, and also present in many living beings. The experimental results demonstrate that the most stable isomer of theobromine-water was formed and detected in supersonic expansions. This isomer very well resembles the structure of the dimers between nucleobases and water, offering similar values of binding energy. A comparison between the results obtained for theobromine-water with those reported in the literature for monohydrates of nucleobases is also offered.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

Characteristics of the excited triplet states of thiolated guanosine derivatives and singlet oxygen generation

Thioguanine is sensitive to UVA light and generates singlet molecular oxygen (1O2*) when exposed to UVA. Three thioguanosine derivatives, 2',3',5'-tri-O-acetyl-6-thioguanosine (ta6TGuo), 2',3',5'-tri-O-acetyl-8-thioguanosine (ta8TGuo), and 2',3',5'-tri-O-acetyl-6,8-dithioguanosine (taDTGuo) were explored photophysically and photochemically. Nanosecond transient absorption and time-resolved near-infrared emission measurements were carried out to investigate the characteristics of their excited triplet states in acetonitrile solution. The quantum yield of intersystem crossing (ΦISC), the intrinsic decay rate constant (k0), the quenching rate constant by 3O2 (kq) and the self-quenching rate constant (kSQ) of their triplet states were all determined. From the precise analysis of the quantum yield of 1O2* generation (ΦΔ) against the concentration of dissolved molecular oxygen, the fraction of the triplet states quenched by dissolved oxygen which gives rise to 1O2* formation (SΔ) was successfully obtained with high accuracy. The ΦΔ values at low oxygen concentrations reveal that these thioguanosines, particularly taDTGuo, can still effectively generate 1O2* at low molecular oxygen concentrations like carcinomatous microenvironments. These findings indicate that taDTGuo would perform well as a potential agent for photo-induced cancer therapies.

How nature covers its bases

The response of nucleobases to UV radiation depends on structure in subtle ways, as revealed by gas-phase experiments.

Excited-State Dynamics in O6-Methylguanosine: Impact of O6-Methylation on the Relaxation Mechanism of Guanine Monomers

Absorption of ultraviolet radiation by DNA bases results in ultrafast internal conversion to the ground state, which minimizes photodamage. However, exogenous and endogenous alkylating agents present in the cellular environment can methylate the nucleobases in DNA. In particular, methylation of guanosine at the O⁶ position in DNA leads to the formation of the O⁶-methylguanosine adduct, which may alter the photostability of DNA. This contribution demonstrates that O⁶-methylation of guanosine red shifts its ground-state absorption spectrum and slows down the rate of internal conversion to the ground state by ∼40-fold in aqueous solution. The 40-fold decrease in the rate of excited-state decay increases the probability of photodamage within cellular DNA. It is proposed that the longer decay lifetime corresponds to relaxation of the excited-state population in O⁶-methylguanosine along a C6-puckered reaction coordinate in the ¹ππ*(L_a) potential energy surface that runs parallel to an ultrafast internal conversion pathway along a C2-puckered coordinate.

Excited State Dynamics of 6-Thioguanine

Here we present the excited state dynamics of jet-cooled 6-thioguanine (6-TG), using resonance-enhanced multiphoton ionization (REMPI), IR-UV double resonance spectroscopy, and pump-probe spectroscopy in the nanosecond and picosecond time domains. We report data on two thiol tautomers, which appear to have different excited state dynamics. These decay to a dark state, possibly a triplet state, with rates depending on tautomer form and on excitation wavelength, with the fastest rate on the order of 10¹⁰ s^-1. We also compare 6-TG with 9-enolguanine, for which we observed decay to a dark state with a 2 orders of magnitude smaller rate. At increased excitation energy (∼+500 cm^-1) an additional pathway appears for the predominant thiol tautomer. Moreover, the excited state dynamics for 6-TG thiols is different from that recently predicted for thiones.

Direct Analysis of Xanthine Stimulants in Archaeological Vessels by Laser Desorption Resonance Enhanced Multiphoton Ionization

Resonance enhanced multiphoton ionization spectroscopy (REMPI) generates simultaneous vibronic spectroscopy and fragment free mass spectrometry to identify molecules within a complex matrix. We combined laser desorption with REMPI spectroscopy to study organic residues within pottery sherds from Maya vessels (600-900 CE) and Mississippian vessels (1100-1200 CE), successfully detecting three molecular markers, caffeine, theobromine, and theophylline, associated with the use of cacao. This analytical approach provides a high molecular specificity, based on both wavelength and mass identification. At the same time, the high detection limit allows for direct laser desorption from sherd scrapings, avoiding the need for extracting organic constituents from the sherd matrix.

Capturing Transient Endoperoxide in the Singlet Oxygen Oxidation of Guanine

The chemistry of singlet O2 toward the guanine base of DNA is highly relevant to DNA lesion, mutation, cell death, and pathological conditions. This oxidative damage is initiated by the formation of a transient endoperoxide through the Diels-Alder cycloaddition of singlet O2 to the guanine imidazole ring. However, no endoperoxide formation was directly detected in native guanine or guanosine, even at -100 °C. Herein, gas-phase ion-molecule scattering mass spectrometry was utilized to capture unstable endoperoxides in the collisions of hydrated guanine ions (protonated or deprotonated) with singlet O2 at ambient temperature. Corroborated by results from potential energy surface exploration, kinetic modeling, and dynamics simulations, various aspects of endoperoxide formation and transformation (including its dependence on guanine ionization and hydration states, as well as on collision energy) were determined. This work has pieced together reaction mechanisms, kinetics, and dynamics data concerning the early stage of singlet O2 induced guanine oxidation, which is missing from conventional condensed-phase studies.

First-principles anharmonic quantum calculations for peptide spectroscopy: VSCF calculations and comparison with experiments

First-principles quantum calculations for anharmonic vibrational spectroscopy of three protected dipeptides are carried out and compared with experimental data.

Gas-Phase IR Spectroscopy of Nucleobases

IR spectroscopy of nucleobases in the gas phase reflects simultaneous advances in both experimental and computational techniques. Important properties, such as excited state dynamics, depend in subtle ways on structure variations, which can be followed by their infrared signatures. Isomer specific spectroscopy is a particularly powerful tool for studying the effects of nucleobase tautomeric form and base pair hydrogen-bonding patterns.

Gas-Phase Conformations and Energetics of Protonated 2'-Deoxyguanosine and Guanosine: IRMPD Action Spectroscopy and Theoretical Studies

The gas-phase structures of protonated 2'-deoxyguanosine, [dGuo+H](+), and its RNA analogue protonated guanosine, [Guo+H](+), are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical electronic structure calculations. IRMPD action spectra are measured over the range extending from ∼550 to 1900 cm(-1) using the FELIX free electron laser and from ∼2800 to 3800 cm(-1) using an optical parametric oscillator/amplifier (OPO/OPA) laser system. The measured IRMPD spectra of [dGuo+H](+) and [Guo+H](+) are compared to each other and to B3LYP/6-311+G(d,p) linear IR spectra predicted for the stable low-energy conformations computed for these species to determine the most favorable site of protonation, identify the structures accessed in the experiments, and elucidate the influence of the 2'-hydroxyl substituent on the structures and the IRMPD spectral features. Theoretical energetics and the measured IRMPD spectra find that N7 protonation is preferred for both [dGuo+H](+) and [Guo+H](+), whereas O6 and N3 protonated conformers are found to be much less stable. The 2'-hydroxyl substituent does not exert a significant influence on the structures and relative stabilities of the stable low-energy conformations of [dGuo+H](+) versus [Guo+H](+) but does provide additional opportunities for hydrogen bonding such that more low-energy structures are found for [Guo+H](+). [dGuo+H](+) and [Guo+H](+) share very parallel IRMPD spectral features in the FELIX and OPO regions, whereas the effect of the 2'-hydroxyl substituent is primarily seen in the relative intensities of the measured IR bands. The measured OPO/OPA spectral signatures, primarily reflecting the IR features associated with the O-H and N-H stretches, provide complementary information to that of the FELIX region and enable the conformers that arise from different protonation sites to be more readily distinguished. Insight gained from this and parallel studies of other DNA and RNA nucleosides and nucleotides should help better elucidate the effects of the local environment on the overall structures of DNA and RNA.

New Method for Double-Resonance Spectroscopy in a Cold Quadrupole Ion Trap and Its Application to UV-UV Hole-Burning Spectroscopy of Protonated Adenine Dimer

A novel method for double-resonance spectroscopy in a cold quadrupole ion trap is presented, which utilizes dipolar resonant excitation of fragment ions in the quadrupole ion trap. Photofragments by a burn laser are removed by applying an auxiliary RF to the trap, and a probe laser detects the depletion of photofragments by the burn laser. By scanning the wavelength of the burn laser, conformation-specific UV spectrum of a cold ion is obtained. This simple and powerful method is applicable to any type of double-resonance spectroscopy in a cold quadrupole ion trap and was applied to UV-UV hole-burning spectroscopy of protonated adenine dimer. It was found that protonated adenine dimer has multiple conformers/tautomers, each with multiple excited states with drastically different excited state dynamics.

The solvent effect and identification of a weakly emissive state in nonradiative dynamics of guanine nucleosides and nucleotides--a combined femtosecond broadband time-resolved fluorescence and transient absorption study

A combined method of femtosecond broadband time-resolved fluorescence (fs-TRF) and transient absorption (fs-TA) was employed to investigate the excited state dynamics of 2'-deoxyguanosine (dG) and 2'-deoxyguanosine 5'-monophosphate (dGMP). Comparative fs-TRF and fs-TA measurements were conducted on dG and dGMP in neutral water, deuterated water, and methanol with excitation wavelengths of 245, 267 and 285 nm. Very similar results were observed with dG and dGMP. The data provide compelling evidence for the co-existence of two nonradiative pathways. One is the generally recognized Laππ* mediated channel, the other involves an unprecedented weakly emissive state which plays a significant role in the overall deactivation processes. The Laππ* channel features biphasic dynamics with time constants (τ1/τ2) of ~0.2/0.8 ps in water and ~0.25/1.0 ps in methanol. The biphasic decay arises due to a partial transfer with τ1 of the Laππ* population to the newly identified state followed by conversion in τ2 of the remaining Laππ* molecules into the electronic ground state. The channel mediated by the weakly emissive species shows solvent-dependent dynamics with time constants (τ3) of ~2.0 ps in water, ~2.3 ps in deuterated water, and ~4.1 ps in methanol. The species features absorption at UV wavelengths (~300-400 nm) and exhibits deeply red-shifted fluorescence (λmax ~ 520 nm) with polarization direction varied markedly from that of the Laππ* but close to the Lbππ*. This species acts as an effective quenching state to the radiative decay of the brightly emissive Laππ* and Lbππ*. It sets in promptly (<~50 fs) after the photoexcitation and is further populated through nonadiabatic coupling with the Laππ*. The overall involvement of this state is facilitated with excitation at high energy and is favoured in methanol over water. According to the spectral character and the solvent effect in particular the kinetic isotope effect, the species is tentatively associated to the πσ* state with charge transfer (CT) character which is considered to be preferentially stabilized by hydrogen-bonding between the guanine amino and surrounding solvent molecules. The result of this study leads to a dramatically different picture of guanine deactivation. It demonstrates a crucial role of the solvent in shaping the nonradiative dynamics of guanine nucleosides and nucleotides. The data presented are important for understanding the detailed photophysics of not only the monomeric guanine but also DNA assemblies that contain guanine in base pairs or have a guanine tetrad as the structural motif.

Coupled-cluster and density functional theory studies of the electronic 0–0 transitions of the DNA bases

The 0–0 transitions of the electronic excitation spectra of the lowest tautomers of the four nucleotide (DNA) bases have been studied using linear-response approximate coupled-cluster singles and doubles (CC2) calculations.

ABSORPTION SPECTRA OF NUCLEIC ACID BASES IN WATER ENVIRONMENT: INSIGHTS INTO FROM COMBINED QM/MM AND CLUSTER-CONTINUUM MODEL CALCULATIONS

Electronic spectra of uracil, thymine, adenine, guanine, and cytosine in the gas phase and aqueous solution have been studied by extensive time-dependent density functional calculations. Calculations show that the Quantum mechanics/molecular mechanics (QM/MM) geometry optimization based on the molecular dynamics (MD) equilibrated configuration can locate an optimal solvated cluster for the base solvation, and the combined QM/MM and cluster-continuum computational protocol is capable of handling the solvent effect on the excited states of nucleic acid bases and providing realistic absorption spectra in water environment with relatively low computational costs. Generally, the vertical excitation energies in aqueous solution by PCM/TD-X3LYP calculations show excellent agreement with the experimental observations and the maximum deviation is less than 0.2 eV. The present results reveal that the hydrogen bond network around the excited-state base and its dipole moment change may remarkably modify the absorption spectra of nucleic acid bases in aqueous solution.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β(μ)) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree-Fock, Møller-Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning's correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β(μ) (9H)/β(μ) (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β(μ) values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

THEORETICAL STUDIES ON PHOTOIONIZATION OF GUANINE TAUTOMERS AND INTERCONVERSION OF CATION RADICALS

The G3MP2B3 and P3 methods have been used to calculate the adiabatic and vertical ionization potentials (IPs) of the eight most stable tautomers of guanine. The calculated energy discrepancy between adiabatic and vertical IPs are in good agreement with the changes in geometry from neutral ground state to stable cation radicals. The geometries of imino-oxo form tautomers have no obvious change in the ionization process, which results in less energy discrepancy between vertical and adiabatic IPs. In the ionization process, the geometries of the amino-oxo and amino-hydroxy form tautomers change from nonplanar to planar structures. Hence the amino-oxo and amino-hydroxy form tautomers have larger energy discrepancy between vertical and adiabatic IPs. Further studies on the interconversion of the cation radicals shed further light on the transition process between the cation radicals and the main pathways are the hydrogen migrations and internal rotations of hydroxy (OH) and imino (NH) groups. The barriers of hydrogen rotations are lower than those of hydrogen migrations. Furthermore, the barriers for the hydrogen migrations between two rings are higher, which are about 3.0 eV.

A new tandem mass spectrometer for photofragment spectroscopy of cold, gas-phase molecular ions

We present here the design of a new tandem mass spectrometer that combines an electrospray ion source with a cryogenically cooled ion trap for spectroscopic studies of cold, gas-phase ions. The ability to generate large ions in the gas phase without fragmentation, cool them to approximately 10 K in an ion trap, and perform photofragment spectroscopy opens up new possibilities for spectroscopic characterization of large biomolecular ions. The incorporation of an ion funnel, together with a number of small enhancements, significantly improves the sensitivity, signal stability, and ease of use compared with the previous instrument built in our laboratory.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Characterizing radiation-induced oxidation of DNA by way of the monohydrated guanine-cytosine radical cation

The interaction of one water molecule with the guanine-cytosine radical cation has been studied with ab initio and density functional methods in order to help elucidate the nature of oxidized aqueous DNA. The theoretical spin density of [GC]*(+) reveals that the radical center is localized on guanine. The adiabatic ionization potential lowers from 7.63 to 6.71 eV in concurrence with the formation of the Watson-Crick base pair and hydration by one water molecule. A natural bond orbital analysis of partial charges shows that approximately 80% of the positive charge persists on guanine upon hydration and formation of the Watson-Crick base pair with cytosine. Hydration energies were computed with second-order Z-averaged perturbation theory (ZAPT2) using the aug-cc-pVDZ basis set at 11 stationary points on the B3LYP/DZP++ potential energy surface. The hydration energy at the global minimum is 14.2 kcal mol(-1). The lowest energy structures correspond to hydration near the glycosidic bond sites. Structural changes in the Watson-Crick base pair are predominantly seen for monohydration in the groove regions of double-helix DNA.

Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.

Tautomeric equilibrium, stability, and hydrogen bonding in 2'-deoxyguanosine monophosphate complexed with Mg2+

The tautomeric equilibrium and hydrogen bonding in nucleotide 2'-deoxyguanosine monophosphate that interacts with hydrated Mg2+ cation (4H2O.Mg[dGMP]) were studied at the MP2/cc-pVDZ//B3LYP/cc-pVDZ and B3LYP/aug-cc-pVTZ//B3LYP/cc-pVDZ levels of theory. The Mg2+ ion forms two inner-shell contacts with the nucleotide, similar to small phosphorylated molecules under physiological conditions. The presence of the phosphate group and the hydrated magnesium cation leads to a change in guanine tautomeric equilibrium of 4H2O.Mg[dGMP] in comparison to free guanine. The influence of the phosphate group and the magnesium cation on tautomeric equilibrium is larger in the anti conformation where the P=O-->Mg and Mg<--N7 coordinate bonds are formed. The canonical oxo form of guanine is more stable (by 6-8 kcal/mol) than the O6-hydroxo form in anti conformation. Thus, the interaction with Mg2+ ion is capable of further suppressing the likelihood of a spontaneous transient formation of the rare tautomer. In the syn conformation of 4H2O.Mg[dGMP], the interaction of the guanine nucleobase with the phosphate group and the magnesium cation is not as strong as in the anti conformation, and the relative stability of guanine tautomers is close to those in free guanine.

Electronic Spectra, Excited State Structures and Interactions of Nucleic Acid Bases and Base Assemblies: A Review

A comprehensive review of recent theoretical and experimental advances in the singlet electronic transitions, excited state structures and dynamics of nucleic acid bases (NABs) and base assemblies are presented. It is well known that NABs absorb ultraviolet radiation, but the absorbed energy is efficiently dissipated in the form of ultrafast internal conversion processes believed to occur in the subpicosecond time scale and, therefore, enabling NABs highly photostable. It is not known how much evolutionary role was played in evolving these molecules and the ultimate selection by nature as genetic materials, but it is well accepted that survival-of-fittest prevails. Recently, significant efforts have been continuously paid to understand the mechanism of electronic excitation deactivation, but universally acceptable mechanism is still elusive. However, recent investigations reveal that electronic excited state geometries of DNA bases are usually nonplanar and this structural nonplanarity may facilitate nonradiative deactivation. Investigation of excited state structures is challenging and, therefore, it is not surprising that despite the impressive theoretical and computational advances, this research area is still hampered by the methodological and computational limitations. Further, stacking has significant influence on the emission properties of molecules. The 2-aminopurine, a fluorescent adenine derivative frequently used in studying DNA dynamics, shows significant attenuations in fluorescence quantum yield when incorporated in the DNA. Theoretical and computational bottlenecks limit a thorough theoretical understanding of effect of stacking interactions on the excited state dynamics of NABs. Despite these limitations the investigations of excited state properties are progressing in the right direction and our better understanding of excited state structure and dynamics of NABs and nucleic acids may help to design preventive strategy for radiation induced illness and photostable materials.

Potential energy surfaces of the microhydrated guanine...cytosine base pair and its methylated analogue

A complete scan of the potential and free-energy surfaces of monohydrated and dihydrated guanine...cytosine and 9-methylguanine...1-methylcytosine base pairs was realized by the molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The most stable and populated structures localized were further fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method with a large basis set. A systematic study of microhydration of these systems using a high-level correlated ab initio approach is presented for the first time. The different behavior of guanine...cytosine and adenine...thymine complexes is also discussed. These studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for better understanding of the influence of the solvent on the stability of the structure of DNA.