Twist‐dependent stacking energy of base‐pair steps in B‐DNA geometry: A density functional theory approach

Stacking energy of all the 10 unique DNA base‐pair steps (bp step) are calculated using density functional theory within the ultrasoft pseudopotential plane wave method and local density approximation for the exchange‐correlation functional. We have studied the dependence of stacking energy on twist angle, an aspect found difficult to explain using classical theory. We have found that the twist angle for different bp steps at stacking energy minimum matches extremely well with the values of average twist obtained from B‐DNA crystal structure data. This indicates that the use of a proper quantum chemical method to calculate the π‐π electronic interactions may explain stacking energy without incorporating hydrophobic interaction through solvent or effect of backbone through pseudobond. From the twist angle‐dependent stacking energy profile, we have also generated the probability distributions of twist for all the bp steps and calculated the variance of the distribution. Our calculated variances show similar trend to that of the experimental data for which sufficient numbers of data are available. The TA, AT, and CG doublets show large variances among the 10 possible bp steps, indicating their maximum flexibility. This might be the case of unusual deformation observed at the TATA‐box while binding to TBP protein. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008