A model for saccular cerebral aneurysm growth by collagen fibre remodelling

Quantifying asymmetry of anterior cerebral arteries as a predictor of anterior communicating artery complex aneurysm

OBJECTIVES

The aim of this study was to establish an anatomical index for early prediction of the risk of development of aneurysms in anterior communicating arterial complex (AcomAC). The asymmetric diameter of one anterior cerebral artery (ACA) to other could alter haemodynamics and may contribute to formation of aneurysms in AcomAC and be a reliable predictor of the risk of development of aneurysms.

DESIGN AND SETTING

This is a retrospective, observational and quantitative study, which used cerebral computed tomography angiography (CCTA) scans in South Australia.

PARTICIPANTS

CCTA scans of 166 adult patients of both sexes were studied.

MAIN OUTCOME MEASURES

The internal diameters of the proximal segments of ACAs (A1s) were measured. Position and presence or absence of aneurysms in AcomAC were determined. The ratio of A1 diameters was taken as a measure of A1 asymmetry.

RESULTS

The ratio of diameters of A1s correlated with the occurrence of AcomAC aneurysms. The risk of development of aneurysms in AcomAC was much greater (80%, OR=47.3) when one A1 segment's radius was at least 50% larger (ie, 2.25 times cross-sectional area) than the other.

CONCLUSION

The general information on asymmetric A1 has been published previously. The present findings have significant contribution since the A1s asymmetry ratios have been categorised in ascending order and matched with the presence of AcomAC aneurysms. The asymmetry ratio of the A1 is a good predictor for the development of AcomAC aneurysms. Reconstruction of the asymmetric A1 could be done if the technology gets advanced.

Epicardial prestrained confinement and residual stresses: a newly observed heart ventricle confinement interface

The heart epicardial layer, with elastin as the dominant component, has not been well investigated, specifically on how it contributes to ventricular biomechanics. In this study, we revealed and quantitatively assessed the overall status of prestraining and residual stresses exerted by the epicardial layer on the heart left ventricle (LV). During porcine heart wall dissection, we discovered that bi-layered LV surface strips, consisting of an epicardial layer and cardiac muscle, always curled towards the epicardial side due to epicardial residual stresses. We hence developed a curling angle characterization technique to intuitively and qualitatively reveal the location-dependency and direction-dependency of epicardial residual stresses. Moreover, by combining prestrain measurement and biaxial mechanical testing, we were able to quantify the epicardial prestrains and residual stresses on the unpressurized intact LV. To investigate the potential mechanical effect of epicardial prestraining, a finite-element (FE) model has been constructed, and we demonstrate that it is the prestraining of the epicardial layer, not the epicardial layer alone, providing an additional resistance mechanism during LV diastolic expansion and ventricular wall protection by reducing myocardial stress. In short, our study on healthy, native porcine hearts has revealed an important phenomenon-the epicardial layer, rich in elastin, acts like a prestrained 'balloon' that wraps around the heart and functions as an extra confinement and protection interface. The obtained knowledge fills a gap in ventricular biomechanics and will help design novel biomimicking materials or prosthetic devices to target the maintenance/recreation of this ventricle confinement interface.

Computational modeling of cardiac growth and remodeling in pressure overloaded hearts-Linking microstructure to organ phenotype

Cardiac growth and remodeling (G&R) refers to structural changes in myocardial tissue in response to chronic alterations in loading conditions. One such condition is pressure overload where elevated wall stresses stimulate the growth in cardiomyocyte thickness, associated with a phenotype of concentric hypertrophy at the organ scale, and promote fibrosis. The initial hypertrophic response can be considered adaptive and beneficial by favoring myocyte survival, but over time if pressure overload conditions persist, maladaptive mechanisms favoring cell death and fibrosis start to dominate, ultimately mediating the transition towards an overt heart failure phenotype. The underlying mechanisms linking biological factors at the myocyte level to biomechanical factors at the systemic and organ level remain poorly understood. Computational models of G&R show high promise as a unique framework for providing a quantitative link between myocardial stresses and strains at the organ scale to biological regulatory processes at the cellular level which govern the hypertrophic response. However, microstructurally motivated, rigorously validated computational models of G&R are still in their infancy. This article provides an overview of the current state-of-the-art of computational models to study cardiac G&R. The microstructure and mechanosensing/mechanotransduction within cells of the myocardium is discussed and quantitative data from previous experimental and clinical studies is summarized. We conclude with a discussion of major challenges and possible directions of future research that can advance the current state of cardiac G&R computational modeling. STATEMENT OF SIGNIFICANCE: The mechanistic links between organ-scale biomechanics and biological factors at the cellular size scale remain poorly understood as these are largely elusive to investigations using experimental methodology alone. Computational G&R models show high promise to establish quantitative links which allow more mechanistic insight into adaptation mechanisms and may be used as a tool for stratifying the state and predict the progression of disease in the clinic. This review provides a comprehensive overview of research in this domain including a summary of experimental data. Thus, this study may serve as a basis for the further development of more advanced G&R models which are suitable for making clinical predictions on disease progression or for testing hypotheses on pathogenic mechanisms using in-silico models.

