On arterial fiber dispersion and auxetic effect

There are two polar contemporary approaches to the constitutive modeling of arterial wall with anisotropy induced by collagen fibers. The first one is based on the angular integration (AI) of the strain energy on a unit sphere for the analytically defined fiber dispersion. The second one is based on the introduction of the generalized structure tensors (GST). AI approach is very involved computationally while GST approach requires somewhat complicated procedure for the exclusion of compressed fibers.We present some middle ground models, which are based on the use of 16 and 8 structure tensors. These models are moderately involved computationally and they allow excluding compressed fibers easily. We use the proposed models to study the role of the fiber dispersion in the constitutive modeling of the arterial wall. Particularly, we study the auxetic effect which can appear in anisotropic materials. The effect means thickening of the tissue in the direction perpendicular to its stretching. Such an effect was not observed in experiments while some simple anisotropic models do predict it. We show that more accurate account of the fiber dispersion suppresses the auxetic effect in a qualitative agreement with experimental observations.     

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.     

Micromechanical modeling of the epimysium of the skeletal muscles

A micromechanical model has been developed to investigate the mechanical properties of the epimysium. In the present model, the collagen fibers in the epimysium are embedded randomly in the ground substance. Two parallel wavy collagen fibers and the surrounding ground substance are used as the repeat unit (unit cell), and the epimysium is considered as an aggregate of unit cells. Each unit cell is distributed in the epimysium with some different angle to the muscle fiber direction. The model allows the progressive straightening of the collagen fiber as well as the effects of fiber reorientation. The predictions of the model compare favorably against experiment. The effects of the collagen fiber volume fraction, collagen fiber waviness at the rest length and the mechanical properties of the collagen fibers and the ground substance are analyzed. This model allows the analysis of mechanical behavior of most soft tissues if appropriate experimental data are available.     

Improved prediction of the collagen fiber architecture in the aortic heart valve

Living tissues show an adaptive response to mechanical loading by changing their internal structure and morphology. Understanding this response is essential for successful tissue engineering of load-bearing structures, such as the aortic valve. In this study, mechanically induced remodeling of the collagen architecture in the aortic valve was investigated. It was hypothesized that, in uniaxially loaded regions, the fibers aligned with the tensile principal stretch direction. For biaxial loading conditions, on the other hand, it was assumed that the collagen fibers aligned with directions situated between the principal stretch directions. This hypothesis has already been applied successfully to study collagen remodeling in arteries. The predicted fiber architecture represented a branching network and resembled the macroscopically visible collagen bundles in the native leaflet. In addition, the complex biaxial mechanical behavior of the native valve could be simulated qualitatively with the predicted fiber directions. The results of the present model might be used to gain further insight into the response of tissue engineered constructs during mechanical conditioning.     

Characterization of collagen fiber architecture in the canine diaphragmatic central tendon.

The diaphragmatic central tendon (DCT), a collagenous soft tissue membrane, acts as a mechanical buffer between the costal and crural muscles. Its direction of mechanical anisotropy has been shown to correspond to the collagen fiber preferred directions. These preferred directions were determined by gross histological examination, and were thus qualitative. In this work we quantified the collagen fiber architecture throughout the DCT using small angle light scattering (SALS). Helium-Neon laser light was passed through tendon specimens and the resultant scattered light distribution, which characterized the local collagen fiber architecture, was recorded with a linear array of five photodiodes. Throughout the DCT two distinct collagen fiber populations were consistently found. For each population three parameters were determined: 1) the preferred directions of collagen fibers, 2) the volume fraction (Vf) of fibers, 3) OI, an orientation index, which ranges from 0 percent for a random network to 100 percent for a perfectly oriented network. Vector maps were used to display results from 1) and 2), and showed a primary group (G1) going from the crural to costal muscles and a secondary one (G2) running perpendicular to G1. Comparisons of Vf between G1 and G2 showed that G1 contained about three times as many fibers as G2, a ratio similar to that found for the degree of mechanical anisotropy. OI were found to be about 60 percent, indicating a high degree of orientation, with no significant regional or population differences (p less than 0.05). These quantitative results suggest that throughout the DCT the degree of mechanical anisotropy is controlled exclusively by Vf.     

Remodelling of continuously distributed collagen fibres in soft connective tissues

