Architectural trends in the human normal and bicuspid aortic valve leaflet and its relevance to valve disease

An Inverse Modeling Approach to Estimate Three-Dimensional Aortic Valve Interstitial Cell Stress Fiber Force Levels

Within the aortic valve (AV) leaflet exists a population of interstitial cells (AVICs) that maintain the constituent tissues by extracellular matrix (ECM) secretion, degradation, and remodeling. AVICs can transition from a quiescent, fibroblast-like phenotype to an activated, myofibroblast phenotype in response to growth or disease. AVIC dysfunction has been implicated in AV disease processes, yet our understanding of AVIC function remains quite limited. A major characteristic of the AVIC phenotype is its contractile state, driven by contractile forces generated by the underlying stress fibers (SF). However, direct assessment of the AVIC SF contractile state and structure within physiologically mimicking three-dimensional environments remains technically challenging, as the size of single SFs are below the resolution of light microscopy. Therefore, in the present study, we developed a three-dimensional (3D) computational approach of AVICs embedded in 3D hydrogels to estimate their SF local orientations and contractile forces. One challenge with this approach is that AVICs will remodel the hydrogel, so that the gel moduli will vary spatially. We thus utilized our previous approach (Khang et al. 2023, "Estimation of Aortic Valve Interstitial Cell-Induced 3D Remodeling of Poly (Ethylene Glycol) Hydrogel Environments Using an Inverse Finite Element Approach," Acta Biomater., 160, pp. 123-133) to define local hydrogel mechanical properties. The AVIC SF model incorporated known cytosol and nucleus mechanical behaviors, with the cell membrane assumed to be perfectly bonded to the surrounding hydrogel. The AVIC SFs were first modeled as locally unidirectional hyperelastic fibers with a contractile force component. An adjoint-based inverse modeling approach was developed to estimate local SF orientation and contractile force. Substantial heterogeneity in SF force and orientations were observed, with the greatest levels of SF alignment and contractile forces occurring in AVIC protrusions. The addition of a dispersed SF orientation to the modeling approach did not substantially alter these findings. To the best of our knowledge, we report the first fully 3D computational contractile cell models which can predict locally varying stress fiber orientation and contractile force levels.

Calcium Fracture and Device Over Expansion in Transcatheter Aortic Valve Replacement for Bicuspid Aortic Valves

Transcatheter aortic valve replacement (TAVR) in patients with bicuspid aortic valve disease (BAV) has potential risks of under expansion and non-circularity which may compromise long-term durability. This study aims to investigate calcium fracture and balloon over expansion in balloon-expandable TAVs on the stent deformation with the aid of simulation. BAV patients treated with the SAPIEN 3 Ultra with pre- and post-TAVR CTs were analyzed (n = 8). Simulations of the stent deployment were performed (1) with baseline simulation allowing calcium fracture, (2) without allowable calcium fracture and (3) with balloon over expansion (1 mm larger diameter). When compared to post CT, baseline simulations had minimal error in expansion (2.5% waist difference) and circularity (3.0% waist aspect ratio difference). When compared to baseline, calcium fracture had insignificant impact on the expansion (- 0.5% average waist difference) and circularity (- 1.6% average waist aspect ratio difference). Over expansion had significantly larger expansion compared to baseline (15.4% average waist difference) but had insignificant impact on the circularity (- 0.5% waist aspect ratio difference). We conclude that stent deformation can be predicted with minimal error, calcium fracture has small differences on the final stent deformation except in extreme calcified cases, and balloon over expansion expands the waist closer to nominal values.

Functional differences in human aortic valve interstitial cells from patients with varying calcific aortic valve disease

Calcific aortic valve disease (CAVD) is characterized by progressive stiffening of aortic valve (AV) tissues, inducing stenosis and insufficiency. Bicuspid aortic valve (BAV) is a common congenital defect in which the AV has two leaflets rather than three, with BAV patients developing CAVD decades years earlier than in the general population. Current treatment for CAVD remains surgical replacement with its continued durability problems, as there are no pharmaceutical therapies or other alternative treatments available. Before such therapeutic approaches can be developed, a deeper understanding of CAVD disease mechanisms is clearly required. It is known that AV interstitial cells (AVICs) maintain the AV extracellular matrix and are typically quiescent in the normal state, transitioning into an activated, myofibroblast-like state during periods of growth or disease. One proposed mechanism of CAVD is the subsequent transition of AVICs into an osteoblast-like phenotype. A sensitive indicator of AVIC phenotypic state is enhanced basal contractility (tonus), so that AVICs from diseased AV will exhibit a higher basal tonus level. The goals of the present study were thus to assess the hypothesis that different human CAVD states lead to different biophysical AVIC states. To accomplish this, we characterized AVIC basal tonus behaviors from diseased human AV tissues embedded in 3D hydrogels. Established methods were utilized to track AVIC-induced gel displacements and shape changes after the application of Cytochalasin D (an actin polymerization inhibitor) to depolymerize the AVIC stress fibers. Results indicated that human diseased AVICs from the non-calcified region of TAVs were significantly more activated than AVICs from the corresponding calcified region. In addition, AVICs from the raphe region of BAVs were more activated than from the non-raphe region. Interestingly, we observed significantly greater basal tonus levels in females compared to males. Furthermore, the overall AVIC shape changes after Cytochalasin suggested that AVICs from TAVs and BAVs develop different stress fiber architectures. These findings are the first evidence of sex-specific differences in basal tonus state in human AVICs in varying disease states. Future studies are underway to quantify stress fiber mechanical behaviors to further elucidate CAVD disease mechanisms.

