Elastic regimes of subisostatic athermal fiber networks

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.

Nanocolloidal hydrogel mimics the structure and nonlinear mechanical properties of biological fibrous networks

Fibrous networks formed by biological polymers such as collagen or fibrin exhibit nonlinear mechanical behavior. They undergo strong stiffening in response to weak shear and elongational strains, but soften under compressional strain, in striking difference with the response to the deformation of flexible-strand networks formed by molecules. The nonlinear properties of fibrous networks are attributed to the mechanical asymmetry of the constituent filaments, for which a stretching modulus is significantly larger than the bending modulus. Studies of the nonlinear mechanical behavior are generally performed on hydrogels formed by biological polymers, which offers limited control over network architecture. Here, we report an engineered covalently cross-linked nanofibrillar hydrogel derived from cellulose nanocrystals and gelatin. The variation in hydrogel composition provided a broad-range change in its shear modulus. The hydrogel exhibited both shear-stiffening and compression-induced softening, in agreement with the predictions of the affine model. The threshold nonlinear stress and strain were universal for the hydrogels with different compositions, which suggested that nonlinear mechanical properties are general for networks formed by rigid filaments. The experimental results were in agreement with an affine model describing deformation of the network formed by rigid filaments. Our results lend insight into the structural features that govern the nonlinear biomechanics of fibrous networks and provide a platform for future studies of the biological impact of nonlinear mechanical properties.

Mechanical criticality of fiber networks at a finite temperature

At zero temperature, spring networks with connectivity below Maxwell's isostatic threshold undergo a mechanical phase transition from a floppy state at small strains to a rigid state for applied shear strain above a critical strain threshold. Disordered networks in the floppy mechanical regime can be stabilized by entropic effects at finite temperature. We develop a scaling theory for this mechanical phase transition at finite temperature, yielding relationships between various scaling exponents. Using Monte Carlo simulations, we verify these scaling relations and identify anomalous entropic elasticity with sublinear T dependence in the linear elastic regime. While our results are consistent with prior studies of phase behavior near the isostatic point, the present work also makes predictions relevant to the broad class of disordered thermal semiflexible polymer networks for which the connectivity generally lies far below the isostatic threshold.

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.

Stiffening mechanisms in stochastic athermal fiber networks

Stochastic athermal networks composed of fibers that deform axially and in bending strain stiffen much faster than thermal networks of axial elements, such as elastomers. Here we investigate the physical origin of stiffening in athermal network materials. To this end, we use models of stochastic networks subjected to uniaxial deformation and identify the emergence of two subnetworks, the stress path subnetwork (SPSN) and the bending support subnetwork (BSSN), which carry most of the axial and bending energies, respectively. The BSSN controls lateral contraction and modulates the organization of the SPSN during deformation. The SPSN is preferentially oriented in the loading direction, while the BSSN's preferential orientation is orthogonal to the SPSN. In nonaffine networks stiffening is exponential, while in close-to-affine networks it is quadratic. The difference is due to a much more modest lateral contraction in the approximately affine case and to a stiffer BSSN. Exponential stiffening emerges from the interplay of the axial and bending deformation modes at the scale of individual or small groups of fibers undergoing large deformations and being subjected to the constraint of rigid cross-links, and it is not necessarily a result of complex interactions involving many connected fibers. An apparent third regime of quadratic stiffening may be evidenced in nonaffinely deforming networks provided the nominal stress is observed. This occurs at large stretches, when the BSSN contribution of stiffening vanishes. However, this regime is not present if the Cauchy stress is used, in which case stiffening is exponential throughout the entire deformation. These results shed light on the physical nature of stiffening in a broad class of materials including connective tissue, the extracellular matrix, nonwovens, felt, and other athermal network materials.

Toughness of Network Materials: Structural Parameters Controlling Damage Accumulation

