A microstructurally motivated model of arterial wall mechanics with mechanobiological implications

Finite Element Simulation of Opening Angle Response of Porcine Aortas Using Layer Specific GAG Distributions in One and Two Layered Solid Matrices

PURPOSE

Recent studies have identified an effect of glycosaminoglycans (GAG) on residual stresses in the aorta, underscoring the need to better understand their biomechanical roles.

METHODS

Aortic ring models for each of the ascending, arch and descending thoracic regions of the porcine thoracic aorta were created in FEBioStudio, using a framework that incorporates the Donnan osmotic swelling in a porous solid matrix. The distribution of fixed charge densities (FCD) through the thickness of the tissue was prescribed as calculated from experimentally quantified sulfated GAG mural distributions. Material parameters for the solid matrix, modeled using a Holmes-Mow constitutive law, were optimized using data from biaxial tensile tests. In addition to modelling the solid matrix as one layer, two layers were considered to capture the differences between the intima-media and the adventitia, for which various stiffness ratios were explored.

RESULTS

As the stiffness of the adventitia with respect to that of the media increased, the simulated opening angle increased. The opening angle also decreased from the ascending to the descending thoracic region in both one- and two-layered solid matrices models. The simulated results were compared against the experimental contribution of GAG to the opening angle, as previously quantified via enzymatic GAG-depletion. When using one layer for the solid matrix, the errors between the simulated opening angles and the experimental contribution of GAG to the opening angle were respectively 28%, 15% and 23% in the ascending, arch and descending thoracic regions. When using two layers for the solid matrix, the smallest errors in the ascending and arch regions were 21% and 5% when the intima-media was modelled as 10 times stiffer, and as twice stiffer than the adventitia, respectively, and 23% in the descending thoracic regions when the intima-media and adventitia shared similar mechanical properties.

CONCLUSIONS

Overall, this study demonstrates that GAG partially contribute to circumferential residual stress, and that GAG swelling is one of several regulators of the opening angle. The minor discrepancies between simulated and experimental opening angles imply that the contribution of GAG extends beyond mere swelling, aligning with previous experimental indications of their interaction with ECM fibers in determining the opening angle.

Mapping the unicellular transcriptome of the ascending thoracic aorta to changes in mechanosensing and mechanoadaptation during aging

Aortic stiffening is an inevitable manifestation of chronological aging, yet the mechano-molecular programs that orchestrate region- and layer-specific adaptations along the length and through the wall of the aorta are incompletely defined. Here, we show that the decline in passive cyclic distensibility is more pronounced in the ascending thoracic aorta (ATA) compared to distal segments of the aorta and that collagen content increases in both the medial and adventitial compartments of the ATA during aging. The single-cell RNA sequencing of aged ATA tissues reveals altered cellular senescence, remodeling, and inflammatory responses accompanied by enrichment of T-lymphocytes and rarefaction of vascular smooth muscle cells, compared to young samples. T lymphocyte clusters accumulate in the adventitia, while the activation of mechanosensitive Piezo-1 enhances vasoconstriction and contributes to the overall functional decline of ATA tissues. These results portray the immuno-mechanical aging of the ATA as a process that culminates in a stiffer conduit permissive to the accrual of multi-gerogenic signals priming to disease development.

AT1b receptors contribute to regional disparities in angiotensin II mediated aortic remodelling in mice

The renin-angiotensin system plays a key role in regulating blood pressure, which has motivated many investigations of associated mouse models of hypertensive arterial remodelling. Such studies typically focus on histological and cell biological changes, not wall mechanics. This study explores tissue-level ramifications of chronic angiotensin II infusion in wild-type (WT) and type 1b angiotensin II (AngII) receptor null (Agtr1b -/-) mice. Biaxial biomechanical and immunohistological changes were quantified and compared in the thoracic and abdominal aorta in these mice following 14 and 28 days of angiotensin II infusion. Preliminary results showed that changes were largely independent of sex. Associated thickening and stiffening of the aortic wall in male mice differed significantly between thoracic and abdominal regions and between genotypes. Notwithstanding multiple biomechanical changes in both WT and Agtr1b -/- mice, AngII infusion caused distinctive wall thickening and inflammation in the descending thoracic aorta of WT, but not Agtr1b -/-, mice. Our study underscores the importance of exploring differential roles of receptor-dependent angiotensin II signalling along the aorta and its influence on distinct cell types involved in regional histomechanical remodelling. Disrupting the AT1b receptor primarily affected inflammatory cell responses and smooth muscle contractility, suggesting potential therapeutic targets.

Vascular Remodeling During Late-Gestation Pregnancy: An In-Vitro Assessment of the Murine Ascending Thoracic Aorta

Maternal mortality due to cardiovascular disease is a rising concern in the U.S. Pregnancy triggers changes in the circulatory system, potentially influencing the structure of the central vasculature. Evidence suggests a link between a woman's pregnancy history and future cardiovascular health, but our understanding remains limited. To fill this gap, we examined the passive mechanics of the murine ascending thoracic aorta during late gestation. By performing biaxial mechanical testing on the ascending aorta, we were able to characterize the mechanical properties of both control and late-gestation tissues. By examining mechanical, structural, and geometric properties, we confirmed that remodeling of the aortic wall occurred. Morphological and mechanical properties of the tissue indicated an outward expansion of the tissue, as reflected in changes in wall thickness (∼12% increase) and luminal diameter (∼6% increase) at its physiologically loaded state in the pregnant group. With these geometric adaptations and despite increased hemodynamic loads, pregnancy did not induce significant changes in the tensile wall stress at the similar physiological pressure levels of the pregnant and control tissues. The alterations also included reduced intrinsic stiffness in the circumferential direction (∼18%) and reduced structural stiffness (∼26%) in the pregnant group. The observed vascular remodeling maintained the elastic stored energy of the aortic wall under systolic loads, indicating preservation of vascular function. Data from our study of pregnancy-related vascular remodeling will provide valuable insights for future investigations of maternal cardiovascular health.

Cell signaling and tissue remodeling in the pulmonary autograft after the Ross procedure: A computational study

In the Ross procedure, a patient's pulmonary valve is transplanted in the aortic position. Despite advantages of this surgery, reoperation is still needed in many cases due to excessive dilatation of the pulmonary autograft. To further understand the failure mechanisms, we propose a multiscale model predicting adaptive processes in the autograft at the cell and tissue scale. The cell-scale model consists of a network model, that includes important signaling pathways and relations between relevant transcription factors and their target genes. The resulting gene activity leads to changes in the mechanical properties of the tissue, modeled as a constrained mixture of collagen, elastin and smooth muscle. The multiscale model is calibrated with findings from experiments in which seven sheep underwent the Ross procedure. The model is then validated against a different set of sheep experiments, for which a qualitative agreement between model and experiment is found. Model outcomes at the cell scale, including the activity of genes and transcription factors, also match experimentally obtained transcriptomics data.

Effect of rapid cooling, frozen storage, and thawing on the passive viscoelastic properties and structure of the rat aorta

