Heart Valve Biomechanics: The Frontiers of Modeling Modalities and the Expansive Capabilities of Ex Vivo Heart Simulation

The field of heart valve biomechanics is a rapidly expanding, highly clinically relevant area of research. While most valvular pathologies are rooted in biomechanical changes, the technologies for studying these pathologies and identifying treatments have largely been limited. Nonetheless, significant advancements are underway to better understand the biomechanics of heart valves, pathologies, and interventional therapeutics, and these advancements have largely been driven by crucial in silico, ex vivo, and in vivo modeling technologies. These modalities represent cutting-edge abilities for generating novel insights regarding native, disease, and repair physiologies, and each has unique advantages and limitations for advancing study in this field. In particular, novel ex vivo modeling technologies represent an especially promising class of translatable research that leverages the advantages from both in silico and in vivo modeling to provide deep quantitative and qualitative insights on valvular biomechanics. The frontiers of this work are being discovered by innovative research groups that have used creative, interdisciplinary approaches toward recapitulating in vivo physiology, changing the landscape of clinical understanding and practice for cardiovascular surgery and medicine.

Influence of aneurysmal aortic root geometry on mechanical stress to the aortic valve leaflet

AIMS

While mechanical stress caused by blood flow, e.g. wall shear stress (WSS), and related parameters, e.g. oscillatory shear index (OSI), are increasingly being recognized as key moderators of various cardiovascular diseases, studies on valves have been limited because of a lack of appropriate imaging modalities. We investigated the influence of aortic root geometry on WSS and OSI on the aortic valve (AV) leaflet.

METHODS AND RESULTS

We applied our novel approach of intraoperative epi-aortic echocardiogram to measure the haemodynamic parameters of WSS and OSI on the AV leaflet. Thirty-six patients were included, which included those who underwent valve-sparing aortic root replacement (VSARR) with no significant aortic regurgitation (n = 17) and coronary artery bypass graft (CABG) with normal AV (n = 19). At baseline, those who underwent VSARR had a higher systolic WSS (0.52 ± 0.12 vs. 0.32 ± 0.08 Pa, respectively, P < 0.001) and a higher OSI (0.37 ± 0.06 vs. 0.29 ± 0.04, respectively, P < 0.001) on the aortic side of the AV leaflet than those who underwent CABG. Multivariate regression analysis revealed that the size of the sinus of Valsalva had a significant association with WSS and OSI. Following VSARR, WSS and OSI values decreased significantly compared with the baseline values (WSS: 0.29 ± 0.12 Pa, P < 0.001; OSI: 0.26 ± 0.09, P < 0.001), and became comparable to the values in those who underwent CABG (WSS, P = 0.42; OSI, P = 0.15).

CONCLUSIONS

Mechanical stress on the AV gets altered in correlation with the size of the aortic root. An aneurysmal aortic root may expose the leaflet to abnormal fluid dynamics. The VSARR procedure appeared to reduce these abnormalities.

Thinner biological tissues induce leaflet flutter in aortic heart valve replacements

Valvular heart disease has recently become an increasing public health concern due to the high prevalence of valve degeneration in aging populations. For patients with severely impacted aortic valves that require replacement, catheter-based bioprosthetic valve deployment offers a minimally invasive treatment option that eliminates many of the risks associated with surgical valve replacement. Although recent percutaneous device advancements have incorporated thinner, more flexible biological tissues to streamline safer deployment through catheters, the impact of such tissues in the complex, mechanically demanding, and highly dynamic valvular system remains poorly understood. The present work utilized a validated computational fluid-structure interaction approach to isolate the behavior of thinner, more compliant aortic valve tissues in a physiologically realistic system. This computational study identified and quantified significant leaflet flutter induced by the use of thinner tissues that initiated blood flow disturbances and oscillatory leaflet strains. The aortic flow and valvular dynamics associated with these thinner valvular tissues have not been previously identified and provide essential information that can significantly advance fundamental knowledge about the cardiac system and support future medical device innovation. Considering the risks associated with such observed flutter phenomena, including blood damage and accelerated leaflet deterioration, this study demonstrates the potentially serious impact of introducing thinner, more flexible tissues into the cardiac system.

