Cytosine Unstacking and Strand Slippage at an Insertion–Deletion Mutation Sequence in an Overhang-Containing DNA Duplex

[Image: see text] Base unstacking in template strands, when accompanied by strand slippage, can result in deletion mutations during strand extension by nucleic acid polymerases. In a GCCC mutation hot-spot sequence, which was previously identified to have a 50% probability of causing such mutations during DNA replication by a Y-family polymerase, a single-base deletion mutation could result from such unstacking of any one of its three template cytosines. In this study, the intrinsic energetic differences in unstacking among these three cytosines in a solvated DNA duplex overhang model were examined using umbrella sampling molecular dynamics simulations. The free energy profiles obtained show that cytosine unstacking grows progressively more unfavorable as one moves inside the duplex from the 5′-end of the overhang template strand. Spontaneous strand slippage occurs in response to such base unstacking in the direction of both the major and minor grooves for all three cytosines. Unrestrained simulations run from three distinct strand-slipped states and one non-strand-slipped state suggest that a more duplexlike environment can help stabilize strand slippage. The possible underlying reasons and biological implications of these observations are discussed in the context of nucleic acid replication active site dynamics.