The effective DNA-DNA interaction force is calculated by computer simulations with explicit tetravalent counterions and monovalent salt. For overcharged DNA molecules, the interaction force shows a double-minimum structure. The positions and depths of these minima are regulated by the counterion density in the bulk. Using two-dimensional lattice sum and free energy perturbation theories, the coexisting phases for DNA bundles are calculated. A DNA-condensation and redissolution transition and a stable mesocrystal with an intermediate lattice constant for high counterion concentration are obtained.
The distance-resolved effective interaction potential between two parallel DNA molecules is calculated by computer simulations with explicit tetravalent counterions and monovalent salt. Adding counterions first yields an attractive minimum in the potential at short distances which then disappears in favor of a shallower minimum at larger separations. The resulting phase diagram includes a DNA-condensation and redissolution transition and a stable mesocrystal with an intermediate lattice constant for high counterion concentration.
The problem of DNA-DNA interaction mediated by divalent counterions is studied using computer simulation. Although divalent counterions cannot condense free DNA molecules in solution, we show that if DNA configurational entropy is restricted, divalent counterions can cause DNA reentrant condensation similar to that caused by tri- or tetra-valent counterions. DNA-DNA interaction is strongly repulsive at small or large counterion concentration and is negligible or slightly attractive for a concentration in between. Implications of our results to experiments of DNA ejection from bacteriophages are discussed. The quantitative result serves to understand electrostatic effects in other experiments involving DNA and divalent counterions.
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.
We develop a robust coarse-grained model for single and double stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduces the experimentally measured salt and sequence dependence of the size of single stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin, and could be described by a standard polymer model. We also show that the persistence length of double stranded DNA is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.
Topology affects physical and biological properties of DNA and impacts fundamental cellular processes, such as gene expression, genome replication, chromosome structure and segregation. In all organisms DNA topology is carefully modulated and the supercoiling degree of defined genome regions may change according to physiological and environmental conditions. Elucidation of structural properties of DNA molecules with different topology may thus help to better understand genome functions. Whereas a number of structural studies have been published on highly negatively supercoiled DNA molecules, only preliminary observations of highly positively supercoiled are available, and a description of DNA structural properties over the full range of supercoiling degree is lacking. Atomic Force Microscopy (AFM) is a powerful tool to study DNA structure at single molecule level. We here report a comprehensive analysis by AFM of DNA plasmid molecules with defined supercoiling degree, covering the full spectrum of biologically relevant topologies, under different observation conditions. Our data, supported by statistical and biochemical analyses, revealed striking differences in the behavior of positive and negative plasmid molecules.