The effects of strain history on aortic valve interstitial cell activation in a 3D hydrogel environment

Aortic valves (AVs) undergo unique stretch histories that include high rates and magnitudes. While major differences in deformation patterns have been observed between normal and congenitally defective bicuspid aortic valves (BAVs), the relation to underlying mechanisms of rapid disease onset in BAV patients remains unknown. To evaluate how the variations in stretch history affect AV interstitial cell (AVIC) activation, high-throughput methods were developed to impart varied cyclical biaxial stretch histories into 3D poly(ethylene) glycol hydrogels seeded with AVICs for 48 h. Specifically, a physiologically mimicking stretch history was compared to two stretch histories with varied peak stretch and stretch rate. Post-conditioned AVICs were imaged for nuclear shape, alpha smooth muscle actin (αSMA) and vimentin (VMN) polymerization, and small mothers against decapentaplegic homologs 2 and 3 (SMAD 2/3) nuclear activity. The results indicated that bulk gel deformations were accurately transduced to the AVICs. Lower peak stretches lead to increased αSMA polymerization. In contrast, VMN polymerization was a function of stretch rate, with SMAD 2/3 nuclear localization and nuclear shape also trending toward stretch rate dependency. Lower than physiological levels of stretch rate led to higher SMAD 2/3 activity, higher VMN polymerization around the nucleus, and lower nuclear elongation. αSMA polymerization did not correlate with VMN polymerization, SMAD 2/3 activity, nor nuclear shape. These results suggest that a negative feedback loop may form between SMAD 2/3, VMN, and nuclear shape to maintain AVIC homeostatic nuclear deformations, which is dependent on stretch rate. These novel results suggest that AVIC mechanobiological responses are sensitive to stretch history and provide insight into the mechanisms of AV disease.

Neural Network Approaches for Soft Biological Tissue and Organ Simulations

Given the functional complexities of soft tissues and organs, it is clear that computational simulations are critical in their understanding and for the rational basis for the development of therapies and replacements. A key aspect of such simulations is accounting for their complex, nonlinear, anisotropic mechanical behaviors. While soft tissue material models have developed to the point of high fidelity, in-silico implementation is typically done using the finite element (FE) method, which remains impractically slow for translational clinical time frames. As a potential path toward addressing the development of high fidelity simulations capable of performing in clinically relevant time frames, we review the use of neural networks (NN) for soft tissue and organ simulation using two approaches. In the first approach, we show how a NN can learn the responses for a detailed meso-structural soft tissue material model. The NN material model not only reproduced the full anisotropic mechanical responses but also demonstrated a considerable efficiency improvement, as it was trained over a range of realizable fibrous structures. In the second approach, we go a step further with the use of a physics-based surrogate model to directly learn the displacement field solution without the need for raw training data or FE simulation datasets. In this approach we utilize a finite element mesh to define the domain and perform the necessary integrations, but not the finite element method (FEM) itself. We demonstrate with this approach, termed neural network finite element (NNFE), results in a trained NNFE model with excellent agreement with the corresponding "ground truth" FE solutions over the entire physiological deformation range on a cuboidal myocardium specimen. More importantly, the NNFE approach provided a significantly decreased computational time for a range of finite element mesh sizes. Specifically, as the FE mesh size increased from 2744 to 175,615 elements, the NNFE computational time increased from 0.1108 s to 0.1393 s, while the "ground truth" FE model increased from 4.541 s to 719.9 s, with the same effective accuracy. These results suggest that NNFE run times are significantly reduced compared with the traditional large-deformation-based finite element solution methods. We then show how a nonuniform rational B-splines (NURBS)-based approach can be directly integrated into the NNFE approach as a means to handle real organ geometries. While these and related approaches are in their early stages, they offer a method to perform complex organ-level simulations in clinically relevant time frames without compromising accuracy.

Biomechanics of Transcatheter Aortic Valve Replacement Complications and Computational Predictive Modeling

Transcatheter aortic valve replacement (TAVR) is a rapidly growing field enabling replacement of diseased aortic valves without the need for open heart surgery. However, due to the nature of the procedure and nonremoval of the diseased tissue, there are rates of complications ranging from tissue rupture and coronary obstruction to paravalvular leak, valve thrombosis, and permanent pacemaker implantation. In recent years, computational modeling has shown a great deal of promise in its capabilities to understand the biomechanical implications of TAVR as well as help preoperatively predict risks inherent to device-patient-specific anatomy biomechanical interaction. This includes intricate replication of stent and leaflet designs and tested and validated simulated deployments with structural and fluid mechanical simulations. This review outlines current biomechanical understanding of device-related complications from TAVR and related predictive strategies using computational modeling. An outlook on future modeling strategies highlighting reduced order modeling which could significantly reduce the high time and cost that are required for computational prediction of TAVR outcomes is presented in this review paper. A summary of current commercial/in-development software is presented in the final section.

