
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.

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.
