Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA model for the dynamics of extensible semiflexible polymers
182
0
0.0
(
0
)
Authors
Gerard T. Barkema
Ask ChatGPT about the research
No Arabic abstract
We present a model for semiflexible polymers in Hamiltonian formulation which interpolates between a Rouse chain and worm-like chain. Both models are realized as limits for the parameters. The model parameters can also be chosen to match the experimental force-extension curve for double-stranded DNA. Near the ground state of the Hamiltonian, the eigenvalues for the longitudinal (stretching) and the transversal (bending) modes of a chain with N springs, indexed by p, scale as lambda_lp ~ (p/N)^2 and lambda_tp ~ p^2(p-1)^2/N^4 respectively for small p. We also show that the associated decay times tau_p ~ (N/p)^4 will not be observed if they exceed the orientational time scale tau_r ~ N^3 for an equally-long rigid rod, as the driven decay is then washed out by diffusive motion.
rate research
Read More
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
We study the dynamical properties of semiflexible polymers with a recently introduced bead-spring model. We focus on double-stranded DNA. The two parameters of the model, $T^*$ and $ u$, are chosen to match its experimental force-extension curve. The bead-spring Hamiltonian is approximated in the first order by the Hessian that is quadratic in the bead positions. The eigenmodels of the Hessian provide the longitudinal (stretching) and transverse (bending) eigenmodes of the polymer, and the corresponding eigenvalues match well with the established phenomenology of semiflexible polymers. Using the longitudinal and transverse eigenmodes, we obtain analytical expressions of (i) the autocorrelation function of the end-to-end vector, (ii) the autocorrelation function of a bond (i.e., a spring, or a tangent) vector at the middle of the chain, and (iii) the mean-square displacement of a tagged bead in the middle of the chain, as sum over the contributions from the modes. We also perform simulations with the full dynamics of the model. The simulations yield numerical values of the correlation functions (i-iii) that agree very well with the analytical expressions for the linearized dynamics. We also study the mean-square displacement of the longitudinal component of the end-to-end vector that showcases strong nonlinear effects in the polymer dynamics, and we identify at least an effective $t^{7/8}$ power-law regime in its time-dependence. Nevertheless, in comparison to the full mean-square displacement of the end-to-end vector the nonlinear effects remain small at all times --- it is in this sense we state that our results demonstrate that the linearized dynamics suffices for dsDNA fragments that are shorter than or comparable to the persistence length. Our results are consistent with those of the wormlike chain (WLC) model, the commonly used descriptive tool of semiflexible polymers.
We explore the effect of an attractive interaction between parallel-aligned polymers, which are perpendicularly grafted on a substrate. Such an attractive interaction could be due to, e.g., reversible cross-links. The competition between permanent grafting favoring a homogeneous state of the polymer brush and the attraction, which tends to induce in-plane collapse of the aligned polymers, gives rise to an instability of the homogeneous phase to a bundled state. In this latter state the in-plane translational symmetry is spontaneously broken and the density is modulated with a finite wavelength, which is set by the length scale of transverse fluctuations of the grafted polymers. We analyze the instability for two models of aligned polymers: directed polymers with a line tension and weakly bending chains with a bending stiffness.
It has become clear in recent years that the simple uniform wormlike chain model needs to be modified in order to account for more complex behavior which has been observed experimentally in some important biopolymers. For example, the large flexibility of short ds-DNA has been attributed to kink or hinge defects. In this paper, we calculate analytically, within the weak bending approximation, the force-extension relation of a wormlike chain with a permanent hinge defect along its contour. The defect is characterized by its bending energy (which can be zero, in the completely flexible case) and its position along the polymer contour. Besides the bending rigidity of the chain, these are the only parameters which describe our model. We show that a hinge defect causes a significant increase in the differential tensile compliance of a pre-stressed chain. In the small force limit, a hinge defect significantly increases the entropic elasticity. Our results apply to any pair of semiflexible segments connected by a hinge. As such, they may also be relevant to cytoskeletal filaments (F-actin, microtubules), where one may treat the cross-link connecting two filaments as a hinge defect.
Semiflexible polymer models are widely used as a paradigm to understand structural phases in biomolecules including folding of proteins. Since stable knots are not so common in real proteins, the existence of stable knots in semiflexible polymers has not been explored much. Here, via extensive replica exchange Monte Carlo simulation we investigate the same for a bead-stick and a bead-spring homopolymer model that covers the whole range from flexible to stiff. We establish the fact that the presence of stable knotted phases in the phase diagram is dependent on the ratio $r_b/r_{rm{min}}$ where $r_b$ is the equilibrium bond length and $r_{rm{min}}$ is the distance for the strongest nonbonded contacts. Our results provide evidence for both models that if the ratio $r_b/r_{rm{min}}$ is outside a small window around unity then depending on the bending stiffness one always encounters stable knotted phases along with the usual frozen and bent-like structures at low temperatures. These findings prompt us to conclude that knots are generic stable phases in semiflexible polymers.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
