A poroviscoelastic description of fibrin gels

Elasticity of whole blood clots measured via Volume Controlled Cavity Expansion

Measuring and understanding the mechanical properties of blood clots can provide insights into disease progression and the effectiveness of potential treatments. However, several limitations hinder the use of standard mechanical testing methods to measure the response of soft biological tissues, like blood clots. These tissues can be difficult to mount, and are inhomogeneous, irregular in shape, scarce, and valuable. To remedy this, we employ in this work Volume Controlled Cavity Expansion (VCCE), a technique that was recently developed, to measure local mechanical properties of soft materials in their natural environment. Through highly controlled volume expansion of a water bubble at the tip of an injection needle, paired with simultaneous measurement of the resisting pressure, we obtain a local signature of whole blood clot mechanical response. Comparing this data with predictive theoretical models, we find that a 1-term Ogden model is sufficient to capture the nonlinear elastic response observed in our experiments and produces shear modulus values that are comparable to values reported in the literature. Moreover, we find that bovine whole blood stored at 4 °C for greater than 2 days exhibits a statistically significant shift in the shear modulus from 2.53 ± 0.44 kPa on day 2 (N = 13) to 1.23 ± 0.18 kPa on day 3 (N = 14). In contrast to previously reported results, our samples did not exhibit viscoelastic rate sensitivity within strain rates ranging from 0.22 - 21.1 s-1. By surveying existing data on whole blood clots for comparison, we show that this technique provides highly repeatable and reliable results, hence we propose the more widespread adoption of VCCE as a path forward to building a better understanding of the mechanics of soft biological materials.

Dynamic remodeling of fiber networks with stiff inclusions under compressive loading

The ability of tissues to sustain and withstand mechanical stress is critical to tissue development and healthy tissue maintenance. The mechanical properties of tissues are typically considered to be dominated by the fibrous extracellular matrix (ECM) component of tissues. Fiber network mechanics can capture certain mechanical features of tissues, such as shear strain stiffening, but is insufficient in describing the compressive response of certain tissues and blood clots that are rich in extracellular matrix. To understand the mechanical response of tissues, we employ a contemporary mechanical model, a fibrous network of fibrin embedded with inert bead inclusions that preserve the volume-conserving constraints of cells in tissues. Combining bulk mechanical rheology and a custom imaging device, we show that the presence of inclusions alters the local dynamic remodeling of the networks undergoing uniaxial compressive strains and demonstrate non-affine correlated motion within a fiber-bead network, predicted to stretch fibers in the network and lead to the ability of the network to stiffen under compression, a key feature of real tissues. These findings have important implications for understanding how local structural properties of cells and ECM fibers impact the bulk mechanical response of real tissues. STATEMENT OF SIGNIFICANCE: To understand why real tissue stiffens under compression, we study a model tissue system which also stiffens: a fibrin network embedded with stiff beads. We design a device that images compression of both fiber and fiber-bead networks. Distinct from previous imaging studies, this setup can dynamically capture network deformation on scales larger than single fibers. From the videos, we see fluid flow and network remodeling are both critical to compression behavior. The fiber-bead network has faster fluid flow, reduced network recovery, and correlated motion during network relaxation. We hypothesize that the beads hinder network relaxation of stretched fibers, thus retaining the applied stress and exhibiting stiffening. Our findings reveal important details for modeling tissue mechanics towards optimizing healthcare.

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.

Structure, mechanical properties, and modeling of cyclically compressed pulmonary emboli

