Rheology of semiflexible bundle networks with transient linkers

Non-Maxwellian viscoelastic stress relaxations in soft matter

Viscoelastic stress relaxation is a basic characteristic of soft matter systems such as colloids, gels, and biological networks. Although the Maxwell model of linear viscoelasticity provides a classical description of stress relaxation, it is often not sufficient for capturing the complex relaxation dynamics of soft matter. In this Tutorial, we introduce and discuss the physics of non-Maxwellian linear stress relaxation as observed in soft materials, the ascribed origins of this effect in different systems, and appropriate models that can be used to capture this relaxation behavior. We provide a basic toolkit that can assist the understanding and modeling of the mechanical relaxation of soft materials for diverse applications.

Organization and dynamics of cross-linked actin filaments in confined environments

The organization of the actin cytoskeleton is impacted by the interplay between physical confinement, features of cross-linking proteins, and deformations of semiflexible actin filaments. Some cross-linking proteins preferentially bind filaments in parallel, although others bind more indiscriminately. However, a quantitative understanding of how the mode of binding influences the assembly of actin networks in confined environments is lacking. Here we employ coarse-grained computer simulations to study the dynamics and organization of semiflexible actin filaments in confined regions upon the addition of cross-linkers. We characterize how the emergent behavior is influenced by the system shape, the number and type of cross-linking proteins, and the length of filaments. Structures include isolated clusters of filaments, highly connected filament bundles, and networks of interconnected bundles and loops. Elongation of one dimension of the system promotes the formation of long bundles that align with the elongated axis. Dynamics are governed by rapid cross-linking into aggregates, followed by a slower change in their shape and connectivity. Cross-linking decreases the average bending energy of short or sparsely connected filaments by suppressing shape fluctuations. However, it increases the average bending energy in highly connected networks because filament bundles become deformed, and small numbers of filaments exhibit long-lived, highly unfavorable configurations. Indiscriminate cross-linking promotes the formation of high-energy configurations due to the increased likelihood of unfavorable, difficult-to-relax configurations at early times. Taken together, this work demonstrates physical mechanisms by which cross-linker binding and physical confinement impact the emergent behavior of actin networks, which is relevant both in cells and in synthetic environments.

Structural Features and Nonlinear Rheology of Self-Assembled Networks of Cross-Linked Semiflexible Polymers

Disordered networks of semiflexible filaments are common support structures in biology. Familiar examples include fibrous matrices in blood clots, bacterial biofilms, and essential components of cells and tissues of plants, animals, and fungi. Despite the ubiquity of these networks in biomaterials, we have only a limited understanding of the relationship between their structural features and their highly strain-sensitive mechanical properties. In this work, we perform simulations of three-dimensional networks produced by the irreversible formation of cross-links between linker-decorated semiflexible filaments. We characterize the structure of networks formed by a simple diffusion-dependent assembly process and measure their associated steady-state rheological features at finite temperature over a range of applied prestrains that encompass the strain-stiffening transition. We quantify the dependence of network connectivity on cross-linker availability and detail the associated connectivity dependence of both linear elasticity and nonlinear strain-stiffening behavior, drawing comparisons with prior experimental measurements of the cross-linker concentration-dependent elasticity of actin gels.

Interplay between Brownian motion and cross-linking kinetics controls bundling dynamics in actin networks

Morphology changes in cross-linked actin networks are important in cell motility, division, and cargo transport. Here, we study the transition from a weakly cross-linked network of actin filaments to a heavily cross-linked network of actin bundles through microscopic Brownian dynamics simulations. We show that this transition occurs in two stages: first, a composite bundle network of small and highly aligned bundles evolves from cross-linking of individual filaments and, second, small bundles coalesce into the clustered bundle state. We demonstrate that Brownian motion speeds up the first stage of this process at a faster rate than the second. We quantify the time to reach the composite bundle state and show that it strongly increases as the mesh size increases only when the concentration of cross-links is small and that it remains roughly constant if we decrease the relative ratio of cross-linkers as we increase the actin concentration. Finally, we examine the dependence of the bundling timescale on filament length, finding that shorter filaments bundle faster because they diffuse faster.

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.

