A theoretical study of structures and electron affinities of radical anions of guanine-cytosine, adenine-thymine, and hypoxanthine-cytosine base pairs

CEBA: A new heterobifunctional reagent for plasmid DNA functionalization by click chemistry

Here, we report the synthesis of 3,6,9-trioxaundecan-1-{4-[(2-Chloroethyl)Ethylamino)]-Benzylamino},11-Azide (CEBA). CEBA alkylates the N7 of guanine of DNA thanks its chloroethyl group and can be coupled by a strain-promoted azide-alkyne cycloaddition to an alkynylated molecule. The optimization of the alkylation level of pDNA reveals that the expression of the encoded gene is preserved when it is randomly modified with at most 1 CEBA molecule per 150 bp. We show that the azido group of CEBA allows the linkage via click chemistry of CEBA-pDNA with a fluorophore or a peptide containing a dibenzocyclooctyne (DBCO) function. This new heterobifunctional reagent opens new ways to equip pDNA easily with signal molecules including peptides and nucleic acids without side products providing great interest for non-viral gene therapy.

Quantifying charge transfer energies at donor-acceptor interfaces in small-molecule solar cells with constrained DFTB and spectroscopic methods

Charge transfer states around the donor-acceptor interface in an organic solar cell determine the device performance in terms of the open circuit voltage. In the present work, we propose a computational scheme based on constrained density functional tight binding theory (c-DFTB) to assess the energy of the lowest charge transfer (CT) state in such systems. A comparison of the c-DFTB scheme with Hartree-Fock based configuration interaction of singles (CIS) and with time-dependent density functional theory (TD-DFT) using the hybrid functional B3LYP reveals that CIS and c-DFTB reproduce the correct Coulomb asymptotics between cationic donor and anionic acceptor configurations, whereas TD-DFT gives a qualitatively wrong excitation energy. Together with an embedding scheme accounting for the polarizable medium, this c-DFTB scheme is applied to several donor-acceptor combinations used in molecular solar cells. The external quantum efficiency of photovoltaic cells based on zinc phthalocyanine-C60 blends reveals a CT band remaining much narrower than the density of states of acceptor HOMO and donor LUMO, an observation which can be interpreted in a natural way in terms of Marcus transfer theory. A detailed comparison with c-DFTB calculations reveals an energy difference of 0.32 eV between calculated and observed absorption from the electronic ground state into the CT state. In a blend of a functionalized thiophene and C60, the photoluminescence spectra differ significantly from neat films, allowing again an assignment to CT states. The proposed computational scheme reproduces the observed trends of the observed open circuit voltages in photovoltaic devices relying on several donor-acceptor blends, finding an offset of 1.16 eV on average. This value is similar as in polymer-fullerene photovoltaic systems where it amounts to about 0.9 eV, indicating that the photophysics of CT states in molecular donor-acceptor blends and in polymer-fullerene blends are governed by the same mechanisms.

Photoelectron spectroscopy and density functional theory studies on the uridine homodimer radical anions

We report the photoelectron spectrum (PES) of the homogeneous dimer anion radical of uridine, (rU)(2)(●-). It features a broad band consisting of an onset of ∼1.2 eV and a maximum at the electron binding energy (EBE) ranging from 2.0 to 2.5 eV. Calculations performed at the B3LYP∕6-31++G∗∗ level of theory suggest that the PES is dominated by dimeric radical anions in which one uridine nucleoside, hosting the excess charge on the base moiety, forms hydrogen bonds via its O8 atom with hydroxyl of the other neutral nucleoside's ribose. The calculated adiabatic electron affinities (AEAGs) and vertical detachment energies (VDEs) of the most stable homodimers show an excellent agreement with the experimental values. The anionic complexes consisting of two intermolecular uracil-uracil hydrogen bonds appeared to be substantially less stable than the uracil-ribose dimers. Despite the fact that uracil-uracil anionic homodimers are additionally stabilized by barrier-free electron-induced proton transfer, their relative thermodynamic stabilities and the calculated VDEs suggest that they do not contribute to the experimental PES spectrum of (rU)(2)(●-).

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.

