A finite element model of stress-mediated vascular adaptation: application to abdominal aortic aneurysms

Despite rapid expansion of our knowledge of vascular adaptation, developing patient-specific models of diseased arteries is still an open problem. In this study, we extend existing finite element models of stress-mediated growth and remodelling of arteries to incorporate a medical image-based geometry of a healthy aorta and, then, simulate abdominal aortic aneurysm. Degradation of elastin initiates a local dilatation of the aorta while stress-mediated turnover of collagen and smooth muscle compensates the loss of elastin. Stress distributions and expansion rates during the aneurysm growth are studied for multiple spatial distribution functions of elastin degradation and kinetic parameters. Temporal variations of the degradation function are also investigated with either direct time-dependent degradation or stretch-induced degradation as possible biochemical and biomechanical mechanisms for elastin degradation. The results show that this computational model has the capability to capture the complexities of aneurysm progression due to variations of geometry, extent of damage and stress-mediated turnover as a step towards patient-specific modelling.     

Progressive structural and biomechanical changes in elastin degraded aorta

Aortic aneurysm is an important clinical condition characterized by common structural changes such as the degradation of elastin, loss of smooth muscle cells, and increased deposition of fibrillary collagen. With the goal of investigating the relationship between the mechanical behavior and the structural/biochemical composition of an artery, this study used a simple chemical degradation model of aneurysm and investigated the progressive changes in mechanical properties. Porcine thoracic aortas were digested in a mild solution of purified elastase (5 U/mL) for 6, 12, 24, 48, and 96 h. Initial size measurements show that disruption of the elastin structure leads to increased artery dilation in the absence of periodic loading. The mechanical properties of the digested arteries, measured with a biaxial tensile testing device, progress through four distinct stages termed (1) initial-softening, (2) elastomer-like, (3) extensible-but-stiff, and (4) collagen-scaffold-like. While stages 1, 3, and 4 are expected as a result of elastin degradation, the S-shaped stress versus strain behavior of the aorta resulting from enzyme digestion has not been reported previously. Our results suggest that gradual changes in the structure of elastin in the artery can lead to a progression through different mechanical properties and thus reveal the potential existence of an important transition stage that could contribute to artery dilation during aneurysm formation.     

Modelling evolution and the evolving mechanical environment of saccular cerebral aneurysms

A fluid–solid-growth (FSG) model of saccular cerebral aneurysm evolution is developed. It utilises a realistic two-layered structural model of the internal carotid artery and explicitly accounts for the degradation of the elastinous constituents and growth and remodelling (G&R) of the collagen fabric. Aneurysm inception is prescribed: a localised degradation of elastin results in a perturbation in the arterial geometry; the collagen fabric adapts, and the artery achieves a new homeostatic configuration. The perturbation to the geometry creates an altered haemodynamic environment. Subsequent degradation of elastin is explicitly linked to low wall shear stress (WSS) in a confined region of the arterial domain. A sidewall saccular aneurysm develops, the collagen fabric adapts and the aneurysm stabilises in size. A quasi-static analysis is performed to determine the geometry at diastolic pressure. This enables the cyclic stretching of the tissue to be quantified, and we propose a novel index to quantify the degree of biaxial stretching of the tissue. Whilst growth is linked to low WSS from a steady (systolic) flow analysis, a pulsatile flow analysis is performed to compare steady and pulsatile flow parameters during evolution. This model illustrates the evolving mechanical environment for an idealised saccular cerebral aneurysm developing on a cylindrical parent artery and provides the guidance to more sophisticated FSG models of aneurysm evolution which link G&R to the local mechanical stimuli of vascular cells.     

Growth and remodeling in a thick-walled artery model: effects of spatial variations in wall constituents

A mathematical model is presented for growth and remodeling of arteries. The model is a thick-walled tube composed of a constrained mixture of smooth muscle cells, elastin and collagen. Material properties and radial and axial distributions of each constituent are prescribed according to previously published data. The analysis includes stress-dependent growth and contractility of the muscle and turnover of collagen fibers. Simulations were conducted for homeostatic conditions and for the temporal response following sudden hypertension. Numerical pressure-radius relations and opening angles (residual stress) show reasonable agreement with published experimental results. In particular, for realistic material and structural properties, the model predicts measured variations in opening angles along the length of the aorta with reasonable accuracy. These results provide a better understanding of the determinants of residual stress in arteries and could lend insight into the importance of constituent distributions in both natural and tissue-engineered blood vessels.     

3D geometric reconstruction of thoracic aortic aneurysms

Background  