Atypical Expression of Smooth Muscle Markers and Co-activators and Their Regulation in Rheumatic Aortic and Calcified Bicuspid Valves

Objective

We have previously reported that human calcified aortic cusps have abundant expression of smooth muscle (SM) markers and co-activators. We hypothesised that cells in bicuspid aortic valve (BAV) cusps and those affected by rheumatic heart valve (RHV) disease may follow a similar phenotypic transition into smooth muscle cells, a process that could be regulated by transforming growth factors (TGFs).

Aims

Cusps from eight patients with BAV and seven patients with RHV were analysed for early and late SM markers and regulators of SM gene expression by immunocytochemistry and compared to healthy aortic valves from 12 unused heart valve donors. The ability of TGFs to induce these markers in valve endothelial cells (VECs) on two substrates was assessed.

Results

In total, 7 out of 8 BAVs and all the RHVs showed an increased and atypical expression of early and late SM markers α-SMA, calponin, SM22 and SM-myosin. The SM marker co-activators were aberrantly expressed in six of the BAV and six of the RHV, in a similar regional pattern to the expression of SM markers. Additionally, regions of VECs, and endothelial cells lining the vessels within the cusps were found to be positive for SM markers and co-activators in three BAV and six RHV. Both BAVs and RHVs were significantly thickened and HIF1α expression was prominent in four BAVs and one RHV. The ability of TGFβs to induce the expression of SM markers and myocardin was greater in VECs cultured on fibronectin than on gelatin. Fibronectin was shown to be upregulated in BAVs and RHVs, within the cusps as well as in the basement membrane.

Conclusion

Bicuspid aortic valves and RHVs expressed increased numbers of SM marker-positive VICs and VECs. Concomittantly, these cells expressed MRTF-A and myocardin, key regulators of SM gene expression. TGFβ1 was able to preferentially upregulate SM markers and myocardin in VECs on fibronectin, and fibronectin was found to be upregulated in BAVs and RHVs. These findings suggest a role of VEC as a source of cells that express SM cell markers in BAVs and RHVs. The similarity between SM marker expression in BAVs and RHVs with our previous study with cusps from patients with aortic stenosis suggests the existance of a common pathological pathway between these different pathologies.

Patient-Specific Quantification of Normal and Bicuspid Aortic Valve Leaflet Deformations from Clinically Derived Images

