Elasticity of a semiflexible filament with a discontinuous tension due to a cross-link or a molecular motor

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.

Elasticity of connected semiflexible quadrilaterals

Using the positional-orientational propagator of a semiflexible filament in the weakly bending regime, we analytically calculate the probability densities associated with the fluctuating tip and the corners of a grafted system of connected quadrilaterals. We calculate closed analytic expressions for the probability densities within the framework of the worm-like chain model, which are valid in the weakly bending regime. The probability densities give the physical quantities related to the elasticity of the system such as the force-extension relation in the fixed extension ensemble, the Poisson's ratio and the average of the force exerted to a confining stiff planar wall by the fluctuating tip of the system. Our analysis reveals that the force-extension relations depend on the contour length of the system (material content), the bending stiffness (chemical nature), the geometrical angle and the number of the quadrilaterals, while the Poisson's ratio depends only on the geometrical angle and the number of the quadrilaterals, and is thus a purely geometric property of the system.

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.

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.

Bimodal probability density characterizes the elastic behavior of a semiflexible polymer in 2D under compression

We explore the elastic behavior of a wormlike chain under compression in terms of exact solutions for the associated probability densities. Strikingly, the probability density for the end-to-end distance projected along the applied force exhibits a bimodal shape in the vicinity of the critical Euler buckling force of an elastic rod, reminiscent of the smeared discontinuous phase transition of a finite system. These two modes reflect the almost stretched and the S-shaped configuration of a clamped polymer induced by the compression. Moreover, we find a bimodal shape of the probability density for the transverse fluctuations of the free end of a cantilevered polymer as fingerprint of its semiflexibility. In contrast to clamped polymers, free polymers display a circularly symmetric probability density and their distributions are identical for compression and stretching forces.

Tensile elasticity of semiflexible polymers with hinge defects

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

Exact solution for the force-extension relation of a semiflexible polymer under compression

Exact solutions for the elastic and thermodynamic properties for the wormlike chain model are elaborated in terms of Mathieu functions. The smearing of the classical Euler buckling instability for clamped polymers is analyzed for the force-extension relation. Interestingly, at strong compression forces the thermal fluctuations lead to larger elongations than for the elastic rod. The susceptibility defined as the derivative of the force-extension relation displays a prominent maximum at a force that approaches the critical Euler buckling force as the persistence length is increased. We also evaluate the excess entropy and heat capacity induced by the compression and find that they vary nonmonotonically with the load. These findings are corroborated by pseudo-Brownian simulations.