Stress-Swelling Finite Element Modeling of Cervical Response With Homeostatic Collagen Fiber Distributions

During pregnancy, the cervix experiences significant mechanical property change due to tissue swelling, and to ongoing changes in the collagen content. In this paper, we model how these two effects contribute to cervical deformation as the pressure load on top of the cervix increases. The cervix and its surrounding supporting ligaments are taken into consideration in the resulting mechanical analysis. The cervix itself is treated as a multilayered tube-like structure, with layer-specific collagen orientation. The cervical tissue in each layer is treated in terms of a collagen constituent that remodels with time within a ground substance matrix that experiences swelling. The load and swelling are taken to change sufficiently slowly so that the collagen properties at any instant can be regarded as being in a state of homeostasis. Among other things, the simulations show how the luminal cross-sectional area varies along its length as a function of pressure and swelling. In general, an increase in pressure causes an overall shortening of the lumen while an increase in swelling has the opposite effect.

Integration of polarized spatial frequency domain imaging (pSFDI) with a biaxial mechanical testing system for quantification of load-dependent collagen architecture in soft collagenous tissues

Collagen fiber networks provide the structural strength of tissues, such as tendons, skin and arteries. Quantifying the fiber architecture in response to mechanical loads is essential towards a better understanding of the tissue-level mechanical behaviors, especially in assessing disease-driven functional changes. To enable novel investigations into these load-dependent fiber structures, a polarized spatial frequency domain imaging (pSFDI) device was developed and, for the first time, integrated with a biaxial mechanical testing system. The integrated instrument is capable of a wide-field quantification of the fiber orientation and the degree of optical anisotropy (DOA), representing the local degree of fiber alignment. The opto-mechanical instrument''s performance was assessed through uniaxial loading on tendon tissues with known collagen fiber microstructures. Our results revealed that the bulk fiber orientation angle of the tendon tissue changed minimally with loading (median ± 0.5*IQR of 52.7° ± 3.3° and 51.9° ± 3.3° under 0 and 3% longitudinal strains, respectively), whereas on a micro-scale, the fibers became better aligned with the direction of loading: the DOA (mean ± SD) increased from 0.149 ± 0.032 to 0.198 ± 0.056 under 0 and 3% longitudinal strains, respectively, p < 0.001. The integrated instrument was further applied to study two representative mitral valve anterior leaflet (MVAL) tissues subjected to various biaxial loads. The fiber orientations within these representative MVAL tissue specimens demonstrated noticeable heterogeneity, with the local fiber orientations dependent upon the sample, the spatial and transmural locations, and the applied loading. Our results also showed that fibers were generally better aligned under equibiaxial (DOA = 0.089 ± 0.036) and circumferentially-dominant loading (DOA = 0.086 ± 0.037) than under the radially-dominant loading (DOA = 0.077 ± 0.034), indicating circumferential predisposition. These novel findings exemplify a deeper understanding of the load-dependent collagen fiber microstructures obtained through the use of the integrated opto-mechanical instrument. STATEMENT OF SIGNIFICANCE: In this study, a novel quantitative opto-mechanical system was developed by combining a polarized Spatial Frequency Domain Imaging (pSFDI) device with a biaxial mechanical tester. The integrated system was used to quantify the load-dependent collagen fiber microstructures in representative tendon and mitral valve anterior leaflet (MVAL) tissues. Our results revealed that MVAL's fiber architectures exhibited load-dependent spatial and transmural heterogeneities, suggesting further microstructural complexity than previously reported in heart valve tissues. These novel findings were possible through the system's ability to, for the first time, capture the load-dependent collagen architecture in the mitral valve anterior leaflet tissue over a wide field of view (e.g., 10 × 10 mm for the MVAL tissue specimens). Such capabilities afford unique future opportunities to improve patient outcomes through concurrent mechanical and microstructural assessments of healthy and diseased tissues in conditions such as heart valve regurgitation and calcification.

Prior Distributions of Material Parameters for Bayesian Calibration of Growth and Remodeling Computational Model of Abdominal Aortic Wall

For the accurate prediction of the vascular disease progression, there is a crucial need for developing a systematic tool aimed toward patient-specific modeling. Considering the interpatient variations, a prior distribution of model parameters has a strong influence on computational results for arterial mechanics. One crucial step toward patient-specific computational modeling is to identify parameters of prior distributions that reflect existing knowledge. In this paper, we present a new systematic method to estimate the prior distribution for the parameters of a constrained mixture model using previous biaxial tests of healthy abdominal aortas (AAs). We investigate the correlation between the estimated parameters for each constituent and the patient's age and gender; however, the results indicate that the parameters are correlated with age only. The parameters are classified into two groups: Group-I in which the parameters ce, ck1, ck2, cm2,Ghc, and ϕe are correlated with age, and Group-II in which the parameters cm1, Ghm, G1e, G2e, and α are not correlated with age. For the parameters in Group-I, we used regression associated with age via linear or inverse relations, in which their prior distributions provide conditional distributions with confidence intervals. For Group-II, the parameter estimated values were subjected to multiple transformations and chosen if the transformed data had a better fit to the normal distribution than the original. This information improves the prior distribution of a subject-specific model by specifying parameters that are correlated with age and their transformed distributions. Therefore, this study is a necessary first step in our group's approach toward a Bayesian calibration of an aortic model. The results from this study will be used as the prior information necessary for the initialization of Bayesian calibration of a computational model for future applications.