Pulmonary embolism occurs when blood flow to a part of the lungs is blocked by a venous thrombus that has traveled from the lower limbs. Little is known about the mechanical behavior of emboli under compressive forces from the surrounding musculature and blood pressure. We measured the stress-strain responses of human pulmonary emboli under cyclic compression, and showed that emboli exhibit a hysteretic stress-strain curve. The fibrin fibers and red blood cells (RBCs) are damaged during the compression process, causing irreversible changes in the structure of the emboli. We showed using electron and confocal microscopy that bundling of fibrin fibers occurs due to compression, and damage is accumulated as more cycles are applied. The stress-strain curves depend on embolus structure, such that variations in composition give quantitatively different responses. Emboli with a high fibrin component demonstrate higher normal stress compared to emboli that have a high RBC component. We compared the compression response of emboli to that of whole blood clots containing various volume fractions of RBCs, and found that RBCs rupture at a certain critical stress. We describe the hysteretic response characteristic of foams, using a model of phase transitions in which the compressed foam is segregated into coexisting rarefied and densified phases whose fractions change during compression. Our model takes account of the rupture of RBCs in the compressed emboli and stresses due to fluid flow through their small pores. Our results can help in classifying emboli as rich in fibrin or rich in red blood cells, and can help in understanding what responses to expect when stresses are applied to thrombi in vivo.

Development of an Injectable Nitric Oxide Releasing Poly(ethylene) Glycol-Fibrin Adhesive Hydrogel

Fibrin microparticles were incorporated into poly(ethylene) glycol (PEG)-fibrinogen hydrogels to create an injectable, composite that could serve as a wound healing support and vehicle to deliver therapeutic factors for tissue engineering. Nitric oxide (NO), a therapeutic agent in wound healing, was loaded into fibrin microparticles by blending S-Nitroso-N-acetyl penicillamine (SNAP) with a fibrinogen solution. The incorporation of microparticles affected swelling behavior and improved tissue adhesivity of composite hydrogels. Controlled NO release was induced via photolytic and thermal activation, and modulated by weight percent of particles incorporated. These NO-releasing composites were non-cytotoxic in culture. Cells maintained morphology, viability, and proliferative character. Fibrin microparticles loaded with SNAP and incorporated into a PEG-fibrinogen matrix, creates a novel injectable composite hydrogel that offers improved tissue adhesivity and inducible NO-release for use as a regenerative support for wound healing and tissue engineering applications.

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.

The biomechanics of subsynovial connective tissue in health and its role in carpal tunnel syndrome

Carpal Tunnel Syndrome (CTS) is the most common surgically treated problem in the hand. Aside from the neuropathy itself, the most common findings are fibrosis of the subsynovial connective tissue (SSCT) and increased intra carpal tunnel pressure. Normally, the SSCT is a multilayer tissue interspersed among the carpal tendons and nerve. As the tendons move, successive SSCT layers are recruited, forming a gliding unit and providing a limit to differential movement. Exceeding this limit, damages the SSCT as has been shown in both cadavers and animal models. This damage leads to a non-inflammatory response with progressive fibrosis and nerve ischemia leaving the SSCT more susceptible to injury. Although the direct consequences for patients are not fully understood, ultrasound research shows that this fibrosis restricts median nerve displacement during tendon loading. This article aims to provide insights into the mechanical properties of SSCT described so far and place it in the context of CTS pathophysiology. A theoretical damage model concerning the SSCT is proposed showing a chain of events and vicious cycles that could lead to the nerve compression as it is found in CTS. Although not complete, this model could explain the pathophysiological pathway of idiopathic CTS.

Phase transitions during compression and decompression of clots from platelet-poor plasma, platelet-rich plasma and whole blood

Blood clots are required to stem bleeding and are subject to a variety of stresses, but they can also block blood vessels and cause heart attacks and ischemic strokes. We measured the compressive response of human platelet-poor plasma (PPP) clots, platelet-rich plasma (PRP) clots and whole blood clots and correlated these measurements with confocal and scanning electron microscopy to track changes in clot structure. Stress-strain curves revealed four characteristic regions, for compression-decompression: (1) linear elastic region; (2) upper plateau or softening region; (3) non-linear elastic region or re-stretching of the network; (4) lower plateau in which dissociation of some newly made connections occurs. Our experiments revealed that compression proceeds by the passage of a phase boundary through the clot separating rarefied and densified phases. This observation motivates a model of fibrin mechanics based on the continuum theory of phase transitions, which accounts for the pre-stress caused by platelets, the adhesion of fibrin fibers in the densified phase, the compression of red blood cells (RBCs), and the pumping of liquids through the clot during compression/decompression. Our experiments and theory provide insights into the mechanical behavior of blood clots that could have implications clinically and in the design of fibrin-based biomaterials.

