Semiflexible polymers: dependence on ensemble and boundary orientations

Forced extension of a wormlike chain in the Gibbs and Helmholtz ensembles

A semiflexible polymer can be stretched by either applying a force to it or by fixing the positions of its endpoints. The two approaches generally yield different results and correspond to experiments performed in either the Gibbs or Helmholtz statistical ensembles. Here, we derive the Helmholtz force-extension relationship for the commonly used wormlike-chain model in the strongly stretched regime. By analyzing it in comparison with the Gibbs ensemble result, we show that equivalence between the two relationships is achieved only in the long-chain thermodynamic limit.

Molecular Dynamics Simulation of a Feather-Boa Model of a Bacterial Chromosome

The chromosome of a bacterium consists of a mega-base pair-long circular DNA, which self-organizes within the micron-sized bacterial cell volume, compacting itself by three orders of magnitude. Unlike eukaryotic chromosomes, it lacks a nuclear membrane and freely floats in the cytosol confined by the cell membrane. It is believed that strong confinement, cross-linking by associated proteins, and molecular crowding all contribute to determine chromosome size and morphology. Modelling the chromosome simply as a circular polymer decorated with closed side loops in a cylindrical confining volume has been shown to already recapture some of the salient properties observed experimentally. Here we describe how a computer simulation can be set up to study structure and dynamics of bacterial chromosomes using this model.

Grafted Semiflexible Nunchucks with a Magnetic Bead Attached to the Free End

Semiflexible nunchucks are block copolymers, which consist of two long blocks of high bending stiffness jointed together by a short block of low bending stiffness. Semiflexible nunchucks that consist of two DNA nanorods jointed by a short segment of double-stranded (ds) DNA and confined in two dimensions have been used in recent experiments by Fygenson and coworkers as a tool to magnify the bending fluctuations of the linking dsDNA, which in turn are used to deduce the persistence length of dsDNA. In a recent theoretical analysis, we showed that in a semiflexible nunchuck with one end grafted, the fluctuations of the position of the free end that is transverse to the grafting direction exhibit a pronounced bimodality, provided that the bending stiffness of the hinge is not very large. In this article, we theoretically analyse a grafted semiflexible nunchuck with a magnetic bead attached to its free end. We show that a transverse magnetic field induces an asymmetry in the bimodal distribution of the transverse fluctuations of the free end. This asymmetry is very sensitive to interactions with a magnetic field and, in principle, could be used in magnetometry (the measurement of a magnetic field or the magnetic moment of the bead). We also investigate how the response of the bimodal distribution of the transverse fluctuations of the free end to a magnetic field depends on the bending stiffness of the nunchuck hinge. In addition, we analyse the closely related systems of a single filament and two filaments jointed at a kink point with one end grafted and the other end attached to a magnetic bead.

Orientational Fluctuations and Bimodality in Semiflexible Nunchucks

Semiflexible nunchucks are block copolymers consisting of two long blocks with high bending rigidity jointed by a short block of lower bending stiffness. Recently, the DNA nanotube nunchuck was introduced as a simple nanoinstrument that mechanically magnifies the bending angle of short double-stranded (ds) DNA and allows its measurement in a straightforward way [Fygenson et al., Nano Lett. 2020, 20, 2, 1388-1395]. It comprises two long DNA nanotubes linked by a dsDNA segment, which acts as a hinge. The semiflexible nunchuck geometry also appears in dsDNA with a hinge defect (e.g., a quenched denaturation bubble or a nick), and in end-linked stiff filaments. In this article, we theoretically investigate various aspects of the conformations and the tensile elasticity of semiflexible nunchucks. We analytically calculate the distribution of bending fluctuations of a wormlike chain (WLC) consisting of three blocks with different bending stiffness. For a system of two weakly bending WLCs end-jointed by a rigid kink, with one end grafted, we calculate the distribution of positional fluctuations of the free end. For a system of two weakly bending WLCs end-jointed by a hinge modeled as harmonic bending spring, with one end grafted, we calculate the positional fluctuations of the free end. We show that, under certain conditions, there is a pronounced bimodality in the transverse fluctuations of the free end. For a semiflexible nunchuck under tension, under certain conditions, there is bimodality in the extension as a function of the hinge position. We also show how steric repulsion affects the bending fluctuations of a rigid-rod nunchuck.

A semiflexible polymer in a gliding assay: reentrant transition, role of turnover and activity

We consider a model of an extensible semiflexible filament moving in two dimensions on a motility assay of motor proteins represented explicitly as active harmonic linkers. Their heads bind stochastically to polymer segments within a capture radius, and extend along the filament in a directed fashion before detaching. Both the extension and detachment rates are load-dependent and generate an active drive on the filament. The filament undergoes a first order phase transition from the open chain to spiral conformation and shows a reentrant behavior in both the active extension and the turnover, defined as the ratio of attachment-detachment rates. Associated with the phase transition, the size and shape of the polymer change non-monotonically, and the relevant autocorrelation functions display a double-exponential decay. The corresponding correlation times show a maximum signifying the dominance of spirals. The orientational dynamics captures the rotation of spirals, and its correlation time decays with activity as a power law.