Aneurysm strength can decrease under calcification

Aneurysms are abnormal dilatations of vessels in the vascular system that are prone to rupture. Prediction of the aneurysm rupture is a challenging and unsolved problem. Various factors can lead to the aneurysm rupture and, in the present study, we examine the effect of calcification on the aneurysm strength by using micromechanical modeling. The calcified tissue is considered as a composite material in which hard calcium particles are embedded in a hyperelastic soft matrix. Three experimentally calibrated constitutive models incorporating a failure description are used for the matrix representation. Two constitutive models describe the aneurysmal arterial wall and the third one - the intraluminal thrombus. The stiffness and strength of the calcified tissue are simulated in uniaxial tension under the varying amount of calcification, i.e. the relative volume of the hard inclusion within the periodic unit cell. In addition, the triaxiality of the stress state, which can be a trigger for the cavitation instability, is tracked. Results of the micromechanical simulation show an increase of the stiffness and a possible decrease of the strength of the calcified tissue as compared to the non-calcified one. The obtained results suggest that calcification (i.e. the presence of hard particles) can significantly affect the stiffness and strength of soft tissue. The development of refined experimental techniques that will allow for the accurate quantitative assessment of calcification is desirable.

Computational Growth and Remodeling of Abdominal Aortic Aneurysms Constrained by the Spine

Abdominal aortic aneurysms (AAAs) evolve over time, and the vertebral column, which acts as an external barrier, affects their biomechanical properties. Mechanical interaction between AAAs and the spine is believed to alter the geometry, wall stress distribution, and blood flow, although the degree of this interaction may depend on AAAs specific configurations. In this study, we use a growth and remodeling (G&R) model, which is able to trace alterations of the geometry, thus allowing us to computationally investigate the effect of the spine for progression of the AAA. Medical image-based geometry of an aorta is constructed along with the spine surface, which is incorporated into the computational model as a cloud of points. The G&R simulation is initiated by local elastin degradation with different spatial distributions. The AAA-spine interaction is accounted for using a penalty method when the AAA surface meets the spine surface. The simulation results show that, while the radial growth of the AAA wall is prevented on the posterior side due to the spine acting as a constraint, the AAA expands faster on the anterior side, leading to higher curvature and asymmetry in the AAA configuration compared to the simulation excluding the spine. Accordingly, the AAA wall stress increases on the lateral, posterolateral, and the shoulder regions of the anterior side due to the AAA-spine contact. In addition, more collagen is deposited on the regions with a maximum diameter. We show that an image-based computational G&R model not only enhances the prediction of the geometry, wall stress, and strength distributions of AAAs but also provides a framework to account for the interactions between an enlarging AAA and the spine for a better rupture potential assessment and management of AAA patients.

Biomechanical Properties of Human Ascending Thoracic Aortic Dissections

Thoracic aortic dissections are associated with a significant risk of morbidity and mortality, and currently challenge our understanding of the biomechanical factors leading to their initiation and propagation. We quantified the biaxial mechanical properties of human type A dissections (n = 16) and modeled the stress–strain data using a microstructurally motivated form of strain energy function. Our results show significantly higher stiffness for dissected tissues as compared to control aorta without arterial disease. Higher stiffness of dissected tissues did not, however, correlate with greater aortic diameter measured prior to surgery nor were there any age dependent differences in the tissue properties.

Cavitation instability as a trigger of aneurysm rupture

Aneurysm formation and growth is accompanied by microstructural alterations in the arterial wall. Particularly, the loss of elastin may lead to tissue disintegration and appearance of voids or cavities at the micron scale. Unstable growth and coalescence of voids may be a predecessor and trigger for the onset of macroscopic cracks. In the present work, we analyze the instability of membrane (2D) and bulk (3D) voids under hydrostatic tension by using two experimentally calibrated constitutive models of abdominal aortic aneurysm enhanced with energy limiters. The limiters provide the saturation value for the strain energy, which indicates the maximum energy that can be stored and dissipated by an infinitesimal material volume. We find that the unstable growth of voids can start when the critical stress is considerably less than the aneurysm strength. Moreover, this critical stress may even approach the arterial wall stress in the physiological range. This finding suggests that cavitation instability can be a rational indicator of the aneurysm rupture.

Investigation of material modeling in fluid-structure interaction analysis of an idealized three-layered abdominal aorta: aneurysm initiation and fully developed aneurysms