The thoracic aortic aneurysm (TAA) is a pathology that involves an expansion of the aortic diameter in the thoracic aorta, leading to risk of rupture. Recent studies have suggested that internal wall stress, which is affected by TAA geometry and the presence or absence of thrombus, is a more reliable predictor of rupture than the maximum diameter, the current clinical criterion. Accurate reconstruction of TAA geometry is a crucial step in patient-specific stress calculations.     

Aorta elastin turnover in normal and hypercholesterolemic Japanese quail

The turnover and degradation of mature elastin from the aortae of Japanese quail were estimated following with l-[U-¹⁴C]lysine by measuring the changes in specific activity of l-[U-¹⁴C]lysine and ¹⁴C-labelled desmosine and isodesmosine (crosslinking amino acids derived from lysyl residues) in elastin over a 39-week period. Only 5% of the variation in radioactivity could be attributed to changes in time. Therefore, it was concluded that the best estimates of mature elastin turnover are only quantifiable in years. Dietary cholesterol in amounts sifficient to induce plaque formation and fragmentation of the elastic lamina in the aorta did not significantly influence turnover time. It would appear that once the total pool of elastin in aorta is stabilized as mature fibers it is not subject to proteolysis or resynthesis of sufficient magnitude to result in measurable turnover.     

Validation of CFD simulations of cerebral aneurysms with implication of geometric variations

BACKGROUND: Computational fluid dynamics (CFD) simulations using medical-image-based anatomical vascular geometry are now gaining clinical relevance. This study aimed at validating the CFD methodology for studying cerebral aneurysms by using particle image velocimetry (PIV) measurements, with a focus on the effects of small geometric variations in aneurysm models on the flow dynamics obtained with CFD. METHOD OF APPROACH: An experimental phantom was fabricated out of silicone elastomer to best mimic a spherical aneurysm model. PIV measurements were obtained from the phantom and compared with the CFD results from an ideal spherical aneurysm model (S1). These measurements were also compared with CFD results, based on the geometry reconstructed from three-dimensional images of the experimental phantom. We further performed CFD analysis on two geometric variations, S2 and S3, of the phantom to investigate the effects of small geometric variations on the aneurysmal flow field. Results. We found poor agreement between the CFD results from the ideal spherical aneurysm model and the PIV measurements from the phantom, including inconsistent secondary flow patterns. The CFD results based on the actual phantom geometry, however, matched well with the PIV measurements. CFD of models S2 and S3 produced qualitatively similar flow fields to that of the phantom but quantitatively significant changes in key hemodynamic parameters such as vorticity, positive circulation, and wall shear stress. CONCLUSION: CFD simulation results can closely match experimental measurements as long as both are performed on the same model geometry. Small geometric variations on the aneurysm model can significantly alter the flow-field and key hemodynamic parameters. Since medical images are subjected to geometric uncertainties, image-based patient-specific CFD results must be carefully scrutinized before providing clinical feedback.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

3D reconstruction and manufacture of real abdominal aortic aneurysms: from CT scan to silicone model

Abdominal aortic aneurysm (AAA) can be defined as a permanent and irreversible dilation of the infrarenal aorta. AAAs are often considered to be an aorta with a diameter 1.5 times the normal infrarenal aorta diameter. This paper describes a technique to manufacture realistic silicone AAA models for use with experimental studies. This paper is concerned with the reconstruction and manufacturing process of patient-specific AAAs. 3D reconstruction from computed tomography scan data allows the AAA to be created. Mould sets are then designed for these AAA models utilizing computer aided designcomputer aided manufacture techniques and combined with the injection-moulding method. Silicone rubber forms the basis of the resulting AAA model. Assessment of wall thickness and overall percentage difference from the final silicone model to that of the computer-generated model was performed. In these realistic AAA models, wall thickness was found to vary by an average of 9.21%. The percentage difference in wall thickness recorded can be attributed to the contraction of the casting wax and the expansion of the silicone during model manufacture. This method may be used in conjunction with wall stress studies using the photoelastic method or in fluid dynamic studies using a laser-Doppler anemometry. In conclusion, these patient-specific rubber AAA models can be used in experimental investigations, but should be assessed for wall thickness variability once manufactured.     

Intrasac pressure changes and vascular remodeling after endovascular repair of abdominal aortic aneurysms: review and biomechanical model simulation

