A constitutive model for vascular tissue that integrates fibril, fiber and continuum levels with application to the isotropic and passive properties of the infrarenal aorta

A fundamental understanding of the mechanical properties of the extracellular matrix (ECM) is critically important to quantify the amount of macroscopic stress and/or strain transmitted to the cellular level of vascular tissue. Structural constitutive models integrate histological and mechanical information, and hence, allocate stress and strain to the different microstructural components of the vascular wall. The present work proposes a novel multi-scale structural constitutive model of passive vascular tissue, where collagen fibers are assembled by proteoglycan (PG) cross-linked collagen fibrils and reinforce an otherwise isotropic matrix material. Multiplicative kinematics account for the straightening and stretching of collagen fibrils, and an orientation density function captures the spatial organization of collagen fibers in the tissue. Mechanical and structural assumptions at the collagen fibril level define a piece-wise analytical stress-stretch response of collagen fibers, which in turn is integrated over the unit sphere to constitute the tissue's macroscopic mechanical properties. The proposed model displays the salient macroscopic features of vascular tissue, and employs the material and structural parameters of clear physical meaning. Likewise, the constitutive concept renders a highly efficient multi-scale structural approach that allows for the numerical analysis at the organ level. Model parameters were estimated from isotropic mean-population data of the normal and aneurysmatic aortic wall and used to predict in-vivo stress states of patient-specific vascular geometries, thought to demonstrate the robustness of the particular Finite Element (FE) implementation. The collagen fibril level of the multi-scale constitutive formulation provided an interface to integrate vascular wall biology and to account for collagen turnover.     

A novel constitutive model for passive right ventricular myocardium: evidence for myofiber–collagen fiber mechanical coupling

The function of right ventricle (RV) is recognized to play a key role in the development of many cardiopulmonary disorders, such as pulmonary arterial hypertension (PAH). Given the strong link between tissue structure and mechanical behavior, there remains a need for a myocardial constitutive model that accurately accounts for right ventricular myocardium architecture. Moreover, most available myocardial constitutive models approach myocardium at the length scale of mean fiber orientation and do not explicitly account for different fibrous constituents and possible interactions among them. In the present work, we developed a fiber-level constitutive model for the passive mechanical behavior of the right ventricular free wall (RVFW). The model explicitly separates the mechanical contributions of myofiber and collagen fiber ensembles, and accounts for the mechanical interactions between them. To obtain model parameters for the healthy passive RVFW, the model was informed by transmural orientation distribution measurements of myo- and collagen fibers and was fit to the mechanical testing data, where both sets of data were obtained from recent experimental studies on non-contractile, but viable, murine RVFW specimens. Results supported the hypothesis that in the low-strain regime, the behavior of the RVFW is governed by myofiber response alone, which does not demonstrate any coupling between different myofiber ensembles. At higher strains, the collagen fibers and their interactions with myofibers begin to gradually contribute and dominate the behavior as recruitment proceeds. Due to the use of viable myocardial tissue, the contribution of myofibers was significant at all strains with the predicted tensile modulus of \(\sim \)32 kPa. This was in contrast to earlier reports (Horowitz et al. 1988) where the contribution of myofibers was found to be insignificant. Also, we found that the interaction between myo- and collagen fibers was greatest under equibiaxial strain, with its contribution to the total stress not exceeding 20 %. The present model can be applied to organ-level computational models of right ventricular dysfunction for efficient diagnosis and evaluation of pulmonary hypertension disorder.     

Multiscale, structure-based modeling for the elastic mechanical behavior of arterial walls

Passive elastic behavior of arterial wall remains difficult to model. Although phenomenological and structural models exist, the question of how the three-dimensional network structure of the collagen in the artery determines its mechanical properties is still open. A model is presented that incorporates a collagen network as well as the noncollagenous material that comprise the artery. The collagen architecture is represented as a network of interconnected fibers, and a neo-Hookean constitutive equation is used to describe the contribution of the noncollagenous matrix. The model is multiscale in that volume-averaging theory is applied to the collagen network, and it is structural in that parameters of the microstructure of the collagen network were considered instead of a macroscopic constitutive law. The computational results provided a good fit to published experimental data for decellularized porcine carotid arteries. The model predicted increased circumferential compliance for increased axial stretch, consistent with previously published reports, and a relatively small sensitivity to open angle. Even at large extensions, the model predicted that the noncollagenous matrix would be in compression, preventing collapse of the collagen network. The incorporation of fiber-fiber interactions led to an accurate model of artery wall behavior with relatively few parameters. The counterintuitive result that the noncollagenous component is in compression during extension and inflation of the tissue suggests that the collagen is important even at small strains, with the noncollagenous components supporting the network, but not resisting the load directly. More accurate representation of the microstructure of the artery wall is needed to explore this issue further.     