Directed force propagation in semiflexible networks

We consider the propagation of tension along specific filaments of a semiflexible filament network in response to the application of a point force using a combination of numerical simulations and analytic theory. We find the distribution of force within the network is highly heterogeneous, with a small number of fibers supporting a significant fraction of the applied load over distances of multiple mesh sizes surrounding the point of force application. We suggest that these structures may be thought of as tensile force chains, whose structure we explore via simulation. We develop self-consistent calculations of the point-force response function and introduce a transfer matrix approach to explore the decay of tension (into bending) energy and the branching of tensile force chains in the network.

Viscoelasticity of 3D actin networks dictated by the mechanochemical characteristics of cross-linkers

In this study, we report a computational investigation on how the mechanochemical characteristics of crosslinking molecules influence the viscoelasticity of three dimensional F-actin networks, an issue of key interest in analyzing the behavior of living cells and biological gels. In particular, it was found that the continuous breakage and rebinding of cross-linkers result in a locally peaked loss modulus in the rheology spectrum of the network, reflecting the fact that maximum energy dissipation is achieved when the driving frequency of the applied oscillating shear becomes comparable to the dissociation/association rate of crosslinking molecules. In addition, we showed that when subjected to constant rate of shear, an actin network can exhibit either strain hardening or softening depending on the ratio between the loading rate and unbinding speed of cross-linkers. A criterion for predicting the transition from softening to hardening was also obtained, in agreement with recent experiments. Finally, significant structural evolution was found to occur in random networks undergoing mechanical "training" (i.e. under a constant applied shear stress over a period of time), eventually leading to a pronounced anisotropic response of the network afterward which again is consistent with experimental observations.

Nonlinear stress relaxation of transiently crosslinked biopolymer networks

A long-standing puzzle in the rheology of living cells is the origin of the experimentally observed long-time stress relaxation. The mechanics of the cell is largely dictated by the cytoskeleton, which is a biopolymer network consisting of transient crosslinkers, allowing for stress relaxation over time. Moreover, these networks are internally stressed due to the presence of molecular motors. In this work we propose a theoretical model that uses a mode-dependent mobility to describe the stress relaxation of such prestressed transient networks. Our theoretical predictions agree favorably with experimental data of reconstituted cytoskeletal networks and may provide an explanation for the slow stress relaxation observed in cells.

Braiding Dynamics in Semiflexible Filament Bundles under Oscillatory Forcing

We examine the nonequilibrium production of topological defects-braids-in semiflexible filament bundles under cycles of compression and tension. During these cycles, the period of compression facilitates the thermally activated pair production of braid/anti-braid pairs, which then may separate when the bundle is under tension. As a result, appropriately tuned alternating periods of compression and extension should lead to the proliferation of braid defects in a bundle so that the linear density of these pairs far exceeds that expected in the thermal equilibrium. Secondly, we examine the slow extension of braided bundles under tension, showing that their end-to-end length creeps nonmonotonically under a fixed force due to braid deformation and the motion of the braid pair along the bundle. We conclude with a few speculations regarding experiments on semiflexible filament bundles and their networks.

A computational framework for modeling cell-matrix interactions in soft biological tissues

Living soft tissues appear to promote the development and maintenance of a preferred mechanical state within a defined tolerance around a so-called set point. This phenomenon is often referred to as mechanical homeostasis. In contradiction to the prominent role of mechanical homeostasis in various (patho)physiological processes, its underlying micromechanical mechanisms acting on the level of individual cells and fibers remain poorly understood, especially how these mechanisms on the microscale lead to what we macroscopically call mechanical homeostasis. Here, we present a novel computational framework based on the finite element method that is constructed bottom up, that is, it models key mechanobiological mechanisms such as actin cytoskeleton contraction and molecular clutch behavior of individual cells interacting with a reconstructed three-dimensional extracellular fiber matrix. The framework reproduces many experimental observations regarding mechanical homeostasis on short time scales (hours), in which the deposition and degradation of extracellular matrix can largely be neglected. This model can serve as a systematic tool for future in silico studies of the origin of the numerous still unexplained experimental observations about mechanical homeostasis.