Stability of the valence anion of cytosine is governed by nucleobases sequence in the double stranded DNA pi-stack: A computational study

The stabilities of the valence anion of cytosine (C(-)) in model trimers of complementary base pairs that possess the B-DNA geometry but differ in base sequence are reported. In order to estimate the energetics of electron attachment to the middle cytosine incorporated in the trimer, a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated cytosine were developed. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. We have demonstrated that contrary to the literature reports, concerning single stranded DNA, the sequence of nucleic bases has a profound effect on the stability of the cytosine valence anion. The anionic 3(')-CCC-5(') complex is the most stable configuration (EA=0.399 eV) and the 3(')-GCG-5(') trimer anion is the most unstable species (EA=-0.193 eV). Moreover, with the energetic correction for the presence of sugar-phosphate backbone all possible double stranded DNA sequences lead to the stable C(-). The predicted electron affinities of the cytosine anion have been compared to the results of analogous studies on the thymine anion published recently [M. Kobyłecka et al., J. Am. Chem. Soc. 130, 15683 (2008)]. The consequences of low-energy barrier proton transfer in the GC anion have been discussed in the context of induced by electrons DNA single strand breaks. The DNA sequences that should dramatically differ in their vulnerability to be damaged by low energy electrons have been proposed.

Valence anion of thymine in the DNA pi-stack

Most of theoretical data on the stability of radical anions supported by nucleic acid bases have been obtained for anions of isolated nucleobases, their nucleosides, or nucleotides. This approach ignores the hallmark forces of DNA, namely, hydrogen bonding and pi-stacking interactions. Since these interactions might be crucial for the electron affinities of nucleobases bound in DNA, we report for the first time on the stability of the thymine valence anion in trimers of complementary bases possessing the regular B-DNA geometry but differing in base sequence. In order to estimate the energetics of electron attachment to a trimer, we developed a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated thymine. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. The two-body interaction energies were corrected for the basis set superposition error, and in benchmark systems, they were extrapolated to the basis set limit and supplemented with correction for higher order correlation terms calculated at the CCSD(T) level. We have demonstrated that the sequence of nucleic bases has a profound effect on the stability of the thymine valence anion: the anionic 5'-CTC-3' (6.0 kcal/mol) sequence is the most stable configuration, and the 5'-GTG-3' (-8.0 kcal/mol) trimer anion is the most unstable species. On the basis of obtained results, one can propose DNA sequences that are different in their vulnerability to damage by low energy electron.

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.

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.

Combined quantum mechanical and molecular mechanical simulations of one- and two-electron reduction potentials of flavin cofactor in water, medium-chain acyl-CoA dehydrogenase, and cholesterol oxidase