Active Brownian particles: mapping to equilibrium polymers and exact computation of moments

It is well known that the path probabilities of Brownian motion correspond to the equilibrium configurational probabilities of flexible Gaussian polymers, while those of active Brownian motion correspond to in-extensible semiflexible polymers. Here we investigate the properties of the equilibrium polymer that corresponds to the trajectories of particles acted on simultaneously by both Brownian and active noise. Through this mapping we can see interesting crossovers in the mechanical properties of the polymer with changing contour length. The polymer end-to-end distribution exhibits Gaussian behaviour for short lengths, which changes to the form of semiflexible filaments at intermediate lengths, to finally go back to a Gaussian form for long contour lengths. By performing a Laplace transform of the governing Fokker-Planck equation of the active Brownian particle, we discuss a direct method to derive exact expressions for all the moments of the relevant dynamical variables, in arbitrary dimensions. These are verified via numerical simulations and used to describe interesting qualitative features such as, for example, dynamical crossovers. Finally we discuss the kurtosis of the ABP's position, which we compute exactly, and show that it can be used to differentiate between active Brownian particles and the active Ornstein-Uhlenbeck process.

Statistical ensemble inequivalence for flexible polymers under confinement in various geometries

The problem of statistical ensemble inequivalence for single polymers has been the subject of intense research. In a recent publication, we show that even though the force-extension relation of a free Gaussian chain exhibits ensemble equivalence, confinement to half-space due to tethering to a planar substrate induces significant inequivalence [S. Dutta and P. Benetatos, Soft Matter, 2018, 14, 6857-6866]. In the present article, we extend that work to the conformational response to confining forces distributed over surfaces. We analyze in both the Helmholtz and the Gibbs ensemble the pressure-volume equation of state of a chain in rectangular, spherical, and cylindrical confinement. We especially consider the case of a directed polymer in a cylinder. We also analyze the case of a tethered chain inside a rectangular box, a sphere, and outside a sphere. In general, confinement causes significant ensemble inequivalence. Remarkably, we recover ensemble equivalence in the limit of squashing confinement. We trace the ensemble inequivalence to the persistence of strong fluctuations. Our work may be useful in the interpretation of single molecule experiments and caging phenomena.

Molecular Dynamics Simulation of a Feather-Boa Model of a Bacterial Chromosome

The chromosome of a bacterium consists of a mega-base pair long circular DNA, which self-organizes within the micron-sized bacterial cell volume, compacting itself by three orders of magnitude. Unlike in eukaryotes, it lacks a nuclear membrane, and freely floats in the cytosol confined by the cell membrane. It is believed that strong confinement, cross-linking by associated proteins, and molecular crowding all contribute to determine chromosome size and morphology. Modeling the chromosome simply as a circular polymer decorated with closed side-loops in a cylindrical confining volume, has been shown to already recapture some of the salient properties observed experimentally. Here, we describe how a computer simulation can be set up to study structure and dynamics of bacterial chromosomes using this model.

Morphological and dynamical properties of semiflexible filaments driven by molecular motors

We consider an explicit model of a semiflexible filament moving in two dimensions on a gliding assay of motor proteins, which attach to and detach from filament segments stochastically, with a detachment rate that depends on the local load experienced. Attached motor proteins move along the filament to one of its ends with a velocity that varies nonlinearly with the motor protein extension. The resultant force on the filament drives it out of equilibrium. The distance from equilibrium is reflected in the end-to-end distribution, modified bending stiffness, and a transition to spiral morphology of the polymer. The local stress dependence of activity results in correlated fluctuations in the speed and direction of the center of mass leading to a series of ballistic-diffusive crossovers in its dynamics.

Inequivalence of fixed-force and fixed-extension statistical ensembles for a flexible polymer tethered to a planar substrate

Recent advances in single macromolecule experiments have sparked interest in the ensemble dependence of force-extension relations. The thermodynamic limit may not be attainable for such systems, which leads to inequivalence of the fixed-force and the fixed-extension ensembles. We consider an ideal Gaussian chain described by the Edwards Hamiltonian with one end tethered to a rigid planar substrate. We analytically calculate the force-extension relation in the two ensembles and we show their inequivalence, which is caused by the confinement of the polymer to half space. The inequivalence is quite remarkable for strong compressional forces. We also perform Monte-Carlo simulations of a tethered wormlike chain with contour length 20 times its persistence length, which corresponds to experiments measuring the conformations of DNA tethered to a wall. The simulations confirm the ensemble inequivalence and qualitatively agree with the analytical predictions of the Gaussian model. Our analysis shows that confinement due to tethering causes ensemble inequivalence, irrespective of the polymer model.

Forced desorption of semiflexible polymers, adsorbed and driven by molecular motors

