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 more » 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. « less
; ; ;
Merks, Roeland M.H.
Award ID(s):
2052515 1646339
Publication Date:
Journal Name:
PLOS Computational Biology
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 amore »single curve regardless of the physical properties or concentration of the cross-linkers.« less
  2. “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 tomore »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.« less
  3. 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 (alternativelymore »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.« less
  4. 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-actinmore »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.« less
  5. Depending on the physical and biochemical properties of actin-binding proteins, actin networks form different types of membrane protrusions at the cell periphery. Actin crosslinkers, which facilitate the interaction of actin filaments with one another, are pivotal in determining the mechanical properties and protrusive behavior of actin networks. Short crosslinkers such as fascin bundle F-actin to form rigid spiky filopodial protrusions. By encapsulation of fascin and actin in giant unilamellar vesicles (GUVs), we show that fascin-actin bundles cause various GUV shape changes by forming bundle networks or straight single bundles depending on GUV size and fascin concentration. We also show that the presence of a long crosslinker, α-actinin, impacts fascin-induced GUV shape changes and significantly impairs the formation of filopodia-like protrusions. Actin bundle-induced GUV shape changes are confirmed by light-induced disassembly of actin bundles leading to the reversal of GUV shape. Our study contributes to advancing the design of shape-changing minimal cells for better characterization of the interaction between lipid bilayer membranes and actin cytoskeleton.