Stretching dynamics of semiflexible polymers

Geometry-induced patterns through mechanochemical coupling

Intracellular protein patterns regulate a variety of vital cellular processes such as cell division and motility, which often involve dynamic cell-shape changes. These changes in cell shape may in turn affect the dynamics of pattern-forming proteins, hence leading to an intricate feedback loop between cell shape and chemical dynamics. While several computational studies have examined the rich resulting dynamics, the underlying mechanisms are not yet fully understood. To elucidate some of these mechanisms, we explore a conceptual model for cell polarity on a dynamic one-dimensional manifold. Using concepts from differential geometry, we derive the equations governing mass-conserving reaction-diffusion systems on time-evolving manifolds. Analyzing these equations mathematically, we show that dynamic shape changes of the membrane can induce pattern-forming instabilities in parts of the membrane, which we refer to as regional instabilities. Deformations of the local membrane geometry can also (regionally) suppress pattern formation and spatially shift already existing patterns. We explain our findings by applying and generalizing the local equilibria theory of mass-conserving reaction-diffusion systems. This allows us to determine a simple onset criterion for geometry-induced pattern-forming instabilities, which is linked to the phase-space structure of the reaction-diffusion system. The feedback loop between membrane shape deformations and reaction-diffusion dynamics then leads to a surprisingly rich phenomenology of patterns, including oscillations, traveling waves, and standing waves, even if these patterns do not occur in systems with a fixed membrane shape. Our paper reveals that the local conformation of the membrane geometry acts as an important dynamical control parameter for pattern formation in mass-conserving reaction-diffusion systems.

A dissipative particle dynamics study of a flexible filament in confined shear flow

In this paper, we investigate the dynamics of a tethered flexible filament due to fluid flow inside a microchannel. We use the finite sized dissipative particle dynamics (FDPD) approach to model this problem. The flexible filament is modeled as a bead-spring system with both extensional and flexural rigidity. The influence of flow rate and bending stiffness on the filament dynamics is studied in terms of the different conformational modes obtained. The competing effects of the hydrodynamic force and elastic force in the presence of Brownian thermal effects of comparable order influence the mode shapes of the filament. The dynamics of the filament motions are further analyzed using proper orthogonal decomposition. An important consequence of the dynamics of the filament is that it causes cross-flow in the micro-channel, which could potentially be exploited in micro-mixing and pumping applications. The cross stream fluid transport is observed to be more pronounced for higher bending stiffness.

Efficient simulation of semiflexible polymers

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 modeled with one bead-spring element per base pair, 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 polymer's 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 ≈1.6, 8, 16, and 16 ps for a dsDNA model chain consisting of 64, 128, 256, and 512 base pairs (i.e., chains of 0.55, 1.11, 2.24, and 4.48 persistence lengths), respectively. Correspondingly, in 1 h, 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 discretized) wormlike chain (WLC) model. Importantly, we demonstrate that at the same level of discretization, i.e., when each discretization element is one base pair long, our algorithm gains about five to six 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 discretized 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 behavior of a dsDNA segment in a shear flow.

Effects of thermal noise on the transitional dynamics of an inextensible elastic filament in stagnation flow

We investigate the dynamics of a single inextensible elastic filament subject to anisotropic friction in a viscous stagnation-point flow, by employing both a continuum model represented by Langevin type stochastic partial differential equations (SPDEs) and a dissipative particle dynamics (DPD) method. Unlike previous works, the filament is free to rotate and the tension along the filament is determined by the local inextensible constraint. The kinematics of the filament is recorded and studied with normal modes analysis. The results show that the filament displays an instability induced by negative tension, which is analogous to Euler buckling of a beam. Symmetry breaking of normal modes dynamics and stretch-coil transitions are observed above the threshold of the buckling instability point. Furthermore, both temporal and spatial noise are amplified resulting from the interaction of thermal fluctuations and nonlinear filament dynamics. Specifically, the spatial noise is amplified with even normal modes being excited due to symmetry breaking, while the temporal noise is amplified with increasing time correlation length and variance.

Dynamics of a semiflexible polymer or polymer ring in shear flow

Polymers exposed to shear flow exhibit a remarkably rich tumbling dynamics. While rigid rods rotate on Jeffery orbits, a flexible polymer stretches and coils up during tumbling. Theoretical results show that in both of these asymptotic regimes the corresponding tumbling frequency f(c) in a linear shear flow of strength γ scales as a power law Wi(2/3) in the Weissenberg number Wi = γτ, where τ is a characteristic time of the polymer's relaxational dynamics. For a flexible polymer these theoretical results are well confirmed by a large body of experimental single molecule studies. However, for the intermediate semiflexible regime, especially relevant for cytoskeletal filaments like F-actin and microtubules, the situation is less clear. While recent experiments on single F-actin filaments are still interpreted within the classical Wi(2/3) scaling law, theoretical results indicated deviations from it. Here we perform extensive Brownian dynamics simulations to explore the tumbling dynamics of semiflexible polymers over a broad range of shear strength and the polymer's persistence length l(p). We find that the Weissenberg number alone does not suffice to fully characterize the tumbling dynamics, and the classical scaling law breaks down. Instead, both the polymer's stiffness and the shear rate are relevant control parameters. Based on our Brownian dynamics simulations we postulate that in the parameter range most relevant for cytoskeletal filaments there is a distinct scaling behavior with f(c) τ* = Wi(3/4)f(c)(x) with Wi = γτ* and the scaling variable x = (l(p)/L)(Wi)(-1/3); here τ* is the time the polymer's center of mass requires to diffuse its own contour length L. Comparing these results with experimental data on F-actin we find that the Wi(3/4) scaling law agrees quantitatively significantly better with the data than the classical Wi(2/3) law. Finally, we extend our results to single ring polymers in shear flow, and find similar results as for linear polymers with slightly different power laws.

