Influence of the dynamic positions of cations on the structure of the DNA minor groove: sequence-dependent effects

Different models for minor groove structures predict that the conformation is essentially fixed by sequence and has an influence on local ion distribution or alternatively that temporal positions of ions around the minor groove can affect the structure if they neutralize cross-strand phosphate charges. Our previous studies show that the minor groove in an AATT dodecamer responds to local sodium ion positions and is narrow when ions neutralize cross-strand phosphate-phosphate charges [J. Am. Chem. Soc. 2000, 122, 10513-10520]. Previous results from a number of laboratories have shown that G-tracts often have a wider minor groove than A-tracts, but they do not indicate whether this is due to reduced flexibility or differences in ion interactions. We have undertaken a molecular dynamics study of a d(TATAGGCCTATA) duplex to answer this question. The results show that the G-tract has the same amplitude of minor groove fluctuations as the A-tract sequence but that it has fewer ion interactions that neutralize cross-strand phosphate charges. These results demonstrate that differences in time-average groove width between A- and G-tracts are due to differences in ion interactions at the minor groove. When ions neutralize the cross-strand phosphates, the minor groove is narrow. When there are no neutralizing ion interactions, the minor groove is wide. The population of structures with no ion interactions is larger with the GGCC than with the AATT duplex, and GGCC has a wider time-average minor groove in agreement with experiment.