Ex Vivo Evidence for the Contribution of Hemodynamic Shear Stress Abnormalities to the Early Pathogenesis of Calcific Bicuspid Aortic Valve Disease

The bicuspid aortic valve (BAV) is the most common congenital cardiac anomaly and is frequently associated with calcific aortic valve disease (CAVD). The most prevalent type-I morphology, which results from left-/right-coronary cusp fusion, generates different hemodynamics than a tricuspid aortic valve (TAV). While valvular calcification has been linked to genetic and atherogenic predispositions, hemodynamic abnormalities are increasingly pointed as potential pathogenic contributors. In particular, the wall shear stress (WSS) produced by blood flow on the leaflets regulates homeostasis in the TAV. In contrast, WSS alterations cause valve dysfunction and disease. While such observations support the existence of synergies between valvular hemodynamics and biology, the role played by BAV WSS in valvular calcification remains unknown. The objective of this study was to isolate the acute effects of native BAV WSS abnormalities on CAVD pathogenesis. Porcine aortic valve leaflets were subjected ex vivo to the native WSS experienced by TAV and type-I BAV leaflets for 48 hours. Immunostaining, immunoblotting and zymography were performed to characterize endothelial activation, pro-inflammatory paracrine signaling, extracellular matrix remodeling and markers involved in valvular interstitial cell activation and osteogenesis. While TAV and non-coronary BAV leaflet WSS essentially maintained valvular homeostasis, fused BAV leaflet WSS promoted fibrosa endothelial activation, paracrine signaling (2.4-fold and 3.7-fold increase in BMP-4 and TGF-β1, respectively, relative to fresh controls), catabolic enzyme secretion (6.3-fold, 16.8-fold, 11.7-fold, 16.7-fold and 5.5-fold increase in MMP-2, MMP-9, cathepsin L, cathepsin S and TIMP-2, respectively) and activity (1.7-fold and 2.4-fold increase in MMP-2 and MMP-9 activity, respectively), and bone matrix synthesis (5-fold increase in osteocalcin). In contrast, BAV WSS did not significantly affect α-SMA and Runx2 expressions and TIMP/MMP ratio. This study demonstrates the key role played by BAV hemodynamic abnormalities in CAVD pathogenesis and suggests the dependence of BAV vulnerability to calcification on the local degree of WSS abnormality.     

The congenital bicuspid aortic valve can experience high-frequency unsteady shear stresses on its leaflet surface

The bicuspid aortic valve (BAV) is a common congenital malformation of the aortic valve (AV) affecting 1% to 2% of the population. The BAV is predisposed to early degenerative calcification of valve leaflets, and BAV patients constitute 50% of AV stenosis patients. Although evidence shows that genetic defects can play a role in calcification of the BAV leaflets, we hypothesize that drastic changes in the mechanical environment of the BAV elicit pathological responses from the valve and might be concurrently responsible for early calcification. An in vitro model of the BAV was constructed by surgically manipulating a native trileaflet porcine AV. The BAV valve model and a trileaflet AV (TAV) model were tested in an in vitro pulsatile flow loop mimicking physiological hemodynamics. Laser Doppler velocimetry was used to make measurements of fluid shear stresses on the leaflet of the valve models using previously established methodologies. Furthermore, particle image velocimetry was used to visualize the flow fields downstream of the valves and in the sinuses. In the BAV model, flow near the leaflets and fluid shear stresses on the leaflets were much more unsteady than for the TAV model, most likely due to the moderate stenosis in the BAV and the skewed forward flow jet that collided with the aorta wall. This additional unsteadiness occurred during mid- to late-systole and was composed of cycle-to-cycle magnitude variability as well as high-frequency fluctuations about the mean shear stress. It has been demonstrated that the BAV geometry can lead to unsteady shear stresses under physiological flow and pressure conditions. Such altered shear stresses could play a role in accelerated calcification in BAVs.     

Three-dimensional macro-scale assessment of regional and temporal wall shear stress characteristics on aortic valve leaflets