A novel in vivo assessment of fluid dynamics on aortic valve leaflet using epi-aortic echocardiogram

BACKGROUND

Mechanical stress caused by blood flow, such as wall shear stress (WSS) and its related parameters, is key moderator of endothelial degeneration. However, an in vivo method to measure WSS on heart valves has not been developed.

METHODS

We developed a novel approach, based on vector flow mapping using intraoperative epi-aortic echocardiogram, to measure WSS and oscillatory shear index (OSI) on the aortic valve. We prospectively enrolled 15 patients with normal valves, who underwent coronary artery bypass graft.

RESULTS

Systolic WSS on the ventricularis (2.40 ± 0.44 Pa [1.45-3.00 Pa]) was higher than systolic WSS on the fibrosa (0.33 ± 0.08 Pa [0.14-0.47 Pa], P < .001) and diastolic WSS on the ventricularis (0.18 ± 0.07 Pa [0.04-0.28 Pa], P < .001). Oscillatory shear index on the fibrosa was higher than on the ventricularis (0.29 ± 0.04 [0.24-0.36] vs 0.05 ± 0.03 [0.01-0.12], P < .001). A pilot study involving two patients with severe aortic regurgitation showed significantly different values in fluid dynamics.

CONCLUSION

Vector flow mapping method using intraoperative epi-aortic echocardiogram is an effective way of measuring WSS and OSI on normal aortic leaflet in vivo, allowing for better understanding of the pathophysiology of aortic valve diseases.

Immersogeometric fluid-structure interaction modeling and simulation of transcatheter aortic valve replacement

The transcatheter aortic valve replacement (TAVR) has emerged as a minimally invasive alternative to surgical treatments of valvular heart disease. TAVR offers many advantages, however, the safe anchoring of the transcatheter heart valve (THV) in the patients anatomy is key to a successful procedure. In this paper, we develop and apply a novel immersogeometric fluid-structure interaction (FSI) framework for the modeling and simulation of the TAVR procedure to study the anchoring ability of the THV. To account for physiological realism, methods are proposed to model and couple the main components of the system, including the arterial wall, blood flow, valve leaflets, skirt, and frame. The THV is first crimped and deployed into an idealized ascending aorta. During the FSI simulation, the radial outward force and friction force between the aortic wall and the THV frame are examined over the entire cardiac cycle. The ratio between these two forces is computed and compared with the experimentally estimated coefficient of friction to study the likelihood of valve migration.

An anisotropic constitutive model for immersogeometric fluid-structure interaction analysis of bioprosthetic heart valves

This paper considers an anisotropic hyperelastic soft tissue model, originally proposed for native valve tissue and referred to herein as the Lee-Sacks model, in an isogeometric thin shell analysis framework that can be readily combined with immersogeometric fluid-structure interaction (FSI) analysis for high-fidelity simulations of bioprosthetic heart valves (BHVs) interacting with blood flow. We find that the Lee-Sacks model is well-suited to reproduce the anisotropic stress-strain behavior of the cross-linked bovine pericardial tissues that are commonly used in BHVs. An automated procedure for parameter selection leads to an instance of the Lee-Sacks model that matches biaxial stress-strain data from the literature more closely, over a wider range of strains, than other soft tissue models. The relative simplicity of the Lee-Sacks model is attractive for computationally-demanding applications such as FSI analysis and we use the model to demonstrate how the presence and direction of material anisotropy affect the FSI dynamics of BHV leaflets.

A framework for designing patient-specific bioprosthetic heart valves using immersogeometric fluid-structure interaction analysis