Biological tissues decay over time after harvesting, which alters their biomechanical properties. This poses logistical challenges for studies investigating passive arterial biomechanics as tissues need to be characterized shortly after excision. Freezing and cryopreservation methods can help alleviate the need for biomechanical testing of fresh tissue in human ex vivo studies. However, these methods tend to eliminate or reduce arterial cell functionality and affect passive biomechanics. Furthermore, their impact on dynamic arterial biomechanics remains unknown despite arterial viscoelastic properties being an integral component contributing to arterial stiffness under in vivo loading conditions. The present study aims to investigate the impact of rapid cooling and subsequent storage at -80 °C on the passive viscoelastic properties of arterial tissue and aid in ascertaining whether this is a suitable method to delay tissue analysis for studies investigating passive arterial biomechanics. Control and frozen abdominal rat aorta segments were quasi-statically and dynamically tested using a biaxial testing set-up. The results were modeled using a constituent-based quasi-linear viscoelastic modeling framework, yielding directional stiffness parameters, individual constituent biomechanical contributions, and a quantification of viscoelastic stiffening under dynamic pressurization conditions. Frozen samples displayed significantly decreased wall thickness, viscoelastic dissipation, viscoelastic stiffening, and significantly decreased circumferential deformation with changes in luminal pressure. Furthermore, frozen samples displayed significantly increased circumferential stiffness, pulse wave velocity, and collagen load bearing. Consequently, these changes should be considered when utilizing this tissue preservation method to delay biomechanical characterization of rat aortic tissue.

Growth and remodeling of the dissected membrane in an idealized dissected aorta model

While transitioning from the acute to chronic phase, the wall of a dissected aorta often expands in diameter and adaptations in thickness and microstructure take place in the dissected membrane. Including the mechanisms, leading to these changes, in a computational model is expected to improve the accuracy of predictions of the long-term complications and optimal treatment timing of dissection patients. An idealized dissected wall was modeled to represent the elastin and collagen production and/or degradation imposed by stress- and inflammation-mediated growth and remodeling, using the homogenized constrained mixture theory. As no optimal growth and remodeling parameters have been defined for aortic dissections, a Latin hypercube sampling with 1000 parameter combinations was assessed for four inflammation patterns, with a varying spatial extent (full/local) and temporal evolution (permanent/transient). The dissected membrane thickening and microstructure was considered together with the diameter expansion over a period of 90 days. The highest success rate was found for the transient inflammation patterns, with about 15% of the samples leading to converged solutions after 90 days. Clinically observed thickening rates were found for 2-4% of the transient inflammation samples, which represented median total diameter expansion rates of about 5 mm/year. The dissected membrane microstructure showed an elastin decrease and, in most cases, a collagen increase. In conclusion, the model with the transient inflammation pattern allowed the reproduction of clinically observed dissected membrane thickening rates, diameter expansion rates and adaptations in microstructure, thus providing guidance in reducing the parameter space in growth and remodeling models of aortic dissections.

An inverse fitting strategy to determine the constrained mixture model parameters: application in patient-specific aorta

The Constrained Mixture Model (CMM) is a novel approach to describe arterial wall mechanics, whose formulation is based on a referential physiological state. The CMM considers the arterial wall as a mixture of load-bearing constituents, each of them with characteristic mass fraction, material properties, and deposition stretch levels from its stress-free state to the in-vivo configuration. Although some reports of this model successfully assess its capabilities, they barely explore experimental approaches to model patient-specific scenarios. In this sense, we propose an iterative fitting procedure of numerical-experimental nature to determine material parameters and deposition stretch values. To this end, the model has been implemented in a finite element framework, and it is calibrated using reported experimental data of descending thoracic aorta. The main results obtained from the proposed procedure consist of a set of material parameters for each constituent. Moreover, a relationship between deposition stretches and residual strain measurements (opening angle and axial stretch) has been numerically proved, establishing a strong consistency between the model and experimental data.

Constituent-based quasi-linear viscoelasticity: a revised quasi-linear modelling framework to capture nonlinear viscoelasticity in arteries

Arteries exhibit fully nonlinear viscoelastic behaviours (i.e. both elastically and viscously nonlinear). While elastically nonlinear arterial models are well established, effective mathematical descriptions of nonlinear viscoelasticity are lacking. Quasi-linear viscoelasticity (QLV) offers a convenient way to mathematically describe viscoelasticity, but its viscous linearity assumption is unsuitable for whole-wall vascular applications. Conversely, application of fully nonlinear viscoelastic models, involving deformation-dependent viscous parameters, to experimental data is impractical and often reduces to identifying specific solutions for each tested loading condition. The present study aims to address this limitation: By applying QLV theory at the wall constituent rather than at the whole-wall level, the deformation-dependent relative contribution of the constituents allows to capture nonlinear viscoelasticity with a unique set of deformation-independent model parameters. Five murine common carotid arteries were subjected to a protocol of quasi-static and harmonic, pseudo-physiological biaxial loading conditions to characterise their viscoelastic behaviour. The arterial wall was modelled as a constrained mixture of an isotropic elastin matrix and four families of collagen fibres. Constituent-based QLV was implemented by assigning different relaxation functions to collagen- and elastin-borne parts of the wall stress. Nonlinearity in viscoelasticity was assessed via the pressure dependency of the dynamic-to-quasi-static stiffness ratio. The experimentally measured ratio increased with pressure, from 1.03 [Formula: see text] 0.03 (mean [Formula: see text] standard deviation) at 80-40 mmHg to 1.58 [Formula: see text] 0.22 at 160-120 mmHg. Constituent-based QLV captured well this trend by attributing the wall viscosity predominantly to collagen fibres, whose recruitment starts at physiological pressures. In conclusion, constituent-based QLV offers a practical and effective solution to model arterial viscoelasticity.

Computational modeling reveals inflammation-driven dilatation of the pulmonary autograft in aortic position

The pulmonary autograft in the Ross procedure, where the aortic valve is replaced by the patient's own pulmonary valve, is prone to failure due to dilatation. This is likely caused by tissue degradation and maladaptation, triggered by the higher experienced mechanical loads in aortic position. In order to further grasp the causes of dilatation, this study presents a model for tissue growth and remodeling of the pulmonary autograft, using the homogenized constrained mixture theory and equations for immuno- and mechano-mediated mass turnover. The model outcomes, compared to experimental data from an animal model of the pulmonary autograft in aortic position, show that inflammation likely plays an important role in the mass turnover of the tissue constituents and therefore in the autograft dilatation over time. We show a better match and prediction of long-term outcomes assuming immuno-mediated mass turnover, and show that there is no linear correlation between the stress-state of the material and mass production. Therefore, not only mechanobiological homeostatic adaption should be taken into account in the development of growth and remodeling models for arterial tissue in similar applications, but also inflammatory processes.

Simple Models of Complex Mechanics for Improved Hypertension Care: Learning to De-stiffen Arteries

Arteries can stiffen via different mechanisms due to the distending effects of blood pressure, the extracellular (ECM) and vascular smooth muscle cells (VSMC). This short review discusses how these simple models can be applied to the complex biomechanics of arteries to gain physiological insight into why an individual's arteries are stiff and identify new therapeutic strategies. In the Multi-Ethnic Study of Atherosclerosis, the important question of whether arteries stiffen with aging due to load-dependent or structural stiffening was investigated. Structural stiffening was consistently observed with aging, but load-dependent stiffening was highly variable. Importantly, the high load-dependent stiffness was associated with future cardiovascular disease events, but structural stiffness was not. Clinical studies in older, hypertensive adults surprisingly show that decreasing vascular smooth muscle tone can cause clinically significant increases in arterial stiffness. To understand this paradox, the author developed a model simple enough for clinical data but with biologically relevant extracellular matrix (ECM) and vascular smooth muscle cell (VSMC) stiffness parameters. The effect of VSMC tone on arterial stiffness depends on the ECM-VSMC stiffness ratio. Future research is needed to develop a framework that incorporates both the blood pressure dependence of arterial stiffness and the VSMC-ECM interaction on hemodynamics. This could result in personalized arterial stiffness treatments and improved CVD outcomes. The subtitle of this review is "Learning to De-Stiffen Arteries" because our results have so far only shown that we can acutely make arteries stiffer. We are optimistic though that the findings and the analytic techniques covered here will be one of the many steps along the path of the arterial stiffness research community learning how to de-stiffen arteries.

How to implement constrained mixture growth and remodeling algorithms for soft biological tissues

