DFT study of polymorphism of the DNA double helix at the level of dinucleoside monophosphates

We apply DFT calculations to deoxydinucleoside monophosphates (dDMPs) which represent minimal fragments of the DNA chain to study the molecular basis of stability of the DNA duplex, the origin of its polymorphism and conformational heterogeneity. In this work, we continue our previous studies of dDMPs where we detected internal energy minima corresponding to the “classical” B conformation (BI-form), which is the dominant form in the crystals of oligonucleotide duplexes. We obtained BI local energy minima for all existing base sequences of dDMPs. In the present study, we extend our analysis to other families of DNA conformations, successfully identifying A, BI, and BII energy minima for all dDMP sequences. These conformations demonstrate distinct differences in sugar ring puckering, but similar sequence-dependent base arrangements. Internal energies of BI and BII conformers are close to each other for nearly all the base sequences. The dGpdG, dTpdG, and dCpdA dDMPs slightly favor the BII conformation, which agrees with these sequences being more frequently experimentally encountered in the BII form. We have found BII-like structures of dDMPs for the base sequences both existing in crystals in BII conformation and those not yet encountered in crystals till now. On the other hand, we failed to obtain dDMP energy minima corresponding to the Z family of DNA conformations, thus giving us the ground to conclude that these conformations are stabilized in both crystals and solutions by external factors, presumably by interactions with various components of the media. Overall the accumulated computational data demonstrate that the A, BI, and BII families of DNA conformations originate from the corresponding local energy minimum conformations of dDMPs, thus determining structural stability of a single DNA strand during the processes of unwinding and rewinding of DNA. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2548–2559, 2010

[1]  Heinz Sklenar,et al.  Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps. , 2004, Biophysical journal.

[2]  Jirí Cerný,et al.  Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations , 2007, J. Comput. Chem..

[3]  F. J. Luque,et al.  Dynamics of B-DNA on the microsecond time scale. , 2007, Journal of the American Chemical Society.

[4]  David A. Dixon,et al.  A local density functional study of the structure and vibrational frequencies of molecular transition-metal compounds , 1992 .

[5]  F. Matthias Bickelhaupt,et al.  Chemistry with ADF , 2001, J. Comput. Chem..

[6]  Igor V. Filippov,et al.  PROSIT: Pseudo-Rotational Online Service and Interactive Tool, Applied to a Conformational Survey of Nucleosides and Nucleotides , 2004, J. Chem. Inf. Model..

[7]  Daniel Svozil,et al.  Geometrical and electronic structure variability of the sugar-phosphate backbone in nucleic acids. , 2008, The journal of physical chemistry. B.

[8]  H. Berman,et al.  Geometric Parameters in Nucleic Acids: Sugar and Phosphate Constituents , 1996 .

[9]  S. Grimme,et al.  Density functional theory including dispersion corrections for intermolecular interactions in a large benchmark set of biologically relevant molecules. , 2006, Physical chemistry chemical physics : PCCP.

[10]  Stefan Grimme,et al.  Accurate description of van der Waals complexes by density functional theory including empirical corrections , 2004, J. Comput. Chem..

[11]  K. Nakao,et al.  Efficiency of the MO method using a Slater-type basis set and non-local density functional formalism for describing DNA base stacking energy , 1999 .

[12]  Dennis R. Salahub,et al.  Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation , 1992 .

[13]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[14]  A. R. Srinivasan,et al.  The nucleic acid database. A comprehensive relational database of three-dimensional structures of nucleic acids. , 1992, Biophysical journal.

[15]  Victor M. Anisimov,et al.  DFT study of minimal fragments of nucleic acid single chain for explication of sequence dependence of DNA duplex conformation , 2009 .

[16]  Andrzej Leś,et al.  DFT Study of B-like Conformations of Deoxydinucleoside Monophosphates Containing Gua and/or Cyt and their Complexes with Na+ Cation , 2008, Journal of Biomolecular Structure and Dynamics.

[17]  Daniel Svozil,et al.  DNA conformations and their sequence preferences , 2008, Nucleic acids research.

[18]  M. Swart,et al.  Performance of various density functionals for the hydrogen bonds in DNA base pairs , 2006 .

[19]  Pavel Hobza,et al.  Accurate interaction energies of hydrogen-bonded nucleic acid base pairs. , 2004, Journal of the American Chemical Society.

[20]  Noriyuki Kurita,et al.  Adjustment of Perdew-Wang exchange functional for describing van der Waals and DNA base-stacking interactions , 2003 .

[21]  Hans Peter Lüthi,et al.  Interaction energies of van der Waals and hydrogen bonded systems calculated using density functional theory: Assessing the PW91 model , 2001 .

[22]  M. Swart,et al.  π-π stacking tackled with density functional theory , 2007, Journal of molecular modeling.