Nonlinear Mechanics of Athermal Branched Biopolymer Networks

Effects of local incompressibility on the rheology of composite biopolymer networks

Fibrous networks such as collagen are common in biological systems. Recent theoretical and experimental efforts have shed light on the mechanics of single component networks. Most real biopolymer networks, however, are composites made of elements with different rigidity. For instance, the extracellular matrix in mammalian tissues consists of stiff collagen fibers in a background matrix of flexible polymers such as hyaluronic acid (HA). The interplay between different biopolymer components in such composite networks remains unclear. In this work, we use 2D coarse-grained models to study the nonlinear strain-stiffening behavior of composites. We introduce a local volume constraint to model the incompressibility of HA. We also perform rheology experiments on composites of collagen with HA. Theoretically and experimentally, we demonstrate that the linear shear modulus of composite networks can be increased by approximately an order of magnitude above the corresponding moduli of the pure components. Our model shows that this synergistic effect can be understood in terms of the local incompressibility of HA, which acts to suppress density fluctuations of the collagen matrix with which it is entangled.

Effects of coordination and stiffness scale separation in disordered elastic networks

Many fibrous materials are modeled as elastic networks featuring a substantial separation between the stiffness scales that characterize different microscopic deformation modes of the network's constituents. This scale separation has been shown to give rise to emergent complexity in these systems' linear and nonlinear mechanical response. Here we study numerically a simple model featuring said stiffness scale separation in two-dimensions and show that its mechanical response is governed by the competition between the characteristic stiffness of collective nonphononic soft modes of the stiff subsystem, and the characteristic stiffness of the soft interactions. We present and rationalize the behavior of the shear modulus of our complex networks across the unjamming transition at which the stiff subsystem alone loses its macroscopic mechanical rigidity. We further establish a relation in the soft-interaction-dominated regime between the shear modulus, the characteristic frequency of nonphononic vibrational modes, and the mesoscopic correlation length that marks the crossover from a disorder-dominated response to local mechanical perturbations in the near field, to a linear, continuumlike response in the far field. The effects of spatial dimension on the observed scaling behavior are discussed, in addition to the interplay between stiffness scales in strain-stiffened networks, which is relevant to understanding the nonlinear mechanics of non-Brownian fibrous biomatter.

Review paper: The importance of consideration of collagen cross-links in computational models of collagen-based tissues

Collagen as the main protein in Extra Cellular Matrix (ECM) is the main load-bearing component of fibrous tissues. Nanostructure and architecture of collagen fibrils play an important role in mechanical behavior of these tissues. Extensive experimental and theoretical studies have so far been performed to capture these properties, but none of the current models realistically represent the complexity of network mechanics because still less is known about the collagen's inner structure and its effect on the mechanical properties of tissues. The goal of this review article is to emphasize the significance of cross-links in computational modeling of different collagen-based tissues, and to reveal the need for continuum models to consider cross-links properties to better reflect the mechanical behavior observed in experiments. In addition, this study outlines the limitations of current investigations and provides potential suggestions for the future work.

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.

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.

Force Transmission in Disordered Fibre Networks

Cells residing in living tissues apply forces to their immediate surroundings to promote the restructuration of the extracellular matrix fibres and to transmit mechanical signals to other cells. Here we use a minimalist model to study how these forces, applied locally by cell contraction, propagate through the fibrous network in the extracellular matrix. In particular, we characterize how the transmission of forces is influenced by the connectivity of the network and by the bending rigidity of the fibers. For highly connected fiber networks the stresses spread out isotropically around the cell over a distance that first increases with increasing contraction of the cell and then saturates at a characteristic length. For lower connectivity, however, the stress pattern is highly asymmetric and is characterised by force chains that can transmit stresses over very long distances. We hope that our analysis of force transmission in fibrous networks can provide a new avenue for future studies on how the mechanical feedback between the cell and the ECM is coupled with the microscopic environment around the cells.

Light-sheet photonic force optical coherence elastography for high-throughput quantitative 3D micromechanical imaging

