Simulated remodeling of loaded collagen networks via strain-dependent enzymatic degradation and constant-rate fiber growth

Mechanical modeling of the maturation process for tissue-engineered implants: Application to biohybrid heart valves

The development of tissue-engineered cardiovascular implants can improve the lives of large segments of our society who suffer from cardiovascular diseases. Regenerative tissues are fabricated using a process called tissue maturation. Furthermore, it is highly challenging to produce cardiovascular regenerative implants with sufficient mechanical strength to withstand the loading conditions within the human body. Therefore, biohybrid implants for which the regenerative tissue is reinforced by standard reinforcement material (e.g. textile or 3d printed scaffold) can be an interesting solution. In silico models can significantly contribute to characterizing, designing, and optimizing biohybrid implants. The first step towards this goal is to develop a computational model for the maturation process of tissue-engineered implants. This paper focuses on the mechanical modeling of textile-reinforced tissue-engineered cardiovascular implants. First, an energy-based approach is proposed to compute the collagen evolution during the maturation process. Then, the concept of structural tensors is applied to model the anisotropic behavior of the extracellular matrix and the textile scaffold. Next, the newly developed material model is embedded into a special solid-shell finite element formulation with reduced integration. Finally, our framework is used to compute two structural problems: a pressurized shell construct and a tubular-shaped heart valve. The results show the ability of the model to predict collagen growth in response to the boundary conditions applied during the maturation process. Consequently, the model can predict the implant's mechanical response, such as the deformation and stresses of the implant.

The non-affine fiber network solver: A multiscale fiber network material model for finite-element analysis

Multiscale mechanical models in biomaterials research have largely relied on simplifying the microstructure in order to make large-scale simulations tractable. The microscale simplifications often rely on approximations of the constituent distributions and assumptions on the deformation of the constituents. Of particular interest in biomechanics are fiber embedded materials, where simplified fiber distributions and assumed affinity in the fiber deformation greatly influence the mechanical behavior. The consequences of these assumptions are problematic when dealing with microscale mechanical phenomena such as cellular mechanotransduction in growth and remodeling, and fiber-level failure events during tissue failure. In this work, we propose a technique for coupling non-affine network models to finite element solvers, allowing for simulation of discrete microstructural phenomena within macroscopically complex geometries. The developed plugin is readily available as an open-source library for use with the bio-focused finite element software FEBio, and the description of the implementation allows for the adaptation to other finite element solvers.

Hybrid Discrete-Continuum Multiscale Model of Tissue Growth and Remodeling

Tissue growth and remodeling (G&R) is often central to disease etiology and progression, so understanding G&R is essential for understanding disease and developing effective therapies. While the state-of-the-art in this regard is animal and cellular models, recent advances in computational tools offer another avenue to investigate G&R. A major challenge for computational models is bridging from the cellular scale (at which changes are actually occurring) to the macroscopic, geometric-scale (at which physiological consequences arise). Thus, many computational models simplify one scale or another in the name of computational tractability. In this work, we develop a discrete-continuum modeling scheme for analyzing G&R, in which we apply changes directly to the discrete cell and extracellular matrix (ECM) architecture and pass those changes up to a finite-element macroscale geometry. We demonstrate the use of the model in three case-study scenarios: the media of a thick-walled artery, and the media and adventitia of a thick-walled artery, and chronic dissection of an arterial wall. We analyze each case in terms of the new and insightful data that can be gathered from this technique, and we compare our results from this model to several others. STATEMENT OF SIGNIFICANCE: This work is significant in that it provides a framework for combining discrete, microstructural- and cellular-scale models to the growth and remodeling of large tissue structures (such as the aorta). It is a significant advance in that it couples the microscopic remodeling with an existing macroscopic finite element model, making it relatively easy to use for a wide range of conceptual models. It has the potential to improve understanding of many growth and remodeling processes, such as organ formation during development and aneurysm formation, growth, and rupture.

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.

Current Insights on the Diverse Structures and Functions in Bacterial Collagen-like Proteins

The dearth of knowledge on the diverse structures and functions in bacterial collagen-like proteins is in stark contrast to the deep grasp of structures and functions in mammalian collagen, the ubiquitous triple-helical scleroprotein that plays a central role in tissue architecture, extracellular matrix organization, and signal transduction. To fill and highlight existing gaps due to the general paucity of data on bacterial CLPs, we comprehensively reviewed the latest insight into their functional and structural diversity from multiple perspectives of biology, computational simulations, and materials engineering. The origins and discovery of bacterial CLPs were explored. Their genetic distribution and molecular architecture were analyzed, and their structural and functional diversity in various bacterial genera was examined. The principal roles of computational techniques in understanding bacterial CLPs' structural stability, mechanical properties, and biological functions were also considered. This review serves to drive further interest and development of bacterial CLPs, not only for addressing fundamental biological problems in collagen but also for engineering novel biomaterials. Hence, both biology and materials communities will greatly benefit from intensified research into the diverse structures and functions in bacterial collagen-like proteins.

