Divalent cations and the electrostatic potential around DNA: Monte Carlo and Poisson-Boltzmann calculations

Area of Science:

• Biophysics
• Computational chemistry
• Molecular modeling

Background:

• Counterion condensation theory explains ion distribution around charged polymers like DNA.
• Understanding ion interactions is crucial for DNA structure, function, and stability.
• Divalent ions play a significant role in DNA condensation and recognition.

Purpose of the Study:

• To test counterion condensation theory predictions for divalent ions against Monte Carlo (MC) simulations.
• To investigate monovalent-divalent ion competition at the DNA surface.
• To evaluate the accuracy of Poisson-Boltzmann (PB) calculations for predicting divalent ion concentrations near DNA.

Main Methods:

• All-atom Monte Carlo (MC) simulations of DNA.
• Varying the partial molar volume of divalent ions to study competition.
• Comparison of MC results with Poisson-Boltzmann (PB) calculations using various dielectric continuum models.

Main Results:

• MC simulations validated counterion condensation theory predictions for divalent ions.
• Standard PB calculations underestimated divalent ion surface densities by 25-30% compared to MC.
• Errors from PB's mean-field approximation were smaller than those from dielectric continuum models.

Conclusions:

• Counterion condensation theory provides a reasonable framework for divalent ion behavior around DNA.
• Poisson-Boltzmann calculations require careful consideration of dielectric models for accurate divalent ion predictions.
• MC simulations offer a more accurate approach for quantifying ion concentrations at the DNA surface.