A novel modelling and simulation approach for the hindered mobility of charged particles in biological hydrogels

This article presents a novel computational model to study the selective filtering of biological hydrogels due to the surface charge and size of diffusing particles. It is the first model that includes the random three-dimensional fibre orientation and connectivity of the biopolymer network and that accounts for elastic deformations of the fibres by means of beam theory. As a key component of the model, novel formulations are proposed both for the electrostatic and repulsive steric interactions between a spherical particle and a beam. In addition to providing a thorough validation of the model, the presented computational studies yield new insights into the underlying mechanisms of hindered particle mobility, especially regarding the influence of the aforementioned aspects that are unique to this model. It is found that the precise distribution of fibre and thus charge agglomerations in the network have a crucial influence on the mobility of oppositely charged particles and gives rise to distinct motion patterns. Considering the high practical significance for instance with respect to targeted drug release or infection defence, the provided proof of concept motivates further advances of the model towards a truly predictive computational tool that allows a case- and patient-specific assessment for real (biological) systems.

Topological defects produce kinks in biopolymer filament bundles

Bundles of stiff filaments are ubiquitous in the living world, found both in the cytoskeleton and in the extracellular medium. These bundles are typically held together by smaller cross-linking molecules. We demonstrate, analytically, numerically, and experimentally, that such bundles can be kinked, that is, have localized regions of high curvature that are long-lived metastable states. We propose three possible mechanisms of kink stabilization: a difference in trapped length of the filament segments between two cross-links, a dislocation where the endpoint of a filament occurs within the bundle, and the braiding of the filaments in the bundle. At a high concentration of cross-links, the last two effects lead to the topologically protected kinked states. Finally, we explore, numerically and analytically, the transition of the metastable kinked state to the stable straight bundle.

Modeling of Cell-Mediated Self-Assembled Colloidal Scaffolds

A critical step in tissue engineering is the design and synthesis of 3D biocompatible matrices (scaffolds) to support and guide the proliferation of cells and tissue growth. The most existing techniques rely on the processing of scaffolds under controlled conditions and then implanting them in vivo, with questions related to biocompatibility and implantation that are still challenging. As an alternative, it was proposed to assemble the scaffolds in loco through the self-organization of colloidal particles mediated by cells. To overcome the difficulty to test experimentally all the relevant parameters, we propose the use of large-scale numerical simulation as a tool to reach useful predictive information and to interpret experimental results. Thus, in this study, we combine experiments, particle-based simulations, and mean-field calculations to show that, in general, the size of the self-assembled scaffold scales with the cell-to-particle ratio. However, we have found an optimal value of this ratio, for which the size of the scaffold is maximal when the cell-cell adhesion is suppressed. These results suggest that the size and structure of the self-assembled scaffolds may be designed by tuning the adhesion between cells in the colloidal suspension.

Smoluchowski equations for linker-mediated irreversible aggregation

We developed a generalized Smoluchowski framework to study linker-mediated aggregation, where linkers and particles are explicitly taken into account. We assume that the bonds between linkers and particles are irreversible, and that clustering occurs through limited diffusion aggregation. The kernel is chosen by analogy with single-component diffusive aggregation but the clusters are distinguished by their number of particles and linkers. We found that the dynamics depends on three relevant factors, all tunable experimentally: (i) the ratio of the diffusion coefficients of particles and linkers; (ii) the relative number of particles and linkers; and (iii) the maximum number of linkers that may bond to a single particle. To solve the Smoluchoski equations analytically we employ a scaling hypothesis that renders the fraction of bondable sites of a cluster independent of the size of the cluster, at each instant. We perform numerical simulations of the corresponding lattice model to test this hypothesis. We obtain results for the asymptotic limit, and the time evolution of the bonding probabilities and the size distribution of the clusters. These findings are in agreement with experimental results reported in the literature and shed light on unexplained experimental observations.

On the mechanical response of the actomyosin cortex during cell indentations

