Structural insights by molecular dynamics simulations into differential repair efficiency for ethano-A versus etheno-A adducts by the human alkylpurine-DNA N-glycosylase

1,N-6-ethenoadenine adducts (epsilonA) are formed by known environmental carcinogens and found to be removed by human alkylpurine-DNA N-glycosylase (APNG). 1,N-6-ethanoadenine (EA) adducts differ from epsilonA by change of a double bond to a single bond in the 5-member exocyclic ring and are formed by chloroethyl nitrosoureas, which are used in cancer therapy. In this work, using purified recombinant human APNG, we show that EA is a substrate for the enzyme. However, the excision efficiency of EA was 65-fold lower than that of epsilonA. Molecular dynamics simulation produced similar structural motifs for epsilonA and EA when incorporated into a DNA duplex, suggesting that there are no specific conformational features in the DNA duplex which can account for the differences in repair efficiency. However, when EA was modeled into the APNG active site, based on the APNG/epsilonA-DNA crystallographic coordinates, in structures produced by 2 ns molecular dynamics simulation, we observed weakening in the stacking interaction between EA and aromatic side chains of the key amino acids in the active site. In contrast, the planar epsilonA is better stacked at the enzyme active site. We propose that the observed destabilization of the EA adduct at the active site, such as reduced stacking interactions, could account for the biochemically observed weaker recognition of EA by APNG as compared to epsilonA.