These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
2 files

Importance of Model Size in Quantum Mechanical Studies of DNA Intercalation

submitted on 03.10.2019, 18:56 and posted on 07.10.2019, 21:39 by Drew P. Harding, Laura J. Kingsley, Glen Spraggon, Steven Wheeler
There is currently a dearth of effective computational tools to design nucleobase-targeting small molecules and molecular mechanics force-fields for nucleobases lag behind their protein-focused counterparts. While quantum chemical methods can provide reliable interaction energies for small molecule-nucleobase interactions, these come at a steep computational cost. As a first step toward refining available tools for predicting small molecule-nucleobase interactions, we assessed the convergence of DFT-computed interaction energies with increasing binding site model size. We find that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in non-covalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.




Email Address of Submitting Author


University of Georgia



ORCID For Submitting Author


Declaration of Conflict of Interest

No conflicts of interest


Read the published paper

in Journal of Computational Chemistry