Molecular dynamics simulations of A • T-rich oligomers:: Sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A . T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na+ intrusions in the first layer of hydration of the minor groove, with long residence times, up to similar to3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove. (C) 2003 Wiley Periodicals, Inc.