Abstract:The computational study of DNA and its interaction with ligands is a highly relevant area of research, with significant consequences for developing new therapeutic strategies. However, the computational description of such large and complex systems requires considering interactions of different types. All these considerations imply a real challenge for computational chemistry. Using quantum methods for the entire system requires significant computational resources, with improvements in parallelization and optimization of theoretical strategies. Computational methods, such as LS-DFT and DLPNO-CCSD(T), may allow performing ab initio QM calculations, including explicitly the electronic structure for large biological systems, at a reasonable computing time. In this work, we study the interaction of small molecules and cations with DNA (duplex-DNA and G-quadruplexes), comparing different computational methods: a linear-scaling DFT (LS-DFT) at LMKLL/DZDP level of theory, semi-empirical methods (PM6-DH2 and PM7), mixed QM/MM, and DLPNO-CCSD(T). Our goal is to demonstrate the adequacy of LS-DFT to treat the different types of interactions present in DNA-dependent systems. We show that LMKLL/DZDP using SIESTA can yield very accurate geometries and energetics in all the different systems considered in this work: duplex DNA (dDNA), phenanthroline intercalating dDNA, G-quadruplexes, and Metal-G-tetrads considering alkaline metals of different sizes. As far as we know, this is the first time that full G-quadruplex geometry optimizations have been carried out using a DFT method thanks to its linear-scaling capabilities. Moreover, we show that LS-DFT provides high-quality structures, and some semi-empirical Hamiltonian can also yield suitable geometries. However, DLPNO-CCSD(T) and LS-DFT are the only methods that accurately describe interaction energies for all the systems considered in our study.

Enhanced sampling molecular dynamics simulations correctly predict the diverse activities of a series of stiff-stilbene G-quadruplex DNA ligands

Ligand Binding to Dynamically Populated G‐Quadruplex DNA

Stabilizing Parallel G-Quadruplex Dna by A New Class of Ligands: Two Non-Planar Alkaloids Through Interaction in Lateral Grooves

Single-molecule imaging reveals a common mechanism shared by G-quadruplex–resolving helicases

Semi-empirical and Linear-Scaling DFT Methods to Characterize duplex DNA and G-quadruplexes in Presence of Interacting Small Molecules

Molecular dynamics simulations of G-quadruplexes: The basic principles and their application to folding and ligand binding

Molecular dynamics simulations reveal the parallel stranded d(GGGA) 3 GGG DNA quadruplex folds via multiple paths from a coil-like ensemble

G4-Ligand-Conjugated Oligonucleotides Mediate Selective Binding and Stabilization of Individual G4 DNA Structures

Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad—a common structural unit of G-quadruplex DNA

Evaluating Molecular Docking Software for Small Molecule Binding to G-Quadruplex DNA

Complex Biophysical and Computational Analyses of G‐Quadruplex Ligands: The Porphyrin Stacks Back

Direct genome-wide identification of G-quadruplex structures by whole-genome resequencing

Probing G-quadruplex-ligand binding using DNA intrinsic fluorescence

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism

Single-molecule visualization of DNA G-quadruplex formation in live cells

Single-Molecule Analysis of the Improved Variants of the G-Quadruplex Recognition Protein G4P

Optimized Virtual Screening Workflow for the Identification of Novel G-Quadruplex Ligands

Unfolding and Conformational Variations of Thrombin-Binding DNA Aptamers: Synthesis, Circular Dichroism and Molecular Dynamics Simulations.

Structural Dynamics Of Human Telomeric G-Quadruplex Loops Studied By Molecular Dynamics Simulations

Quantitative analysis and prediction of G-quadruplex forming sequences in double-stranded DNA

The Transition from Unfolded to Folded G-Quadruplex DNA Analyzed and Interpreted by Two-Dimensional Infrared Spectroscopy