Dynamic simulation pericardial bioprosthetic heart valve function

While providing nearly trouble-free function for 10-12 years, current bioprosthetic heart valves (BHV) continue to suffer from limited long-term durability. This is usually a result of leaflet calcification and/or structural degeneration, which may be related to regions of stress concentration associated with complex leaflet deformations. In the current work, a dynamic three-dimensional finite element analysis of a pericardial BHV was performed with a recently developed FE implementation of the generalized nonlinear anisotropic Fung-type elastic constitutive model for pericardial BHV tissues (W. Sun and M.S. Sacks, 2005, [Biomech. Model. Mechanobiol., 4(2-3), pp. 190-199]). The pericardial BHV was subjected to time-varying physiological pressure loading to compute the deformation and stress distribution during the opening phase of the valve function. A dynamic sequence of the displacements revealed that the free edge of the leaflet reached the fully open position earlier and the belly region followed. Asymmetry was observed in the resulting displacement and stress distribution due to the fiber direction and the anisotropic characteristics of the Fung-type elastic constitutive material model. The computed stress distribution indicated relatively high magnitudes near the free edge of the leaflet with local bending deformation and subsequently at the leaflet attachment boundary. The maximum computed von Mises stress during the opening phase was 33.8 kPa. The dynamic analysis indicated that the free edge regions of the leaflets were subjected to significant flexural deformation that may potentially lead to structural degeneration after millions of cycles of valve function. The regions subjected to time varying flexural deformation and high stresses of the present study also correspond to regions of tissue valve calcification and structural failure reported from explanted valves. In addition, the present simulation also demonstrated the importance of including the bending component together with the in-plane material behavior of the leaflets towards physiologically realistic deformation of the leaflets. Dynamic simulations with experimentally determined leaflet material specification can be potentially used to modify the valve towards an optimal design to minimize regions of stress concentration and structural failure.     

Simulated bioprosthetic heart valve deformation under quasi-static loading

For more than 40 years, the replacement of diseased natural heart valves with prosthetic devices has dramatically extended the quality and length of the lives of millions of patients worldwide. However, bioprosthetic heart valves (BHV) continue to fail due to structural failure resulting from poor tissue durability and faulty design. Clearly, an in-depth understanding of the biomechanical behavior of BHV at both the tissue and functional prosthesis levels is essential to improving BHV design and to reduce rates of failure. In this study, we simulated quasi-static BHV leaflet deformation under 40, 80, and 120 mm Hg quasi-static transvalvular pressures. A Fung-elastic material model was used that incorporated material parameters and axes derived from actual leaflet biaxial tests and measured leaflet collagen fiber structure. Rigorous experimental validation of predicted leaflet strain field was used to validate the model results. An overall maximum discrepancy of 2.36% strain between the finite element (FE) results and experiment measurements was obtained, indicating good agreement between computed and measured major principal strains. Parametric studies utilizing the material parameter set from one leaflet for all three leaflets resulted in substantial variations in leaflet stress and strain distributions. This result suggests that utilization of actual leaflet material properties is essential for accurate BHV FE simulations. The present study also underscores the need for rigorous experimentation and accurate constitutive models in simulating BHV function and design.     

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.     

Mapping the spatial variation of mitral valve elastic properties using air-pulse optical coherence elastography

The mitral valve is a highly heterogeneous tissue composed of two leaflets, anterior and posterior, whose unique composition and regional differences in material properties are essential to overall valve function. While mitral valve mechanics have been studied for many decades, traditional testing methods limit the spatial resolution of measurements and can be destructive. Optical coherence elastography (OCE) is an emerging method for measuring viscoelastic properties of tissues in a noninvasive, nondestructive manner. In this study, we employed air-pulse OCE to measure the spatial variation in mitral valve elastic properties with micro-scale resolution at 1 mm increments along the radial length of the leaflets. We analyzed differences between the leaflets, as well as between regions of the valve. We found that the anterior leaflet has a higher elastic wave velocity, which is reported as a surrogate for stiffness, than the posterior leaflet, most notably at the annular edge of the sample. In addition, we found a spatial elastic gradient in the anterior leaflet, where the annular edge was found to have a greater elastic wave velocity than the free edge. This gradient was less pronounced in the posterior leaflet. These patterns were confirmed using established uniaxial tensile testing methods. Overall, the anterior leaflet was stiffer and had greater heterogeneity in its mechanical properties than the posterior leaflet. This study measures differences between the two mitral leaflets with greater resolution than previously feasible and demonstrates a method that may be suitable for assessing valve mechanics following repair or during the engineering of synthetic valve replacements.     