A quantitative method to determine the orientation of collagen fibers in the dermis.

We have developed a quantitative microscopic method to determine changes in the orientation of collagen fibers in the dermis resulting from mechanical stress. The method is based on the use of picrosirius red-stained cryostat sections of piglet skin in which collagen fibers reflect light strongly when epipolarization microscopy is used. Digital images of sections were converted into binary images that were analyzed quantitatively on the basis of the length of the collagen fibers in the plane of the section as a measure for the orientation of the fibers. The length of the fibers was expressed in pixels and the mean length of the 10 longest fibers in the image was taken as the parameter for the orientation of the fibers. To test the procedure in an experimental setting, we used skin after 0 and 30 min of skin stretching. The orientation of the fibers in sections of control skin differed significantly from the orientation of fibers in sections of skin that was stretched mechanically for 30 min [76 +/- 15 (n=5) vs 132 +/- 36 (n=5)]. The method described here is a relatively simple way to determine (changes in) the orientation of individual collagen fibers in connective tissue and can also be applied for analysis of the orientation of any other structural element in tissues so long as a representative binary image can be created.     

A unified multiscale mechanical model for soft collagenous tissues with regular fiber arrangement

In this paper the mechanical response of soft collagenous tissues with regular fiber arrangement (RSCTs) is described by means of a nanoscale model and a two-step micro–macro homogenization technique. The non-linear collagen constitutive behavior is modeled at the nanoscale by a novel approach accounting for entropic mechanisms as well as stretching effects occurring in collagen molecules. Crimped fibers are reduced to equivalent straight ones at the microscale and the constitutive response of RSCTs at the macroscale is formulated by homogenizing a fiber reinforced material. This approach has been applied to different RSCTs (tendon, periodontal ligament and aortic media), resulting effective and accurate as proved by the excellent agreement with available experimental data. The model is based on few parameters, directly related to histological and morphological evidences and whose sensitivity has been widely investigated. Applications to simulation of some physiopathological mechanisms are also proposed, providing confirmation of clinical evidences and quantitative indications helpful for clinical practice.     

Inferring spatial variations of microstructural properties from macroscopic mechanical response

Disease alters tissue microstructure, which in turn affects the macroscopic mechanical properties of tissue. In elasticity imaging, the macroscopic response is measured and is used to infer the spatial distribution of the elastic constitutive parameters. When an empirical constitutive model is used, these parameters cannot be linked to the microstructure. However, when the constitutive model is derived from a microstructural representation of the material, it allows for the possibility of inferring the local averages of the spatial distribution of the microstructural parameters. This idea forms the basis of this study. In particular, we first derive a constitutive model by homogenizing the mechanical response of a network of elastic, tortuous fibers. Thereafter, we use this model in an inverse problem to determine the spatial distribution of the microstructural parameters. We solve the inverse problem as a constrained minimization problem and develop efficient methods for solving it. We apply these methods to displacement fields obtained by deforming gelatin–agar co-gels and determine the spatial distribution of agar concentration and fiber tortuosity, thereby demonstrating that it is possible to image local averages of microstructural parameters from macroscopic measurements of deformation.     

Collagen fiber alignment does not explain mechanical anisotropy in fibroblast populated collagen gels

Many load-bearing soft tissues exhibit mechanical anisotropy. In order to understand the behavior of natural tissues and to create tissue engineered replacements, quantitative relationships must be developed between the tissue structures and their mechanical behavior. We used a novel collagen gel system to test the hypothesis that collagen fiber alignment is the primary mechanism for the mechanical anisotropy we have reported in structurally anisotropic gels. Loading constraints applied during culture were used to control the structural organization of the collagen fibers of fibroblast populated collagen gels. Gels constrained uniaxially during culture developed fiber alignment and a high degree of mechanical anisotropy, while gels constrained biaxially remained isotropic with randomly distributed collagen fibers. We hypothesized that the mechanical anisotropy that developed in these gels was due primarily to collagen fiber orientation. We tested this hypothesis using two mathematical models that incorporated measured collagen fiber orientations: a structural continuum model that assumes affine fiber kinematics and a network model that allows for nonaffine fiber kinematics. Collagen fiber mechanical properties were determined by fitting biaxial mechanical test data from isotropic collagen gels. The fiber properties of each isotropic gel were then used to predict the biaxial mechanical behavior of paired anisotropic gels. Both models accurately described the isotropic collagen gel behavior. However, the structural continuum model dramatically underestimated the level of mechanical anisotropy in aligned collagen gels despite incorporation of measured fiber orientations; when estimated remodeling-induced changes in collagen fiber length were included, the continuum model slightly overestimated mechanical anisotropy. The network model provided the closest match to experimental data from aligned collagen gels, but still did not fully explain the observed mechanics. Two different modeling approaches showed that the level of collagen fiber alignment in our uniaxially constrained gels cannot explain the high degree of mechanical anisotropy observed in these gels. Our modeling results suggest that remodeling-induced redistribution of collagen fiber lengths, nonaffine fiber kinematics, or some combination of these effects must also be considered in order to explain the dramatic mechanical anisotropy observed in this collagen gel model system.     

