No Arabic abstract
We propose a new statistical mechanics model for the melting transition of DNA. Base pairing and stacking are treated as separate degrees of freedom, and the interplay between pairing and stacking is described by a set of local rules which mimic the geometrical constraints in the real molecule. This microscopic mechanism intrinsically accounts for the cooperativity related to the free energy penalty of bubble nucleation. The model describes both the unpairing and unstacking parts of the spectroscopically determined experimental melting curves. Furthermore, the model explains the observed temperature dependence of the effective thermodynamic parameters used in models of the nearest neighbor (NN) type. We compute the partition function for the model through the transfer matrix formalism, which we also generalize to include non local chain entropy terms. This part introduces a new parametrization of the Yeramian-like transfer matrix approach to the Poland-Scheraga description of DNA melting. The model is exactly solvable in the homogeneous thermodynamic limit, and we calculate all observables without use of the grand partition function. As is well known, models of this class have a first order or continuous phase transition at the temperature of complete strand separation depending on the value of the exponent of the bubble entropy.
We study the melting of a double stranded DNA in the presence of stretching forces, via 3D Monte-Carlo simulations, exactly solvable models and heuristic arguments. The resulting force-temperature phase diagram is dramatically different for the cases where the force is applied to only one strand or to both. Different assumptions on the monomer size of single and double stranded DNA lead to opposite conclusions as to whether DNA melts or not as it overstretches.
The flexibility and stiffness of small DNA play a fundamental role ranging from several biophysical processes to nano-technological applications. Here, we estimate the mechanical properties of short double-stranded DNA (dsDNA) having length ranging from 12 base-pairs (bps) to 56 bps, paranemic crossover (PX) DNA, and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grain models $-$ Martini and oxDNA. To calculate the persistence length ($L_p$) and the stretch modulus ($gamma$) of the dsDNA, we incorporate the worm-like chain and elastic rod model, while for DNT, we implement our previously developed theoretical framework. We compare and contrast all the results with previously reported all-atom molecular dynamics (MD) simulation and experimental results. The mechanical properties of dsDNA ($L_p$ $sim$ 50nm, $gamma sim$ 800-1500 pN), PX DNA ($gamma sim$ 1600-2000 pN) and DNTs ($L_p sim 1-10 mu$m, $gamma sim$ 6000-8000 pN) estimated using Martini soft elastic network and oxDNA are in very good agreement with the all-atom MD and experimental values, while the stiff elastic network Martini reproduces order of magnitude higher values of $L_p$ and $gamma$. The high flexibility of small dsDNA is also depicted in our calculations. However, Martini models proved inadequate to capture the salt concentration effects on the mechanical properties with increasing salt molarity. OxDNA captures the salt concentration effect on small dsDNA mechanics. But it is found to be ineffective to reproduce the salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time evolved PX DNA and DNT structures from the oxDNA models are comparable to the all-atom MD simulated structures. Our findings provide a route to study the mechanical properties of DNA nanostructures with increased time and length scales and has a remarkable implication in the context of DNA nanotechnology.
We discuss the distribution of ions around highly charged PEs when there is competition between monovalent and multivalent ions, pointing out that in this case the number of condensed ions is sensitive to short-range interactions, salt, and model-dependent approximations. This sensitivity is discussed in the context of recent experiments on DNA aggregation, induced by multivalent counterions such as spermine and spermidine.
The amount and type of self-entanglement of DNA filaments is significantly affected by spatial confinement, which is ubiquitous in biological systems. Motivated by recent advancements in single DNA molecule experiments based on nanofluidic devices, and by the introduction of algorithms capable of detecting knots in open chains, we investigate numerically the entanglement of linear, open DNA chains confined inside nano-slits. The results regard the abundance, type and length of occurring knots and are compared with recent findings for DNA inside nano-channels. In both cases, the width of the confining region, D, spans the 30nm- 1mu m range and the confined DNA chains are 1 to 4mu m long. It is found that the knotting probability is maximum for slit widths in the 70-100nm range. However, over the considered DNA contour lengths, the maximum incidence of knots remains below 20%, while for channel confinement it tops 50%. Further differences of the entanglement are seen for the average contour length of the knotted region which drops significantly below D ~100nm for channel-confinement, while it stays approximately constant for slit-like confinement. These properties ought to reverberate in different kinetic properties of linear DNA depending on confinement and could be detectable experimentally or exploitable in nano-technological applications.
We study DNA self-assembly and DNA computation using a coarse-grained DNA model within the directional dynamic bonding framework {[}C. Svaneborg, Comp. Phys. Comm. 183, 1793 (2012){]}. In our model, a single nucleotide or domain is represented by a single interaction site. Complementary sites can reversibly hybridize and dehybridize during a simulation. This bond dynamics induces a dynamics of the angular and dihedral bonds, that model the collective effects of chemical structure on the hybridization dynamics. We use the DNA model to perform simulations of the self-assembly kinetics of DNA tetrahedra, an icosahedron, as well as strand displacement operations used in DNA computation.