Buckling a Semiflexible Polymer Chain under Compression

Nucleosome Array Deformation in Chromatin is Sustained by Bending, Twisting and Kinking of Linker DNA

Chromatin in the nucleus undergoes mechanical stresses from different sources during the various stages of cell life. Here a trinucleosome array is used as the minimal model to study the mechanical response to applied stress at the molecular level. By using large-scale, all-atom steered-molecular dynamics simulations, we show that the largest part of mechanical stress in compression is accommodated by the DNA linkers joining pairs of nucleosomes, which store the elastic energy accumulated by the applied force. Different mechanical instabilities (Euler bending, Brazier kinking, twist-bending) can deform the DNA canonical structure, as a function of the increasing force load. An important role of the histone tails in assisting the DNA deformation is highlighted. The overall response of the smallest chromatin fragment to compressive stress leaves the nucleosome assembly with a substantial plastic deformation and localised defects, which can have a potential impact on DNA transcription, downstream signaling pathways, the regulation of gene expression, and DNA repair.

Compression-induced buckling of a semiflexible filament in two and three dimensions

The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. For filaments in the low-temperature limit and under a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression. We determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates this transition. We find the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble, the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show the theory is highly sensitive to fluctuations in length in two dimensions, and that the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts.

Characterizing the free-energy landscapes of DNA origamis

We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.

Instability caused swimming of ferromagnetic filaments in pulsed field

Magnetic filaments driven by external magnetic field are an interesting topic of research in-terms of the possible bio-medical applications. In this paper, we investigated the applicability of using ferromagnetic filaments as micro swimmers both experimentally and numerically. It was found that applying a pulse wave field profile with a duty cycle of 30[Formula: see text] induced experimentally observable swimming, which is similar to the breast stroke of micro algae. Good agreement with numerical simulations was found. Moreover, for stable continuous swimming, an initial filament shape is required to avoid transition to the structurally preferred non-swimming S-like mode.

Compression and Stretching of Confined Linear and Ring Polymers by Applying Force

We use Langevin dynamics to study the deformations of linear and ring polymers in different confinements by applying compression and stretching forces on their two sides. Our results show that the compression deformations are the results of an interplay among of polymer rigidity, degree of confinement, and force applied. When the applied force is beyond the threshold required for the buckling transition, the semiflexible chain under the strong confinement firstly buckles; then comes helical deformation. However, under the same force loading, the semiflexible chain under the weaker confinement exhibits buckling instability and shrinks from the folded ends/sides until it becomes three-folded structures. This happens because the strong confinement not only strongly reduces the buckling wavelength, but also increases the critical buckling force threshold. For the weakly confined polymers, in compression process, the flexible linear polymer collapses into condensed states under a small external force, whereas the ring polymer only shows slight shrinkage, due to the excluded volume interactions of two strands in the crowded states. These results are essential for understanding the deformations of the ring biomacromolecules and polymer chains in mechanical compression or driven transport.

Loops versus lines and the compression stiffening of cells

Both animal and plant tissue exhibit a nonlinear rheological phenomenon known as compression stiffening, or an increase in moduli with increasing uniaxial compressive strain. Does such a phenomenon exist in single cells, which are the building blocks of tissues? One expects an individual cell to compression soften since the semiflexible biopolymer-based cytoskeletal network maintains the mechanical integrity of the cell and in vitro semiflexible biopolymer networks typically compression soften. To the contrary, we find that mouse embryonic fibroblasts (mEFs) compression stiffen under uniaxial compression via atomic force microscopy studies. To understand this finding, we uncover several potential mechanisms for compression stiffening. First, we study a single semiflexible polymer loop modeling the actomyosin cortex enclosing a viscous medium modeled as an incompressible fluid. Second, we study a two-dimensional semiflexible polymer/fiber network interspersed with area-conserving loops, which are a proxy for vesicles and fluid-based organelles. Third, we study two-dimensional fiber networks with angular-constraining crosslinks, i.e. semiflexible loops on the mesh scale. In the latter two cases, the loops act as geometric constraints on the fiber network to help stiffen it via increased angular interactions. We find that the single semiflexible polymer loop model agrees well with the experimental cell compression stiffening finding until approximately 35% compressive strain after which bulk fiber network effects may contribute. We also find for the fiber network with area-conserving loops model that the stress-strain curves are sensitive to the packing fraction and size distribution of the area-conserving loops, thereby creating a mechanical fingerprint across different cell types. Finally, we make comparisons between this model and experiments on fibrin networks interlaced with beads as well as discuss implications for single cell compression stiffening at the tissue scale.

Force-displacement relations at compression of dsDNA macromolecules

The elasticity of dsDNA molecules is investigated by Monte Carlo simulations based on a coarse-grained model of DNA. The force-displacement (f-r) curves are computed under the constraints of the constant force (Gibbs) or the constant length (Helmholtz) ensemble. Particular attention was paid to the compressional (negative) and weak tensile forces. It was confirmed that simulations using the vector Gibbs ensemble fail to represent the compression behavior of polymers. Simulations using the scalar Gibbs protocol resulted in a qualitatively correct compressional response of DNA provided that the quadratic averages of displacements were employed. Furthermore, a well-known shortcoming of the popular Marko-Siggia relation for DNA elasticity at weak tensile forces is elucidated. Conversely, the function f-r from the simulation at the constant length constraint, as well as the new closed-form expressions, provides a realistic depiction of the DNA elasticity over the wide range of negative and positive forces. Merely a qualitative resemblance of the compression functions f-r predicted by the employed approaches supports the notion that the elastic response of DNA molecules may be greatly affected by the specifics of the experimental setups and the kind of averaging of the measured variable.

Elastic behavior of a semiflexible polymer in 3D subject to compression and stretching forces

We elucidate the elastic behavior of a wormlike chain in 3D under compression and provide exact solutions for the experimentally accessible force-extension relation in terms of generalized spheroidal wave functions. In striking contrast to the classical Euler buckling instability, the force-extension relation of a clamped semiflexible polymer exhibits a smooth crossover from an almost stretched to a buckled configuration. In particular, the associated susceptibility, which measures the strength of the response of the polymer to the applied force, displays a prominent peak in the vicinity of the critical Euler buckling force. For increasing persistence length, the force-extension relation and the susceptibility of semiflexible polymers approach the behavior of a classical rod, whereas thermal fluctuations permit more flexible polymers to resist the applied force. Furthermore, we find that semiflexible polymers confined to 2D can oppose the applied force more strongly than in 3D.

Unfolding Kinetics of a Wormlike Chain under Elongational Flow

A simple theory of the unfolding kinetics of a semi-flexible polymer chain is presented in terms of a Kramers type picture for the energy of elongation. The hydrodynamic interactions are discussed in terms of slender body theory. It turns out that the elongation of the chain is basically linear in time and independent of the viscosity. The former prediction agrees with experiments on the stretching dynamics of DNA under planar elongational flow. Nevertheless, the theory overestimates the experimental rate by a significant amount for reasons that are unclear.