Interstitial flow, pressure and residual stress in the aging carotid artery model in FEBio.

Vascular smooth muscle cells (VSMCs) are subject to interstitial flow-induced shear stress, which is a critical parameter in cardiovascular disease progression. Transmural pressure loading and residual stresses alter the hydraulic conductivity of the arterial layers and modulate the interstitial fluid flux through the arterial wall. In this paper, a biphasic multilayer model of a common carotid artery (CCA) with anisotropic fiber-reinforced soft tissue and strain-dependent permeability is developed in FEBio software. After the verification of the numerical predictions, age-related arterial thickening and stiffening effects on arterial deformation and interstitial flow are computed under physiological geometry and physical parameters. We found that circumferential residual stress shifts outward in each layer and its gradient increases up to 6 times with aging. Internally pressurized CCA displays nonlinear deformation. In the aged artery, the circumferential stress becomes greater on the media layer (82-158 kPa) and lower on the intima and adventitia (19-23 kPa and 25-28 kPa, respectively). The radial compression of the intima reduces the total hydraulic conductivity by 48% in the young and 16% in the aged arterial walls. Consequently, the average radial interstitial flux increases with pressure by 14% in the young and 91% in the aged arteries. Accordingly, the flow shear stress experienced by the VSMCs becomes more significant for aged arteries, which may accelerate cardiovascular disease progression compared to young arteries.

Myocardial Biomechanics and the Consequent Differentially Expressed Genes of the Left Atrial Ligation Chick Embryonic Model of Hypoplastic Left Heart Syndrome.

Left atrial ligation (LAL) of the chick embryonic heart is a model of the hypoplastic left heart syndrome (HLHS) where a purely mechanical intervention without genetic or pharmacological manipulation is employed to initiate cardiac malformation. It is thus a key model for understanding the biomechanical origins of HLHS. However, its myocardial mechanics and subsequent gene expressions are not well-understood. We performed finite element (FE) modeling and single-cell RNA sequencing to address this. 4D high-frequency ultrasound imaging of chick embryonic hearts at HH25 (ED 4.5) were obtained for both LAL and control. Motion tracking was performed to quantify strains. Image-based FE modeling was conducted, using the direction of the smallest strain eigenvector as the orientations of contractions, the Guccione active tension model and a Fung-type transversely isotropic passive stiffness model that was determined via micro-pipette aspiration. Single-cell RNA sequencing of left ventricle (LV) heart tissues was performed for normal and LAL embryos at HH30 (ED 6.5) and differentially expressed genes (DEG) were identified.After LAL, LV thickness increased by 33%, strains in the myofiber direction increased by 42%, while stresses in the myofiber direction decreased by 50%. These were likely related to the reduction in ventricular preload and underloading of the LV due to LAL. RNA-seq data revealed potentially related DEG in myocytes, including mechano-sensing genes (Cadherins, NOTCH1, etc.), myosin contractility genes (MLCK, MLCP, etc.), calcium signaling genes (PI3K, PMCA, etc.), and genes related to fibrosis and fibroelastosis (TGF-ß, BMP, etc.). We elucidated the changes to the myocardial biomechanics brought by LAL and the corresponding changes to myocyte gene expressions. These data may be useful in identifying the mechanobiological pathways of HLHS.

Embryonic aortic arch material properties obtained by optical coherence tomography-guided micropipette aspiration.

It is challenging to determine the in vivo material properties of a very soft, mesoscale arterial vesselsof size âˆ¼ 80 to 120 µm diameter. This information is essential to understand the early embryonic cardiovascular development featuring rapidly evolving dynamic microstructure. Previous research efforts to describe the properties of the embryonic great vessels are very limited. Our objective is to measure the local material properties of pharyngeal aortic arch tissue of the chick-embryo during the early Hamburger-Hamilton (HH) stages, HH18 and HH24. Integrating the micropipette aspiration technique with optical coherence tomography (OCT) imaging, a clear vision of the aspirated arch geometry is achieved for an inner pipette radius of Rp = 25 µm. The aspiration of this region is performed through a calibrated negatively pressurized micro-pipette. A computational finite element model is developed to model the nonlinear behaviour of the arch structure by considering the geometry-dependent constraints. Numerical estimations of the nonlinear material parameters for aortic arch samples are presented. The exponential material nonlinearity parameter (a) of aortic arch tissue increases statistically significantly from a = 0.068 ± 0.013 at HH18 to a = 0.260 ± 0.014 at HH24 (p = 0.0286). As such, the aspirated tissue length decreases from 53 µm at HH18 to 34 µm at HH24. The calculated NeoHookean shear modulus increases from 51 Pa at HH18 to 93 Pa at HH24 which indicates a statistically significant stiffness increase. These changes are due to the dynamic changes of collagen and elastin content in the media layer of the vessel during development.

