Modeling the Conformational Preference of the Carbon-Bonded Covalent Adduct Formed upon Exposure of 2′-Deoxyguanosine to Ochratoxin A

Abstract

The conformational flexibility of the C8-linked guanine adduct formed from attachment of ochratoxin A (OTA) was analyzed using a systematic computational approach and models ranging from the nucleobase to the adducted DNA helix. A focus was placed on the influence of the C8-modification of 2'-deoxyguanosine (dG) on the preferred relative arrangement of the nucleobase and the C8-substituent and, more importantly, the anti/syn conformational preference with respect to the glycosidic bond. Although OTA is twisted with respect to the base in the nucleobase model, addition of the deoxyribose sugar induces a further twist and restricts rotation about the C-C linkage due to close contacts between OTA and the sugar. The nucleoside model preferentially adpots a syn orientation (by 10-20 kJ mol(-1) depending on the OTA conformation) due to the presence of an O5'-H···N3 interaction. However, when this hydrogen bond is eliminated, which better mimics the DNA environment, a small (<5 kJ mol(-1)) anti/syn energy difference is predicted. Inclusion of the 5'-monophosphate group leads to an up to 20 kJ mol(-1) preference for the syn (nucleotide) conformation due to stabilizing base-phosphate interactions involving the amino group of guanine. Nevertheless, MD simulations and free energy analysis predict that both syn- and anti-conformations of OTB-dG are equally stable in helices when paired opposite cytosine. These results indicate that the adduct will likely adopt a syn conformation in an isolated nucleoside and nucleotide, while a mixture of syn and anti conformations will be observed in DNA duplexes. Since the syn conformation could stabilize base mismatches upon DNA replication or Z-DNA structures with varied biological outcomes, future computational and experimental work should elucidate the consequences of the conformational preference of this potentially harmful DNA lesion.