On the unusual stability of valence anions of thymine based on very rare tautomers: A computational study

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).     

Photoelectron spectroscopy of adiabatically bound valence anions of rare tautomers of the nucleic acid bases

Anionic states of nucleic acid bases (NABs) are involved in DNA damage by low-energy electrons and in charge transfer through DNA. Previous gas phase studies of free, unsolvated NAB parent anions probed mostly dipole-bound states, which are not present in condensed phase environments. Recently, we demonstrated that very rare tautomers of uracil (U), cytosine (C), adenine (A), and guanine (G), which are obtained from canonical tautomers through N-to-C proton transfers, support valence anionic states. Here we report the photoelectron spectrum of the final member of the NABs series: the valence state of the thymine (T) anion. Additionally, we summarized the work of all five NABs. All of the newfound anionic tautomers of the NABs may be formed via dissociative electron attachment followed by hydrogen atom reattachment to a carbon atom. Furthermore, these unusual tautomers may affect the structure and properties of DNA and RNA exposed to low-energy electrons. The new valence states observed here, unlike dipole bound states, could exist in condensed phases and may be relevant to radiobiological damage.     

Stabilization of very rare tautomers of 1-methylcytosine by an excess electron

We characterized valence anionic states of 1-methylcytosine using various electronic structure methods. We found that the most stable valence anion is related to neither the canonical amino-oxo nor a rare imino-oxo tautomer, in which a proton is transferred from the N4 to N3 atom. Instead, it is related to an imino-oxo tautomer, in which the C5 atom is protonated. This anion is characterized by an electron vertical detachment energy (VDE) of 2.12 eV and it is more stable than the anion based on the canonical tautomer by 1.0 kcal/mol. The latter is characterized by a VDE of 0.31 eV. Another unusual low-lying imino-oxo tautomer with a VDE of 3.60 eV has the C6 atom protonated and is 3.6 kcal/mol less stable than the anion of the canonical tautomer. All these anionic states are adiabatically unbound with respect to the canonical amino-oxo neutral, with the instability of 5.8 kcal/mol for the most stable valence anion. The mechanism of formation of anionic tautomers with carbon atoms protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to the C5 or C6 atom. The six-member ring structure of anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Indeed the neutral systems collapse without a barrier to a linear or a bicyclo structure, which might be viewed as lesions to DNA or RNA. Within the PCM hydration model, the anions become adiabatically bound with respect to the corresponding neutrals, and the two most stable tautomers have a carbon atom protonated.     

Molecular structures and energetics of the (TiO2)n (n = 1-4) clusters and their anions

The (TiO2)n clusters and their anions for n = 1-4 have been studied with coupled cluster theory [CCSD(T)] and density functional theory (DFT). For n > 1, numerous conformations are located for both the neutral and anionic clusters, and their relative energies are calculated at both the DFT and CCSD(T) levels. The CCSD(T) energies are extrapolated to the complete basis set limit for the monomer and dimer and calculated up to the triple-zeta level for the trimer and tetramer. The adiabatic and vertical electron detachment energies of the anionic clusters to the ground and first excited states of the neutral clusters are calculated at both levels and compared with the experimental results. The comparison allows for the definitive assignment of the ground-state structures of the anionic clusters. Anions of the dimer and tetramer are found to have very closely lying conformations within 2 kcal/mol at the CCSD(T) level, whereas that of the trimer does not. In addition, accurate clustering energies and heats of formation are calculated for the neutral clusters and compared with the available experimental data. Estimates of the titanium-oxygen bond energies show that they are stronger than the group VIB transition metal-oxygen bonds except for tungsten. The atomization energies of these clusters display much stronger basis set dependence than the clustering energies. This allows the calculation of more accurate heats of formation for larger clusters on the basis of calculated clustering energies.     

Solvation free energies of molecules. The most stable anionic tautomers of uracil