A mechanical model is presented to analyze the mechanics and dynamics of the cell cortex during indentation. We investigate the impact of active contraction on the cross-linked actin network for different probe sizes and indentation rates. The essential molecular mechanisms of filament stretching, cross-linking and motor activity, are represented by an active and viscous mechanical continuum. The filaments behave as worm-like chains linked either by passive rigid linkers or by myosin motors. In the first example, the effects of probe size and loading rate are evaluated using the model for an idealized rounded cell shape in which properties are based on the results of parallel-plate rheometry available in the literature. Extreme cases of probe size and indentation rate are taken into account. Afterward, AFM experiments were done by engaging smooth muscle cells with both sharp and spherical probes. By inverse analysis with finite element software, our simulations mimicking the experimental conditions show the model is capable of fitting the AFM data. The results provide spatiotemporal dependence on the size and rate of the mechanical stimuli. The model captures the general features of the cell response. It characterizes the actomyosin cortex as an active solid at short timescales and as a fluid at longer timescales by showing (1) higher levels of contraction in the zones of high curvature; (2) larger indentation forces as the probe size increases; and (3) increase in the apparent modulus with the indentation depth but no dependence on the rate of the mechanical stimuli. The methodology presented in this work can be used to address and predict microstructural dependence on the force generation of living cells, which can contribute to understanding the broad spectrum of results in cell experiments.

Varying crosslinking motifs drive the mesoscale mechanics of actin-microtubule composites

The cytoskeleton precisely tunes its mechanics by altering interactions between semiflexible actin filaments, rigid microtubules, and crosslinking proteins. We use optical tweezers microrheology and confocal microscopy to characterize how varying crosslinking motifs impact the mesoscale mechanics and mobility of actin-microtubule composites. We show that, upon subtle changes in crosslinking patterns, composites can exhibit two distinct classes of force response - primarily elastic versus more viscous. For example, a composite in which actin and microtubules are crosslinked to each other but not to themselves is markedly more elastic than one in which both filaments are independently crosslinked. Notably, this distinction only emerges at mesoscopic scales in response to nonlinear forcing, whereas varying crosslinking motifs have little impact on the microscale mechanics and mobility. Our unexpected scale-dependent results not only inform the physics underlying key cytoskeleton processes and structures, but, more generally, provide valuable perspective to materials engineering endeavors focused on polymer composites.

Crosslinking actin networks produces compressive force

Actin has been shown to be essential for clathrin-mediated endocytosis in yeast. However, actin polymerization alone is likely insufficient to produce enough force to deform the membrane against the huge turgor pressure of yeast cells. In this paper, we used Brownian dynamics simulations to demonstrate that crosslinking of a meshwork of nonpolymerizing actin filaments is able to produce compressive forces. We show that the force can be up to several thousand pico-Newtons if the crosslinker has a high stiffness. The force decays over time as a result of crosslinker turnover, and is a result of converting chemical binding energy into elastic energy.

Origin of Slow Stress Relaxation in the Cytoskeleton

Dynamically cross-linked semiflexible biopolymers such as the actin cytoskeleton govern the mechanical behavior of living cells. Semiflexible biopolymers nonlinearly stiffen in response to mechanical loads, whereas the cross-linker dynamics allow for stress relaxation over time. Here we show, through rheology and theoretical modeling, that the combined nonlinearity in time and stress leads to an unexpectedly slow stress relaxation, similar to the dynamics of disordered systems close to the glass transition. Our work suggests that transient cross-linking combined with internal stress can explain prior reports of soft glassy rheology of cells, in which the shear modulus increases weakly with frequency.

Optimal number of linkers per monomer in linker-mediated aggregation

We study the dynamics of diffusion-limited irreversible aggregation of monomers, where bonds are mediated by linkers. We combine kinetic Monte Carlo simulations of a lattice model with a mean-field theory to study the dynamics when the diffusion of aggregates is negligible and only monomers diffuse. We find two values of the number of linkers per monomer which maximize the size of the largest aggregate. We explain the existence of the two maxima based on the distribution of linkers per monomer. This observation is well described by a simple mean-field model. We also show that a relevant parameter is the ratio of the diffusion coefficients of monomers and linkers. In particular, when this ratio is close to ten, the two maxima merge at a single maximum.

