Contribution of the intrinsic mechanical energy of the phosphodiester linkage to the relative stability of the A, BI, and BII forms of duplex DNA.

Abstract

Canonical forms of duplex DNA are known to sample well-defined regions of the alpha, beta, gamma, epsilon, and zeta dihedral angles that define the conformation of the phosphodiester linkage in the backbone of oligonucleotides. While extensive studies of base composition and base sequence dependent effects on the sampling of the A, B1, and BII canonical forms of duplex DNA have been presented, our understanding of the intrinsic contribution of the five dihedral degrees of freedom associated with the phosphodiester linkage to the conformational properties of duplex DNA is still limited. To better understand this contribution, ab initio quantum mechanical (QM) calculations were performed on a model compound representative of the phosphodiester backbone to systematically sample the energetics about the alpha, beta, gamma, epsilon, and zeta dihedral angles relevant to the conformational properties of duplex DNA. Low-energy regions of dihedral potential energy surfaces are shown to correlate with the regions of dihedral space sampled in experimental crystal structures of the canonical forms of DNA, validating the utility of the model compound and emphasizing the contribution of the intrinsic mechanical properties of the phosphodiester backbone to the conformational properties of duplex DNA. Those contributions include the relative stability of the A, BI, and BII conformations of duplex DNA, where the gas-phase energetics favor the BI form over the A and BII forms. In addition, subtle features of the potential energy surfaces mimic changes in the probability distributions of alpha, beta, gamma, epsilon, and zeta dihedral angles in A, BI, and BII forms of DNA as well as with conformations sampled in single-stranded DNA. These results show that the intrinsic mechanical properties of the phosphodiester backbone make a significant contribution to conformational properties of duplex DNA observed in the condensed phase and allow for the prediction that single-stranded DNA primarily samples folded conformations thereby possibly lowering the entropic barrier to the formation of duplex DNA.