No Arabic abstract
The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling $G$ between twisting and bending, known to arise from the groove asymmetry of the DNA double helix. Here, the effect of $G$ on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness $C_text{eff}$ for small twist-bend coupling. We find that the bare $G$ is screened by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with $G=0$, but with long-wavelength bending and twisting rigidities that are shifted by $G$-dependent amounts. Using results for torsional and bending rigidities for freely-fluctuating DNA, we show how our perturbative results can be extended to a non-perturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo triad and molecular dynamics oxDNA models, characterized by different degrees of coarse-graining, validating the perturbative and non-perturbative analyses. While our theory is in generally-good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely-used coarse-grained models.
By combining analytical and numerical calculations, we investigate the minimal-energy shape of short DNA loops of approximately $100$ base pairs (bp). We show that in these loops the excess twist density oscillates as a response to an imposed bending stress, as recently found in DNA minicircles and observed in nucleosomal DNA. These twist oscillations, here referred to as twist waves, are due to the coupling between twist and bending deformations, which in turn originates from the asymmetry between DNA major and minor grooves. We introduce a simple analytical variational shape, that reproduces the exact loop energy up to the fourth significant digit, and is in very good agreement with shapes obtained from coarse-grained simulations. We, finally, analyze the loop dynamics at room temperature, and show that the twist waves are robust against thermal fluctuations. They perform a normal diffusive motion, whose origin is briefly discussed.
Recent work indicates that twist-bend coupling plays an important role in DNA micromechanics. Here we investigate its effect on bent DNA. We provide an analytical solution of the minimum-energy shape of circular DNA, showing that twist-bend coupling induces sinusoidal twist waves. This solution is in excellent agreement with both coarse-grained simulations of minicircles and nucleosomal DNA data, which is bent and wrapped around histone proteins in a superhelical conformation. Our analysis shows that the observed twist oscillation in nucleosomal DNA, so far attributed to the interaction with the histone proteins, is an intrinsic feature of free bent DNA, and should be observable in other protein-DNA complexes.
We propose a statistical mechanics model for DNA melting in which base stacking and pairing are explicitly introduced as distinct degrees of freedom. Unlike previous approaches, this model describes thermal denaturation of DNA secondary structure in the whole experimentally accessible temperature range. Base pairing is described through a zipper model, base stacking through an Ising model. We present experimental data on the unstacking transition, obtained exploiting the observation that at moderately low pH this transition is moved down to experimentally accessible temperatures. These measurements confirm that the Ising model approach is indeed a good description of base stacking. On the other hand, comparison with the experiments points to the limitations of the simple zipper model description of base pairing.
Recent magnetic tweezers experiments have reported systematic deviations of the twist response of double-stranded DNA from the predictions of the twistable worm-like chain model. Here we show, by means of analytical results and computer simulations, that these discrepancies can be resolved if a coupling between twist and bend is introduced. We obtain an estimate of 40 $pm$ 10 nm for the twist-bend coupling constant. Our simulations are in good agreement with high-resolution, magnetic-tweezers torque data. Although the existence of twist-bend coupling was predicted long ago (Marko and Siggia, Macromolecules 27, 981 (1994)), its effects on the mechanical properties of DNA have been so far largely unexplored. We expect that this coupling plays an important role in several aspects of DNA statics and dynamics.
We consider self-propelled particles undergoing run-and-tumble dynamics (as exhibited by E. coli) in one dimension. Building on previous analyses at drift-diffusion level for the one-particle density, we add both interactions and noise, enabling discussion of domain formation by self-trapping, and other collective phenomena. Mapping onto detailed-balance systems is possible in certain cases.