New insights into the structure of abasic DNA from molecular dynamics simulations

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5′-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson–Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20–30° at the abasic site. The Watson–Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.