On collagen fiber morphoelasticity and homeostatic remodeling tone

A variety of biochemical and physical processes participate in the creation and maintenance of collagen in biological tissue. Under mechanical stimuli these collagen fibers undergo continuous processes of morphoelastic change. The model presented here is motivated by experimental reports of stretch-stabilization of the collagen fibers to enzymatic degradation. The fiber structure is modeled in terms of a fiber density evolution that is regulated by means of a fixed creation rate and a mechano-sensitive dissolution rate. The theory accounts for the possibly different natural configurations of the fiber unit constituents and the ground substance matrix. It also generalizes previous theoretical descriptions so as to account for finite survival times of the individual fiber units. Special consideration is given to steady state fiber-remodeling processes in which fiber creation and dissolution are in balance. Fiber assembly processes that involve prestretching the fiber constituents yield a homeostatic stress response with a characteristic fiber tone. Fiber density returns to homeostasis after mechanical disruption when sufficient time has passed.

Stress-Swelling Finite Element Modeling of Cervical Response With Homeostatic Collagen Fiber Distributions

During pregnancy, the cervix experiences significant mechanical property change due to tissue swelling, and to ongoing changes in the collagen content. In this paper, we model how these two effects contribute to cervical deformation as the pressure load on top of the cervix increases. The cervix and its surrounding supporting ligaments are taken into consideration in the resulting mechanical analysis. The cervix itself is treated as a multilayered tube-like structure, with layer-specific collagen orientation. The cervical tissue in each layer is treated in terms of a collagen constituent that remodels with time within a ground substance matrix that experiences swelling. The load and swelling are taken to change sufficiently slowly so that the collagen properties at any instant can be regarded as being in a state of homeostasis. Among other things, the simulations show how the luminal cross-sectional area varies along its length as a function of pressure and swelling. In general, an increase in pressure causes an overall shortening of the lumen while an increase in swelling has the opposite effect.

The Balance between Actomyosin Contractility and Microtubule Polymerization Regulates Hierarchical Protrusions That Govern Efficient Fibroblast-Collagen Interactions

Fibroblasts undergo a critical transformation from an initially inactive state to a morphologically different and contractile state after several hours of being embedded within a physiologically relevant three-dimensional (3D) fibrous collagen-based extracellular matrix (ECM). However, little is known about the critical mechanisms by which fibroblasts adapt themselves and their microenvironment in the earliest stage of cell-matrix interaction. Here, we identified the mechanisms by which fibroblasts interact with their 3D collagen fibrous matrices in the early stages of cell-matrix interaction and showed that fibroblasts use energetically efficient hierarchical micro/nano-scaled protrusions in these stages as the primary means for the transformation and adaptation. We found that actomyosin contractility in these protrusions in the early stages of cell-matrix interaction restricts the growth of microtubules by applying compressive forces on them. Our results show that actomyosin contractility and microtubules work in concert in the early stages of cell-matrix interaction to adapt fibroblasts and their microenvironment to one another. These early stage interactions result in responses to disruption of the microtubule network and/or actomyosin contractility that are opposite to well-known responses to late-stage disruption and reveal insight into the ways that cells adapt themselves and their ECM recursively.

A Mesh-Free Approach to Incorporate Complex Anisotropic and Heterogeneous Material Properties into Eye-Specific Finite Element Models

Commercial finite element modeling packages do not have the tools necessary to effectively incorporate the complex anisotropic and heterogeneous material properties typical of the biological tissues of the eye. We propose a mesh-free approach to incorporate realistic material properties into finite element models of individual human eyes. The method is based on the idea that material parameters can be estimated or measured at so called control points, which are arbitrary and independent of the finite element mesh. The mesh-free approach approximates the heterogeneous material parameters at the Gauss points of each finite element while the boundary value problem is solved using the standard finite element method. The proposed method was applied to an eye-specific model a human posterior pole and optic nerve head. We demonstrate that the method can be used to effectively incorporate experimental measurements of the lamina cribrosa micro-structure into the eye-specific model. It was convenient to define characteristic material orientations at the anterior and posterior scleral surface based on the eye-specific geometry of each sclera. The mesh-free approach was effective in approximating these characteristic material directions with smooth transitions across the sclera. For the first time, the method enabled the incorporation of the complex collagen architecture of the peripapillary sclera into an eye-specific model including the recently discovered meridional fibers at the anterior surface and the depth dependent width of circumferential fibers around the scleral canal. The model results suggest that disregarding the meridional fiber region may lead to an underestimation of local strain concentrations in the retina. The proposed approach should simplify future studies that aim to investigate collagen remodeling in the sclera and optic nerve head or in other biological tissues with similar challenges.

