Source author record

Enrico Carlon

Enrico Carlon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Equilibrium Fluctuations of DNA Plectonemes

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.

preprint2022arXiv

Mechanical properties of Nucleic Acids and the non-local Twistable Wormlike Chain model

Mechanical properties of nucleic acids play an important role in many biological processes which often involve physical deformations of these molecules. At sufficiently long length scales (say above $\sim 20-30$ base pairs) the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites gets coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model which clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables which describe the twist and bending of the molecule. Differently from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior which matches experimental data.

preprint2022arXiv

Rigid Base Biasing in Molecular Dynamics enables enhanced sampling of DNA conformations

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.

preprint2016arXiv

Torque-Induced Rotational Dynamics in Polymers: Torsional Blobs and Thinning

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.

preprint2014arXiv

Models of DNA denaturation dynamics: universal properties

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 $τ\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.

preprint2014arXiv

Rotational dynamics of entangled polymers

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.

preprint2014arXiv

Unwinding Dynamics of a Helically Wrapped Polymer

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.

preprint2013arXiv

Unwinding relaxation dynamics of polymers

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.

preprint2012arXiv

Fractional Brownian motion and the critical dynamics of zipping polymers

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 $τ\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.

preprint2012arXiv

Probing Hybridization parameters from microarray experiments: nearest neighbor model and beyond

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.

preprint2011arXiv

Equilibrium winding angle of a polymer around a bar

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 $<θ^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 $<θ^2>\sim (\ln L)^{2α}$, with $α=0.75(1)$. The ratio $γ= <θ^4>/<θ^2>^2=3.74(5)$ is incompatible with the gaussian value $γ=3$, but consistent with the observation that the tail of the probability distribution function $p(θ)$ is found to decrease slower than a gaussian function. These findings are at odds with the existing first-order $ε$-expansion results.

preprint2011arXiv

Nonequilibrium effects in DNA microarrays: a multiplatform study

It has recently been shown that in some DNA microarrays the time needed to reach thermal equilibrium may largely exceed the typical experimental time, which is about 15h in standard protocols (Hooyberghs et al. Phys. Rev. E 81, 012901 (2010)). In this paper we discuss how this breakdown of thermodynamic equilibrium could be detected in microarray experiments without resorting to real time hybridization data, which are difficult to implement in standard experimental conditions. The method is based on the analysis of the distribution of fluorescence intensities I from different spots for probes carrying base mismatches. In thermal equilibrium and at sufficiently low concentrations, log I is expected to be linearly related to the hybridization free energy $ΔG$ with a slope equal to $1/RT_{exp}$, where $T_{exp}$ is the experimental temperature and R is the gas constant. The breakdown of equilibrium results in the deviation from this law. A model for hybridization kinetics explaining the observed experimental behavior is discussed, the so-called 3-state model. It predicts that deviations from equilibrium yield a proportionality of $\log I$ to $ΔG/RT_{eff}$. Here, $T_{eff}$ is an effective temperature, higher than the experimental one. This behavior is indeed observed in some experiments on Agilent arrays. We analyze experimental data from two other microarray platforms and discuss, on the basis of the results, the attainment of equilibrium in these cases. Interestingly, the same 3-state model predicts a (dynamical) saturation of the signal at values below the expected one at equilibrium.

preprint2010arXiv

Elastic Lattice Polymers

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 $ρ_m = 1 - <l>/m$ scales asymptotically for large $m$ as $ρ_m=ρ_\infty(1-θ/m + ...)$, where $θ$ is the polymer entropic exponent, so that $θ$ can be determined from the analysis of $ρ_m$. We perform simulations for elastic lattice polymer loops with various sizes and knots, in which we measure $ρ_m$. The resulting estimates support the hypothesis that the exponent $θ$ 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.

preprint2002arXiv

A simple model of DNA denaturation and mutually avoiding walks statistics

Recently Garel, Monthus and Orland (Europhys. Lett. v 55, 132 (2001)) considered a model of DNA denaturation in which excluded volume effects within each strand are neglected, while mutual avoidance is included. Using an approximate scheme they found a first order denaturation. We show that a first order transition for this model follows from exact results for the statistics of two mutually avoiding random walks, whose reunion exponent is c > 2, both in two and three dimensions. Analytical estimates of c due to the interactions with other denaturated loops, as well as numerical calculations, indicate that the transition is even sharper than in models where excluded volume effects are fully incorporated. The probability distribution of distances between homologous base pairs decays as a power law at the transition.

preprint1999arXiv

Density Matrix Renormalization Group and Reaction-Diffusion Processes

The density matrix renormalization group (DMRG) is applied to some one-dimensional reaction-diffusion models in the vicinity of and at their critical point. The stochastic time evolution for these models is given in terms of a non-symmetric ``quantum Hamiltonian'', which is diagonalized using the DMRG method for open chains of moderate lengths (up to about 60 sites). The numerical diagonalization methods for non-symmetric matrices are reviewed. Different choices for an appropriate density matrix in the non-symmetric DMRG are discussed. Accurate estimates of the steady-state critical points and exponents can then be found from finite-size scaling through standard finite-lattice extrapolation methods. This is exemplified by studying the leading relaxation time and the density profiles of diffusion-annihilation and of a branching-fusing model in the directed percolation universality class.