Abstract
Understanding looping probabilities, including the particular case of ring closure or cyclization, of fluctuating polymers (e.g., DNA) is important in many applications in molecular biology and chemistry. In a continuum limit the configuration of a polymer is a curve in the group $\mathit{SE}$(3) of rigid body displacements, whose energy can be modeled via the Cosserat theory of elastic rods. Cosserat rods are a more detailed version of the classic wormlike-chain (WLC) model, which we show to be more appropriate in short-length scale, or stiff, regimes, where the contributions of extension and shear deformations are not negligible and lead to noteworthy high values for the cyclization probabilities (or $J$-factors). We therefore observe that the Cosserat framework is a candidate for gaining a better understanding of the enhanced cyclization of short DNA molecules reported in various experiments, which is not satisfactorily explained by WLC-type models. Characterizing the stochastic fluctuations about minimizers of the energy by means of Laplace expansions in a (real) path integral formulation, we develop efficient analytical approximations for the two cases of full looping, in which both end-to-end relative translation and rotation are prescribed, and of marginal looping probabilities, where only end-to-end translation is prescribed. For isotropic Cosserat rods, certain looping boundary value problems admit nonisolated families of critical points of the energy due to an associated continuous symmetry. For the first time, taking inspiration from (imaginary) path integral techniques, a quantum mechanical probabilistic treatment of Goldstone modes in statistical rod mechanics sheds light on $J$-factor computations for isotropic rods in the semiclassical context. All the results are achieved exploiting appropriate Jacobi fields arising from Gaussian path integrals and show good agreement when compared with intense Monte Carlo simulations for the target examples.
Highlights
It is widely known that polymers involved in biological and chemical processes are anything but static objects
The present study is general and is applicable to various end-to-end statistics, we focus on the computation of ring-closure or cyclization probabilities for elastic rods, targeting three significant aspects
The questions we are trying to answer would be interpreted as follows: what is a good estimate of the probability of the end monomers coming into contact with each other? How is the latter value changing if we impose an orientation constraint on the binding site? How does the shape of the cross section affect the statistics? And, what happens if we deviate from the standard inextensible and unshearable model and incorporate shear and extension as possible deformations?
Summary
It is widely known that polymers involved in biological and chemical processes are anything but static objects. This is a prerogative of the Cosserat, more general framework, where the centerline displacement and the cross-sectional rotation are considered as independent variables We show that these additional degrees of freedom are crucial in the analysis of polymer chains in short-length scale, or stiff, regimes, in both the full and marginal cases, where the system exploits extension and shear deformations for minimizing the overall elastic energy, in the face of an increasingly penalizing bending contribution. This allows the cyclization probability density to take high values even when the WLC model (and Kirchhoff) is vanishing exponentially.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.