Compressive instabilities enable cell-induced extreme densification patterns in the fibrous extracellular matrix: Discrete model predictions

We present a new model and extensive computations that explain the dramatic remodelling undergone by a fibrous collagen extracellular matrix (ECM), when subjected to contractile mechanical forces from embedded cells or cell clusters. This remodelling creates complex patterns, comprising multiple narrow localised bands of severe densification and fiber alignment, extending far into the ECM, often joining distant cells or cell clusters (such as tumours). Most previous models cannot capture this behaviour, as they assume stable mechanical fiber response with stress an increasing function of fiber stretch, and a restriction to small displacements. Our fully nonlinear network model distinguishes between two types of single-fiber nonlinearity: fibers that undergo stable (supercritical) buckling (as in previous work) versus fibers that suffer unstable (subcritical) buckling collapse. The model allows unrestricted, arbitrarily large displacements (geometric nonlinearity). Our assumptions on single-fiber instability are supported by recent simulations and experiments on buckling of individual beams with a hierarchical microstructure, such as collagen fibers. We use simple scenarios to illustrate, for the first time, two distinct compressive-instability mechanisms at work in our model: unstable buckling collapse of single fibers, and snap-through of multiple-fiber groups. The latter is possible even when single fibers are stable. Through simulations of large fiber networks, we show how these instabilities lead to spatially extended patterns of densification, fiber alignment and ECM remodelling induced by cell contraction. Our model is simple, but describes a very complex, multi-stable energy landscape, using sophisticated numerical optimisation methods that overcome the difficulties caused by instabilities in large systems. Our work opens up new ways of understanding the unique biomechanics of fibrous-network ECM, by fully accounting for nonlinearity and associated loss of stability in fiber networks. Our results provide new insights on tumour invasion and metastasis.

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.

Micromechanics of fibrous scaffolds and their stiffness sensing by cells

Mechanical properties of the tissue engineering scaffolds are known to play a crucial role in cell response. Therefore, an understanding of the cell-scaffold interactions is of high importance. Here, we have utilized discrete fiber network model to quantitatively study the micromechanics of fibrous scaffolds with different fiber arrangements and cross-linking densities. We observe that localized forces on the scaffold result in its anisotropic deformation even for isotropic fiber arrangements. We also see an exponential decay of the displacement field with distance from the location of applied force. This nature of the decay allows us to estimate the characteristic length for force transmission in fibrous scaffolds. Furthermore, we also looked at the stiffness sensing of fibrous scaffolds by individual cells and its dependence on the cellular sensing mechanism. For this, we considered two conditions- stress-controlled, and strain-controlled application of forces by a cell. With fixed strain, we find that the stiffness sensed by a cell is proportional to the scaffold's 'macroscopic' elastic modulus. However, under fixed stress application by the cell, the stiffness sensed by the cell also depends on the cell's own stiffness. In fact, the stiffness values for the same scaffold sensed by the stiff and soft cells can differ from each other by an order of magnitude. The insights from this work will help in designing tissue engineering scaffolds for applications where mechanical stimuli are a critical factor.

Probing soft fibrous materials by indentation

Indentation is often used to measure the stiffness of soft materials whose main structural component is a network of filaments, such as the cellular cytoskeleton, connective tissue, gels, and the extracellular matrix. For elastic materials, the typical procedure requires fitting the experimental force-displacement curve with the Hertz model, which predicts that f=kδ1.5 and k is proportional to the reduced modulus of the indented material, E/(1-ν2). Here we show using explicit models of fiber networks that the Hertz model applies to indentation in network materials provided the indenter radius is larger than approximately 12lc, where lc is the mean segment length of the network. Using smaller indenters leads to a relation between force and indentation displacement of the form f=kδq, where q is observed to increase with decreasing indenter radius. Using the Hertz model to interpret results of indentations in network materials using small indenters leads to an inferred modulus smaller than the real modulus of the material. The origin of this departure from the classical Hertz model is investigated. A compacted, stiff network region develops under the indenter, effectively increasing the indenter size and modifying its shape. This modification is marginal when large indenters are used. However, when the indenter radius is small, the effect of the compacted layer is pronounced as it changes the indenter profile from spherical towards conical. This entails an increase of exponent q above the value of 1.5 corresponding to spherical indenters. STATEMENT OF SIGNIFICANCE: The article presents a study of indentation in network biomaterials and demonstrates a size effect which precludes the use of the Hertz model to infer the elastic constants of the material. The size effect occurs once the indenter radius is smaller than approximately 12 times the mean segment length of the network. This result provides guidelines for the selection of indentation conditions that guarantee the applicability of the Hertz model. At the same time, the finding may be used to infer the mean segment length of the network based on indentations with indenters of various sizes. Hence, the method can be used to evaluate this structural parameter which is not easily accessible in experiments.

