No Arabic abstract
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.
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.
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 means of computer simulations of a coarse-grained DNA model we show that the DNA hairpin zippering dynamics is anomalous, i.e. the characteristic time T scales non-linearly with N, the hairpin length: T ~ N^a with a>1. This is in sharp contrast with the prediction of the zipper model for which T ~ N. We show that the anomalous dynamics originates from an increase in the friction during zippering due to the tension built in the closing strands. From a simple polymer model we get a = 1+ nu = 1.59 with nu the Flory exponent, a result which is in agreement with the simulations. We discuss transition path times data where such effects should be detected.
Test experiments of hybridization in DNA microarrays show systematic deviations from the equilibrium isotherms. We argue that these deviations are due to the presence of a partially hybridized long-lived state, which we include in a kinetic model. Experiments confirm the model predictions for the intensity vs. free energy behavior. The existence of slow relaxation phenomena has important consequences for the specificity of microarrays as devices for the detection of a target sequence from a complex mixture of nucleic acids.