The aortic valve (AV) achieves unidirectional blood flow between the left ventricle and the aorta. Although hemodynamic stresses have been shown to regulate valvular biology, the native wall shear stress (WSS) experienced by AV leaflets remains largely unknown. The objective of this study was to quantify computationally the macro-scale leaflet WSS environment using fluid–structure interaction modeling. An arbitrary Lagrangian–Eulerian approach was implemented to predict valvular flow and leaflet dynamics in a three-dimensional AV geometry subjected to physiologic transvalvular pressure. Local WSS characteristics were quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI) and temporal shear gradient (TSG). The dominant radial WSS predicted on the leaflets exhibited high amplitude and unidirectionality on the ventricularis (TSM>7.50 dyn/cm², OSI < 0.17, TSG>325.54 dyn/cm² s) but low amplitude and bidirectionality on the fibrosa (TSM < 2.73 dyn/cm², OSI>0.38, TSG < 191.17 dyn/cm² s). The radial WSS component computed in the leaflet base, belly and tip demonstrated strong regional variability (ventricularis TSM: 7.50–22.32 dyn/cm², fibrosa TSM: 1.26–2.73 dyn/cm²). While the circumferential WSS exhibited similar spatially dependent magnitude (ventricularis TSM: 1.41–3.40 dyn/cm², fibrosa TSM: 0.42–0.76 dyn/cm²) and side-specific amplitude (ventricularis TSG: 101.73–184.43 dyn/cm² s, fibrosa TSG: 41.92–54.10 dyn/cm² s), its temporal variations were consistently bidirectional (OSI>0.25). This study provides new insights into the role played by leaflet–blood flow interactions in valvular function and critical hemodynamic stress data for the assessment of the hemodynamic theory of AV disease.     

A multiscale computational comparison of the bicuspid and tricuspid aortic valves in relation to calcific aortic stenosis

Patients with bicuspid aortic valve (BAV) are more likely to develop a calcific aortic stenosis (CAS), as well as a number of other ailments, as compared to their cohorts with normal tricuspid aortic valves (TAV). It is currently unknown whether the increase in risk of CAS is caused by the geometric differences between the tricuspid and bicuspid valves or whether the increase in risk is caused by the same underlying factors that produce the geometric difference. CAS progression is understood to be a multiscale process, mediated at the cell level. In this study, we employ multiscale finite-element simulations of the valves. We isolate the effect of one geometric factor, the number of cusps, in order to explore its effect on multiscale valve mechanics, particularly in relation to CAS. The BAV and TAV are modeled by a set of simulations describing the cell, tissue, and organ length scales. These simulations are linked across the length scales to create a coherent multiscale model. At each scale, the models are three-dimensional, dynamic, and incorporate accurate nonlinear constitutive models of the valve leaflet tissue. We compare results between the TAV and BAV at each length scale. At the cell-scale, our region of interest is the location where calcification develops, near the aortic-facing surface of the leaflet. Our simulations show the observed differences between the tricuspid and bicuspid valves at the organ scale: the bicuspid valve shows greater flexure in the solid phase and stronger jet formation in the fluid phase relative to the tricuspid. At the cell-scale, however, we show that the region of interest is shielded against strain by the wrinkling of the fibrosa. Thus, the cellular deformations are not significantly different between the TAV and BAV in the calcification-prone region. This result supports the assertion that the difference in calcification observed in the BAV versus TAV may be due primarily to factors other than the simple geometric difference between the two valves.     

Experimental measurement of dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. Although exact causes and mechanisms of AV calcification are unclear, previous studies suggest that mechanical forces play a role. Since calcium deposits occur almost exclusively on the aortic surfaces of AV leaflets, it has been hypothesized that adverse patterns of fluid shear stress on the aortic surface of AV leaflets promote calcification. The current study characterizes AV leaflet aortic surface fluid shear stresses using Laser Doppler velocimetry and an in vitro pulsatile flow loop. The valve model used was a native porcine valve mounted on a suturing ring and preserved using 0.15% glutaraldehyde solution. This valve model was inserted in a mounting chamber with sinus geometries, which is made of clear acrylic to provide optical access for measurements. To understand the effects of hemodynamics on fluid shear stress, shear stress was measured across a range of conditions: varying stroke volumes at the same heart rate and varying heart rates at the same stroke volume. Systolic shear stress magnitude was found to be much higher than diastolic shear stress magnitude due to the stronger flow in the sinuses during systole, reaching up to 20 dyn/cm² at mid-systole. Upon increasing stroke volume, fluid shear stresses increased due to stronger sinus fluid motion. Upon increasing heart rate, fluid shear stresses decreased due to reduced systolic duration that restricted the formation of strong sinus flow. Significant changes in the shear stress waveform were observed at 90 beats/min, most likely due to altered leaflet dynamics at this higher heart rate. Overall, this study represents the most well-resolved shear stress measurements to date across a range of conditions on the aortic side of the AV. The data presented can be used for further investigation to understand AV biological response to shear stresses.     

