Grand-canonical simulation of DNA condensation with two salts, affect of divalent counterion size

TL;DR

This study uses grand-canonical Monte Carlo simulations to explore how the size of divalent counterions influences DNA condensation behavior in restricted environments, revealing size-dependent effects on DNA aggregation.

Contribution

It introduces a simple parameter, counterion size, to explain variations in DNA condensation and details a new simulation method for systems with two salts.

Findings

Counterion size affects DNA condensation behavior.

Smaller or larger counterions induce different DNA aggregation states.

The simulation method handles two-salt systems effectively.

Abstract

The problem of DNA$-$DNA interaction mediated by divalent counterions is studied using a generalized Grand-canonical Monte-Carlo simulation for a system of two salts. The effect of the divalent counterion size on the condensation behavior of the DNA bundle is investigated. Experimentally, it is known that multivalent counterions have strong effect on the DNA condensation phenomenon. While tri- and tetra-valent counterions are shown to easily condense free DNA molecules in solution into toroidal bundles, the situation with divalent counterions are not as clear cut. Some divalent counterions like Mg$^{+2}$ are not able to condense free DNA molecules in solution, while some like Mn$^{+2}$ can condense them into disorder bundles. In restricted environment such as in two dimensional system or inside viral capsid, Mg$^{+2}$ can have strong effect and able to condense them, but the condensation varies qualitatively with different system, different coions. It has been suggested that divalent counterions can induce attraction between DNA molecules but the strength of the attraction is not strong enough to condense free DNA in solution. However, if the configuration entropy of DNA is restricted, these attractions are enough to cause appreciable effects. The variations among different divalent salts might be due to the hydration effect of the divalent counterions. In this paper, we try to understand this variation using a very simple parameter, the size of the divalent counterions. We investigate how divalent counterions with different sizes can leads to varying qualitative behavior of DNA condensation in restricted environments. Additionally a Grand canonical Monte-Carlo method for simulation of systems with two different salts is presented in detail.