Researcher profile

Enrico Carlon

Enrico Carlon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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.

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 &#34;elastic&#34; 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.

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&#39;&#39;, 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.