Many materials have a network of fibers as their main structural component and are referred to as network materials. Their strength and toughness are important in both engineering and biology. In this work we consider stochastic model fiber networks without pre-existing cracks and study their rupture mechanism. These materials soften as the crosslinks or fibers fail and exhibit either brittle failure immediately after the peak stress, or a more gradual, ductile rupture in the post peak regime. We observe that ductile failure takes place at constant energy release rate defined in the absence of pre-existing cracks as the strain derivative of the specific energy released. The network parameters controlling the energy release rate are identified and discussed in relation to the Lake-Thomas theory which applies to crack growth situations. We also observe a ductile to brittle failure transition as the network becomes more affine and relate the embrittlement to the reduction of mechanical heterogeneity of the network. Further, we confirm previous reports that the network strength scales linearly with the bond strength and with the crosslink density. The present results extend the Lake-Thomas theory to networks without pre-existing cracks which fail by the gradual accumulation of distributed damage and contribute to the development of a physical picture of failure in stochastic network materials.

Nonaffine Deformation of Semiflexible Polymer and Fiber Networks

Networks of semiflexible or stiff polymers such as most biopolymers are known to deform inhomogeneously when sheared. The effects of such nonaffine deformation have been shown to be much stronger than for flexible polymers. To date, our understanding of nonaffinity in such systems is limited to simulations or specific 2D models of athermal fibers. Here, we present an effective medium theory for nonaffine deformation of semiflexible polymer and fiber networks, which is general to both 2D and 3D and in both thermal and athermal limits. The predictions of this model are in good agreement with both prior computational and experimental results for linear elasticity. Moreover, the framework we introduce can be extended to address nonlinear elasticity and network dynamics.

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.

Multiscale architecture: Mechanics of composite cytoskeletal networks

Different types of biological cells respond differently to mechanical stresses, and these responses are mainly governed by the cytoskeleton. The main components of this biopolymer network are actin filaments, microtubules, and intermediate filaments, whose mechanical and dynamic properties are highly distinct, thus opening up a large mechanical parameter space. Aside from experiments on whole, living cells, "bottom-up" approaches, utilizing purified, reconstituted protein systems, tremendously help to shed light on the complex mechanics of cytoskeletal networks. Such experiments are relevant in at least three aspects: (i) from a fundamental point of view, cytoskeletal networks provide a perfect model system for polymer physics; (ii) in materials science and "synthetic cell" approaches, one goal is to fully understand properties of cellular materials and reconstitute them in synthetic systems; (iii) many diseases are associated with cell mechanics, so a thorough understanding of the underlying phenomena may help solving pressing biomedical questions. In this review, we discuss the work on networks consisting of one, two, or all three types of filaments, entangled or cross-linked, and consider active elements such as molecular motors and dynamically growing filaments. Interestingly, tuning the interactions among the different filament types results in emergent network properties. We discuss current experimental challenges, such as the comparability of different studies, and recent methodological advances concerning the quantification of attractive forces between filaments and their influence on network mechanics.

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.

Stretchy and disordered: Toward understanding fracture in soft network materials via mesoscopic computer simulations

Soft network materials exist in numerous forms ranging from polymer networks, such as elastomers, to fiber networks, such as collagen. In addition, in colloidal gels, an underlying network structure can be identified, and several metamaterials and textiles can be considered network materials as well. Many of these materials share a highly disordered microstructure and can undergo large deformations before damage becomes visible at the macroscopic level. Despite their widespread presence, we still lack a clear picture of how the network structure controls the fracture processes of these soft materials. In this Perspective, we will focus on progress and open questions concerning fracture at the mesoscopic scale, in which the network architecture is clearly resolved, but neither the material-specific atomistic features nor the macroscopic sample geometries are considered. We will describe concepts regarding the network elastic response that have been established in recent years and turn out to be pre-requisites to understand the fracture response. We will mostly consider simulation studies, where the influence of specific network features on the material mechanics can be cleanly assessed. Rather than focusing on specific systems, we will discuss future challenges that should be addressed to gain new fundamental insights that would be relevant across several examples of soft network materials.

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.

