Abstract
We have used molecular dynamics simulations to study the structure and dynamics of a range of DNA duplexes containing the 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapydG) lesion that can result from oxidative damage at guanine. Compared to the corresponding undamaged DNA duplexes, FapydG-containing duplexes show little gross structural changes-the damaged base remains stacked in to the DNA double helix and retains hydrogen bonds to its cytosine partner. However, the experimentally observed reduction in DNA stability that accompanies lesion formation can be explained by a careful energetic analysis of the simulation data. Irrespective of the nature of the base pairs on either side of the lesion site, conversion of a guanine to a FapydG base results in increased dynamical flexibility in the base (but not in the DNA as a whole) that significantly weakens its hydrogen-bonding interactions. Surprisingly, the stacking interactions with its neighbours are not greatly altered. The formamido group adopts a non-planar conformation that can interact significantly and in a sequence-dependent manner with its 3'-neighbour. We conclude that the recognition of FapydG lesions by the repair protein formamidopyrimidine-DNA glycosylase probably does not involve the protein capturing an already-extrahelical FapydG base, but rather it relies on detecting alterations to the DNA structure and flexibility created by the lesion site.