An experimental and modeling study of the viscoelastic behavior of collagen gel

Stress relaxation in network materials: the contribution of the network

Stress relaxation in network materials with permanent crosslinks is due to the transport of fluid within the network (poroelasticity), the viscoelasticity of the matrix and the viscoelasticity of the network. While relaxation associated with the matrix was studied extensively, the contribution of the network remains unexplored. In this work we consider two and three-dimensional stochastic fiber networks with viscoelastic fibers and explore the dependence of stress relaxation on network structure. We observe that relaxation has two regimes - an initial exponential regime, followed by a stretched exponential regime - similar to the situation in other disordered materials. The stretch exponent is a function of density, fiber diameter and the network structure, and has a minimum at the transition between the affine and non-affine regimes of network behavior. The relaxation time constant of the first, exponential regime is similar to the relaxation time constant of individual fibers and is independent of network density and fiber diameter. The relaxation time constant of the second, stretched exponential regime is a weak function of network parameters. The stretched exponential emerges from the heterogeneity of relaxation dynamics on scales comparable with the mesh size, with higher heterogeneity leading to smaller stretch exponents. In composite networks of fibers whose relaxation time constant is selected from a distribution with set mean, the stretch exponent decreases with increasing the coefficient of variation of the fiber time constant distribution. As opposed to thermal glass formers and colloids, in these athermal systems the dynamic heterogeneity is introduced by the network structure and does not evolve during relaxation. While in thermal systems the control parameter is the temperature, in this athermal case the control parameter is a non-dimensional structural parameter which describes the degree of non-affinity of the network.

Local tissue heterogeneity may modulate neuronal responses via altered axon strain fields: insights about innervated joint capsules from a computational model

In innervated collagenous tissues, tissue scale loading may contribute to joint pain by transmitting force through collagen fibers to the embedded mechanosensitive axons. However, the highly heterogeneous collagen structures of native tissues make understanding this relationship challenging. Recently, collagen gels with embedded axons were stretched and the resulting axon signals were measured, but these experiments were unable to measure the local axon strain fields. Computational discrete fiber network models can directly determine axon strain fields due to tissue scale loading. Therefore, this study used a discrete fiber network model to identify how heterogeneous collagen networks (networks with multiple collagen fiber densities) change axon strain due to tissue scale loading. In this model, a composite cylinder (axon) was embedded in a Delaunay network (collagen). Homogeneous networks with a single collagen volume fraction and two types of heterogeneous networks with either a sparse center or dense center were created. Measurements of fiber forces show higher magnitude forces in sparse regions of heterogeneous networks and uniform force distributions in homogeneous networks. The average axon strain in the sparse center networks decreases when compared to homogeneous networks with similar collagen volume fractions. In dense center networks, the average axon strain increases compared to homogeneous networks. The top 1% of axon strains are unaffected by network heterogeneity. Based on these results, the interaction of tissue scale loading, collagen network heterogeneity, and axon strains in native musculoskeletal tissues should be considered when investigating the source of joint pain.

A viscoelastic model for human myocardium

Understanding the biomechanics of the heart in health and disease plays an important role in the diagnosis and treatment of heart failure. The use of computational biomechanical models for therapy assessment is paving the way for personalized treatment, and relies on accurate constitutive equations mapping strain to stress. Current state-of-the art constitutive equations account for the nonlinear anisotropic stress-strain response of cardiac muscle using hyperelasticity theory. While providing a solid foundation for understanding the biomechanics of heart tissue, most current laws neglect viscoelastic phenomena observed experimentally. Utilizing experimental data from human myocardium and knowledge of the hierarchical structure of heart muscle, we present a fractional nonlinear anisotropic viscoelastic constitutive model. The model is shown to replicate biaxial stretch, triaxial cyclic shear and triaxial stress relaxation experiments (mean error ∼7.68%), showing improvements compared to its hyperelastic (mean error ∼24%) counterparts. Model sensitivity, fidelity and parameter uniqueness are demonstrated. The model is also compared to rate-dependent biaxial stretch as well as different modes of biaxial stretch, illustrating extensibility of the model to a range of loading phenomena. STATEMENT OF SIGNIFICANCE: The viscoelastic response of human heart tissues has yet to be integrated into common constitutive models describing cardiac mechanics. In this work, a fractional viscoelastic modeling approach is introduced based on the hierarchical structure of heart tissue. From these foundations, the current state-of-the-art biomechanical models of the heart muscle are transformed using fractional viscoelasticity, replicating passive muscle function across multiple experimental tests. Comparisons are drawn with current models to highlight the improvements of this approach and predictive responses show strong qualitative agreement with experimental data. The proposed model presents the first constitutive model aimed at capturing viscoelastic nonlinear response across multiple testing regimes, providing a platform for better understanding the biomechanics of myocardial tissue in health and disease.

Liver Bioreactor Design Issues of Fluid Flow and Zonation, Fibrosis, and Mechanics: A Computational Perspective

