Molecular mechanics studies on poly(purine).poly(pyrimidine) sequences in DNA: polymorphism and local variability

Energy minimization has been carried out on three poly(purine).poly(pyrimidine) sequences--d(G)10.d(C)10, d(A)10.d(T)10, and d(AG)5.d(CT)5--using the molecular mechanics program AMBER (Assisted Model Building and Energy Refinement). In order to extensively scan the conformational space available, five different helical models were studied, three of them being right-handed helices while the other two were left helical. For all three sequences the right-handed A- and B-type helices are energetically slightly preferred over the left helices, but the energy difference between the various right-handed helices is only marginal. A detailed analysis has been carried out to characterize the local structural variability in the refined structures, both in terms of torsion angles as well as other parameters such as base-pair tilt, wedge roll, and wedge tilt, etc. All three sequences exhibit similar structural features for a particular form, but both the forms A and B show significant deviations from fiber models. In particular, the A-form structures have higher unit rise (2.7 A), and lower unit twist (31 degrees) and base-pair tilt (12 degrees), compared to the fiber model, which has corresponding values of 2.56 A, 32.7 degrees, and 20 degrees, respectively. All these changes indicate that the refined models are closer to the A-form structure observed in crystals of oligonucleotides. In the refined B-for models, the helical parameters are close to the fiber B-form, although the torsion angles show considerable variations. None of the three sequences examined, including the d(A)n.d(T)n sequence, show any pronounced curvature for the B-form structure.