Mechanism of the OH Radical Addition to Adenine from Quantum-Chemistry Determinations of Reaction Paths and Spectroscopic Tracking of the Intermediates.

The OH radical is a well-known mediator in the oxidation of biological structures like DNA. Over the past decades, the precise events taking place after reaction of DNA nucleobases with OH radical have been widely investigated by the scientific community. Thirty years after the proposal of the main routes for the reaction of •OH with adenine ( Vieira , A. ; Steenken , S. J. Am. Chem. Soc. 1990 , 112 , 6986 - 6994 ), the present work demonstrates that the OH radical addition to C4 position is a minor pathway. Instead, the dehydration process is mediated by the A5OH adduct. Conclusions are based on density functional theory calculations for the ground-state reactivity and highly accurate multiconfigurational computations for the excited states of the radical intermediates. The methodology has been also used to study the mechanism giving rise to the mutagens 8-oxoA and FAPyA. Taking into account the agreement between the experimental data and the theoretical results, it is concluded that addition to the C5 and C8 positions accounts for at least ∼44.5% of the total •OH reaction in water solution. Finally, the current findings suggest that hydrophobicity in the DNA/RNA surroundings facilitates the formation of 8-oxoA and FAPyA.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2.

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

Excitation of nucleobases from a computational perspective I: reaction paths.