Flavin adenine dinucleotide (FAD) is a common cofactor in redox proteins, and its reduction potentials are controlled by the protein environment. This regulation is mainly responsible for the versatile catalytic functions of flavoenzymes. In this article, we report computations of the reduction potentials of FAD in medium-chain acyl-CoA dehydrogenase (MCAD) and cholesterol oxidase (CHOX). In addition, the reduction potentials of lumiflavin in aqueous solution have also been computed. Using molecular dynamics and free-energy perturbation techniques, we obtained the free-energy changes for two-electron/two-proton as well as one-electron/one-proton addition steps. We employed a combined quantum mechanical and molecular mechanical (QM/MM) potential, in which the flavin ring was represented by the self-consistent-charge density functional tight-binding (SCC-DFTB) method, while the rest of the enzyme-solvent system was treated by classical force fields. The computed two-electron/two-proton reduction potentials for lumiflavin and the two enzyme-bound FADs are in reasonable agreement with experimental data. The calculations also yielded the pKa values for the one-electron reduced semiquinone (FH*) and the fully reduced hydroquinone (FH2) forms. The pKa of the FAD semiquinone in CHOX was found to be around 4, which is 4 units lower than that in the enzyme-free state and 2 units lower than that in MCAD; this supports the notion that oxidases have a greater ability than dehydrogenases to stabilize anionic semiquinones. In MCAD, the flavin ring interacts with four hydrophobic residues and has a significantly bent structure, even in the oxidized state. The present study shows that this bending of the flavin imparts a significant destabilization (approximately 5 kcal/mol) to the oxidized state. The reduction potential of lumiflavin was also computed using DFT (M06-L and B3LYP functionals with 6-31+G(d,p) basis set) with the SM6 continuum solvation model, and the results are in good agreement with results from explicit free-energy simulations, which supports the conclusion that the SCC-DFTB/MM computation is reasonably accurate for both 1e(-)/1H+ and 2e(-)/2H+ reduction processes. These results suggest that the first coupled electron-proton addition is stepwise for both the free and the two enzyme-bound flavins. In contrast, the second coupled electron-proton addition is also stepwise for the free flavin but is likely to be concerted when the flavin is bound to either the dehydrogenase or the oxidase enzyme.

Combined quantum mechanical and molecular mechanical simulations of one- and two-electron reduction potentials of flavin cofactor in water, medium-chain acyl-CoA dehydrogenase, and cholesterol oxidase

Flavin adenine dinucleotide (FAD) is a common cofactor in redox proteins, and its reduction potentials are controlled by the protein environment. This regulation is mainly responsible for the versatile catalytic functions of flavoenzymes. In this article, we report computations of the reduction potentials of FAD in medium-chain acyl-CoA dehydrogenase (MCAD) and cholesterol oxidase (CHOX). In addition, the reduction potentials of lumiflavin in aqueous solution have also been computed. Using molecular dynamics and free-energy perturbation techniques, we obtained the free-energy changes for two-electron/two-proton as well as one-electron/one-proton addition steps. We employed a combined quantum mechanical and molecular mechanical (QM/MM) potential, in which the flavin ring was represented by the self-consistent-charge density functional tight-binding (SCC-DFTB) method, while the rest of the enzyme-solvent system was treated by classical force fields. The computed two-electron/two-proton reduction potentials for lumiflavin and the two enzyme-bound FADs are in reasonable agreement with experimental data. The calculations also yielded the pKa values for the one-electron reduced semiquinone (FH*) and the fully reduced hydroquinone (FH2) forms. The pKa of the FAD semiquinone in CHOX was found to be around 4, which is 4 units lower than that in the enzyme-free state and 2 units lower than that in MCAD; this supports the notion that oxidases have a greater ability than dehydrogenases to stabilize anionic semiquinones. In MCAD, the flavin ring interacts with four hydrophobic residues and has a significantly bent structure, even in the oxidized state. The present study shows that this bending of the flavin imparts a significant destabilization (approximately 5 kcal/mol) to the oxidized state. The reduction potential of lumiflavin was also computed using DFT (M06-L and B3LYP functionals with 6-31+G(d,p) basis set) with the SM6 continuum solvation model, and the results are in good agreement with results from explicit free-energy simulations, which supports the conclusion that the SCC-DFTB/MM computation is reasonably accurate for both 1e(-)/1H+ and 2e(-)/2H+ reduction processes. These results suggest that the first coupled electron-proton addition is stepwise for both the free and the two enzyme-bound flavins. In contrast, the second coupled electron-proton addition is also stepwise for the free flavin but is likely to be concerted when the flavin is bound to either the dehydrogenase or the oxidase enzyme.

Density Functional Study toward Understanding Dehydrogenation of the Adenine−Thymine Base Pair and Its Anion