In this paper, we review existing clinical research data on post-endovascular repair (EVAR) intrasac pressure and relation with abdominal aortic aneurysm (AAA) size changes. Based on the review, we hypothesize that intrasac pressure has a significant impact on post-EVAR AAA size changes, and post-EVAR remodeling depends also on how the pressure has changed over a period of time. The previously developed model of an AAA based on a constrained mixture approach is extended to include vascular adaptation after EVAR using an idealized geometry. Computational simulation shows that the same mechanism of collagen stress-mediated remodeling in AAA expansion induces the aneurysm wall to shrink in a reduced sac-pressure after post-EVAR. Computational simulation suggests that the intrasac pressure of 60 mm?Hg is a critical value. At this value, the AAA remains stable, while values above cause the AAA to expand and values below cause the AAA to shrink. There are, however, variations between individuals due to different cellular sensitivities in stress-mediated adaptation. Computer simulation also indicates that an initial decrease in intrasac pressure helps the AAA shrink even if the pressure increases after some time. The presented study suggests that biomechanics has a major effect on initial adaptation after EVAR and also illustrates the utility of a computational model of vascular growth and remodeling in predicting diameter changes during the progression and after the treatment of AAAs.     

Computational study of growth and remodelling in the aortic arch

Opening angles (OAs) are associated with growth and remodelling in arteries. One curiosity has been the relatively large OAs found in the aortic arch of some animals. Here, we use computational models to explore the reasons behind this phenomenon. The artery is assumed to contain a smooth muscle/collagen phase and an elastin phase. In the models, growth and remodelling of smooth muscle/collagen depends on wall stress and fluid shear stress. Remodelling of elastin, which normally turns over very slowly, is neglected. The results indicate that OAs generally increase with longitudinal curvature (torus model), earlier elastin production during development, and decreased wall stiffness. Correlating these results with available experimental data suggests that all of these effects may contribute to the large OAs in the aortic arch. The models also suggest that the slow turnover rate of elastin limits longitudinal growth. These results should promote increased understanding of the causes of residual stress in arteries.     

A volumetric model for growth of arterial walls with arbitrary geometry and loads

Stress and deformation in arterial wall tissue are factors which may influence significantly its response and evolution. In this work we develop models based on nonlinear elasticity and finite element numerical solutions for the mechanical behaviour and the remodelling of the soft tissue of arteries, including anisotropy induced by the presence of collagen fibres. Remodelling and growth in particular constitute important features in order to interpret stenosis and atherosclerosis. The main object of this work is to model accurately volumetric growth, induced by fluid shear stress in the intima and local wall stress in arteries with patient-specific geometry and loads. The model is implemented in a nonlinear finite element setting which may be applied to realistic 3D geometries obtained from in vivo measurements. The capabilities of this method are demonstrated in several examples. Firstly a stenotic process on an idealised geometry induced by a non-uniform shear stress distribution is considered. Following the growth of a right coronary artery from an in vivo reconstructed geometry is presented. Finally, experimental measurements for growth under hypertension for rat carotid arteries are modelled.     

Identification of in vivo material and geometric parameters of a human aorta: toward patient-specific modeling of abdominal aortic aneurysm

Recent advances in computational modeling of vascular adaptations and the need for their extension to patient-specific modeling have introduced new challenges to the path toward abdominal aortic aneurysm modeling. First, the fundamental assumption in adaptation models, namely the existence of vascular homeostasis in normal vessels, is not easy to implement in a vessel model built from medical images. Second, subjecting the vessel wall model to the normal pressure often makes the configuration deviate from the original geometry obtained from medical images. To address those technical challenges, in this work, we propose a two-step optimization approach; first, we estimate constitutive parameters of a healthy human aorta intrinsic to the material by using biaxial test data and a weighted nonlinear least-squares parameter estimation method; second, we estimate the distributions of wall thickness and anisotropy using a 2-D parameterization of the vessel wall surface and a global approximation scheme integrated within an optimization routine. A direct search method is implemented to solve the optimization problem. The numerical optimization method results in a considerable improvement in both satisfying homeostatic condition and minimizing the deviation of geometry from the original shape based on in vivo images. Finally, the utility of the proposed technique for patient-specific modeling is demonstrated in a simulation of an abdominal aortic aneurysm enlargement.     

The canine tail artery as a model for cerebral aneurysm studies