Deficient signaling via Alk2 (Acvr1) leads to bicuspid aortic valve development

Bicuspid aortic valve (BAV) is the most common congenital cardiac anomaly in humans. Despite recent advances, the molecular basis of BAV development is poorly understood. Previously it has been shown that mutations in the Notch1 gene lead to BAV and valve calcification both in human and mice, and mice deficient in Gata5 or its downstream target Nos3 have been shown to display BAVs. Here we show that tissue-specific deletion of the gene encoding Activin Receptor Type I (Alk2 or Acvr1) in the cushion mesenchyme results in formation of aortic valve defects including BAV. These defects are largely due to a failure of normal development of the embryonic aortic valve leaflet precursor cushions in the outflow tract resulting in either a fused right- and non-coronary leaflet, or the presence of only a very small, rudimentary non-coronary leaflet. The surviving adult mutant mice display aortic stenosis with high frequency and occasional aortic valve insufficiency. The thickened aortic valve leaflets in such animals do not show changes in Bmp signaling activity, while Map kinase pathways are activated. Although dysfunction correlated with some pro-osteogenic differences in gene expression, neither calcification nor inflammation were detected in aortic valves of Alk2 mutants with stenosis. We conclude that signaling via Alk2 is required for appropriate aortic valve development in utero, and that defects in this process lead to indirect secondary complications later in life.     

Defining the Role of Fluid Shear Stress in the Expression of Early Signaling Markers for Calcific Aortic Valve Disease

Calcific aortic valve disease (CAVD) is an active process presumably triggered by interplays between cardiovascular risk factors, molecular signaling networks and hemodynamic cues. While earlier studies demonstrated that alterations in fluid shear stress (FSS) on the fibrosa could trigger inflammation, the mechanisms of CAVD pathogenesis secondary to side-specific FSS abnormalities are poorly understood. This knowledge could be critical to the elucidation of key CAVD risk factors such as congenital valve defects, aging and hypertension, which are known to generate FSS disturbances. The objective of this study was to characterize ex vivo the contribution of isolated and combined abnormalities in FSS magnitude and frequency to early valvular pathogenesis. The ventricularis and fibrosa of porcine aortic valve leaflets were exposed simultaneously to different combinations of sub-physiologic/physiologic/supra-physiologic levels of FSS magnitude and frequency for 24, 48 and 72 hours in a double cone-and-plate device. Endothelial activation and paracrine signaling were investigated by measuring cell-adhesion molecule (ICAM-1, VCAM-1) and cytokine (BMP-4, TGF-β1) expressions, respectively. Extracellular matrix (ECM) degradation was characterized by measuring the expression and activity of the proteases MMP-2, MMP-9, cathepsin L and cathepsin S. The effect of the FSS treatment yielding the most significant pathological response was examined over a 72-hour period to characterize the time-dependence of FSS mechano-transduction. While cytokine expression was stimulated under elevated FSS magnitude at normal frequency, ECM degradation was stimulated under both elevated FSS magnitude at normal frequency and physiologic FSS magnitude at abnormal frequency. In contrast, combined FSS magnitude and frequency abnormalities essentially maintained valvular homeostasis. The pathological response under supra-physiologic FSS magnitude peaked at 48 hours but was then maintained until the 72-hour time point. This study confirms the sensitivity of valve leaflets to both FSS magnitude and frequency and suggests the ability of supra-physiologic FSS levels or abnormal FSS frequencies to initiate CAVD mechanisms.     

