A second life for the crystallographic structure of Berenil-dodecanucleotide complex: a computational revisitation thirty years after its publication

This study revisits the pioneering work of Professor Neidle, and co-workers, on the crystal structure of complexes formed between groove binders and DNA sequences. The original research revealed a DNA-ligand complex consisting of a dodecanucleotide bound with Berenil [1,3-bis(4′-amidinophenyl)-triazene] an anti-trypanocidal drug. This article aims to delve deeper into the structural dynamics of this system, showcasing the role played by water molecules in stabilizing the interaction between the ligand and the DNA. With this work, we reevaluate the findings from the original crystallographic study by employing modern molecular dynamics techniques, including Supervised Molecular Dynamics (SuMD) for generating binding trajectories, Thermal Titration Molecular Dynamics for assessing unbinding events, and AquaMMapS to identify regions occupied by stationary water molecules. The study addresses a minor and a major groove binding mode and assesses their strength and specificity using TTMD simulations, generating unbinding trajectories. This comprehensive approach integrates the understanding of the interaction of this DNA-ligand complex, which originated with the valuable work of Professor Neidle, resulting in an in-depth insight into the pivotal role of water molecules with this DNA, a behavior detected and extendable even to other nucleic acid complexes.