Extracellular matrix remodelling plays an essential role in tissue engineering of load-bearing structures. The goal of this study is to model changes in collagen fibre content and orientation in soft connective tissues due to mechanical stimuli. A theory is presented describing the mechanical condition within the tissue and accounting for the effects of collagen fibre alignment and changes in fibre content. A fibre orientation tensor is defined to represent the continuous distribution of collagen fibre directions. A constitutive model is introduced to relate the fibre configuration to the macroscopic stress within the material. The constitutive model is extended with a structural parameter, the fibre volume fraction, to account for the amount of fibres present within the material. It is hypothesised that collagen fibre reorientation is induced by macroscopic deformations and the amount of collagen fibres is assumed to increase with the mean fibre stretch. The capabilities of the model are demonstrated by considering remodelling within a biaxially stretched cube. The model is then applied to analyse remodelling within a closed stented aortic heart valve. The computed preferred fibre orientation runs from commissure to commissure and resembles the fibre directions in the native aortic valve.     

Structural modeling reveals microstructure-strength relationship for human ascending thoracic aorta

High lethality of aortic dissection necessitates accurate predictive metrics for dissection risk assessment. The not infrequent incidence of dissection at aortic diameters <5.5 cm, the current threshold guideline for surgical intervention (Nishimura et al., 2014), indicates an unmet need for improved evidence-based risk stratification metrics. Meeting this need requires a fundamental understanding of the structural mechanisms responsible for dissection evolution within the vessel wall. We present a structural model of the repeating lamellar structure of the aortic media comprised of elastic lamellae and collagen fiber networks, the primary load-bearing components of the vessel wall. This model was used to assess the role of these structural features in determining in-plane tissue strength, which governs dissection initiation from an intimal tear. Ascending aortic tissue specimens from three clinically-relevant patient populations were considered: non-aneurysmal aorta from patients with morphologically normal tricuspid aortic valve (CTRL), aneurysmal aorta from patients with tricuspid aortic valve (TAV), and aneurysmal aorta from patients with bicuspid aortic valve (BAV). Multiphoton imaging derived collagen fiber organization for each patient cohort was explicitly incorporated in our model. Model parameters were calibrated using experimentally-measured uniaxial tensile strength data in the circumferential direction for each cohort, while the model was validated by contrasting simulated tissue strength against experimentally-measured strength in the longitudinal direction. Orientation distribution, controlling the fraction of loaded collagen fibers at a given stretch, was identified as a key feature governing anisotropic tissue strength for all patient cohorts.     

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.     

Remodelling of the angular collagen fiber distribution in cardiovascular tissues

Understanding collagen fiber remodelling is desired to optimize the mechanical conditioning protocols in tissue-engineering of load-bearing cardiovascular structures. Mathematical models offer strong possibilities to gain insight into the mechanisms and mechanical stimuli involved in these remodelling processes. In this study, a framework is proposed to investigate remodelling of angular collagen fiber distribution in cardiovascular tissues. A structurally based model for collagenous cardiovascular tissues is extended with remodelling laws for the collagen architecture, and the model is subsequently applied to the arterial wall and aortic valve. For the arterial wall, the model predicts the presence of two helically arranged families of collagen fibers. A branching, diverging hammock-type fiber architecture is predicted for the aortic valve. It is expected that the proposed model may be of great potential for the design of improved tissue engineering protocols and may give further insight into the pathophysiology of cardiovascular diseases.     

A multiaxial constitutive law for mammalian left ventricular myocardium in steady-state barium contracture or tetanus.

The constitutive law of the material comprising any structure is essential for mechanical analysis since this law enables calculation of the stresses from the deformations and vice versa. To date, there is no constitutive law for actively contracting myocardial tissue. Using 2,3-butanedione monoxime to protect the myocardium from mechanical trauma, we subjected thin midwall slices of rabbit myocardium to multiaxial stretching first in the passive state and then during steady-state barium contracture or during tetani in ryanodine-loaded tissue. Assuming transverse isotropy in both the passive and active conditions, we used our previously described methods (Humphrey et al., 1990a) to obtain both passive and active constitutive laws. The major results of this study are: (1) This is the first multiaxial constitutive law for actively contracting mammalian myocardium. (2) The functional forms of the constitutive law for barium contracture and ryanodine-induced tetani are the same but differ from those in the passive state. Hence, one cannot simply substitute differing values for the coefficients of the passive law to describe the active tissue properties. (3) There are significant stresses developed in the cross-fiber direction (more than 40 percent of those in the fiber direction) that cannot be attributed to either deformation effects or nonparallel muscle fibers. These results provide the foundation for future mechanical analyses of the heart.     

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.     