The Relationship Between Scleral Strain Change and Differential Cumulative Intraocular Pressure Exposure in the Nonhuman Primate Chronic Ocular Hypertension Model

Purpose

To determine the relationship between peripapillary scleral strain change and cumulative differential IOP exposure in nonhuman primates (NHPs) with unilateral chronic ocular hypertension.

Methods

Posterior scleral shells from 6 bilaterally normal and 10 unilateral chronic ocular hypertension NHPs were pressurized from 5 to 45 mm Hg, and the resulting full-field, three-dimensional, scleral surface deformations were acquired using laser speckle interferometry. Scleral tensile strain (local tissue deformation) was calculated by analytical differentiation of the displacement field; zero strain was assumed at 5 mm Hg. Maximum principal strain was used to represent the scleral strain, and strains were averaged over a 15°-wide (∼3.6-mm) circumpapillary region adjacent to the ONH. The relative difference in mean strain was calculated between fellow eyes and compared with the differential cumulative IOP exposure within NHPs during the study period. The relationship between the relative difference in scleral strain and the differential cumulative IOP exposure in fellow eyes was assessed using an F test and quadratic regression model.

Results

Relative differential scleral tensile strain was significantly associated with differential cumulative IOP exposure in contralateral eyes in the chronic ocular hypertension NHPs, with the bilaterally normal NHPs showing no significant strain difference between fellow eyes. The sclera in the chronic ocular hypertension eyes was more compliant than in their fellow eyes at low levels of differential cumulative IOP exposure, but stiffer at larger differential IOPs (P < 0.0001).

Conclusions

These cross-sectional findings suggest that longitudinal IOP-induced changes in scleral mechanical behavior are dependent on the magnitude of differential cumulative IOP exposure.

Concentration-Dependent Effects of Fibroblast-Like Synoviocytes on Collagen Gel Multiscale Biomechanics and Neuronal Signaling: Implications for Modeling Human Ligamentous Tissues

Abnormal loading of a joint's ligamentous capsule causes pain by activating the capsule's nociceptive afferent fibers, which reside in the capsule's collagenous matrix alongside fibroblast-like synoviocytes (FLS) and transmit pain to the dorsal root ganglia (DRG). This study integrated FLS into a DRG-collagen gel model to better mimic the anatomy and physiology of human joint capsules; using this new model, the effect of FLS on multiscale biomechanics and cell physiology under load was investigated. Primary FLS cells were co-cultured with DRGs at low or high concentrations, to simulate variable anatomical FLS densities, and failed in tension. Given their roles in collagen degradation and nociception, matrix-metalloproteinase (MMP-1) and neuronal expression of the neurotransmitter substance P were probed after gel failure. The amount of FLS did not alter (p > 0.3) the gel failure force, displacement, or stiffness. FLS doubled regional strains at both low (p < 0.01) and high (p = 0.01) concentrations. For high FLS, the collagen network showed more reorganization at failure (p < 0.01). Although total MMP-1 and neuronal substance P were the same regardless of FLS concentration before loading, protein expression of both increased after failure, but only in low FLS gels (p ≤ 0.02). The concentration-dependent effect of FLS on microstructure and cellular responses implies that capsule regions with different FLS densities experience variable microenvironments. This study presents a novel DRG-FLS co-culture collagen gel system that provides a platform for investigating the complex biomechanics and physiology of human joint capsules, and is the first relating DRG and FLS interactions between each other and their surrounding collagen network.

A computational multi-scale approach to investigate mechanically-induced changes in tricuspid valve anterior leaflet microstructure

