Benchmark RI-MP2 database of nucleic acid base trimers: performance of different density functional models for prediction of structures and binding energies
- 9 July 2007
- journal article
- research article
- Published by Royal Society of Chemistry (RSC) in Physical Chemistry Chemical Physics
- Vol. 9 (36) , 5000-5008
Abstract
A new database of nucleic acid base trimers has been developed that includes 141 geometries and stabilization energies obtained at the RI-MP2 level of theory with the TZVPP basis set. Compared to previously compiled biologically oriented databases, this new construct includes considerably more complicated structures; the various intermolecular interactions in the trimers are quite heterogeneous and in particular include simultaneous hydrogen bonding and stacking interactions, which is similar to the situation in actual biopolymers. Validation against these benchmark data is therefore a more demanding task for approximate models, since correct descriptions of all energy terms are unlikely to be accomplished by fortuitous cancellations of systematic errors. The density functionals TPSS (both with and without an empirical dispersion term), PWB6K, M05-2X, and BH&H, and the self-consistent charge density functional tight binding method augmented with an empirical dispersion term (SCC-DFTB-D) were assessed for their abilities accurately to compute structures and energies. The best reproduction of the BSSE corrected RI-MP2 stabilization energies was achieved by the TPSS functional (TZVPP basis set) combined with empirical dispersion; removal of the dispersion correction leads to significantly degraded performance. The M05-2X and PWB6K functionals performed very well in reproducing the RI-MP2 geometries, but showed a systematic moderate underestimation of the magnitude of base stacking interactions. The SCC-DFTB-D method predicts geometries in fair agreement with RI-MP2; given its computational efficiency it represents a good option for initial scanning of analogous biopolymeric potential energy surfaces. BH&H gives geometries of comparable quality to the other functionals but significantly overestimates interaction energies other than stacking.Keywords
This publication has 51 references indexed in Scilit:
- Van der Waals Interactions in Aromatic Systems: Structure and Energetics of Dimers and Trimers of PyridineChemphyschem, 2005
- The X3LYP extended density functional accurately describes H-bonding but fails completely for stackingPhysical Chemistry Chemical Physics, 2005
- Helium dimer dispersion forces and correlation potentials in density functional theoryThe Journal of Chemical Physics, 2002
- RI-MP2 calculations with extended basis sets—a promising tool for study of H-bonded and stacked DNA base pairsPhysical Chemistry Chemical Physics, 2001
- Ab Initio Calculation of Nonbonded Interactions: Are We There Yet?The Journal of Physical Chemistry A, 2000
- The Nature of the Hydrogen Bond in DNA Base Pairs: The Role of Charge Transfer and Resonance AssistanceChemistry – A European Journal, 1999
- Charge Transfer and Environment Effects Responsible for Characteristics of DNA Base PairingAngewandte Chemie International Edition, 1999
- RI-MP2: first derivatives and global consistencyTheoretical Chemistry Accounts, 1997
- Density functional theory and molecular clustersJournal of Computational Chemistry, 1995
- Use of approximate integrals in ab initio theory. An application in MP2 energy calculationsChemical Physics Letters, 1993