Microhydration of cytosine and its radical anion: Cytosine∙(H2O)n (n=1–5)

Spectroscopic observation and ultrafast coherent vibrational dynamics of the aqueous phenylalanine radical

The phenylalanine radical (Phe˙) has been proposed to mediate biological electron transport (ET) and exhibit long-lived electronic coherences following attosecond photoionization. However, the coupling of ultrafast structural reorganization to the oxidation/ionization of biomolecules such as phenylalanine remains unexplored. Moreover, studies of ET involving Phe˙ are hindered by its hitherto unobserved electronic spectrum. Here, we report the spectroscopic observation and coherent vibrational dynamics of aqueous Phe˙, prepared by sub-6 fs photodetachment of phenylalaninate anions. Sub-picosecond transient absorption spectroscopy reveals the ultraviolet absorption signature of Phe˙. Ultrafast structural reorganization drives coherent vibrational motion involving nine fundamental frequencies and one overtone. DFT calculations rationalize the absence of the decarboxylation reaction, a photodegradation pathway previously identified for Phe˙. Our findings guide the interpretation of future attosecond experiments aimed at elucidating coherent electron motion in photoionized aqueous biomolecules and pave way for the spectroscopic identification of Phe˙ in studies of biological ET.

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.

Electron and hole interactions with P, Z, and P:Z and the formation of mutagenic products by proton transfer reactions

Z would act as an electron acceptor and P would capture a hole in the unnatural DNA. The latter process would produce mutagenic products via a proton transfer reaction.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Mechanism for Fluorescence Quenching of Tryptophan by Oxamate and Pyruvate: Conjugation and Solvation-Induced Photoinduced Electron Transfer

Oxamate and pyruvate are isoelectronic molecules. They both quench tryptophan fluorescence with Stern-Volmer constants of 16 and 20 M^-1, respectively, which are comparable to that of arcrylamide, a commonly used probe for protein structure. On the other hand, it is well known that neither the carboxylate group of these molecules nor the amide group is a good quencher. To find the mechanism of the quenching by oxamate and pyruvate, density functional theory computations with a polarizable continuum model, solvation based on density, and explicit waters, were performed. Results indicate that both molecules can be an electron acceptor via photoinduced electron transfer. There are two requirements. First, the carboxylate and amide moieties must be in direct contact to bring about noticeable quenching. The conjugation between the amide (or the keto) group and the carboxylate group leads to a lower π* orbital, which is the lowest unoccupied molecular orbital (LUMO), and can then accept an electron from the excited tryptophan. Second, since oxamate and pyruvate ions have high electron density, hydrogen bonds with waters, which can be simulated by an explicit water model, are essential. Their LUMO energies are strongly influenced by water in aqueous solution. The above findings demonstrate how tryptophan fluorescence gets quenched in aqueous solution. The findings may be important in dealing with those problems where frontier orbitals are considered, especially with molecules having high electron density.

Bifurcate localization modes of excess electron in aqueous Ca(2+)amide solution revealed by ab initio molecular dynamics simulation: towards hydrated electron versus hydrated amide anion

In this work, we conduct ab initio molecular dynamics simulations on the localization dynamics of an excess electron (EE) in acetamide/Ca(2+) aqueous solutions with three different interaction modes of Ca(2+) with acetamide: tight contact, solvent-shared state, and separated interaction. The simulated results reveal that an EE could exhibit two different localization behaviors in these acetamide/Ca(2+) aqueous solutions depending on different amideCa(2+) interactions featuring different contact distances. For the tight contact and solvent-shared state of amideCa(2+) solutions, vertically injected diffuse EEs follow different mechanisms with different dynamics, forming a cavity-shaped hydrated electron or a hydrated amide anion, respectively. Meanwhile, for the separated state, only one localization pattern of a vertically injected diffuse EE towards the formation of hydrated amide anion is observed. The hindrance of hydrated Ca(2+) and the attraction of the hydrated amide group originating from its polarity and low energy π* orbital are the main driving forces. Additionally, different EE localization modes have different effects on the interaction between the amide group and Ca(2+) in turn. This work provides an important basis for further understanding the mechanisms and dynamics of localizations/transfers of radiation-produced EEs and associated EE-induced lesions and damage to biological species in real biological environments or other aqueous solutions.

Microhydration effects on the structures and electrophilic properties of cytidine

Adiabatic electron affinities (AEAs) for cytidine hydrates with up to four water molecules.

Quantum chemical MP2 results on some hydrates of cytosine: binding sites, energies and the first hydration shell