The clinical benefit of patient-specific modeling of heart valve disease remains an unrealized goal, often a result of our limited understanding of the in vivo milieu. This is particularly true in assessing bicuspid aortic valve (BAV) disease, the most common cardiac congenital defect in humans, which leads to premature and severe aortic stenosis or insufficiency (AS/AI). However, assessment of BAV risk for AS/AI on a patient-specific basis is hampered by the substantial degree of anatomic and functional variations that remain largely unknown. The present study was undertaken to utilize a noninvasive computational pipeline ( https://doi.org/10.1002/cnm.3142 ) that directly yields local heart valve leaflet deformation information using patient-specific real-time three-dimensional echocardiographic imaging (rt-3DE) data. Imaging data was collected for patients with normal tricuspid aortic valve (TAV, [Formula: see text]) and those with BAV ([Formula: see text] with fused left and right coronary leaflets and [Formula: see text] with fused right and non-coronary leaflets), from which the medial surface of each leaflet was extracted. The resulting deformation analysis resulted in, for the first time, quantified differences between the in vivo functional deformations of the TAV and BAV leaflets. Our approach was able to capture the complex, heterogeneous surface deformation fields in both TAV and BAV leaflets. We were able to identify and quantify differences in stretch patterns between leaflet types, and found in particular that stretches experienced by BAV leaflets during closure differ from those of TAV leaflets in terms of both heterogeneity as well as overall magnitude. Deformation is a key parameter in the clinical assessment of valvular function, and serves as a direct means to determine regional variations in structure and function. This study is an essential step toward patient-specific assessment of BAV based on correlating leaflet deformation and AS/AI progression, as it provides a means for assessing patient-specific stretch patterns.

Parameterization, geometric modeling, and isogeometric analysis of tricuspid valves

Approximately 1.6 million patients in the United States are affected by tricuspid valve regurgitation, which occurs when the tricuspid valve does not close properly to prevent backward blood flow into the right atrium. Despite its critical role in proper cardiac function, the tricuspid valve has received limited research attention compared to the mitral and aortic valves on the left side of the heart. As a result, proper valvular function and the pathologies that may cause dysfunction remain poorly understood. To promote further investigations of the biomechanical behavior and response of the tricuspid valve, this work establishes a parameter-based approach that provides a template for tricuspid valve modeling and simulation. The proposed tricuspid valve parameterization presents a comprehensive description of the leaflets and the complex chordae tendineae for capturing the typical three-cusp structural deformation observed from medical data. This simulation framework develops a practical procedure for modeling tricuspid valves and offers a robust, flexible approach to analyze the performance and effectiveness of various valve configurations using isogeometric analysis. The proposed methods also establish a baseline to examine the tricuspid valve's structural deformation, perform future investigations of native valve configurations under healthy and disease conditions, and optimize prosthetic valve designs.

Transcatheter Heart Valves: A Biomaterials Perspective

Heart valve disease is prevalent throughout the world, and the number of heart valve replacements is expected to increase rapidly in the coming years. Transcatheter heart valve replacement (THVR) provides a safe and minimally invasive means for heart valve replacement in high-risk patients. The latest clinical data demonstrates that THVR is a practical solution for low-risk patients. Despite these promising results, there is no long-term (>20 years) durability data on transcatheter heart valves (THVs), raising concerns about material degeneration and long-term performance. This review presents a detailed account of the materials development for THVRs. It provides a brief overview of THVR, the native valve properties, the criteria for an ideal THV, and how these devices are tested. A comprehensive review of materials and their applications in THVR, including how these materials are fabricated, prepared, and assembled into THVs is presented, followed by a discussion of current and future THVR biomaterial trends. The field of THVR is proliferating, and this review serves as a guide for understanding the development of THVs from a materials science and engineering perspective.

A Fluid-Structure Interaction Study of Different Bicuspid Aortic Valve Phenotypes Throughout the Cardiac Cycle

The bicuspid aortic valve (BAV) is a congenital malformation of the aortic valve with a variety of structural features. The current research on BAV mainly focuses on the systolic phase, while ignoring the diastolic hemodynamic characteristics and valve mechanics. The purpose of this study is to compare the differences in hemodynamics and mechanical properties of BAV with different phenotypes throughout the cardiac cycle by means of numerical simulation. Based on physiological anatomy, we established an idealized tricuspid aortic valve (TAV) model and six phenotypes of BAV models (including Type 0 a-p, Type 0 lat, Type 1 L-R, Type 1 N-L, Type 1 R-N, and Type 2), and simulated the dynamic changes of the aortic valve during the cardiac cycle using the fluid-structure interaction method. The morphology of the leaflets, hemodynamic parameters, flow patterns, and strain were analyzed. Compared with TAV, the cardiac output and effective orifice area of different BAV phenotypes decreased certain degree, along with the peak velocity and mean pressure difference increased both. Among all BAV models, Type 2 exhibited the worst hemodynamic performance. During the systole, obvious asymmetric flow field was observed in BAV aorta, which was related to the orientation of BAV. Higher strain was generated in diastole for BAV models. The findings of this study suggests specific differences in the hemodynamic characteristics and valve mechanics of different BAV phenotypes, including different severity of stenosis, flow patterns, and leaflet strain, which may be critical for prediction of other subsequent aortic diseases and differential treatment strategy for certain BAV phenotype.

Isogeometric finite element-based simulation of the aortic heart valve: Integration of neural network structural material model and structural tensor fiber architecture representations

The functional complexity of native and replacement aortic heart valves (AVs) is well known, incorporating such physical phenomenons as time-varying non-linear anisotropic soft tissue mechanical behavior, geometric non-linearity, complex multi-surface time varying contact, and fluid-structure interactions to name a few. It is thus clear that computational simulations are critical in understanding AV function and for the rational basis for design of their replacements. However, such approaches continued to be limited by ad-hoc approaches for incorporating tissue fibrous structure, high-fidelity material models, and valve geometry. To this end, we developed an integrated tri-leaflet valve pipeline built upon an isogeometric analysis framework. A high-order structural tensor (HOST)-based method was developed for efficient storage and mapping the two-dimensional fiber structural data onto the valvular 3D geometry. We then developed a neural network (NN) material model that learned the responses of a detailed meso-structural model for exogenously cross-linked planar soft tissues. The NN material model not only reproduced the full anisotropic mechanical responses but also demonstrated a considerable efficiency improvement, as it was trained over a range of realizable fibrous structures. Results of parametric simulations were then performed, as well as population-based bicuspid AV fiber structure, that demonstrated the efficiency and robustness of the present approach. In summary, the present approach that integrates HOST and NN material model provides an efficient computational analysis framework with increased physical and functional realism for the simulation of native and replacement tri-leaflet heart valves.

Patient-Specific Bicuspid Aortic Valve Biomechanics: A Magnetic Resonance Imaging Integrated Fluid-Structure Interaction Approach

Congenital bicuspid aortic valve (BAV) consists of two fused cusps and represents a major risk factor for calcific valvular stenosis. Herein, a fully coupled fluid-structure interaction (FSI) BAV model was developed from patient-specific magnetic resonance imaging (MRI) and compared against in vivo 4-dimensional flow MRI (4D Flow). FSI simulation compared well with 4D Flow, confirming direction and magnitude of the flow jet impinging onto the aortic wall as well as location and extension of secondary flows and vortices developing at systole: the systolic flow jet originating from an elliptical 1.6 cm² orifice reached a peak velocity of 252.2 cm/s, 0.6% lower than 4D Flow, progressively impinging on the ascending aorta convexity. The FSI model predicted a peak flow rate of 22.4 L/min, 6.7% higher than 4D Flow, and provided BAV leaflets mechanical and flow-induced shear stresses, not directly attainable from MRI. At systole, the ventricular side of the non-fused leaflet revealed the highest wall shear stress (WSS) average magnitude, up to 14.6 Pa along the free margin, with WSS progressively decreasing towards the belly. During diastole, the aortic side of the fused leaflet exhibited the highest diastolic maximum principal stress, up to 322 kPa within the attachment region. Systematic comparison with ground-truth non-invasive MRI can improve the computational model ability to reproduce native BAV hemodynamics and biomechanical response more realistically, and shed light on their role in BAV patients' risk for developing complications; this approach may further contribute to the validation of advanced FSI simulations designed to assess BAV biomechanics.

The role of stress concentration in calcified bicuspid aortic valve

Calcific aortic valve disease (CAVD) is the most common valvular heart disease in the aging population, and is now believed to be a slow, progressive, yet actively regulated process. The disease progression can be divided into two phases: initiation phase associated with lipid deposition and inflammation response, and the later propagation phase with active calcification growth. It has been hypothesized that elevated mechanical stress plays a major role in both phases of disease progression. In order to identify a direct link between leaflet stress and calcification development, we performed patient-specific finite-element (FE) analyses of six bicuspid aortic valves (BAV), where the leaflets, raphe and calcifications were all considered. The results showed that during the initiation phase, calcium buildup is likely to occur along the leaflet-root attachment curve (ATC), and the commissures, which are subject to the most drastic changes in stress during the cardiac cycle. During the propagation phase, the presence of calcification would lead to local stress concentration along its boundary, hence further calcification growth. Three patterns of calcification formation were identified on BAV leaflets: 'radial', which extended radially from ATC into the leaflet belly region; 'commissure to commissure', which extended circumferentially along the coaptation; and 'raphe', which located in the vicinity of the raphe. Furthermore, we found a strong correlation between regions with a high risk of calcium buildup and regions with elevated mechanical stress. The high-risk regions predicted at diastole on the non-calcified leaflet from FE models agreed reasonably well with the in vivo calcification locations, which indicates that patient-specific FE modelling could help us to evaluate the potential risk of calcification formation in the early stage of CAVD.

Effect of Residual and Transformation Choice on Computational Aspects of Biomechanical Parameter Estimation of Soft Tissues

Several nonlinear and anisotropic constitutive models have been proposed to describe the biomechanical properties of soft tissues, and reliably estimating the unknown parameters in these models using experimental data is an important step towards developing predictive capabilities. However, the effect of parameter estimation technique on the resulting biomechanical parameters remains under-analyzed. Standard off-the-shelf techniques can produce unreliable results where the parameters are not uniquely identified and can vary with the initial guess. In this study, a thorough analysis of parameter estimation techniques on the resulting properties for four multi-parameter invariant-based constitutive models is presented. It was found that linear transformations have no effect on parameter estimation for the presented cases, and nonlinear transforms are necessary for any improvement. A distinct focus is put on the issue of non-convergence, and we propose simple modifications that not only improve the speed of convergence but also avoid convergence to a wrong solution. The proposed modifications are straightforward to implement and can avoid severe problems in the biomechanical analysis. The results also show that including the fiber angle as an unknown in the parameter estimation makes it extremely challenging, where almost all of the formulations and models fail to converge to the true solution. Therefore, until this issue is resolved, a non-mechanical—such as optical—technique for determining the fiber angle is required in conjunction with the planar biaxial test for a robust biomechanical analysis.

Mechanics of the Tricuspid Valve-From Clinical Diagnosis/Treatment, In-Vivo and In-Vitro Investigations, to Patient-Specific Biomechanical Modeling

Proper tricuspid valve (TV) function is essential to unidirectional blood flow through the right side of the heart. Alterations to the tricuspid valvular components, such as the TV annulus, may lead to functional tricuspid regurgitation (FTR), where the valve is unable to prevent undesired backflow of blood from the right ventricle into the right atrium during systole. Various treatment options are currently available for FTR; however, research for the tricuspid heart valve, functional tricuspid regurgitation, and the relevant treatment methodologies are limited due to the pervasive expectation among cardiac surgeons and cardiologists that FTR will naturally regress after repair of left-sided heart valve lesions. Recent studies have focused on (i) understanding the function of the TV and the initiation or progression of FTR using both in-vivo and in-vitro methods, (ii) quantifying the biomechanical properties of the tricuspid valve apparatus as well as its surrounding heart tissue, and (iii) performing computational modeling of the TV to provide new insight into its biomechanical and physiological function. This review paper focuses on these advances and summarizes recent research relevant to the TV within the scope of FTR. Moreover, this review also provides future perspectives and extensions critical to enhancing the current understanding of the functioning and remodeling tricuspid valve in both the healthy and pathophysiological states.

Serotonin receptor 2B signaling with interstitial cell activation and leaflet remodeling in degenerative mitral regurgitation

AIMS

Mitral valve interstitial cells (MVIC) play an important role in the pathogenesis of degenerative mitral regurgitation (MR) due to mitral valve prolapse (MVP). Numerous clinical studies have observed serotonin (5HT) dysregulation in cardiac valvulopathies; however, the impact of 5HT-mediated signaling on MVIC activation and leaflet remodeling in MVP have been investigated to a limited extent. Here we test the hypothesis that 5HT receptors (5HTRs) signaling contributes to MVP pathophysiology.

METHODS AND RESULTS

Diseased human MV leaflets were obtained during cardiac surgery for MVP; normal MV leaflets were obtained from heart transplants. MV RNA was used for microarray analysis of MVP patients versus control, highlighting genes that indicate the involvement of 5HTR pathways and extracellular matrix remodeling in MVP. Human MV leaflets were also studied in vitro and ex vivo with biomechanical testing to assess remodeling in the presence of a 5HTR2B antagonist (LY272015). MVP leaflets from Cavalier King Charles Spaniels were used as a naturally acquired in vivo model of MVP. These canine MVP leaflets (N=5/group) showed 5HTR2B upregulation. This study also utilized CB57.1ML/6 mice in order to determine the effect of Angiotensin II infusion on MV remodeling. Histological analysis showed that MV thickening due to chronic Angiotensin II remodeling is mitigated by a 5HTR2B antagonist (LY272015) but not by 5HTR2A inhibitors.

CONCLUSION

In humans, MVP is associated with an upregulation in 5HTR2B expression and increased 5HT receptor signaling in the leaflets. Antagonism of 5HTR2B mitigates MVIC activation in vitro and MV remodeling in vivo. These observations support the view that 5HTR signaling is involved not only in previously reported 5HT-related valvulopathies, but it is also involved in the pathological remodeling of MVP.

Computational methods for the aortic heart valve and its replacements

INTRODUCTION

Replacement with a prosthetic device remains a major treatment option for the patients suffering from heart valve disease, with prevalence growing resulting from an ageing population. While the most popular replacement heart valve continues to be the bioprosthetic heart valve (BHV), its durability remains limited. There is thus a continued need to develop a general understanding of the underlying mechanisms limiting BHV durability to facilitate development of a more durable prosthesis. In this regard, computational models can play a pivotal role as they can evaluate our understanding of the underlying mechanisms and be used to optimize designs that may not always be intuitive. Areas covered: This review covers recent progress in computational models for the simulation of BHV, with a focus on aortic valve (AV) replacement. Recent contributions in valve geometry, leaflet material models, novel methods for numerical simulation, and applications to BHV optimization are discussed. This information should serve not only to infer reliable and dependable BHV function, but also to establish guidelines and insight for the design of future prosthetic valves by analyzing the influence of design, hemodynamics and tissue mechanics. Expert commentary: The paradigm of predictive modeling of heart valve prosthesis are becoming a reality which can simultaneously improve clinical outcomes and reduce costs. It can also lead to patient-specific valve design.

An improved parameter estimation and comparison for soft tissue constitutive models containing an exponential function

Motivated by the well-known result that stiffness of soft tissue is proportional to the stress, many of the constitutive laws for soft tissues contain an exponential function. In this work, we analyze properties of the exponential function and how it affects the estimation and comparison of elastic parameters for soft tissues. In particular, we find that as a consequence of the exponential function there are lines of high covariance in the elastic parameter space. As a result, one can have widely varying mechanical parameters defining the tissue stiffness but similar effective stress–strain responses. Drawing from elementary algebra, we propose simple changes in the norm and the parameter space, which significantly improve the convergence of parameter estimation and robustness in the presence of noise. More importantly, we demonstrate that these changes improve the conditioning of the problem and provide a more robust solution in the case of heterogeneous material by reducing the chances of getting trapped in a local minima. Based upon the new insight, we also propose a transformed parameter space which will allow for rational parameter comparison and avoid misleading conclusions regarding soft tissue mechanics.

Computational comparison of regional stress and deformation characteristics in tricuspid and bicuspid aortic valve leaflets

The bicuspid aortic valve (BAV) is the most common congenital valvular defect and a major risk factor for secondary calcific aortic valve disease. While hemodynamics is presumed to be a potential contributor to this complication, the validation of this theory has been hampered by the limited knowledge of the mechanical stress abnormalities experienced by BAV leaflets and their dependence on the heterogeneous BAV fusion patterns. The objective of this study was to compare computationally the regional and temporal fluid wall shear stress (WSS) and structural deformation characteristics in tricuspid aortic valve (TAV), type-0, and type-I BAV leaflets. Arbitrary Lagrangian-Eulerian fluid-structure interaction models were designed to simulate the flow and leaflet dynamics in idealized TAV, type-0, and type-I BAV geometries subjected to physiologic transvalvular pressure. The regional leaflet mechanics was quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI), temporal shear gradient (TSG), and stretch. The simulations identified regions of WSS overloads and increased WSS bidirectionality (174% increase in temporal shear magnitude, 0.10 increase in OSI on type-0 leaflets) in BAV leaflets relative to TAV leaflets. BAV leaflets also experienced larger radial deformations than TAV leaflets (4% increase in type-0 BAV leaflets). Type-I BAV leaflets exhibited contrasted WSS environments marked by WSS overloads on the non-coronary leaflet and sub-physiologic WSS levels on the fused leaflet. This study provides important insights into the mechanical characteristics of BAV leaflets, which may further our understanding of the role played by hemodynamic forces in BAV disease. Copyright © 2016 John Wiley & Sons, Ltd.