A novel Fontan Y-graft for interrupted inferior vena cava and azygos continuation.

OBJECTIVES: To evaluate the hemodynamicdynamic advantage of a new Fontan surgical template that is intended for complex single-ventricle patients with interrupted inferior vena cava-azygos and hemi-azygos continuation. The new technique has emerged from a comprehensive pre-surgical simulation campaign conducted to facilitate a balanced hepatic flow and somatic Fontan pathway growth after Kawashima procedure. METHODS: For 9 patients, aged 2 to18 years, majority having poor preoperative oxygen saturation, a pre-surgical computational fluid dynamics customization is conducted. Both the traditional Fontan pathways and the proposed novel Y-graft templates are considered. Numerical model was validated against in vivo phase-contrast magnetic resonance imaging data and in vitro experiments. RESULTS: The proposed template is selected and executed for 6 out of the 9 patients based on its predicted superior hemodynamic performance. Pre-surgical simulations performed for this cohort indicated that flow from the hepatic veins (HEP) do not reach to the desired lung. The novel Y-graft template, customized via a right- or left-sided displacement of the total cavopulmonary connection anastomosis location resulted a drastic increase in HEP flow to the desired lung. Orientation of HEP to azygos direct shunt is found to be important as it can alter the flow pattern from 38% in the caudally located direct shunt to 3% in the cranial configuration with significantly reversed flow. The postoperative measurements prove that oxygen saturation increased significantly (P-value = 0.00009) to normal levels in 1 year follow-up. CONCLUSIONS: The new Y-graft template, if customized for the individual patient, is a viable alternative to the traditional surgical pathways. This template addresses the competing hemodynamic design factors of low physiological venous pressure, high postoperative oxygen saturation, low energy loss and balanced hepatic growth factor distribution possibly assuring adequate lung development. DATE AND NUMBER OF IRB APPROVAL: 25 October 2019, 280011928-604.01.01.

Computational modeling of vascular growth in patient-specific pulmonary arterial patch reconstructions.

Recent progress in vascular growth mechanics has involved the use of computational algorithms to address clinical problems with the use of three-dimensional patient specific geometries. The objective of this study is to establish a predictive computational model for the volumetric growth of pulmonary arterial (PA) tissue following complex cardiovascular patch reconstructive surgeries for congenital heart disease patients. For the first time in the literature, the growth mechanics and performance of artificial cardiovascular patches in contact with the growing PA tissue domain is established. An elastic-growing material model was developed in the open source FEBio software suite to first examine the surgical patch reconstruction process for an idealized main PA anatomy as a benchmark model and then for the patient-specific PA of a newborn. Following patch reconstruction, high levels of stress and strain are compensated by growth on the arterial tissue. As this growth progresses, the arterial tissue is predicted to stiffen to limit elastic deformations. We simulated this arterial growth up to the age of 18 years, when somatic growth plateaus. Our research findings show that the non-growing patch material remains in a low strain state throughout the simulation timeline, while experiencing high stress hot-spots. Arterial tissue growth along the surgical stitch lines is triggered mainly due to PA geometry and blood pressure, rather than due to material property differences in the artificial and native tissue. Thus, non-uniform growth patterns are observed along the arterial tissue proximal to the sutured boundaries. This computational approach is effective for the pre-surgical planning of complex patch surgeries to quantify the unbalanced growth of native arteries and artificial non-growing materials to develop optimal patch biomechanics for improved postoperative outcomes.

Microstructure of early embryonic aortic arch and its reversibility following mechanically altered hemodynamic load release.

In the embryonic heart, blood flow is distributed through a bilaterally paired artery system composed of the aortic arches (AAs). The purpose of this study is to establish an understanding of the governing mechanism of microstructural maturation of the AA matrix and its reversibility, toward the desired macroscopic vessel lumen diameter and thickness for healthy, abnormal, and in ovo repaired abnormal mechanical loading. While matrix-remodeling mechanisms were significantly different for normal versus conotruncal banding (CTB), both led to an increase in vessel lumen. Correlated with right-sided flow increase at Hamburger & Hamilton stages 21, intermittent load switching between collagen I and III with elastin and collagen-IV defines the normal process. However, decreases in collagen I, elastin, vascular endothelial growth factor-A, and fibrillin-1 in CTB were recovered almost fully following the CTB-release model, primarily because of the pressure load changes. The complex temporal changes in matrix proteins are illustrated through a predictive finite-element model based on elastin and collagen load-sharing mechanism to achieve lumen area increase and thickness increase resulting from wall shear stress and tissue strain, respectively. The effect of embryonic timing in cardiac interventions on AA microstructure was established where abnormal mechanical loading was selectively restored at the key stage of development. Recovery of the normal mechanical loading via early fetal intervention resulted in delayed microstructural maturation. Temporal elastin increase, correlated with wall shear stress, is required for continuous lumen area growth.NEW & NOTEWORTHY The present study undertakes comparative analyses of the mechanistic differences of the arterial matrix microstructure and dynamics in the three fundamental processes of control, conotruncal banded, and released conotruncal band in avian embryo. Among other findings, this study provides specific evidence on the restorative role of elastin during the early lumen growth process. During vascular development, a novel intermittent load-switching mechanism between elastin and collagen, triggered by a step increase in wall shear stress, governs the chronic vessel lumen cross-sectional area increase. Mimicking the fetal cardiovascular interventions currently performed in humans, the early release of the abnormal mechanical load rescues the arterial microstructure with time lag.

