An official website of the United States government
Motivated by the observation of the storage of excess elastic free energy -- prestress in cross linked semiflexible filament networks, we consider the problem of the conformational statistics of a single semiflexible polymer in a quenched random potential. The random potential, which represents the effect of cross linking to other filaments is assumed to have a finite correlation length and mean strength. We examine the statistical distribution of curvature in the limit that the filaments are much shorter than their thermal persistence length. We compare our theoretical predictions to finite element Brownian dynamics simulations. Lastly we comment on the validity of replica field techniques in addressing these questions.
more »
« less
Bending fluctuations in semiflexible, inextensible, slender filaments in Stokes flow: Toward a spectral discretization
Semiflexible slender filaments are ubiquitous in nature and cell biology, including in the cytoskeleton, where reorganization of actin filaments allows the cell to move and divide. Most methods for simulating semiflexible inextensible fibers/polymers are based on discrete (bead-link or blob-link) models, which become prohibitively expensive in the slender limit when hydrodynamics is accounted for. In this paper, we develop a novel coarse-grained approach for simulating fluctuating slender filaments with hydrodynamic interactions. Our approach is tailored to relatively stiff fibers whose persistence length is comparable to or larger than their length and is based on three major contributions. First, we discretize the filament centerline using a coarse non-uniform Chebyshev grid, on which we formulate a discrete constrained Gibbs–Boltzmann (GB) equilibrium distribution and overdamped Langevin equation for the evolution of unit-length tangent vectors. Second, we define the hydrodynamic mobility at each point on the filament as an integral of the Rotne–Prager–Yamakawa kernel along the centerline and apply a spectrally accurate “slender-body” quadrature to accurately resolve the hydrodynamics. Third, we propose a novel midpoint temporal integrator, which can correctly capture the Ito drift terms that arise in the overdamped Langevin equation. For two separate examples, we verify that the equilibrium distribution for the Chebyshev grid is a good approximation of the blob-link one and that our temporal integrator for overdamped Langevin dynamics samples the equilibrium GB distribution for sufficiently small time step sizes. We also study the dynamics of relaxation of an initially straight filament and find that as few as 12 Chebyshev nodes provide a good approximation to the dynamics while allowing a time step size two orders of magnitude larger than a resolved blob-link simulation. We conclude by applying our approach to a suspension of cross-linked semiflexible fibers (neglecting hydrodynamic interactions between fibers), where we study how semiflexible fluctuations affect bundling dynamics. We find that semiflexible filaments bundle faster than rigid filaments even when the persistence length is large, but show that semiflexible bending fluctuations only further accelerate agglomeration when the persistence length and fiber length are of the same order.
more »
« less
- PAR ID:
- 10437391
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 15
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)In cells, cytoskeletal filament networks are responsible for cell movement, growth, and division. Filaments in the cytoskeleton are driven and organized by crosslinking molecular motors. In reconstituted cytoskeletal systems, motor activity is responsible for far-from-equilibrium phenomena such as active stress, self-organized flow, and spontaneous nematic defect generation. How microscopic interactions between motors and filaments lead to larger-scale dynamics remains incompletely understood. To build from motor–filament interactions to predict bulk behavior of cytoskeletal systems, more computationally efficient techniques for modeling motor–filament interactions are needed. Here, we derive a coarse-graining hierarchy of explicit and continuum models for crosslinking motors that bind to and walk on filament pairs. We compare the steady-state motor distribution and motor-induced filament motion for the different models and analyze their computational cost. All three models agree well in the limit of fast motor binding kinetics. Evolving a truncated moment expansion of motor density speeds the computation by 103–106 compared to the explicit or continuous-density simulations, suggesting an approach for more efficient simulation of large networks. These tools facilitate further study of motor–filament networks on micrometer to millimeter length scales.more » « less
-
Merks, Roeland M.H. (Ed.)Cross-linked actin networks are the primary component of the cell cytoskeleton and have been the subject of numerous experimental and modeling studies. While these studies have demonstrated that the networks are viscoelastic materials, evolving from elastic solids on short timescales to viscous fluids on long ones, questions remain about the duration of each asymptotic regime, the role of the surrounding fluid, and the behavior of the networks on intermediate timescales. Here we perform detailed simulations of passively cross-linked non-Brownian actin networks to quantify the principal timescales involved in the elastoviscous behavior, study the role of nonlocal hydrodynamic interactions, and parameterize continuum models from discrete stochastic simulations. To do this, we extend our recent computational framework for semiflexible filament suspensions, which is based on nonlocal slender body theory, to actin networks with dynamic cross linkers and finite filament lifetime. We introduce a model where the cross linkers are elastic springs with sticky ends stochastically binding to and unbinding from the elastic filaments, which randomly turn over at a characteristic rate. We show that, depending on the parameters, the network evolves to a steady state morphology that is either an isotropic actin mesh or a mesh with embedded actin bundles. For different degrees of bundling, we numerically apply small-amplitude oscillatory shear deformation to extract three timescales from networks of hundreds of filaments and cross linkers. We analyze the dependence of these timescales, which range from the order of hundredths of a second to the actin turnover time of several seconds, on the dynamic nature of the links, solvent viscosity, and filament bending stiffness. We show that the network is mostly elastic on the short time scale, with the elasticity coming mainly from the cross links, and viscous on the long time scale, with the effective viscosity originating primarily from stretching and breaking of the cross links. We show that the influence of nonlocal hydrodynamic interactions depends on the network morphology: for homogeneous meshworks, nonlocal hydrodynamics gives only a small correction to the viscous behavior, but for bundled networks it both hinders the formation of bundles and significantly lowers the resistance to shear once bundles are formed. We use our results to construct three-timescale generalized Maxwell models of the networks.more » « less
-
We propose a novel integral model describing the motion of both flexible and rigid slender fibers in viscous flow and develop a numerical method for simulating dynamics of curved rigid fibers. The model is derived from nonlocal slender body theory (SBT), which approximates flow near the fiber using singular solutions of the Stokes equations integrated along the fiber centerline. In contrast to other models based on (singular) SBT, our model yields a smooth integral kernel which incorporates the (possibly varying) fiber radius naturally. The integral operator is provably negative definite in a nonphysical idealized geometry, as expected from the partial differential equation theory. This is numerically verified in physically relevant geometries. We discuss the convergence and stability of a numerical method for solving the integral equation. The accuracy of the model and method is verified against known models for ellipsoids. Finally, we develop an algorithm for computing dynamics of rigid fibers with complex geometries in the case where the fiber density is much greater than that of the fluid, for example, in turbulent gas-fiber suspensions.more » « less
-
null (Ed.)Many-body interactions in systems of active matter can cause particles to move collectively and self-organize into dynamic structures with long-range order. In cells, the self-assembly of cytoskeletal filaments is critical for cellular motility, structure, intracellular transport, and division. Semiflexible cytoskeletal filaments driven by polymerization or motor-protein interactions on a two-dimensional substrate, such as the cell cortex, can induce filament bending and curvature leading to interesting collective behavior. For example, the bacterial cell-division filament FtsZ is known to have intrinsic curvature that causes it to self-organize into rings and vortices, and recent experiments reconstituting the collective motion of microtubules driven by motor proteins on a surface have observed chiral symmetry breaking of the collective behavior due to motor-induced curvature of the filaments. Previous work on the self-organization of driven filament systems have not studied the effects of curvature and filament structure on collective behavior. In this work, we present Brownian dynamics simulation results of driven semiflexible filaments with intrinsic curvature and investigate how the interplay between filament rigidity and radius of curvature can tune the self-organization behavior in homochiral systems and heterochiral mixtures. We find a curvature-induced reorganization from polar flocks to self-sorted chiral clusters, which is modified by filament flexibility. This transition changes filament transport from ballistic to diffusive at long timescales.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0144242
