Reconstructing vascular homeostasis by growth-based prestretch and optimal fiber deposition
Introduction
Mathematical modeling of the mechanical behavior of vascular tissue requires consideration of many important factors. While the mechanical characterization and constituent models of arterial tissue have been well developed, a challenging factor is the consideration of residual stress in the material in vivo. This residual stress determines the overall stress state of the material in the loaded in vivo configuration, as well as vascular homeostasis, which represents the preferred mechanical configuration in which vascular constituents operate (Fung, 1991). Thus, while proper constituent models are necessary to capture the detailed mechanical behavior of vascular tissue, the inclusion of residual stress is also an essential factor for predicting biomechanical behavior in vivo and vascular remodeling (Holzapfel and Ogden, 2010).
Another important factor in modeling arterial tissue is the consideration of fiber orientation. Arteries are generally composed of three layers (intima, media, and adventitia) and the major mechanical constituents are elastin, collagen, and smooth muscle fibers (Humphrey, 2008). Elastin is effective under small strain, while collagen bears the majority of loading when deformation is large (Roach and Burton, 1957), and smooth muscle fibers provide additional and vasoactive support (Fung, 1993). To model these components, different constitutive relations have been proposed, including Fung-type models (Chuong and Fung, 1983) and formulations proposed by Holzapfel et al. (2000). Notably, the arterial wall is often considered to be a two-layer fiber-reinforced material. A non-collagenous ground matrix, including elastin, is typically described as an isotropic neo-Hookean material (Dorrington and McCrum, 1977). The collagen and smooth muscle impart anisotropic behavior and are often described by a Fung-type response (Holzapfel et al., 2000; Bellini et al., 2014). To characterize the deformation of the various constituents, a constrained mixture model (Humphrey and Rajagopal, 2002; Cyron et al., 2016) is often employed, particularly in vascular growth and remodeling studies (Figueroa et al., 2009a; Zeinali-Davarani et al., 2011; Wu and Shadden, 2015). Fiber orientations have traditionally been prescribed from statistical analyses of histological data (Holzapfel, 2006) and rule-based approaches (Augustin et al., 2014), with more recent works (Gasser et al., 2006; Holzapfel and Ogden, 2017) considering the dispersion of fiber orientation. Nonetheless, prescription of fiber orientation is nontrivial in realistic branched vascular models.
The in vivo stress state of vascular tissue relies critically on residual stress and fiber orientation. Namely, the homeostatic state is hypothesized to be defined by relatively uniform stress distribution in the transmural direction of each layer in the vessel wall, and that non-uniform residual stress helps to maintain this condition (Fung, 1991). In idealized cylindrical vascular geometries, the residual stress field is often specified based on an opening angle (Alastrué et al., 2007), however, this approach is often inadequate in realistic vascular geometries. Pierce et al. (2015) prescribed residual stress by specifying the deformation gradient to map the stress-free configuration to the actual mixture configuration of the vessel wall. Bellini et al. (2014) incorporated a prestretch ratio into the total deformation to model the effect of residual stresses and showed that the stress distribution tends to be more uniform when residual stress is included. Holzapfel et al. investigated the 3D behavior of the residual stress experimentally (Holzapfel et al., 2007) and developed a theoretical framework to incorporate residual stress in different vascular layers (Holzapfel and Ogden, 2010). Previous works have often required parameter information that can be difficult to obtain or is of challenging applicability when realistic geometries are considered. Likewise, the residual stress distribution is rarely considered when the vascular homeostatic state needs to be reconstructed in non-idealized geometries, see, e.g., (Maes et al., 2019; Mousavi et al., 2017). To the best of our knowledge, this is the first time that a residual stress field is reconstructed in general vascular geometries with anomalies such as bifurcations.
We propose a two-stage approach to the reconstruction of the vascular homeostatic state. The first stage aims to generate an appropriate residual stress level for a given configuration of the vessel, and the second stage aims to generate optimal fiber alignments based on the stress field. Inspired by kinematic G&R models (Rodriguez et al., 1994), (Taber, 1998), and (Taber and Humphrey, 2001), we introduce a prestretch tensor, obtained through an iterative process to match a prescribed homeostatic stress distribution at physiologic pressure. The deposition angle of the collagen fibers is then defined by the solution of an optimization problem, which enables the vessel to sustain the stress with minimal amount of biomass, while stochasticity is incorporated into this process to account for fiber dispersion. We propose supervised learning to predict the prestretch in vascular bifurcations or other geometric anomalies, where theory-based methods struggle due to topological singularities. Based on this overall framework, the numerical stability of the generation of the residual stress field and fiber directions is improved for realistic vascular geometries. Beyond addressing the vascular homeostatic state reconstruction problem, the method proposed here provides a starting point for general biomechanical analyses, such as vascular tissue mechanics or fluid-structure interaction simulations, where residual stress and proper fiber alignment set the stage for the mechanical response of the tissue.
Section snippets
Constrained mixture model
The vessel wall is considered a constrained mixture (Humphrey and Rajagopal, 2002) of collagen, elastin, and smooth muscle fibers, where all constituents deform as a continuum. The turnover of collagen and smooth muscle fibers is a process of degradation of existing material and production of new material, described mathematically aswhere k denotes different collagen or smooth muscle fiber families; Mk(t) is the local mass density of constituent k in the
Results
The framework proposed above was implemented to first generate a residual stress field based on the evolution equation (15) and then to iteratively generate collagen fiber directions based on the optimality conditions (28) along with deposition randomness (29). The code was developed in Python, and the finite element solution was implemented using FEniCS (Logg et al., 2012). All computations were performed using linear tetrahedral elements with an average edge size smaller than 0.7 mm. The
Discussion
A computational framework to generate the vascular homeostatic state has been proposed, which includes the generation of appropriate prestretch distribution and the specification of physiological vascular fiber deposition directions. A benefit of this framework is that in order to generate the residual stress distribution only a homeostatic stress value (or distribution) needs to be specified. In particular, the specification of the natural configurations of the different constituents and the
Conclusions
A framework was presented to reconstruct the vascular homeostatic state by developing a method to compute the residual stress distribution and fiber directions. The residual stress is implicitly incorporated by iteratively computing a growth prestretch based on a prescribed homeostatic stress distribution. Vascular fiber directions are further incorporated based on the solution of a physiologically-motivated optimization problem. The methodology was applied to an idealized cylindrical geometry
Author contributions
Jiacheng Wu: Conceptualization, Methodology, Software, Visualization, Formal Analysis, Writing - Original Draft. Christoph Augustin: Conceptualization, Methodology, Visualization, Formal Analysis, Writing - Original Draft. Shawn Shadden: Conceptualization, Writing - Review and Editing, Supervision, Project administration, Funding acquisition.
Declaration of competing interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
This work was supported in part by the American Heart Association (Award No. 18EIA33900046) to SCS. Additionally, this research has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska–Curie action H2020-MSCA-IF-2016 InsiliCardio, GA No. 750835 and under the ERA-NET co-fund action No. 680969 (ERA-CVD SICVALVES) funded by the Austrian Science Fund (FWF), Grant I 4652-B to CMA.
References (71)
- et al.
A computational method to assess the in vivo stresses and unloaded configuration of patient-specific blood vessels
J. Comput. Appl. Math.
(2013) - et al.
An alternate formulation of blood vessel mechanics and the meaning of the in vivo property
J. Biomech.
(1995) - et al.
Methods to compute 3D residual stress distributions in hyperelastic tubes with application to arterial walls
Int. J. Eng. Sci.
(2010) - et al.
Stress-modulated collagen fiber remodeling in a human carotid bifurcation
J. Theor. Biol.
(2007) - et al.
On non-physical response in models for fiber-reinforced hyperelastic materials
Int. J. Solid Struct.
(2010) Determination of material models for arterial walls from uniaxial extension tests and histological structure
J. Theor. Biol.
(2006)- et al.
On the tension-compression switch in soft fibrous solids
Eur. J. Mech. Solid.
(2015) - et al.
Constrained mixture modeling affects material parameter identification from planar biaxial tests
J. Mech. Behav. Biomed. Mater.
(2019) - et al.
A method for incorporating three-dimensional residual stretches/stresses into patient-specific finite element simulations of arteries
J. Mech. Behav. Biomed. Mater.
(2015) - et al.
Stress-dependent finite growth in soft elastic tissues
J. Biomech.
(1994)
On the physical assumptions underlying the volumetric-isochoric split and the case of anisotropy
Eur. J. Mech. Solid.
Strain energy density function and uniform strain hypothesis for arterial mechanics
J. Biomech.
Experimental measurements of elastic properties of media and adventitia of bovine carotid arteries
J. Biomech.
Medical image-based simulation of abdominal aortic aneurysm growth
Mech. Res. Commun.
Numerical framework for patient-specific computational modelling of vascular tissue
Int. J. Numer. Methods Biomed. Eng.
Assessing the use of the “opening angle method” to enforce residual stresses in patient-specific arteries
Ann. Biomed. Eng.
The vascular modeling toolkit
Classical and all-floating feti methods for the simulation of arterial tissues
Int. J. Numer. Methods Eng.
Competition between radial expansion and thickening in the enlargement of an intracranial saccular aneurysm
J. Elasticity
A theoretical model of enlarging intracranial fusiform aneurysms
J. Biomech. Eng.
A microstructurally motivated model of arterial wall mechanics with mechanobiological implications
Ann. Biomed. Eng.
Three-dimensional stress distribution in arteries
J. Biomech. Eng.
Residual stress in arteries
A homogenized constrained mixture (and mechanical analog) model for growth and remodeling of soft tissue
Biomech. Model. Mechanobiol.
Preferred fiber orientations in healthy arteries and veins understood from netting analysis
Math. Mech. Solid
Growth and remodeling of load-bearing biological soft tissues
Meccanica
Computational Modeling in Biomechanics
Elastin as a rubber
Biopolymers
The simple tension problem at large volumetric strains computed from finite hyperelastic material laws
Acta Mech.
Biomechanical phenotyping of central arteries in health and disease: advantages of and methods for murine models
Ann. Biomed. Eng.
A computational framework for fluid–solid-growth modeling in cardiovascular simulations
Comput. Methods Appl. Mech. Eng.
A computational framework for fluid–solid-growth modeling in cardiovascular simulations
Comput. Methods Appl. Mech. Eng.
NIH public access
Methods
Three-dimensional collagen organization of human brain arteries at different transmural pressures
J. Vasc. Res.
Thermodynamic relations for high elastic materials
Trans. Faraday Soc.
Cited by (4)
Drivers of vascular growth and remodeling: A computational framework to promote benign adaptation in the Ross procedure
2023, Journal of the Mechanical Behavior of Biomedical MaterialsRobust and efficient fixed-point algorithm for the inverse elastostatic problem to identify myocardial passive material parameters and the unloaded reference configuration
2022, Journal of Computational PhysicsCitation Excerpt :Heterogeneity in fiber stretch in an end-diastolic state impairs the Frank-Starling mechanism, as shown in a recent modeling study [8]. In future studies, the presented workflow will be extended by using ideas based on growth and remodeling [27,88] to address this issue. Third, we presented the methodology only for single-chamber LV models, while solid- and electromechanical whole-organ simulations of the heart are becoming feasible [9,11,72,78].
Fibrillar Collagen: A Review of the Mechanical Modeling of Strain-Mediated Enzymatic Turnover
2021, Applied Mechanics Reviews
- 1
Authors contributed equally.