Molecular‐mechanics studies on d(CGCGAATTCGCG)2 and dA12·dT12: An illustration of the coupling between sugar repuckering and DNA twisting

Abstract

AbstractMolecular‐mechanics calculations have been carried out on the base‐paired deoxy dodecanucleoside undecaphosphates d(CGCGAATTCGCG)2 and d(A12)·d(T12). These refinements were carried out using the model‐built Arnott B‐DNA geometry as initial coordinates (with a helix repaeat of 10.0 residues/turn), as well as helix repeats ranging from 9 to 12 residues/turn. There was some variation in the optimum calculated helix repeat, depending on the dielectric model, the presence or absence of counterions, and the method used for inclusion for nonbonded interactions; the most interesting general result of these calculations was the coupling between furanose sugar puckering and twist. This coupling was observed for all models. With a helix repeat of 9.0 residues/turn, all sugars remain C(2′)endo after refinement; as the helix repear increases to 12.0 residues/turn, the number of sugars repuckering to O(1′)endo and C(3′)endo increases also. With our most rigorous model (i.e., a model with no cutoff distance for nonbonded interactions) and a helix repeat of 10.0 residues/turn, we find a greater tendency for pyrimidine than purine repuckering in d(CGCGAATTCGCG)2, in agreement with the x‐ray structural data of Drew et al. [(1981) Proc. Natl. Acad. Sci. USA 78, 2179–2185].We also carried out a number of calculations in which we “forced” one of two deoxy sugars to repucker or one of the C3′‐O3′‐P‐O5′ (ω) torsion angles to change from gauche− to trans using dihedral angle constraints. After the constraints were removed, some of these structures “reverted” to the sugar pucker of the initial structures, while others remained repuckered. In all cases, the energies for repuckered structures after refinement were very similar to energies of the initial structure. Experiments and theory suggest that local conformational fluctuations play an essential role in nmr relaxation of 31P and 13C atoms in double‐helical DNA. The results of our previous calculations on hexanucleoside phosphates and the calculations presented there are consistent with an important contribution to nmr relaxation processes of conformational changes in the torsion angle ω′ from gauche− to trans and deoxy sugar repuckering from C(2′)endo to C(3′)endo. Specifically, the calculations presented here indicate a very flexible phosphate backbone in helixes having an intermediate helix repeat of 10 to 11 residues/turn. These helixes may accommodate sugars of variable pucker without significantly disrupting base–base hydrogen‐bonding and stacking interactions. All of the variant structures are similar in energy, suggesting that conversion between them can occur on a nanosecond time scale, as observed in nmr relaxation experiments.