Different material models for an idealized three-layered abdominal aorta are compared using computational techniques to study aneurysm initiation and fully developed aneurysms. The computational model includes fluid-structure interaction (FSI) between the blood vessel and the blood. In order to model aneurysm initiation, the medial region was degenerated to mimic the medial loss occurring in the inception of an aneurysm. Various cases are considered in order to understand their effects on the initiation of an abdominal aortic aneurysm. The layers of the blood vessel were modeled using either linear elastic materials or Mooney-Rivlin (otherwise known as hyperelastic) type materials. The degenerated medial region was also modeled in either linear elastic or hyperelastic-type materials and assumed to be in the shape of an arc with a thin width or a circular ring with different widths. The blood viscosity effect was also considered in the initiation mechanism. In addition, dynamic analysis of the blood vessel was performed without interaction with the blood flow by applying time-dependent pressure inside the lumen in a three-layered abdominal aorta. The stresses, strains, and displacements were compared for a healthy aorta, an initiated aneurysm and a fully developed aneurysm. The study shows that the material modeling of the vessel has a sizable effect on aneurysm initiation and fully developed aneurysms. Different material modeling of degeneration regions also affects the stress-strain response of aneurysm initiation. Additionally, the structural analysis without considering FSI (called noFSI) overestimates the peak von Mises stress by 52% at the interfaces of the layers.

Modeling rupture of growing aneurysms

Growth and rupture of aneurysms are driven by micro-structural alterations of the arterial wall yet precise mechanisms underlying the process remain to be uncovered. In the present work we examine a scenario when the aneurysm evolution is dominated by turnover of collagen fibers. In the latter case it is natural to hypothesize that rupture of individual fibers (or their bonds) causes the overall aneurysm rupture. We examine this hypothesis in computer simulations of growing aneurysms in which constitutive equations describe both collagen evolution and failure. Failure is enforced in constitutive equations by limiting strain energy that can be accumulated in a fiber. Within the proposed theoretical framework we find a range of parameters that lead to the aneurysm rupture. We conclude in a qualitative agreement with clinical observations that some aneurysms will rupture while others will not.

Physical factors effecting cerebral aneurysm pathophysiology

Many factors that are either blood-, wall-, or hemodynamics-borne have been associated with the initiation, growth, and rupture of intracranial aneurysms. The distribution of cerebral aneurysms around the bifurcations of the circle of Willis has provided the impetus for numerous studies trying to link hemodynamic factors (flow impingement, pressure, and/or wall shear stress) to aneurysm pathophysiology. The focus of this review is to provide a broad overview of such hemodynamic associations as well as the subsumed aspects of vascular anatomy and wall structure. Hemodynamic factors seem to be correlated to the distribution of aneurysms on the intracranial arterial tree and complex, slow flow patterns seem to be associated with aneurysm growth and rupture. However, both the prevalence of aneurysms in the general population and the incidence of ruptures in the aneurysm population are extremely low. This suggests that hemodynamic factors and purely mechanical explanations by themselves may serve as necessary, but never as necessary and sufficient conditions of this disease's causation. The ultimate cause is not yet known, but it is likely an additive or multiplicative effect of a handful of biochemical and biomechanical factors.

An efficient framework for optimization and parameter sensitivity analysis in arterial growth and remodeling computations

Computational models for vascular growth and remodeling (G&R) are used to predict the long-term response of vessels to changes in pressure, flow, and other mechanical loading conditions. Accurate predictions of these responses are essential for understanding numerous disease processes. Such models require reliable inputs of numerous parameters, including material properties and growth rates, which are often experimentally derived, and inherently uncertain. While earlier methods have used a brute force approach, systematic uncertainty quantification in G&R models promises to provide much better information. In this work, we introduce an efficient framework for uncertainty quantification and optimal parameter selection, and illustrate it via several examples. First, an adaptive sparse grid stochastic collocation scheme is implemented in an established G&R solver to quantify parameter sensitivities, and near-linear scaling with the number of parameters is demonstrated. This non-intrusive and parallelizable algorithm is compared with standard sampling algorithms such as Monte-Carlo. Second, we determine optimal arterial wall material properties by applying robust optimization. We couple the G&R simulator with an adaptive sparse grid collocation approach and a derivative-free optimization algorithm. We show that an artery can achieve optimal homeostatic conditions over a range of alterations in pressure and flow; robustness of the solution is enforced by including uncertainty in loading conditions in the objective function. We then show that homeostatic intramural and wall shear stress is maintained for a wide range of material properties, though the time it takes to achieve this state varies. We also show that the intramural stress is robust and lies within 5% of its mean value for realistic variability of the material parameters. We observe that prestretch of elastin and collagen are most critical to maintaining homeostasis, while values of the material properties are most critical in determining response time. Finally, we outline several challenges to the G&R community for future work. We suggest that these tools provide the first systematic and efficient framework to quantify uncertainties and optimally identify G&R model parameters.

Turnover of fibrillar collagen in soft biological tissue with application to the expansion of abdominal aortic aneurysms