Inflation of an artery leading to aneurysm formation and rupture

Formation and rupture of aneurysms due to the inflation of an artery with collagen fibers distributed in two preferred directions, subjected to internal pressure and axial stretch are examined within the framework of nonlinear elasticity. A two layer tube model with a fiber-reinforced composite based incompressible anisotropic hyperelastic constitutive material is employed to model the stress-strain behavior of the artery wall with distributed collagen fibers. The artery wall takes up a uniform inflation deformation, and there are no aneurysms in the artery under the normal condition. But an aneurysm may be formed in arteries when the stiffness of the fibers is decreased to a certain value or the direction of the fibers is changed to a certain degree towards the circumferential direction. The aneurysm may expand to much large extent and become complex in shape. One portion of the aneurysm becomes highly distended as a bubble while the rest remains lightly inflated. The rupture of the aneurysm is discussed along with the distribution of stresses. Critical pressures and the rupture pressures are given for different collagen fiber orientations or stiffness. Furthermore, the stability of the solutions is discussed to explain the formation of aneurysm.     

Efficient computational simulation of actin stress fiber remodeling

Understanding collagen and stress fiber remodeling is essential for the development of engineered tissues with good functionality. These processes are complex, highly interrelated, and occur over different time scales. As a result, excessive computational costs are required to computationally predict the final organization of these fibers in response to dynamic mechanical conditions. In this study, an analytical approximation of a stress fiber remodeling evolution law was derived. A comparison of the developed technique with the direct numerical integration of the evolution law showed relatively small differences in results, and the proposed method is one to two orders of magnitude faster.     

A closed-form structural model of planar fibrous tissue mechanics

Structural models of tissue mechanics, in which the tissue is represented as a sum or integral of fiber contributions for a distribution of fiber orientations, are a popular tool to represent the complex mechanical behavior of soft tissues. A significant practical challenge, however, is evaluation of the integral that defines the stress. Numerical integration is accurate but computationally demanding, posing an impediment to incorporation of structural models into large-scale finite-element simulations. In this paper, a closed-form analytic evaluation of the integral is derived for fibers distributed according to a von Mises distribution and an exponential fiber stress–strain law.     

Hyperelastic modeling of the combined effects of tissue swelling and deformation-related collagen renewal in fibrous soft tissue

A continuum mechanics constitutive model is presented for the interaction between swelling and collagen remodeling in biological soft tissue. The model is inherently two-way: swelling stretches the collagen fibers which affects their rate of degradation—the remodeled fibrous microarchitecture provides selective directional stiffening that causes the swollen tissue to expand more in the unreinforced directions. The constitutive model specifically treats stretch-stabilization wherein the rate of enzymatic-induced degradation of collagen is a decreasing function of fiber stretch. New collagen replacement takes place in a generally swollen environment, and this synthesis is tracked as a function of time by means of a time integration scheme that accounts for the historical sequence of collagen recreation. The model allows for the specification of the collagen pre-stretch at the time of first synthesis, thus allowing for the consideration of either initially limp replacement fiber or initially pre-tensioned replacement fiber. Loading and swelling that occurs on time scales that are commensurate with the natural time scales for fiber degradation and replacement lead to the consideration of time-integral constitutive equations. Loading and swelling that take place on time scales that are very different from that of the remodeling time scales provide a simplified treatment in which there are definite notions of a short-time instantaneous response and also a large-time approach to a steady-state condition of homeostasis.     

Constitutive modeling of mouse carotid arteries using experimentally measured microstructural parameters

