学术文献库

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.