A better understanding of the inherent properties of vascular tissue to adapt to its mechanical environment is crucial to improve the predictability of biomechanical simulations. Fibrillar collagen in the vascular wall plays a central role in tissue adaptation owing to its relatively short lifetime. Pathological alterations of collagen turnover may fail to result in homeostasis and could be responsible for abdominal aortic aneurysm (AAA) growth at later stages of the disease. For this reason our previously reported multiscale constitutive framework (Martufi, G. & Gasser, T. C. 2011 J. Biomech. 44, 2544-2550 (doi:10.1016/j.jbiomech.2011.07.015)) has been enriched by a collagen turnover model. Specifically, the framework's collagen fibril level allowed a sound integration of vascular wall biology, and the impact of collagen turnover on the macroscopic properties of AAAs was studied. To this end, model parameters were taken from the literature and/or estimated from clinical follow-up data of AAAs (on average 50.7 mm-large). Likewise, the in vivo stretch of the AAA wall was set, such that 10 per cent of collagen fibres were engaged. Results showed that the stretch spectrum, at which collagen fibrils are deposed, is the most influential parameter, i.e. it determines whether the vascular geometry grows, shrinks or remains stable over time. Most importantly, collagen turnover also had a remarkable impact on the macroscopic stress field. It avoided high stress gradients across the vessel wall, thus predicted a physiologically reasonable stress field. Although the constitutive model could be successfully calibrated to match the growth of small AAAs, a rigorous validation against experimental data is crucial to further explore the model's descriptive and predictive capabilities.

SENSITIVITY ANALYSIS OF COMPUTATIONAL STRUCTURAL DYNAMICS IN A CEREBRAL ANEURYSM MODEL TO WALL THICKNESS AND MODEL

Intracranial saccular aneurysms tend to be thin walled and stiffer compared with a normal artery. The current work describes computational structural dynamics (CSD) in an anatomically realistic model of a cerebral aneurysm located in the ophthalmic region, using different wall thickness, model data for the artery and aneurysm, and geometry size. The model was obtained from three-dimensional rotational angiography image data. The wall was assumed three-dimensional hyperelastic solid with different thickness in the artery and in the aneurysm regions. The effects of carotid siphon length are reported. The CSD was solved with the finite elements package ADINA. The predictions of stress and strain on the aneurysm wall were compared.

Constrained Mixture Models as Tools for Testing Competing Hypotheses in Arterial Biomechanics: A Brief Survey

Hypothesis testing via numerical models has emerged as a powerful tool which permits the verification of theoretical frameworks against canonical experimental and clinical observations. Cleverly designed computational experiments also inspire new methodologies by elucidating important biological processes and restricting parametric spaces. Constrained mixture models of arterial growth and remodeling (G&R) can facilitate the design of computational experiments which can bypass technical limitations in the laboratory, by considering illustrative special cases. The resulting data may then inform the design of focused experimental techniques and the development of improved theories. This work is a survey of computational hypothesis-testing studies, which exploit the unique abilities offered by the constrained mixture theory of arterial G&R. Specifically, we explore the core hypotheses integrated in these models, review their basic mathematical conceptualizations, and recapitulate their most salient and illuminating findings. We then assess how a decade's worth of constrained mixture models have contributed to a lucid, emerging picture of G&R mechanisms.

COMPUTATIONAL MODELING OF FORMATION OF A CEREBRAL ANEURYSM UNDER THE INFLUENCE OF SMOOTH MUSCLE CELL RELAXATION

The mechanics of cerebral aneurysm pathogenesis, evolution and rupture are not yet well understood. This paper presents a numerical analysis of the formation of a saccular cerebral aneurysm in for the first time in a 3D model of the basilar artery bifurcation under normal and hypertensive blood pressure. Due to the excessive endothelium derived nitric oxide produced in high wall shear stress, we assumed that smooth muscle cell relaxation is the origin of the aneurysm formation. Arterial wall remodeling under constant tension was considered to be the other mechanism of disease evolution. The wall was constructed from two elastic and hyperelastic isotropic regions. The flow was considered steady, laminar, Newtonian, and incompressible. The fully coupled fluid and structure models were solved with the finite elements package ADINA 8.5. The wall shear stress, effective stress and deformation distributions under normal and hypertensive blood pressure were compared to a healthy bifurcation. The model shows that although the malfunction of the endothelial cell layer and the corresponding smooth muscle cell-related loss of vascular tone is important to the inception of the disease; A saccular aneurysm may not be formed by this mechanism alone, and also requires the fiber-related arterial wall remodeling for further development.

On the biomechanics and mechanobiology of growing skin

