Search

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs. all-atom models using stiffness matrices in Fourier space ($q$-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version cgDNA+. Our analysis shows that while cgDNA+ follows closely the $q$-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables which are rather strong in the $q \to 0$ (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects which manifest themselves at longer length scales. In particular, cgDNA performs poorly on the twist stiffness with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist is strongly length scale dependent: DNA is torsionally soft at a few base pair distances, but becomes more rigid at distances of a few dozens base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated., Comment: 9 pages, 6 figures

All-atom simulations have become increasingly popular to study conformational and dynamical properties of nucleic acids as they are accurate and provide high spatial and time resolutions. This high resolution however comes at a heavy computational cost and within the time scales of simulations nucleic acids weakly fluctuate around their ideal structure exploring a limited set of conformations. We introduce the RBB-NA algorithm which is capable of controlling rigid base parameters in all-atom simulations of Nucleic Acids. With suitable biasing potentials this algorithm can "force" a DNA or RNA molecule to assume specific values of the six rotational (tilt, roll, twist, buckle, propeller, opening) and/or the six translational parameters (shift, slide, rise, shear, stretch, stagger). The algorithm enables the use of advanced sampling techniques to probe the structure and dynamics of locally strongly deformed Nucleic Acids. We illustrate its performance showing some examples in which DNA is strongly twisted, bent or locally buckled. In these examples RBB-NA reproduces well the unconstrained simulations data and other known features of DNA mechanics, but it also allows one to explore the anharmonic behavior characterizing the mechanics of nucleic acids in the high deformation regime., Comment: 12 pages, 6 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

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

We derive an analytical expression for the transition path time (TPT) distribution for a one-dimensional particle crossing a parabolic barrier. The solution is expressed in terms of the eigenfunctions and eigenvalues of the associated Fokker-Planck equation. The particle performs an anomalous dynamics generated by a power-law memory kernel, which includes memoryless Markovian dynamics as a limiting case. Our result takes into account absorbing boundary conditions, extending existing results obtained for free boundaries. We show that TPT distributions obtained from numerical simulations are in excellent agreement with analytical results, while the typically employed free boundary conditions lead to a systematic overestimation of the barrier height. These findings may be useful in the analysis of experimental results on transition path times. A web tool to perform this analysis is freely available., Comment: 7 pages; 6 figures

The capture and translocation of biomolecules through nanometer-scale pores are processes with a potential large number of applications, and hence they have been intensively studied in the recent years. The aim of this paper is to review existing models of the capture process by a nanopore, together with some recent experimental data of short single- and double-stranded DNA captured by Cytolysin A (ClyA) nanopore. ClyA is a transmembrane protein of bacterial origin which has been recently engineered through site-specific mutations, to allow the translocation of double- and single-stranded DNA. A comparison between theoretical estimations and experiments suggests that for both cases the capture is a reaction-limited process. This is corroborated by the observed salt dependence of the capture rate, which we find to be in quantitative agreement with the theoretical predictions., Comment: Published in JPCM Special Issue "Transport in Narrow Channels"

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

We develop a theoretical description of the critical zipping dynamics of a self-folding polymer. We use tension propagation theory and the formalism of the generalized Langevin equation applied to a polymer that contains two complementary parts which can bind to each other. At the critical temperature, the (un)zipping is unbiased and the two strands open and close as a zipper. The number of closed base pairs $n(t)$ displays a subdiffusive motion characterized by a variance growing as $\langle \Delta n^2(t) \rangle \sim t^\alpha$ with $\alpha < 1$ at long times. Our theory provides an estimate of both the asymptotic anomalous exponent $\alpha$ and of the subleading correction term, which are both in excellent agreement with numerical simulations. The results indicate that the tension propagation theory captures the relevant features of the dynamics and shed some new insights on related polymer problems characterized by anomalous dynamical behavior., Comment: 8 pages, 3 figures, submitted

By using the blob theory and computer simulations, we investigate the properties of a linear polymer performing a stationary rotational motion around a long impenetrable rod. In particular, in the simulations the rotation is induced by a torque applied to the end of the polymer that is tethered to the rod. Three different regimes are found, in close analogy with the case of polymers pulled by a constant force at one end. For low torques the polymer rotates maintaining its equilibrium conformation. At intermediate torques the polymer assumes a trumpet shape, being composed by blobs of increasing size. At even larger torques the polymer is partially wrapped around the rod. We derive several scaling relations between various quantities as angular velocity, elongation and torque. The analytical predictions match the simulation data well. Interestingly, we find a "thinning" regime where the torque has a very weak (logarithmic) dependence on the angular velocity. We discuss the origin of this behavior, which has no counterpart in polymers pulled by an applied force., Comment: 30 pages, 8 figures, 1 TOC figure; video abstract at https://youtu.be/LwicoSkh3mI

Some recent results on the rotational dynamics of polymers are reviewed and extended. We focus here on the relaxation of a polymer, either flexible or semiflexible, initially wrapped around a rigid rod. We also study the steady polymer rotation generated by a constant torque on the rod. The interplay of frictional and entropic forces leads to a complex dynamical behavior characterized by non-trivial universal exponents. The results are based on extensive simulations of polymers undergoing Rouse dynamics and on an analytical approach using force balance and scaling arguments. The analytical results are in general in good agreement with the simulations, showing how a simplified approach can correctly capture the complex dynamical behavior of rotating polymers., Comment: 13 pages; 7 figures; proceedings of the International Workshop on "Brownian Motion in Confined Geometries", Max Planck Institute for the Physics of Complex Systems in Dresden from 17 - 21 March 2014; to appear in EPJ-ST

We study the rotational dynamics of a flexible polymer initially wrapped around a rigid rod and unwinding from it. This dynamics is of interest in several problems in biology and constitutes a fundamental instance of polymer relaxation from a state of minimal entropy. We investigate the dynamics of several quantities such as the total and local winding angles and metric quantities. The results of simulations performed in two and three dimensions, with and without self-avoidance, are explained by a theory based on scaling arguments and on a balance between frictional and entropic forces. The early stage of the dynamics is particularly rich, being characterized by three coexisting phases., Comment: 9 pages, 9 figures

We briefly review some of the models used to describe DNA denaturation dynamics, focusing on the value of the dynamical exponent $z$, which governs the scaling of the characteristic time $\tau\sim L^z$ as a function of the sequence length $L$. The models contain different degrees of simplifications, in particular sometimes they do not include a description for helical entanglement: we discuss how this aspect influences the value of $z$, which ranges from $z=0$ to $z \approx 3.3$. Connections with experiments are also mentioned.

The relaxation dynamics of a polymer wound around a fixed obstacle constitutes a fundamental instance of polymer with twist and torque and it is of relevance also for DNA denaturation dynamics. We investigate it by simulations and Langevin equation analysis. The latter predicts a relaxation time scaling as a power of the polymer length times a logarithmic correction related to the equilibrium fluctuations of the winding angle. The numerical data support this result and show that at short times the winding angle decreases as a power-law. This is also in agreement with the Langevin equation provided a winding-dependent friction is used, suggesting that such reduced description of the system captures the basic features of the problem., Comment: 4 pages, 5 figures. Accepted for publication in Phys. Rev. Lett

In this article it is shown how optimized and dedicated microarray experiments can be used to study the thermodynamics of DNA hybridization for a large number of different conformations in a highly parallel fashion. In particular, free energy penalties for mismatches are obtained in two independent ways and are shown to be correlated with values from melting experiments in solution reported in the literature. The additivity principle, which is at the basis of the nearest-neighbor model, and according to which the penalty for two isolated mismatches is equal to the sum of the independent penalties, is thoroughly tested. Additivity is shown to break down for a mismatch distance below 5 nt. The behavior of mismatches in the vicinity of the helix edges, and the behavior of tandem mismatches are also investigated. Finally, some thermodynamic outlying sequences are observed and highlighted. These sequences contain combinations of GA mismatches. The analysis of the microarray data reported in this article provides new insights on the DNA hybridization parameters and can help to increase the accuracy of hybridization-based technologies., Comment: 13 pages, 11 figures, 1 table, Supplementary Data available in Appendix

We consider two complementary polymer strands of length $L$ attached by a common end monomer. The two strands bind through complementary monomers and at low temperatures form a double stranded conformation (zipping), while at high temperature they dissociate (unzipping). This is a simple model of DNA (or RNA) hairpin formation. Here we investigate the dynamics of the strands at the equilibrium critical temperature $T=T_c$ using Monte Carlo Rouse dynamics. We find that the dynamics is anomalous, with a characteristic time scaling as $\tau \sim L^{2.26(2)}$, exceeding the Rouse time $\sim L^{2.18}$. We investigate the probability distribution function, the velocity autocorrelation function, the survival probability and boundary behaviour of the underlying stochastic process. These quantities scale as expected from a fractional Brownian motion with a Hurst exponent $H=0.44(1)$. We discuss similarities and differences with unbiased polymer translocation., Comment: 7 pages, 8 figures

The winding angle probability distribution of a planar self-avoiding walk has been known exactly since a long time: it has a gaussian shape with a variance growing as $<\theta^2>\sim \ln L$. For the three-dimensional case of a walk winding around a bar, the same scaling is suggested, based on a first-order epsilon-expansion. We tested this three-dimensional case by means of Monte Carlo simulations up to length $L\approx25\,000$ and using exact enumeration data for sizes $L\le20$. We find that the variance of the winding angle scales as $<\theta^2>\sim (\ln L)^{2\alpha}$, with $\alpha=0.75(1)$. The ratio $\gamma = <\theta^4>/<\theta^2>^2=3.74(5)$ is incompatible with the gaussian value $\gamma =3$, but consistent with the observation that the tail of the probability distribution function $p(\theta)$ is found to decrease slower than a gaussian function. These findings are at odds with the existing first-order $\epsilon$-expansion results., Comment: 18 pages, 12 figures, 1 table

We investigate by Monte Carlo simulations the zipping and unzipping dynamics of two polymers connected by one end and subject to an attractive interaction between complementary monomers. In zipping, the polymers are quenched from a high temperature equilibrium configuration to a low temperature state, so that the two strands zip up by closing up a "Y"-fork. In unzipping, the polymers are brought from a low temperature double stranded configuration to high temperatures, so that the two strands separate. Simulations show that the unzipping time, $\tau_u$, scales as a function of the polymer length as $\tau_u \sim L$, while the zipping is characterized by anomalous dynamics $\tau_z \sim L^\alpha$ with $\alpha = 1.37(2)$. This exponent is in good agreement with simulation results and theoretical predictions for the scaling of the translocation time of a forced polymer passing through a narrow pore. We find that the exponent $\alpha$ is robust against variations of parameters and temperature, whereas the scaling of $\tau_z$ as a function of the driving force shows the existence of two different regimes: the weak forcing ($\tau_z \sim 1/F$) and strong forcing ($\tau_z$ independent of $F$) regimes. The crossover region is possibly characterized by a non-trivial scaling in $F$, matching the prediction of recent theories of polymer translocation. Although the geometrical setup is different, zipping and translocation share thus the same type of anomalous dynamics. Systems where this dynamics could be experimentally investigated are DNA (or RNA) hairpins: our results imply an anomalous dynamics for the hairpins closing times, but not for the opening times., Comment: 15 pages, 9 figures

We study a model of "elastic" lattice polymer in which a fixed number of monomers $m$ is hosted by a self-avoiding walk with fluctuating length $l$. We show that the stored length density $\rho_m = 1 - /m$ scales asymptotically for large $m$ as $\rho_m=\rho_\infty(1-\theta/m + ...)$, where $\theta$ is the polymer entropic exponent, so that $\theta$ can be determined from the analysis of $\rho_m$. We perform simulations for elastic lattice polymer loops with various sizes and knots, in which we measure $\rho_m$. The resulting estimates support the hypothesis that the exponent $\theta$ is determined only by the number of prime knots and not by their type. However, if knots are present, we observe strong corrections to scaling, which help to understand how an entropic competition between knots is affected by the finite length of the chain., Comment: 10 pages

We investigate the Rubinstein-Duke model for polymer reptation by means of density-matrix renormalization group techniques both in absence and presence of a driving field. In the former case the renewal time \tau and the diffusion coefficient D are calculated for chains up to N=150 reptons and their scaling behavior in N is analyzed. Both quantities scale as powers of N: $\tau \sim N^z$ and $D \sim 1/N^x$ with the asymptotic exponents z=3 and x=2, in agreement with the reptation theory. For an intermediate range of lengths, however, the data are well-fitted by some effective exponents whose values are quite sensitive to the dynamics of the end reptons. We find 2.7