The occluded canine tail artery, which comes off in the same plane as the aortoiliac junction, has been used as a flow model for cerebral aneurysms. These experiments were designed to determine if it is a realistic distensible model of human intracranial aneurysms. Distensibility studies were done on the aorta, and the iliac and tail arteries of four dogs. From these pressure-volume studies, tension-strain curves, elastances, and collagen slack were obtained. The tail artery is stiffer longitudinally and more distensible circumferentially than the other vessels. The iliac arteries and the aorta are not significantly different. The elastance of elastin and collagen is lower in the tail artery, and the collagen is more wavy circumferentially. Longitudinally, the collagen slack is least for the tail artery, and the elastance of elastin is not different in all three vessels. The number of elastin layers in the iliac and tail arteries seen in cross section is not significantly different, but the aorta is different from both these vessels. In another four dogs the aorta proximal to the trifurcation was cannulated and infused with saline to increase pressure. India ink marks were put on the surface to measure changes in length. Photographs were taken at intervals of 10 mmHg(1 mmHg = 133.3 Pa). This was done with the vessels tethered and untethered in the body and then taken out and studied with the same method in vitro. Arteries tethered in the body expanded circumferentially more than longitudinally. The tail artery becomes less distensible if untethered in the body and therefore acts more like an aneurysm.(ABSTRACT TRUNCATED AT 250 WORDS)     

Linked mechanical and biological aspects of remodeling in mouse pulmonary arteries with hypoxia-induced hypertension

Supravalvular aortic stenosis (SVAS) is associated with decreased elastin and altered arterial mechanics. Mice with a single deletion in the elastin gene (ELN(+/-)) are models for SVAS. Previous studies have shown that elastin haploinsufficiency in these mice causes hypertension, decreased arterial compliance, and changes in arterial wall structure. Despite these differences, ELN(+/-) mice have a normal life span, suggesting that the arteries remodel and adapt to the decreased amount of elastin. To test this hypothesis, we performed in vitro mechanical tests on abdominal aorta, ascending aorta, and left common carotid artery from ELN(+/-) and wild-type (C57BL/6J) mice. We compared the circumferential and longitudinal stress-stretch relationships and residual strains. The circumferential stress-stretch relationship is similar between genotypes and changes <3% with longitudinal stretch at lengths within 10% of the in vivo value. At mean arterial pressure, the circumferential stress in the ascending aorta is higher in ELN(+/-) than in wild type. Although arterial pressures are higher, the increased number of elastic lamellae in ELN(+/-) arteries results in similar tension/lamellae compared with wild type. The longitudinal stress-stretch relationship is similar between genotypes for most arteries. Compared with wild type, the in vivo longitudinal stretch is lower in ELN(+/-) abdominal and carotid arteries and the circumferential residual strain is higher in ELN(+/-) ascending aorta. The increased circumferential residual strain brings the transmural strain distribution in ELN(+/-) ascending aorta close to wild-type values. The mechanical behavior of ELN(+/-) arteries is likely due to the reduced elastin content combined with adaptive remodeling during vascular development.     

Mechanics, mechanobiology, and modeling of human abdominal aorta and aneurysms

Biomechanical factors play fundamental roles in the natural history of abdominal aortic aneurysms (AAAs) and their responses to treatment. Advances during the past two decades have increased our understanding of the mechanics and biology of the human abdominal aorta and AAAs, yet there remains a pressing need for considerable new data and resulting patient-specific computational models that can better describe the current status of a lesion and better predict the evolution of lesion geometry, composition, and material properties and thereby improve interventional planning. In this paper, we briefly review data on the structure and function of the human abdominal aorta and aneurysmal wall, past models of the mechanics, and recent growth and remodeling models. We conclude by identifying open problems that we hope will motivate studies to improve our computational modeling and thus general understanding of AAAs.     

Role of elastin anisotropy in structural strain energy functions of arterial tissue

The vascular wall exhibits nonlinear anisotropic mechanical properties. The identification of a strain energy function (SEF) is the preferred method to describe its complex nonlinear elastic properties. Earlier constituent-based SEF models, where elastin is modeled as an isotropic material, failed in describing accurately the tissue response to inflation–extension loading. We hypothesized that these shortcomings are partly due to unaccounted anisotropic properties of elastin. We performed inflation–extension tests on common carotid of rabbits before and after enzymatic degradation of elastin and applied constituent-based SEFs, with both an isotropic and an anisotropic elastin part, on the experimental data. We used transmission electron microscopy (TEM) and serial block-face scanning electron microscopy (SBFSEM) to provide direct structural evidence of the assumed anisotropy. In intact arteries, the SEF including anisotropic elastin with one family of fibers in the circumferential direction fitted better the inflation–extension data than the isotropic SEF. This was supported by TEM and SBFSEM imaging, which showed interlamellar elastin fibers in the circumferential direction. In elastin-degraded arteries, both SEFs succeeded equally well in predicting anisotropic wall behavior. In elastase-treated arteries fitted with the anisotropic SEF for elastin, collagen engaged later than in intact arteries. We conclude that constituent-based models with an anisotropic elastin part characterize more accurately the mechanical properties of the arterial wall when compared to models with simply an isotropic elastin. Microstructural imaging based on electron microscopy techniques provided evidence for elastin anisotropy. Finally, the model suggests a later and less abrupt collagen engagement after elastase treatment.     