Frustrated binding of biopolymer crosslinkers

Transiently crosslinked actin filament networks allow cells to combine elastic rigidity with the ability to deform viscoelastically. Theoretical models of semiflexible polymer networks predict that the crosslinker unbinding rate governs the timescale beyond which viscoelastic flow occurs. However a direct comparison between network and crosslinker dynamics is lacking. Here we measure the network's stress relaxation timescale using rheology and the lifetime of bound crosslinkers using fluorescence recovery after photobleaching (FRAP). Intriguingly, we observe that the crosslinker unbinding rate measured by FRAP is more than an order of magnitude slower than the rate measured by rheology. We rationalize this difference with a three-state model where crosslinkers are bound to either 0, 1 or 2 filaments, which allows us to extract crosslinker transition rates that are otherwise difficult to access. We find that the unbinding rate of singly bound crosslinkers is nearly two orders of magnitude slower than for doubly bound ones. We attribute the increased unbinding rate of doubly bound crosslinkers to the high stiffness of biopolymers, which frustrates crosslinker binding.

Linear rheology of reversibly cross-linked biopolymer networks

We suggest a simple model for reversible cross-links, binding, and unbinding to/from a network of semiflexible polymers. The resulting frequency dependent response of the network to an applied shear is calculated via Brownian dynamics simulations. It is shown to be rather complex with the time scale of the linkers competing with the excitations of the network. If the lifetime of the linkers is the longest time scale, as is indeed the case in most biological networks, then a distinct low frequency peak of the loss modulus develops. The storage modulus shows a corresponding decay from its plateau value, which for irreversible cross-linkers extends all the way to the static limit. This additional relaxation mechanism can be controlled by the relative weight of reversible and irreversible linkers.

Mechanical hysteresis in actin networks

Understanding the response of complex materials to external force is central to fields ranging from materials science to biology. Here, we describe a novel type of mechanical adaptation in cross-linked networks of F-actin, a ubiquitous protein found in eukaryotic cells. We show that shear stress changes the network's nonlinear mechanical response even long after that stress is removed. The duration, magnitude and direction of forcing history all change this mechanical response. While the mechanical hysteresis is long-lived, it can be simply erased by force application in the opposite direction. We further show that the observed mechanical adaptation is consistent with stress-dependent changes in the nematic order of the constituent filaments. Thus, this mechanical hysteresis arises from the changes in non-linear response that originates from stress-induced changes to filament orientation. This demonstrates that F-actin networks can exhibit analog read-write mechanical hysteretic properties, which can be used for adaptation to mechanical stimuli.

Filament turnover tunes both force generation and dissipation to control long-range flows in a model actomyosin cortex

Actomyosin-based cortical flow is a fundamental engine for cellular morphogenesis. Cortical flows are generated by cross-linked networks of actin filaments and myosin motors, in which active stress produced by motor activity is opposed by passive resistance to network deformation. Continuous flow requires local remodeling through crosslink unbinding and and/or filament disassembly. But how local remodeling tunes stress production and dissipation, and how this in turn shapes long range flow, remains poorly understood. Here, we study a computational model for a cross-linked network with active motors based on minimal requirements for production and dissipation of contractile stress: Asymmetric filament compliance, spatial heterogeneity of motor activity, reversible cross-links and filament turnover. We characterize how the production and dissipation of network stress depend, individually, on cross-link dynamics and filament turnover, and how these dependencies combine to determine overall rates of cortical flow. Our analysis predicts that filament turnover is required to maintain active stress against external resistance and steady state flow in response to external stress. Steady state stress increases with filament lifetime up to a characteristic time τm, then decreases with lifetime above τm. Effective viscosity increases with filament lifetime up to a characteristic time τc, and then becomes independent of filament lifetime and sharply dependent on crosslink dynamics. These individual dependencies of active stress and effective viscosity define multiple regimes of steady state flow. In particular our model predicts that when filament lifetimes are shorter than both τc and τm, the dependencies of effective viscosity and steady state stress on filament turnover cancel one another, such that flow speed is insensitive to filament turnover, and shows a simple dependence on motor activity and crosslink dynamics. These results provide a framework for understanding how animal cells tune cortical flow through local control of network remodeling.