A detailed quantum chemical investigation was undertaken to obtain the structure and energetics of cytosine hydrates Cyt·nH2O, with n = 1 to 7. The MP2(fc)/aug-cc-pVDZ level was used as the standard, with some DFT (B3LYP) and coupled cluster calculations, as well as calculations with the aug-cc-pVTZ basis set added for comparison. In a systematic search for microhydrated forms of cytosine, we have found that several structures have not yet been reported in the literature. The energies of different isomers, as well as binding energies are compared. When predicting the stability of a complex, we suggest using a scheme where the water molecules are extracted from a finite model of bulk water. Finally, based on energetic data, we suggest a rational definition of the first hydration shell; with this definition, it contains just six water molecules.

Protection of DNA against low-energy electrons by amino acids: a first-principles molecular dynamics study

The physical shielding and chemical stabilizing of the low energy electron by glycine around the thymine through proton transfer.

Structures, electrophilic properties, and hydrogen bonds of cytidine, uridine, and their radical anions: Microhydration effects

Structures, electrophilic properties, and hydrogen bonds of the neutral and anionic monohydrated nucleoside, (cytidine)H2O, and (uridine)H2O have been systematically investigated using density functional theory. Various water-binding sites were predicted by explicitly considering the optimized monohydrated structures. Meanwhile, predictions of electron affinities and vertical detachment energies were also carried out to investigate their electrophilic properties. By examining the singly occupied molecular orbital and natural population analysis, we found the excess negative charge is localized on the cytidine and uridine moiety in anionic monohydrates. This may be the reason why the strength of hydrogen bonding undergoes an obvious change upon the extra electron attachment. Based on the electron density (ρ) and reduced density gradient (RDG), we present an approach to map and analyze the weak interaction (especially hydrogen bond) in monohydrated cytidine and uridine. The scatter plots of RDG versus ρ allow us to identify the different type interactions. Meanwhile, the maps of the gradient isosurfaces show a rich visualization of hydrogen bond, van der Waals interaction, and steric effect.

Electron attachment to the guanine-cytosine nucleic acid base pair and the effects of monohydration and proton transfer

The guanine-cytosine (GC) radical anion and its interaction with a single water molecule is studied using ab initio and density functional methods. Z-averaged second-order perturbation theory (ZAPT2) was applied to GC radical anion for the first time. Predicted spin densities show that the radical character is localized on cytosine. The Watson-Crick monohydrated GC anion is compared to neutral GC·H2O, as well as to the proton-transferred analogue on the basis of structural and energetic properties. In all three systems, local minima are identified that correspond to water positioned in the major and minor grooves of macromolecular DNA. On the anionic surface, two novel structures have water positioned above or below the GC plane. On the neutral and anionic surfaces, the global minimum can be described as water interacting with the minor groove. These structures are predicted to have hydration energies of 9.7 and 11.8 kcal mol(-1), respectively. Upon interbase proton-transfer (PT), the anionic global minimum has water positioned in the major groove, and the hydration energy increases to 13.4 kcal mol(-1). PT GC·H2O(•-) has distonic character; the radical character resides on cytosine, while the negative charge is localized on guanine. The effects of proton transfer are further investigated through the computed adiabatic electron affinities (AEA) of GC and monohydrated GC, and the vertical detachment energies (VDE) of the corresponding anions. Monohydration increases the AEAs and VDEs by only 0.1 eV, while proton-transfer increases the VDEs substantially (0.8 eV). The molecular charge distribution of monohydrated guanine-cytosine radical anion depends heavily on interbase proton transfer.

Excess electron localization in solvated DNA bases

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Electron attachment to a hydrated DNA duplex: the dinucleoside phosphate deoxyguanylyl-3',5'-deoxycytidine

The minimal essential section of DNA helices, the dinucleoside phosphate deoxyguanylyl-3',5'-deoxycytidine dimer octahydrate, [dGpdC](2), has been constructed, fully optimized, and analyzed by using quantum chemical methods at the B3LYP/6-31+G(d,p) level of theory. Study of the electrons attached to [dGpdC](2) reveals that DNA double strands are capable of capturing low-energy electrons and forming electronically stable radical anions. The relatively large vertical electron affinity (VEA) predicted for [dGpdC](2) (0.38 eV) indicates that the cytosine bases are good electron captors in DNA double strands. The structure, charge distribution, and molecular orbital analysis for the fully optimized radical anion [dGpdC](2)(·-) suggest that the extra electron tends to be redistributed to one of the cytosine base moieties, in an electronically stable structure (with adiabatic electron affinity (AEA) 1.14 eV and vertical detachment energy (VDE) 2.20 eV). The structural features of the optimized radical anion [dGpdC](2)(·-) also suggest the probability of interstrand proton transfer. The interstrand proton transfer leads to a distonic radical anion [d(G-H)pdC:d(C+H)pdG](·-), which contains one deprotonated guanine anion and one protonated cytosine radical. This distonic radical anion is predicted to be more stable than [dGpdC](2)(·-). Therefore, experimental evidence for electron attachment to the DNA double helices should be related to [d(G-H)pdC:d(C+H)pdG](·-) complexes, for which the VDE might be as high as 2.7 eV (in dry conditions) to 3.3 eV (in fully hydrated conditions). Effects of the polarizable medium have been found to be important for increasing the electron capture ability of the dGpdC dimer. The ultimate AEA value for cytosine in DNA duplexes is predicted to be 2.03 eV in aqueous solution.