Nonlinear anisotropic stress analysis of anatomically realistic cerebral aneurysms

BACKGROUND: Static deformation analysis and estimation of wall stress distribution of patient-specific cerebral aneurysms can provide useful insights into the disease process and rupture. METHOD OF APPROACH: The three-dimensional geometry of saccular cerebral aneurysms from 27 patients (18 unruptured and nine ruptured) was reconstructed based on computer tomography angiography images. The aneurysm wall tissue was modeled using a nonlinear, anisotropic, hyperelastic material model (Fung-type) which was incorporated in a user subroutine in ABAQUS. Effective material fiber orientations were assumed to align with principal surface curvatures. Static deformation of the aneurysm models were simulated assuming uniform wall thickness and internal pressure load of 100 mm Hg. RESULTS: The numerical analysis technique was validated by quantitative comparisons to results in the literature. For the patient-specific models, in-plane stresses in the aneurysm wall along both the stiff and weak fiber directions showed significant regional variations with the former being higher. The spatial maximum of stress ranged from as low as 0.30 MPa in a small aneurysm to as high as 1.06 MPa in a giant aneurysm. The patterns of distribution of stress, strain, and surface curvature were found to be similar. Sensitivity analyses showed that the computed stress is mesh independent and not very sensitive to reasonable perturbations in model parameters, and the curvature-based criteria for fiber orientations tend to minimize the total elastic strain energy in the aneurysms wall. Within this small study population, there were no statistically significant differences in the spatial means and maximums of stress and strain values between the ruptured and unruptured groups. However, the ratios between the stress components in the stiff and weak fiber directions were significantly higher in the ruptured group than those in the unruptured group. CONCLUSIONS: A methodology for nonlinear, anisotropic static deformation analysis of geometrically realistic aneurysms was developed, which can be used for a more accurate estimation of the stresses and strains than previous methods and to facilitate prospective studies on the role of stress in aneurysm rupture.     

Elastin as a self-organizing biomaterial: use of recombinantly expressed human elastin polypeptides as a model for investigations of structure and self-assembly of elastin

Elastin is the major extracellular matrix protein of large arteries such as the aorta, imparting characteristics of extensibility and elastic recoil. Once laid down in tissues, polymeric elastin is not subject to turnover, but is able to sustain its mechanical resilience through thousands of millions of cycles of extension and recoil. Elastin consists of ca. 36 domains with alternating hydrophobic and cross-linking characteristics. It has been suggested that these hydrophobic domains, predominantly containing glycine, proline, leucine and valine, often occurring in tandemly repeated sequences, are responsible for the ability of elastin to align monomeric chains for covalent cross-linking. We have shown that small, recombinantly expressed polypeptides based on sequences of human elastin contain sufficient information to self-organize into fibrillar structures and promote the formation of lysine-derived cross-links. These cross-linked polypeptides can also be fabricated into membrane structures that have solubility and mechanical properties reminiscent of native insoluble elastin. Understanding the basis of the self-organizational ability of elastin-based polypeptides may provide important clues for the general design of self-assembling biomaterials.     

Improving the Efficiency of Abdominal Aortic Aneurysm Wall Stress Computations

An abdominal aortic aneurysm is a pathological dilation of the abdominal aorta, which carries a high mortality rate if ruptured. The most commonly used surrogate marker of rupture risk is the maximal transverse diameter of the aneurysm. More recent studies suggest that wall stress from models of patient-specific aneurysm geometries extracted, for instance, from computed tomography images may be a more accurate predictor of rupture risk and an important factor in AAA size progression. However, quantification of wall stress is typically computationally intensive and time-consuming, mainly due to the nonlinear mechanical behavior of the abdominal aortic aneurysm walls. These difficulties have limited the potential of computational models in clinical practice. To facilitate computation of wall stresses, we propose to use a linear approach that ensures equilibrium of wall stresses in the aneurysms. This proposed linear model approach is easy to implement and eliminates the burden of nonlinear computations. To assess the accuracy of our proposed approach to compute wall stresses, results from idealized and patient-specific model simulations were compared to those obtained using conventional approaches and to those of a hypothetical, reference abdominal aortic aneurysm model. For the reference model, wall mechanical properties and the initial unloaded and unstressed configuration were assumed to be known, and the resulting wall stresses were used as reference for comparison. Our proposed linear approach accurately approximates wall stresses for varying model geometries and wall material properties. Our findings suggest that the proposed linear approach could be used as an effective, efficient, easy-to-use clinical tool to estimate patient-specific wall stresses.