Skin displays an impressive functional plasticity, which allows it to adapt gradually to environmental changes. Tissue expansion takes advantage of this adaptation, and induces a controlled in situ skin growth for defect correction in plastic and reconstructive surgery. Stretches beyond the skin's physiological limit invoke several mechanotransduction pathways, which increase mitotic activity and collagen synthesis, ultimately resulting in a net gain in skin surface area. However, the interplay between mechanics and biology during tissue expansion remains unquantified. Here, we present a continuum model for skin growth that summarizes the underlying mechanotransduction pathways collectively in a single phenomenological variable, the strain-driven area growth. We illustrate the governing equations for growing biological membranes, and demonstrate their computational solution within a nonlinear finite element setting. In displacement-controlled equi-biaxial extension tests, the model accurately predicts the experimentally observed histological, mechanical, and structural features of growing skin, both qualitatively and quantitatively. Acute and chronic elastic uniaxial stretches are 25% and 10%, compared to 36% and 10% reported in the literature. Acute and chronic thickness changes are -28% and -12%, compared to -22% and -7% reported in the literature. Chronic fractional weight gain is 3.3, compared to 2.7 for wet weight and 3.3 for dry weight reported in the literature. In two clinical cases of skin expansion in pediatric forehead reconstruction, the model captures the clinically observed mechanical and structural responses, both acutely and chronically. Our results demonstrate that the field theories of continuum mechanics can reliably predict the mechanical manipulation of thin biological membranes by controlling their mechanotransduction pathways through mechanical overstretch. We anticipate that the proposed skin growth model can be generalized to arbitrary biological membranes, and that it can serve as a valuable tool to virtually manipulate living tissues, simply by means of changes in the mechanical environment.

Effects of dispersion of fiber orientation on the mechanical property of the arterial wall

Effects of dispersion of the fiber orientation on the mechanical property of the arterial wall in health and disease subjected to the combined internal pressure and axial loading are examined within the framework of the finite deformation hyper-elasticity theory. Taking into account the residual stress, a two layer thick-walled circular cylindrical tube model with the fiber-reinforced incompressible composite hyper-elastic material is employed. Stress-radius curves and stress distributions of the arterial wall are given in health and disease considering dispersion of the fiber orientation. With instability/bifurcation analysis, it is found that not only the uniform inflation of the arterial wall in health, but also the instability or bifurcation problem for the arterial wall in disease may be described by this model.

Perspectives on biological growth and remodeling

The continuum mechanical treatment of biological growth and remodeling has attracted considerable attention over the past fifteen years. Many aspects of these problems are now well-understood, yet there remain areas in need of significant development from the standpoint of experiments, theory, and computation. In this perspective paper we review the state of the field and highlight open questions, challenges, and avenues for further development.

Numerical investigation of the hemodynamics in anatomically realistic lateral cerebral aneurysms

Hemodynamically induced stress plays an important role in the progression and rupture of cerebral aneurysms. The current work describes computational fluid dynamics (CFD) simulations in anatomically realistic models of cerebral aneurysms. Twenty lateral aneurysms models were investigated. The models were obtained from three-dimensional rotational angiographic imaging data and CFD were studied under the same physiologically representative waveform of inflow. The flow was assumed to be laminar, non-Newtonian, and incompressible. The CFD models were solved with the finite elements package ADINA. Predictions of velocity field and wall shear stress (WSS) on the aneurysms were compared for the different cases. Linear correlations between the WSS on the aneurysm fundus at peak systole for lateral aneurysms with an area index were found.

Modeling of fibroblast-controlled strengthening and remodeling of uniaxially constrained collagen gels

A theoretical model for the remodeling of collagen gels is proposed. The collagen fabric is modeled as a network of collagen fibers, which in turn are composed of collagen fibrils. In the model, the strengthening of collagen fabric is accomplished by fibroblasts, which continuously recruit and attach more collagen fibrils to existing collagen fibers. The fibroblasts also accomplish a reorientation of collagen fibers. Fibroblasts are assumed to reorient collagen fibers toward the direction of maximum material stiffness. The proposed model is applied to experiments in which fibroblasts were inserted into a collagen gel. The model is able to predict the force-strain curves for the experimental collagen gels, and the final distribution of collagen fibers also agrees qualitatively with the experiments.

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.

Constitutive modelling of arteries

This review article is concerned with the mathematical modelling of the mechanical properties of the soft biological tissues that constitute the walls of arteries. Many important aspects of the mechanical behaviour of arterial tissue can be treated on the basis of elasticity theory, and the focus of the article is therefore on the constitutive modelling of the anisotropic and highly nonlinear elastic properties of the artery wall. The discussion focuses primarily on developments over the last decade based on the theory of deformation invariants, in particular invariants that in part capture structural aspects of the tissue, specifically the orientation of collagen fibres, the dispersion in the orientation, and the associated anisotropy of the material properties. The main features of the relevant theory are summarized briefly and particular forms of the elastic strain-energy function are discussed and then applied to an artery considered as a thick-walled circular cylindrical tube in order to illustrate its extension–inflation behaviour. The wide range of applications of the constitutive modelling framework to artery walls in both health and disease and to the other fibrous soft tissues is discussed in detail. Since the main modelling effort in the literature has been on the passive response of arteries, this is also the concern of the major part of this article. A section is nevertheless devoted to reviewing the limited literature within the continuum mechanics framework on the active response of artery walls, i.e. the mechanical behaviour associated with the activation of smooth muscle, a very important but also very challenging topic that requires substantial further development. A final section provides a brief summary of the current state of arterial wall mechanical modelling and points to key areas that need further modelling effort in order to improve understanding of the biomechanics and mechanobiology of arteries and other soft tissues, from the molecular, to the cellular, tissue and organ levels.

Influence of medial collagen organization and axial in situ stretch on saccular cerebral aneurysm growth