Dielectric properties of fully hydrated nucleotides in the terahertz frequency range

We use terahertz time domain spectroscopy (THz-TDS) to determine the complex frequency-dependent dielectric response of all four nucleotides at different dilute concentrations. In addition, the suspension model's ability to extract the dielectric response of just the nucleotide with the hydration shell epsilon(b) excluding the dielectric information relating to the bulk will be verified. The suspension model enables us to make the determination that the nucleotides have influences on the water molecules out to the fourth hydration shell. We use a two Debye relaxation fit model for water, all concentrations and all epsilon(b) values. We observed how the nucleotides affect the relaxation parameters in relation to that of pure bulk water. With this information, we notice a transition between purines and pyrimidines, where one is a hydrogen-bond network structure building type material with a low concentration increment and the other is a structure breaking type material with a low concentration decrement. Due to conductivity measurements, we determine that kinetic depolarization is a negligible affect compared to that of dielectric saturation, which we find to dominate where a decrement is found.

Electron attachment to hydrated oligonucleotide dimers: guanylyl-3',5'-cytidine and cytidylyl-3',5'-guanosine

The dinucleoside phosphate deoxycytidylyl-3',5'-deoxyguanosine (dCpdG) and deoxyguanylyl-3',5'-deoxycytidine (dGpdC) systems are among the largest to be studied by reliable theoretical methods. Exploring electron attachment to these subunits of DNA single strands provides significant progress toward definitive predictions of the electron affinities of DNA single strands. The adiabatic electron affinities of the oligonucleotides are found to be sequence dependent. Deoxycytidine (dC) on the 5' end, dCpdG, has larger adiabatic electron affinity (AEA, 0.90 eV) than dC on the 3' end of the oligomer (dGpdC, 0.66 eV). The geometric features, molecular orbital analyses, and charge distribution studies for the radical anions of the cytidine-containing oligonucleotides demonstrate that the excess electron in these anionic systems is dominantly located on the cytosine nucleobase moiety. The pi-stacking interaction between nucleobases G and C seems unlikely to improve the electron-capturing ability of the oligonucleotide dimers. The influence of the neighboring base on the electron-capturing ability of cytosine should be attributed to the intensified proton accepting-donating interaction between the bases. The present investigation demonstrates that the vertical detachment energies (VDEs) of the radical anions of the oligonucleotides dGpdC and dCpdG are significantly larger than those of the corresponding nucleotides. Consequently, reactions with low activation barriers, such as those for O-C sigma bond and N-glycosidic bond breakage, might be expected for the radical anions of the guanosine-cytosine mixed oligonucleotides.

Elementary lesions in DNA subunits: electron, hydrogen atom, proton, and hydride transfers

