Isomers of the uracil dimer:: An ab initio benchmark study

被引:26
作者
Frey, Jann A. [1 ]
Mueller, Andreas [1 ]
Losada, Martin [1 ]
Leutwyler, Samuel [1 ]
机构
[1] Univ Bern, Dept Chem & Biochem, CH-3012 Bern, Switzerland
关键词
D O I
10.1021/jp0683162
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
Benchmark ab initio calculations at the correlated level are reported for ten isomers of the uracil dimer (U center dot U): six are doubly N-H center dot center dot center dot O hydrogen bonded, three have a C-H center dot center dot center dot O and an N-H center dot center dot center dot O hydrogen bond, and one is doubly C-H center dot center dot center dot O hydrogen bonded. Their structures were optimized at the correlated level by using second-order Moller-Plesset perturbation theory (MP2), resolution of identity MP2 (RIMP2), and the binding energies D-e calculated with the coupled-cluster method with singles, doubles, and iterative triples, CCSD(T). The MP2 and RIMP2 binding energies D-e are extrapolated to the complete basis set (CBS) limit, using the aug-cc-pVXZ (X = D, T, Q) basis sets, giving binding energies accurate to +/- 0.07 kcal/mol. With one exception, the correlation energy contributions at the CCSD(T) level increase the binding energies, although the changes are small, +0.03 to -0.27 kcal/mol (or 0.1% to 2.2%). The most stable isomer is the doubly N1-H center dot center dot center dot O hydrogen-bonded HB4 form, with D-e[CCSD(T)]= -19.04 kcal/mol. The biologically relevant HB2 dimer has D-e[CCSD(T)] = -12.64 kcal/mol, and the HB5 dimer that constitutes the main structural motif of the uracil crystal has -13.20 kcal/mol. The "Calcutta" dimer, which occurs in an RNA hexamer, is among the weakest isomers, with D-e[CCSD(T)] = -9.81 kcal/mol. We compare to the binding energies calculated with the B3LYP, PBE, and PW91 density functionals; the PW91/6-311++G(d,p) binding energies agree with the CBS benchmark values, to within < 2%. A useful single-molecule descriptor for the strengths of the individual hydrogen bonds can be derived from the gas-phase acidity Delta E-0(A-H) of the N-H or C-H donor groups and the gas-phase proton affinity PA(0)(B) of the CO groups of the uracil monomer. The calculated hydrogen bond energies D-e(infinity) correlate well with the difference between gas-phase acidity and basicity, Delta E-0(A-H) - PA(0)(B).
引用
收藏
页码:3534 / 3542
页数:9
相关论文
共 40 条
[1]   ELECTRONIC-STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS - THE PROGRAM SYSTEM TURBOMOLE [J].
AHLRICHS, R ;
BAR, M ;
HASER, M ;
HORN, H ;
KOLMEL, C .
CHEMICAL PHYSICS LETTERS, 1989, 162 (03) :165-169
[2]  
BAEYENDS KJ, 1995, NAT STRUCT BIOL, V2, P52
[3]   CALCULATION OF SMALL MOLECULAR INTERACTIONS BY DIFFERENCES OF SEPARATE TOTAL ENERGIES - SOME PROCEDURES WITH REDUCED ERRORS [J].
BOYS, SF ;
BERNARDI, F .
MOLECULAR PHYSICS, 1970, 19 (04) :553-&
[4]   THE ELECTRONIC-SPECTRA OF THE PYRIMIDINE-BASES URACIL AND THYMINE IN A SUPERSONIC MOLECULAR-BEAM [J].
BRADY, BB ;
PETEANU, LA ;
LEVY, DH .
CHEMICAL PHYSICS LETTERS, 1988, 147 (06) :538-543
[5]  
CARRABINE JA, 1971, BIOCHEMISTRY-US, V10, P292
[6]   Functional insights from the structure of the 30S ribosomal subunit and its interactions with antibiotics [J].
Carter, AP ;
Clemons, WM ;
Brodersen, DE ;
Morgan-Warren, RJ ;
Wimberly, BT ;
Ramakrishnan, V .
NATURE, 2000, 407 (6802) :340-348
[7]  
Chandra AK, 1999, J PHYS CHEM A, V103, P8853
[8]   Theoretical study of the acidity and basicity of the cytosine tautomers and their 1:1 complexes with water [J].
Chandra, AK ;
Michalska, D ;
Wysokinsky, R ;
Zeegers-Huyskens, T .
JOURNAL OF PHYSICAL CHEMISTRY A, 2004, 108 (44) :9593-9600
[9]   CODON-ANTICODON PAIRING - WOBBLE HYPOTHESIS [J].
CRICK, FHC .
JOURNAL OF MOLECULAR BIOLOGY, 1966, 19 (02) :548-&
[10]  
DABKOWSKA I, 2004, J CHEM PHYS, V122