Physically-based structural modeling of a typical regenerative tissue analog bridges material macroscale continuum and cellular microscale discreteness and elucidates the hierarchical characteristics of cell-matrix interaction

This paper presents a comprehensive physically-based structural modelling for the passive and active biomechanical processes in a typical engineered tissue - namely, cell-compacted collagen gel. First, it introduces a sinusoidal curve analog for quantifying the mechanical response of the collagen fibrils and a probability distribution function of the characteristic crimp ratio for taking into account the fibrillar geometric entropic effect. The constitutive framework based on these structural characteristics precisely reproduces the nonlinearity, the viscoelasticity, and fairly captures the Poisson effect exhibiting in the macroscale tensile tests; which, therefore, substantially validates the structural modelling for the analysis of the cell-gel interaction during collagen gel compaction. Second, a deterministic molecular clutch model specific to the interaction between the cell pseudopodium and the collagen network is developed, which emphasizes the dependence of traction force on clutch number altering with the retrograde flow velocity, actin polymeric velocity, and the deformation of the stretched fibril. The modelling reveals the hierarchical features of cellular substrate sensing, i.e. a biphasic traction force response to substrate elasticity begins at the level of individual fibrils and develops into the second biphasic sensing by means of the fibrillar number integration at the whole-cell level. Singular in crossing the realms of continuum and discrete mechanics, the methodologies developed in this study for modelling the filamentous materials and cell-fibril interaction deliver deep insight into the temporospatially dynamic 3D cell-matrix interaction, and are able to bridge the cellular microscale and material macroscale in the exploration of related topics in mechanobiology.

Anomalous mechanics of Zn2+-modified fibrin networks

Fibrin is the main component of blood clots. The mechanical properties of fibrin are therefore of critical importance in successful hemostasis. One of the divalent cations released by platelets during hemostasis is Zn2+; however, its effect on the network structure of fibrin gels and on the resultant mechanical properties remains poorly understood. Here, by combining mechanical measurements with three-dimensional confocal microscopy imaging, we show that Zn2+ can tune the fibrin network structure and alter its mechanical properties. In the presence of Zn2+, fibrin protofibrils form large bundles that cause a coarsening of the fibrin network due to an increase in fiber diameter and reduction of the total fiber length. We further show that the protofibrils in these bundles are loosely coupled to one another, which results in a decrease of the elastic modulus with increasing Zn2+ concentrations. We explore the elastic properties of these networks at both low and high stress: At low stress, the elasticity originates from pulling the thermal slack out of the network, and this is consistent with the thermal bending of the fibers. By contrast, at high stress, the elasticity exhibits a common master curve consistent with the stretching of individual protofibrils. These results show that the mechanics of a fibrin network are closely correlated with its microscopic structure and inform our understanding of the structure and physical mechanisms leading to defective or excessive clot stiffness.

Strength of stochastic fibrous materials under multiaxial loading

Many biological and engineering materials are made from fibers organized in the form of a stochastic crosslinked network, and the mechanics of the network controls the behavior of the material. In this work we investigate the strength of stochastic networks without pre-existing damage which fail due to crosslink rupture. Athermal networks ranging from approximately affine to strongly non-affine are subjected to multiaxial loading and the strength is evaluated using numerical models. It is observed that once the stress is normalized by the strength measured in uniaxial tension, the failure surface becomes approximately independent of network parameters. This extends the relation between strength and network parameters previously established in (S. Deogekar, M. R. Islam, R. C. Picu, Parameters controlling the strength of stochastic fibrous materials, Int. J. Solids Struct., 2019, 168, 194-202) to the multiaxial case. The failure surface depends on both first two invariants of the stress. Strongly non-affine networks behave somewhat different from the affine networks under loadings close to the hydrostatic and pure shear loading modes, while the difference disappears in the first quadrant of the principal stress space. The results are compared with experimental data from the literature.