Anionic states of nucleic acid bases are suspected to play a role in the radiation damage processes of DNA. Our recent studies suggested that the excess electron attachment to the nucleic acid bases can stabilize some rare tautomers, i.e. imine-enamine tautomers and other tautomers with a proton being transferred from nitrogen sites to carbon sites (with respect to the canonical tautomer). So far, these new anionic tautomers have been characterized by the gas-phase electronic structure calculations and photoelectron spectroscopy experiments. In the current contribution we explore the effect of water solvation on the stability of the new anionic tautomers of uracil. The accurate free energies of solvation are calculated in a two step approach. The major contribution was calculated using the classical free-energy perturbation adiabatic-charging approach, where it is assumed that the solvated molecule has the charge distribution given by the polarizable continuum model. In the second step the free energy of solvation is refined by taking into account the real, average solvent charge distribution. This is done using our accelerated QM/MM simulations, where the QM energy of the solute is calculated in the mean potential averaged over many MD steps. We found that in water solution three of the recently identified anionic tautomers are 6.5-3.6 kcal mol(-1) more stable than the anion of the canonical tautomer.     

Expanded-size bases in naturally sized DNA: evaluation of steric effects in Watson-Crick pairing

We describe physicochemical properties in DNA of altered-size nucleobases that retain Watson-Crick analogous hydrogen-bonding ability. Size-expanded analogues of adenine and thymine (xA and xT, respectively, which are expanded by benzo-fusion) were incorporated into natural DNA oligonucleotides, and their effects on helix stability were measured. Base stacking studies revealed that the two stretched analogues stack much more strongly than do their naturally sized counterparts. In contrast to this, pairing studies showed that single substitutions of the new bases are destabilizing to the natural helix as compared to A or T in standard A-T pairs in the same context, unless multiple adjacent substitutions are used. Interestingly, the size-expanded bases displayed selective recognition of the hydrogen-bonding complementary partners, suggesting that Watson-Crick analogous pairs were still formed despite local backbone strain. In an attempt to compensate for the added size of the expanded adenine, we tested a formamide deoxynucleoside, which Leonard proposed as a shortened thymine analogue (F(o)). Data showed, however, that this compound adopts a conformation unfavorable for pairing. On the basis of the combined thermodynamic data, we estimate the energetic cost of the 2.4 A stretching of an isolated base pair in DNA at ca. +1 to 2 kcal/mol. Notably, during the pairing studies, the two size-expanded nucleobases were found to display significant changes in fluorescence on formation of stacked versus unstacked structures, suggesting possible applications in probing nucleic acid structures and biochemical mechanisms.     

Optimized Slater-type basis sets for the elements 1-118

Seven different types of Slater type basis sets for the elements H (Z = 1) up to E118 (Z = 118), ranging from a double zeta valence quality up to a quadruple zeta valence quality, are tested in their performance in neutral atomic and diatomic oxide calculations. The exponents of the Slater type functions are optimized for the use in (scalar relativistic) zeroth-order regular approximated (ZORA) equations. Atomic tests reveal that, on average, the absolute basis set error of 0.03 kcal/mol in the density functional calculation of the valence spinor energies of the neutral atoms with the largest all electron basis set of quadruple zeta quality is lower than the average absolute difference of 0.16 kcal/mol in these valence spinor energies if one compares the results of ZORA equation with those of the fully relativistic Dirac equation. This average absolute basis set error increases to about 1 kcal/mol for the all electron basis sets of triple zeta valence quality, and to approximately 4 kcal/mol for the all electron basis sets of double zeta quality. The molecular tests reveal that, on average, the calculated atomization energies of 118 neutral diatomic oxides MO, where the nuclear charge Z of M ranges from Z = 1-118, with the all electron basis sets of triple zeta quality with two polarization functions added are within 1-2 kcal/mol of the benchmark results with the much larger all electron basis sets, which are of quadruple zeta valence quality with four polarization functions added. The accuracy is reduced to about 4-5 kcal/mol if only one polarization function is used in the triple zeta basis sets, and further reduced to approximately 20 kcal/mol if the all electron basis sets of double zeta quality are used. The inclusion of g-type STOs to the large benchmark basis sets had an effect of less than 1 kcal/mol in the calculation of the atomization energies of the group 2 and group 14 diatomic oxides. The basis sets that are optimized for calculations using the frozen core approximation (frozen core basis sets) have a restricted basis set in the core region compared to the all electron basis sets. On average, the use of these frozen core basis sets give atomic basis set errors that are approximately twice as large as the corresponding all electron basis set errors and molecular atomization energies that are close to the corresponding all electron results. Only if spin-orbit coupling is included in the frozen core calculations larger errors are found, especially for the heavier elements, due to the additional approximation that is made that the basis functions are orthogonalized on scalar relativistic core orbitals.     

