Abstract
Quantum mechanical (QM) and hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations of a recently reported dinuclear mercury(II)-mediated base pair were performed aiming to analyse its intramolecular bonding pattern, its stability, and to obtain clues on the mechanism of the incorporation of mercury(II) into the DNA. The dynamic distance constraint was employed to find initial structures, control the dissociation process in an unbiased fashion and to determine the free energy required. A strong influence of the exocyclic carbonyl or amino groups of neighbouring base pairs on both the bonding pattern and the mechanism of incorporation was observed. During the dissociation simulation, an amino group of an adenine moiety of the adjacent base pair acts as a turnstile to rotate the mercury(II) ion out of the DNA core region. The calculations provide an important insight into the mechanism of formation of this dinuclear metal-mediated base pair and indicate that the exact location of a transition metal ion in a metal-mediated base pair may be more ambiguous than derived from simple model building.