Tissue engineering, with the goal of repairing or replacing damaged tissue and organs, has continued to make dramatic science-based advances since its origins in the late 1980’s and early 1990’s. Such advances are always multi-disciplinary in nature, from basic biology and chemistry through physics and mathematics to various engineering and computer fields. This review will focus its attention on two topics critical for tissue engineering liver development: (a) fluid flow, zonation, and drug screening, and (b) biomechanics, tissue stiffness, and fibrosis, all within the context of 3D structures. First, a general overview of various bioreactor designs developed to investigate fluid transport and tissue biomechanics is given. This includes a mention of computational fluid dynamic methods used to optimize and validate these designs. Thereafter, the perspective provided by computer simulations of flow, reactive transport, and biomechanics responses at the scale of the liver lobule and liver tissue is outlined, in addition to how bioreactor-measured properties can be utilized in these models. Here, the fundamental issues of tortuosity and upscaling are highlighted, as well as the role of disease and fibrosis in these issues. Some idealized simulations of the effects of fibrosis on lobule drug transport and mechanics responses are provided to further illustrate these concepts. This review concludes with an outline of some practical applications of tissue engineering advances and how efficient computational upscaling techniques, such as dual continuum modeling, might be used to quantify the transition of bioreactor results to the full liver scale.

Poroelastic Modeling of Highly Hydrated Collagen Hydrogels: Experimental Results vs. Numerical Simulation With Custom and Commercial Finite Element Solvers

This study presents a comparison between the performances of two Finite Element (FE) solvers for the modeling of the poroelastic behavior of highly hydrated collagen hydrogels. Characterization of collagen hydrogels has been a widespread challenge since this is one of the most used natural biomaterials for Tissue Engineering (TE) applications. V-Biomech® is a free custom FE solver oriented to soft tissue modeling, while Abaqus® is a general-purpose commercial FE package which is widely used for biomechanics computational modeling. Poroelastic simulations with both solvers were compared to two experimental protocols performed by Busby et al. (2013) and Chandran and Barocas (2004), also using different implementations of the frequently used Neo-Hookean hyperelastic model. The average differences between solvers outputs were under 5% throughout the different tests and hydrogel properties. Thus, differences were small enough to be considered negligible and within the variability found experimentally from one sample to another. This work demonstrates that constitutive modeling of soft tissues, such as collagen hydrogels can be achieved with either V-Biomech or Abaqus standard options (without user-subroutine), which is important for the biomechanics and biomaterials research community. V-Biomech has shown its potential for the validation of biomechanical characterization of soft tissues, while Abaqus' versatility is useful for the modeling and analysis of TE applications where other complex phenomena may also need to be captured.

Viscoelasticity of spinal cord and meningeal tissues

Compared to the outer dura mater, the mechanical behavior of spinal pia and arachnoid meningeal layers has received very little attention in the literature. This is despite experimental evidence of their importance with respect to the overall spinal cord stiffness and recovery following compression. Accordingly, inclusion of the mechanical contribution of the pia and arachnoid maters would improve the predictive accuracy of finite element models of the spine, especially in the distribution of stresses and strain through the cord's cross-section. However, to-date, only linearly elastic moduli for what has been previously identified as spinal pia mater is available in the literature. This study is the first to quantitatively compare the viscoelastic behavior of isolated spinal pia-arachnoid-complex, neural tissue of the spinal cord parenchyma, and intact construct of the two. The results show that while it only makes up 5.5% of the overall cross-sectional area, the thin membranes of the innermost meninges significantly affect both the elastic and viscous response of the intact construct. Without the contribution of the pia and arachnoid maters, the spinal cord has very little inherent stiffness and experiences significant relaxation when strained. The ability of the fitted non-linear viscoelastic material models of each condition to predict independent data within experimental variability supports their implementation into future finite element computational studies of the spine.

STATEMENT OF SIGNIFICANCE

The neural tissue of the spinal cord is surrounded by three fibrous layers called meninges which are important in the behavior of the overall spinal-cord-meningeal construct. While the mechanical properties of the outermost layer have been reported, the pia mater and arachnoid mater have received considerably less attention. This study is the first to directly compare the behavior of the isolated neural tissue of the cord, the isolated pia-arachnoid complex, and the construct of these individual components. The results show that, despite being very thin, the inner meninges significantly affect the elastic and time-dependent response of the spinal cord, which may have important implications for studies of spinal cord injury.

Micromechanical study of the load transfer in a polycaprolactone-collagen hybrid scaffold when subjected to unconfined and confined compression