Biological soft tissues are constantly adapting to their mechanical environment and remodel to restore certain mechanobiological homeostatic conditions. These effects can be modeled using the constrained mixture theory, that assumes degradation of material over time and the gradual replacement of extant material by newly deposited material. While this theory presents an elegant way to grasp phenomena of growth and remodeling in soft biological tissues, implementation difficulties may arise. Therefore, we give a detailed overview of the mathematical description of the constrained mixture theory and its homogenized equivalent, and present practical suggestions to numerically implement the theories. These implementations are thoroughly tested with multiple example growth and remodeling models. Results show a good correspondence between both theories, with the homogenized theory favored in terms of time efficiency. Results of a step time convergence study show the importance of choosing a small enough time step, especially when using the classical theory.

A novel constitutive model considering the role of elastic lamellae' structural heterogeneity in homogenizing transmural stress distribution in arteries

Understanding how the homeostatic stress state can be reached in arterial tissues can provide new insights into vascular physiology. Even though the function of maintaining homeostasis is often linked to the concentric layers of medial elastic lamellae, how the lamellae are capable of evenly distributing the stress transmurally remains to be understood. The recent microstructural study by Yu et al. (2018 J. R. Soc. Interface 15, 20180492) revealed that, circumferentially, lamellar layers closer to the lumen are wavier than the ones further away from it and, thus, experience more unfolding when subjected to blood pressure. Motivated by this peculiar finding, the current study, for the first time, proposes a novel approach to model elastic lamellae and such structural heterogeneity using the extensible worm-like chain model. When implemented into the material description of the conventional two-layer artery model, in which adventitial collagen is modelled using the inextensible worm-like chain model, it is demonstrated that structural heterogeneity in elastic lamellae plays an important role in dictating transmural stress distribution and, therefore, the homeostasis of the arterial wall.

Mechanical characterization of human umbilical arteries by thick-walled models: Enhanced vascular compliance by removing an abluminal lining

The decellularized human umbilical artery (HUA) is considered as a promising option for small-diameter, tissue-engineered vascular grafts (TEVGs). Our previous study showed that the HUA bears a thin, watertight lining on its outermost abluminal surface. Removal of this abluminal lining layer improves efficacy of the perfusion-assisted decellularization of the HUA and increases its compliance. As stress across the wall is believed to affect growth and remodeling of the TEVG, it is imperative to mechanically characterize the HUA using thick-walled models. Combining inflation experiments and computational methods, we investigate the mechanical properties of the HUA before and after the abluminal lining removal to examine the HUA's wall mechanics. The inflation tests of five HUAs were performed to obtain the mechanical and geometrical response of the vessel wall before and after the lining layer removal. Using nonlinear hyperelastic models, the same responses are obtained computationally using the thick-walled models. The experimental data are incorporated into the computational models to estimate the mechanical and orientation parameters of the fibers and isotropic matrix of different layers in the HUAs. The parameter fitting of both thick-walled models (before and after the abluminal lining removal) results in most of the R-squared values for measuring the goodness of fitting to be over 0.90 for all samples. The compliance of the HUA increases from a mean value of 2.60% per 100 mmHg before the removal of the lining to a mean value of 4.21% per 100 mmHg after the removal. The results reveal that, although the abluminal lining is thin, it is stiff and capable of supporting majority of the high luminal pressure, and that the inner layer is far less stressed than the abluminal lining. Computational simulations also show that removal of the abluminal lining increases the circumferential wall stress by up to 280 kPa under the in vivo luminal pressure. The integrated computational and experimental approaches provide more accurate estimates of the material behaviors of HUAs employed in grafts and, in turn, the study enhances our understanding of interactions between the graft and the native vessel on vascular growth and remodeling.

G2Φnet: Relating genotype and biomechanical phenotype of tissues with deep learning

Many genetic mutations adversely affect the structure and function of load-bearing soft tissues, with clinical sequelae often responsible for disability or death. Parallel advances in genetics and histomechanical characterization provide significant insight into these conditions, but there remains a pressing need to integrate such information. We present a novel genotype-to-biomechanical phenotype neural network (G2Φnet) for characterizing and classifying biomechanical properties of soft tissues, which serve as important functional readouts of tissue health or disease. We illustrate the utility of our approach by inferring the nonlinear, genotype-dependent constitutive behavior of the aorta for four mouse models involving defects or deficiencies in extracellular constituents. We show that G2Φnet can infer the biomechanical response while simultaneously ascribing the associated genotype by utilizing limited, noisy, and unstructured experimental data. More broadly, G2Φnet provides a powerful method and a paradigm shift for correlating genotype and biomechanical phenotype quantitatively, promising a better understanding of their interplay in biological tissues.

Shear stress and intravascular pressure effects on vascular dynamics: two-phase blood flow in elastic microvessels accounting for the passive stresses

We study the steady hemodynamics in physiological elastic microvessels proposing an advanced fluid-structure interaction model. The arteriolar tissue is modeled as a two-layer fiber-reinforced hyperelastic material representing its Media and Adventitia layers. The constitutive model employed (Holzapfel et al. in J Elast 61:1-48, 2000) is parametrized via available data on stress-strain experiments for arterioles. The model is completed by simulating the blood/plasma flow in the lumen, using the thixotropic elasto-viscoplastic model in its core, and the linear Phan-Thien and Tanner viscoelastic model in its annular part. The Cell-Free Layer (CFL) and the Fåhraeus and Fåhraeus-Lindqvist effects are considered via analytical expressions based on experimental data (Giannokostas et al. in Materials (Basel) 14:367, 2021b). The coupling between tissue deformation and blood flow is achieved through the experimentally verified pressure-shear hypothesis (Pries et al. Circ Res 77:1017-1023, 1995). Our calculations confirm that the increase in the reference inner radius produces larger expansion. Also, by increasing the intraluminal pressure, the thinning of the walls is more pronounced and it may reach 40% of the initial thickness. Comparing our predictions with those in rigid-wall microtubes, we conclude that apart from the vital importance of vasodilation, there is an up to 25% reduction in wall shear stress. The passive vasodilation contributes to the decrease in the tissue stress fields and affects the hemodynamic features such as the CFL thickness, reducing the plasma layer when blood flows in vessels with elastic walls, in quantitative agreement with previous experiments. Our calculations verify the correctness of the pressure-shear hypothesis but not that of the Laplace law.

Smooth muscle tone alters arterial stiffness: the importance of the extracellular matrix to vascular smooth muscle stiffness ratio

BACKGROUND

Recent studies show that vascular smooth muscle (VSM) is more important to elastic artery mechanics than previously believed. It remains unclear whether increased VSM tone increases or decreases arterial stiffness.

METHODS AND RESULTS

We developed a novel arterial mechanics model based on pressure-diameter relationships that incorporates the contributions of extracellular matrix (ECM) and VSM to arterial stiffness measures. This model is advantageous because it simple enough to use with limited clinical data but has biologically relevant parameters which include ECM stiffness, VSM stiffness, and VSM tone. The model was used to retrospectively analyze the effects of nitroglycerin-induced vasodilation in four clinical studies. Stiffness parameters were modeled for five arterial regions including both elastic and muscular arteries. The model describes complex experimental data with changing VSM tone and blood pressure. Our analysis found that when ECM is less stiff than VSM, increasing VSM tone increases arterial stiffness. The opposite is seen when ECM is stiffer than VSM, increasing VSM tone decreases stiffness. Our results also suggest that VSM tone is a compensatory mechanism for elevated ECM stiffness in hypertensive individuals.

CONCLUSION