Numerous studies have suggested that medical image derived computational mechanics models could be developed to reduce mortality and morbidity due to cardiovascular diseases by allowing for patient-specific surgical planning and customized medical device design. In this work, we present a novel framework for designing prosthetic heart valves using a parametric design platform and immersogeometric fluid-structure interaction (FSI) analysis. We parameterize the leaflet geometry using several key design parameters. This allows for generating various perturbations of the leaflet design for the patient-specific aortic root reconstructed from the medical image data. Each design is analyzed using our hybrid arbitrary Lagrangian-Eulerian/immersogeometric FSI methodology, which allows us to efficiently simulate the coupling of the deforming aortic root, the parametrically designed prosthetic valves, and the surrounding blood flow under physiological conditions. A parametric study is performed to investigate the influence of the geometry on heart valve performance, indicated by the effective orifice area and the coaptation area. Finally, the FSI simulation result of a design that balances effective orifice area and coaptation area reasonably well is compared with patient-specific phase contrast magnetic resonance imaging data to demonstrate the qualitative similarity of the flow patterns in the ascending aorta.

Computational Modeling for Fluid–Structure Interaction of Bioprosthetic Heart Valve with Different Suture Density: Comparison with Dynamic Structure Simulation

In this paper, a parametric geometry model based on elliptic and conic surfaces was developed for bioprosthetic heart valve (BHV) simulation. The valve material was modeled by a hyperelastic nonlinear anisotropic solid model. Different suture densities could be substituted by various bonded points between artery vessel and the leaflets as boundary conditions in the computational modeling. Besides these two assumptions that dynamic structure (DS) and fluid–structure interaction (FSI) both shared, the latter need incompressible viscous Newton fluid model to depict bloodstream passing through the BHV. Immersed boundary (IB) method was introduced to solve the FSI simulation. In addition, the DS analysis applied transvalvular pressure on the valve while FSI had left ventricular pressure on fluid inlet as initials. There was inconsistency between the moments of the maximum deformation and the maximum loading in both simulations. Although a similar trend of deformation lagging the load was viewed, the extent of delay in FSI was much smaller compared with that in DS simulation. The deformed profiles in cross-sectional views were shown in one picture to illustrate different dynamic responses caused by distinct assumptions. Percent of open area at the moments when the maximum deformation occurred was defined to show which calculation achieved better approximation for precise hemodynamics. Fixed point was given as boundaries between BHV and artery in the modeling part. Calculations showed that the more the fixed points in this bonded contact, the lower the principal stress was. The maximum shear stress showed a different trend. It had a different trend. Stress concentration in the conjunction area made it high-risk to be teared. Different suture densities had significant impaction in FSI simulations. With that analysis our work achieved a more comprehensive simulation to describe true hemodynamics of a BHV implanted in artery. The artery vessel had particular dynamic response under such assumptions, gradient existed in the maximum principal stress distribution diagram, from inner wall through which blood passing to the outer wall. Results showed a large suture density was suggested in BHV implantation.

The effects of positioning of transcatheter aortic valves on fluid dynamics of the aortic root

Transcatheter aortic valve implantation is a novel treatment for severe aortic valve stenosis. Due to the recent use of this technology and the procedural variability, there is very little data that quantify the hemodynamic consequences of variations in valve placement. Changes in aortic wall stresses and fluid retention in the sinuses of Valsalva can have a significant effect on the clinical response a patient has to the procedure. By comprehensively characterizing complex flow in the sinuses of Valsalva using digital particle image velocimetry and an advanced heart-flow simulator, various positions of a deployed transcatheter valve with respect to a bioprosthetic aortic valve (valve-in-valve) were tested in vitro. Displacements of the transcatheter valve were axial and directed below the simulated native valve annulus. It was determined that for both blood residence time and aortic Reynolds stresses, it is optimal to have the annulus of the transcatheter valve deployed as close to the aortic valve annulus as possible.

Dynamic and fluid-structure interaction simulations of bioprosthetic heart valves using parametric design with T-splines and Fung-type material models