How fast does a signal propagate through proteins?

As the molecular basis of signal propagation in the cell, proteins are regulated by perturbations, such as mechanical forces or ligand binding. The question arises how fast such a signal propagates through the protein molecular scaffold. As a first step, we have investigated numerically the dynamics of force propagation through a single (Ala) protein following a sudden increase in the stretching forces applied to its end termini. The force propagates along the backbone into the center of the chain on the picosecond scale. Both conformational and tension dynamics are found in good agreement with a coarse-grained theory of force propagation through semiflexible polymers. The speed of force propagation of 50Å ps⁻¹ derived from these simulations is likely to determine an upper speed limit of mechanical signal transfer in allosteric proteins or molecular machines.

Propagation of tension along a polymer chain

We study the propagation of tension caused by an external force along a long polymeric molecule in two different settings, namely along a free polymer in three-dimensional (3D) space being pulled from one end, and along a prestretched circular polymer, confined in a narrow circular tube. We show that in both cases, the tension propagation is governed by a diffusion equation, and in particular, the tension front propagates as t(1/2) along the contour of the chain. The results are confirmed numerically, and by molecular dynamics simulations in the case of the 3D polymer. We also compare our results with the previously suggested ones for the translocation setting, and discuss why tension propagation is slower in that case.

Longitudinal response of confined semiflexible polymers

The longitudinal response of single semiflexible polymers to sudden changes in externally applied forces is known to be controlled by the propagation and relaxation of backbone tension. Under many experimental circumstances, realized, for example, in nanofluidic devices or in polymeric networks or solutions, these polymers are effectively confined in a channel- or tubelike geometry. By means of heuristic scaling laws and rigorous analytical theory, we analyze the tension dynamics of confined semiflexible polymers for various generic experimental setups. It turns out that in contrast to the well-known linear response, the influence of confinement on the nonlinear dynamics can largely be described as that of an effective prestress. We also study the free relaxation of an initially confined chain, finding a surprising superlinear ~t(9/8) growth law for the change in end-to-end distance at short times.

Tension dynamics and viscoelasticity of extensible wormlike chains

The dynamic response of prestressed semiflexible biopolymers is characterized by the propagation and relaxation of tension, which arises due to the near inextensibility of a stiff backbone. It is coupled to the dynamics of contour length stored in thermal undulations but also to the local relaxation of elongational strain. We present a systematic theory of tension dynamics for stiff yet extensible wormlike chains. Our results show that even moderate prestress gives rise to distinct Rouse-like extensibility signatures in the high-frequency viscoelastic response.

Freely relaxing polymers remember how they were straightened

The relaxation of initially straight semiflexible polymers has been discussed mainly with respect to the longest relaxation time. The biologically relevant nonequilibrium dynamics on shorter times is comparatively poorly understood, partly because "initially straight" can be realized in manifold ways. Combining Brownian dynamics simulations and systematic theory, we demonstrate how different experimental preparations give rise to specific short-time and universal long-time dynamics. We also discuss boundary effects and the onset of the stretch-coil transition.

Power-law creep behavior of a semiflexible chain

Rheological properties of adherent cells are essential for their physiological functions, and microrheological measurements on living cells have shown that their viscoelastic responses follow a weak power law over a wide range of time scales. This power law is also influenced by mechanical prestress borne by the cytoskeleton, suggesting that cytoskeletal prestress determines the cell's viscoelasticity, but the biophysical origins of this behavior are largely unknown. We have recently developed a stochastic two-dimensional model of an elastically joined chain that links the power-law rheology to the prestress. Here we use a similar approach to study the creep response of a prestressed three-dimensional elastically jointed chain as a viscoelastic model of semiflexible polymers that comprise the prestressed cytoskeletal lattice. Using a Monte Carlo based algorithm, we show that numerical simulations of the chain's creep behavior closely correspond to the behavior observed experimentally in living cells. The power-law creep behavior results from a finite-speed propagation of free energy from the chain's end points toward the center of the chain in response to an externally applied stretching force. The property that links the power law to the prestress is the chain's stiffening with increasing prestress, which originates from entropic and enthalpic contributions. These results indicate that the essential features of cellular rheology can be explained by the viscoelastic behaviors of individual semiflexible polymers of the cytoskeleton.