Based on retrospective analysis of four clinical studies, we propose a simple hypothesis for the role of VSM tone on arterial stiffness: increased VSM tone increases arterial stiffness when VSM is stiffer than ECM and decreases arterial stiffness when ECM is stiffer than VSM. This hypothesis and the methods used in this study could have important implications for understanding arterial physiology in both hypertension and cardiovascular disease and deserve further exploration.

About prestretch in homogenized constrained mixture models simulating growth and remodeling in patient-specific aortic geometries

Evolution of mechanical and structural properties in the Ascending Thoracic Aorta (ATA) is the results of complex mechanobiological processes. In this work, we address some numerical challenges in order to elaborate computational models of these processes. For that, we extend the state of the art of homogenized constrained mixture (hCM) models. In these models, prestretches are assigned to the mixed constituents in order to ensure local mechanical equilibrium macroscopically, and to maintain a homeostatic level of tension in collagen fibers microscopically. Although the initial prestretches were assumed as homogeneous in idealized straight tubes, more elaborate prestretch distributions need to be considered for curved geometrical models such as patient-specific ATA. Therefore, we introduce prestretches having a three-dimensional gradient across the ATA geometry in the homeostatic reference state. We test different schemes with the objective to ensure stable growth and remodeling (G&R) simulations on patient-specific curved vessels. In these simulations, aneurysm progression is triggered by tissue changes in the constituents such as mass degradation of intramural elastin. The results show that the initial prestretches are not only critical for the stability of numerical simulations, but they also affect the G&R response. Eventually, we submit that initial conditions required for G&R simulations need to be identified regionally for ensuring realistic patient-specific predictions of aneurysm progression.

Intramural Distributions of GAGs and Collagen vs. Opening Angle of the Intact Porcine Aortic Wall

The heterogeneity and contribution of collagen and elastin to residual stresses have been thoroughly studied, but more recently, glycosaminoglycans (GAGs) also emerged as potential regulators. In this study, the opening angle of aortic rings (an indicator of circumferential residual stresses) and the mural distributions of sulfated GAGs (sGAG), collagen, and elastin were quantified in the ascending, aortic arch and descending thoracic regions of 5- to 6-month-old pigs. The opening angle correlated positively with the aortic ring's mean radius and thickness, with good and moderate correlations respectively. The correlations between the sGAG, collagen, elastin, and collagen:sGAG ratio and the opening angle were evaluated to identify aortic compositional factors that could play roles in regulating circumferential residual stresses. The total collagen:sGAG ratio displayed the strongest correlation with the opening angle (r = - 0.715, p < 0.001), followed by the total sGAG content which demonstrated a good correlation (r = 0.623, p < 0.001). Additionally, the intramural gradients of collagen, sGAG and collagen:sGAG correlated moderately with the opening angle. We propose that, in addition to the individual role sGAG play through their content and intramural gradient, the interaction between collagen and sGAG should be considered when evaluating circumferential residual stresses in the aorta.

Growth and remodeling in the pulmonary autograft: Computational evaluation using kinematic growth models and constrained mixture theory

Computational investigations of how soft tissues grow and remodel are gaining more and more interest and several growth and remodeling theories have been developed. Roughly, two main groups of theories for soft tissues can be distinguished: kinematic-based growth theory and theories based on constrained mixture theory. Our goal was to apply these two theories on the same experimental data. Within the experiment, a pulmonary artery was exposed to systemic conditions. The change in diameter was followed-up over time. A mechanical and microstructural analysis of native pulmonary artery and pulmonary autograft was conducted. Whereas the kinematic-based growth theory is able to accurately capture the growth of the tissue, it does not account for the mechanobiological processes causing this growth. The constrained mixture theory takes into account the mechanobiological processes including removal, deposition and adaptation of all structural constituents, allowing us to simulate a changing microstructure and mechanical behavior.

Progressive Microstructural Deterioration Dictates Evolving Biomechanical Dysfunction in the Marfan Aorta

Medial deterioration leading to thoracic aortic aneurysms arises from multiple causes, chief among them mutations to the gene that encodes fibrillin-1 and leads to Marfan syndrome. Fibrillin-1 microfibrils associate with elastin to form elastic fibers, which are essential structural, functional, and instructional components of the normal aortic wall. Compromised elastic fibers adversely impact overall structural integrity and alter smooth muscle cell phenotype. Despite significant progress in characterizing clinical, histopathological, and mechanical aspects of fibrillin-1 related aortopathies, a direct correlation between the progression of microstructural defects and the associated mechanical properties that dictate aortic functionality remains wanting. In this paper, age-matched wild-type, Fbn1 C1041G/+, and Fbn1 mgR/mgR mouse models were selected to represent three stages of increasing severity of the Marfan aortic phenotype. Ex vivo multiphoton imaging and biaxial mechanical testing of the ascending and descending thoracic aorta under physiological loading conditions demonstrated that elastic fiber defects, collagen fiber remodeling, and cell reorganization increase with increasing dilatation. Three-dimensional microstructural characterization further revealed radial patterns of medial degeneration that become more uniform with increasing dilatation while correlating strongly with increased circumferential material stiffness and decreased elastic energy storage, both of which comprise aortic functionality.

Uncertainty quantification in subject-specific estimation of local vessel mechanical properties

Quantitative estimation of local mechanical properties remains critically important in the ongoing effort to elucidate how blood vessels establish, maintain, or lose mechanical homeostasis. Recent advances based on panoramic digital image correlation (pDIC) have made high-fidelity 3D reconstructions of small-animal (e.g., murine) vessels possible when imaged in a variety of quasi-statically loaded configurations. While we have previously developed and validated inverse modeling approaches to translate pDIC-measured surface deformations into biomechanical metrics of interest, our workflow did not heretofore include a methodology to quantify uncertainties associated with local point estimates of mechanical properties. This limitation has compromised our ability to infer biomechanical properties on a subject-specific basis, such as whether stiffness differs significantly between multiple material locations on the same vessel or whether stiffness differs significantly between multiple vessels at a corresponding material location. In the present study, we have integrated a novel uncertainty quantification and propagation pipeline within our inverse modeling approach, relying on empirical and analytic Bayesian techniques. To demonstrate the approach, we present illustrative results for the ascending thoracic aorta from three mouse models, quantifying uncertainties in constitutive model parameters as well as circumferential and axial tangent stiffness. Our extended workflow not only allows parameter uncertainties to be systematically reported, but also facilitates both subject-specific and group-level statistical analyses of the mechanics of the vessel wall.

Investigation of Fiber-Driven Mechanical Behavior of Human and Porcine Bladder Tissue Tested Under Identical Conditions

The urinary bladder is a highly dynamic organ that undergoes large deformations several times per day. Mechanical characteristics of the tissue are crucial in determining the function and dysfunction of the organ. Yet, literature reporting on the mechanical properties of human bladder tissue is scarce and, at times, contradictory. In this study, we focused on mechanically testing tissue from both human and pig bladders using identical protocols to validate the use of pigs as a model for the human bladder. Furthermore, we tested the effect of two treatments on tissue mechanical properties. Namely, elastase to digest elastin fibers, and oxybutynin to reduce smooth muscle cell spasticity. Additionally, mechanical properties based on the anatomical direction of testing were evaluated. We implemented two different material models to aid in the interpretation of the experimental results. We found that human tissue behaves similarly to pig tissue at high deformations (collagen-dominated behavior) while we detected differences between the species at low deformations (amorphous matrix-dominated behavior). Our results also suggest that elastin could play a role in determining the behavior of the fiber network. Finally, we confirmed the anisotropy of the tissue, which reached higher stresses in the transverse direction when compared to the longitudinal direction.

From Uniaxial Testing of Isolated Layers to a Tri-Layered Arterial Wall: A Novel Constitutive Modelling Framework