Nonlinear Mechanical Properties of Prestressed Branched Fibrous Networks

Random fiber networks constitute the solid skeleton of many biological materials such as the cytoskeleton of cells and extracellular matrix of soft tissues. These random networks show unique mechanical properties such as nonlinear shear strain-stiffening and strain softening when subjected to preextension and precompression, respectively. In this study, we perform numerical simulations to characterize the influence of axial prestress on the nonlinear mechanical response of random network structures as a function of their micromechanical and geometrical properties. We build our numerical network models using the microstructure of disordered hexagonal lattices and quantify their nonlinear shear response as a function of uniaxial prestress strain. We consider three different material models for individual fibers and fully characterize their influence on the mechanical response of prestressed networks. Moreover, we investigate both the influence of geometric disorder keeping the network connectivity constant and the influence of the randomness in the stiffness of individual fibers keeping their mean stiffness constant. The effects of network connectivity and bending rigidity of fibers are also determined. Several important conclusions are made, including that the tensile and compressive prestress strains, respectively, increase and decrease the initial network shear stiffness but have no effect on the maximal shear modulus. We discuss the findings in terms of microstructural properties such as the local strain energy distribution.

Strain Stiffening and Negative Normal Force of Agarose Hydrogels

Inspired by the specific strain stiffening and negative normal force phenomena in several biological networks, herein, we show strain stiffening and negative normal force in agarose hydrogels. We use both pre-strain and strain amplitude sweep protocols in dynamic rheological measurements where the gel slip was suppressed by the in situ gelation in the cross-hatched parallel plate rheometer geometry. Within the stiffening region, we show the scaling relation for the differential modulus K ∝ σ1, where σ is stress. The strain at the onset of stiffening is almost constant throughout the concentration range. The gels show negative apparent normal stress difference when sheared as a result of the gel contraction. The pore size of the hydrogel is large enough to allow water to move with respect to the network to balance the pressure difference caused by the hoop stress. The rheological analysis together with scanning electron microscopy suggests that the agarose gels can be described using subisostatic athermal network models where the connectivity dictates the stiffening behavior. Therefore, the simple agarose gels appear to capture several of the viscoelastic properties, which were previously thought to be characteristic to biological protein macromolecules.

Mechanical properties of subisostatic random networks composed of nonlinear fibers

Fibrous protein networks provide structural integrity to different biological materials such as soft tissues. These networks display an unusual exponential strain-stiffening behavior when subjected to mechanical loads. This nonlinear strain-stiffening behavior has so far been explained in terms of the network microstructure and the flexibility of constituting fibers. Here, we conduct a comprehensive computational study to characterize the importance of material properties of individual fibers in the overall nonlinear mechanical response of random fiber networks. To this end, we consider three nonlinear material models, ranging from an almost linear form to a highly nonlinear one, for the fibers of subisostatic disordered networks. We characterize the amount of strain-stiffening as a function of bending rigidity of the fibers, the amount of nonlinearity of the fibers, and the connectivity of random networks. We find that networks composed of highly nonlinear fibers exhibit much more strain-stiffening than networks made up of linear fibers. Furthermore, the local strain distribution becomes more homogenous as the amount of nonlinearity in the material models increases. Increasing the network connectivity signifies the importance of the nonlinear material response of individual fibers in the overall mechanical behavior of networks. The constitutive behavior of fibers plays an important role in defining the failure response of networks particularly in the damage initiation and evolution. These important findings for how the mechanical response of individual fibers affects the overall mechanical properties of random networks could find applications in designing new biomimetic materials and could help scientists better understand the mechanical properties of biological materials.

Nonlinear Poisson Effect Governed by a Mechanical Critical Transition