Scaffolds are used in diverse tissue engineering applications as hosts for cell proliferation and extracellular matrix formation. One of the most used tissue engineering materials is collagen, which is well known to be a natural biomaterial, also frequently used as cell substrate, given its natural abundance and intrinsic biocompatibility. This study aims to evaluate how the macroscopic biomechanical stimuli applied on a construct made of polycaprolactone scaffold embedded in a collagen substrate translate into microscopic stimuli at the cell level. Eight poro-hyperelastic finite element models of 3D printed hybrid scaffolds from the same batch were created, along with an equivalent model of the idealized geometry of that scaffold. When applying an 8% confined compression at the macroscopic level, local fluid flow of up to 20 [Formula: see text]m/s and octahedral strain levels mostly under 20% were calculated in the collagen substrate. Conversely unconfined compression induced fluid flow of up to 10 [Formula: see text]m/s and octahedral strain from 10 to 35%. No relevant differences were found amongst the scaffold-specific models. Following the mechanoregulation theory based on Prendergast et al. (J Biomech 30:539-548, 1997. https://doi.org/10.1016/S0021-9290(96)00140-6 ), those results suggest that mainly cartilage or fibrous tissue formation would be expected to occur under unconfined or confined compression, respectively. This in silico study helps to quantify the microscopic stimuli that are present within the collagen substrate and that will affect cell response under in vitro bioreactor mechanical stimulation or even after implantation.

Predicting cell viability within tissue scaffolds under equiaxial strain: multi-scale finite element model of collagen-cardiomyocytes constructs

Successful tissue engineering and regenerative therapy necessitate having extensive knowledge about mechanical milieu in engineered tissues and the resident cells. In this study, we have merged two powerful analysis tools, namely finite element analysis and stochastic analysis, to understand the mechanical strain within the tissue scaffold and residing cells and to predict the cell viability upon applying mechanical strains. A continuum-based multi-length scale finite element model (FEM) was created to simulate the physiologically relevant equiaxial strain exposure on cell-embedded tissue scaffold and to calculate strain transferred to the tissue scaffold (macro-scale) and residing cells (micro-scale) upon various equiaxial strains. The data from FEM were used to predict cell viability under various equiaxial strain magnitudes using stochastic damage criterion analysis. The model validation was conducted through mechanically straining the cardiomyocyte-encapsulated collagen constructs using a custom-built mechanical loading platform (EQUicycler). FEM quantified the strain gradients over the radial and longitudinal direction of the scaffolds and the cells residing in different areas of interest. With the use of the experimental viability data, stochastic damage criterion, and the average cellular strains obtained from multi-length scale models, cellular viability was predicted and successfully validated. This methodology can provide a great tool to characterize the mechanical stimulation of bioreactors used in tissue engineering applications in providing quantification of mechanical strain and predicting cellular viability variations due to applied mechanical strain.

Combined numerical and experimental biomechanical characterization of soft collagen hydrogel substrate

This work presents a combined experimental-numerical framework for the biomechanical characterization of highly hydrated collagen hydrogels, namely with 0.20, 0.30 and 0.40% (by weight) of collagen concentration. Collagen is the most abundant protein in the extracellular matrix of animals and humans. Its intrinsic biocompatibility makes collagen a promising substrate for embedding cells within a highly hydrated environment mimicking natural soft tissues. Cell behaviour is greatly influenced by the mechanical properties of the surrounding matrix, but the biomechanical characterization of collagen hydrogels has been challenging up to now, since they present non-linear poro-viscoelastic properties. Combining the stiffness outcomes from rheological experiments with relevant literature data on collagen permeability, poroelastic finite element (FE) models were developed. Comparison between experimental confined compression tests available in the literature and analogous FE stress relaxation curves showed a close agreement throughout the tests. This framework allowed establishing that the dynamic shear modulus of the collagen hydrogels is between 0.0097 ± 0.018 kPa for the 0.20% concentration and 0.0601 ± 0.044 kPa for the 0.40% concentration. The Poisson's ratio values for such conditions lie within the range of 0.495-0.485 for 0.20% and 0.480-0.470 for 0.40%, respectively, showing that rheology is sensitive enough to detect these small changes in collagen concentration and thus allowing to link rheology results with the confined compression tests. In conclusion, this integrated approach allows for accurate constitutive modelling of collagen hydrogels. This framework sets the grounds for the characterization of related hydrogels and to the use of this collagen parameterization in more complex multiscale models.

Three-dimensional morphogenesis of MDCK cells induced by cellular contractile forces on a viscous substrate

Substrate physical properties are essential for many physiological events such as embryonic development and 3D tissue formation. Physical properties of the extracellular matrix such as viscoelasticity and geometrical constraints are understood as factors that affect cell behaviour. In this study, we focused on the relationship between epithelial cell 3D morphogenesis and the substrate viscosity. We observed that Madin-Darby Canine Kidney (MDCK) cells formed 3D structures on a viscous substrate (Matrigel). The structures appear as a tulip hat. We then changed the substrate viscosity by genipin (GP) treatment. GP is a cross-linker of amino groups. Cells cultured on GP-treated-matrigel changed their 3D morphology in a substrate viscosity-dependent manner. Furthermore, to elucidate the spatial distribution of the cellular contractile force, localization of mono-phosphorylated and di-phosphorylated myosin regulatory light chain (P-MRLCs) was visualized by immunofluorescence. P-MRLCs localized along the periphery of epithelial sheets. Treatment with Y-27632, a Rho-kinase inhibitor, blocked the P-MRLCs localization at the edge of epithelial sheets and halted 3D morphogenesis. Our results indicate that the substrate viscosity, the substrate deformation, and the cellular contractile forces induced by P-MRLCs play crucial roles in 3D morphogenesis.