We formulate and characterize a model to describe the dynamics of semiflexible polymers in the presence of activity due to motor proteins attached irreversibly to a substrate, and a transverse pulling force acting on one end of the filament. The stochastic binding-unbinding of the motor proteins and their ability to move along the polymer generate active forces. As the pulling force reaches a threshold value, the polymer eventually desorbs from the substrate. Performing underdamped Langevin dynamics simulation of the polymer, and with stochastic motor activity, we obtain desorption phase diagrams. The correlation time for fluctuations in the desorbed fraction increases as one approaches complete desorption, captured quantitatively by a power law spectral density. We present theoretical analysis of the phase diagram using mean field approximations in the weakly bending limit of the polymer and performing linear stability analysis. This predicts an increase in the desorption force with the polymer bending rigidity, active velocity and processivity of the motor proteins to capture the main features of the simulation results.

Disordered, stretched, and semiflexible biopolymers in two dimensions

We study the effects of intrinsic sequence-dependent curvature for a two-dimensional semiflexible biopolymer with short-range correlation in intrinsic curvatures. We show exactly that when not subjected to any external force, such a system is equivalent to a system with a well-defined intrinsic curvature and a proper renormalized persistence length. We find the exact expression for the distribution function of the equivalent system. However, we show that such an equivalent system does not always exist for the polymer subjected to an external force. We find that under an external force, the effect of sequence disorder depends upon the averaging order, the degree of disorder, and the experimental conditions, such as the boundary conditions. Furthermore, a short to moderate length biopolymer may be much softer or has a smaller apparent persistent length than what would be expected from the "equivalent system." Moreover, under a strong stretching force and for a long biopolymer, the sequence disorder is immaterial for elasticity. Finally, the effect of sequence disorder may depend upon the quantity considered.

End-to-end distribution for a wormlike chain in arbitrary dimensions

We construct an efficient methodology for calculating wormlike chain statistics in arbitrary D dimensions over all chain rigidities, from fully rigid to completely flexible. The structure of our exact analytical solution for the end-to-end distribution function for a wormlike chain in arbitrary D dimensions in Fourier-Laplace space (i.e., Fourier-transformed end position and Laplace-transformed chain length) adopts the form of an infinite continued fraction, which is advantageous for its compact structure and stability for numerical implementation. We then proceed to present a step-by-step methodology for performing the Fourier-Laplace inversion in order to make full use of our results in general applications. Asymptotic methods for evaluating the Laplace inversion (power-law expansion and Rayleigh-Schrödinger perturbation theory) are employed in order to improve the accuracy of the numerical inversions of the end-to-end distribution function in real space. We adapt our results to the evaluation of the single-chain structure factor, rendering simple, closed-form expressions that facilitate comparison with scattering experiments. Using our techniques, the accuracy of the end-to-end distribution function is enhanced up to the limit of the machine precision. We demonstrate the utility of our methodology with realizations of the chain statistics, giving a general methodology that can be applied to a wide range of biophysical problems.

Unique elastic properties of the spectrin tetramer as revealed by multiscale coarse-grained modeling

The force-extension profile of tetrameric spectrin is determined by using multiscale computer simulation. Fluctuation results of atomistic simulations of double spectrin repeat units (DSRU) are used to systematically build a coarse-grained (CG) model for the tetrameric form of spectrin. It is found that the spectrin tetramer can be modeled as a soft polymer with a unique flat force-extension profile over the range of biologically important lengths. It is also concluded that in the cytoskeletal network of the red blood cell the tetramer is in an "overcompressed" state. These findings are in contrast to the commonly used models of spectrin tetramer elasticity, namely the "entropic spring" polymer models. From these results, it is concluded that stable intact helical linker regions are needed to maintain the soft elasticity of the spectrin tetramer.

Elasticity of two-dimensional filaments with constant spontaneous curvature

We study the mechanical property of a two-dimensional filament with constant spontaneous curvature and under uniaxial applied force. We derive the equation that governs the stable shape of the filament and obtain analytical solutions for the equation. We find that for a long filament with positive initial azimuth angle (the azimuth angle is the angle between x axis and the tangent of the filament) and under large stretching force, the azimuth angle is a two-valued function of the arclength, decreases first, and then increases with increasing arclength. Otherwise, the azimuth angle is a monotonic function of arclength. At finite temperature, we derive the differential equation that governs the partition function and find exact solution of the partition function for a filament free of force. We obtain closed-form expressions on the force-extension relation for a filament under low force and for a long filament under strong stretching force. Our results show that for a biopolymer with moderate length and not too small spontaneous curvature, the effect of the spontaneous curvature cannot be replaced by a simple renormalization of the persistence length in the wormlike chain model.

Bimodal distribution function of a three-dimensional wormlike chain with a fixed orientation of one end

We study the distribution function of a three-dimensional wormlike chain with a fixed orientation of one chain end using the exact representation of the distribution function in terms of the Green's function of the quantum rigid rotator in a homogeneous external field. The transverse one-dimensional distribution function of the free chain end displays a bimodal shape in the intermediate range of chain lengths (1.3L{p},...,3.5L{p}). We also present analytical results for short and long chains, which are in complete agreement with the results of previous studies obtained using different methods.