Biopolym. Cell. 2011; 27(1):74-81.
Molecular Biophysics
Analysis of 2-deoxy-D-ribofuranose molecule conformational capacity with the quantum-mechanical density functional method
1Nikolaienko T. Yu., 1Bulavin L. A., 2, 3Hovorun D. M.
  1. Taras Shevchenko National University of Kyiv
    64/13, Volodymyrska Str., Kyiv, Ukraine, 01601
  2. Institute of Molecular Biology and Genetics, NAS of Ukraine
    150, Akademika Zabolotnoho Str., Kyiv, Ukraine, 03680
  3. Institute of High Technologies,
    Taras Shevchenko National University of Kyiv
    2, korp.5, Pr. Akademika Hlushkova, Kyiv, Ukraine, 03022

Abstract

Aim. To reveal the exhaustive conformational family for α- and β-anomeric forms of the 2-deoxy-D-ribofuranose molecule and to elucidate whether this molecule can be considered as the representative DNA sugar model. Methods. The molecule conformations were described with four torsion angles (β, γ, ε and χ') and the sugar pseudorotation angle P. In 810 initial geometries for each of the anomeric forms all torsion angles (β, γ, ε and χ') were set to each of three values: +60, –60 and 180, while the furanose ring pseudorotation angle was set to each of 10 different values. The geometry optimization had been carried out with the quantum-mechanical density functional theory (DFT) method with the B3LYP exchange-correlation functional and standard 6-31G(d, p) basis set. Single point electron energies in the optimized geometries were calculated at the MP2/6-311++G (d, p) theory level while the MP2/cc-pVQZ theory level was used to calculate the electron energies of energetically most favorable conformations for α- and β-anomeric forms. Normal vibrations anharmonicity was described by the spectroscopic second-order perturbation theory (PT2). Results. The exhaustive conformational families for both anomeric forms of the 2-deoxy-D- ribofuranose molecule were obtained. For α-anomer it consists of 89 structures and contains 107 structures for beta-anomer. It is shown, that in gas-phase mixture of α- and β-anomers in thermodynamic equilibrium the former would dominate with the 82:18 % ratio. In each conformation the normal vibrations anharmonicity has been characterized. Conclusions. The presence of B- and AII-DNA-like conformations in conformational families of both anomeric forms is established. It is shown, that the anharmonic correction can lead to significant (more than kBT under standard conditions) changes in the vibrational part of the molecule’s Gibbs energy.
Keywords: 2-deoxy-D-ribofuranose, conformational analysis, ab initio calculations

References

[1] Saenger W. Principles of nucleic acid structure. New York: Springer, 1984; 556 p.
[2] Kurita N., Danilov V. I., Anisimov V. M. The structure of Watson-Crick DNA base pairs obtained by MP2 optimization Chem. Phys. Lett 2005 404, N 1–3 P. 164–170.
[3] Bulavin LA, Nikolaienko TYu, Hovorun DM. Phosphoric acid molecule structural softness: ab initio quantum-mechanical study. Dopovidi Nats Akad Nauk Ukrainy. 2009;(10):80-7.
[4] Perahia D., Pullman B., Saran A. Molecular orbital calculations on the conformation of nucleic acids and their constituents. IX. The geometry of the phosphate group: Key to the conformation of polynucleotides? Biochim. Biophys. Acta 1974 340, N 3 P. 299–313.
[5] Furberg S. The crystal structure of phosphoric acid Acta Chim. Scand 1955 9, N 10 P. 1557–1566.
[6] Zhurakivsky RO, Hovorun DM. Conformational properties of 1'-deoxyribose, the model sugar residue of nucleosides: the DFT quantum-mechanical investigation. Dopovidi Nats Akad Nauk Ukrainy. 2008;(3):167-76.
[7] Zhurakivs'kyi RO, Yurenko EP, Hovorun DM. Conformational properties of 1',2'-deoxyribose - the model sugar residue of 2'-deoxyribonucleosides: results of a nonempiric quantum mechanical study. Dopovidi Nats Akad Nauk Ukrainy. 2006; (8):207-13.
[8] Brown N. S., Bicknell R. Thymidine phosphorylase, 2-deoxyD-ribose and angiogenesis Biochem. J 1998 334, pt 1 P. 1–8.
[9] Cloran F., Carmichael I., Serianni A. S. 2-Deoxy--D-erythro-pentofuranose: hydroxymethyl group conformation and substituent effects on molecular structure, ring geometry, and NMR spin-spin coupling constants from quantum chemical calculations J. Am. Chem. Soc 2001 123, N 20 P. 4781–4791.
[10] Frisch M. J., Trucks G. W., Schlegel H. B., Scuseria G. E., Robb M. A., Cheeseman J. R., Gaussian 03, Revision C.02 Wallingford CT: Gaussian Inc., 2004.
[11] Foresman J. B., Frisch A. Exploring chemistry with electronic structure methods Pittsburgh: Gaussian Inc., 1996.
[12] Boryskina O. P., Tkachenko M. Yu., Shestopalova A. V. Variability of DNA structure and protein-nucleic acid reconginition Biopolym. Cell 2010 26, N 5 P. 360–372.
[13] Barone V. Anharmonic vibrational properties by a fully automated second-order perturbative approach J. Chem. Phys 2005 122, N 1 P. 014108, 1–10.
[14] Barone V. Vibrational zero-point energies and thermodynamic functions beyond the harmonic approximation J. Chem. Phys 2004 120, N 7 P. 3059–3065.
[15] Irikura K. K., Johnson R. D., Kacker R. N. Uncertainties in scaling factors for ab initio vibrational frequencies J. Phys. Chem. A 2005 109, N 37 P. 8430–8437.
[16] Merrick J. P., Moran D., Radom L. An evaluation of harmonic vibrational frequency scale factors J. Phys. Chem. A 2007 111, N 45 P. 11683–11700.