The main intrinsic photochemical events in nucleobases can be described on theoretical grounds within the realm of non-adiabatic computational photochemistry. From a static standpoint, the photochemical reaction path approach (PRPA), through the computation of the respective minimum energy path (MEP), can be regarded as the most suitable strategy in order to explore the electronically excited isolated nucleobases. Unfortunately, the PRPA does not appear widely in the studies reported in the last decade. The main ultrafast decay observed experimentally for the gas-phase excited nucleobases is related to the computed barrierless MEPs from the bright excited state connecting the initial Franck-Condon region and a conical intersection involving the ground state. At the highest level of theory currently available (CASPT2//CASPT2), the lowest excited (1)(ππ*) hypersurface for cytosine has a shallow minimum along the MEP deactivation pathway. In any case, the internal conversion processes in all the natural nucleobases are attained by means of interstate crossings, a self-protection mechanism that prevents the occurrence of photoinduced damage of nucleobases by ultraviolet radiation. Many alternative and secondary paths have been proposed in the literature, which ultimately provide a rich and constructive interplay between experimentally and theoretically oriented research.

Complete-active-space second-order perturbation theory (CASPT2//CASSCF) study of the dissociative electron attachment in canonical DNA nucleobases caused by low-energy electrons (0-3 eV).

Low-energy (0-3 eV) ballistic electrons originated during the irradiation of biological material can interact with DNA/RNA nucleobases yielding transient-anion species which undergo decompositions. Since the discovery that these reactions can eventually lead to strand breaking of the DNA chains, great efforts have been dedicated to their study. The main fragmentation at the 0-3 eV energy range is the ejection of a hydrogen atom from the specific nitrogen positions. In the present study, the methodological approach introduced in a previous work on uracil [I. González-Ramírez et al., J. Chem. Theory Comput. 8, 2769-2776 (2012)] is employed to study the DNA canonical nucleobases fragmentations of N-H bonds induced by low-energy electrons. The approach is based on minimum energy path and linear interpolation of internal coordinates computations along the N-H dissociation channels carried out at the complete-active-space self-consistent field//complete-active-space second-order perturbation theory level. On the basis of the calculated theoretical quantities, new assignations for the adenine and cytosine anion yield curves are provided. In addition, the π1 (-) and π2 (-) states of the pyrimidine nucleobases are expected to produce the temporary anions at electron energies close to 1 and 2 eV, respectively. Finally, the present theoretical results do not allow to discard neither the dipole-bound nor the valence-bound mechanisms in the range of energies explored, suggesting that both possibilities may coexist in the experiments carried out with the isolated nucleobases.

Tuning the photophysical properties of anti-B18H22: efficient intersystem crossing between excited singlet and triplet states in new 4,4'-(HS)2-anti-B18H20.

The tuning of the photophysical properties of the highly fluorescent boron hydride cluster anti-B18H22 (1), by straightforward chemical substitution to produce 4,4'-(HS)2-anti-B18H20 (2), facilitates intersystem crossing from excited singlet states to a triplet manifold. This subsequently enhances O2((1)Δg) singlet oxygen production from a quantum yield of ΦΔ âˆ¼ 0.008 in 1 to 0.59 in 2. This paper describes the synthesis and full structural characterization of the new compound 4,4'-(HS)2-anti-B18H20 (2) and uses UV-vis spectroscopy coupled with density functional theory (DFT) and ab initio computational studies to delineate and explain its photophysical properties.

Communication: Electronic UV-Vis transient spectra of the ∙OH reaction products of uracil, thymine, cytosine, and 5,6-dihydrouracil by using the complete active space self-consistent field second-order perturbation (CASPT2//CASSCF) theory.

Addition of ∙OH radicals to pyrimidine nucleobases is a common reaction in DNA/RNA damage by reactive oxygen species. Among several experimental techniques, transient absorption spectroscopy has been during the last decades used to characterize such compounds. Discrepancies have however appeared in the assignment of the adduct or adducts responsible for the reported transient absorption UV-Vis spectra. In order to get an accurate assignment of the transient spectra and a unified description of the absorption properties of the ∙OH reaction products of pyrimidines, a systematic complete active space self-consistent field second-order perturbation (CASPT2//CASSCF) theory study has been carried out on the uracil, thymine, and cytosine ∙OH addition adducts, as well as on the 5,6-dihydrouracil hydrogen abstraction products. With the obtained findings, the C5OH contributions to the lowest-energy band can be finally discarded. Instead, a bright (2)(π2) state of the C6OH adducts is determined to be the main responsible in all compounds for the absorption band in the Vis range.

Ab initio determination of the ionization potentials of water clusters (H2O)n (n = 2-6).

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic ionization potentials of several water clusters: dimer, trimer, tetramer, pentamer, hexamer book, hexamer ring, hexamer cage, and hexamer prism. The present results establish reference values at a level not reported before for these systems, calibrating different computational strategies and helping to discard less reliable theoretical and experimental data. The systematic study with the increasing size of the water cluster allows obtaining some clues on the structure and reductive properties of liquid water.

On the photophysics and photochemistry of the water dimer.

The photochemistry of the water dimer irradiated by UV light is studied by means of the complete active space perturbation theory//complete active space self-consistent field (CASPT2//CASSCF) method and accurate computational approaches like as minimum energy paths. Both electronic structure computations and ab initio molecular dynamics simulations are carried out. The results obtained show small shifts relative to a single water molecule on the vertical excitation energies of the dimer due to the hydrogen bond placed between the water donor (W(D)) and the water acceptor (W(A)). A red-shift and a blue-shift are predicted for the W(D) and W(A), respectively, supporting previous theoretical and experimental results. The photoinduced chemistry of the water dimer is described as a process occurring between two single water molecules in which the effect of the hydrogen bond plays a minor role. Thus, the photoinduced decay routes correspond to two photodissociation processes, one for each water molecule. The proposed mechanism for the decay channels of the lowest-lying excited states of the system is established as the photochemical production of a hydrogen-bonded H(2)O...HO species plus a hydrogen H atom.

A distance-dependent parameterization of the extended Hubbard model for conjugated and aromatic hydrocarbons derived from stretched ethene.

The Hubbard model, which is widely used in physics but is mostly unfamiliar to chemists, provides an attractive yet simple model for chemistry beyond the self consistent field molecular orbital approximation. The Hubbard model adds an effective electron-electron repulsion when two electrons occupy the same atomic orbital to the familiar Hückel Hamiltonian. Thus it breaks the degeneracy between excited singlet and triplet states and allows an explicit treatment of electron correlation. We show how to evaluate the parameters of the model from high-level ab initio calculations on two-atom fragments and then to transfer the parameters to large molecules and polymers where accurate ab initio calculations are difficult or impossible. The recently developed MS-RASPT2 method is used to generate accurate potential energy curves for ethene as a function of carbon-carbon bond length, which are used to parameterize the model for conjugated hydrocarbons. Test applications to several conjugated/aromatic molecules show that even though the model is very simple, it is capable of reasonably accurate predictions for bond lengths, and predicts molecular excitation energies in reasonable agreement with those from the MS-RASPT2 method.

The role of adenine excimers in the photophysics of oligonucleotides.

Energies and structures of different arrangements of the stacked adenine homodimer have been computed at the ab initio CASPT2 level of theory in isolation and in an aqueous environment. Adenine dimers are shown to form excimer singlet states with different degrees of stacking and interaction. A model for a 2-fold decay dynamics of adenine oligomers can be supported in which, after initial excitation in the middle UV range, unstacked or slightly stacked pairs of nucleobases will relax by an ultrafast internal conversion to the ground state, localizing the excitation in the monomer and through the corresponding conical intersection with the ground state. On the other hand, long-lifetime intrastrand stacked excimer singlet states will be formed in different conformations, including neutral and charge transfer dimers, which originate the red-shifted emission observed in the oligonucleotide chains and that will evolve toward the same monomer decay channel after surmounting an energy barrier. By computing the transient absorption spectra for the different structures considered and their relative stability in vacuo and in water, it is concluded that in the adenine homodimers the maximum-overlap face-to-face orientations are the most stable excimer conformations observed in experiment.

Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching.

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophore's optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Determination of the electron-detachment energies of 2'-deoxyguanosine 5'-monophosphate anion: influence of the conformation.

The vertical electron-detachment energies (VDEs) of the singly charged 2'-deoxyguanosine 5'-monophosphate anion (dGMP-) are determined by using the multiconfigurational second-order perturbation CASPT2 method at the MP2 ground-state equilibrium geometry of relevant conformers. The origin of the unique low-energy band in the gas phase photoelectron spectrum of dGMP-, with maximum at around 5.05 eV, is unambiguously assigned to electron detachment from the highest occupied molecular orbital of pi-character belonging to guanine fragment of a syn conformation. The presence of a short H-bond linking the 2-amino and phosphate groups, the guanine moiety acting as proton donor, is precisely responsible for the pronounced decrease of the computed VDE with respect to that obtained in other conformations. As a whole, the present research supports the nucleobase as the site with the lowest ionization potential in negatively charged (deprotonated) nucleotides at the most stable conformations as well as for B-DNA-like type arrangements, in agreement with experimental evidence.

Photophysics of 1-aminonaphthalene: a theoretical and time-resolved experimental study.

The photophysics of 1-aminonaphthalene (1-napthylamine, AMN) has been investigated on the basis of a constructive experimental-theoretical interplay derived from time-resolved measurements and high-level quantum-chemical ab initio CASPT2//CASSCF calculations. Transient ionization signals at femtosecond resolution were collected for AMN cold isolated molecules following excitation from the vibrationless ground level to a number of vibrational states (within the pump resolution) in the lowest accessible excited state and further multiphoton ionization probing at 500, 800, and 1300 nm. Theory predicts two pipi* states, (1)L(b) and (1)L(a), as the lowest singlet electronic excitations, with adiabatic transitions from S(0) at 3.50 and 3.69 eV, respectively. Since the associated oscillator strength for the lowest transition is exceedingly small, the (1)L(b) state is not expected to become populated significantly and the (1)L(a) state appears as the main protagonist of the AMN photophysics. Though calculations foresee a surface crossing between (1)L(a) and the lower (1)L(b) states, no dynamical signature of it is observed in the time-dependent measurements. In the relaxation of (1)L(a), the radiant emission competes with the intersystem crossing and internal conversion channels. The rates of these mechanisms have been determined at different excitation energies. The internal conversion is mediated by a (1)L(a)/S(0) conical intersection located 0.7 eV above the (1)L(a) minimum. The relaxation of a higher-lying singlet excited state, observed above 40 000 cm(-1) (4.96 eV) and calculated at 5.18 eV, has been also explored.

A three-state model for the photophysics of guanine.

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated 1(pi pi* L(a)) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L(a))CI. On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, 1(pi pi* L(b)) and, in particular, 1(n(O) pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.

Molecular basis of DNA photodimerization: intrinsic production of cyclobutane cytosine dimers.

Based on CASPT2 results, the present contribution establishes for the first time that cytosine photodimer formation (C< >C) is mediated along the triplet and singlet manifold by a singlet-triplet crossing, (T1/S0)X, and by a conical intersection, (S1/S0)CI, respectively. The former can be accessed in a barrierless way from a great variety of photochemical avenues and exhibits a covalent single bond between the ethene C6-C6' carbon atoms of each monomer. The efficiency of the stepwise triplet mechanism, however, would be modulated by the effectiveness of the intersystem crossing mechanism. The results provide the grounds for the understanding of the potential photogenotoxicity of endogenous and exogenous compounds via triplet-triplet sensitization, with a lower bound for cytosine oligonucleotides predicted to be 2.70 eV, and give support to the traditional view of the primary role of triplet excited states in the photochemistry of DNA, a well-known source of photoproducts in solution under triplet photosensitization conditions. The function played by singlet excimers (excited dimers) to explain both the red-shifted fluorescence and photoreaction is highlighted. A rationale on the pronounced wavelength dependence of the observed fluorescence is offered. Geometrical arrangements at the time of light irradiation close to, but energetically above, (S1/S0)CI are suggested as reactive orientations that become prone to produce C< >C directly, with no energy barrier. Because of the outstanding intrinsic ability of cytosine to form stable relaxed excimers, the system located near the bound relaxed excimer has to accumulate enough vibrational energy to surmount a small barrier of 0.2 eV to reach (S1/S0)CI, making the overall process to proceed at a slower relative rate as compared to other compounds such as thymine, which is not susceptible of forming so stable excimers.

Photoreactivity of furocoumarins and DNA in PUVA therapy: formation of psoralen-thymine adducts.

The mechanism of the [2+2] cycloaddition photoreaction of psoralen and a DNA nucleobase, thymine, cornerstone of the furocoumarin-based PUVA (psoralen+UVA radiation) phototherapy, has been studied by the quantum-chemical multiconfigurational CASPT2 method. Triplet- and singlet-mediated mono- and diadduct formations have been determined to take place via singlet-triplet crossings and conical intersections, correlated with the initially promoted triplet or singlet states in different possible reactive orientations. Pyrone-side monoadducts are suggested to be formed in the triplet manifold of the system, and to be less prone to yield diadducts because of the properties of the monoadduct lowest triplet state and the minor accessibility of its excited singlet states. Furan-side monoadducts are better produced in the singlet manifold after reaching a conical intersection with the ground state of the system. From there, the absorption of a second photon would in this case trigger the formation of the diadduct. The proposed mechanisms enable rationalizing the phototherapeutic behavior of several furocoumarins.

Origin of the absorption maxima of the photoactive yellow protein resolved via ab initio multiconfigurational methods.

We discuss the role of the protein in controlling the absorption spectra of photoactive yellow protein (PYP), the archetype xanthopsin photoreceptor, using quantum mechanics/molecular mechanics (QM/MM) methods based on ab initio multireference perturbation theory, combined with molecular dynamics (MD) simulations. It is shown that in order to get results in agreement with the experimental data, it is necessary to use a model that allows for a proper relaxation of the whole system and treats the states involved in the electronic spectrum in a balanced way, avoiding biased results due to the effect of nonrepresentative electrostatic interactions on the chromophore.

Theoretical insight into the intrinsic ultrafast formation of cyclobutane pyrimidine dimers in UV-irradiated DNA: thymine versus cytosine.

The higher formation yields measured in the ultrafast photoinduced formation of cyclobutane thymine dimers (T<>T) with respect to those of cytosine (C<>C) are explained, on the basis of ab initio CASPT2 results, by the existence in thymine of more reactive orientations and a less efficient photoreversibility, whereas in cytosine the funnel toward the photolesion becomes competitive with that mediating the internal conversion of the excited-cytosine monomer.

Ab initio determination of the electron affinities of DNA and RNA nucleobases.

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic electron affinities of the five canonical DNA and RNA nucleobases: uracil, thymine, cytosine, adenine, and guanine. The present results aim for the accurate determination of the intrinsic electron acceptor properties of the isolated nucleic acid bases as described by their electron affinities, establishing an overall set of theoretical reference values at a level not reported before and helping to rule out less reliable theoretical and experimental data and to calibrate theoretical strategies.

On the intrinsic population of the lowest triplet state of thymine.

The population of the lowest triplet state of thymine after near-UV irradiation has been established, on the basis of CASPT2//CASSCF quantum chemical calculations, to take place via three distinct intersystem crossing mechanisms from the initially populated singlet bright 1pipi* state. Two singlet-triplet crossings have been found along the minimum-energy path for ultrafast decay of the singlet state at 4.8 and 4.0 eV, involving the lowest 3npi* and 3pipi* states, respectively. Large spin-orbit coupling elements predict efficient intersystem crossing processes in both cases. Another mechanism involving energy transfer from the lowest 1npi* state with much larger spin-orbit coupling terms can also be proposed. The wavelength dependence measured for the triplet quantum yield of pyrimidine nucleobases is explained by the location and accessibility of the singlet-triplet crossing regions.