Conformational Flexibility of C8-Phenoxyl-2′-deoxyguanosine Nucleotide Adducts

• Published in: The journal of physical chemistry. B 2010-04, Vol.114 (12), p.4373-4382
• Main Authors: Millen, Andrea L, Manderville, Richard A, Wetmore, Stacey D
• Format: Article
• Language: English
• Subjects: B: Biophysical Chemistry; Deoxyguanosine - chemistry; DNA - chemistry; Nucleic Acid Conformation; Protein Conformation
• ISSN: 1520-6106; 1520-5207
• DOI: 10.1021/jp911993f

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.