When DNA is damaged by ionizing radiation, the genes in a cell may acquire mutations or the cell could die. The smallest known DNA-damaging unit is an electron, often low-energy secondary electrons. Additional electrons and transfers involving hydrogen atoms, protons, and hydride anions can damage DNA subunits, including individual nucleobases and nucleoside pairs. Researchers would like to better understand the molecular mechanisms involved in DNA damage from ionizing radiation. In this Account, we highlight our theoretical investigations of the molecular mechanisms of DNA damage using quantum mechanical models. Our investigations use robust theoretical methods with computations conducted in the gas phase and with solution models. We calculate adiabatic electron affinities (AEAs), which describe the energetics of electronic attachment to closed-shell DNA subunits, for the free bases, nucleosides, nucleotides, base pairs, and single and double DNA strand units. Electron affinities for free nucleobases yield the order uracil > thymine > cytosine > guanine > adenine and the same order for the DNA nucleosides, mononucleotides, and nucleoside 3',5'-diphosphates. AEA values increase steadily with the size and complexity of the system because of greater hydration, glycosylation, nucleotide formation, and base pairing. We predict and experimental results partially confirm that most of the more complex and hydrated species are observable as radical anions. Our modeling studies indicate that depyrimidination reactions of radical anion nucleosides release cytosine more often than thymine. Recent experimental results support those findings. Our theoretical studies of DNA base-pair radical anions predict increases in electron affinity accompanying H bonding and solvation. Electron addition facilitates some proton transfers in these pairs, which results in strongly perturbed pairing configurations. Of all nucleobase moieties within the more complex radical anion systems, thymine is best able to retain a negative charge. Charge and spin are well-separated in some of these systems. Radical species derived via hydrogen abstraction from DNA subunits yield large AEA values because they form closed-shell anions. Our studies predict single-strand breaks following H abstraction from nucleotides. Some H-abstraction processes in the DNA base pairs lead to severe distortions in pairing configuration based on our calculations. This body of systematic theoretical studies provides realistic descriptions of some events that lead to elementary DNA lesions, while providing rationalizations for many observed phenomena. Such approaches can lead to the design of new experiments, which would contribute to our understanding of the chemical physics of nucleic acids.

Formation of spectral intermediate G-C and A-T anion complex in duplex DNA studied by pulse radiolysis

The dynamics of electron adducts of 2'-deoxynucleotides and oligonucelotides (ODNs) were measured spectroscopically by nanosecond pulse radiolysis. The radical anions of the nucleotides were produced within 10 ns by the reaction of hydrated electrons (e(aq)(-)) and were protonated to form the corresponding neutral radicals. At pH 7.0, the radical anion of deoxythymidine (dT(*-)) was protonated to form the neutral radical dT(H)(*) in the time range of microseconds. The rate constant for the protonation was determined as 1.8 x 10(10) M(-1) s(-1). In contrast, the neutral radical of dC(H)(*) was formed immediately after the pulse, suggesting that the protonation occurs within 10 ns. The transient spectra of excess electrons of the double-stranded ODNs 5'-TAATTTAATAT-3' (AT) and 5'-CGGCCCGGCGC-3' (GC) differed from those of pyrimidine radicals (C and T) and their composite. In contrast, the spectra of the electron adducts of the single-stranded ODNs GC and AT exhibited characteristics of C and T, respectively. These results suggest that, in duplex ODNs, the spectral intermediates of G-C and A-T anions complex were formed. On the microsecond time scale, the subsequent changes in absorbance of the ODN AT had a first-order rate constant of 4 x 10(4) s(-1), reflecting the protonation of T.

Effects of Microsolvation on the Adenine−Uracil Base Pair and Its Radical Anion: Adenine−Uracil Mono- and Dihydrates

Microhydration effects upon the adenine-uracil (AU) base pair and its radical anion have been investigated by explicitly considering various structures of their mono- and dihydrates at the B3LYP/DZP++ level of theory. For the neutral AU base pair, 5 structures were found for the monohydrate and 14 structures for the dihydrate. In the lowest-energy structures of the neutral mono- and dihydrates, one and two water molecules bind to the AU base pair through a cyclic hydrogen bond via the N(9)-H and N(3) atoms of the adenine moiety, while the lowest-lying anionic mono- and dihydrates have a water molecule which is involved in noncyclic hydrogen bonding via the O4 atom of the uracil unit. Both the vertical detachment energy (VDE) and adiabatic electron affinity (AEA) of the AU base pair are predicted to increase upon hydration. While the VDE and AEA of the unhydrated AU pair are 0.96 and 0.40 eV, respectively, the corresponding predictions for the lowest-lying anionic dihydrates are 1.36 and 0.75 eV, respectively. Because uracil has a greater electron affinity than adenine, an excess electron attached to the AU base pair occupies the pi* orbital of the uracil moiety. When the uracil moiety participates in hydrogen bonding as a hydrogen bond acceptor (e.g., the N(6)-H(6a)...O(4) hydrogen bond between the adenine and uracil bases and the O(w)-H(w)...N and O(w)-H(w)...O hydrogen bonds between the AU pair and the water molecules), the transfer of the negative charge density from the uracil moiety to either the adenine or water molecules efficiently stabilizes the system. In addition, anionic structures which have C-H...O(w) contacts are energetically more favorable than those with N-H...O(w) hydrogen bonds, because the C-H...O(w) contacts do not allow the unfavorable electron density donation from the water to the uracil moiety. This delocalization effect makes the energetic ordering for the anionic hydrates very different from that for the corresponding neutrals.