DNA polymorphism and local variation in base-pair orientation: a theoretical rationale.

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B----A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.