Density functional theory was used to study the thermodynamics and kinetics for the glycosidic bond cleavage in deoxyuridine. Two reaction pathways were characterized for the unimolecular decomposition in vacuo. However, these processes are associated with large reaction barriers and highly endothermic reaction energies, which is in agreement with experiments that suggest a (water) nucleophile is required for the nonenzymatic glycosidic bond cleavage. Two (S(N)1 and S(N)2) reaction pathways were characterized for direct hydrolysis of the glycosidic bond by a single water molecule; however, both pathways also involve very large barriers. Activation of the water nucleophile via partial proton abstraction steadily decreases the barrier and leads to a more exothermic reaction energy as the proton affinity of the molecule interacting with water increases. Indeed, our data suggests that the barrier heights and reaction energies range from that for hydrolysis by water to that for hydrolysis by the hydroxyl anion, which represents the extreme of (full) water activation (deprotonation). Hydrogen bonds between small molecules (hydrogen fluoride, water, or ammonia) and the nucleobase were found to further decrease the barrier and overall reaction energy but not to the extent that the same hydrogen-bonding interactions increase the acidity of the nucleobase. Our results suggest that the nature of the nucleophile plays a more important role in reducing the barrier to glycosidic bond cleavage than the nature of the small molecule bound, and models with more than one hydrogen fluoride molecule interacting with the nucleobase provide further support for this conclusion. Our results lead to a greater fundamental understanding of the effects of the nucleophile, activation of the nucleophile, and interactions with the nucleobase for this important biological reaction.

The structural and spectral properties of (ortho and para) C8-aryl-purine adducts formed from carbon attachment by phenolic toxins were investigated through DFT calculations and UV-vis absorbance and emission studies. The global minima of both the deoxyadenosine (dA) and deoxyguanosine (dG) adducts adopted a syn conformation about the glycosidic bond due to the presence of an O5'-H...N3 hydrogen bond, where the anti minima are 20-30 kJ mol-1 higher in energy. While the nucleobase adducts are planar, the presence of the deoxyribose sugar induces a twist about the carbon-carbon bond connecting the phenol and nucleobase rings. ortho-Phenolic adducts are less twisted than the corresponding para adducts due to stabilization provided by an intramolecular O-H...N7 bond. Solvation calculations, in combination with UV-vis studies, demonstrate that the structural preference is solvent dependent, where solvents with hydrogen-bonding abilities disrupt the intramolecular O-H...N7 hydrogen bond such that a greater degree of twist is observed, and less polar solvents stabilize the planar structure. Indeed, the ratio of twisted to planar conformers is estimated to be as large as 50:50 in some aprotic solvents. Thus, the combined experimental and computational approach has provided a greater understanding of the structure of the ortho- and para-dA and dG C-bonded phenoxyl adducts as the first step to understanding the biological consequences of this form of DNA damage.

Previous computational work suggests that isolated C8-phenoxyl-2'-deoxyguanosine nucleoside adducts preferentially adopt a syn orientation about the glycosidic bond, which is the first step in the mechanism by which many bulky C8 adducts exert their mutagenic effects. Since it is not clear whether these results can be directly extrapolated to the preferred conformation in DNA helices, approaches that more accurately reflect the physiological environment were used in the present study to understand the anti/syn preference of the ortho and para C8-phenoxyl-2'-deoxyguanosine adducts. Using nucleoside models and methods (B3LYP) similar to those previously implemented, we determine that the syn conformer is less stable than previously predicted when geometries relevant to B-DNA are considered. This indicates that the conformational energy trend is model dependent and stresses the importance of considering models that better mimic the DNA environment when determining the conformational preference of damaged bases. Therefore, we expanded our computational model to include the 5'-monophosphate group. Since the correct anti/syn energy trend for 2'-deoxyguanosine (dG) 5'-monophosphate has only been found using very specific computational models and prior knowledge of the biologically relevant nucleotide conformation, which is unavailable for most damaged systems, we initially benchmark our computational approach by studying the natural nucleotide. Despite the wide use of gas-phase optimizations in the current literature, only through the implementation of solvation-phase optimizations, as well as the use of a counterion model for the phosphate backbone, is the correct anti/syn energy trend predicted. Indeed, this is the first time in the literature that a biologically relevant syn structure is characterized for dG using methods suitable for studying bulky DNA adducts. Subsequently, our newly identified approach for DNA lesions was used to study C8-phenoxyl DNA adducts. In contrast to previously published results, we predict that the ortho and para adducts may adopt both the anti and syn conformations in DNA helices. These results have implications for the base-pairing properties and mutagenicity of these adducts, which must be further considered in future work.