Electron affinity of the guanine-cytosine base pair and structural perturbations upon anion formation

The adiabatic electron affinity (AEA) for the Watson-Crick guanine-cytosine (GC) DNA base pair is predicted using a range of density functional methods with double- and triple-zeta plus polarization plus diffuse (DZP++ and TZ2P++) basis sets in an effort to bracket the true electron affinity. The methods used have been calibrated against a comprehensive tabulation of experimental electron affinities (Chem.Rev. 2002, 102, 231). Optimized structures for GC and the GC anion are compared to the neutral and anionic forms of the individual bases as well as Rich's 1976 X-ray structure for sodium guanylyl-3',5'-cytidine nonahydrate, GpC.9H(2)O. Structural distortions and natural population (NPA) charge distributions of the GC anion indicate that the unpaired electron is localized primarily on the cytosine moiety. Unlike treatments using second-order perturbation theory (MP2), density functional theory consistently predicts a substantial positive adiabatic electron affinity for the GC pair (e.g., TZ2P++/B3LYP: +0.48 eV). The stabilization of C(-) via three hydrogen bonds to guanine is sufficient to facilitate adiabatic binding of an electron to GC and is also consistent with the positive experimental electron affinities obtained by photoelectron spectroscopy of cytosine anions incrementally microsolvated with water molecules. The pairing (dissociation) energy for GC(-) (35.6 kcal/mol) is determined with inclusion of electron correlation and shows the anion to have greater thermodynamic stability; the pairing energy for neutral GC (TZ2P++/B3LYP 23.9 kcal/mol) compares favorably to previous MP2/6-31G (23.4 kcal/mol) results and a debated experiment (21.0 kcal/mol).     

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.     

Barrier-free proton transfer in the valence anion of 2'-deoxyadenosine-5'-monophosphate. II. A computational study

