On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level.

The geometries and interaction energies of stacked and hydrogen-bonded uracil dimers and a stacked adeninecdots, three dots, centeredthymine pair were studied by means of high-level quantum chemical calculations. Specifically, standard as well as counterpoise-corrected optimizations were performed at second-order Moller-Plesset (MP2) and coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)] levels with various basis sets up to the complete basis set limit. The results can be summarized as follows: (i) standard geometry optimization with small basis set (e.g., 6-31G(*)) provides fairly reasonable intermolecular separation; (ii) geometry optimization with extended basis sets at the MP2 level underestimates the intermolecular distances compared to the reference CCSD(T) results, whereas the MP2/cc-pVTZ counterpoise-corrected optimization agrees well with the reference geometries and, therefore, is recommended as a next step for improving MP2/cc-pVTZ geometries; (iii) the stabilization energy of stacked nucleic acids base pairs depends considerably on the method used for geometry optimization, so the use of reliable geometries, such as counterpoise-corrected MP2/cc-pVTZ ones, is recommended; (iv) the density functional theory methods fail completely in locating the energy minima for stacked structures and when the geometries from MP2 calculations are used, the resulting stabilization energies are strongly underestimated; (v) the self-consistent charges-density functional tight binding method, with inclusion of the empirical dispersion energy, accurately reproduces interaction energies and geometries of dispersion-bonded (stacked) complexes; this method can thus be recommended for prescanning the potential energy surfaces of van der Waals complexes.

[1]  Jirí Cerný,et al.  The X3LYP extended density functional accurately describes H-bonding but fails completely for stacking. , 2005, Physical chemistry chemical physics : PCCP.

[2]  Pavel Hobza,et al.  Unexpectedly strong energy stabilization inside the hydrophobic core of small protein rubredoxin mediated by aromatic residues: correlated ab initio quantum chemical calculations. , 2005, Journal of the American Chemical Society.

[3]  P. Wawrzyniak,et al.  Theoretical study of the complex between formic acid and argon , 2004 .

[4]  Pavel Hobza,et al.  True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment. , 2003, Journal of the American Chemical Society.

[5]  Pavel Hobza,et al.  On the convergence of the (ΔECCSD(T)−ΔEMP2) term for complexes with multiple H-bonds , 2002 .

[6]  Pavel Hobza,et al.  Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations. , 2002, Journal of the American Chemical Society.

[7]  K. M. Ervin,et al.  Dynamics of the gas-phase reactions of chloride ion with fluoromethane: high excess translational activation energy for an endothermic S(N)2 reaction. , 2002, Journal of the American Chemical Society.

[8]  Thomas Frauenheim,et al.  Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment , 2001 .

[9]  P. Hobza,et al.  Interaction Energies of Hydrogen-Bonded Formamide Dimer, Formamidine Dimer, and Selected DNA Base Pairs Obtained with Large Basis Sets of Atomic Orbitals , 2000 .

[10]  P Hobza,et al.  Noncovalent interactions: a challenge for experiment and theory. , 2000, Chemical reviews.

[11]  P Hobza,et al.  Structure, energetics, and dynamics of the nucleic Acid base pairs: nonempirical ab initio calculations. , 1999, Chemical reviews.

[12]  Trygve Helgaker,et al.  Basis-set convergence of the energy in molecular Hartree–Fock calculations , 1999 .

[13]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[14]  Donald G. Truhlar,et al.  Basis-set extrapolation , 1998 .

[15]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

[16]  Miquel Duran,et al.  How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? , 1996 .

[17]  Pavel Hobza,et al.  Density functional theory and molecular clusters , 1995, J. Comput. Chem..

[18]  Thomas A. Halgren,et al.  The representation of van der Waals (vdW) interactions in molecular mechanics force fields: potential form, combination rules, and vdW parameters , 1992 .

[19]  David Feller,et al.  Application of systematic sequences of wave functions to the water dimer , 1992 .

[20]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[21]  Pavel Hobza,et al.  Intermolecular interactions between medium-sized systems. Nonempirical and empirical calculations of interaction energies. Successes and failures , 1988 .

[22]  L. Kroon-Batenburg,et al.  The use of a moment-optimized DZP basis set for describing the interaction in the water dimer , 1985 .

[23]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[24]  Pavel Hobza,et al.  Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies , 2003 .

[25]  Pavel Hobza,et al.  Reliable theoretical treatment of molecular clusters: Counterpoise-corrected potential energy surface and anharmonic vibrational frequencies of the water dimer , 1999 .