Biaxial stress relaxation of semilunar heart valve leaflets during simulated collagen catabolism: Effects of collagenase concentration and equibiaxial strain state

Heart valve leaflet collagen turnover and remodeling are innate to physiological homeostasis; valvular interstitial cells routinely catabolize damaged collagen and affect repair. Moreover, evidence indicates that leaflets can adapt to altered physiological (e.g. pregnancy) and pathological (e.g. hypertension) mechanical load states, tuning collagen structure and composition to changes in pressure and flow. However, while valvular interstitial cell-secreted matrix metalloproteinases are considered the primary effectors of collagen catabolism, the mechanisms by which damaged collagen fibers are selectively degraded remain unclear. Growing evidence suggests that the collagen fiber strain state plays a key role, with the strain-dependent configuration of the collagen molecules either masking or presenting proteolytic sites, thereby protecting or accelerating collagen proteolysis. In this study, the effects of equibiaxial strain state on collagen catabolism were investigated in porcine aortic valve and pulmonary valve tissues. Bacterial collagenase (0.2 and 0.5 mg/mL) was utilized to simulate endogenous matrix metalloproteinases, and biaxial stress relaxation and biochemical collagen concentration served as functional and compositional measures of collagen catabolism, respectively. At a collagenase concentration of 0.5 mg/mL, increasing the equibiaxial strain imposed during stress relaxation (0%, 37.5%, and 50%) yielded significantly lower median collagen concentrations in the aortic valve (p = 0.0231) and pulmonary valve (p = 0.0183), suggesting that relatively large strain magnitudes may enhance collagen catabolism. Collagen concentration decreases were paralleled by trends of accelerated normalized stress relaxation rate with equibiaxial strain in aortic valve tissues. Collectively, these in vitro results indicate that biaxial strain state is capable of affecting the susceptibility of valvular collagens to catabolism, providing a basis for further investigation of how such phenomena may manifest at different strain magnitudes or in vivo.