Under extensional strain, fiber networks can exhibit an anomalously large and nonlinear Poisson effect accompanied by a dramatic transverse contraction and volume reduction for applied strains as small as a few percent. We demonstrate that this phenomenon is controlled by a collective mechanical phase transition that occurs at a critical uniaxial strain that depends on network connectivity. This transition is punctuated by an anomalous peak in the apparent Poisson's ratio and other critical signatures such as diverging nonaffine strain fluctuations.

Parameters controlling the strength of stochastic fibrous materials

Many materials of everyday use are fibrous and their strength is important in most applications. In this work we study the dependence of the strength of random fiber networks on structural parameters such as the network density, cross-link density, fiber tortuosity, and the strength of the inter-fiber cross-links. Athermal networks of cellular and fibrous type are considered. We conclude that the network strength scales linearly with the cross-link number density and with the cross-link strength for a broad range of network parameters, and for both types of networks considered. Network strength is independent of fiber material properties and of fiber tortuosity. This information can be used to design fiber networks for specified strength and, generally, to understand the mechanical behavior of fibrous materials.

Strain-Stiffening of Agarose Gels

Strain-stiffening is one of the characteristic properties of biological hydrogels and extracellular matrices, where the stiffness increases upon increased deformation. Whereas strain-stiffening is ubiquitous in protein-based materials, it has been less observed for polysaccharide and synthetic polymer gels. Here we show that agarose, that is, a common linear polysaccharide, forms helical fibrillar bundles upon cooling from aqueous solution. The hydrogels with these semiflexible fibrils show pronounced strain-stiffening. However, to reveal strain-stiffening, suppressing wall slippage turned as untrivial. Upon exploring different sample preparation techniques and rheological architectures, the cross-hatched parallel plate geometries and in situ gelation in the rheometer successfully prevented the slippage and resolved the strain-stiffening behavior. Combining with microscopy, we conclude that strain-stiffening is due to the semiflexible nature of the agarose fibrils and their geometrical connectivity, which is below the central-force isostatic critical connectivity. The biocompatibility and the observed strain-stiffening suggest the potential of agarose hydrogels in biomedical applications.

Stress-stabilized subisostatic fiber networks in a ropelike limit

The mechanics of disordered fibrous networks such as those that make up the extracellular matrix are strongly dependent on the local connectivity or coordination number. For biopolymer networks this coordination number is typically between 3 and 4. Such networks are sub-isostatic and linearly unstable to deformation with only central force interactions, but exhibit a mechanical phase transition between floppy and rigid states under strain. The introduction of weak bending interactions stabilizes these networks and suppresses the critical signatures of this transition. We show that applying external stress can also stabilize subisostatic networks with only tensile central force interactions, i.e., a ropelike potential. Moreover, we find that the linear shear modulus shows a power-law scaling with the external normal stress, with a non-mean-field exponent. For networks with finite bending rigidity, we find that the critical stain shifts to lower values under prestress.

Strong triaxial coupling and anomalous Poisson effect in collagen networks

While cells within tissues generate and sense 3D states of strain, the current understanding of the mechanics of fibrous extracellular matrices (ECMs) stems mainly from uniaxial, biaxial, and shear tests. Here, we demonstrate that the multiaxial deformations of fiber networks in 3D cannot be inferred solely based on these tests. The interdependence of the three principal strains gives rise to anomalous ratios of biaxial to uniaxial stiffness between 8 and 9 and apparent Poisson's ratios larger than 1. These observations are explained using a microstructural network model and a coarse-grained constitutive framework that predicts the network Poisson effect and stress-strain responses in uniaxial, biaxial, and triaxial modes of deformation as a function of the microstructural properties of the network, including fiber mechanics and pore size of the network. Using this theoretical approach, we found that accounting for the Poisson effect leads to a 100-fold increase in the perceived elastic stiffness of thin collagen samples in extension tests, reconciling the seemingly disparate measurements of the stiffness of collagen networks using different methods. We applied our framework to study the formation of fiber tracts induced by cellular forces. In vitro experiments with low-density networks showed that the anomalous Poisson effect facilitates higher densification of fibrous tracts, associated with the invasion of cancerous acinar cells. The approach developed here can be used to model the evolving mechanics of ECM during cancer invasion and fibrosis.

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.