Quantitative characterisation of micro-scale mechanical properties of the extracellular matrix (ECM) and dynamic cell-ECM interactions can significantly enhance fundamental discoveries and their translational potential in the rapidly growing field of mechanobiology. However, quantitative 3D imaging of ECM mechanics with cellular-scale resolution and dynamic monitoring of cell-mediated changes to pericellular viscoelasticity remain a challenge for existing mechanical characterisation methods. Here, we present light-sheet photonic force optical coherence elastography (LS-pfOCE) to address this need by leveraging a light-sheet for parallelised, non-invasive, and localised mechanical loading. We demonstrate the capabilities of LS-pfOCE by imaging the micromechanical heterogeneity of fibrous collagen matrices and perform live-cell imaging of cell-mediated ECM micromechanical dynamics. By providing access to 4D spatiotemporal variations in the micromechanical properties of 3D biopolymer constructs and engineered cellular systems, LS-pfOCE has the potential to drive new discoveries in mechanobiology and contribute to the development of novel biomechanics-based clinical diagnostics and therapies.

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.

Energetic rigidity. I. A unifying theory of mechanical stability

Rigidity regulates the integrity and function of many physical and biological systems. This is the first of two papers on the origin of rigidity, wherein we propose that "energetic rigidity," in which all nontrivial deformations raise the energy of a structure, is a more useful notion of rigidity in practice than two more commonly used rigidity tests: Maxwell-Calladine constraint counting (first-order rigidity) and second-order rigidity. We find that constraint counting robustly predicts energetic rigidity only when the system has no states of self-stress. When the system has states of self-stress, we show that second-order rigidity can imply energetic rigidity in systems that are not considered rigid based on constraint counting, and is even more reliable than shear modulus. We also show that there may be systems for which neither first- nor second-order rigidity imply energetic rigidity. The formalism of energetic rigidity unifies our understanding of mechanical stability and also suggests new avenues for material design.

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.

Jamming and arrest of cell motion in biological tissues

Collective cell motility is crucial to many biological processes including morphogenesis, wound healing, and cancer invasion. Recently, the biology and biophysics communities have begun to use the term 'cell jamming' to describe the collective arrest of cell motion in tissues. Although this term is widely used, the underlying mechanisms are varied. In this review, we highlight three independent mechanisms that can potentially drive arrest of cell motion - crowding, tension-driven rigidity, and reduction of fluctuations - and propose a framework that connects all three. Because multiple mechanisms may be operating simultaneously, this emphasizes that experiments should strive to identify which mechanism dominates in a given situation. We also discuss how specific cell-scale and molecular-scale biological processes, such as cell-cell and cell-substrate interactions, control aspects of these underlying physical mechanisms.

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.

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.

Elastic Anisotropy Governs the Range of Cell-Induced Displacements

The unique nonlinear mechanics of the fibrous extracellular matrix (ECM) facilitates long-range cell-cell mechanical communications that would be impossible for linear elastic substrates. Past research has described the contribution of two separated effects on the range of force transmission, including ECM elastic nonlinearity and fiber alignment. However, the relation between these different effects is unclear, and how they combine to dictate force transmission range is still elusive. Here, we combine discrete fiber simulations with continuum modeling to study the decay of displacements induced by a contractile cell in fibrous networks. We demonstrate that fiber nonlinearity and fiber reorientation both contribute to the strain-induced elastic anisotropy of the cell's local environment. This elastic anisotropy is a "lumped" parameter that governs the slow decay of displacements, and it depends on the magnitude of applied strain, either an external tension or an internal contraction, as a model of the cell. Furthermore, we show that accounting for artificially prescribed elastic anisotropy dictates the decay of displacements induced by a contracting cell. Our findings unify previous single effects into a mechanical theory that explains force transmission in fibrous networks. This work may provide insights into biological processes that involve communication of distant cells mediated by the ECM, such as those occurring in morphogenesis, wound healing, angiogenesis, and cancer metastasis. It may also provide design parameters for biomaterials to control force transmission between cells as a way to guide morphogenesis in tissue engineering.

Stiffness Sensing by Cells

Physical stimuli are essential for the function of eukaryotic cells, and changes in physical signals are important elements in normal tissue development as well as in disease initiation and progression. The complexity of physical stimuli and the cellular signals they initiate are as complex as those triggered by chemical signals. One of the most important, and the focus of this review, is the effect of substrate mechanical properties on cell structure and function. The past decade has produced a nearly exponentially increasing number of mechanobiological studies to define how substrate stiffness alters cell biology using both purified systems and intact tissues. Here we attempt to identify common features of mechanosensing in different systems while also highlighting the numerous informative exceptions to what in early studies appeared to be simple rules by which cells respond to mechanical stresses.

