Search

A stretched DNA molecule which is also under- or overwound, undergoes a buckling transition forming intertwined looped domains called plectonemes. Here we develop a simple theory that extends the two-phase model of stretched supercoiled DNA, allowing for the coexistence of multiple plectonemic domains by including positional and length distribution entropies. Such a multiplectoneme phase is favored in long DNA molecules in which the gain of positional entropy compensates for the cost of nucleating a plectoneme along a stretched DNA segment. Despite its simplicity, the developed theory is shown to be in excellent agreement with Monte Carlo simulations of the twistable wormlike chain model. The theory predicts more plectonemes than experimentally observed, which we attribute to the limited resolution of experimental data. Since plectonemes are detected through fluorescence signals, those shorter than the observable threshold are likely missed., Comment: 15 pages, 11 figures

Coarse-grained models have played an important role in the study of the behavior of DNA at length scales beyond a few hundred base pairs. Traditionally, these models have relied on structurally featureless and sequence-independent approaches, such as the twistable worm-like chain. However, research over the past decade has illuminated the substantial impact of DNA sequence even at the kilo-base pair scale. Several robust sequence-dependent models have emerged, capturing intricacies at the base pair-step level. Here we introduce an analytical framework for coarse-graining such models to the 2 to 40-base pair scale while preserving essential structural and dynamical features. These faithful coarse-grained parametrizations enable efficient sampling of large molecules. Rather than providing a fully parametrized model, we present the methodology and software necessary for mapping any base pair-step model to the desired level of coarse-graining. Finally, we provide application examples of our method, including estimates of the persistence length and effective torsional stiffness of DNA in a setup mimicking a freely orbiting tweezer, as well as simulations of intrinsically helical DNA., Comment: 23 pages, 9 figures

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e. over- or under-wounded. They are ubiquitous in cellular DNA and their physical properties have attracted significant interest both from the experimental and modeling side. In this work, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched and of a plectonemic phase. Several different mechanisms are found to contribute to extension fluctuations, characterized by the extension variance. We find the dominant contribution to the variance to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the Twistable Wormlike Chain and analyze the fluctuation of various quantities, which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations., Comment: 13 pages, 8 figures

While analyzing the persistence length of DNA atomistic simulations Wu et al. [J. Chem. Phys. 142, 125103 (2015)] introduced an empirical formula to account for the observed length-dependence. In particular they found that the persistence length increases with the distance. Here, we derive the formula by Wu et al. using a non-local twistable wormlike chain which introduces couplings between distal sites. Finally, we show that the same formula can account for the length-scale dependence of the torsional persistence length and is, in fact, applicable to any kind of polymer model with non-local couplings., Comment: 3 pages, 1 Figure

The mechanical properties of DNA are typically described by elastic theories with purely local couplings (on-site models). We discuss and analyze coarse-grained (oxDNA) and all-atom simulations, which indicate that in DNA distal sites are coupled. Hence, off-site models provide a more realistic description of the mechanics of the double helix. We show that off-site interactions are responsible for a length scale dependence of the elasticity, and we develop an analytical framework to estimate bending and torsional persistence lengths in models including these interactions. Our simulations indicate that off-site couplings are particularly strong for certain degrees of freedom, while they are very weak for others. If stiffness parameters obtained from DNA data are used, the theory predicts large length scale dependent effects for torsional fluctuations and a modest effect in bending fluctuations, which is in agreement with experiments., Comment: 16 pages, 9 figures; video abstract at https://youtu.be/8FUdFw8IbuY

To study the elastic properties of rod-like DNA nanostructures, we perform long simulations of these structure using the oxDNA coarse-grained model. By analysing the fluctuations in these trajectories we obtain estimates of the bend and twist persistence lengths, and the underlying bend and twist elastic moduli and couplings between them. Only on length scales beyond those associated with the spacings between the interhelix crossovers do the bending fluctuations behave like those of a worm-like chain. The obtained bending persistence lengths are much larger than that for double-stranded DNA and increase non-linearly with the number of helices, whereas the twist moduli increase approximately linearly. To within the numerical error in our data, the twist-bend coupling constants are of order zero. That the bending persistence lengths we obtain are generally somewhat higher than in experiment probably reflects both that the simulated origami have no assembly defects and that the oxDNA extensional modulus for double-stranded DNA is too large.

By combining analytical results and simulations of various coarse-grained models we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternance of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, not the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants., Comment: 11 pages, 7 figures

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., Comment: 16 pages, 6 figures

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.

It is well-established that many physical properties of DNA at sufficiently long length scales can be understood by means of simple polymer models. One of the most widely used elasticity models for DNA is the twistable worm-like chain (TWLC), which describes the double helix as a continuous elastic rod with bending and torsional stiffness. An extension of the TWLC, which has recently received some attention, is the model by Marko and Siggia, who introduced an additional twist-bend coupling, expected to arise from the groove asymmetry. By performing computer simulations of two available versions of oxDNA, a coarse-grained model of nucleic acids, we investigate the microscopic origin of twist-bend coupling. We show that this interaction is negligible in the oxDNA version with symmetric grooves, while it appears in the oxDNA version with asymmetric grooves. Our analysis is based on the calculation of the covariance matrix of equilibrium deformations, from which the stiffness parameters are obtained. The estimated twist-bend coupling coefficient from oxDNA simulations is $G=30\pm1$~nm. The groove asymmetry induces a novel twist length scale and an associated renormalized twist stiffness $\kappa_{\rm t} \approx 80$~nm, which is different from the intrinsic torsional stiffness $C \approx 110$~nm. This naturally explains the large variations on experimental estimates of the intrinsic stiffness performed in the past., Comment: 20 pages, 8 figures