Force percolation of contractile active gels

Living systems provide a paradigmatic example of active soft matter. Cells and tissues comprise viscoelastic materials that exert forces and can actively change shape. This strikingly autonomous behavior is powered by the cytoskeleton, an active gel of semiflexible filaments, crosslinks, and molecular motors inside cells. Although individual motors are only a few nm in size and exert minute forces of a few pN, cells spatially integrate the activity of an ensemble of motors to produce larger contractile forces (∼nN and greater) on cellular, tissue, and organismal length scales. Here we review experimental and theoretical studies on contractile active gels composed of actin filaments and myosin motors. Unlike other active soft matter systems, which tend to form ordered patterns, actin-myosin systems exhibit a generic tendency to contract. Experimental studies of reconstituted actin-myosin model systems have long suggested that a mechanical interplay between motor activity and the network's connectivity governs this contractile behavior. Recent theoretical models indicate that this interplay can be understood in terms of percolation models, extended to include effects of motor activity on the network connectivity. Based on concepts from percolation theory, we propose a state diagram that unites a large body of experimental observations. This framework provides valuable insights into the mechanisms that drive cellular shape changes and also provides design principles for synthetic active materials.

Intact Telopeptides Enhance Interactions between Collagens

Collagen is the fundamental structural component of a wide range of connective tissues and of the extracellular matrix. It undergoes self-assembly from individual triple-helical proteins into well-ordered fibrils, a process that is key to tissue development and homeostasis, and to processes such as wound healing. Nucleation of this assembly is known to be slowed considerably by pepsin removal of short nonhelical regions that flank collagen's triple helix, known as telopeptides. Using optical tweezers to perform microrheology measurements, we explored the changes in viscoelasticity of solutions of collagen with and without intact telopeptides. Our experiments reveal that intact telopeptides contribute a significant frequency-dependent enhancement of the complex shear modulus. An analytical model of polymers associating to establish chemical equilibrium among higher-order species shows trends in G' and G″ consistent with our experimental observations, including a concentration-dependent crossover in G″/c around 300 Hz. This work suggests that telopeptides facilitate transient intermolecular interactions between collagen proteins, even in the acidic conditions used here.

Multiscale model predicts increasing focal adhesion size with decreasing stiffness in fibrous matrices

We describe a multiscale model that incorporates force-dependent mechanical plasticity induced by interfiber cross-link breakage and stiffness-dependent cellular contractility to predict focal adhesion (FA) growth and mechanosensing in fibrous extracellular matrices (ECMs). The model predicts that FA size depends on both the stiffness of ECM and the density of ligands available to form adhesions. Although these two quantities are independent in commonly used hydrogels, contractile cells break cross-links in soft fibrous matrices leading to recruitment of fibers, which increases the ligand density in the vicinity of cells. Consequently, although the size of focal adhesions increases with ECM stiffness in nonfibrous and elastic hydrogels, plasticity of fibrous networks leads to a departure from the well-described positive correlation between stiffness and FA size. We predict a phase diagram that describes nonmonotonic behavior of FA in the space spanned by ECM stiffness and recruitment index, which describes the ability of cells to break cross-links and recruit fibers. The predicted decrease in FA size with increasing ECM stiffness is in excellent agreement with recent observations of cell spreading on electrospun fiber networks with tunable cross-link strengths and mechanics. Our model provides a framework to analyze cell mechanosensing in nonlinear and inelastic ECMs.

The dynamics of filament assembly define cytoskeletal network morphology

The actin cytoskeleton is a key component in the machinery of eukaryotic cells, and it self-assembles out of equilibrium into a wide variety of biologically crucial structures. Although the molecular mechanisms involved are well characterized, the physical principles governing the spatial arrangement of actin filaments are not understood. Here we propose that the dynamics of actin network assembly from growing filaments results from a competition between diffusion, bundling and steric hindrance, and is responsible for the range of observed morphologies. Our model and simulations thus predict an abrupt dynamical transition between homogeneous and strongly bundled networks as a function of the actin polymerization rate. This suggests that cells may effect dramatic changes to their internal architecture through minute modifications of their nonequilibrium dynamics. Our results are consistent with available experimental data.