Effect of hyaluronic acid on microscale deformations of collagen gels

As fibrous collagen is the most abundant protein in mammalian tissues, gels of collagen fibers have been extensively used as an extracellular matrix scaffold to study how cells sense and respond to cues from their microenvironment. Other components of native tissues, such as glycosaminoglycans like hyaluronic acid, can affect cell behavior in part by changing the mechanical properties of the collagen gel. Prior studies have quantified the effects of hyaluronic acid on the mechanical properties of collagen gels in experiments of uniform shear or compression at the macroscale. However, there remains a lack of experimental studies of how hyaluronic acid changes the mechanical properties of collagen gels at the scale of a cell. Here, we studied how addition of hyaluronic acid to gels of collagen fibers affects the local field of displacements in response to contractile loads applied on length scales similar to those of a contracting cell. Using spherical poly(N-isopropylacrylamide) particles, which contract when heated, we induced displacement in gels of collagen and collagen with hyaluronic acid. Displacement fields were quantified using a combination of confocal microscopy and digital image correlation. Results showed that hyaluronic acid suppressed the distance over which displacements propagated, suggesting that it caused the network to become more linear. Additionally, hyaluronic acid had no statistical effect on heterogeneity of the displacement fields, but it did make the gels more elastic by substantially reducing the magnitude of permanent deformations. Lastly, we examined the effect of hyaluronic acid on fiber remodeling due to localized forces and found that hyaluronic acid partially - but not fully - inhibited remodeling. This result is consistent with prior studies suggesting that fiber remodeling is associated with a phase transition resulting from an instability caused by nonlinearity of the collagen gel.

Random Fiber Network Loaded by a Point Force

This article presents the displacement field produced by a point force acting on an athermal random fiber network (the Green function for the network). The problem is defined within the limits of linear elasticity, and the field is obtained numerically for nonaffine networks characterized by various parameter sets. The classical Green function solution applies at distances from the point force larger than a threshold which is independent of the network parameters in the range studied. At smaller distances, the nonlocal nature of fiber interactions modifies the solution.

Effect of matrix heterogeneity on cell mechanosensing

Cells sense mechanical signals within the extracellular matrix, the most familiar being stiffness, but matrix stiffness cannot be simply described by a single value. Randomness in matrix structure causes stiffness at the scale of a cell to vary by more than an order of magnitude. Additionally, the extracellular matrix contains ducts, blood vessels, and, in cancer or fibrosis, regions with abnormally high stiffness. These different features could alter the stiffness sensed by a cell, but it is unclear whether the change in stiffness is large enough to overcome the noise caused by heterogeneity due to the random fibrous structure. Here we used a combination of experiments and modeling to determine the extent to which matrix heterogeneity disrupts the potential for cell sensing of a locally stiff feature in the matrix. Results showed that, at the scale of a single cell, spatial heterogeneity in local stiffness was larger than the increase in stiffness due to a stiff feature. The heterogeneity was reduced only for large length scales compared to the fiber length. Experiments verified this conclusion, showing spheroids of cells, which were large compared to the average fiber length, spreading preferentially toward stiff inclusions. Hence, the propagation of mechanical cues through the matrix depends on length scale, with single cells being able to sense only the stiffness of the nearby fibers and multicellular structures, such as tumors, also sensing the stiffness of distant matrix features.

Collective forces of tumor spheroids in three-dimensional biopolymer networks

We describe a method for quantifying the contractile forces that tumor spheroids collectively exert on highly nonlinear three-dimensional collagen networks. While three-dimensional traction force microscopy for single cells in a nonlinear matrix is computationally complex due to the variable cell shape, here we exploit the spherical symmetry of tumor spheroids to derive a scale-invariant relationship between spheroid contractility and the surrounding matrix deformations. This relationship allows us to directly translate the magnitude of matrix deformations to the total contractility of arbitrarily sized spheroids. We show that our method is accurate up to strains of 50% and remains valid even for irregularly shaped tissue samples when considering only the deformations in the far field. Finally, we demonstrate that collective forces of tumor spheroids reflect the contractility of individual cells for up to 1 hr after seeding, while collective forces on longer timescales are guided by mechanical feedback from the extracellular matrix.

Modulus of Fibrous Collagen at the Length Scale of a Cell