Normal stress anisotropy and marginal stability in athermal elastic networks

Hydrogels of semiflexible biopolymers such as collagen have been shown to contract axially under shear strain, in contrast to the axial dilation observed for most elastic materials. Recent work has shown that this behavior can be understood in terms of the porous, two-component nature and consequent time-dependent compressibility of hydrogels. The apparent normal stress measured by a torsional rheometer reflects only the tensile contribution of the axial component σzz on long (compressible) timescales, crossing over to the first normal stress difference, N1 = σxx - σzz at short (incompressible) times. While the behavior of N1 is well understood for isotropic viscoelastic materials undergoing affine shear deformation, biopolymer networks are often anisotropic and deform nonaffinely. Here, we numerically study the normal stresses that arise under shear in subisostatic, athermal semiflexible polymer networks. We show that such systems exhibit strong deviations from affine behavior and that these anomalies are controlled by a rigidity transition as a function of strain.

The Role of Network Architecture in Collagen Mechanics

Collagen forms fibrous networks that reinforce tissues and provide an extracellular matrix for cells. These networks exhibit remarkable strain-stiffening properties that tailor the mechanical functions of tissues and regulate cell behavior. Recent models explain this nonlinear behavior as an intrinsic feature of disordered networks of stiff fibers. Here, we experimentally validate this theoretical framework by measuring the elastic properties of collagen networks over a wide range of self-assembly conditions. We show that the model allows us to quantitatively relate both the linear and nonlinear elastic behavior of collagen networks to their underlying architecture. Specifically, we identify the local coordination number (or connectivity) 〈z〉 as a key architectural parameter that governs the elastic response of collagen. The network elastic response reveals that 〈z〉 decreases from 3.5 to 3 as the polymerization temperature is raised from 26 to 37°C while being weakly dependent on concentration. We furthermore infer a Young's modulus of 1.1 MPa for the collagen fibrils from the linear modulus. Scanning electron microscopy confirms that 〈z〉 is between three and four but is unable to detect the subtle changes in 〈z〉 with polymerization conditions that rheology is sensitive to. Finally, we show that, consistent with the model, the initial stress-stiffening response of collagen networks is controlled by the negative normal stress that builds up under shear. Our work provides a predictive framework to facilitate future studies of the regulatory effect of extracellular matrix molecules on collagen mechanics. Moreover, our findings can aid mechanobiological studies of wound healing, fibrosis, and cancer metastasis, which require collagen matrices with tunable mechanical properties.

Hindrances to precise recovery of cellular forces in fibrous biopolymer networks

How cells move through the three-dimensional extracellular matrix (ECM) is of increasing interest in attempts to understand important biological processes such as cancer metastasis. Just as in motion on flat surfaces, it is expected that experimental measurements of cell-generated forces will provide valuable information for uncovering the mechanisms of cell migration. However, the recovery of forces in fibrous biopolymer networks may suffer from large errors. Here, within the framework of lattice-based models, we explore possible issues in force recovery by solving the inverse problem: how can one determine the forces cells exert to their surroundings from the deformation of the ECM? Our results indicate that irregular cell traction patterns, the uncertainty of local fiber stiffness, the non-affine nature of ECM deformations and inadequate knowledge of network topology will all prevent the precise force determination. At the end, we discuss possible ways of overcoming these difficulties.

Programming the mechanics of cohesive fiber networks by compression

