Crack Propagation Versus Fiber Alignment in Collagen Gels: Experiments and Multiscale Simulation

Collagen gel contraction assays: From modelling wound healing to quantifying cellular interactions with three-dimensional extracellular matrices

Cells respond to and actively remodel the extracellular matrix (ECM). The dynamic and bidirectional interaction between cells and ECM, especially their mechanical interactions, has been found to play an essential role in triggering a series of complex biochemical and biomechanical signal pathways and in regulating cellular functions and behaviours. The collagen gel contraction assay (CGCA) is a widely used method to investigate cell-ECM interactions in 3D environments and provides a mechanically associated readout reflecting 3D cellular contractility. In this review, we summarize various versions of CGCA, with an emphasis on recent high-throughput and low-consumption CGCA techniques. More importantly, we focus on the technique of force monitoring during the contraction of collagen gel, which provides a quantitative characterization of the overall forces generated by all the resident cells in the collagen hydrogel. Accordingly, we present recent biological applications of the CGCA, which have expanded from the initial wound healing model to other studies concerning cell-ECM interactions, including fibrosis, cancer, tissue repair and the preparation of biomimetic microtissues.

A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue

The heterogeneous, nonlinear, anisotropic material behavior of biological tissues makes precise definition of an accurate constitutive model difficult. One possible solution to this issue would be to define microstructural elements and perform fully coupled multiscale simulation. However, for complex geometries and loading scenarios, the computational costs of such simulations can be prohibitive. Ideally then, we should seek a method that contains microstructural detail, but leverages the speed of classical continuum-based finite-element (FE) modeling. In this work, we demonstrate the use of the Holzapfel-Gasser-Ogden (HGO) model [1, 2] to fit the behavior of microstructural network models. We show that Delaunay microstructural networks can be fit to the HGO strain energy function by calculating fiber network strain energy and average fiber stretch ratio. We then use the HGO constitutive model in a FE framework to improve the speed of our hybrid model, and demonstrate that this method, combined with a material property update scheme, can match a full multiscale simulation. This method gives us flexibility in defining complex FE simulations that would be impossible, or at least prohibitively time consuming, in multiscale simulation, while still accounting for microstructural heterogeneity.

Connectivity and plasticity determine collagen network fracture

Collagen forms the structural scaffold of connective tissues in all mammals. Tissues are remarkably resistant against mechanical deformations because collagen molecules hierarchically self-assemble in fibrous networks that stiffen with increasing strain. Nevertheless, collagen networks do fracture when tissues are overloaded or subject to pathological conditions such as aneurysms. Prior studies of the role of collagen in tissue fracture have mainly focused on tendons, which contain highly aligned bundles of collagen. By contrast, little is known about fracture of the orientationally more disordered collagen networks present in many other tissues such as skin and cartilage. Here, we combine shear rheology of reconstituted collagen networks with computer simulations to investigate the primary determinants of fracture in disordered collagen networks. We show that the fracture strain is controlled by the coordination number of the network junctions, with less connected networks fracturing at larger strains. The hierarchical structure of collagen fine-tunes the fracture strain by providing structural plasticity at the network and fiber level. Our findings imply that low connectivity and plasticity provide protective mechanisms against network fracture that can optimize the strength of biological tissues.

Mechanics of a two-fiber model with one nested fiber network, as applied to the collagen-fibrin system

The mechanical behavior of collagen-fibrin (col-fib) co-gels is both scientifically interesting and clinically relevant. Collagen-fibrin networks are a staple of tissue engineering research, but the mechanical consequences of changes in co-gel composition have remained difficult to predict or even explain. We previously observed fundamental differences in failure behavior between collagen-rich and fibrin-rich co-gels, suggesting an essential change in how the two components interact as the co-gel's composition changes. In this work, we explored the hypothesis that the co-gel behavior is due to a lack of percolation by the dilute component. We generated a series of computational models based on interpenetrating fiber networks. In these models, the major network component percolated the model space but the minor component did not, instead occupying a small island embedded within the larger network. Each component was assigned properties based on a fit of single-component gel data. Island size was varied to match the relative concentrations of the two components. The model predicted that networks rich in collagen, the stiffer component, would roughly match pure-collagen gel behavior with little additional stress due to the fibrin, as seen experimentally. For fibrin-rich gels, however, the model predicted a smooth increase in the overall network strength with added collagen, as seen experimentally but not consistent with an additive parallel model. We thus conclude that incomplete percolation by the low-concentration component of a co-gel is a major determinant of its macroscopic properties, especially if the low-concentration component is the stiffer component.

STATEMENT OF SIGNIFICANCE

Models for the behavior of fibrous networks have useful applications in many different fields, including polymer science, textiles, and tissue engineering. In addition to being important structural components in soft tissues and blood clots, these protein networks can serve as scaffolds for bioartificial tissues. Thus, their mechanical behavior, especially in co-gels, is both interesting from a materials science standpoint and significant with regard to tissue engineering.

Multiscale mechanics of the cervical facet capsular ligament, with particular emphasis on anomalous fiber realignment prior to tissue failure

The facet capsular ligaments encapsulate the bilateral spinal facet joints and are common sources of painful injury due to afferent innervation. These ligaments exhibit architectural complexity, which is suspected to contribute to the experimentally observed lack of co-localization between macroscopic strain and microstructural tissue damage. The heterogeneous and multiscale nature of this ligament, combined with challenges in experimentally measuring its microscale mechanics, hinders the ability to understand sensory mechanisms under normal or injurious loading. Therefore, image-based, subject-specific, multiscale finite-element models were constructed to predict the mechanical responses of the human cervical facet capsular ligament under uniaxial tensile stretch. The models precisely simulated the force-displacement responses for all samples ([Formula: see text]) and showed promise in predicting the magnitude and location of peak regional strains at two different displacements. Yet, there was a loss of agreement between the model and experiment in terms of fiber organization at large tissue stretch, possibly due to a lack of accounting for tissue failure. The mean fiber stretch ratio predicted by the models was found to be significantly higher in regions that exhibited anomalous fiber realignment experimentally than in regions with normal realignment ([Formula: see text]). The development of microstructural abnormalities was associated with the predicted fiber-level stretch ([Formula: see text]), but not with the elemental maximum principal stress or maximum principal strain by logistic regression. The multiscale models elucidate a potential mechanical basis for predicting injury-prone tissue domains and for defining the relationships between macroscopic ligament stretch and microscale pathophysiology in the subfailure regime.