Docking studies on DNA-ligand interactions: building and application of a protocol to identify the binding mode

Despite DNA being an important target for several drugs, most of the docking programs are validated only for proteins and their ligands. In this paper, we used AutoDock 4.0 to perform self-dockings and cross dockings between two DNA ligands (a minor groove binder and an intercalator) and four distinct receptors: 1) crystallographic DNA without intercalation gap; 2) crystallographic DNA with intercalation gap; 3) canonical B-DNA; and 4) modified B-DNA with intercalation gap. Besides being efficient in self-dockings, AutoDock is capable of correctly identifying two of the main DNA binding modes with the condition that the target possesses an artificial intercalation gap. Therefore, we suggest a default protocol to identify DNA binding modes which uses a modified canonical DNA (with gap) as receptor. This protocol was applied to dock two different Troger bases to DNA and the predicted binding modes agree with those suggested, yet not established, by experimental data. We also applied the protocol to dock aflatoxin B(1) exo-8,9-epoxide, and the results are in complete agreement with experimental data from the literature. We propose that this approach can be used to investigate other ligands whose binding mode to DNA remains unknown, yielding a suitable starting point for further theoretical studies such as molecular dynamics simulations.