A model for saccular cerebral aneurysm growth, proposed by Kroon and Holzapfel (2007, "A Model for Saccular Cerebral Aneurysm Growth in a Human Middle Cerebral Artery," J. Theor. Biol., 247, pp. 775-787; 2008, "Modeling of Saccular Aneurysm Growth in a Human Middle Cerebral Artery," ASME J. Biomech. Eng., 130, p. 051012), is further investigated. A human middle cerebral artery is modeled as a two-layer cylinder where the layers correspond to the media and the adventitia. The immediate loss of media in the location of the aneurysm is taken to be responsible for the initiation of the aneurysm growth. The aneurysm is regarded as a development of the adventitia, which is composed of several distinct layers of collagen fibers perfectly aligned in specified directions. The collagen fibers are the only load-bearing constituent in the aneurysm wall; their production and degradation depend on the stretch of the wall and are responsible for the aneurysm growth. The anisotropy of the surrounding media was modeled using the strain-energy function proposed by Holzapfel et al. (2000, "A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models," J. Elast., 61, pp. 1-48), which is valid for an elastic material with two families of fibers. It was shown that the inclusion of fibers in the media reduced the maximum principal Cauchy stress and the maximum shear stress in the aneurysm wall. The thickness increase in the aneurysm wall due to material growth was also decreased. Varying the fiber angle in the media from a circumferential direction to a deviation of 10 deg from the circumferential direction did, however, only show a little effect. Altering the axial in situ stretch of the artery had a much larger effect in terms of the steady-state shape of the aneurysm and the resulting stresses in the aneurysm wall. The peak values of the maximum principal stress and the thickness increase both became significantly higher for larger axial stretches.

Open Problems in Computational Vascular Biomechanics: Hemodynamics and Arterial Wall Mechanics

The vasculature consists of a complex network of vessels ranging from large arteries to arterioles, capillaries, venules, and veins. This network is vital for the supply of oxygen and nutrients to tissues and the removal of carbon dioxide and waste products from tissues. Because of its primary role as a pressure-driven chemomechanical transport system, it should not be surprising that mechanics plays a vital role in the development and maintenance of the normal vasculature as well as in the progression and treatment of vascular disease. This review highlights some past successes of vascular biomechanics, but emphasizes the need for research that synthesizes complementary advances in molecular biology, biomechanics, medical imaging, computational methods, and computing power for purposes of increasing our understanding of vascular physiology and pathophysiology as well as improving the design of medical devices and clinical interventions, including surgical procedures. That is, computational mechanics has great promise to contribute to the continued improvement of vascular health.

Coupling the Hemodynamic Environment to the Evolution of Cerebral Aneurysms: Computational Framework and Numerical Examples

The physiological mechanisms that give rise to the inception and development of a cerebral aneurysm are accepted to involve the interplay between the local mechanical forces acting on the arterial wall and the biological processes occurring at the cellular level. In fact, the wall shear stresses (WSSs) that act on the endothelial cells are thought to play a pivotal role. A computational framework is proposed to explore the link between the evolution of a cerebral aneurysm and the influence of hemodynamic stimuli that act on the endothelial cells. An aneurysm evolution model, which utilizes a realistic microstructural model of the arterial wall, is combined with detailed 3D hemodynamic solutions. The evolution of the blood flow within the developing aneurysm determines the distributions of the WSS and the spatial WSS gradient (WSSG) that act on the endothelial cell layer of the tissue. Two illustrative examples are considered: Degradation of the elastinous constituents is driven by deviations of WSS or the WSSG from normotensive values. This model provides the basis to further explore the etiology of aneurysmal disease.

Coupling hemodynamics with vascular wall mechanics and mechanobiology to understand intracranial aneurysms

Arteries exhibit a remarkable ability to adapt in response to sustained alterations in hemodynamic loading, to heal in response to injuries, and to compensate in response to diverse disease conditions. Nevertheless, such compensatory adaptations are limited and many vascular disorders, if untreated, lead to significant morbidity or mortality. Parallel advances in vascular biology, medical imaging, biomechanics, and computational methods promise to provide increased insight into many arterial diseases, including intracranial aneurysms. In particular, although it may be possible to identify useful clinical correlations between either the blood flow patterns within or the shape of aneurysms and their rupture-potential, our ultimate goal should be to couple studies of hemodynamics with those of wall mechanics and the underlying mechanobiology so that we can understand better the mechanisms by which aneurysms arise, enlarge, and rupture and thereby identify better methods of treatment. This paper presents one such approach to fluid-solid-growth (FSG) modeling of intracranial aneurysms.

Modelling the mechanical response of elastin for arterial tissue