The dynamics of actin cytoskeleton is essential to the function of living cells. Here, Foffano et al. describe a nonequilibrium filament model to mimic the formation of cytoskeleton and pinpoint the key role played by the actin entanglement during the transition from homogeneous to bundled networks.

Discontinuous bundling transition in semiflexible polymer networks induced by Casimir interactions

Fluctuation-induced interactions are an important organizing principle in a variety of soft matter systems. We investigate the role of fluctuation-based or thermal Casimir interactions between cross linkers in a semiflexible network. One finds that, by integrating out the polymer degrees of freedom, there is an attractive logarithmic potential between nearest-neighbor cross linkers in a bundle, with a significantly weaker next-nearest-neighbor interaction. Here we show that a one-dimensional gas of these strongly interacting linkers in equilibrium with a source of unbound ones admits a discontinuous phase transition between a sparsely and a densely bound bundle. This discontinuous transition induced by the long-ranged nature of the Casimir interaction allows for a similarly abrupt structural transition in semiflexible filament networks between a low cross linker density isotropic phase and a higher cross link density bundle network. We support these calculations with the results of finite element Brownian dynamics simulations of semiflexible filaments and transient cross linkers.

Semiflexible Biopolymers in Bundled Arrangements

Bundles and networks of semiflexible biopolymers are key elements in cells, lending them mechanical integrity while also enabling dynamic functions. Networks have been the subject of many studies, revealing a variety of fundamental characteristics often determined via bulk measurements. Although bundles are equally important in biological systems, they have garnered much less scientific attention since they have to be probed on the mesoscopic scale. Here, we review theoretical as well as experimental approaches, which mainly employ the naturally occurring biopolymer actin, to highlight the principles behind these structures on the single bundle level.

Viscoelasticity of reversibly crosslinked networks of semiflexible polymers

We present a theoretical framework for the linear and nonlinear viscoelastic properties of reversibly crosslinked networks of semiflexible polymers. In contrast to affine models where network strain couples to the polymer end-to-end distance, in our model strain rather serves to locally distort the network structure. This induces bending modes in the polymer filaments, the properties of which are slaved to the surrounding network structure. Specifically, we investigate the frequency-dependent linear rheology, in particular in combination with crosslink binding-unbinding processes. We also develop schematic extensions to describe the nonlinear response during creep measurements as well as during constant strain-rate ramps.

Response of biopolymer networks governed by the physical properties of cross-linking molecules

In this study, we examine how the physical properties of cross-linking molecules affect the bulk response of bio-filament networks, an outstanding question in the study of biological gels and the cytoskeleton. We show that the stress-strain relationship of such networks typically undergoes linear increase - strain hardening - stress serration - total fracture transitions due to the interplay between the bending and stretching of individual filaments and the deformation and breakage of cross-linkers. Interestingly, the apparent network modulus is found to scale with the linear and rotational stiffness of the crosslinks to a power exponent of 0.78 and 0.13, respectively. In addition, the network fracture energy will reach its minimum at intermediate rotational compliance values, reflecting the fact that most of the strain energy will be stored in the distorted filaments with rigid cross-linkers while the imposed deformation will be "evenly" distributed among significantly more crosslinking molecules with high rotational compliance.

Weak reversible cross links may decrease the strength of aligned fiber bundles

Reversible cross-linking is an effective strategy to specifically tailor the mechanical properties of polymeric materials that can be found in a variety of biological as well as man-made materials. Using a simple model in this paper the influence of weak, reversible cross-links on the mechanical properties of aligned fiber bundles is investigated. Special emphasis in this analysis is put on the strength of the investigated structures. Using Monte Carlo methods two topologies of cross-links exceeding the strength of the covalent backbone are studied. Most surprisingly only two cross-links are sufficient to break the backbone of a multi chain system, resulting in a reduced strength of the material. The found effect crucially depends on the ratio of inter- to intra-chain cross-links and, thus, on the grafting density that determines this ratio.