The tricuspid valve is an atrioventricular valve that prevents blood backflow from the right ventricle into the right atrium during ventricular contractions. It is important to study mechanically induced microstructural alterations in the tricuspid valve leaflets, as this aids both in understanding valvular diseases and in the development of new engineered tissue replacements. The structure and composition of the extracellular matrix (ECM) fiber networks are closely tied to an overall biomechanical function of the tricuspid valve. In this study, we conducted experiments and implemented a multiscale modeling approach to predict ECM microstructural changes to tissue-level mechanical responses in a controlled loading environment. In particular, we characterized a sample of a porcine anterior leaflet at a macroscale using a biaxial mechanical testing method. We then generated a three-dimensional finite element model, to which computational representations of corresponding fiber networks were incorporated based on properties of the microstructural architecture obtained from small angle light scattering. Using five different biaxial boundary conditions, we performed iterative simulations to obtain model parameters with an overall R² value of 0.93. We observed that mechanical loading could markedly alter the underlying ECM architecture. For example, a relatively isotropic fiber network (with an anisotropy index value α of 28%) became noticeably more anisotropic (with an α of 40%) when it underwent mechanical loading. We also observed that the mechanical strain was distributed in a different manner at the ECM/fiber level as compared to the tissue level. The approach presented in this study has the potential to be implemented in pathophysiologically altered biomechanical and structural conditions and to bring insights into the mechanobiology of the tricuspid valve. STATEMENT OF SIGNIFICANCE: Quantifying abnormal cellar/ECM-level deformation of tricuspid valve leaflets subjected to a modified loading environment is of great importance, as it is believed to be linked to valvular remodeling responses. For example, developing surgical procedures or engineered tissue replacements that maintain/mimic ECM-level mechanical homeostasis could lead to more durable outcomes. To quantify leaflet deformation, we built a multiscale framework encompassing the contributions of disorganized ECM components and organized fibers, which can predict the behavior of the tricuspid valve leaflets under physiological loading conditions both at the tissue level and at the ECM level. In addition to future in-depth studies of tricuspid valve pathologies, our model can be used to characterize tissues in other valves of the heart.

Predicting and understanding collagen remodeling in human native heart valves during early development

The hemodynamic functionality of heart valves strongly depends on the distribution of collagen fibers, which are their main load-bearing constituents. It is known that collagen networks remodel in response to mechanical stimuli. Yet, the complex interplay between external load and collagen remodeling is poorly understood. In this study, we adopted a computational approach to simulate collagen remodeling occurring in native fetal and pediatric heart valves. The computational model accounted for several biological phenomena: cellular (re)orientation in response to both mechanical stimuli and topographical cues provided by collagen fibers; collagen deposition and traction forces along the main cellular direction; collagen degradation decreasing with stretch; and cell-mediated collagen prestretch. Importantly, the computational results were well in agreement with previous experimental data for all simulated heart valves. Simulations performed by varying some of the computational parameters suggest that cellular forces and (re)orientation in response to mechanical stimuli may be fundamental mechanisms for the emergence of the circumferential collagen alignment usually observed in native heart valves. On the other hand, the tendency of cells to coalign with collagen fibers is essential to maintain and reinforce that circumferential alignment during development. STATEMENT OF SIGNIFICANCE: The hemodynamic functionality of heart valves is strongly influenced by the alignment of load-bearing collagen fibers. Currently, the mechanisms that are responsible for the development of the circumferential collagen alignment in native heart valves are not fully understood. In the present study, cell-mediated remodeling of native human heart valves during early development was computationally simulated to understand the impact of individual mechanisms on collagen alignment. Our simulations successfully predicted the degree of collagen alignment observed in native fetal and pediatric semilunar valves. The computational results suggest that the circumferential collagen alignment arises from cell traction and cellular (re)orientation in response to mechanical stimuli, and with increasing age is reinforced by the tendency of cells to co-align with pre-existing collagen fibers.

Fiber Network Models Predict Enhanced Cell Mechanosensing on Fibrous Gels

The propagation of mechanical signals through nonlinear fibrous tissues is much more extensive than through continuous synthetic hydrogels. Results from recent studies indicate that increased mechanical propagation arises from the fibrous nature of the material rather than the strain-stiffening property. The relative importance of different parameters of the fibrous network structure to this propagation, however, remains unclear. In this work, we directly compared the mechanical response of substrates of varying thickness subjected to a constant cell traction force using either a nonfibrous strain-stiffening continuum-based model or a volume-averaged fiber network model consisting of two different types of fiber network structures: one with low fiber connectivity (growth networks) and one with high fiber connectivity (Delaunay networks). The growth network fiber models predicted a greater propagation of substrate displacements through the model and a greater sensitivity to gel thickness compared to the more connected Delaunay networks and the nonlinear continuum model. Detailed analysis of the results indicates that rotational freedom of the fibers in a network with low fiber connectivity is critically important for enhanced, long-range mechanosensing. Our findings demonstrate the utility of multiscale models in predicting cells mechanosensing on fibrous gels, and they provide a more complete understanding of how cell traction forces propagate through fibrous tissues, which has implications for the design of engineered tissues and the stem cell niche.