A structural constitutive model for collagenous cardiovascular tissues incorporating the angular fiber distribution

Accurate constitutive models are required to gain further insight into the mechanical behavior of cardiovascular tissues. In this study, a structural constitutive framework for cardiovascular tissues is introduced that accounts for the angular distribution of collagen fibers. To demonstrate its capabilities, the model is applied to study the biaxial behavior of the arterial wall and the aortic valve. The pressure-radius relationships of the arterial wall accurately describe experimentally observed sigma-shaped curves. In addition, the nonlinear and anisotropic mechanical properties of the aortic valve can be analyzed with the proposed model. We expect that the current model offers strong possibilities to further investigate the complex mechanical behavior of cardiovascular tissues, including their response to mechanical stimuli.     

A physically motivated constitutive model for cell-mediated compaction and collagen remodeling in soft tissues

Collagen is the main load-bearing component of many soft tissues and has a large influence on the mechanical behavior of tissues when exposed to mechanical loading. Therefore, it is important to increase our understanding of collagen remodeling in soft tissues to understand the mechanisms behind pathologies and to control the development of the collagen network in engineered tissues. In the present study, a constitutive model was developed by coupling a recently developed model describing the orientation and contractile stresses exerted by cells in response to mechanical stimuli to physically motivated collagen remodeling laws. In addition, cell-mediated contraction of the collagen fibers was included as a mechanism for tissue compaction. The model appeared to be successful in predicting a range of experimental observations, which are (1) the change in transition stretch of periosteum after remodeling at different applied stretches, (2) the compaction and alignment of collagen fibers in tissue-engineered strips, (3) the fiber alignment in cruciform gels with different arm widths, and (4) the alignment of collagen fibers in engineered vascular grafts. Moreover, by changing the boundary conditions, the model was able to predict a helical architecture in the vascular graft without assuming the presence of two helical fiber families a priori. Ultimately, this model may help to increase our understanding of collagen remodeling in physiological and pathological conditions, and it may provide a tool for determining the optimal experimental conditions for obtaining native-like collagen architectures in engineered tissues.     

A failure model for ligaments

We propose a failure model for ligament which assumes that sequential uncrimping and stretching of collagen fibers is responsible for the mechanical response of ligament. We further assume that the fibers rupture sequentially and in a brittle, strain-limited manner. The model was fit to stress strain curves obtained from medial collateral ligaments of New Zealand White rabbits from two age groups (4 and 7 months). The model indicated that collagen modulus values ranged from 300 to 680 MPa and that fiber failure strains ranged from 6 to 22%. The model provides a convenient means of describing the elastic and failure response of ligament using four structurally based parameters.     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Characterisation of Anisotropic and Non-linear Behaviour of Human Skin In Vivo

The aim of this work is to characterise the in-plane mechanical behaviour of human skin in vivo. For this purpose the structural skin model proposed by Lanir [1] is employed and a mixed numerical-experimental method is developed. The numerical-experimental method is based on the confrontation of measured data from an experiment, with calculated data from a finite element model, eventually leading to the determination of some of the parameters of a constitutive model, in the present case Lanir's Skin Model. Since collagen, the main constituent of skin, dominates the anisotropic and non-linear behaviour of skin, the parameters of Lanir's Skin Model concerning the mechanical behaviour of the collagen fibres are estimated. In vivo experiments were carried out on the volar forearm. During the experiments, reaction forces and the displacement field at different states of deformation are measured. Both data sets are used for the determination of the parameters.     

Experimentally validated microstructural 3D constitutive model of coronary arterial media

Accurate modeling of arterial response to physiological or pathological loads may shed light on the processes leading to initiation and progression of a number of vascular diseases and may serve as a tool for prediction and diagnosis. In this study, a microstructure based hyperelastic constitutive model is developed for passive media of porcine coronary arteries. The most general model contains 12 independent parameters representing the three-dimensional inner fibrous structure of the media and includes the effects of residual stresses and osmotic swelling. Parameter estimation and model validation were based on mechanical data of porcine left anterior descending (LAD) media under radial inflation, axial extension, and twist tests. The results show that a reduced four parameter model is sufficient to reliably predict the passive mechanical properties. These parameters represent the stiffness and the helical orientation of each lamellae fiber and the stiffness of the interlamellar struts interconnecting these lamellae. Other structural features, such as orientational distribution of helical fibers and anisotropy of the interlamellar network, as well as possible transmural distribution of structural features, were found to have little effect on the global media mechanical response. It is shown that the model provides good predictions of the LAD media twist response based on parameters estimated from only biaxial tests of inflation and extension. In addition, good predictive capabilities are demonstrated for the model behavior at high axial stretch ratio based on data of law stretches.     

