An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.

Institute of Biophysics Academy of Sciences of the Czech Republic and National Center for Biomolecular Research Královopolská 135 612 65 Brno Czech Republic

References provided by Crossref.org

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

Molecular dynamics simulations and their application to four-stranded DNA

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E