Stress-stabilized subisostatic fiber networks in a ropelike limit

Unexpected softening of a fibrous matrix by contracting inclusions

Cells respond to the stiffness of their surrounding environment, but quantifying the stiffness of a fibrous matrix at the scale of a cell is complicated, due to the effects of nonlinearity and complex force transmission pathways resulting from randomness in fiber density and connections. While it is known that forces produced by individual contractile cells can stiffen the matrix, it remains unclear how simultaneous contraction of multiple cells in a fibrous matrix alters the stiffness at the scale of a cell. Here, we used computational modeling and experiments to quantify the stiffness of a random fibrous matrix embedded with multiple contracting inclusions, which mimicked the contractile forces of a cell. The results showed that when the matrix was free to contract as a result of the forces produced by the inclusions, the matrix softened rather than stiffened, which was surprising given that the contracting inclusions applied tensile forces to the matrix. Using the computational model, we identified that the underlying cause of the softening was that the majority of the fibers were under a local state of axial compression, causing buckling. We verified that this buckling-induced matrix softening was sufficient for cells to sense and respond by altering their morphology and force generation. Our findings reveal that the localized forces induced by cells do not always stiffen the matrix; rather, softening can occur in instances wherein the matrix can contract in response to the cell-generated forces. This study opens up new possibilities to investigate whether cell-induced softening contributes to maintenance of homeostatic conditions or progression of disease. STATEMENT OF SIGNIFICANCE: Mechanical interactions between cells and the surrounding matrix strongly influence cellular functions. Cell-induced forces can alter matrix properties, and much prior literature in this area focused on the influence of individual contracting cells. Cells in tissues are rarely solitary; rather, they are interspersed with neighboring cells throughout the matrix. As a result, the mechanics are complicated, leaving it unclear how the multiple contracting cells affect matrix stiffness. Here, we show that multiple contracting inclusions within a fibrous matrix can cause softening that in turn affects cell sensing and response. Our findings provide new directions to determine impacts of cell-induced softening on maintenance of tissue or progression of disease.

Effective medium theory for mechanical phase transitions of fiber networks

Networks of stiff fibers govern the elasticity of biological structures such as the extracellular matrix of collagen. These networks are known to stiffen nonlinearly under shear or extensional strain. Recently, it has been shown that such stiffening is governed by a strain-controlled athermal but critical phase transition, from a floppy phase below the critical strain to a rigid phase above the critical strain. While this phase transition has been extensively studied numerically and experimentally, a complete analytical theory for this transition remains elusive. Here, we present an effective medium theory (EMT) for this mechanical phase transition of fiber networks. We extend a previous EMT appropriate for linear elasticity to incorporate nonlinear effects via an anharmonic Hamiltonian. The mean-field predictions of this theory, including the critical exponents, scaling relations and non-affine fluctuations qualitatively agree with previous experimental and numerical results.

Effect of connectivity on the elasticity of athermal network materials

Network materials with stochastic structure are ubiquitous in biology and engineering, which drives the current interest in establishing relations between their structure and mechanical behavior. In this work we focus on the effect of connectivity defined by the number of fibers emerging from a crosslink, z, and compare networks with identical (z-homogeneous) and distinct (z-heterogeneous) z at the crosslinks. We observe that the functional form of strain stiffening is z-independent, and that the central z-dependent parameter is the small strain stiffness, E₀. We confirm previous results indicating that the functional form of E₀(z) is a power function with 3 regimes and observe that this applies to a broad range of z. However, the scaling exponents are different in the z-homogeneous and z-heterogeneous cases. We confirm that increasing z across the Maxwell's central force isostatic point leads to a transition from bending to axial energy storage. However, we observe that this does not necessarily imply that deformation becomes affine in the large z limit. In fact, networks of fibers with low bending stiffness retain a relaxation mode based on the rotational degree of freedom of the crosslinks which allows E₀ in the large z limit to be smaller than the affine model prediction. We also conclude that in the z-heterogeneous case, the mean connectivity z̄ is sufficient to evaluate the effect of connectivity on E₀ and that higher moments of the distribution of z are less important.

Mechanics of fiber networks under a bulk strain