Beam finite-element model of a molecular motor for the simulation of active fibre networks

Molecular motors are proteins that excessively increase the efficiency of subcellular transport processes. They allow for cell division, nutrient transport and even macroscopic muscle movement. In order to understand the effect of motors in large biopolymer networks, e.g. the cytoskeleton, we require a suitable model of a molecular motor. In this contribution, we present such a model based on a geometrically exact beam finite-element formulation. We discuss the numerical model of a non-processive motor such as myosin II, which interacts with actin filaments. Based on experimental data and inspired by the theoretical understanding offered by the power-stroke model and the swinging-cross-bridge model, we parametrize our numerical model in order to achieve the effect that a physiological motor has on its cargo. To this end, we introduce the mechanical and mathematical foundations of the model, then discuss its calibration, prove its usefulness by conducting finite-element simulations of actin-myosin motility assays and assess the influence of motors on the rheology of semi-flexible biopolymer networks.

Microstructure and rheology of microfibril-polymer networks

By using an adsorbing polymer in combination with mechanical de-agglomeration, the microstructure and rheological properties of networks of microfibrils could be controlled. By the addition of sodium carboxymethyl cellulose during de-agglomeration of networks of bacterial cellulose, the microstructure could be changed from an inhomogeneous network with bundles of microfibrils and voids to a more homogeneous spread and alignment of the particles. As a result the macroscopic rheological properties were altered. Although still elastic and gel-like in nature, the elasticity and viscous behavior of the network as a function of microfibril concentration is altered. The microstructure is thus changed by changing the surface properties of the building blocks leading to a direct influence on the materials macroscopic behavior.

Computational analysis of morphologies and phase transitions of cross-linked, semi-flexible polymer networks

The eukaryotic cytoskeleton is a protein fibre network mainly consisting of the semi-flexible biopolymer F-actin, microtubules and intermediate filaments. It is well known to exhibit a pronounced structural polymorphism, which enables intracellular processes such as cell adhesion, cell motility and cell division. We present a computational study on cross-linked networks of semi-flexible polymers, which offers a detailed analysis of the network structure and phase transitions from one morphology to another. We elaborate the morphological differences, their mechanical implications and the order of the observed phase transitions. Finally, we present a perspective on how the information gained in our simulations can be exploited in order to build both flexible and accurate, microstructurally informed, homogenized constitutive models of the cytoskeleton.

Bond breaking dynamics in semiflexible networks under load

We examine the bond-breaking dynamics of transiently cross-linked semiflexible networks using a single filament model in which that filament is peeled from an array of cross-linkers. We examine the effect of quenched disorder in the placement of the linkers along the filament and the effect of stochastic bond-breaking (assuming Bell model unbinding kinetics) on the dynamics of filament cross-linker dissociation and the statistics of ripping events. We find that bond forces decay exponentially away from the point of loading and that bond breaking proceeds sequentially down the linker array from the point of loading in a series of stochastic ripping events. We compare these theoretical predictions to the observed trajectories of large beads in a cross-linked microtubule network and identify the observed jumps of the bead with the linker rupture events predicted by the single filament model.

A single charge in the actin binding domain of fascin can independently tune the linear and non-linear response of an actin bundle network

Actin binding proteins (ABPs) not only set the structure of actin filament assemblies but also mediate the frequency-dependent viscoelastic moduli of cross-linked and bundled actin networks. Point mutations in the actin binding domain of those ABPs can tune the association and dissociation dynamics of the actin/ABP bond and thus modulate the network mechanics both in the linear and non-linear response regime. We here demonstrate how the exchange of a single charged amino acid in the actin binding domain of the ABP fascin triggers such a modulation of the network rheology. Whereas the overall structure of the bundle networks is conserved, the transition point from strain-hardening to strain-weakening sensitively depends on the cross-linker off-rate and the applied shear rate. Our experimental results are consistent both with numerical simulations of a cross-linked bundle network and a theoretical description of the bundle network mechanics which is based on non-affine bending deformations and force-dependent cross-link dynamics.