Rigidity and auxeticity transitions in networks with strong bond-bending interactions

A widely studied model for gels or biopolymeric fibrous materials are networks with central force interactions, such as Hookean springs. Less commonly studied are materials whose mechanics are dominated by non-central force interactions such as bond-bending potentials. Inspired by recent experimental advancements in designing colloidal gels with tunable interactions, we study the micro- and macroscopic elasticity of two-dimensional planar graphs with strong bond-bending potentials, in addition to weak central forces. We introduce a theoretical framework that allows us to directly investigate the limit in which the ratio of characteristic central-force to bending stiffnesses vanishes. In this limit we show that a generic isostatic point exists at [Formula: see text], coinciding with the isostatic point of frames with central-force interactions in two dimensions. We further demonstrate the emergence of a stiffening transition when the coordination is increased towards the isostatic point, which shares similarities with the strain-induced stiffening transition observed in biopolymeric fibrous materials, and coincides with an auxeticity transition above which the material's Poisson's ratio approaches -1 when bond-bending interactions dominate.

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.

Scaling Theory for Mechanical Critical Behavior in Fiber Networks

As a function of connectivity, spring networks exhibit a critical transition between floppy and rigid phases at an isostatic threshold. For connectivity below this threshold, fiber networks were recently shown theoretically to exhibit a rigidity transition with corresponding critical signatures as a function of strain. Experimental collagen networks were also shown to be consistent with these predictions. We develop a scaling theory for this strain-controlled transition. Using a real-space renormalization approach, we determine relations between the critical exponents governing the transition, which we verify for the strain-controlled transition using numerical simulations of both triangular lattice-based and packing-derived 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.

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.

Combined experimental and computational characterization of crosslinked collagen-based hydrogels

Collagen hydrogels are widely used for in-vitro experiments and tissue engineering applications. Their use has been extended due to their biocompatibility with cells and their capacity to mimic biological tissues; nevertheless their mechanical properties are not always optimal for these purposes. Hydrogels are formed by a network of polymer filaments embedded on an aqueous substrate and their mechanical properties are mainly defined by the filament network architecture and the individual filament properties. To increase properties of native collagen, such as stiffness or strain-stiffening, these networks can be modified by adding crosslinking agents that alter the network architecture, increasing the unions between filaments. In this work, we have investigated the effect of one crosslinking agent, transglutaminase, in collagen hydrogels with varying collagen concentration. We have observed a linear dependency of the gel rigidity on the collagen concentration. Moreover, the addition of transglutaminase has induced an earlier strain-stiffening of the collagen gels. In addition, to better understand the mechanical implications of collagen concentration and crosslinkers inclusion, we have adapted an existing computational model, based on the worm-like chain model (WLC), to reproduce the mechanical behavior of the collagen gels. With this model we can estimate the parameters of the biopolymer networks without more sophisticated techniques, such as image processing or network reconstruction, or, inversely, predict the mechanical properties of a defined collagen network.

Probing Gelation and Rheological Behavior of a Self-Assembled Molecular Gel

Molecular gels have been investigated over the last few decades; however, mechanical behavior of these self-assembled gels is not well understood, particularly how these materials fail at large strain. Here, we report the gelation and rheological behavior of a molecular gel formed by self-assembly of a low molecular weight gelator (LMWG), di-Fmoc-l-lysine, in 1-propanol/water mixture. Gels were prepared by solvent-triggered technique, and gelation was tracked using Fourier transform infrared (FTIR) spectroscopy and shear rheology. FTIR spectroscopy captures the formation of hydrogen bonding between the gelator molecules, and the change in IR spectra during the gelation process correlates with the gelation kinetics results captured by rheology. Self-assembly of gelator molecules leads to a fiber-like structure, and these long fibers topologically interact to form a gel-like material. Stretched-exponential function can capture the stress-relaxation data. Stress-relaxation time for these gels have been found to be long owing to long fiber dimensions, and the stretching exponent value of 1/3 indicates polydispersity in fiber dimensions. Cavitation rheology captures fracture-like behavior of these gels, and critical energy release rate has been estimated to be of the order 0.1 J/m². Our results provide new understanding of the rheological behavior of molecular gels and their structural origin.

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.