Biopolymer networks are common in biological systems from the cytoskeleton of individual cells to collagen in the extracellular matrix. The mechanics of these systems under applied strain can be explained in some cases by a phase transition from soft to rigid states. For collagen networks, it has been shown that this transition is critical in nature and it is predicted to exhibit diverging fluctuations near a critical strain that depends on the network's connectivity and structure. Whereas prior work focused mostly on shear deformation that is more accessible experimentally, here we study the mechanics of such networks under an applied bulk or isotropic extension. We confirm that the bulk modulus of subisostatic fiber networks exhibits similar critical behavior as a function of bulk strain. We find different nonmean-field exponents for bulk as opposed to shear. We also confirm a similar hyperscaling relation to what was previously found for shear.

Rigidity transitions in zero-temperature polygons

We study geometrical clues of a rigidity transition due to the emergence of a system-spanning state of self-stress in underconstrained systems of individual polygons and spring networks constructed from such polygons. When a polygon with harmonic bond edges and an area spring constraint is subject to an expansive strain, we observe that convexity of the polygon is a necessary condition for such a self-stress. We prove that the cyclic configuration of the polygon is a sufficient condition for the self-stress. This correspondence of geometry and rigidity is akin to the straightening of a one dimensional chain of springs to rigidify it. We predict the onset of the rigidity transition and estimate the transition strain using purely geometrical methods. These findings help determine the rigidity of an area-preserving polygon just by knowing its geometry. Since two-dimensional spring networks can be considered as a network of polygons, we look for similar geometric features in underconstrained spring networks under isotropic expansive strain. We observe that all polygons attain convexity at the rigidity transition such that the fraction of convex, but not cyclic, polygons predicts the onset of the rigidity transition. Acyclic polygons in the network correlate with larger tensions, forming effective force chains.

Biomechanical origins of inherent tension in fibrin networks

Blood clots form at the site of vascular injury to seal the wound and prevent bleeding. Clots are in tension as they perform their biological functions and withstand hydrodynamic forces of blood flow, vessel wall fluctuations, extravascular muscle contraction and other forces. There are several mechanisms that generate tension in a blood clot, of which the most well-known is the contraction/retraction caused by activated platelets. Here we show through experiments and modeling that clot tension is generated by the polymerization of fibrin. Our mathematical model is built on the hypothesis that the shape of fibrin monomers having two-fold symmetry and off-axis binding sites is ultimately the source of inherent tension in individual fibers and the clot. As the diameter of a fiber grows during polymerization the fibrin monomers must suffer axial twisting deformation so that they remain in register to form the half-staggered arrangement characteristic of fibrin protofibrils. This deformation results in a pre-strain that causes fiber and network tension. Our results for the pre-strain in single fibrin fibers is in agreement with experiments that measured it by cutting fibers and measuring their relaxed length. We connect the mechanics of a fiber to that of the network using the 8-chain model of polymer elasticity. By combining this with a continuum model of swellable elastomers we can compute the evolution of tension in a constrained fibrin gel. The temporal evolution and tensile stresses predicted by this model are in qualitative agreement with experimental measurements of the inherent tension of fibrin clots polymerized between two fixed rheometer plates. These experiments also revealed that increasing thrombin concentration leads to increasing internal tension in the fibrin network. Our model may be extended to account for other mechanisms that generate pre-strains in individual fibers and cause tension in three-dimensional proteinaceous polymeric networks.

Stiffening of under-constrained spring networks under isotropic strain

Disordered spring networks are a useful paradigm to examine macroscopic mechanical properties of amorphous materials. Here, we study the elastic behavior of under-constrained spring networks, i.e. networks with more degrees of freedom than springs. While such networks are usually floppy, they can be rigidified by applying external strain. Recently, an analytical formalism has been developed to predict the scaling behavior of the elastic network properties close to this rigidity transition. Here we numerically show that these predictions apply to many different classes of spring networks, including phantom triangular, Delaunay, Voronoi, and honeycomb networks. The analytical predictions further imply that the shear modulus G scales linearly with isotropic stress T close to the rigidity transition. However, this seems to be at odds with recent numerical studies suggesting an exponent between G and T that is smaller than one for some network classes. Using increased numerical precision and shear stabilization, we demonstrate here that close to the transition a linear scaling, G ∼ T, holds independent of the network class. Finally, we show that our results are not or only weakly affected by finite-size effects, depending on the network class.

Materials science and mechanosensitivity of living matter