Mathematical model of chronic pancreatitis

Chronic pancreatitis (CP) is a progressive inflammatory disease of the pancreas, leading to its fibrotic destruction. There are currently no drugs that can stop or slow the progression of the disease. The etiology of the disease is multifactorial, whereas recurrent attacks of acute pancreatitis are thought to precede the development of CP. A better understanding of the pathology of CP is needed to facilitate improved diagnosis and treatment strategies for this disease. The present paper develops a mathematical model of CP based on a dynamic network that includes macrophages, pancreatic stellate cells, and prominent cytokines that are present at high levels in the CP microenvironment. The model is represented by a system of partial differential equations. The model is used to explore in silico potential drugs that could slow the progression of the disease, for example infliximab (anti-TNF-[Formula: see text]) and tocilizumab or siltuximab (anti-IL-6/IL-6R).

Mathematical modeling of liver fibrosis

Fibrosis is the formation of excessive fibrous connective tissue in an organ or tissue, which occurs in reparative process or in response to inflammation. Fibrotic diseases are characterized by abnormal excessive deposition of fibrous proteins, such as collagen, and the disease is most commonly progressive, leading to organ disfunction and failure. Although fibrotic diseases evolve in a similar way in all organs, differences may occur as a result of structure and function of the specific organ. In liver fibrosis, the gold standard for diagnosis and monitoring the progression of the disease is biopsy, which is invasive and cannot be repeated frequently. For this reason there is currently a great interest in identifying non-invasive biomarkers for liver fibrosis. In this paper, we develop for the first time a mathematical model of liver fibrosis by a system of partial differential equations. We use the model to explore the efficacy of potential and currently used drugs aimed at blocking the progression of liver fibrosis. We also use the model to develop a diagnostic tool based on a combination of two biomarkers.

Tissue Strain Reorganizes Collagen With a Switchlike Response That Regulates Neuronal Extracellular Signal-Regulated Kinase Phosphorylation In Vitro: Implications for Ligamentous Injury and Mechanotransduction

Excessive loading of ligaments can activate the neural afferents that innervate the collagenous tissue, leading to a host of pathologies including pain. An integrated experimental and modeling approach was used to define the responses of neurons and the surrounding collagen fibers to the ligamentous matrix loading and to begin to understand how macroscopic deformation is translated to neuronal loading and signaling. A neuron-collagen construct (NCC) developed to mimic innervation of collagenous tissue underwent tension to strains simulating nonpainful (8%) or painful ligament loading (16%). Both neuronal phosphorylation of extracellular signal-regulated kinase (ERK), which is related to neuroplasticity (R2 ≥ 0.041; p ≤ 0.0171) and neuronal aspect ratio (AR) (R2 ≥ 0.250; p < 0.0001), were significantly correlated with tissue-level strains. As NCC strains increased during a slowly applied loading (1%/s), a "switchlike" fiber realignment response was detected with collagen reorganization occurring only above a transition point of 11.3% strain. A finite-element based discrete fiber network (DFN) model predicted that at bulk strains above the transition point, heterogeneous fiber strains were both tensile and compressive and increased, with strains in some fibers along the loading direction exceeding the applied bulk strain. The transition point identified for changes in collagen fiber realignment was consistent with the measured strain threshold (11.7% with a 95% confidence interval of 10.2-13.4%) for elevating ERK phosphorylation after loading. As with collagen fiber realignment, the greatest degree of neuronal reorientation toward the loading direction was observed at the NCC distraction corresponding to painful loading. Because activation of neuronal ERK occurred only at strains that produced evident collagen fiber realignment, findings suggest that tissue strain-induced changes in the micromechanical environment, especially altered local collagen fiber kinematics, may be associated with mechanotransduction signaling in neurons.

Micromechanical Modeling Study of Mechanical Inhibition of Enzymatic Degradation of Collagen Tissues

This study investigates how the collagen fiber structure influences the enzymatic degradation of collagen tissues. We developed a micromechanical model of a fibrous collagen tissue undergoing enzymatic degradation based on two central hypotheses. The collagen fibers are crimped in the undeformed configuration. Enzymatic degradation is an energy activated process and the activation energy is increased by the axial strain energy density of the fiber. We determined the intrinsic degradation rate and characteristic energy for mechanical inhibition from fibril-level degradation experiments and applied the parameters to predict the effect of the crimped fiber structure and fiber properties on the degradation of bovine cornea and pericardium tissues under controlled tension. We then applied the model to examine the effect of the tissue stress state on the rate of tissue degradation and the anisotropic fiber structures that developed from enzymatic degradation.

