Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences.

Interplay between divalent cations (Mg ²⁺ and Ca ²⁺ ) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT ₃₀ and dA ₃₀ , with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg ²⁺ and Ca ²⁺ , we find that Mg ²⁺ binds to the minor grooves and phosphate groups, whereas Ca ²⁺ binds specifically to the minor groove. Both Mg ²⁺ and Ca ²⁺ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA ₃₀ conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT ₃₀ populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.