In-vivo heterogeneous functional and residual strains in human aortic valve leaflets

Residual and physiological functional strains in soft tissues are known to play an important role in modulating organ stress distributions. Yet, no known comprehensive information on residual strains exist, or non-invasive techniques to quantify in-vivo deformations for the aortic valve (AV) leaflets. Herein we present a completely non-invasive approach for determining heterogeneous strains - both functional and residual - in semilunar valves and apply it to normal human AV leaflets. Transesophageal 3D echocardiographic (3DE) images of the AV were acquired from open-heart transplant patients, with each AV leaflet excised after heart explant and then imaged in a flattened configuration ex-vivo. Using an established spline parameterization of both 3DE segmentations and digitized ex-vivo images (Aggarwal et al., 2014), surface strains were calculated for deformation between the ex-vivo and three in-vivo configurations: fully open, just-coapted, and fully-loaded. Results indicated that leaflet area increased by an average of 20% from the ex-vivo to in-vivo open states, with a highly heterogeneous strain field. The increase in area from open to just-coapted state was the highest at an average of 25%, while that from just-coapted to fully-loaded remained almost unaltered. Going from the ex-vivo to in-vivo mid-systole configurations, the leaflet area near the basal attachment shrank slightly, whereas the free edge expanded by ~10%. This was accompanied by a 10° -20° shear along the circumferential-radial direction. Moreover, the principal stretches aligned approximately with the circumferential and radial directions for all cases, with the highest stretch being along the radial direction. Collectively, these results indicated that even though the AV did not support any measurable pressure gradient in the just-coapted state, the leaflets were significantly pre-strained with respect to the excised state. Furthermore, the collagen fibers of the leaflet were almost fully recruited in the just-coapted state, making the leaflet very stiff with marginal deformation under full pressure. Lastly, the deformation was always higher in the radial direction and lower along the circumferential one, the latter direction made stiffer by the preferential alignment of collagen fibers. These results provide significant insight into the distribution of residual strains and the in-vivo strains encountered during valve opening and closing in AV leaflets, and will form an important component of the tool that can evaluate valve׳s functional properties in a non-invasive manner.