Experimental technique of measuring dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact cause and mechanism of the progression of AV calcification is unknown, although mechanical forces have been known to play a role. It is thus important to characterize the mechanical environment of the AV. In the current study, we establish a methodology of measuring shear stresses experienced by the aortic surface of the AV leaflets using an in vitro valve model and adapting the laser Doppler velocimetry (LDV) technique. The valve model was constructed from a fresh porcine aortic valve, which was trimmed and sutured onto a plastic stented ring, and inserted into an idealized three-lobed sinus acrylic chamber. Valve leaflet location was measured by obtaining the location of highest back-scattered LDV laser light intensity. The technique of performing LDV measurements near to biological surfaces as well as the leaflet locating technique was first validated in two phantom flow systems: (1) steady flow within a straight tube with AV leaflet adhered to the wall, and (2) steady flow within the actual valve model. Dynamic shear stresses were then obtained by applying the techniques on the valve model in a physiologic pulsatile flow loop. Results show that aortic surface shear stresses are low during early systole (<5 dyn/cm2) but elevated to its peak during mid to late systole at about 18-20 dyn/cm2. Low magnitude shear stress (<5 dyn/cm2) was observed during early diastole and dissipated to zero over the diastolic duration. Systolic shear stress was observed to elevate only with the formation of sinus vortex flow. The presented technique can also be used on other in vitro valve models such as congenitally geometrically malformed valves, or to investigate effects of hemodynamics on valve shear stress. Shear stress data can be used for further experiments investigating effects of fluid shear stress on valve biology, for conditioning tissue engineered AV, and to validate numerical simulations.     

Simulations of morphotype-dependent hemodynamics in non-dilated bicuspid aortic valve aortas

Bicuspid aortic valves (BAVs) generate flow abnormalities that may promote aortopathy. While positive helix fraction (PHF) index, flow angle (θ), flow displacement (d) and wall shear stress (WSS) exhibit abnormalities in dilated BAV aortas, it is unclear whether those anomalies stem from the abnormal valve anatomy or the dilated aorta. Therefore, the objective of this study was to quantify the early impact of different BAV morphotypes on aorta hemodynamics prior to dilation. Fluid-structure interaction models were designed to quantify standard peak-systolic flow metrics and temporal WSS characteristics in a realistic non-dilated aorta connected to functional tricuspid aortic valve (TAV) and type-I BAVs. While BAVs generated increased helicity (PHF>0.68) in the middle ascending aorta (AA), larger systolic flow skewness (θ>11.2°) and displacement (d>6.8 mm) relative to the TAV (PHF=0.51; θ<5.5°; d<3.3 mm), no distinct pattern was observed between morphotypes. In contrast, WSS magnitude and directionality abnormalities were BAV morphotype- and site-dependent. Type-I BAVs subjected the AA convexity to peak-systolic WSS overloads (up to 1014% difference vs. TAV). While all BAVs increased WSS unidirectionality on the proximal AA relative to the TAV, the most significant abnormality was achieved by the BAV with left-right-coronary cusp fusion on the wall convexity (up to 0.26 decrease in oscillatory shear index vs. TAV). The results indicate the existence of strong hemodynamic abnormalities in non-dilated type-I BAV AAs, their colocalization with sites vulnerable to dilation and the superior specificity of WSS metrics over global hemodynamic metrics to the valve anatomy.     

Helical flows and asymmetry of blood jet in dilated ascending aorta with normally functioning bicuspid valve

Bicuspid aortic valve (BAV) is associated with aortic dilatation and aneurysm. Several studies evidenced an eccentric systolic flow in ascending aorta associated with increased wall shear stresses (WSS) and the occurrence of an helical systolic flow. This study seeks to elucidate the connections between jet asymmetry and helical flow in patients with normally functioning BAV and dilated ascending aorta. We performed a computational parametric study by varying, for a patient-specific geometry, the valve area and the flow rate entering the aorta and drawing also a tricuspid valve (TAV). We considered also phase-contrast magnetic resonance imaging of four BAV and TAV patients. Measurement of normalized flow asymmetry index, systolic WSS and of a new index (positive helix fraction, PHF) quantifying the presence of a single a single helical flow were performed. In our computation, BAV cases featured higher values of all indices with respect to TAV in both numerical and imaged-based results. Moreover, all indices increased with decreasing valve area and/or with increasing flow rate. This allowed to separate the BAV and TAV cases with respect to the jet asymmetry, WSS localization and helical flow. Interestingly, these results were obtained without modeling the leaflets.     

Unraveling divergent gene expression profiles in bicuspid and tricuspid aortic valve patients with thoracic aortic dilatation: the ASAP study