Fibrous networks are ideal functional materials since they provide mechanical rigidity at low weight. Here, we demonstrate that fibrous networks of the blood clotting protein fibrin undergo a strong and irreversible increase in their mechanical rigidity in response to uniaxial compression. This rigidification can be precisely controlled by the level of applied compressive strain, providing a means to program the network rigidity without having to change its composition. To identify the underlying mechanism we measure single fiber-fiber interactions using optical tweezers. We further develop a minimal computational model of cohesive fiber networks that shows that stiffening arises due to the formation of new bonds in the compressed state, which develop tensile stress when the network is re-expanded. The model predicts that the network stiffness after a compression cycle obeys a power-law dependence on tensile stress, which we confirm experimentally. This finding provides new insights into how biological tissues can adapt themselves independently of any cellular processes, offering new perspectives to inspire the design of reprogrammable materials.

Stress-induced plasticity of dynamic collagen networks

The structure and mechanics of tissues is constantly perturbed by endogenous forces originated from cells, and at the same time regulate many important cellular functions such as migration, differentiation, and growth. Here we show that 3D collagen gels, major components of connective tissues and extracellular matrix (ECM), are significantly and irreversibly remodeled by cellular traction forces, as well as by macroscopic strains. To understand this ECM plasticity, we develop a computational model that takes into account the sliding and merging of ECM fibers. We have confirmed the model predictions with experiment. Our results suggest the profound impacts of cellular traction forces on their host ECM during development and cancer progression, and suggest indirect mechanical channels of cell-cell communications in 3D fibrous matrices.The structure and mechanics of tissues is constantly perturbed by endogenous forces originated from cells. Here the authors show that 3D collagen gels, major components of connective tissues and extracellular matrix, are significantly and irreversibly remodelled by cellular traction forces and by macroscopic strains.

Nonlinear Actin Deformations Lead to Network Stiffening, Yielding, and Nonuniform Stress Propagation

We use optical tweezers microrheology and fluorescence microscopy to apply nonlinear microscale strains to entangled and cross-linked actin networks, and measure the resulting stress and actin filament deformations. We couple nonlinear stress response and relaxation to the velocities and displacements of individual fluorescent-labeled actin segments, at varying times throughout the strain and varying distances from the strain path, to determine the underlying molecular dynamics that give rise to the debated nonlinear response and stress propagation of cross-linked and entangled actin networks at the microscale. We show that initial stress stiffening arises from acceleration of strained filaments due to molecular extension along the strain, while softening and yielding is coupled to filament deceleration, halting, and recoil. We also demonstrate a surprising nonmonotonic dependence of filament deformation on cross-linker concentration. Namely, networks with no cross-links or substantial cross-links both exhibit fast initial filament velocities and reduced molecular recoil while intermediate cross-linker concentrations display reduced velocities and increased recoil. We show that these collective results are due to a balance of network elasticity and force-induced cross-linker unbinding and rebinding. We further show that cross-links dominate entanglement dynamics when the length between cross-linkers becomes smaller than the length between entanglements. In accord with recent simulations, we demonstrate that post-strain stress can be long-lived in cross-linked networks by distributing stress to a small fraction of highly strained connected filaments that span the network and sustain the load, thereby allowing the rest of the network to recoil and relax.

A novel microstructural interpretation for the biomechanics of mouse skin derived from multiscale characterization

Skin is a complex, multi-layered organ, with important functions in the protection of the body. The dermis provides structural support to the epidermal barrier, and thus has attracted a large number of mechanical studies. As the dermis is made of a mixture of stiff fibres embedded in a soft non-fibrillar matrix, it is classically considered that its mechanical response is based on an initial alignment of the fibres, followed by the stretching of the aligned fibres. Using a recently developed set-up combining multiphoton microscopy with mechanical assay, we imaged the fibres network evolution during dermis stretching. These observations, combined with a wide set of mechanical tests, allowed us to challenge the classical microstructural interpretation of the mechanical properties of the dermis: we observed a continuous alignment of the collagen fibres along the stretching. All our results can be explained if each fibre contributes by a given stress to the global response. This plastic response is likely due to inner sliding inside each fibre. The non-linear mechanical response is due to structural effects of the fibres network in interaction with the surrounding non-linear matrix. This multiscale interpretation explains our results on genetically-modified mice with a simple alteration of the dermis microstructure.

