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
Author(s) / Creator(s):
; ; ;
Merks, Roeland M.H.
Date Published:
Journal Name:
PLOS Computational Biology
Page Range / eLocation ID:
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. Abstract

    Actin, an important component of eukaryotic cell cytoskeleton, regulates cell shape and transport. The morphology and biochemical properties of actin filaments are determined by their structure and protein-protein contacts. Crowded environments can organize filaments into bundles, but less is known about how they affect F-actin structure. This study used 2D IR spectroscopy and spectral calculations to examine how crowding and bundling impact the secondary structure and local environments in filaments and weakly or strongly bundled networks. The results reveal that bundling induces changes in actin’s secondary structure, leading to a decrease inβ-sheet and an increase in loop conformations. Strongly bundled networks exhibit a decrease in backbone solvent exposure, with less perturbedα-helices and nearly “locked”β-sheets. Similarly, the loops become less hydrated but maintain a dynamic environment. These findings highlight the role of loop structure in actin network morphology and stability under morphology control by PEG.

    more » « less
  5. “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