This paper builds on a recently developed immersogeometric fluid-structure interaction (FSI) methodology for bioprosthetic heart valve (BHV) modeling and simulation. It enhances the proposed framework in the areas of geometry design and constitutive modeling. With these enhancements, BHV FSI simulations may be performed with greater levels of automation, robustness and physical realism. In addition, the paper presents a comparison between FSI analysis and standalone structural dynamics simulation driven by prescribed transvalvular pressure, the latter being a more common modeling choice for this class of problems. The FSI computation achieved better physiological realism in predicting the valve leaflet deformation than its standalone structural dynamics counterpart.

Coronary Flow Impacts Aortic Leaflet Mechanics and Aortic Sinus Hemodynamics

Mechanical stresses on aortic valve leaflets are well-known mediators for initiating processes leading to calcific aortic valve disease. Given that non-coronary leaflets calcify first, it may be hypothesized that coronary flow originating from the ostia significantly influences aortic leaflet mechanics and sinus hemodynamics. High resolution time-resolved particle image velocimetry (PIV) measurements were conducted to map the spatiotemporal characteristics of aortic sinus blood flow and leaflet motion with and without physiological coronary flow in a well-controlled in vitro setup. The in vitro setup consists of a porcine aortic valve mounted in a physiological aorta sinus chamber with dynamically controlled coronary resistance to emulate physiological coronary flow. Results were analyzed using qualitative streak plots illustrating the spatiotemporal complexity of blood flow patterns, and quantitative velocity vector and shear stress contour plots to show differences in the mechanical environments between the coronary and non-coronary sinuses. It is shown that the presence of coronary flow pulls the classical sinus vorticity deeper into the sinus and increases flow velocity near the leaflet base. This creates a beneficial increase in shear stress and washout near the leaflet that is not seen in the non-coronary sinus. Further, leaflet opens approximately 10% farther into the sinus with coronary flow case indicating superior valve opening area. The presence of coronary flow significantly improves leaflet mechanics and sinus hemodynamics in a manner that would reduce low wall shear stress conditions while improving washout at the base of the leaflet.

Cardiac valve cells and their microenvironment--insights from in vitro studies

During every heartbeat, cardiac valves open and close coordinately to control the unidirectional flow of blood. In this dynamically challenging environment, resident valve cells actively maintain homeostasis, but the signalling between cells and their microenvironment is complex. When homeostasis is disrupted and the valve opening obstructed, haemodynamic profiles can be altered and lead to impaired cardiac function. Currently, late stages of cardiac valve diseases are treated surgically, because no drug therapies exist to reverse or halt disease progression. Consequently, investigators have sought to understand the molecular and cellular mechanisms of valvular diseases using in vitro cell culture systems and biomaterial scaffolds that can mimic the extracellular microenvironment. In this Review, we describe how signals in the extracellular matrix regulate valve cell function. We propose that the cellular context is a critical factor when studying the molecular basis of valvular diseases in vitro, and one should consider how the surrounding matrix might influence cell signalling and functional outcomes in the valve. Investigators need to build a systems-level understanding of the complex signalling network involved in valve regulation, to facilitate drug target identification and promote in situ or ex vivo heart valve regeneration.

Fluid-structure interaction analysis of bioprosthetic heart valves: Significance of arterial wall deformation

We propose a framework that combines variational immersed-boundary and arbitrary Lagrangian-Eulerian (ALE) methods for fluid-structure interaction (FSI) simulation of a bioprosthetic heart valve implanted in an artery that is allowed to deform in the model. We find that the variational immersed-boundary method for FSI remains robust and effective for heart valve analysis when the background fluid mesh undergoes deformations corresponding to the expansion and contraction of the elastic artery. Furthermore, the computations presented in this work show that the arterial wall deformation contributes significantly to the realism of the simulation results, leading to flow rates and valve motions that more closely resemble those observed in practice.

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.

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(2), and reversed in direction at the end of systole for 15-25 ms, and reached a significant negative magnitude of 40-51 dyn/cm(2). Shear stresses from the theoretical model applied to in vivo data showed that shear stresses peaked at 77-92 dyn/cm(2) 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(2).