Biomechanical Behavior of Bioprosthetic Heart Valve Heterograft Tissues: Characterization, Simulation, and Performance

The use of replacement heart valves continues to grow due to the increased prevalence of valvular heart disease resulting from an ageing population. Since bioprosthetic heart valves (BHVs) continue to be the preferred replacement valve, there continues to be a strong need to develop better and more reliable BHVs through and improved the general understanding of BHV failure mechanisms. The major technological hurdle for the lifespan of the BHV implant continues to be the durability of the constituent leaflet biomaterials, which if improved can lead to substantial clinical impact. In order to develop improved solutions for BHV biomaterials, it is critical to have a better understanding of the inherent biomechanical behaviors of the leaflet biomaterials, including chemical treatment technologies, the impact of repetitive mechanical loading, and the inherent failure modes. This review seeks to provide a comprehensive overview of these issues, with a focus on developing insight on the mechanisms of BHV function and failure. Additionally, this review provides a detailed summary of the computational biomechanical simulations that have been used to inform and develop a higher level of understanding of BHV tissues and their failure modes. Collectively, this information should serve as a tool not only to infer reliable and dependable prosthesis function, but also to instigate and facilitate the design of future bioprosthetic valves and clinically impact cardiology.

On the effects of leaflet microstructure and constitutive model on the closing behavior of the mitral valve