Augmenting Surgery via Multi-scale Modeling and Translational Systems Biology in the Era of Precision Medicine: A Multidisciplinary Perspective

In this era of tremendous technological capabilities and increased focus on improving clinical outcomes, decreasing costs, and increasing precision, there is a need for a more quantitative approach to the field of surgery. Multiscale computational modeling has the potential to bridge the gap to the emerging paradigms of Precision Medicine and Translational Systems Biology, in which quantitative metrics and data guide patient care through improved stratification, diagnosis, and therapy. Achievements by multiple groups have demonstrated the potential for (1) multiscale computational modeling, at a biological level, of diseases treated with surgery and the surgical procedure process at the level of the individual and the population; along with (2) patient-specific, computationally-enabled surgical planning, delivery, and guidance and robotically-augmented manipulation. In this perspective article, we discuss these concepts, and cite emerging examples from the fields of trauma, wound healing, and cardiac surgery.

Cell-matrix interaction during strain-dependent remodelling of simulated collagen networks

The importance of tissue remodelling is widely accepted, but the mechanism by which the remodelling process occurs remains poorly understood. At the tissue scale, the concept of tensional homeostasis, in which there exists a target stress for a cell and remodelling functions to move the cell stress towards that target, is an important foundation for much theoretical work. We present here a theoretical model of a cell in parallel with a network to study what factors of the remodelling process help the cell move towards mechanical stability. The cell-network system was deformed and kept at constant stress. Remodelling was modelled by simulating strain-dependent degradation of collagen fibres and four different cases of collagen addition. The model did not lead to complete tensional homeostasis in the range of conditions studied, but it showed how different expressions for deposition and removal of collagen in a fibre network can interact to modulate the cell's ability to shield itself from an imposed stress by remodelling the surroundings. This study also showed how delicate the balance between deposition and removal rates is and how sensitive the remodelling process is to small changes in the remodelling rules.

A Combined In Vitro Imaging and Multi-Scale Modeling System for Studying the Role of Cell Matrix Interactions in Cutaneous Wound Healing

Many cell types remodel the extracellular matrix of the tissues they inhabit in response to a wide range of environmental stimuli, including mechanical cues. Such is the case in dermal wound healing, where fibroblast migrate into and remodel the provisional fibrin matrix in a complex manner that depends in part on the local mechanical environment and the evolving multi-scale mechanical interactions of the system. In this study, we report on the development of an image-based multi-scale mechanical model that predicts the short-term (24 hours), structural reorganization of a fibrin gel by fibroblasts. These predictive models are based on an in vitro experimental system where clusters of fibroblasts (i.e., explants) were spatially arranged into a triangular geometry onto the surface of fibrin gels that were subjected to either Fixed or Free in-plane mechanical constraints. Experimentally, regional differences in short-term structural remodeling and cell migration were observed for the two gel boundary conditions. A pilot experiment indicated that these small differences in the short-term remodeling of the fibrin gel translate into substantial differences in long-term (4 weeks) remodeling, particularly in terms of collagen production. The multi-scale models were able to predict some regional differences in remodeling and qualitatively similar reorganization patterns for the two boundary conditions. However, other aspects of the model, such as the magnitudes and rates of deformation of gel, did not match the experiments. These discrepancies between model and experiment provide fertile ground for challenging model assumptions and devising new experiments to enhance our understanding of how this multi-scale system functions. These efforts will ultimately improve the predictions of the remodeling process, particularly as it relates to dermal wound healing and the reduction of patient scarring. Such models could be used to recommend patient-specific mechanical-based treatment dependent on parameters such as wound geometry, location, age, and health.

Remodeling of fibrous extracellular matrices by contractile cells: predictions from discrete fiber network simulations

Contractile forces exerted on the surrounding extracellular matrix (ECM) lead to the alignment and stretching of constituent fibers within the vicinity of cells. As a consequence, the matrix reorganizes to form thick bundles of aligned fibers that enable force transmission over distances larger than the size of the cells. Contractile force-mediated remodeling of ECM fibers has bearing on a number of physiologic and pathophysiologic phenomena. In this work, we present a computational model to capture cell-mediated remodeling within fibrous matrices using finite element-based discrete fiber network simulations. The model is shown to accurately capture collagen alignment, heterogeneous deformations, and long-range force transmission observed experimentally. The zone of mechanical influence surrounding a single contractile cell and the interaction between two cells are predicted from the strain-induced alignment of fibers. Through parametric studies, the effect of cell contractility and cell shape anisotropy on matrix remodeling and force transmission are quantified and summarized in a phase diagram. For highly contractile and elongated cells, we find a sensing distance that is ten times the cell size, in agreement with experimental observations.