The dehydrogenated radicals and anions of Watson-Crick adenine-thymine (A-T) base pair have been investigated by the B3LYP/DZP++ approach. Calculations show that the dehydrogenated radicals and anions have relatively high stabilities compared with the single base adenine and thymine. The electron attachment to the A-T base pair and its derivatives significantly modifies the hydrogen bond interactions and results in remarkable structural changes. As for the dehydrogenated A-T radicals, they have relatively high electron affinities and different dehydrogenation properties with respect to their constituent elements. The relatively low-cost hydrogen eliminations correspond to the (N9)-H (adenine) and (N1)-H (thymine) bonds cleavage. Both dehydrogenation processes have Gibbs free energies of reaction DeltaG degrees of 13.4 and 17.2 kcal mol-1, respectively. The solvent water exhibits significant effect on electron attachment and dehydrogenation properties of the A-T base pair and its derivatives. In the dehydrogenating process, the anionic A-T fragment gradually changes its electronic configuration from pi* to sigma* state, like the single bases adenine and thymine.

Interaction of Guanine, Its Anions, and Radicals with Lysine in Different Charge States

Modification in DNA or protein structure can severely affect DNA-protein interactions and the functioning of biological systems. Some new insights into radiation-induced effects of guanine-lysine interactions have been obtained here by theoretical investigations. Geometries of zwitterionic and non-zwitterionic lysine in different charge states (neutral, radical cation, and protonated cation) were optimized employing the B3LYP/6-31G** and B3LYP/AUG-cc-pVDZ levels of hybrid density functional theory (DFT) and using the second-order Møller-Plesset perturbation theory along with the 6-31G** basis set. In the case of neutral lysine in the gas phase, no zwitterionic structure was obtained. The non-zwitterionic structures of lysine in radical and protonated cationic forms are appreciably more stable than the corresponding zwitterionic structures in the gas phase as obtained at all levels of theory employed here. Binding of guanine and different dehydrogenated guanine radicals with lysine in different charge states was studied at the B3LYP/6-31G** level of DFT. When guanine makes a complex with the lysine radical cation, large amounts of spin and positive charge densities are transferred from the lysine radical cation to guanine and the guanine is thus converted from its normal form to the radical cationic form. Complexation of the lysine radical cation with the H1-hydrogen-abstracted guanine radical leads to CO2 liberation and proton transfer from lysine. These results are compared with the available experimental ones.

Does the A·T or G·C Base-Pair Possess Enhanced Stability? Quantifying the Effects of CH···O Interactions and Secondary Interactions on Base-Pair Stability Using a Phenomenological Analysis and ab Initio Calculations

An empirically based relationship between overall complex stability (-DeltaG degrees ) and various possible component interactions is developed to probe the question of whether the A.T/U and G.C base-pairs exhibit enhanced stability relative to similarly hydrogen-bonded complexes. This phenomenological approach suggests ca. 2-2.5 kcal mol-1 in additional stability for A.T owing to a group interaction containing a CH...O contact. Pairing geometry and the role of the CH...O interaction in the A.T base-pair were also probed using MP2/6-31+G(d,p) calculations and a double mutant cycle. The ab initio studies indicated that Hoogsteen geometry is preferred over Watson-Crick geometry in A.T by ca. 1 kcal mol-1. Factors that might contribute to the preference for Hoogsteen geometry are a shorter CH...O contact, a favorable alignment of dipoles, and greater distances between secondary repulsive sites. The CH...O interaction was also investigated in model complexes of adenine with ketene and isocyanic acid. The ab initio calculations support the result of the phenomenological approach that the A.T base-pair does have enhanced stability relative to hydrogen-bonded complexes with just N-H...N and N-H...O hydrogen bonds.

The guanine cation radical: investigation of deprotonation states by ESR and DFT