The propensity of four representative conformations of 2(')-deoxyadenosine-5(')-monophosphate (5(')-dAMPH) to bind an excess electron has been studied at the B3LYP6-31++G(d,p) level. While isolated canonical adenine does not support stable valence anions in the gas phase, all considered neutral conformations of 5(')-dAMPH form adiabatically stable anions. The type of an anionic 5(')-dAMPH state, i.e., the valence, dipole bound, or mixed (valence/dipole bound), depends on the internal hydrogen bond(s) pattern exhibited by a particular tautomer. The most stable anion results from an electron attachment to the neutral syn-south conformer. The formation of this anion is associated with a barrier-free proton transfer triggered by electron attachment and the internal rotation around the C4(')-C5(') bond. The adiabatic electron affinity of the a_south-syn anion is 1.19 eV, while its vertical detachment energy is 1.89 eV. Our results are compared with the photoelectron spectrum (PES) of 5(')-dAMPH(-) measured recently by Stokes et al., [J. Chem. Phys. 128, 044314 (2008)]. The computational VDE obtained for the most stable anionic structure matches well with the experimental electron binding energy region of maximum intensity. A further understanding of DNA damage might require experimental and computational studies on the systems in which purine nucleotides are engaged in hydrogen bonding.     

Electron affinities of the DNA and RNA bases.

Adiabatic electron affinities (AEAs) for the DNA and RNA bases are predicted by using a range of density functionals with a double-zeta plus polarization plus diffuse (DZP++) basis set in an effort to bracket the true EAs. Although the AEAs exhibit moderate fluctuations with respect to the choice of functional, systematic trends show that the covalent uracil (U) and thymine (T) anions are bound by 0.05-0.25 eV while the adenine (A) anion is clearly unbound. The computed AEAs for cytosine (C) and guanine (G) oscillate between small positive and negative values for the three most reliable functional combinations (BP86, B3LYP, and BLYP), and it remains unclear if either covalent anion is bound. AEAs with B3LYP/TZ2P++ single points are 0.19 (U), 0.16 (T), 0.07 (G), -0.02 (C), and -0.17 eV (A). Favorable comparisons are made to experimental estimates extrapolated from photoelectron spectra data for the complexes of the nucleobases with water. However, experimental values scaled from liquid-phase reduction potentials are shown to overestimate the AEAs by as much as 1.5 eV. Because the uracil and thymine covalent EAs are in energy ranges near those of their dipole-bound counterparts, preparation and precise experimental measurement of the thermodynamically stable covalent anions may prove challenging.     

Successive attachment of electrons to protonated Guanine: (G+H)* radicals and (G+H)- anions

The structures, energetics, and vibrational frequencies of nine hydrogenated 9H-keto-guanine radicals (G+H)(*) and closed-shell anions (G+H)(-) are predicted using the carefully calibrated (Chem. Rev. 2002, 102, 231) B3LYP density functional method in conjunction with a DZP++ basis set. These radical and anionic species come from consecutive electron attachment to the corresponding protonated (G+H)(+) cations in low pH environments. The (G+H)(+) cations are studied using the same level of theory. The proton affinity (PA) of guanine computed in this research (228.1 kcal/mol) is within 0.7 kcal/mol of the latest experiment value. The radicals range over 41 kcal/mol in relative energy, with radical r1, in which H is attached at the C8 site of guanine, having the lowest energy. The lowest energy anion is a2, derived by hydride ion attachment at the C2 site of guanine. No stable N2-site hydride should exist in the gas phase. Structure a9 was predicted to be dissociative in this research. The theoretical adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies were computed, with AEAs ranging from 0.07 to 3.12 eV for the nine radicals.     

Effects of electron attachment on C5'-O5' and C1'-N1 bond cleavages of pyrimidine nucleotides: A theoretical study

Sugar-base C(1')-N(1) and phosphate-sugar C(5')-O(5') bond breakings of 2'-deoxycytidine-5'-monophosphates (dCMP) and 2'-deoxythymidine-5'- monophosphates (dTMP) and their radical anions have been explored theoretically at the B3LYP/DZP++ level of theory. Calculations show that the low-energy electrons attachment to the pyrimidine nucleotides results in remarkable structural and chemical bonding changes. Predicted Gibbs free energies of reaction DeltaG for the C(5')-O(5') bond dissociation process of the radical anions are -14.6 and -11.5 kcal mol(-1), respectively, and such dissociation processes may be intrinsically spontaneous in the gas phase. Furthermore, the C(5')-O(5') bond cleavage processes of the anionic dCMP and dTMP were predicted to have activation energies of 6.9 and 8.0 kcal mol(-1) in the gas phase, respectively, much lower than the barriers for the C(1')-N(1) bond breaking process, showing that the C-O bond dissociation in DNA single strand breaks is a dominant process as observed experimentally.     

DNA nucleosides and their radical anions: molecular structures and electron affinities

The deoxyribonucleosides have been studied to determine the properties of combinations of 2-deoxyribose with each of the isolated DNA bases for both neutral and anionic species. We have used a carefully calibrated theoretical method [Chem. Rev. 2002, 102, 231], employing the B3LYP hybrid Hartree-Fock/DFT functional with the DZP++ basis set. Predictions are made of the geometric parameters, adiabatic electron affinities, charge distributions based on natural population analysis, and decomposition enthalpy for the neutral and anionic forms of the four 2'-deoxyribonucleosides in DNA: 2'-deoxyriboadenosine (dA), 2'-deoxyribocytidine (dC), 2'-deoxyriboguanosine (dG), and 2'-deoxyribothymidine (dT). Geometric changes in the anions show that the glycosidic bond exhibits little change with excess charge for the guanosine but significant shortening for the adenosine and for the pyrimidines. The zero-point corrected adiabatic electron affinities in eV for each of the 2'-deoxyribonucleosides are as follows: 0.06, dA; 0.09, dG; 0.33, dC; and 0.44, dT. These values are uniformly greater than those of the corresponding isolated bases (-0.28, A; -0.07, G; 0.03, C; 0.20, T) and the isolated 2-deoxyribose (-0.38) at the same level of theory. The vertical detachment energies of dT and dC are substantial, 0.72 and 0.94 eV, and these anions should be observable. A high VDE, 0.91 eV, is also found for dA but its anion is unlikely to be stable due to the small AEA of 0.06 eV. The high VDE reflects the fact that the molecular structures of the anions and the corresponding neutral species are quite different. Valence character is displayed for the SOMOs of dA, dC, and dT, while some component of diffuse character is visible in the SOMO of dG. Further analysis of electronic changes upon electron attachment include an examination of the NPA charges, which show that in the neutral 2'-deoxyribonucleosides the sum of NPA charges for every base is the same, -0.28 with the sum of 2-deoxyribose charges being positive, +0.28. In the anions, the trend in charge division varies based on the nature of the excess electron in the anions. Thermodynamically, the overall enthalpy change for the reaction of water with the neutral nucleosides to give bases and ribose is approximately zero. The analogous decomposition is exothermic by 8 to 11 kcal mol-1 for the anions, indicating possible challenges for anionic gas-phase nucleoside exploration in the presence of water.     

Photoinduced reductive repair of thymine glycol: implications for excess electron transfer through DNA containing modified bases

Photoinduced reduction of thymine glycol in oligodeoxynucleotides was investigated using either a reduced form of flavin adenine dinucleotide (FADH(-)) as an intermolecular electron donor or covalently linked phenothiazine (PTZ) as an intramolecular electron donor. Intermolecular electron donation from photoexcited flavin (FADH(-)) to free thymidine glycol generated thymidine in high yield, along with a small amount of 6-hydroxy-5,6-dihydrothymidine. In the case of photoreduction of 4-mer long single-stranded oligodeoxynucleotides containing thymine glycol by *FADH(-), the restoration yield of thymine was varied depending on the sequence of oligodeoxynucleotides. Time-resolved spectroscopic study on the photoreduction by laser-excited N,N-dimethylaniline (DMA) suggested elimination of a hydroxyl ion from the radical anion of thymidine glycol with a rate constant of approximately 10(4) s(-1) generates 6-hydroxy-5,6-dihydrothymidine (6-HOT(*)) as a key intermediate, followed by further reduction of 6-HOT(*) to thymidine or 6-hydroxy-5,6-dihydrothymdine (6-HOT). On the other hand, an excess electron injected into double-stranded DNA containing thymine glycol was not trapped at the lesion but was further transported along the duplex. Considering redox properties of the nucleobases and PTZ, competitive excess electron trapping at pyrimidine bases (thymine, T and cytosine, C) which leads to protonation of the radical anion (T(-)(*), C(-)(*)) or rapid back electron transfer to the radical cation of PTZ (PTZ(+)(*)), is presumably faster than elimination of the hydroxyl ion from the radical anion of thymine glycol in DNA.     

Do Photolyases Need To Provide Considerable Activation Energy for the Splitting of Cyclobutane Pyrimidine Dimer Radical Anions?

cis-syn Cyclobutane pyrimidine dimers, major UV-induced DNA lesions, are efficiently repaired by DNA photolyases. The key step of the repair reaction is a light-driven electron transfer from the FADH(-) cofactor to the dimer; the resulting radical anion splits spontaneously. Whether the splitting reaction requires considerable activation energy is still under dispute. Recent reports show that the splitting reaction of a dimer radical anion has a significant activation barrier (0.45 eV), and so photolyases have to provide considerable energy. However, these results contradict observations that cis-syn dimer radical anions split into monomers at -196 degrees C, and that the full process of DNA photoreactivation was fast (1.5-2 ns). To investigate the activation energies of dimer radical anions, three model compounds 1-3 were prepared. These include a covalently linked cyclobutane thymine dimer and a tryptophan residue (1) or a flavin unit (3), and the covalently linked uracil dimer and tryptophan (2). Their properties of photosensitised splitting of the dimer units by tryptophan or flavin unit were investigated over a large temperature range, -196 to 70 degrees C. The activation energies were obtained from the temperature dependency of splitting reactions for 1 and 2, 1.9 kJ mol(-1) and 0.9 kJ mol(-1) for the thymine and uracil dimer radical anions, respectively. These values are much lower than that obtained for E. coli photolyase (0.45 eV), and are surmountable at -196 degrees C. The activation energies provide support for previous observations that repair efficiencies for uracil dimers are higher than thymine dimers, both in enzymatic and model systems. The mechanisms of highly efficient enzymatic DNA repair are discussed.     

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

Calibration and applications of the ΔMP2 method for calculating core electron binding energies

We calibrated a method for the evaluation of core electron binding energies, based on the energy differences between the cation and neutral molecule evaluated at the level of M?ller-Plesset perturbation theory. The central feature of the method is the use of a mixed basis set: a large all-electron basis set is used for the atom whose core electron is removed, while the model core potential basis set is employed for all remaining atoms. Calibration was carried out for 55 molecules and 114 binding energies of 1s core electrons for the atoms C, N, O, and F. The average absolute deviation for all the core electron binding energies is 0.163 eV. The method was applied to the calculation of the core electron binding energies of five nucleic acid bases (uracil, adenine, cytosine, guanine, and thymine) and several of their low-energy tautomers.     

Extension of the PDDG/PM3 and PDDG/MNDO semiempirical molecular orbital methods to the halogens

The new semiempirical methods, PDDG/PM3 and PDDG/MNDO, have been parameterized for halogens. For comparison, the original MNDO and PM3 were also reoptimized for the halogens using the same training set; these modified methods are referred to as MNDO' and PM3'. For 442 halogen-containing molecules, the smallest mean absolute error (MAE) in heats of formation is obtained with PDDG/PM3 (5.6 kcal/mol), followed by PM3' (6.1 kcal/mol), PDDG/MNDO (6.6 kcal/mol), PM3 (8.1 kcal/mol), MNDO' (8.5 kcal/mol), AM1 (11.1 kcal/mol), and MNDO (14.0 kcal/mol). For normal-valent halogen-containing molecules, the PDDG methods also provide improved heats of formation over MNDO/d. Hypervalent compounds were not included in the training set and improvements over the standard NDDO methods with sp basis sets were not obtained. For small haloalkanes, the PDDG methods yield more accurate heats of formation than are obtained from density functional theory (DFT) with the B3LYP and B3PW91 functionals using large basis sets. PDDG/PM3 and PM3' also give improved binding energies over the standard NDDO methods for complexes involving halide anions, and they are competitive with B3LYP/6-311++G(d,p) results including thermal corrections. Among the semiempirical methods studied, PDDG/PM3 also generates the best agreement with high-level ab initio G2 and CCSD(T) intrinsic activation energies for S(N)2 reactions involving methyl halides and halide anions. Finally, the MAEs in ionization potentials, dipole moments, and molecular geometries show that the parameter sets for the PDDG and reoptimized NDDO methods reduce the MAEs in heats of formation without compromising the other important QM observables.     

The 2'-deoxyadenosine-5'-phosphate anion, the analogous radical, and the different hydrogen-abstracted radical anions: molecular structures and effects on DNA damage

The 2'-deoxyadenosine-5'-phosphate (5'-dAMP) anion and its related radicals have been studied by reliably calibrated theoretical approaches. This study reveals important physical characteristics of 5'-dAMP radical related processes. One-electron oxidation of the 5'-dAMP anion is found on both the phosphoryl group and the adenine base with electron detachment energies close to that of phosphate. Partial removal of electron density from the adenine fragment leads to an extended pi system which includes the amine group of the adenine. Although the radical-centered carbon increases the extent of bonding with its adjacent atoms, it usually weakens the chemical bonds between the atoms at the alpha- and beta-positions. This tendency should be important in predicting the reactivity of the sugar-based radicals. The overall stability sequence of the H-abstracted 5'-dAMP anionic radicals is consistent with the analogous results for the H-abstracted neutral radicals of the adenosine nucleoside: aliphatic radicals > aromatic radicals. The negatively charged phosphoryl group attached to atom C(5)' of the ribose does not change this energetic sequence. All the H-abstraction produced 5'-dAMP radical anions are distonic radical anions. Studies have shown that the charge-radical-separating feature of the distonic radical anions is biologically relevant. This result should be important in understanding the reactive properties of these H-abstraction-produced anion radicals.