Living systems are composed of molecules that are synthesized by cells that use energy sources within their surroundings to create fascinating materials that have mechanical properties optimized for their biological function. Their functionality is a ubiquitous aspect of our lives. We use wood to construct furniture, bacterial colonies to modify the texture of dairy products and other foods, intestines as violin strings, bladders in bagpipes, and so on. The mechanical properties of these biological materials differ from those of other simpler synthetic elastomers, glasses, and crystals. Reproducing their mechanical properties synthetically or from first principles is still often unattainable. The challenge is that biomaterials often exist far from equilibrium, either in a kinetically arrested state or in an energy consuming active state that is not yet possible to reproduce de novo. Also, the design principles that form biological materials often result in nonlinear responses of stress to strain, or force to displacement, and theoretical models to explain these nonlinear effects are in relatively early stages of development compared to the predictive models for rubberlike elastomers or metals. In this Review, we summarize some of the most common and striking mechanical features of biological materials and make comparisons among animal, plant, fungal, and bacterial systems. We also summarize some of the mechanisms by which living systems develop forces that shape biological matter and examine newly discovered mechanisms by which cells sense and respond to the forces they generate themselves, which are resisted by their environment, or that are exerted upon them by their environment. Within this framework, we discuss examples of how physical methods are being applied to cell biology and bioengineering.

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.

Shear-induced phase transition and critical exponents in three-dimensional fiber networks

When subject to applied strain, fiber networks exhibit nonlinear elastic stiffening. Recent theory and experiments have shown that this phenomenon is controlled by an underlying mechanical phase transition that is critical in nature. Growing simulation evidence points to non-mean-field behavior for this transition and a hyperscaling relation has been proposed to relate the corresponding critical exponents. Here, we report simulations on two distinct network structures in three dimensions. By performing a finite-size scaling analysis, we test hyperscaling and identify various critical exponents. From the apparent validity of hyperscaling, as well as the non-mean-field exponents we observe, our results suggest that the upper critical dimension for the strain-controlled phase transition is above three, in contrast to the jamming transition that represents another athermal, mechanical phase transition.

Long-range mechanical signaling in biological systems

Cells can respond to signals generated by other cells that are remarkably far away. Studies from at least the 1920's showed that cells move toward each other when the distance between them is on the order of a millimeter, which is many times the cell diameter. Chemical signals generated by molecules diffusing from the cell surface would move too slowly and dissipate too fast to account for these effects, suggesting that they might be physical rather than biochemical. The non-linear elastic responses of sparsely connected networks of stiff or semiflexible filament such as those that form the extracellular matrix (ECM) and the cytoskeleton have unusual properties that suggest multiple mechanisms for long-range signaling in biological tissues. These include not only direct force transmission, but also highly non-uniform local deformations, and force-generated changes in fiber alignment and density. Defining how fibrous networks respond to cell-generated forces can help design new methods to characterize abnormal tissues and can guide development of improved biomimetic materials.

The vimentin cytoskeleton: when polymer physics meets cell biology

The proper functions of tissues depend on the ability of cells to withstand stress and maintain shape. Central to this process is the cytoskeleton, comprised of three polymeric networks: F-actin, microtubules, and intermediate filaments (IFs). IF proteins are among the most abundant cytoskeletal proteins in cells; yet they remain some of the least understood. Their structure and function deviate from those of their cytoskeletal partners, F-actin and microtubules. IF networks show a unique combination of extensibility, flexibility and toughness that confers mechanical resilience to the cell. Vimentin is an IF protein expressed in mesenchymal cells. This review highlights exciting new results on the physical biology of vimentin intermediate filaments and their role in allowing whole cells and tissues to cope with stress.

Finite size effects in critical fiber networks

Fibrous networks such as collagen are common in physiological systems. One important function of these networks is to provide mechanical stability for cells and tissues. At physiological levels of connectivity, such networks would be mechanically unstable with only central-force interactions. While networks can be stabilized by bending interactions, it has also been shown that they exhibit a critical transition from floppy to rigid as a function of applied strain. Beyond a certain strain threshold, it is predicted that underconstrained networks with only central-force interactions exhibit a discontinuity in the shear modulus. We study the finite-size scaling behavior of this transition and identify both the mechanical discontinuity and critical exponents in the thermodynamic limit. We find both non-mean-field behavior and evidence for a hyperscaling relation for the critical exponents, for which the network stiffness is analogous to the heat capacity for thermal phase transitions. Further evidence for this is also found in the self-averaging properties of fiber networks.

A minimal-length approach unifies rigidity in underconstrained materials

We present an approach to understand geometric-incompatibility-induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length [Formula: see text], determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate [Formula: see text] from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.