Spatiotemporal remodeling of embryonic aortic arch: stress distribution, microstructure, and vascular growth in silico.

The microstructure for mature vessels has been investigated in detail, while there is limited information about the embryonic stages, in spite of their importance in the prognosis of congenital heart defects. It is hypothesized that the embryonic vasculature represents a disorganized but dynamic soft tissue, which rapidly evolves toward a specialized multi-cellular vascular structure under mechanical loading. Here the microstructural evolution process of the embryonic pharyngeal aortic arch structure was simulated using an in ovo validated long-term growth and remodeling computational model, implemented as an in-house FEBio plug-in. Optical coherence tomography-guided servo-null pressure measurements are assigned as boundary conditions through the critical embryonic stages. The accumulation of key microstructural constituents was recorded through zoom confocal microscopy for all six embryonic arch arteries simultaneously. The total amount and the radial variation slope of the collagen along the arch wall thickness in different arch types and for different embryonic times, with different dimension scales, were normalized and compared statistically. The arch growth model shows that the stress levels around the lumen boundary increase from [Formula: see text] (Stage 18) to a level higher than [Formula: see text] (Stage 24), depending on matrix constituent production rates, while the homeostatic strain level is kept constant. The statistical tests show that although the total collagen levels differentiate among bilateral positions of the same arch, the shape coefficient of the matrix microstructural gradient changes with embryonic time, proving radial localization, in accordance with numerical model results. In vivo cell number (DAPI) and vascular endothelial growth factor distributions followed similar trends.

Computational Pre-surgical Planning of Arterial Patch Reconstruction: Parametric Limits and In Vitro Validation.

Surgical treatment of congenital heart disease (CHD) involves complex vascular reconstructions utilizing artificial and native surgical materials. A successful surgical reconstruction achieves an optimal hemodynamic profile through the graft in spite of the complex post-operative vessel growth pattern and the altered pressure loading. This paper proposes a new in silico patient-specific pre-surgical planning framework for patch reconstruction and investigates its computational feasibility. The proposed protocol is applied to the patch repair of main pulmonary artery (MPA) stenosis in the Tetralogy of Fallot CHD template. The effects of stenosis grade, the three-dimensional (3D) shape of the surgical incision and material properties of the artificial patch are investigated. The release of residual stresses due to the surgical incision and the extra opening of the incision gap for patch implantation are simulated through a quasi-static finite-element vascular model with shell elements. Implantation of different unloaded patch shapes is simulated. The patched PA configuration is pressurized to the physiological post-operative blood pressure levels of 25 and 45 mmHg and the consequent post-operative stress distributions and patched artery shapes are computed. Stress-strain data obtained in-house, through the biaxial tensile tests for the mechanical properties of common surgical patch materials, Dacron, Polytetrafluoroethylene, human pericardium and porcine xenopericardium, are employed to represent the mechanical behavior of the patch material. Finite-element model is experimentally validated through the actual patch surgery reconstructions performed on the 3D printed anatomical stenosis replicas. The post-operative recovery of the initially narrowed lumen area and post-op tortuosity are quantified for all modeled cases. A computational fluid dynamics solver is used to evaluate post-operative pressure drop through the patch-reconstructed outflow tract. According to our findings, the shorter incisions made at the throat result in relatively low local peak stress values compared to other patch design alternatives. Longer cut and double patch cases are the most effective in repairing the initial stenosis level. After the patch insertion, the pressure drop in the artery due to blood flow decreases from 9.8 to 1.35 mmHg in the conventional surgical configuration. These results are in line with the clinical experience where a pressure gradient at or above 50 mmHg through the MPA can be an indication to intervene. The main strength of the proposed pre-surgical planning framework is its capability to predict the intra-operative and post-operative 3D vascular shape changes due to intramural pressure, cut length and configuration, for both artificial and native patch materials.