Contributions of Glycosaminoglycans to Collagen Fiber Recruitment in Constitutive Modeling of Arterial Mechanics

The contribution of glycosaminoglycans (GAGs) to the biological and mechanical functions of biological tissue has emerged as an important area of research. GAGs provide structural basis for the organization and assembly of extracellular matrix (ECM). The mechanics of tissue with low GAG content can be indirectly affected by the interaction of GAGs with collagen fibers, which have long been known to be one of the primary contributors to soft tissue mechanics. Our earlier study showed that enzymatic GAG depletion results in straighter collagen fibers that are recruited at lower levels of stretch, and a corresponding shift in earlier arterial stiffening (Mattson et al., 2016). In this study, the effect of GAGs on collagen fiber recruitment was studied through a structure-based constitutive model. The model incorporates structural information, such as fiber orientation distribution, content, and recruitment of medial elastin, medial collagen, and adventitial collagen fibers. The model was first used to study planar biaxial tensile stress-stretch behavior of porcine descending thoracic aorta. Changes in elastin and collagen fiber orientation distribution, and collagen fiber recruitment were then incorporated into the model in order to predict the stress-stretch behavior of GAG depleted tissue. Our study shows that incorporating early collagen fiber recruitment into the model predicts the stress-stretch response of GAG depleted tissue reasonably well (rms = 0.141); considering further changes of fiber orientation distribution does not improve the predicting capability (rms = 0.149). Our study suggests an important role of GAGs in arterial mechanics that should be considered in developing constitutive models.     

A Constitutive Model for a Maturing Fibrin Network

Blood clot formation is crucial to maintain normal physiological conditions but at the same time involved in many diseases. The mechanical properties of the blood clot are important for its functioning but complicated due to the many processes involved. The main structural component of the blood clot is fibrin, a fibrous network that forms within the blood clot, thereby increasing its mechanical rigidity. A constitutive model for the maturing fibrin network is developed that captures the evolving mechanical properties. The model describes the fibrin network as a network of fibers that become thicker in time. Model parameters are related to the structural properties of the network, being the fiber length, bending stiffness, and mass-length ratio. Results are compared with rheometry experiments in which the network maturation is followed in time for various loading frequencies and fibrinogen concentrations. Three parameters are used to capture the mechanical behavior including the mass-length ratio. This parameter agrees with values determined using turbidimetry experiments and is subsequently used to derive the number of protofibrils and fiber radius. The strength of the model is that it describes the mechanical properties of the maturing fibrin network based on it structural quantities. At the same time the model is relatively simple, which makes it suitable for advanced numerical simulations of blood clot formation during flow in blood vessels.     

A Constitutive Model for a Maturing Fibrin Network

Blood clot formation is crucial to maintain normal physiological conditions but at the same time involved in many diseases. The mechanical properties of the blood clot are important for its functioning but complicated due to the many processes involved. The main structural component of the blood clot is fibrin, a fibrous network that forms within the blood clot, thereby increasing its mechanical rigidity. A constitutive model for the maturing fibrin network is developed that captures the evolving mechanical properties. The model describes the fibrin network as a network of fibers that become thicker in time. Model parameters are related to the structural properties of the network, being the fiber length, bending stiffness, and mass-length ratio. Results are compared with rheometry experiments in which the network maturation is followed in time for various loading frequencies and fibrinogen concentrations. Three parameters are used to capture the mechanical behavior including the mass-length ratio. This parameter agrees with values determined using turbidimetry experiments and is subsequently used to derive the number of protofibrils and fiber radius. The strength of the model is that it describes the mechanical properties of the maturing fibrin network based on it structural quantities. At the same time the model is relatively simple, which makes it suitable for advanced numerical simulations of blood clot formation during flow in blood vessels.     

Stretching type II collagen with optical tweezers

Collagens are the most abundant proteins of vertebrates and they provide mechanical and supportive functions in a wide range of connective tissues. Knowledge of the mechanical properties of single collagen molecules is essential in studying the self-assembly of collagen, the interaction between cells and extracellular matrix, the etiology of tissue degeneration and mechanism of regeneration, and the relationship between the structures and mechanical properties of tissues. Here we stretched single type II collagen molecules in neutral pH solution using optical tweezers. The molecular parameters of collagen were obtained by fitting force-extension curves into worm-like chain elasticity model. The molecule length of type II collagen monomer was 295.8 nm. The persistence length of type II collagen monomer was 11.2 nm. These observations indicate that collagen molecules are flexible rather than rigid rod molecules at neutral pH solution.