This work reports ESR studies that identify the favored site of deprotonation of the guanine cation radical (G*+) in an aqueous medium at 77 K. Using ESR and UV-visible spectroscopy, one-electron oxidized guanine is investigated in frozen aqueous D2O solutions of 2'-deoxyguanosine (dGuo) at low temperatures at various pHs at which the guanine cation radical, G*+ (pH 3-5), singly deprotonated species, G(-H)* (pH 7-9), and doubly deprotonated species, G(-2H)*- (pH > 11), are found. C-8-deuteration of dGuo to give 8-D-dGuo removes the major proton hyperfine coupling at C-8. This isolates the anisotropic nitrogen couplings for each of the three species and aids our analyses. These anisotropic nitrogen couplings were assigned to specific nitrogen sites by use of 15N-substituted derivatives at N1, N2, and N3 atoms in dGuo. Both ESR and UV-visible spectra are reported for each of the species: G*+, G(-H)*, and G(-2H)*-. The experimental anisotropic ESR hyperfine couplings are compared to those obtained from DFT calculations for the various tautomers of G(-H)*. Using the B3LYP/6-31G(d) method, the geometries and energies of G*+ and its singly deprotonated state in its two tautomeric forms, G(N1-H)* and G(N2-H)*, were investigated. In a nonhydrated state, G(N2-H)* is found to be more stable than G(N1-H)*, but on hydration with seven water molecules G(N1-H)* is found to be more stable than G(N2-H)*. The theoretically calculated hyperfine coupling constants (HFCCs) of G*+, G(N1-H)*, and G(-2H)*- match the experimentally observed HFCCs best on hydration with seven or more waters. For G(-2H)*-, the hyperfine coupling constant (HFCC) at the exocyclic nitrogen atom (N2) is especially sensitive to the number of hydrating water molecules; good agreement with experiment is not obtained until nine or 10 waters of hydration are included.

Binding of Gold Clusters with DNA Base Pairs: A Density Functional Study of Neutral and Anionic GC−Aun and AT−Aun (n = 4, 8) Complexes

Binding of clusters of gold atoms (Au) with the guanine-cytosine (GC) and adenine-thymine (AT) Watson-Crick DNA base pairs was studied using the density functional theory (DFT). Geometries of the neutral GC-Au(n) and AT-Au(n) and the corresponding anionic (GC-Au(n))(-1) and (AT-Au(n))(-1) (n = 4, 8) complexes were fully optimized in different electronic states, that is, singlet and triplet states for the neutral complexes and doublet and quartet states for the anionic complexes, using the B3LYP density functional method. The 6-31+G basis set was used for all atoms except gold. For gold atoms, the Los Alamos effective core potential (ECP) basis set LanL2DZ was employed. Vibrational frequency calculations were performed to ensure that the optimized structures corresponded to potential energy surface minima. The gold clusters around the neutral GC and AT base pairs have a T-shaped structure, which satisfactorily resemble those observed experimentally and in other theoretical studies. However, in anionic GC and AT base pairs, the gold clusters have extended zigzag and T-shaped structures. We found that guanine and adenine have high affinity for Au clusters, with their N3 and N7 sites being preferentially involved in binding with the same. The calculated adiabatic electron affinities (AEAs) of the GC-Au(n)complexes (n = 4, 8) were found to be much larger than those of the isolated base pairs.

The deprotonated guanine-cytosine base pair

