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.
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
submitted on 03.10.2019 and posted on 07.10.2019by 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.