STATEMENT OF SIGNIFICANCE

The objective of this paper is to measure and mathematically model the compression behavior of various human blood clots. We show by a combination of confocal and scanning electron microscopy that compression proceeds by the passage of a front through the sample that separates a densified region of the clot from a rarefied region, and that the compression/decompression response is reversible with hysteresis. These observations form the basis of a model for the compression response of clots based on the continuum theory of phase transitions. Our studies may reveal how clot rheology under large compression in vivo due to muscle contraction, platelet retraction and hydrodynamic flow varies under various pathophysiological conditions and could inform the design of fibrin based biomaterials.

Delineation of the mechanisms of tendon gliding resistance within the carpal tunnel

BACKGROUND

Forceful, high-velocity, and repetitive manual hand tasks contribute to the onset of carpal tunnel syndrome. This study aimed to isolate and identify mechanisms that contribute to tendon gliding resistance in the carpal tunnel.

METHODS

Eight human cadaver hands (four pairs) were used. Tendon gliding resistance (force, energy, and stiffness) was measured under different conditions: with intact and with divided subsynovial connective tissue, at 2mm/s and 60mm/s tendon excursion velocity, and with and without relaxation time before tendon excursion.

RESULTS

Subsynovial connective tissue stretching substantially contributed to increased gliding resistance force and energy during higher tendon excursion velocities, and subsynovial connective tissue stiffening was observed. Poroelastic properties of the tendon (and possibly the subsynovial connective tissue) also appear to be involved because relaxation time significantly increased gliding resistance force and energy (P<0.01), and the difference in energy and force between high- and low-velocity tendon excursions increased with relaxation time (P=0.01 and P<0.01). Lastly, without relaxation time, no difference in force and energy was observed (P=0.06 and P=0.60), suggesting contact friction.

INTERPRETATION

These findings are consistent with the hypothesis that the mechanics of tendon motion within the carpal tunnel are affected by the integrity of the subsynovial connective tissue. While not tested here, in carpal tunnel syndrome this tissue is known to be the fibrotic, thickened, and less-fluid-permeable. An extrapolation of our findings suggests that these changes in the subsynovial connective tissue of carpal tunnel syndrome patients could increase contact friction and carpal tunnel pressure.

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.

Hydraulic fracture during epithelial stretching

The origin of fracture in epithelial cell sheets subject to stretch is commonly attributed to excess tension in the cells' cytoskeleton, in the plasma membrane, or in cell-cell contacts. Here, we demonstrate that for a variety of synthetic and physiological hydrogel substrates the formation of epithelial cracks is caused by tissue stretching independently of epithelial tension. We show that the origin of the cracks is hydraulic; they result from a transient pressure build-up in the substrate during stretch and compression manoeuvres. After pressure equilibration, cracks heal readily through actomyosin-dependent mechanisms. The observed phenomenology is captured by the theory of poroelasticity, which predicts the size and healing dynamics of epithelial cracks as a function of the stiffness, geometry and composition of the hydrogel substrate. Our findings demonstrate that epithelial integrity is determined in a tension-independent manner by the coupling between tissue stretching and matrix hydraulics.

Deciphering mechanical regulation of chondrogenesis in fibrin-polyurethane composite scaffolds enriched with human mesenchymal stem cells: a dual computational and experimental approach

