skip to main content


Title: Simulations of dynamically cross-linked actin networks: Morphology, rheology, and hydrodynamic interactions
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
Award ID(s):
2052515 1646339
NSF-PAR ID:
10338788
Author(s) / Creator(s):
; ; ;
Editor(s):
Merks, Roeland M.H.
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
17
Issue:
12
ISSN:
1553-7358
Page Range / eLocation ID:
e1009240
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Cells dynamically control their material properties through remodeling of the actin cytoskeleton, an assembly of cross-linked networks and bundles formed from the biopolymer actin. We recently found that cross-linked networks of actin filaments reconstituted in vitro can exhibit adaptive behavior and thus serve as a model system to understand the underlying mechanisms of mechanical adaptation of the cytoskeleton. In these networks, training, in the form of applied shear stress, can induce asymmetry in the nonlinear elasticity. Here, we explore control over this mechanical hysteresis by tuning the concentration and mechanical properties of cross-linking proteins in both experimental and simulated networks. We find that this effect depends on two conditions: the initial network must exhibit nonlinear strain stiffening, and filaments in the network must be able to reorient during training. Hysteresis depends strongly and non-monotonically on cross-linker concentration, with a peak at moderate concentrations. In contrast, at low concentrations, where the network does not strain stiffen, or at high concentrations, where filaments are less able to rearrange, there is little response to training. Additionally, we investigate the effect of changing cross-linker properties and find that longer or more flexible cross-linkers enhance hysteresis. Remarkably plotting hysteresis against alignment after training yields a single curve regardless of the physical properties or concentration of the cross-linkers. 
    more » « less
  2. The shape of most animal cells is controlled by the actin cortex, a thin network of dynamic actin filaments (F-actin) situated just beneath the plasma membrane. The cortex is held far from equilibrium by both active stresses and polymer turnover: Molecular motors drive deformations required for cell morphogenesis, while actin-filament disassembly dynamics relax stress and facilitate cortical remodeling. While many aspects of actin-cortex mechanics are well characterized, a mechanistic understanding of how nonequilibrium actin turnover contributes to stress relaxation is still lacking. To address this, we developed a reconstituted in vitro system of entangled F-actin, wherein the steady-state length and turnover rate of F-actin are controlled by the actin regulatory proteins cofilin, profilin, and formin, which sever, recycle, and assemble filaments, respectively. Cofilin-mediated severing accelerates the turnover and spatial reorganization of F-actin, without significant changes to filament length. We demonstrate that cofilin-mediated severing is a single-timescale mode of stress relaxation that tunes the low-frequency viscosity over two orders of magnitude. These findings serve as the foundation for understanding the mechanics of more physiological F-actin networks with turnover and inform an updated microscopic model of single-filament turnover. They also demonstrate that polymer activity, in the form of ATP hydrolysis on F-actin coupled to nucleotide-dependent cofilin binding, is sufficient to generate a form of active matter wherein asymmetric filament disassembly preserves filament number despite sustained severing. 
    more » « less
  3. 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
  4. “Viscosity is the most ubiquitous dissipative mechanical behavior” (Maugin, 1999). Despite its ubiquity, even for those systems where the mechanisms causing viscous and other forms of dissipation are known there are only a few quantitative models that extract the macroscopic rheological response from these microscopic mechanisms. One such mechanism is the stochastic breaking and forming of bonds which is present in polymer networks with transient cross-links, strong inter-layer bonding between graphene sheets, and sliding dry friction. In this paper we utilize a simple yet flexible model to show analytically how stochastic bonds can induce an array of rheological behaviors at the macroscale. We find that varying the bond interactions induces a Maxwell-type macroscopic material behavior with Newtonian viscosity, shear thinning, shear thickening, or solid like friction when subjected to shear at constant rates. When bond rupture is independent of the force applied, Newtonian viscosity is the predominant behavior. When bond breaking is accelerated by the applied force, a shear thinning response becomes most prevalent. Further connections of the macroscopic response to the interaction potential and rates of bonding and unbonding are illustrated through phase diagrams and analysis of limiting cases. Finally, we apply this model to polymer networks and to experimental data on “solid bridges” in polydisperse granular media. We imagine possible applications to material design through engineering bonds with specific interactions to bring about a desired macroscopic behavior. 
    more » « less
  5. null (Ed.)
    The rapid equilibrium fluctuations of water molecules are intimately connected to the rheological response; molecular motions resetting the local structure and stresses seen as flow and volume changes. In the case of water or hydrogen bonding liquids generally, the relationship is a non-trivial consideration due to strong directional interactions complicating theoretical models and necessitating clear observation of the timescale and nautre of the associated equilibrium motions. Recent work has illustrated a coincidence of timescales for short range sub-picosecond motions and the implied timescale for the shear viscosity response in liquid water. Here, neutron and light scattering methods are used to experimentally illustrate the timescale of bulk viscosity and provide a description of the associated molecular relaxation. Brillouin scattering has been used to establish the timescale of bulk viscosity; and borrowing the Maxwell approach, the ratio of the bulk viscosity, ζ , to the bulk modulus, K , yields a relaxation time, τ B , which emerges on the order of 1–2 ps in the 280 K to 303 K temperature range. Inelastic neutron scattering is subsequently used to describe the motions of water and heavy water at the molecular scale, providing both coherent and incoherent scattering data. A rotational (alternatively described as localized) motion of water protons on the 1–2 ps timescale is apparent in the incoherent scattering spectra of water, while the coherent spectra from D 2 O on the length scale of the first sharp diffraction peak, describing the microscopic density fluctuations of water, confirms the relaxation of water structure at a comparable timescale of 1–2 ps. The coincidence of these three timescales provides a mechanistic description of the bulk viscous response, with the local structure resetting due to rotational/localized motions on the order of 1–2 ps, approximately three times slower than the relaxations associated with shear viscosity. In this way we show that the shear viscous response is most closely associated with changes in water network connectivity, while the bulk viscous response is associated with local density fluctuations. 
    more » « less