We compare two constitutive models proposed to model the elastinous constituents of an artery. Holzapfel and Weizsäcker [1998. Biomechanical behavior of the arterial wall and its numerical characterization. Comput. Biol. Med. 28, 377-392] attribute a neo-Hookean response, i.e. Psi=c(I(1)-3)), to the elastin whilst Zulliger et al. [2004a. A strain energy function for arteries accounting for wall composition and structure. J. Biomech. 37, 989-1000] propose Psi=c(I(1)-3)(3/2). We analyse these constitutive models for two specific cases: (i) uniaxial extension of an elastinous sheet; (ii) inflation of a cylindrical elastinous membrane. For case (i) we illustrate the functional relationships between: (a) the Cauchy stress (CS) and the Green-Lagrange (GL) strain; (b) the tangent modulus (gradient of the CS-GL strain curve) and linearised strain. The predicted mechanical responses are compared with recent uniaxial extension tests on elastin [Gundiah, N., Ratcliffe, M.B., Pruitt, L.A., 2007. Determination of strain energy function for arterial elastin: experiments using histology and mechanical tests. J. Biomech. 40, 586-594; Lillie, M.A., Gosline, J.M., 2007a. Limits to the durability of arterial elastic tissue. Biomaterials 28, 2021-2031; 2007b. Mechanical properties of elastin along the thoracic aorta in the pig. J. Biomech. 40, 2214-2221]. The neo-Hookean model accurately predicts the mechanical response of a single elastin fibre. However, it is unable to accurately capture the mechanical response of arterial elastin, e.g. the initial toe region of arterial elastin (if it exists) or the gradual increase in modulus of arterial elastin that occurs as it is stretched. The alternative constitutive model (n=32) yields a nonlinear mechanical response that departs from recent uniaxial test data mentioned above, for the same stretch range. For case (ii) we illustrate the pressure-circumferential stretch relationships and the gradients of the pressure-circumferential stretch curves: significant qualitative differences are observed. For the neo-Hookean model, the gradient decreases rapidly to zero, however, for n=32, the gradient decreases more gradually to a constant value. We conclude that whilst the neo-Hookean model has limitations, it appears to capture more accurately the mechanical response of elastin.

VASCULAR MECHANICS, MECHANOBIOLOGY, AND REMODELING

Arteries exhibit a remarkable ability to adapt in response to sustained alterations in hemodynamic loading as well as in response to disease, injury, and clinical treatment. A better understanding of such adaptations will be aided greatly by formulating, testing, and refining appropriate theoretical frameworks for modeling the biomechanics and associated mechanobiology. The goal of this brief review is to highlight some recent developments in the use of a constrained mixture theory of arterial growth and remodeling, with particular attention to the requisite constitutive relations, and to highlight future directions of needed research.

A theoretical model for fibroblast-controlled growth of saccular cerebral aneurysms

A new theoretical model for the growth of saccular cerebral aneurysms is proposed by extending the recent constitutive framework of Kroon and Holzapfel [2007a. A model for saccular cerebral aneurysm growth by collagen fibre remodelling. J. Theor. Biol. 247, 775-787]. The continuous turnover of collagen is taken to be the driving mechanism in aneurysmal growth. The collagen production rate depends on the magnitude of the cyclic deformation of fibroblasts, caused by the pulsating blood pressure during the cardiac cycle. The volume density of fibroblasts in the aneurysmal tissue is taken to be constant throughout the growth process. The growth model is assessed by considering the inflation of an axisymmetric membranous piece of aneurysmal tissue, with material characteristics representative of a cerebral aneurysm. The diastolic and systolic states of the aneurysm are computed, together with its load-free state. It turns out that the value of collagen pre-stretch, that determines growth speed and stability of the aneurysm, is of pivotal importance. The model is able to predict aneurysms with typical berry-like shapes observed clinically, and the predicted wall stresses correlate well with the experimentally obtained ultimate stresses of this type of tissue. The model predicts that aneurysms should fail when reaching a size of about 1.2-3.6mm, which is smaller than what has been clinically observed. With some refinements, the model may, however, be used to predict future growth of diagnosed aneurysms.

Modeling of Saccular Aneurysm Growth in a Human Middle Cerebral Artery

Saccular aneurysm growth in a human middle cerebral artery is modeled. The aneurysm growth model was presented in a companion paper by Kroon and Holzapfel (“A Model for Saccular Cerebral Aneurysm Growth by Collagen Fibre Remodelling,” J. Theor. Biol., in press) and was assessed there for axisymmetric growth. The aneurysm growth model is now evaluated for a more realistic setting. The middle cerebral artery is modeled as a two-layered cylinder, where the layers correspond to the media and the adventitia. An instant loss of the media in a region of the artery wall initiates the growth of the saccular aneurysm. The aneurysm wall is assumed to be a development of the adventitia of the original healthy artery, and collagen is assumed to be the only load-bearing constituent in the adventitia and in the aneurysm wall. The collagen is organized in a number of distinct layers where fibers in a specific layer are perfectly aligned in a certain fiber direction. The production of new collagen is taken to depend on the stretching of the aneurysm wall, and the continuous remodeling of the collagen fibers is responsible for the aneurysm growth. The general behavior of the growth model is investigated and also the influence of the structural organization of the collagen fabric. The analysis underlines the fact that the material behavior of aneurysmal tissue cannot be expected to be isotropic. The model predictions agree well with clinical and experimental results, for example, in terms of aneurysm size and shape, wall stress levels, and wall thickness.