Thoracic aortic aneurysm (TAA) is a common complication in patients with a bicuspid aortic valve (BAV), the most frequent congenital heart disorder. For unknown reasons TAA occurs at a younger age, with a higher frequency in BAV patients than in patients with a tricuspid aortic valve (TAV), resulting in an increased risk for aortic dissection and rupture. To investigate the increased TAA incidence in BAV patients, we obtained tissue biopsy samples from nondilated and dilated aortas of 131 BAV and TAV patients. Global gene expression profiles were analyzed from controls and from aortic intima-media and adventitia of patients (in total 345 samples). Of the genes found to be differentially expressed with dilation, only a few (<4%) were differentially expressed in both BAV and TAV patients. With the use of gene set enrichment analysis, the cell adhesion and extracellular region gene ontology sets were identified as common features of TAA in both BAV and TAV patients. Immune response genes were observed to be particularly overexpressed in the aortic media of dilated TAV samples. The divergent gene expression profiles indicate that there are fundamental differences in TAA etiology in BAV and TAV patients. Immune response activation solely in the aortic media of TAV patients suggests that inflammation is involved in TAA formation in TAV but not in BAV patients. Conversely, genes were identified that were only differentially expressed with dilation in BAV patients. The result has bearing on future clinical studies in which separate analysis of BAV and TAV patients is recommended.     

Experimental measurement of dynamic fluid shear stress on the ventricular surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact causes and mechanisms of AV calcification are unclear, although previous studies suggest that mechanical forces play a role. It has been clinically demonstrated that calcification preferentially occurs on the aortic surface of the AV. This is hypothesized to be due to differences in the mechanical environments on the two sides of the valve. It is thus necessary to characterize fluid shear forces acting on both sides of the leaflet to test this hypothesis. The current study is one of two studies characterizing dynamic shear stress on both sides of the AV leaflets. In the current study, shear stresses on the ventricular surface of the AV leaflets were measured experimentally on two prosthetic AV models with transparent leaflets in an in vitro pulsatile flow loop using two-component Laser Doppler Velocimetry (LDV). Experimental measurements were utilized to validate a theoretical model of AV ventricular surface shear stress based on the Womersley profile in a straight tube, with corrections for the opening angle of the valve leaflets. This theoretical model was applied to in vivo data based on MRI-derived volumetric flow rates and valve dimension obtained from the literature. Experimental results showed that ventricular surface shear stress was dominated by the streamwise component. The systolic shear stress waveform resembled a half-sinusoid during systole and peaks at 64–71 dyn/cm², and reversed in direction at the end of systole for 15–25?ms, and reached a significant negative magnitude of 40–51 dyn/cm². Shear stresses from the theoretical model applied to in vivo data showed that shear stresses peaked at 77–92 dyn/cm² and reversed in direction for substantial period of time (108–110?ms) during late systole with peak negative shear stress of 35–38 dyn/cm².     

Normal and abnormal development of the aortic wall and valve: correlation with clinical entities

Dilation of the wall of the thoracic aorta can be found in patients with a tricuspid (TAV) as well as a bicuspid aortic valve (BAV) with and without a syndromic component. BAV is the most common congenital cardiovascular malformation, with a population prevalence of 0.5–2 %. The clinical course is often characterised by aneurysm formation and in some cases dissection. The non-dilated aortic wall is less well differentiated in all BAV as compared with TAV, thereby conferring inherent developmental susceptibility. Furthermore, a turbulent flow, caused by the inappropriate opening of the bicuspid valve, could accelerate the degenerative process in the aortic wall. However, not all patients with bicuspidy develop clinical complications during their life. We postulate that the increased vulnerability for aortic complications in a subset of patients with BAV is caused by a defect in the early development of the aorta and aortic valve. This review discusses histological and molecular genetic aspects of the normal and abnormal development of the aortic wall and semilunar valves. Aortopathy associated with BAV could be the result of a shared developmental defect during embryogenesis.     

Supra-annular Valve-in-Valve implantation reduces blood stasis on the transcatheter aortic valve leaflets

