RI-MP2 calculations with extended basis sets—a promising tool for study of H-bonded and stacked DNA base pairs

Abstract

The interaction energies of 9-methyladenine···1-methylthymine H-bonded and 9-methyladenine···9-methylguanine stacked pairs were evaluated at the MP2 and resolution of the identity MP2 (RI-MP2) levels. The interaction energies are almost identical for both methods. The RI-MP2 method is about one order of magnitude faster than the exact MP2 method and, therefore, this method is well suited for the study of large models of biological importance. The basis set dependence of both base pairs was studied and reasonable relative energy characteristics were obtained only if at least the valence double-ζ-polarisation SVP basis set with modified (diffuse) polarisation functions on non-hydrogen and hydrogen atoms was used. However, reliable absolute stabilisation energies of both base pairs were only obtained if at least the aug-SVP basis set is applied. Stack interaction energies and three- and four-body energy terms for selected base-pair steps of B-DNA were calculated for the first time using the RI-MP2 method and the aug-SVP basis set.