Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

Institute of Biophysics Academy of Sciences of the Czech Republic Kralovopolska 135 612 65 Brno Czech Republic

References provided by Crossref.org

Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions

Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies

QM/MM studies of hairpin ribozyme self-cleavage suggest the feasibility of multiple competing reaction mechanisms

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

Stretched DNA investigated using molecular-dynamics and quantum-mechanical calculations

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Molecular dynamics simulations and their application to four-stranded DNA

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

Molecular dynamics simulations of sarcin-ricin rRNA motif