Changes in the local mechanical environment and tissue mechanical properties affect the biological activity of cells and play a key role in a variety of diseases, such as cancer, arthritis, nephropathy, and cardiovascular disease. Constitutive relations have long been used to predict the local mechanical environment within biological tissues and to investigate the relationship between biological responses and mechanical stimuli. Recent constitutive relations for soft tissues consider both material and structural properties by incorporating parameters that describe microstructural organization, such as fiber distributions, fiber angles, fiber crimping, and constituent volume fractions. The recently developed technique of imaging the microstructure of a single artery as it undergoes multiple deformations provides quantitative structural data that can reduce the number of estimated parameters by using parameters that are truly experimentally intractable. Here, we employed nonlinear multiphoton microscopy to quantify collagen fiber organization in mouse carotid arteries and incorporated measured fiber distribution data into structurally motivated constitutive relations. Microscopy results demonstrate that collagen fibers deform in an affine manner over physiologically relevant deformations. The incorporation of measured fiber angle distributions into constitutive relations improves the model's predictive accuracy and does not significantly reduce the goodness of fit. The use of measured structural parameters rather than estimated structural parameters promises to improve the predictive capabilities of the local mechanical environment, and to extend the utility of intravital imaging methods for estimating the mechanical behavior of tissues using in situ structural information.     

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.     

The influence of fiber dispersion on the mechanical response of aortic tissues in health and disease: a computational study

Changes in the structural components of aortic tissues have been shown to play a significant role in the pathogenesis of aortic degeneration. Therefore, reliable stress analyses require a suitable and meaningful constitutive model that captures micro-structural changes. As recent data show, in-plane and out-of-plane collagen fiber dispersions vary significantly between healthy and aneurysmatic aortic walls. The aim of this study is to computationally investigate the influence of fiber dispersion on the mechanical response of aortic tissues in health and disease. In particular, the influence of three different fiber dispersions is studied: (i) non-rotationally symmetric dispersion, the most realistic assumption for aortic tissues; (ii) transversely isotropic dispersion, a special case; (iii) perfectly aligned fibers (no dispersion in either plane), another special case. Explicit expressions for the stress and elasticity tensors as needed for the implementation in a finite element code are provided. Three representative numerical examples are studied: planar biaxial extension, inflation of residually stressed and pre-stretched aortic segments and inflation of an idealized abdominal aortic aneurysm (AAA) geometry. For the AAA geometry the case of isotropic dispersion is additionally analyzed. Documented structural and mechanical parameters are taken from human aortas (healthy media/adventitia and AAA). The influence of fiber dispersions upon magnitudes and distributions of stresses and deformations are presented and analyzed. Stresses vary significantly, especially in the AAA case, where material stiffening is significantly influenced by fiber dispersion. The results highlight the need to incorporate the structural differences into finite element simulations to obtain more accurate stress predictions. Additionally, results show the capability of one constitutive model to represent different scenarios of aortic micro-structures allowing future studies of collagen reorientation during disease progression.     

On the Presence of Affine Fibril and Fiber Kinematics in the Mitral Valve Anterior Leaflet

In this study, we evaluated the hypothesis that the constituent fibers follow an affine deformation kinematic model for planar collagenous tissues. Results from two experimental datasets were utilized, taken at two scales (nanometer and micrometer), using mitral valve anterior leaflet (MVAL) tissues as the representative tissue. We simulated MVAL collagen fiber network as an ensemble of undulated fibers under a generalized two-dimensional deformation state, by representing the collagen fibrils based on a planar sinusoidally shaped geometric model. The proposed approach accounted for collagen fibril amplitude, crimp period, and rotation with applied macroscopic tissue-level deformation. When compared to the small angle x-ray scattering measurements, the model fit the data well, with an r² = 0.976. This important finding suggests that, at the homogenized tissue-level scale of ∼1 mm, the collagen fiber network in the MVAL deforms according to an affine kinematics model. Moreover, with respect to understanding its function, affine kinematics suggests that the constituent fibers are largely noninteracting and deform in accordance with the bulk tissue. It also suggests that the collagen fibrils are tightly bounded and deform as a single fiber-level unit. This greatly simplifies the modeling efforts at the tissue and organ levels, because affine kinematics allows a straightforward connection between the macroscopic and local fiber strains. It also suggests that the collagen and elastin fiber networks act independently of each other, with the collagen and elastin forming long fiber networks that allow for free rotations. Such freedom of rotation can greatly facilitate the observed high degree of mechanical anisotropy in the MVAL and other heart valves, which is essential to heart valve function. These apparently novel findings support modeling efforts directed toward improving our fundamental understanding of tissue biomechanics in healthy and diseased conditions.