Bioprosthetic heart valve leaflet motion monitored by dual camera stereo photogrammetry

Dual camera stereo photogrammetry (DCSP) was applied to investigate the leaflet motion of bioprosthetic heart valves (BHVs) in a physiologic pulse flow loop (PFL). A 25-mm bovine pericardial valve was installed in the aortic valve position of the PFL, which was operated at a pulse rate of 70 beats/min and a cardiac output of 5 l/min. The systolic/diastolic aortic pressure was maintained at 120/80 mmHg to mimic the physiologic load experienced by the aortic valve. The leaflet of the test valve was marked with 80 India ink dots to form a fan-shaped matrix. From the acquired image sequences, 3-D coordinates of the marker matrix were derived and hence the surface contour, local mean and Gaussian curvatures at each opening and closing phase during one cardiac cycle were reconstructed. It is generally believed that the long-term failure rate of BHV is related to the uneven distribution of mechanical stresses occurring in the leaflet material during opening and closing. Unfortunately, a quantitative analysis of the leaflet motion under physiological conditions has not been reported. The newly developed technique permits frame-by-frame mapping of the leaflet surface, which is essential for dynamic analysis of stress-strain behavior in BHV.     

A two-dimensional fluid-structure interaction model of the aortic valve [correction of value

Failure of synthetic heart valves is usually caused by tearing and calcification of the leaflets. Leaflet fiber-reinforcement increases the durability of these valves by unloading the delicate parts of the leaflets, maintaining their physiological functioning. The interaction of the valve with the surrounding fluid is essential when analyzing its functioning. However, the large differences in material properties of fluid and structure and the finite motion of the leaflets complicate blood-valve interaction modeling. This has, so far, obstructed numerical analyses of valves operating under physiological conditions. A two-dimensional fluid-structure interaction model is presented, which allows the Reynolds number to be within the physiological range, using a fictitious domain method based on Lagrange multipliers to couple the two phases. The extension to the three-dimensional case is straightforward. The model has been validated experimentally using laser Doppler anemometry for measuring the fluid flow and digitized high-speed video recordings to visualize the leaflet motion in corresponding geometries. Results show that both the fluid and leaflet behaviour are well predicted for different leaflet thicknesses.     

High resolution three-dimensional strain mapping of bioprosthetic heart valves using digital image correlation

Transcatheter aortic valve replacement (TAVR) is a safe and effective treatment option for patients deemed at high and intermediate risk for surgical aortic valve replacement. Similar to surgical aortic valves (SAVs), transcatheter aortic valves (TAVs) undergo calcification and mechanical wear over time. However, to date, there have been limited publications on the long-term durability of TAV devices. To assess longevity and mechanical strength of TAVs in comparison to surgical bioprosthetic valves, three-dimensional deformation analysis and strain measurement of the leaflets become an inevitable part of the evaluation. The goal of this study was to measure and compare leaflet displacement and strain of two commonly used TAVs in a side-by-side comparison with a commonly used SAV using a high-resolution digital image correlation (DIC) system. 26-mm Edwards SAPIEN 3, 26-mm Medtronic CoreValve, and 25-mm Carpentier-Edwards PERIMOUNT Magna surgical bioprosthesis were examined in a custom-made valve testing apparatus. A time-varying, spatially uniform pressure was applied to the leaflets at different loading rates. GOM ARAMIS® software was used to map leaflet displacement and strain fields during loading and unloading. High displacement regions were found to be at the leaflet belly region of the three bioprosthetic valves. In addition, the frame of the surgical bioprosthesis was found to be remarkably flexible, in contrary to CoreValve and SAPIEN 3 in which the stent was nearly rigid under a similar loading condition. The experimental DIC measurements can be used to characterize the anisotropic materiel behavior of the bioprosthetic heart valve leaflets and validate heart valve computational simulations.     

In vitro assessment of a combined radiofrequency ablation and cryo-anchoring catheter for treatment of mitral valve prolapse

Percutaneous approaches to mitral valve repair are an attractive alternative to surgical repair or replacement. Radiofrequency ablation has the potential to approximate surgical leaflet resection by using resistive heating to reduce leaflet size, and cryogenic temperatures on a percutaneous catheter can potentially be used to reversibly adhere to moving mitral valve leaflets for reliable application of radiofrequency energy. We tested a combined cryo-anchoring and radiofrequency ablation catheter using excised porcine mitral valves placed in a left heart flow loop capable of reproducing physiologic pressure and flow waveforms. Transmitral flow and pressure were monitored during the cryo-anchoring procedure and compared to baseline flow conditions, and the extent of radiofrequency energy delivery to the mitral valve was assessed post-treatment. Long term durability of radiofrequency ablation treatment was assessed using statically treated leaflets placed in a stretch bioreactor for four weeks. Transmitral flow and pressure waveforms were largely unaltered during cryo-anchoring. Parameter fitting to mechanical data from leaflets treated with radiofrequency ablation and cryo-anchoring revealed significant mechanical differences from untreated leaflets, demonstrating successful ablation of mitral valves in a hemodynamic environment. Picrosirius red staining showed clear differences in morphology and collagen birefringence between treated and untreated leaflets. The durability study indicated that statically treated leaflets did not significantly change size or mechanics over four weeks. A cryo-anchoring and radiofrequency ablation catheter can adhere to and ablate mitral valve leaflets in a physiologic hemodynamic environment, providing a possible percutaneous alternative to surgical leaflet resection of mitral valve tissue.     

Biaxial mechanical properties of bovine jugular venous valve leaflet tissues

Venous valve incompetence has been implicated in diseases ranging from chronic venous insufficiency (CVI) to intracranial venous hypertension. However, while the mechanical properties of venous valve leaflet tissues are central to CVI biomechanics and mechanobiology, neither stress–strain curves nor tangent moduli have been reported. Here, equibiaxial tensile mechanical tests were conducted to assess the tangent modulus, strength and anisotropy of venous valve leaflet tissues from bovine jugular veins. Valvular tissues were stretched to 60% strain in both the circumferential and radial directions, and leaflet tissue stress–strain curves were generated for proximal and distal valves (i.e., valves closest and furthest from the right heart, respectively). Toward linking mechanical properties to leaflet microstructure and composition, Masson’s trichrome and Verhoeff–Van Gieson staining and collagen assays were conducted. Results showed: (1) Proximal bovine jugular vein venous valves tended to be bicuspid (i.e., have two leaflets), while distal valves tended to be tricuspid; (2) leaflet tissues from proximal valves exhibited approximately threefold higher peak tangent moduli in the circumferential direction than in the orthogonal radial direction (i.e., proximal valve leaflet tissues were anisotropic; \(p<0.01\)); (3) individual leaflets excised from the same valve apparatus appeared to exhibit different mechanical properties (i.e., intra-valve variability); and (4) leaflets from distal valves exhibited a trend of higher soluble collagen concentrations than proximal ones (i.e., inter-valve variability). To the best of the authors’ knowledge, this is the first study reporting biaxial mechanical properties of venous valve leaflet tissues. These results provide a baseline for studying venous valve incompetence at the tissue level and a quantitative basis for prosthetic venous valve design.     

Functional analysis of bioprosthetic heart valves

Glutaraldehyde-treated bovine pericardium is used successfully as bioprosthetic material in the manufacturing of heart valves leaflets. The mechanical properties of bovine pericardial aortic valve leaflets seem to influence its mechanical behaviour and the failure mechanisms. In this study the effect of orthotropy on tricuspid bioprosthetic aortic valve was analysed, using a three-dimensional finite element model, during the entire cardiac cycle. Multiaxial tensile tests were also performed to determine the anisotropy of pericardium. Seven different models of the same valve were analysed using different values of mechanical characteristics from one leaflet to another, considering pericardium as an orthotropic material. The results showed that even a small difference between values along the two axes of orthotropy can negatively influence leaflets performance as regard both displacement and stress distribution. Leaflets of bovine pericardium bioprostheses could be manufactured to be similar to natural human heart valves reproducing their well-known anisotropy. In this way it could be possible to improve the manufacturing process, durability and function of pericardial bioprosthetic valves.     

Modeling of long-term fatigue damage of soft tissue with stress softening and permanent set effects

One of the major failure modes of bioprosthetic heart valves is non-calcific structural deterioration due to fatigue of the tissue leaflets. Experimental methods to characterize tissue fatigue properties are complex and time-consuming. A constitutive fatigue model that could be calibrated by isolated material tests would be ideal for investigating the effects of more complex loading conditions. However, there is a lack of tissue fatigue damage models in the literature. To address these limitations, in this study, a phenomenological constitutive model was developed to describe the stress softening and permanent set effects of tissue subjected to long-term cyclic loading. The model was used to capture characteristic uniaxial fatigue data for glutaraldehyde-treated bovine pericardium and was then implemented into finite element software. The simulated fatigue response agreed well with the experimental data and thus demonstrates feasibility of this approach.     

A note on the blood analog for in-vitro testing of heart valve bioprostheses

Glycerol solution with the viscosity coefficient similar to that of blood is used in evaluating the performance characteristics of prosthetic heart valves in the laboratory. However, physiological saline solution is used as a test fluid in testing tissue heart valves even though the viscosity coefficient does not match that of human blood. It is commonly believed that glycerol is absorbed by the tissue valves and hence the leaflets become stiff, making the test results invalid. However, in our laboratory a comparison of tissue valves exposed to glycerine solution at various times does not indicate any difference in the leaflet opening characteristics. Hence, it is suggested that glycerine solution be used as a test fluid for the evaluation of tissue valves also.     

Modeling leaflet correction techniques in aortic valve repair: A finite element study

In aortic valve sparing surgery, cusp prolapse is a common cause of residual aortic insufficiency. To correct cusp pathology, native leaflets of the valve frequently require adjustment which can be performed using a variety of described correction techniques, such as central or commissural plication, or resuspension of the leaflet free margin. The practical question then arises of determining which surgical technique provides the best valve performance with the most physiologic coaptation. To answer this question, we created a new finite element model with the ability to simulate physiologic function in normal valves, and aortic insufficiency due to leaflet prolapse in asymmetric, diseased or sub-optimally repaired valves. The existing leaflet correction techniques were simulated in a controlled situation, and the performance of the repaired valve was quantified in terms of maximum leaflets stress, valve orifice area, valve opening and closing characteristics as well as total coaptation area in diastole. On the one hand, the existing leaflet correction techniques were shown not to adversely affect the dynamic properties of the repaired valves. On the other hand, leaflet resuspension appeared as the best technique compared to central or commissural leaflet plication. It was the only method able to achieve symmetric competence and fix an individual leaflet prolapse while simultaneously restoring normal values for mechanical stress, valve orifice area and coaptation area.     

A nonlinear anisotropic model for porcine aortic heart valves

The anisotropic property of porcine aortic valve leaflet has potentially significant effects on its mechanical behaviour and the failure mechanisms. However, due to its complex nature, testing and modelling the anisotropic porcine aortic valves remains a continuing challenge to date. This study has developed a nonlinear anisotropic finite element model for porcine heart valves. The model is based on the uniaxial experimental data of porcine aortic heart valve leaflet and the properties of nonlinear composite material. A finite element code is developed to solve this problem using the 8-node super-parameter nonlinear shells and the update Lagrangian method. The stress distribution and deformation of the porcine aortic valves with either uniform and non-uniform thicknesses in closed phase and loaded condition are calculated. The results showed significant changes in the stress distributions due to the anisotropic property of the leaflets. Compared with the isotropic valve at the same loading condition, it is found that the site of the peak stress of the anisotropic leaflet is different; the maximum longitudinal normal stress is increased, but the maximum transversal normal stress and in-plane shear stress are reduced. We conclude that it is very important to consider the anisotropic property of the porcine heart valves in order to understand the failure mechanism of such valves in vivo.     

Aortic valve leaflet mechanical properties facilitate diastolic valve function

This work was concerned with the numerical simulation of the behaviour of aortic valves whose material can be modelled as non-linear elastic anisotropic. Linear elastic models for the valve leaflets with parameters used in previous studies were compared with hyperelastic models, incorporating leaflet anisotropy with pronounced stiffness in the circumferential direction through a transverse isotropic model. The parameters for the hyperelastic models were obtained from fits to results of orthogonal uniaxial tensile tests on porcine aortic valve leaflets. The computational results indicated the significant impact of transverse isotropy and hyperelastic effects on leaflet mechanics; in particular, increased coaptation with peak values of stress and strain in the elastic limit. The alignment of maximum principal stresses in all models follows approximately the coarse collagen fibre distribution found in aortic valve leaflets. The non-linear elastic leaflets also demonstrated more evenly distributed stress and strain which appears relevant to long-term scaffold stability and mechanotransduction.     

Platelet activation of mechanical versus bioprosthetic heart valves during systole

Thrombus formation is a major concern for recipients of mechanical heart valves (MHVs), which requires them to take anticoagulant drugs for the rest of their lives. Bioprosthetic heart valves (BHVs) do not require life-long anticoagulant therapy but deteriorate after 10–15 years. The thrombus formation is initiated by the platelet activation which is thought to be mainly generated in MHVs by the flow through the hinge and the leakage flow during the diastole. However, our results show that the activation in the bulk flow during the systole phase might play an essential role as well. This is based on our results obtained by comparing the thrombogenic performance of a MHV and a BHV (as control) in terms of shear induced platelet activation under exactly the same conditions. Three different mathematical activation models including linear level of activation, damage accumulation, and Soares model are tested to quantify the platelet activation during systole using the previous simulations of the flow through MHV and BHV in a straight aorta under the same physiologic flow conditions. Results indicate that the platelet activation in the MHV at the beginning of the systole phase is slightly less than the BHV. However, at the end of the systole phase the platelet activation by the bulk flow for the MHV is several folds (1.41, 5.12, and 2.81 for linear level of activation, damage accumulation, and Soares model, respectively) higher than the BHV for all tested platelet activation models.     

Straightening of curved pattern of collagen fibers under load controls aortic valve shape

The network of collagen fibers in the aortic valve leaflet is believed to play an important role in the strength and durability of the valve. However, in addition to its stress-bearing role, such a fiber network has the potential to produce functionally important shape changes in the closed valve under pressure load. We measured the average pattern of the collagen network in porcine aortic valve leaflets after staining for collagen. We then used finite element simulation to explore how this collagen pattern influences the shape of the closed valve. We observed a curved or bent pattern, with collagen fibers angled downward from the commissures toward the center of the leaflet to form a pattern that is concave toward the leaflet free edge. Simulations showed that these curved fiber trajectories straighten under pressure load, leading to functionally important changes in closed valve shape. Relative to a pattern of straight collagen fibers running parallel to the leaflet free edge, the concave pattern of curved fibers produces a closed valve with a 40% increase in central leaflet coaptation height and with decreased leaflet billow, resulting in a more physiological closed valve shape. Furthermore, simulations show that these changes in loaded leaflet shape reflect changes in leaflet curvature due to modulation of in-plane membrane stress resulting from straightening of the curved fibers. This effect appears to play an important role in normal valve function and may have important implications for the design of prosthetic and tissue engineered replacement valves.     

Pericardial heterograft valves: An assessment of leaflet stresses and their implications for heart valve design

Bovine pericardium, stabilized with glutaraldehyde, is used widely in the construction of heart valve substitutes, but the design and construction of valve substitutes from this material are empirically based. Collagenous tissue can support tension, but experimental evidence indicates that flexure-induced compressive stresses can lead to fatigue failure. This study uses experimental results obtained from cyclic uniaxial load tests to predict the type and magnitude of operational stresses which occur in pericardial heterograft leaflets. Both Young's modulus and Poisson's ratio varied with uniaxial loading in pericardium, chemically modified free of tension. Leaflet stresses were analysed by using effective incremental representations of these parameters. In leaflets with unrestricted rotation at the point of attachment to the stent, the mid-plane tensions always exceeded the bending stresses, and no zones of leaflet compression were predicted. In contrast, with totally restricted leaflet rotation induced by clamping (possibly between a male and female frame) the bending stresses were greater than the mid-plane tensions at the hinge line and significant compressive stresses were predicted at this site. If elastic boundary conditions were introduced at the stent (possibly by wrapping the stent in pericardium) then the compressive stresses were reduced as the degree of elasticity was increased. Glutaraldehyde fixation of the pericardium under load produced a stiffer material; higher compressive stresses at the stent and significant increases in total stress were predicted for this tissue. The application of elevated pressure loading also increased the compressive and total stresses in the leaflet. Finally, it was shown that bicuspid leaflets were likely to experience higher stresses than tricuspid leaflets. This simple stress analysis should help valve designers of pericardial heterografts to identify those conditions which lead to tissue compression, high total stress, and ultimately material fatigue.     

Cell-mediated retraction versus hemodynamic loading – A delicate balance in tissue-engineered heart valves

Preclinical studies of tissue-engineered heart valves (TEHVs) showed retraction of the heart valve leaflets as major failure of function mechanism. This retraction is caused by both passive and active cell stress and passive matrix stress. Cell-mediated retraction induces leaflet shortening that may be counteracted by the hemodynamic loading of the leaflets during diastole. To get insight into this stress balance, the amount and duration of stress generation in engineered heart valve tissue and the stress imposed by physiological hemodynamic loading are quantified via an experimental and a computational approach, respectively.     

Comparative study of cellular and extracellular matrix composition of native and tissue engineered heart valves.

Tissue engineering of heart valves utilizes biodegradable or metabolizable scaffolds for remodeling by seeded autologous cells. The aim of this study was to determine and compare extracellular matrix (ECM) formations, cellular phenotypes and cell location of native and tissue engineered (TE) valve leaflets. Ovine carotid arteries, ovine and porcine hearts were obtained from slaughterhouses. Cells were isolated from carotid arteries and dissected ovine, porcine and TE leaflets. TE constructs were fabricated from decellularized porcine pulmonary valves, seeded ovine arterial cells and subsequent 16 days dynamic in vitro culture using a pulsatile bioreactor. Native and TE valves were studied by histology (hematoxylin-eosin, resorcin-fuchsin, Movat pentachrome), NIR femtosecond multiphoton laser scanning microscopy and scanning electron microscopy (SEM). Cells of native and TE tissues were identified and localized by immunohistochemistry. Arterial, valvular and re-isolated TE-construct cells were processed for immunocytochemistry and Western blotting. ECM analysis and SEM revealed characteristical and comparable structures in native and TE leaflets. Most cells in native leaflets stained strongly positive for vimentin. Cells positive to alpha-smooth muscle actin (alpha-SMA), myosin and calponin were only found at the ventricular (inflow) side of ovine aortic and porcine pulmonary valve leaflets. Cells from TE constructs had a strong expression of vimentin, alpha-SMA, myosin, calponin and h-caldesmon throughout the entire leaflet. Comparable ECM formation and endothelial cell lining of native and TE leaflets could be demonstrated. However, immunostaining revealed significant differences between valvular cell phenotypes of native and TE leaflets. These results may be essential for further cardiovascular tissue engineering efforts.