STATEMENT OF SIGNIFICANCE

Soft tissues, as skin, tendon or aorta, are made of extra-cellular matrix, with very few cells embedded inside. The matrix is a mixture of water and biomolecules, which include the collagen fibre network. The role of the collagen is fundamental since the network is supposed to control the tissue mechanical properties and remodeling: the cells attach to the collagen fibres and feel the network deformations. This paper challenges the classical link between fibres organization and mechanical properties. To do so, it uses multiscale observations combined to a large set of mechanical loading. It thus appears that the behaviour at low stretches is mostly controlled by the network structural response, while, at large stretches, the fibre inner-sliding dominate.

Strain-driven criticality underlies nonlinear mechanics of fibrous networks

Networks with only central force interactions are floppy when their average connectivity is below an isostatic threshold. Although such networks are mechanically unstable, they can become rigid when strained. It was recently shown that the transition from floppy to rigid states as a function of simple shear strain is continuous, with hallmark signatures of criticality [Sharma et al., Nature Phys. 12, 584 (2016)1745-247310.1038/nphys3628]. The nonlinear mechanical response of collagen networks was shown to be quantitatively described within the framework of such mechanical critical phenomenon. Here, we provide a more quantitative characterization of critical behavior in subisostatic networks. Using finite-size scaling we demonstrate the divergence of strain fluctuations in the network at well-defined critical strain. We show that the characteristic strain corresponding to the onset of strain stiffening is distinct from but related to this critical strain in a way that depends on critical exponents. We confirm this prediction experimentally for collagen networks. Moreover, we find that the apparent critical exponents are largely independent of the spatial dimensionality. With subisostaticity as the only required condition, strain-driven criticality is expected to be a general feature of biologically relevant fibrous networks.

Elasticity of fibrous networks under uniaxial prestress

We present theoretical and experimental studies of the elastic response of fibrous networks subjected to uniaxial strain. Uniaxial compression or extension is applied to extracellular networks of fibrin and collagen using a shear rheometer with free water in/outflow. Both uniaxial stress and the network shear modulus are measured. Prior work [van Oosten, et al., Sci. Rep., 2015, 6, 19270] has shown softening/stiffening of these networks under compression/extension, together with a nonlinear response to shear, but the origin of such behaviour remains poorly understood. Here, we study how uniaxial strain influences the nonlinear mechanics of fibrous networks. Using a computational network model with bendable and stretchable fibres, we show that the softening/stiffening behaviour can be understood for fixed lateral boundaries in 2D and 3D networks with comparable average connectivities to the experimental extracellular networks. Moreover, we show that the onset of stiffening depends strongly on the imposed uniaxial strain. Our study highlights the importance of both uniaxial strain and boundary conditions in determining the mechanical response of hydrogels.

Nonlinear Mechanics of Athermal Branched Biopolymer Networks

Naturally occurring biopolymers such as collagen and actin form branched fibrous networks. The average connectivity in branched networks is generally below the isostatic threshold at which central force interactions marginally stabilize the network. In the submarginal regime, for connectivity below this threshold, such networks are unstable toward small deformations unless stabilized by additional interactions such as bending. Here we perform a numerical study on the elastic behavior of such networks. We show that the nonlinear mechanics of branched networks is qualitatively similar to that of filamentous networks with freely hinged cross-links. In agreement with a recent theoretical study,1 we find that branched networks also exhibit nonlinear mechanics consistent with athermal critical phenomena controlled by strain. We obtain the critical exponents capturing the nonlinear elastic behavior near the critical point by performing scaling analysis of the stiffening curves. We find that the exponents evolve with the connectivity in the network. We show that the nonlinear mechanics of disordered networks, independent of the detailed microstructure, can be characterized by a strain-driven second-order phase transition, and that the primary quantitative differences among different architectures are in the critical exponents describing the transition.