The extracellular matrix provides macroscale structural support to tissues as well as microscale mechanical cues, like stiffness, to the resident cells. As those cues modulate gene expression, proliferation, differentiation, and motility, quantifying the stiffness that cells sense is crucial to understanding cell behavior. Whereas the macroscopic modulus of a collagen network can be measured in uniform extension or shear, quantifying the local stiffness sensed by a cell remains a challenge due to the inhomogeneous and nonlinear nature of the fiber network at the scale of the cell. To address this challenge, we designed an experimental method to measure the modulus of a network of collagen fibers at this scale. We used spherical particles of an active hydrogel (poly N-isopropylacrylamide) that contract when heated, thereby applying local forces to the collagen matrix and mimicking the contractile forces of a cell. After measuring the particles' bulk modulus and contraction in networks of collagen fibers, we applied a nonlinear model for fibrous materials to compute the modulus of the local region surrounding each particle. We found the modulus at this length scale to be highly heterogeneous, with modulus varying by a factor of 3. In addition, at different values of applied strain, we observed both strain stiffening and strain softening, indicating nonlinearity of the collagen network. Thus, this experimental method quantifies local mechanical properties in a fibrous network at the scale of a cell, while also accounting for inherent nonlinearity.

Force chains in cell-cell mechanical communication

Force chains (FCs) are a key determinant of the micromechanical properties and behaviour of heterogeneous materials, such as granular systems. However, less is known about FCs in fibrous materials, such as the networks composing the extracellular matrix (ECM) of biological systems. Using a finite-element computational model, we simulated the contraction of a single cell and two nearby cells embedded in two-dimensional fibrous elastic networks and analysed the tensile FCs that developed in the ECM. The role of ECM nonlinear elasticity on FC formation was evaluated by considering linear and nonlinear, i.e. exhibiting 'buckling' and/or 'strain-stiffening', stress-strain curves. The effect of the degree of cell contraction and network coordination value was assessed. We found that nonlinear elasticity of the ECM fibres influenced the structure of the FCs, facilitating the transition towards more distinct chains that were less branched and more radially oriented than the chains formed in linear elastic networks. When two neighbouring cells contract, a larger number of FCs bridged between the cells in nonlinear networks, and these chains had a larger effective rigidity than the chains that did not reach a neighbouring cell. These results suggest that FCs function as a route for mechanical communication between distant cells and highlight the contribution of ECM fibre nonlinear elasticity to the formation of FCs.

Heterogeneity and nonaffinity of cell-induced matrix displacements

Cell contractile forces deform and reorganize the surrounding matrix, but the relationship between the forces and the resulting displacements is complicated by the fact that the fibrous structure brings about a complex set of mechanical properties. Many studies have quantified nonlinear and time-dependent properties at macroscopic scales, but it is unclear whether macroscopic properties apply to the scale of a cell, where the matrix is composed of a heterogeneous network of fibers. To address this question, we mimicked the contraction of a cell embedded within a fibrous collagen matrix and quantified the resulting displacements. The data revealed displacements that were heterogeneous and nonaffine. The heterogeneity was reproducible during cyclic loading, and it decreased with decreasing fiber length. Both the experiments and a fiber network model showed that the heterogeneous displacements decayed over distance at a rate no faster than the average displacement field, indicating no transition to homogeneous continuum behavior. Experiments with cells fully embedded in collagen matrices revealed the presence of heterogeneous displacements as well, exposing the dramatic heterogeneity in matrix reorganization that is induced by cells at different positions within the same fibrous matrix.

Nonlinear Elasticity of the ECM Fibers Facilitates Efficient Intercellular Communication

Biological cells embedded in fibrous matrices have been observed to form intercellular bands of dense and aligned fibers through which they mechanically interact over long distances. Such matrix-mediated cellular interactions have been shown to regulate various biological processes. This study aimed to explore the effects of elastic nonlinearity of the fibers contained in the extracellular matrix (ECM) on the transmission of mechanical loads between contracting cells. Based on our biological experiments, we developed a finite-element model of two contracting cells embedded within a fibrous network. The individual fibers were modeled as showing linear elasticity, compression microbuckling, tension stiffening, or both of the latter two. Fiber compression buckling resulted in smaller loads in the ECM, which were primarily directed toward the neighboring cell. These loads decreased with increasing cell-to-cell distance; when cells were >9 cell diameters apart, no such intercellular interaction was observed. Tension stiffening further contributed to directing the loads toward the neighboring cell, though to a smaller extent. The contraction of two neighboring cells resulted in mutual attraction forces, which were considerably increased by tension stiffening and decayed with increasing cell-to-cell distances. Nonlinear elasticity contributed also to the onset of force polarity on the cell boundaries, manifested by larger contractile forces pointing toward the neighboring cell. The density and alignment of the fibers within the intercellular band were greater when fibers buckled under compression, with tension stiffening further contributing to this structural remodeling. Although previous studies have established the role of the ECM nonlinear mechanical behavior in increasing the range of force transmission, our model demonstrates the contribution of nonlinear elasticity of biological gels to directional and efficient mechanical signal transfer between distant cells, and rehighlights the importance of using fibrous gels in experimental settings for facilitating intercellular communication. VIDEO ABSTRACT.