Human thoracic and abdominal aortic aneurysmal tissues: Damage experiments, statistical analysis and constitutive modeling

Development of aortic aneurysms includes significant morphological changes within the tissue: collagen content increases, elastin content reduces and smooth muscle cells degenerate. We seek to quantify the impact of these changes on the passive mechanical response of aneurysms in the supra-physiological loading range via mechanical testing and constitutive modeling. We perform uniaxial extension tests on circumferentially and axially oriented strips from five thoracic (65.6 years ± 13.4, mean ± SD) and eight abdominal (63.9 years ± 11.4) aortic fusiform aneurysms to investigate both continuous and discontinuous softening during supra-physiological loading. We determine the significance of the differences between the fitted model parameters: diseased thoracic versus abdominal tissues, and healthy (Weisbecker et al., J. Mech. Behav. Biomed. Mater. 12, 93-106, 2012) versus diseased tissues. We also test correlations among these parameters and age, Body Mass Index (BMI) and preoperative aneurysm diameter, and investigate histological cuts. Tissue response is anisotropic for all tests and the anisotropic pseudo-elastic damage model fits the data well for both primary loading and discontinuous softening which we interpret as damage. We found statistically relevant differences between model parameters fitted to diseased thoracic versus abdominal tissues, as well as between those fitted to healthy versus diseased tissues. Only BMI correlated with fitted model parameters in abdominal aortic aneurysmal tissues.

Chain registry and load-dependent conformational dynamics of collagen

Degradation of fibrillar collagen is critical for tissue maintenance. Yet, understanding collagen catabolism has been challenging partly due to a lack of atomistic picture for its load-dependent conformational dynamics, as both mechanical load and local unfolding of collagen affect its cleavage by matrix metalloproteinase (MMP). We use molecular dynamics simulation to find the most cleavage-prone arrangement of α chains in a collagen triple helix and find amino acids that modulate stability of the MMP cleavage domain depending on the chain registry within the molecule. The native-like state is mechanically inhomogeneous, where the cleavage site interfaces a stiff region and a locally unfolded and flexible region along the molecule. In contrast, a triple helix made of the stable glycine-proline-hydroxyproline motif is uniformly flexible and is dynamically stabilized by short-lived, low-occupancy hydrogen bonds. These results provide an atomistic basis for the mechanics, conformation, and stability of collagen that affect catabolism.

Scleral permeability varies by mouse strain and is decreased by chronic experimental glaucoma

PURPOSE

To determine differences in scleral permeability, as measured by diffusion of macromolecules, by using fluorescence recovery after photobleaching (FRAP), with reference to differences by mouse strain, scleral region, and the effect of experimental glaucoma.

METHODS

In three mouse strains (B6, CD1, and B6 mice with mutation in collagen 8α2 [Aca23]), we used FRAP to measure the diffusion of fluorescein isothiocyanate-dextran, molecular weight 40 kDa, into a photobleached zone of sclera. Scleral regions near the optic nerve head (peripapillary) and two successively more anterior regions were compared. Sclera from mouse eyes subjected to chronically elevated intraocular pressure after bead injection into the anterior chamber were compared to fellow eye controls. FRAP data were compared against estimated retinal ganglion cell axon loss in glaucomatous eyes.

RESULTS

Diffusion rates of dextran molecules in the sclera were significantly greater in Aca23 and B6 mice than in CD1 mice in a multivariate model adjusted for region and axial length (P < 0.0001). Dextran diffusion significantly decreased in glaucomatous eyes, and the decline increased with greater axon loss (P = 0.0003, multivariable model). Peripapillary scleral permeability was higher in CD1 than B6 and Aca23 mice (P < 0.05, multivariable model, adjusted by Bonferroni).

CONCLUSIONS

Measurement of the diffusion rates of dextran molecules in the sclera showed that glaucoma leads to decreased scleral permeability in all three mouse strains tested. Among mouse strains tested, those that were more susceptible to glaucomatous loss of retinal ganglion cells had a lower scleral permeability at baseline.

A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load

A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.

Multiscale mechanical simulations of cell compacted collagen gels