Recent long-term studies showed an unsatisfactory recurrence rate of severe mitral regurgitation 3-5 years after surgical repair, suggesting that excessive tissue stresses and the resulting strain-induced tissue failure are potential etiological factors controlling the success of surgical repair for treating mitral valve (MV) diseases. We hypothesized that restoring normal MV tissue stresses in MV repair techniques would ultimately lead to improved repair durability through the restoration of MV normal homeostatic state. Therefore, we developed a micro- and macro- anatomically accurate MV finite element model by incorporating actual fiber microstructural architecture and a realistic structure-based constitutive model. We investigated MV closing behaviors, with extensive in vitro data used for validating the proposed model. Comparative and parametric studies were conducted to identify essential model fidelity and information for achieving desirable accuracy. More importantly, for the first time, the interrelationship between the local fiber ensemble behavior and the organ-level MV closing behavior was investigated using a computational simulation. These novel results indicated not only the appropriate parameter ranges, but also the importance of the microstructural tuning (i.e., straightening and re-orientation) of the collagen/elastin fiber networks at the macroscopic tissue level for facilitating the proper coaptation and natural functioning of the MV apparatus under physiological loading at the organ level. The proposed computational model would serve as a logical first step toward our long-term modeling goal-facilitating simulation-guided design of optimal surgical repair strategies for treating diseased MVs with significantly enhanced durability.

An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves

In this paper, we develop a geometrically flexible technique for computational fluid-structure interaction (FSI). The motivating application is the simulation of tri-leaflet bioprosthetic heart valve function over the complete cardiac cycle. Due to the complex motion of the heart valve leaflets, the fluid domain undergoes large deformations, including changes of topology. The proposed method directly analyzes a spline-based surface representation of the structure by immersing it into a non-boundary-fitted discretization of the surrounding fluid domain. This places our method within an emerging class of computational techniques that aim to capture geometry on non-boundary-fitted analysis meshes. We introduce the term "immersogeometric analysis" to identify this paradigm. The framework starts with an augmented Lagrangian formulation for FSI that enforces kinematic constraints with a combination of Lagrange multipliers and penalty forces. For immersed volumetric objects, we formally eliminate the multiplier field by substituting a fluid-structure interface traction, arriving at Nitsche's method for enforcing Dirichlet boundary conditions on object surfaces. For immersed thin shell structures modeled geometrically as surfaces, the tractions from opposite sides cancel due to the continuity of the background fluid solution space, leaving a penalty method. Application to a bioprosthetic heart valve, where there is a large pressure jump across the leaflets, reveals shortcomings of the penalty approach. To counteract steep pressure gradients through the structure without the conditioning problems that accompany strong penalty forces, we resurrect the Lagrange multiplier field. Further, since the fluid discretization is not tailored to the structure geometry, there is a significant error in the approximation of pressure discontinuities across the shell. This error becomes especially troublesome in residual-based stabilized methods for incompressible flow, leading to problematic compressibility at practical levels of refinement. We modify existing stabilized methods to improve performance. To evaluate the accuracy of the proposed methods, we test them on benchmark problems and compare the results with those of established boundary-fitted techniques. Finally, we simulate the coupling of the bioprosthetic heart valve and the surrounding blood flow under physiological conditions, demonstrating the effectiveness of the proposed techniques in practical computations.

The stretch responsive microRNA miR-148a-3p is a novel repressor of IKBKB, NF-κB signaling, and inflammatory gene expression in human aortic valve cells

Bicuspid aortic valves calcify at a significantly higher rate than normal aortic valves, a process that involves increased inflammation. Because we have previously found that bicuspid aortic valve experience greater stretch, we investigated the potential connection between stretch and inflammation in human aortic valve interstitial cells (AVICs). Microarray, quantitative PCR (qPCR), and protein assays performed on AVICs exposed to cyclic stretch showed that stretch was sufficient to increase expression of interleukin and metalloproteinase family members by more than 1.5-fold. Conditioned medium from stretched AVICs was sufficient to activate leukocytes. microRNA sequencing and qPCR experiments demonstrated that miR-148a-3p was repressed in both stretched AVICs (43% repression) and, as a clinical correlate, human bicuspid aortic valves (63% reduction). miR-148a-3p was found to be a novel repressor of IKBKB based on data from qPCR, luciferase, and Western blot experiments. Furthermore, increasing miR-148a-3p levels in AVICs was sufficient to decrease NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) signaling and NF-κB target gene expression. Our data demonstrate that stretch-mediated activation of inflammatory pathways is at least partly the result of stretch-repression of miR-148a-3p and a consequent failure to repress IKBKB. To our knowledge, we are the first to report that cyclic stretch of human AVICs activates inflammatory genes in a tissue-autonomous manner via a microRNA that regulates a central inflammatory pathway.