Fibrin-polyurethane composite scaffolds support chondrogenesis of human mesenchymal stem cells (hMSCs) derived from bone marrow and due to their robust mechanical properties allow mechanical loading in dynamic bioreactors, which has been shown to increase the chondrogenic differentiation of MSCs through the transforming growth factor beta pathway. The aim of this study was to use the finite element method, mechanical testing, and dynamic in vitro cell culture experiments on hMSC-enriched fibrin-polyurethane composite scaffolds to quantitatively decipher the mechanoregulation of chondrogenesis within these constructs. The study identified compressive principal strains as the key regulator of chondrogenesis in the constructs. Although dynamic uniaxial compression did not induce chondrogenesis, multiaxial loading by combined application of dynamic compression and interfacial shear induced significant chondrogenesis at locations where all the three principal strains were compressive and had a minimum magnitude of 10%. In contrast, no direct correlation was identified between the level of pore fluid velocity and chondrogenesis. Due to the high permeability of the constructs, the pore fluid pressures could not be increased sufficiently by mechanical loading, and instead, chondrogenesis was induced by triaxial compressive deformations of the matrix with a minimum magnitude of 10%. Thus, it can be concluded that dynamic triaxial compressive deformations of the matrix is sufficient to induce chondrogenesis in a threshold-dependent manner, even where the pore fluid pressure is negligible.

Fibrin polymerization in blood coagulation-a statistical model

A theoretical model for the growth of fibrin clots is derived. The model is based on a statistical description of the polymerization process underlying the formation of the fibrin polymeric network. The model provides insights regarding the role of various factors, such as thrombin concentration, plasmin concentration, and the local shear rate in the coagulation process. In particular, the effect of these factors on the mechanical properties of the clot is studied. Numerical results are in very good agreement with quantitative and qualitative experimental observations. Importantly, no fitting parameters are used, and all model parameters, such as fibrin persistence length and monomer size, are in accordance with experimental reports.

Multiscale models of thrombogenesis

To restrict the loss of blood follow from the rupture of blood vessels, the human body rapidly forms a clot consisting of platelets and fibrin. However, to prevent pathological clotting within vessels as a result of vessel damage, the response must be regulated. Clots forming within vessels (thrombi) can restrict the flow of blood causing damage to tissues in the flow field. Additionally, fragments dissociating from the primary thrombus (emboli) may lodge and clog vessels in the brain (causing ischemic stroke) or lungs (resulting in pulmonary embolism). Pathologies related to the obstruction of blood flow through the vasculature are the major cause of mortality in the United States. Venous thromboembolic disease alone accounts for 900,000 hospitalizations and 300,000 deaths per year and the incidence will increase as the population ages (Wakefield et al. J Vasc Surg 2009, 49:1620-1623). Thus, understanding the interplay between the many processes involved in thrombus development is of significant biomedical value. In this article, we first review computational models of important subprocesses of hemostasis/thrombosis including coagulation reactions, platelet activation, and fibrin assembly, respectively. We then describe several multiscale models integrating these subprocesses to simulate temporal and spatial development of thrombi. The development of validated computational models and predictive simulations will enable one to explore how the variation of multiple hemostatic factors affects thrombotic risk providing an important new tool for thrombosis research.

The consolidation behavior of silk hydrogels

Hydrogels have mechanical properties and structural features that are similar to load-bearing soft tissues including intervertebral disc and articular cartilage, and can be implanted for tissue restoration or for local release of therapeutic factors. To help predict their performance, mechanical characterization and mathematical modeling are the available methods for use in tissue engineering and drug delivery settings. In this study, confined compression creep tests were performed on silk hydrogels, over a range of concentrations, to examine the phenomenological behavior of the gels under a physiological loading scenario. Based on the observed behavior, we show that the time-dependent response can be explained by a consolidation mechanism, and modeled using Biot's poroelasticity theory. Two observations are in strong support of this modeling framework, namely, the excellent numerical agreement between increasing load step creep data and the linear Terzaghi theory, and the similar values obtained from numerical simulations and direct measurements of the permeability coefficient. The higher concentration gels (8% and 12% w/v) clearly show a strain-stiffening response to creep loading with increasing loads, while the lower concentration gel (4% w/v) does not. A nonlinear elastic constitutive formulation is employed to account for the stiffening. Furthermore, an empirical formulation is used to represent the deformation-dependent permeability.