Ions around DNA: Monte Carlo estimates of distribution with improved electrostatic potentials.

Abstract

AbstractCoulomb's law for the electrostatic interactions between ions is modified when discontinuities in dielectric constant (relative permittivity) occur. In a DNA solution such a discontinuity occurs at the interface between the DNA molecular helix and the surrounding water. We take the modified interaction potentials from a previous report [Macromolecules (1986) 19, 1186–1194] and use them with the Monte Carlo method to find the distribution of univalent and bivalent counterions around the DNA helix in the absence of coions (i.e., no added salt). In comparing the ion distribution with the modified potential to the distribution without, we find that the effects of the modifications to the potentials are considerable. The modifications tend to drive the ions out of the grooves of the helix, especially out of the major groove. This result comes partly from the repulsion exerted on the ions by the low‐permittivity helix and partly from the concentration of the field of the phosphates at the surface of the helix, a concentration that is also caused by the discontinuity in permittivity.