Mechanical testing and constitutive modelling of isolated arterial layers yields insight into the individual layers' mechanical properties, but per se fails to recapitulate the in vivo loading state, neglecting layer-specific residual stresses. The aim of this study was to develop a testing/modelling framework that integrates layer-specific uniaxial testing data into a three-layered model of the arterial wall, thereby enabling study of layer-specific mechanics under realistic (patho)physiological conditions. Circumferentially and axially oriented strips of pig thoracic aortas (n = 10) were tested uniaxially. Individual arterial layers were then isolated from the wall, tested, and their mechanical behaviour modelled using a hyperelastic strain energy function. Subsequently, the three layers were computationally assembled into a single flat-walled sample, deformed into a cylindrical vessel, and subjected to physiological tension-inflation. At the in vivo axial stretch of 1.10 ± 0.03, average circumferential wall stress was 75 ± 9 kPa at 100 mmHg, which almost doubled to 138 ± 15 kPa at 160 mmHg. A ~ 200% stiffening of the adventitia over the 60 mmHg pressure increase shifted layer-specific load-bearing from the media (65 ± 10% → 61 ± 14%) to the adventitia (28 ± 9% → 32 ± 14%). Our approach provides valuable insight into the (patho)physiological mechanical roles of individual arterial layers at different loading states, and can be implemented conveniently using simple, inexpensive and widely available uniaxial testing equipment.

Constrained Mixture Models of Soft Tissue Growth and Remodeling - Twenty Years After

Soft biological tissues compromise diverse cell types and extracellular matrix constituents, each of which can possess individual natural configurations, material properties, and rates of turnover. For this reason, mixture-based models of growth (changes in mass) and remodeling (change in microstructure) are well-suited for studying tissue adaptations, disease progression, and responses to injury or clinical intervention. Such approaches also can be used to design improved tissue engineered constructs to repair, replace, or regenerate tissues. Focusing on blood vessels as archetypes of soft tissues, this paper reviews a constrained mixture theory introduced twenty years ago and explores its usage since by contrasting simulations of diverse vascular conditions. The discussion is framed within the concept of mechanical homeostasis, with consideration of solid-fluid interactions, inflammation, and cell signaling highlighting both past accomplishments and future opportunities as we seek to understand better the evolving composition, geometry, and material behaviors of soft tissues under complex conditions.

A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue

The heterogeneous, nonlinear, anisotropic material behavior of biological tissues makes precise definition of an accurate constitutive model difficult. One possible solution to this issue would be to define microstructural elements and perform fully coupled multiscale simulation. However, for complex geometries and loading scenarios, the computational costs of such simulations can be prohibitive. Ideally then, we should seek a method that contains microstructural detail, but leverages the speed of classical continuum-based finite-element (FE) modeling. In this work, we demonstrate the use of the Holzapfel-Gasser-Ogden (HGO) model [1, 2] to fit the behavior of microstructural network models. We show that Delaunay microstructural networks can be fit to the HGO strain energy function by calculating fiber network strain energy and average fiber stretch ratio. We then use the HGO constitutive model in a FE framework to improve the speed of our hybrid model, and demonstrate that this method, combined with a material property update scheme, can match a full multiscale simulation. This method gives us flexibility in defining complex FE simulations that would be impossible, or at least prohibitively time consuming, in multiscale simulation, while still accounting for microstructural heterogeneity.

Vascular Mechanobiology: Homeostasis, Adaptation, and Disease

Cells of the vascular wall are exquisitely sensitive to changes in their mechanical environment. In healthy vessels, mechanical forces regulate signaling and gene expression to direct the remodeling needed for the vessel wall to maintain optimal function. Major diseases of arteries involve maladaptive remodeling with compromised or lost homeostatic mechanisms. Whereas homeostasis invokes negative feedback loops at multiple scales to mediate mechanobiological stability, disease progression often occurs via positive feedback that generates mechanobiological instabilities. In this review, we focus on the cell biology, wall mechanics, and regulatory pathways associated with arterial health and how changes in these processes lead to disease. We discuss how positive feedback loops arise via biomechanical and biochemical means. We conclude that inflammation plays a central role in overriding homeostatic pathways and suggest future directions for addressing therapeutic needs.

Histological and biomechanical properties of systemic arteries in young and old Warmblood horses

Arterial rupture is a well-recognized cause of sudden death in horses, which mainly affects older horses. The arterial wall is known to stiffen with age, although the underlying age-related histological and biomechanical changes remain unclear. The purpose of this study was to investigate the effect of aging by histological analysis of the arterial wall and examination of the arterial wall biomechanical properties using an inflation-extension test. Entire circular samples of the proximal and distal aorta, cranial and caudal common carotid, external iliac, femoral and median artery were collected from 6 young (6 years) and 14 old horses (≥15 years). Samples of all arteries were histologically examined and intima media thickness as well as area % of elastin, smooth muscle actin and collagen type I and III were determined. Older horses had a significantly larger intima media thickness and a significantly higher area % of smooth muscle actin compared to young horses. Samples of the proximal and distal aorta, the caudal common carotid and the external iliac artery were mechanically assessed using an in-house developed inflation-extension device with ultrasound analysis. Rupture occurred in a minority of arteries (8/78) at high pressures (between 250-300 mmHg), and mostly occurred in older horses (7/8). Pressure-area, pressure-compliance and pressure-distensibility curves were constructed. A significant difference in the pressure-area curves of the distal aorta, common carotid artery and external iliac artery, the pressure-compliance curves of the proximal aorta and carotid artery and the pressure-distensibility curve of the proximal aorta was observed between young and old horses. Results demonstrate an effect of age on the histological and biomechanical properties of the arterial wall, which might explain why arterial rupture occurs more often in older horses.

Excessive adventitial stress drives inflammation-mediated fibrosis in hypertensive aortic remodelling in mice

Hypertension induces significant aortic remodelling, often adaptive but sometimes not. To identify immuno-mechanical mechanisms responsible for differential remodelling, we studied thoracic aortas from 129S6/SvEvTac and C57BL/6 J mice before and after continuous 14-day angiotensin II infusion, which elevated blood pressure similarly in both strains. Histological and biomechanical assessments of excised vessels were similar at baseline, suggesting a common homeostatic set-point for mean wall stress. Histology further revealed near mechano-adaptive remodelling of the hypertensive 129S6/SvEvTac aortas, but a grossly maladaptive remodelling of C57BL/6 J aortas. Bulk RNA sequencing suggested that increased smooth muscle contractile processes promoted mechano-adaptation of 129S6/SvEvTac aortas while immune processes prevented adaptation of C57BL/6 J aortas. Functional studies confirmed an increased vasoconstrictive capacity of the former while immunohistochemistry demonstrated marked increases in inflammatory cells in the latter. We then used multiple computational biomechanical models to test the hypothesis that excessive adventitial wall stress correlates with inflammatory cell infiltration. These models consistently predicted that increased vasoconstriction against an increased pressure coupled with modest deposition of new matrix thickens the wall appropriately, restoring wall stress towards homeostatic consistent with adaptive remodelling. By contrast, insufficient vasoconstriction permits high wall stresses and exuberant inflammation-driven matrix deposition, especially in the adventitia, reflecting compromised homeostasis and gross maladaptation.

A distributed lumped parameter model of blood flow with fluid-structure interaction

A distributed lumped parameter (DLP) model of blood flow was recently developed that can be simulated in minutes while still incorporating complex sources of energy dissipation in blood vessels. The aim of this work was to extend the previous DLP modeling framework to include fluid-structure interactions (DLP-FSI). This was done by using a simple compliance term to calculate pressure that does not increase the simulation complexity of the original DLP models. Verification and validation studies found DLP-FSI simulations had good agreement compared to analytical solutions of the wave equations, experimental measurements of pulsatile flow in elastic tubes, and in vivo MRI measurements of thoracic aortic flow. This new development of DLP-FSI allows for significantly improved computational efficiency of FSI simulations compared to FSI approaches that solve the full 3D conservation of mass and momentum equations while also including the complex sources of energy dissipation occurring in cardiovascular flows that other simplified models neglect.

