DNA Bending [Strahs & Schlick, 2000]

Adenine-rich sequences (A-tracts) occur strategically in regions of curvature [Sinden, 1994] and draw continued interest because the origin and nature of bending (as well as interpretation of crystal and solution data) remain unresolved. To explain dissimilar bend directions in A-tract crystals, versus unique bend directions in solution, we simulated two A-tracts whose structures had also been determined crystallographically: CGCGA6CG and CGCA6GCG.

Our simulations captured a subtle conversion (from initially dissimilar to similar bends equivalent to minor groove compression) [Strahs & Schlick, 2000]. The subtlety of the changes from the crystallographic structures together with development of bend directions as predicted in solution experimental studies suggested that we could reconcile data from various studies through common bend-stabilizing motifs in our simulations [Young et al., 1995; Young & Beveridge, 1998; Sprous et al., 1999; Sherer et al., 1999].

These motifs include large amounts of propeller twisting in AT base pairs, consistent sugar puckering differences between A and T sugars, a minor-groove water spine with high occupancy and long lifetimes stabilized by a narrow minor groove, and large rolling on the 5' side in relatively straight A-tract moieties. We explain [Strahs & Schlick, 2000] part of the differing interpretations by program-dependent [Lu & Olson, 1999, Lu et al., 1999] measurements of curvature. We avoid this sensitivity in our study by defining global bend angles -- global roll and tilt -- that depend on all base pair bends rather than only on external values, as generally done.



Figure 8 from reference [Strahs & Schlick, 2000]. It shows the minor groove water spine of the two A-tract dodecamers as viewed looking into the minor groove. Only water molecules in the primary hydration shell interacting within the A-tract and connecting the A-tract to the base pair on the 3' side are shown.

Hydrogen bond interactions (green mesh surface) between the water molecule (oxygen, orange; hydrogen, white) and the adenine N3 (blue) and thymine O2 (red) atoms are represented.

Next to each water position, the percentage of the 2 ns molecular dynamics simulation during which spine interactions were observed is indicated in red. The average (blue) and maximum (red) lifetimes of each spine position are plotted in the lower panels.

Our calculated average and maximum lifetimes of the water spine molecules are in excellent agreement with various NMR studies. Intriguingly, the geometrical requirements of the water spine stabilize the propeller twist and perhaps the sugar puckering differences between the A and T sugars.




Relevant Links






Click to go back to Research