Awareness of the harmful effects of radiation has increased interest in finding the mechanisms of DNA damage. Radical and anion formation among the DNA base pairs are thought to be important steps in such damage [Collins, G. P. (2003) Sci. Am. 289 (3), 26-27]. Energetic properties and optimized geometries of 10 radicals and their respective anions derived through hydrogen abstraction from the Watson-Crick guanine-cytosine (G-C) base pair have been studied using reliable theoretical methods. The most favorable deprotonated structure (dissociation energy 42 kcal x mol(-1), vertical detachment energy 3.79 eV) ejects the proton analogous to the cytosine glycosidic bond in DNA. This structure is a surprisingly large 12 kcal x mol(-1) lower in energy than any of the other nine deprotonated G-C structures. This system retains the qualitative G-C structure but with the H...O2 distance dramatically reduced from 1.88 to 1.58 A, an extremely short hydrogen bond. The most interesting deprotonated G-C structure is a "reverse wobble" incorporating two N-H...N hydrogen bonds. Three different types of relaxation energies (4.3-54 kcal x mol(-1)) are defined and reported to evaluate the energy released via different mechanisms for the preparation of the deprotonated species. Relative energies, adiabatic electron affinities (ranging from 1.93 to 3.65 eV), and pairing energies are determined to discern which radical will most alter the G-C properties. The most stable deprotonated base pair corresponds to the radical with the largest adiabatic electron affinity, 3.65 eV. This value is an enormous increase over the electron affinity (0.60 eV) of the closed-shell G-C base pair.

Microsolvation Effect, Hydrogen-Bonding Pattern, and Electron Affinity of the Uracil−Water Complexes U−(H2O)n (n = 1, 2, 3)

To achieve a systematic understanding of the influence of microsolvation on the electron accepting behaviors of nucleobases, the reliable theoretical method (B3LYP/DZP++) has been applied to a comprehensive conformational investigation on the uracil-water complexes U-(H(2)O)(n) (n = 1, 2, 3) in both neutral and anionic forms. For the neutral complexes, the conformers of hydration on the O2 of uracil are energetically favored. However, hydration on the O4 atom of uracil is more stable for the radical anions. The electron structure analysis for the H-bonding patterns reveal that the CH...OH(2) type H-bond exists only for di- and trihydrated uracil complexes in which a water dimer or trimer is involved. The electron density structure analysis and the atoms-in-molecules (AIM) analysis for U-(H(2)O)(n) suggest a threshold value of the bond critical point (BCP) density to justify the CH...OH(2) type H-bond; that is, CH...OH(2) could be considered to be a H-bond only when its BCP density value is equal to or larger than 0.010 au. The positive adiabatic electron affinity (AEA) and vertical detachment energy (VDE) values for the uracil-water complexes suggest that these hydrated uracil anions are stable. Moreover, the average AEA and VDE of U-(H(2)O)(n) increase as the number of the hydration waters increases.

Structural and Energetic Characterization of a DNA Nucleoside Pair and Its Anion: Deoxyriboadenosine (dA) − Deoxyribothymidine (dT)

The geometries of the DNA nucleoside pairs between 2'-deoxyriboadenosine (dA) and 2'-deoxyribothymidine (dT) and its anion (dAdT-) were fully optimized using carefully calibrated density functional methods. The addition of an electron to dAdT results in remarkable changes to the two hydrogen bonding distances, the H...O distance decreasing by 0.303 angstroms and the N...H distance increasing by 0.229 angstroms. The electron affinity of the dAdT pair was studied to reveal the correct trends of adiabatic electron affinity (EA(ad)) under the influence of the additional components to the individual bases. The consequence of negative charge in terms of structural variations, energetic changes, and charge distribution were explored. The EA(ad) of dAdT is predicted to be positive (0.60 eV), and it exhibits a substantial increase compared with those of the corresponding bases A and T and the nucleic acid base pair AT. The effects of pairing and the addition of the sugar moiety on the EA(ad) are well described as the summation of the individual influences. The influence of the pairing on the EA is comparable to that of the addition of 2-deoxyribose. The excess charge is mainly located on the thyminyl moiety in the anionic dAdT pair. The positive vertical electron affinity (VEA = 0.20 eV) for dAdT suggests that it is able to form a stable anion through electron attachment. A large vertical detachment energy (VDE = 1.14 eV) has been determined for the anionic dAdT nucleoside pair. Therefore, one may expect that the stable anionic dAdT nucleoside pair should be able to undergo the subsequent glycosidic bond cleavage process.

Marked Influences on the Adenine−Cytosine Base Pairs by Electron Attachment and Ionization