Leaflet thrombosis following transcatheter aortic valve replacement (TAVR) and Valve-in-Valve (ViV) procedures has been increasingly recognized. This study aimed to investigate the effect of positioning of the transcatheter aortic valve (TAV) in ViV setting on the flow dynamics aspect of post-ViV thrombosis by quantifying the blood stasis in the intra-annular and supra-annular settings. To that end, two idealized computational models, representing ViV intra-annular and supra-annular positioning of a TAV were developed in a patient-specific geometry. Three-dimensional flow fields were then obtained via fluid-solid interaction modeling to study the difference in blood residence time (BRT) on the TAV leaflets in the two settings. At the end of diastole, a strip of high BRT ($⩾1.2s$) region was observed on the TAV leaflets in the ViV intra-annular positioning at the fixed boundary where the leaflets are attached to the frame. Such a high BRT region was absent on the TAV leaflets in the supra-annular positioning. The maximum value of BRT on the surface of non-, right, and left coronary leaflets of the TAV in the supra-annular positioning were 53%, 11%, and 27% smaller compared to the intra-annular positioning, respectively. It was concluded that the geometric confinement of TAV by the leaflets of the failed bioprosthetic valve in ViV intra-annular positioning increases the BRT on the leaflets and may act as a permissive factor in valvular thrombosis. The absence of such a geometric confinement in the ViV supra-annular positioning leads to smaller BRT and subsequently less likelihood of leaflet thrombosis.     

Bicuspid aortic valve hemodynamics induces abnormal medial remodeling in the convexity of porcine ascending aortas

The type-I bicuspid aortic valve (BAV), which differs from the normal tricuspid aortic valve (TAV) most commonly by left-right coronary cusp fusion, is frequently associated with secondary aortopathies. While BAV aortic dilation has been linked to a genetic predisposition, hemodynamics has emerged as a potential alternate etiology. However, the link between BAV hemodynamics and aortic medial degeneration has not been established. The objective of this study was to compare the regional wall shear stresses (WSS) in a TAV and BAV ascending aorta (AA) and to isolate ex vivo their respective impact on aortic wall remodeling. The WSS environments generated in the convex region of a TAV and BAV AA were predicted through fluid–structure interaction (FSI) simulations in an aorta model subjected to both valvular flows. Remodeling of porcine aortic tissue exposed to TAV and BAV AA WSS for 48 h in a cone-and-plate bioreactor was investigated via immunostaining, immunoblotting and zymography. FSI simulations revealed the existence of larger and more unidirectional WSS in the BAV than in the TAV AA convexity. Exposure of normal aortic tissue to BAV AA WSS resulted in increased MMP-2 and MMP-9 expressions and MMP-2 activity but similar fibrillin-1 content and microfibril organization relative to the TAV AA WSS treatment. This study confirms the sensitivity of aortic tissue to WSS abnormalities and demonstrates the susceptibility of BAV hemodynamic stresses to focally mediate aortic medial degradation. The results provide compelling support to the important role of hemodynamics in BAV secondary aortopathy.     

Patient-specific modeling of biomechanical interaction in transcatheter aortic valve deployment

The objective of this study was to develop a patient-specific computational model to quantify the biomechanical interaction between the transcatheter aortic valve (TAV) stent and the stenotic aortic valve during TAV intervention. Finite element models of a patient-specific stenotic aortic valve were reconstructed from multi-slice computed tomography (MSCT) scans, and TAV stent deployment into the aortic root was simulated. Three initial aortic root geometries of this patient were analyzed: (a) aortic root geometry directly reconstructed from MSCT scans, (b) aortic root geometry at the rapid right ventricle pacing phase, and (c) aortic root geometry with surrounding myocardial tissue. The simulation results demonstrated that stress, strain, and contact forces of the aortic root model directly reconstructed from MSCT scans were significantly lower than those of the model at the rapid ventricular pacing phase. Moreover, the presence of surrounding myocardium slightly increased the mechanical responses. Peak stresses and strains were observed around the calcified regions in the leaflets, suggesting the calcified leaflets helped secure the stent in position. In addition, these elevated stresses induced during TAV stent deployment indicated a possibility of tissue tearing and breakdown of calcium deposits, which might lead to an increased risk of stroke. The potential of paravalvular leak and occlusion of coronary ostia can be evaluated from simulated post-deployment aortic root geometries. The developed computational models could be a valuable tool for pre-operative planning of TAV intervention and facilitate next generation TAV device design.     

Interlayer micromechanics of the aortic heart valve leaflet