Evolving structure-function relations during aortic maturation and aging revealed by multiphoton microscopy

The evolving microstructure and mechanical properties that promote homeostasis in the aorta are fundamental to age-specific adaptations and disease progression. We combine ex vivo multiphoton microscopy and biaxial biomechanical phenotyping to quantify and correlate layer-specific microstructural parameters, for the primary extracellular matrix components (fibrillar collagen and elastic lamellae) and cells (endothelial, smooth muscle, and adventitial), with mechanical properties of the mouse aorta from weaning through natural aging up to one year. The aging endothelium was characterized by progressive reductions in cell density and altered cellular orientation. The media similarly showed a progressive decrease in smooth muscle cell density and alignment though with inter-lamellar widening from intermediate to older ages, suggesting cell hypertrophy, matrix accumulation, or both. Despite not changing in tissue thickness, the aging adventitia exhibited a marked thickening and straightening of collagen fiber bundles and reduction in cell density, suggestive of age-related remodeling not growth. Multiple microstructural changes correlated with age-related increases in circumferential and axial material stiffness, among other mechanical metrics. Because of the importance of aging as a risk factor for cardiovascular diseases, understanding the normal progression of structural and functional changes is essential when evaluating superimposed disease-related changes as a function of the age of onset.

Potential damage in pulmonary arterial hypertension: An experimental study of pressure-induced damage of pulmonary artery

Pulmonary arterial hypertension (PAH) is associated with elevated pulmonary arterial pressure. PAH prognosis remains poor with a 15% mortality rate within 1 year, even with modern clinical management. Previous clinical studies proposed wall shear stress (WSS) to be an important hemodynamic factor affecting cell mechanotransduction, growth and remodeling, and disease progress in PAH. However, WSS in vivo is typically at most 2.5 Pa and a doubt has been cast whether WSS alone can drive disease progress. Furthermore, our current understanding of PAH pathology largely comes from small animals' studies in which caliber enlargement, a hallmark of PAH in humans, is rarely reported. Therefore, a large-animal experiment on pulmonary arteries (PAs) is needed to validate whether increased pressure can induce enlargement of PAs caliber. In this study, we use an inflation testing device to characterize the mechanical behavior, both nonlinear elastic behavior and irreversible damage of porcine arteries. The parameters of elastic behavior are estimated from the inflation test at a low-pressure range before and after over-pressurization. Then, histological images are qualitatively examined for medial and adventitial layers. This study sheds light on the relevance of pressure-induced damage mechanism in human PAH.

Structural and functional remodeling of the female Apoe-/- mouse aorta due to chronic cigarette smoke exposure

Despite a decline in popularity over the past several decades, cigarette smoking remains a leading cause of cardiovascular morbidity and mortality. Yet, the effects of cigarette smoking on vascular structure and function are largely unknown. To evaluate changes in the mechanical properties of the aorta that occur with chronic smoking, we exposed female apolipoprotein E-deficient mice to mainstream cigarette smoke daily for 24 wk, with room air as control. By the time of euthanasia, cigarette-exposed mice had lower body mass but experienced larger systolic/diastolic blood pressure when compared with controls. Smoking was associated with significant wall thickening, reduced axial stretch, and circumferential material softening of the aorta. Although this contributed to maintaining intrinsic tissue stiffness at control levels despite larger pressure loads, the structural stiffness became significantly larger. Furthermore, the aorta from cigarette-exposed mice exhibited decreased ability to store elastic energy and augment diastolic blood flow. Histological analysis revealed a region-dependent increase in the cross-sectional area due to smoking. Increased smooth muscle and extracellular matrix content led to medial thickening in the ascending aorta, whereas collagen deposition increased the thickness of the descending thoracic and abdominal aorta. Atherosclerotic lesions were larger in exposed vessels and featured a necrotic core overlaid by a thinned fibrous cap and macrophage infiltration, consistent with a vulnerable phenotype. Collectively, our data indicate that cigarette smoking decreases the mechanical functionality of the aorta, inflicts morphometric alterations to distinct segments of the aorta, and accelerates the progression of atherosclerosis.NEW & NOTEWORTHY We studied the effects of chronic cigarette smoking on the structure and function of the aorta in a mouse model of nose-only aerosol inhalation. Our data indicated that exposure to cigarette smoke impairs vascular function by reducing the ability of the aorta to store elastic energy and by decreasing aortic distensibility. Combined with a more vulnerable atherosclerotic phenotype, these findings reveal the biomechanical mechanisms that support the development of cardiovascular disease due to cigarette smoking.

Differential biomechanical responses of elastic and muscular arteries to angiotensin II-induced hypertension

Elastic and muscular arteries are distinguished by their distinct microstructures, biomechanical properties, and smooth muscle cell contractile functions. They also exhibit differential remodeling in aging and hypertension. Although regional differences in biomechanical properties have been compared, few studies have quantified biaxial differences in response to hypertension. Here, we contrast passive and active changes in large elastic and medium- and small-sized muscular arteries in adult mice in response to chronic infusion of angiotensin over 14 days. We found a significant increase in wall thickness, both medial and adventitial, in the descending thoracic aorta that associated with trends of an increased collagen:elastin ratio. There was adventitial thickening in the small-sized mesenteric artery, but also significant changes in elastic lamellar structure and contractility. An increased contractile response to phenylephrine coupled with a reduced vasodilatory response to acetylcholine in the mesenteric artery suggested an increased contractile state in response to hypertension. Overall reductions in the calculated gradients in pulse wave velocity and elastin energy storage capability from elastic-to-muscular arteries suggested a possible transfer of excessive pulsatile energy into the small-sized muscular arteries resulting in significant functional consequences in response to hypertension.

Constitutive modeling using structural information on collagen fiber direction and dispersion in human superficial femoral artery specimens of different ages

Arterial mechanics plays an important role in vascular pathophysiology and repair, and advanced imaging can inform constitutive models of vascular behavior. We have measured the mechanical properties of 14 human superficial femoral arteries (SFAs) (age 12-70, mean 48±19 years) using planar biaxial extension, and determined the preferred collagen fiber direction and dispersion using multiphoton microscopy. The collagen fiber direction and dispersion were evaluated using second-harmonic generation imaging and modeled using bivariate von Mises distributions. The microstructures of elastin and collagen were assessed using two-photon fluorescence imaging and conventional bidirectional histology. The mechanical and structural data were used to describe the SFA mechanical behavior using two- and four-fiber family invariant-based constitutive models. Older SFAs were stiffer and mechanically more nonlinear than younger specimens. In the adventitia, collagen fibers were undulated and diagonally-oriented, while in the media, they were straight and circumferentially-oriented. The media was rich in collagen that surrounded the circumferentially-oriented smooth muscle cells, and the elastin was present primarily in the internal and external elastic laminae. Older SFAs had a more circumferential collagen fiber alignment, a decreased circumferential-radial fiber dispersion, but the same circumferential-longitudinal fiber dispersion as younger specimens. Both the two- and the four-fiber family constitutive models were able to capture the experimental data, and the fits were better for the four-fiber family formulation. Our data provide additional details on the SFA intramural structure and inform structurally-based constitutive models.

Insights on the Pathogenesis of Aneurysm through the Study of Hereditary Aortopathies