Engineered tissues are commonly stretched or compressed (i.e., conditioned) during culture to stimulate extracellular matrix (ECM) production and to improve the mechanical properties of the growing construct. The relationships between mechanical stimulation and ECM remodeling, however, are complex, interdependent, and dynamic. Thus, theoretical models are required for understanding the underlying phenomena so that the conditioning process can be optimized to produce functional engineered tissues. Here, we continue our development of multiscale mechanical models by simulating the effect of cell tractions on developing isometric tension and redistributing forces in the surrounding fibers of a collagen gel embedded with explants. The model predicted patterns of fiber reorganization that were similar to those observed experimentally. Furthermore, the inclusion of cell compaction also changed the distribution of fiber strains in the gel compared to the acellular case, particularly in the regions around the cells where the highest strains were found.

Microscale fiber network alignment affects macroscale failure behavior in simulated collagen tissue analogs

A tissue's microstructure determines its failure properties at larger length scales, however, the specific relationship between microstructure and macroscopic failure in native and engineered soft tissues (such as capsular ligaments, aortic aneurysms, or vascular grafts) has proven elusive. In this study, variations in the microscale fiber alignment in collagen gel tissue analogs were modeled in order to understand their effects on macroscale damage and failure outcomes. The study employed a multiscale finite-element (FE) model for damage and failure in collagen-based materials. The model relied on microstructural representative volume elements (RVEs) that consisted of stochastically-generated networks of discrete type-I collagen fibers. Fiber alignment was varied within RVEs and between layers of RVEs in a macroscopic FE model of a notched dogbone geometry. The macroscale stretch and the microscale response of fibers for each of the differently aligned cases were compared as the dogbone was uniaxially extended to failure. Networks with greater fiber alignment parallel to the direction of extension failed at smaller strains (with a 6-22% reduction in the Green strain at failure), however, at greater grip forces (a 28-60% increase) than networks with fibers aligned perpendicular to the extension. Alternating layers of crisscrossed network alignments (aligned ±45 deg to the direction of extension) failed at smaller strains but at greater grip forces than those created using one fiber alignment type. In summary, variations in microscale structure via fiber alignment produced different macroscale failure trends. To conclude, these findings may be significant in the realm of tissue engineering and in soft tissue biomechanics.

A model of strain-dependent glomerular basement membrane maintenance and its potential ramifications in health and disease

A model is developed and analyzed for type IV collagen turnover in the kidney glomerular basement membrane (GBM), which is the primary structural element in the glomerular capillary wall. The model incorporates strain dependence in both deposition and removal of the GBM, leading to an equilibrium tissue strain at which deposition and removal are balanced. The GBM thickening decreases tissue strain per unit of transcapillary pressure drop according to the law of Laplace, but increases the transcapillary pressure drop required to maintain glomerular filtration. The model results are in agreement with the observed GBM alterations in Alport syndrome and thin basement membrane disease, and the model-predicted linear relation between the inverse capillary radius and inverse capillary thickness at equilibrium is consistent with published data on different mammals. In addition, the model predicts a minimum achievable strain in the GBM based on the geometry, properties, and mechanical environment; that is, an infinitely thick GBM would still experience a finite strain. Although the model assumptions would be invalid for an extremely thick GBM, the minimum achievable strain could be significant in diseases, such as Alport syndrome, characterized by focal GBM thickening. Finally, an examination of reasonable values for the model parameters suggests that the oncotic pressure drop-the osmotic pressure difference between the plasma and the filtrate due to large molecules-plays an important role in setting the GBM strain and, thus, leakage of protein into the urine may be protective against some GBM damage.

Conformational dynamics accompanying the proteolytic degradation of trimeric collagen I by collagenases

Collagenases are the principal enzymes responsible for the degradation of collagens during embryonic development, wound healing, and cancer metastasis. However, the mechanism by which these enzymes disrupt the highly chemically and structurally stable collagen triple helix remains incompletely understood. We used a single-molecule magnetic tweezers assay to characterize the cleavage of heterotrimeric collagen I by both the human collagenase matrix metalloproteinase-1 (MMP-1) and collagenase from Clostridium histolyticum. We observe that the application of 16 pN of force causes an 8-fold increase in collagen proteolysis rates by MMP-1 but does not affect cleavage rates by Clostridium collagenase. Quantitative analysis of these data allows us to infer the structural changes in collagen associated with proteolytic cleavage by both enzymes. Our data support a model in which MMP-1 cuts a transient, stretched conformation of its recognition site. In contrast, our findings suggest that Clostridium collagenase is able to cleave the fully wound collagen triple helix, accounting for its lack of force sensitivity and low sequence specificity. We observe that the cleavage of heterotrimeric collagen is less force sensitive than the proteolysis of a homotrimeric collagen model peptide, consistent with studies suggesting that the MMP-1 recognition site in heterotrimeric collagen I is partially unwound at equilibrium.