While the mechanical behaviors of the fibrosa and ventricularis layers of the aortic valve (AV) leaflet are understood, little information exists on their mechanical interactions mediated by the GAG-rich central spongiosa layer. Parametric simulations of the interlayer interactions of the AV leaflets in flexure utilized a tri-layered finite element (FE) model of circumferentially oriented tissue sections to investigate inter-layer sliding hypothesized to occur. Simulation results indicated that the leaflet tissue functions as a tightly bonded structure when the spongiosa effective modulus was at least 25 % that of the fibrosa and ventricularis layers. Novel studies that directly measured transmural strain in flexure of AV leaflet tissue specimens validated these findings. Interestingly, a smooth transmural strain distribution indicated that the layers of the leaflet indeed act as a bonded unit, consistent with our previous observations (Stella and Sacks in J Biomech Eng 129:757–766, 2007) of a large number of transverse collagen fibers interconnecting the fibrosa and ventricularis layers. Additionally, when the tri-layered FE model was refined to match the transmural deformations, a layer-specific bimodular material model (resulting in four total moduli) accurately matched the transmural strain and moment-curvature relations simultaneously. Collectively, these results provide evidence, contrary to previous assumptions, that the valve layers function as a bonded structure in the low-strain flexure deformation mode. Most likely, this results directly from the transverse collagen fibers that bind the layers together to disable physical sliding and maintain layer residual stresses. Further, the spongiosa may function as a general dampening layer while the AV leaflets deforms as a homogenous structure despite its heterogeneous architecture.     

Simulated elliptical bioprosthetic valve deformation: Implications for asymmetric transcatheter valve deployment

The asymmetric, elliptical shape of a transcatheter aortic valve (TAV), after implantation into a calcified aortic root, has been clinically observed. However, the impact of elliptical TAV configuration on TAV leaflet stress and strain distribution and valve regurgitation is largely unknown. In this study, we developed computational models of elliptical TAVs based on a thin pericardial bioprosthetic valve model recently developed. Finite element and computational fluid dynamics simulations were performed to investigate TAV leaflet structural deformation and central backflow leakage, and compared with those of a nominal symmetric TAV. From the results, we found that for a distorted TAV with an elliptical eccentricity of 0.68, the peak stress increased significantly by 143% compared with the nominal circular TAV. When the eccentricity of an elliptical TAV was larger than 0.5, a central backflow leakage was likely to occur. Also, deployment of a TAV with a major calcified region perpendicular to leaflet coaptation line was likely to cause a larger valve leakage. In conclusion, the computational models of elliptical TAVs developed in this study could improve our understanding of the biomechanics involved in a TAV with an elliptical configuration and facilitate optimal design of next-generation TAV devices.     

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\(0.36-0.87\)) with a strong correlation with the organ length-scale radial strain (\(R^{2}= 0.95\)); conversely, circumferential cell strains spanned a very narrow range (\(0.75-0.88\)) showing no correlation with the circumferential strain at the organ level (\(R^{2}= 0.02\)). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

     

Micromechanics of the fibrosa and the ventricularis in aortic valve leaflets.

The elastic response of aortic valve cusps is a summation of its fibrous components. To investigate the micromechanical function of valve leaflet constituents, we separated the fibrosa and the ventricularis from fresh and glutaraldehyde-fixed leaflets and tested them individually. The ventricularis was stiffer circumferentially than radially (7.41 kPa vs 3.68 kPa, p less than 0.00001) and was more extensible radially (62.7% vs 21.8% strain to high modulus phase, p less than 0.00001). The fibrosa was also stiffer circumferentially than radially (13.02 kPa vs 4.65 kPa, p less than 0.0008), but had uniform extensibility. Glutaraldehyde fixation did not affect the circumferential elastic modulus of the fibrosa, but reduced its radial modulus from 4.65 kPa to 2.32 kPa (p less than 0.0078). The elastic modulus of the ventricularis remained unchanged. Fixation also reduced the extensibility of the ventricularis circumferentially (from 21.8% to 15.2% strain, p less than 0.018), but not radially, and increased the radial extensibility of the fibrosa from 27.7% to 46.1% (p less than 0.0048). These data show that while the ventricularis contains a large amount of elastin, the amount of radially oriented collagen is similar to that of the fibrosa. The fibrosa, by itself, has the same extensibility in both directions (about 23% strain), but can extend much more radially when connected to the rest of the leaflet because it is attached to the ventricularis in a highly folded configuration. The two layers therefore complement each other during aortic valve function, and become detrimentally altered by fixation in glutaraldehyde.