Thoracic aortic aneurysms (TAA) are permanent and localized dilations of the aorta that predispose patients to a life-threatening risk of aortic dissection or rupture. The identification of pathogenic variants that cause hereditary forms of TAA has delineated fundamental molecular processes required to maintain aortic homeostasis. Vascular smooth muscle cells (VSMCs) elaborate and remodel the extracellular matrix (ECM) in response to mechanical and biochemical cues from their environment. Causal variants for hereditary forms of aneurysm compromise the function of gene products involved in the transmission or interpretation of these signals, initiating processes that eventually lead to degeneration and mechanical failure of the vessel. These include mutations that interfere with transduction of stimuli from the matrix to the actin-myosin cytoskeleton through integrins, and those that impair signaling pathways activated by transforming growth factor-β (TGF-β). In this review, we summarize the features of the healthy aortic wall, the major pathways involved in the modulation of VSMC phenotypes, and the basic molecular functions impaired by TAA-associated mutations. We also discuss how the heterogeneity and balance of adaptive and maladaptive responses to the initial genetic insult might contribute to disease.

A Chemomechanobiological Model of the Long-Term Healing Response of Arterial Tissue to a Clamping Injury

Vascular clamping often causes injury to arterial tissue, leading to a cascade of cellular and extracellular events. A reliable in silico prediction of these processes following vascular injury could help us to increase our understanding thereof, and eventually optimize surgical techniques or drug delivery to minimize the amount of long-term damage. However, the complexity and interdependency of these events make translation into constitutive laws and their numerical implementation particularly challenging. We introduce a finite element simulation of arterial clamping taking into account acute endothelial denudation, damage to extracellular matrix, and smooth muscle cell loss. The model captures how this causes tissue inflammation and deviation from mechanical homeostasis, both triggering vascular remodeling. A number of cellular processes are modeled, aiming at restoring this homeostasis, i.e., smooth muscle cell phenotype switching, proliferation, migration, and the production of extracellular matrix. We calibrated these damage and remodeling laws by comparing our numerical results to in vivo experimental data of clamping and healing experiments. In these same experiments, the functional integrity of the tissue was assessed through myograph tests, which were also reproduced in the present study through a novel model for vasodilator and -constrictor dependent smooth muscle contraction. The simulation results show a good agreement with the in vivo experiments. The computational model was then also used to simulate healing beyond the duration of the experiments in order to exploit the benefits of computational model predictions. These results showed a significant sensitivity to model parameters related to smooth muscle cell phenotypes, highlighting the pressing need to further elucidate the biological processes of smooth muscle cell phenotypic switching in the future.

Regulation of SMC traction forces in human aortic thoracic aneurysms

Smooth muscle cells (SMCs) usually express a contractile phenotype in the healthy aorta. However, aortic SMCs have the ability to undergo profound changes in phenotype in response to changes in their extracellular environment, as occurs in ascending thoracic aortic aneurysms (ATAA). Accordingly, there is a pressing need to quantify the mechanobiological effects of these changes at single cell level. To address this need, we applied Traction Force Microscopy (TFM) on 759 cells coming from three primary healthy (AoPrim) human SMC lineages and three primary aneurysmal (AnevPrim) human SMC lineages, from age and gender matched donors. We measured the basal traction forces applied by each of these cells onto compliant hydrogels of different stiffness (4, 8, 12, 25 kPa). Although the range of force generation by SMCs suggested some heterogeneity, we observed that: 1. the traction forces were significantly larger on substrates of larger stiffness; 2. traction forces in AnevPrim were significantly higher than in AoPrim cells. We modelled computationally the dynamic force generation process in SMCs using the motor-clutch model and found that it accounts well for the stiffness-dependent traction forces. The existence of larger traction forces in the AnevPrim SMCs were related to the larger size of cells in these lineages. We conclude that phenotype changes occurring in ATAA, which were previously known to reduce the expression of elongated and contractile SMCs (rendering SMCs less responsive to vasoactive agents), tend also to induce stronger SMCs. Future work aims at understanding the causes of this alteration process in aortic aneurysms.

Adventitial remodeling protects against aortic rupture following late smooth muscle-specific disruption of TGFβ signaling

Altered signaling through transforming growth factor-beta (TGFβ) increases the risk of aortic dissection in patients, which has been confirmed in mouse models. It is well known that altered TGFβ signaling affects matrix turnover, but there has not been a careful examination of associated changes in structure-function relations. In this paper, we present new findings on the rupture potential of the aortic wall following late postnatal smooth muscle cell (SMC)-specific disruption of type I and II TGFβ receptors in a mouse model with demonstrated dissection susceptibility. Using a combination of custom computer-controlled biaxial tests and quantitative histology and immunohistochemistry, we found that loss of TGFβ signaling in SMCs compromises medial properties but induces compensatory changes in the adventitia that preserve wall strength above that which is needed to resist in vivo values of wall stress. These findings emphasize the different structural defects that lead to aortic dissection and rupture - compromised medial integrity and insufficient adventitial strength, respectively. Relative differences in these two defects, in an individual subject at a particular time, likely reflects the considerable phenotypic diversity that is common in clinical presentations of thoracic aortic dissection and rupture. There is, therefore, a need to move beyond examinations of bulk biological assays and wall properties to cell- and layer-specific studies that delineate pathologic and compensatory changes in wall biology and composition, and thus the structural integrity of the aortic wall that can dictate differences between life and death.

From Transcript to Tissue: Multiscale Modeling from Cell Signaling to Matrix Remodeling

Tissue-level biomechanical properties and function derive from underlying cell signaling, which regulates mass deposition, organization, and removal. Here, we couple two existing modeling frameworks to capture associated multiscale interactions-one for vessel-level growth and remodeling and one for cell-level signaling-and illustrate utility by simulating aortic remodeling. At the vessel level, we employ a constrained mixture model describing turnover of individual wall constituents (elastin, intramural cells, and collagen), which has proven useful in predicting diverse adaptations as well as disease progression using phenomenological constitutive relations. Nevertheless, we now seek an improved mechanistic understanding of these processes; we replace phenomenological relations in the mixture model with a logic-based signaling model, which yields a system of ordinary differential equations predicting changes in collagen synthesis, matrix metalloproteinases, and cell proliferation in response to altered intramural stress, wall shear stress, and exogenous angiotensin II. This coupled approach promises improved understanding of the role of cell signaling in achieving tissue homeostasis and allows us to model feedback between vessel mechanics and cell signaling. We verify our model predictions against data from the hypertensive murine infrarenal abdominal aorta as well as results from validated phenomenological models, and consider effects of noisy signaling and heterogeneous cell populations.

Computational Study of Growth and Remodeling in Ascending Thoracic Aortic Aneurysms Considering Variations of Smooth Muscle Cell Basal Tone

In this paper, we investigate the progression of Ascending Thoracic Aortic Aneurysms (ATAA) using a computational model of Growth and Remodeling (G&R) taking into account the composite (elastin, four collagen fiber families and Smooth Muscle Cells—SMCs) and multi-layered (media and adventitia) nature of the aorta. The G&R model, which is based on the homogenized Constrained Mixture theory, is implemented as a UMAT in the Abaqus finite-element package. Each component of the mixture is assigned a strain energy density function: nearly-incompressible neo-Hookean for elastin and Fung-type for collagen and SMCs. Active SMCs tension is additionally considered, through a length-tension relationship having a classic inverted parabola shape, in order to investigate its effects on the progression of ATAA in a patient-specific model. A sensitivity analysis is performed to evaluate the potential impact of variations in the parameters of the length-tension relationships. These variations reflect in variations of SMCs normal tone during ATAA progression, with active stress contributions ranging between 30% (best case scenario) and 0% (worst case scenario) of the total wall circumferential stress. Low SMCs active stress in the worst case scenarios, in fact, affect the rates of collagen deposition by which the elastin loss is gradually compensated by collagen deposition in the simulated ATAA progression, resulting eventually in larger aneurysm diameters. The types of length-tension relationships leading to a drop of SMCs active stress in our simulations reveal a critical condition which could also result in SMCs apoptosis.

