Molecular mechanics and dynamics calculations on (dA)10.cntdot.(dT)10 incorporating distance constraints derived from NMR relaxation measurements

Abstract
Structural constraints derived from proton NMR relaxation measurements on poly(dA).cntdot.poly(dT) in the form of interproton separations and orientation have been combined with molecular mechanics and annealed molecular dynamics calculations to derive a model for the solution-state structure of this molecule. Three different possible starting configurations, including the standard A and B forms of Arnott and Hukins [Arnott, S., and Hukins, D. W. L. (1972) Biochem. Biophys. Res. Commun. 47, 1506-1509] and the heteronomous (H) structure [Arnott, S., Chandrasekaran, R., Hall, I. H., and Puigjaner, L. C. (1983) Nucleic Acids Res. 11, 4141-4155], were examined. Both the B- and H-DNA structures converged to the same B-like structure (approximately C2''-endo conformation on both the A and T sugars, glycosidic bond torsional angle of 63-73.degree.) with the same energies and average helical parameters that gave good fits of the NMR relaxation rates. This model also accounts for the experimental observation [Behling, R. W., and Kearns, D. R. (1986) Biochemistry 25, 3335-3346] that the AH2 proton interacts more strongly with the H1'' sugar proton on the T strand than on the A strand. Although the helix repeat angle (39.degree.) is larger than that for standard B-DNA (36.degree.), this does not result in a significantly smaller minor groove, as monitored by the interstrand P-P separation. Calculations starting with the A-DNA structure lead to a very high energy structure that gave a poorer fit of the NMR data.