The reverse wobble and the reverse Hoogsteen adenine-cytosine mispairs regarding their radical cations and anions are studied with the hybrid three-parameter B3LYP density functional method and 6-31+G(d), 6-311+G(2df,2p) basis sets. Hydrogen bonding mispairs are remarkably influenced by electron attachment and ionization. Only one stronger hydrogen bond N6-H (in adenine)...N3 (in cytosine) exists in the radical pair, while the strengths of two N-H...N hydrogen bonds in the neutral pair are comparable. Geometrical coplanarity is found for the neutral and cationic pairs, in contrast to the anionic pairs in which the cytosine moiety exhibits significant deformation due to electron attachment. Dissociation energies for the neutral and radical pairs are slightly higher than those of the adenine-thymine pairs but much smaller than those of the guanine-cytosine pairs. Valence-bound anions of these two adenine-cytosine pairs are thermodynamically stable by 0.1-0.2 eV with respect to the neutral pairs. On the basis of the comparison between the experimental data of the solvated clusters and the calculated values, these two pairs can be quantitatively equivalent to the clusters in which each base is solvated by five water molecules.

Adiabatic Electron Affinities of the Polyhydrated Adenine−Thymine Base Pair: A Density Functional Study

Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.

AT Base Pair Anions versus (9-Methyl-A)(1-Methyl-T) Base Pair Anions

The anionic base pairs of adenine and thymine, (AT)(-), and 9-methyladenine and 1-methylthymine, (MAMT)(-), have been investigated both theoretically and experimentally in a complementary, synergistic study. Calculations on (AT)(-) found that it had undergone a barrier-free proton transfer (BFPT) similar to that seen in other dimer anion systems and that its structural configuration was neither Watson-Crick (WC) nor Hoogsteen (HS). The vertical detachment energy (VDE) of (AT)(-) was determined by anion photoelectron spectroscopy and found to be in agreement with the VDE value predicted by theory for the BFPT mechanism. An AT pair in DNA is structurally immobilized into the WC configuration, in part, by being bonded to the sugars of the double helix. This circumstance was mimicked by methylating the sites on both A and T where these sugars would have been tied, viz., 9-methyladenine and 1-methylthymine. Calculations found no BFPT in (MAMT)(-) and a resulting (MAMT)(-) configuration that was either HS or WC, with the configurations differing in stability by ca. 2 kcal/mol. The photoelectron spectrum of (MAMT)(-) occurred at a completely different electron binding energy than had (AT)(-). Moreover, the VDE value of (MAMT)(-) was in agreement with that predicted by theory. The configuration of (MAMT)(-) and its lack of electron-induced proton transfer are inter-related. While there may be other pathways for electron-induced DNA alterations, BFPT in the WC/HS configurations of (AT)(-) is not feasible.

A TDDFT Study of the Optical Response of DNA Bases, Base Pairs, and Their Tautomers in the Gas Phase

We present calculations of the optical response of the DNA bases and base pairs both in their normal and tautomeric forms in the gas phase, using time-dependent density functional theory (TDDFT). These calculations are performed in real time within the adiabatic approximation with a basis of local orbitals. Our results for the individual bases are in good agreement with experiment and computationally more demanding calculations of chemical accuracy. The optical response of base pairs indicates that the differences between normal and tautomeric forms in certain cases are significant enough to provide a means of identification.

Electron Binding to Nucleic Acid Bases. Experimental and Theoretical Studies. A Review

An in-depth knowledge of an excess electron binding mechanism to DNA and RNA nucleobases is important for our understanding of radiation damage influence on the biological functions of nucleic acids, as well as for the possible use of DNA molecules as wires in molecular electronic circuits. The of anions created by electron attachment to individual nucleic acid bases is discussed in detail. The principles of the experimental and theoretical approaches to the description of these anions are outlined, and the available results concerning valence- and dipole-bound anions of nucleic acid bases are reviewed. A review with 167 references.