Reconstructing vascular homeostasis by growth-based prestretch and optimal fiber deposition

Computational modeling of cardiovascular biomechanics should generally start from a homeostatic state. This is particularly relevant for image-based modeling, where the reference configuration is the loaded in vivo state obtained from imaging. This state includes residual stress of the vascular constituents, as well as anisotropy from the spatially varying orientation of collagen and smooth muscle fibers. Estimation of the residual stress and fiber orientation fields is a formidable challenge in realistic applications. To help address this challenge, we herein develop a growth based Algorithm to recover a residual stress distribution in vascular domains such that the stress state in the loaded configuration is equal to a prescribed homeostatic stress distribution at physiologic pressure. A stress-driven fiber deposition process is included in the framework, which defines the distribution of the fiber alignments in the vascular homeostatic state based on a minimization procedure. Numerical simulations are conducted to test this two-stage homeostasis generation algorithm in both idealized and non-idealized geometries, yielding results that agree favorably with prior numerical and experimental data.

Electrospun Tissue-Engineered Arterial Graft Thickness Affects Long-Term Composition and Mechanics

Wall stress is often lower in tissue-engineered constructs than in comparable native tissues due to the use of stiff polymeric materials having thicker walls. In this work, we sought to design a murine arterial graft having a more favorable local mechanical environment for the infiltrating cells; we used electrospinning to enclose a compliant inner core of poly(glycerol sebacate) with a stiffer sheath of poly(caprolactone) to reduce the potential for rupture. Two scaffolds were designed that differed in the thickness of the core as previous computational simulations found that circumferential wall stresses could be increased in the core toward native values by increasing the ratio of the core:sheath. Our modified electrospinning protocols reduced swelling of the core upon implantation and eliminated residual stresses in the sheath, both of which had contributed to the occlusion of implanted grafts during pilot studies. For both designs, a subset of implanted grafts occluded due to thrombosis or ruptured due to suspected point defects in the sheath. However, there were design-based differences in collagen content and mechanical behavior during early remodeling of the patent samples, with the thinner-core scaffolds having more collagen and a stiffer behavior after 12 weeks of implantation than the thicker-core scaffolds. By 24 weeks, the thicker-core scaffolds also became stiff, with similar amounts of collagen but increased smooth muscle cell and elastin content. These data suggest that increasing wall stress toward native values may provide a more favorable environment for normal arterial constituents to form despite the overall stiffness of the construct remaining elevated due to the absolute increase in load-bearing constituents.

Vascular adaptation in the presence of external support - A modeling study

Vascular grafts have long been used to replace damaged or diseased vessels with considerable success, but a new approach is emerging where native vessels are merely supported, not replaced. Although external supports have been evaluated in diverse situations - ranging from aneurysmal disease to vein grafts or the Ross operation - optimal supports and procedures remain wanting. In this paper, we present a novel application of a growth and remodeling model well suited for parametrically exploring multiple designs of external supports while accounting for mechanobiological and immunobiological responses of the supported native vessel. These results suggest that a load bearing external support can reduce vessel thickening in response to pressure elevation. Results also suggest that the final adaptive state of the vessel depends on the structural stiffness of the support via a mechano-driven adaptation, although luminal encroachment may be a complication in the presence of chronic inflammation. Finally, the supported vessel can stiffen (structurally and materially) along circumferential and axial directions, which could have implications on overall hemodynamics and thus subsequent vascular remodeling. The proposed framework can provide valuable insights into vascular adaptation in the presence of external support, accelerate rational design, and aid translation of this emerging approach.

Unraveling the multilayer mechanical response of aorta using layer-specific residual stresses and experimental properties

To test the capability of the multilayer model, we used previously published layer-specific experimental data relating to the axial pre-stretch, the opening angle, the fiber distribution obtained by polarized light microscopy measurements, and the uniaxial and biaxial response of the porcine descending and abdominal aorta. We fitted the mechanical behavior of each arterial layer using Gasser, Holzapfel and Ogden strain energy function using the dispersion parameter κ as phenomenological parameter obtained during the fitting procedure or computed from the experimental fiber distribution. A multilayer finite element model of the whole aorta with the dimensions of the circumferential and longitudinal strips were then built using layer-specific material parameters previously fitted. This model was used to capture the whole aorta response under uniaxial and biaxial stress states and to reproduce the response of the whole aorta to internal pressure. Our results show that a model based on a multilayer structure without residual stresses is unable to render the uniaxial and biaxial mechanical response of the aorta (R²=0.6954 and R²=0.8582 for descending thoracic aorta (DTA) and infrarenal abdominal aorta (IAA), respectively). Only an appropriate multilayer model that includes layer-specific residual stresses can reproduce the response of the whole aorta (R²=0.9787 and R²=0.9636 for DTA and IAA respectively). In addition, a multilayer model without residual stresses produces the same stress distribution as a monolayer model without residual stresses where the maximal value of circumferential and longitudinal stresses appears at the inner radius of the intima. Finally, if layer-specific residual stresses are not available, there is less error the stress distribution using a monolayer model with residual stresses that a multilayer model without residual stresses.

Fast, Rate-Independent, Finite Element Implementation of a 3D Constrained Mixture Model of Soft Tissue Growth and Remodeling

Constrained mixture models of soft tissue growth and remodeling can simulate many evolving conditions in health as well as in disease and its treatment, but they can be computationally expensive. In this paper, we derive a new fast, robust finite element implementation based on a concept of mechanobiological equilibrium that yields fully resolved solutions and allows computation of quasi-equilibrated evolutions when imposed perturbations are slow relative to the adaptive process. We demonstrate quadratic convergence and verify the model via comparisons with semi-analytical solutions for arterial mechanics. We further examine the enlargement of aortic aneurysms for which we identify new mechanobiological insights into factors that affect the nearby non-aneurysmal segment as it responds to the changing mechanics within the diseased segment. Because this new 3D approach can be implemented within many existing finite element solvers, constrained mixture models of growth and remodeling can now be used more widely.

Evaluation of the Stress-Growth Hypothesis in Saphenous Vein Perfusion Culture

The great saphenous vein (GSV) has served as a coronary artery bypass graft (CABG) conduit for over 50 years. Despite prevalent use, first-year failure rates remain high compared to arterial autograft options. Amongst other factors, vein graft failure can be attributed to material and mechanical mismatching that lead to apoptosis, inflammation, and intimal-medial hyperplasia. Through the implementation of the continuum mechanical-based theory of "stress-mediated growth and remodeling," we hypothesize that the mechanical properties of porcine GSV grafts can be favorably tuned for CABG applications prior to implantation using a prolonged but gradual transition from venous to arterial loading conditions in an inflammatory and thrombogenic deficient environment. To test this hypothesis, we used a hemodynamic-mimetic perfusion bioreactor to guide remodeling through stepwise incremental changes in pressure and flow over the course of 21-day cultures. Biaxial mechanical testing of vessels pre- and post-remodeling was performed, with results fit to structurally-motivated constitutive models using non-parametric bootstrapping. The theory of "small-on-large" was used to describe appropriate stiffness moduli, while histology and viability assays confirmed microstructural adaptations and vessel viability. Results suggest that stepwise transition from venous-to-arterial conditions results in a partial restoration of circumferential stretch and circumferential, but not axial, stress through vessel dilation and wall thickening in a primarily outward remodeling process. These remodeled tissues also exhibited decreased mechanical isotropy and circumferential, but not axial, stiffening. In contrast, only increases in axial stiffness were observed using culture under venous perfusion conditions and those tissues experienced moderate intimal resorption.