A noninvasive method for the determination of in vivo mitral valve leaflet strains

A method for non-invasive estimation of mitral valve annular regional strains

INTRODUCTION

In this study, a method to assess local deformations along the mitral annulus curve is proposed.

METHODS

It employs the known global geometry (from echocardiography) of the annulus during the cardiac cycle which is approximated with a closed cubic spline to generate a smooth mathematical representation of the annulus at each available time point. A point-wise mapping of the annular geometries between two consecutive times is established by minimizing the global displacements of the annulus. Hence, the displacements of the mitral annulus are determined and used to calculate the regional strains along the annulus line.

RESULTS

Data obtained from sonomicrometric markers are used to test the method detailed herein. The results show that our method can predict annular displacements throughout the cardiac cycle. Strain values computed with this approach are in line with previously experimentally measured strains reported in the literature. Finally, the method is applied and illustrated on an echocardiographic recording of a healthy individual.

CONCLUSION

The numerical method provided can be used to capture the regional annular strains by echocardiography and may help to predict regional dysfunctions in the mitral annulus, providing information on the pathological mechanisms.

FEBio FINESSE: An Open-Source Finite Element Simulation Approach to Estimate In Vivo Heart Valve Strains Using Shape Enforcement

PURPOSE

Finite element simulations are an enticing tool to evaluate heart valve function; however, patient-specific simulations derived from 3D echocardiography are hampered by several technical challenges. The objective of this work is to develop an open-source method to enforce matching between finite element simulations and in vivo image-derived heart valve geometry in the absence of patient-specific material properties, leaflet thickness, and chordae tendineae structures.

METHODS

We evaluate FEBio Finite Element Simulations with Shape Enforcement (FINESSE) using three synthetic test cases considering a range of model complexity. FINESSE is then used to estimate the in vivo valve behavior and leaflet strains for three pediatric patients.

RESULTS

Our results suggest that FINESSE can be used to enforce finite element simulations to match an image-derived surface and estimate the first principal leaflet strains within ± 0.03 strain. Key considerations include: (i) defining the user-defined penalty, (ii) omitting the leaflet commissures to improve simulation convergence, and (iii) emulating the chordae tendineae behavior via prescribed leaflet free edge motion or a chordae emulating force. In all patient-specific cases, FINESSE matched the target surface with median errors of approximately the smallest voxel dimension. Further analysis revealed valve-specific findings, such as the tricuspid valve leaflet strains of a 2-day old patient with HLHS being larger than those of two 13-year old patients.

CONCLUSIONS

FEBio FINESSE can be used to estimate patient-specific in vivo heart valve leaflet strains. The development of this open-source pipeline will enable future studies to begin linking in vivo leaflet mechanics with patient outcomes.

In Vivo Deformation of the Human Basilar Artery

An estimated 6.8 million people in the United States have an unruptured intracranial aneurysms, with approximately 30,000 people suffering from intracranial aneurysms rupture each year. Despite the development of population-based scores to evaluate the risk of rupture, retrospective analyses have suggested the limited usage of these scores in guiding clinical decision-making. With recent advancements in imaging technologies, artery wall motion has emerged as a promising biomarker for the general study of neurovascular mechanics and in assessing the risk of intracranial aneurysms. However, measuring arterial wall deformations in vivo itself poses several challenges, including how to image local wall motion and deriving the anisotropic wall strains over the cardiac cycle. To overcome these difficulties, we first developed a novel in vivo MRI-based imaging method to acquire cardiac gated images of the human basilar artery (BA) over the cardiac cycle. Next, complete BA endoluminal surfaces from each frame were segmented, producing high-resolution point clouds of the endoluminal surfaces. From these point clouds we developed a novel B-spline-based surface representation, then exploited the local support nature of B-splines to determine the local endoluminal surface strains. Results indicated distinct regional and temporal variations in BA wall deformation, highlighting the heterogeneous nature BA function. These included large circumferential strains (up to ∼ 20 % ), and small longitudinal strains, which were often contractile and out of phase with the circumferential strains patterns. Of particular interest was the temporal phase lag in the maximum circumferential perimeter length, which indicated that the BA deforms asynchronously over the cardiac cycle. In summary, the proposed method enabled local deformation analysis, allowing for the successful reproduction of local features of the BA, such as regional principal stretches, areal changes, and pulsatile motion. Integrating the proposed method into existing population-based scores has the potential to improve our understanding of mechanical properties of human BA and enhance clinical decision-making.

Simulated Effects of Acute Left Ventricular Myocardial Infarction on Mitral Regurgitation in an Ovine Model

Ischemic mitral regurgitation (IMR) occurs from incomplete coaptation of the mitral valve (MV) after myocardial infarction (MI), typically worsened by continued remodeling of the left ventricular (LV). The importance of LV remodeling is clear as IMR is induced by the post-MI dual mechanisms of mitral annular dilation and leaflet tethering from papillary muscle (PM) distension via the MV chordae tendineae (MVCT). However, the detailed etiology of IMR remains poorly understood, in large part due to the complex interactions of the MV and the post-MI LV remodeling processes. Given the patient-specific anatomical complexities of the IMR disease processes, simulation-based approaches represent an ideal approach to improve our understanding of this deadly disease. However, development of patient-specific models of left ventricle-mitral valve (LV-MV) interactions in IMR are complicated by the substantial variability and complexity of the MR etiology itself, making it difficult to extract underlying mechanisms from clinical data alone. To address these shortcomings, we developed a detailed ovine LV-MV finite element (FE) model based on extant comprehensive ovine experimental data. First, an extant ovine LV FE model (Sci. Rep. 2021 Jun 29;11(1):13466) was extended to incorporate the MV using a high fidelity ovine in vivo derived MV leaflet geometry. As it is not currently possible to image the MVCT in vivo, a functionally equivalent MVCT network was developed to create the final LV-MV model. Interestingly, in pilot studies, the MV leaflet strains did not agree well with known in vivo MV leaflet strain fields. We then incorporated previously reported MV leaflet prestrains (J. Biomech. Eng. 2023 Nov 1;145(11):111002) in the simulations. The resulting LV-MV model produced excellent agreement with the known in vivo ovine MV leaflet strains and deformed shapes in the normal state. We then simulated the effects of regional acute infarctions of varying sizes and anatomical locations by shutting down the local myocardial contractility. The remaining healthy (noninfarcted) myocardium mechanical behaviors were maintained, but allowed to adjust their active contractile patterns to maintain the prescribed pressure-volume loop behaviors in the acute post-MI state. For all cases studied, the LV-MV simulation demonstrated excellent agreement with known LV and MV in vivo strains and MV regurgitation orifice areas. Infarct location was shown to play a critical role in resultant MV leaflet strain fields. Specifically, extensional deformations of the posterior leaflets occurred in the posterobasal and laterobasal infarcts, while compressive deformations of the anterior leaflet were observed in the anterobasal infarct. Moreover, the simulated posterobasal infarct induced the largest MV regurgitation orifice area, consistent with experimental observations. The present study is the first detailed LV-MV simulation that reveals the important role of MV leaflet prestrain and functionally equivalent MVCT for accurate predictions of LV-MV interactions. Importantly, the current study further underscored simulation-based methods in understanding MV function as an integral part of the LV.

Computational simulation model of transcatheter edge-to-edge mitral valve repair: a proof-of-concept study

AIMS

As transcatheter mitral valve (MV) interventions are expanding and more device types and sizes become available, a tool supporting operators in pre-procedural planning and the clinical decision-making process is highly desirable. We sought to develop a finite element computational simulation model to predict the results of transcatheter edge-to-edge repair (TEER) interventions.

METHODS AND RESULTS

We prospectively enrolled patients with secondary mitral regurgitation (MR) referred for a clinically indicated TEER. The 3D trans-oesophageal echocardiograms performed at the beginning of the procedure were used to perform the simulation. On the 3D dynamic model of the MV that was first obtained, we simulated the clip implantation using the same clip type, size, number, and implantation location that was used during the intervention. The 3D model of the MV obtained after the simulation of the clip implantation was compared with the clinical results obtained at the end of the intervention. We analysed the degree and location of residual MR and the shape and area of the diastolic MV area. We performed computational simulation on five patients. Overall, the simulated models predicted well the degree and location of the residual regurgitant orifice(s) but tended to underestimate the diastolic mitral orifice area.

CONCLUSION

In this proof-of-concept study, we present preliminary results on our algorithm simulating clip implantation in five patients with functional MR. We show promising results regarding the feasibility and accuracy in terms of predicting residual MR and the need to improve the estimation of the diastolic MV area.

Bayesian Optimization-Based Inverse Finite Element Analysis for Atrioventricular Heart Valves

Inverse finite element analysis (iFEA) of the atrioventricular heart valves (AHVs) can provide insights into the in-vivo valvular function, such as in-vivo tissue strains; however, there are several limitations in the current state-of-the-art that iFEA has not been widely employed to predict the in-vivo, patient-specific AHV leaflet mechanical responses. In this exploratory study, we propose the use of Bayesian optimization (BO) to study the AHV functional behaviors in-vivo. We analyzed the efficacy of Bayesian optimization to estimate the isotropic Lee-Sacks material coefficients in three benchmark problems: (i) an inflation test, (ii) a simplified leaflet contact model, and (iii) an idealized AHV model. Then, we applied the developed BO-iFEA framework to predict the leaflet properties for a patient-specific tricuspid valve under a congenital heart defect condition. We found that the BO could accurately construct the objective function surface compared to the one from a [Formula: see text] grid search analysis. Additionally, in all cases the proposed BO-iFEA framework yielded material parameter predictions with average element errors less than 0.02 mm/mm (normalized by the simulation-specific characteristic length). Nonetheless, the solutions were not unique due to the presence of a long-valley minima region in the objective function surfaces. Parameter sets along this valley can yield functionally equivalent outcomes (i.e., closing behavior) and are typically observed in the inverse analysis or parameter estimation for the nonlinear mechanical responses of the AHV. In this study, our key contributions include: (i) a first-of-its-kind demonstration of the BO method used for the AHV iFEA; and (ii) the evaluation of a candidate AHV in-silico modeling approach wherein the chordae could be substituted with equivalent displacement boundary conditions, rendering the better iFEA convergence and a smoother objective surface.

Patient-Specific Quantitative In-Vivo Assessment of Human Mitral Valve Leaflet Strain Before and After MitraClip Repair

PURPOSE

Mitral regurgitation (MR) is a highly prevalent and deadly cardiac disease characterized by improper mitral valve (MV) leaflet coaptation. Among the plethora of available treatment strategies, the MitraClip is an especially safe option, but optimizing its long-term efficacy remains an urgent challenge.

METHODS

We applied our noninvasive image-based strain computation pipeline [1] to intraoperative transesophageal echocardiography datasets taken from ten patients undergoing MitraClip repair, spanning a range of MR etiologies and MitraClip configurations. We then analyzed MV leaflet strains before and after MitraClip implementation to develop a better understanding of (1) the pre-operative state of human regurgitant MV, and (2) the MitraClip's impact on the MV leaflet deformations.

RESULTS

The MV pre-operative strain fields were highly variable, underscoring both the heterogeneity of the MR in the patient population and the need for patient-specific treatment approaches. Similarly, there were no consistent overall post-operative strain patterns, although the average A2 segment radial strain difference between pre- and post-operative states was consistently positive. In contrast, the post-operative strain fields were better correlated to their respective pre-operative strain fields than to the inter-patient post-operative strain fields. This quantitative result implies that the patient specific pre-operative state of the MV guides its post-operative deformation, which suggests that the post-operative state can be predicted using pre-operative data-derived modelling alone.

CONCLUSIONS

The pre-operative MV leaflet strain patterns varied considerably across the range of MR disease states and after MitraClip repair. Despite large inter-patient heterogeneity, the post-operative deformation appears principally dictated by the pre-operative deformation state. This novel finding suggests that though the variation in MR functional state and MitraClip-induced deformation were substantial, the post-operative state can be predicted from the pre-operative data alone. This study suggests that, with use of larger patient cohort and corresponding long-term outcomes, quantitative predictive factors of MitraClip durability can be identified.

Virtual blebbistatin: A robust and rapid software approach to motion artifact removal in optical mapping of cardiomyocytes

Fluorescent reporters of cardiac electrophysiology provide valuable information on heart cell and tissue function. However, motion artifacts caused by cardiac muscle contraction interfere with accurate measurement of fluorescence signals. Although drugs such as blebbistatin can be applied to stop cardiac tissue from contracting by uncoupling calcium-contraction, their usage prevents the study of excitation-contraction coupling and, as we show, impacts cellular structure. We therefore developed a robust method to remove motion computationally from images of contracting cardiac muscle and to map fluorescent reporters of cardiac electrophysiological activity onto images of undeformed tissue. When validated on cardiomyocytes derived from human induced pluripotent stem cells (iPSCs), in both monolayers and engineered tissues, the method enabled efficient and robust reduction of motion artifact. As with pharmacologic approaches using blebbistatin for motion removal, our algorithm improved the accuracy of optical mapping, as demonstrated by spatial maps of calcium transient decay. However, unlike pharmacologic motion removal, our computational approach allowed direct analysis of calcium-contraction coupling. Results revealed calcium-contraction coupling to be more uniform across cells within engineered tissues than across cells in monolayer culture. The algorithm shows promise as a robust and accurate tool for optical mapping studies of excitation-contraction coupling in heart tissue.

Mechanical Behavior and Collagen Structure of Degenerative Mitral Valve Leaflets and a Finite Element Model of Primary Mitral Regurgitation

Degenerative mitral valve disease is the main cause of primary mitral regurgitation with two phenotypes: fibroelastic deficiency (FED) often with localized myxomatous degeneration and diffuse myxomatous degeneration or Barlow's disease. Myxomatous degeneration disrupts the microstructure of the mitral valve leaflets, particularly the collagen fibers, which affects the mechanical behavior of the leaflets. The present study uses biaxial mechanical tests and second harmonic generation microscopy to examine the mechanical behavior of Barlow and FED tissue. Three tissue samples were harvested from a FED patient and one sample is from a Barlow patient. Then we use an appropriate constitutive model by excluding the collagen fibers under compression. Finally, we built an FE model based on the echocardiography of patients diagnosed with FED and Barlow and the characterized material model and collagen fiber orientation. The Barlow sample and the FED sample from the most affected segment showed different mechanical behavior and collagen structure compared to the other two FED samples. The FE model showed very good agreement with echocardiography with 2.02±1.8 mm and 1.05±0.79 mm point-to-mesh distance errors for Barlow and FED patients, respectively. It has also been shown that the exclusion of collagen fibers under compression provides versatility for the material model; it behaves stiff in the belly region, preventing excessive bulging, while it behaves very softly in the commissures to facilitate folding. STATEMENT OF SIGNIFICANCE: None.

Simulation of Mitral Valve Plasticity in Response to Myocardial Infarction

Left ventricular myocardial infarction (MI) has broad and debilitating effects on cardiac function. In many cases, MI leads to ischemic mitral regurgitation (IMR), a condition characterized by incompetency of the mitral valve (MV). IMR has many deleterious effects as well as a high mortality rate. While various clinical treatments for IMR exist, success of these procedures remains limited, in large part because IMR dramatically alters the geometry and function of the MV in ways that are currently not well understood. Previous investigations of post-MI MV remodeling have elucidated that MV tissues have a significant ability to undergo a form of permanent inelastic deformations in the first phase of the post-MI period. These changes appear to be attributable to the altered loading and boundary conditions on the MV itself, as opposed to an independent pathophysiological process. Mechanistically, these results suggest that the MV mostly responds passively to MI during the first 8 weeks post-MI by undergoing a permanent deformation. In the present study, we developed the first computational model of this post-MI MV remodeling process, which we term "mitral valve plasticity." Integrating methodologies and insights from previous studies of in vivo ovine MV function, image-based patient-specific model development, and post-MI MV adaptation, we constructed a representative geometric model of a pre-MI MV. We then performed finite element simulations of the entire MV apparatus under time-dependent boundary conditions and accounting for changes to material properties equivalent to those observed 0-8 weeks post-MI. Our results suggest that during this initial period of adaptation, the MV response to MI can be accurately modeled using a soft tissue plasticity approach, similar to permanent set frameworks that have been applied previously in the context of exogenously crosslinked tissues.

In vivo assessment of mitral valve leaflet remodelling following myocardial infarction

Each year, more than 40,000 people undergo mitral valve (MV) repair surgery domestically to treat regurgitation caused by myocardial infarction (MI). Although continual MV tissue remodelling following repair is believed to be a major contributor to regurgitation recurrence, the effects of the post-MI state on MV remodelling remain poorly understood. This lack of understanding limits our ability to predict the remodelling of the MV both post-MI and post-surgery to facilitate surgical planning. As a necessary first step, the present study was undertaken to noninvasively quantify the effects of MI on MV remodelling in terms of leaflet geometry and deformation. MI was induced in eight adult Dorset sheep, and real-time three-dimensional echocardiographic (rt-3DE) scans were collected pre-MI as well as at 0, 4, and 8 weeks post-MI. A previously validated image-based morphing pipeline was used to register corresponding open- and closed-state scans and extract local in-plane strains throughout the leaflet surface at systole. We determined that MI induced permanent changes in leaflet dimensions in the diastolic configuration, which increased with time to 4 weeks, then stabilised. MI substantially affected the systolic shape of the MV, and the range of stretch experienced by the MV leaflet at peak systole was substantially reduced when referred to the current time-point. Interestingly, when we referred the leaflet strains to the pre-MI configuration, the systolic strains remained very similar throughout the post-MI period. Overall, we observed that post-MI ventricular remodeling induced permanent changes in the MV leaflet shape. This predominantly affected the MV's diastolic configuration, leading in turn to a significant decrease in the range of stretch experienced by the leaflet when referenced to the current diastolic configuration. These findings are consistent with our previous work that demonstrated increased plastic (i.e. non-recoverable) leaflet deformations post-MI, that was completely accounted for by the associated changes in collagen fiber structure. Moreover, we demonstrated through noninvasive methods that the state of the MV leaflet can elucidate the progression and extent of MV adaptation following MI and is thus highly relevant to the design of current and novel patient specific minimally invasive surgical repair strategies.

Effect of transcatheter edge-to-edge repair device position on diastolic hemodynamic parameters: An echocardiography-based simulation study

Background

Transcatheter edge-to-edge repair (TEER) has developed from innovative technology to an established treatment strategy of mitral regurgitation (MR). The risk of iatrogenic mitral stenosis after TEER is, however, a critical factor in the conflict of interest between maximal reduction of MR and minimal impairment of left ventricular filling. We aim to investigate systematically the impact of device position on the post treatment hemodynamic outcome by involving the patient-specific segmentation of the diseased mitral valve.

Materials and methods

Transesophageal echocardiographic image data of ten patients with severe MR (age: 57 ± 8 years, 20% female) were segmented and virtually treated with TEER at three positions by using a position based dynamics approach. Pre- and post-interventional patient geometries were preprocessed for computational fluid dynamics (CFD) and simulated at peak-diastole with patient-specific blood flow boundary conditions. Simulations were performed with boundary conditions mimicking rest and stress. The simulation results were compared with clinical data acquired for a cohort of 21 symptomatic MR patients (age: 79 ± 6 years, 43% female) treated with TEER.

Results

Virtual TEER reduces the mitral valve area (MVA) from 7.5 ± 1.6 to 2.6 ± 0.6 cm2. Central device positioning resulted in a 14% smaller MVA than eccentric device positions. Furthermore, residual MVA is better predictable for central than for eccentric device positions (R 2 = 0.81 vs. R 2 = 0.49). The MVA reduction led to significantly higher maximal diastolic velocities (pre: 0.9 ± 0.2 m/s, post: 2.0 ± 0.5 m/s) and pressure gradients (pre: 1.5 ± 0.6 mmHg, post: 16.3 ± 9 mmHg) in spite of a mean flow rate reduction by 23% due to reduced MR after the treatment. On average, velocities were 12% and pressure gradients were 25% higher with devices in central compared to lateral or medial positions.

Conclusion

Virtual TEER treatment combined with CFD is a promising tool for predicting individual morphometric and hemodynamic outcomes. Such a tool can potentially be used to support clinical decision making, procedure planning, and risk estimation to prevent post-procedural iatrogenic mitral stenosis.

Mitral Annular Elasticity Determines Severity of Regurgitation in Barlow's Disease

AIMS

Barlow's disease (BD) with late systolic mitral regurgitation provides diagnostic and therapeutic challenges. The mechanisms of the regurgitation are still unclear. We hypothesized that the onset and the severity of late systolic regurgitation are determined by annulus dynamics and the mechanical stresses imposed by the left ventricle.

METHODS AND RESULTS

Ten patients with BD and mitral annulus disjunction (MAD) were compared with ten healthy controls. Resting blood pressure (BP) was measured and transthoracic three-dimensional (3D) echocardiography was analyzed using a holographic display that allows tracking and measurements of mitral annulus surface area (ASA) throughout the cardiac cycle. A novel annulus elastance index (dASA/dP) was calculated between aortic valve opening and onset of mitral regurgitation. Severity of MAD was quantified as the disjunction index (mm*degree). Leaflet coaptation area was calculated using a finite element model. Peak systolic ASA in controls and patients were 9.3±0.6 and 21.1±3.1 cm², respectively (p<0.001). In patients ASA increased rapidly during LV ejection and onset of mitral regurgitation coincided closely with peak upslope of annulus area change (dASA/dt). The finite element model showed a close association between rapid annulus displacement and coaptation area deficit in BD. Systolic annulus elastance index (0.058±0.036 cm²/mmHg) correlated strongly with disjunction index (r=0.91, p<0.0001). Moreover, regurgitation volume showed a positive correlation with systolic BP (r=0.80, p<0.01) CONCLUSION: The present pilot study supports the hypothesis that annulus dilatation may accentuate mitral valve regurgitation in patients with Barlow's disease. A novel annulus elastance index may predict the severity of mitral valve regurgitation in selected patients.

Patient-Specific Three-Dimensional Ultrasound Derived Computational Modeling of the Mitral Valve

Several new techniques to repair the mitral valve affected by functional mitral regurgitation are in development. However, due to the heterogeneity of valve lesions between patients, predicting the outcomes of novel treatment approaches is challenging. We present a patient-specific, 3D ultrasound-derived computational model of the mitral valve for procedure planning, that faithfully mimics the pathological valve dynamics. 3D ultrasound images were obtained in three pigs induced with heart failure and which developed functional mitral regurgitation. For each case, images were segmented, and finite element model of mitral valve was constructed. Annular and papillary muscle dynamics were extracted and imposed as kinematic boundary conditions, and the chordae were pre-strained to induce valve tethering. Valve closure was simulated by applying physiologic transvalvular pressure on the leaflets. Agreement between simulation results and truth datasets was confirmed, with accurate location of regurgitation jets and coaptation defects. Inclusion of kinematic patient-specific boundary conditions was necessary to achieve these results, whereas use of idealized boundary conditions deviated from the truth dataset. Due to the impact of boundary conditions on the model, the effect of repair strategies on valve closure varied as well, indicating that our approach of using patient-specific boundary conditions for mitral valve modeling is valid.

Clinical Impact of Computational Heart Valve Models

This paper provides a review of engineering applications and computational methods used to analyze the dynamics of heart valve closures in healthy and diseased states. Computational methods are a cost-effective tool that can be used to evaluate the flow parameters of heart valves. Valve repair and replacement have long-term stability and biocompatibility issues, highlighting the need for a more robust method for resolving valvular disease. For example, while fluid-structure interaction analyses are still scarcely utilized to study aortic valves, computational fluid dynamics is used to assess the effect of different aortic valve morphologies on velocity profiles, flow patterns, helicity, wall shear stress, and oscillatory shear index in the thoracic aorta. It has been analyzed that computational flow dynamic analyses can be integrated with other methods to create a superior, more compatible method of understanding risk and compatibility.

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.

Pre-surgical Prediction of Ischemic Mitral Regurgitation Recurrence Using In Vivo Mitral Valve Leaflet Strains

Ischemic mitral regurgitation (IMR) is a prevalent cardiac disease associated with substantial morbidity and mortality. Contemporary surgical treatments continue to have limited long-term success, in part due to the complex and multi-factorial nature of IMR. There is thus a need to better understand IMR etiology to guide optimal patient specific treatments. Herein, we applied our finite element-based shape-matching technique to non-invasively estimate peak systolic leaflet strains in human mitral valves (MVs) from in-vivo 3D echocardiographic images taken immediately prior to and post-annuloplasty repair. From a total of 21 MVs, we found statistically significant differences in pre-surgical MV size, shape, and deformation patterns between the with and without IMR recurrence patient groups at 6 months post-surgery. Recurrent MVs had significantly less compressive circumferential strains in the anterior commissure region compared to the recurrent MVs (p = 0.0223) and were significantly larger. A logistic regression analysis revealed that average pre-surgical circumferential leaflet strain in the Carpentier A1 region independently predicted 6-month recurrence of IMR (optimal cutoff value - 18%, p = 0.0362). Collectively, these results suggest greater disease progression in the recurrent group and underscore the highly patient-specific nature of IMR. Importantly, the ability to identify such factors pre-surgically could be used to guide optimal treatment methods to reduce post-surgical IMR recurrence.

Uncertainty quantification in subject-specific estimation of local vessel mechanical properties

Quantitative estimation of local mechanical properties remains critically important in the ongoing effort to elucidate how blood vessels establish, maintain, or lose mechanical homeostasis. Recent advances based on panoramic digital image correlation (pDIC) have made high-fidelity 3D reconstructions of small-animal (e.g., murine) vessels possible when imaged in a variety of quasi-statically loaded configurations. While we have previously developed and validated inverse modeling approaches to translate pDIC-measured surface deformations into biomechanical metrics of interest, our workflow did not heretofore include a methodology to quantify uncertainties associated with local point estimates of mechanical properties. This limitation has compromised our ability to infer biomechanical properties on a subject-specific basis, such as whether stiffness differs significantly between multiple material locations on the same vessel or whether stiffness differs significantly between multiple vessels at a corresponding material location. In the present study, we have integrated a novel uncertainty quantification and propagation pipeline within our inverse modeling approach, relying on empirical and analytic Bayesian techniques. To demonstrate the approach, we present illustrative results for the ascending thoracic aorta from three mouse models, quantifying uncertainties in constitutive model parameters as well as circumferential and axial tangent stiffness. Our extended workflow not only allows parameter uncertainties to be systematically reported, but also facilitates both subject-specific and group-level statistical analyses of the mechanics of the vessel wall.

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.

Scaling digital twins from the artisanal to the industrial

Mathematical modeling and simulation are moving from being powerful development and analysis tools towards having increased roles in operational monitoring, control and decision support, in which models of specific entities are continually updated in the form of a digital twin. However, current digital twins are largely the result of bespoke technical solutions that are difficult to scale. We discuss two exemplar applications that motivate challenges and opportunities for scaling digital twins, and that underscore potential barriers to wider adoption of this technology.

In Vivo Image-Based 4D Modeling of Competent and Regurgitant Mitral Valve Dynamics

BACKGROUND

In vivo characterization of mitral valve dynamics relies on image analysis algorithms that accurately reconstruct valve morphology and motion from clinical images. The goal of such algorithms is to provide patient-specific descriptions of both competent and regurgitant mitral valves, which can be used as input to biomechanical analyses and provide insights into the pathophysiology of diseases like ischemic mitral regurgitation (IMR).

OBJECTIVE

The goal is to generate accurate image-based representations of valve dynamics that visually and quantitatively capture normal and pathological valve function.

METHODS

We present a novel framework for 4D segmentation and geometric modeling of the mitral valve in real-time 3D echocardiography (rt-3DE), an imaging modality used for pre-operative surgical planning of mitral interventions. The framework integrates groupwise multi-atlas label fusion and template-based medial modeling with Kalman filtering to generate quantitatively descriptive and temporally consistent models of valve dynamics.

RESULTS

The algorithm is evaluated on rt-3DE data series from 28 patients: 14 with normal mitral valve morphology and 14 with severe IMR. In these 28 data series that total 613 individual 3DE images, each 3D mitral valve segmentation is validated against manual tracing, and temporal consistency between segmentations is demonstrated.

CONCLUSIONS

Automated 4D image analysis allows for reliable non-invasive modeling of the mitral valve over the cardiac cycle for comparison of annular and leaflet dynamics in pathological and normal mitral valves. Future studies can apply this algorithm to cardiovascular mechanics applications, including patient-specific strain estimation, fluid dynamics simulation, inverse finite element analysis, and risk stratification for surgical treatment.

On the role of predicted in vivo mitral valve interstitial cell deformation on its biosynthetic behavior

Ischemic mitral regurgitation (IMR), a frequent complication of myocardial infarction, is characterized by regurgitation of blood from the left ventricle back into the left atrium. Physical interventions via surgery or less-invasive techniques are the only available therapies for IMR, with valve repair via undersized ring annuloplasty (URA) generally preferred over valve replacement. However, recurrence of IMR after URA occurs frequently and is attributed to continued remodeling of the MV and infarct region of the left ventricle. The mitral valve interstitial cells (MVICs) that maintain the tissue integrity of the MV leaflets are highly mechanosensitive, and altered loading post-URA is thought to lead to aberrant MVIC-directed tissue remodeling. Although studies have investigated aspects of mechanically directed VIC activation and remodeling potential, there remains a substantial disconnect between organ-level biomechanics and cell-level phenomena. Herein, we utilized an extant multiscale computational model of the MV that linked MVIC to organ-level MV biomechanical behaviors to simulate changes in MVIC deformation following URA. A planar biaxial bioreactor system was then used to cyclically stretch explanted MV leaflet tissue, emulating the in vivo changes in loading following URA. This simulation-directed experimental investigation revealed that post-URA deformations resulted in decreased MVIC activation and collagen mass fraction. These results are consistent with the hypothesis that URA failures post-IMR are due, in part, to reduced MVIC-mediated maintenance of the MV leaflet tissue resulting from a reduction in physical stimuli required for leaflet tissue homeostasis. Such information can inform the development of novel URA strategies with improved durability.

The time has come to extend the expiration limit of cryopreserved allograft heart valves

Despite the wide choice of commercial heart valve prostheses, cryopreserved semilunar allograft heart valves (C-AHV) are required, and successfully transplanted in selected groups of patients. The expiration limit (EL) criteria have not been defined yet. Most Tissue Establishments (TE) use the EL of 5 years. From physiological, functional, and surgical point of view, the morphology and mechanical properties of aortic and pulmonary roots represent basic features limiting the EL of C-AHV. The aim of this work was to review methods of AHV tissue structural analysis and mechanical testing from the perspective of suitability for EL validation studies. Microscopic structure analysis of great arterial wall and semilunar leaflets tissue should clearly demonstrate cells as well as the extracellular matrix components by highly reproducible and specific histological staining procedures. Quantitative morphometry using stereological grids has proved to be effective, as the exact statistics was feasible. From mechanical testing methods, tensile test was the most suitable. Young's moduli of elasticity, ultimate stress and strain were shown to represent most important AHV tissue mechanical characteristics, suitable for exact statistical analysis. C-AHV are prepared by many different protocols, so as each TE has to work out own EL for C-AHV.

Mitral regurgitation after transcatheter aortic valve replacement

Patients undergoing transcatheter aortic valve replacement (TAVR) might have an associated significant MR that can potentially lead to left ventricular (LV) failure after procedure. Considering the specific alterations in the mitral valve in TAVR scenario and the widespread use of TAVR in recent years, it appears important to know and understand the anatomical, functional and clinical implications to develop adequate strategies for the future. Patients with severe mitral regurgitation (MR) have been generally excluded from randomized clinical trials, making poor the impact that associated MR can have on clinical outcomes after TAVR. Several factors must be considered whose presence influences the severity of MR. For example, the elevated prevalence of coronary disease with consequent ischemic MR may account for LV dilation observed at the end stage of aortic stenosis. Evidence randomized studies and registries suggests that the rate of concomitant moderate-to-severe MR in patients undergoing TAVR oscillates between 2% and 33%, and patients with moderate to severe MR may have hemodynamic frailty with clinical deterioration during mechanical intervention. Short- and long-term outcomes, including cardiac mortality, appear to be influenced by the existence of preoperative moderate-to-severe MR or by the postprocedural worsening of mild MR, generally due to adverse LV remodeling. The incidence and the prognostic effect of concomitant MR in patients undergoing TAVR requires specific attention as might trigger adjunctive strategy treatment which should be carefully evaluated in clinical trials.

KEYWORDS

Mitral regurgitation (MR); mitral valve; transcatheter aortic valve; transcatheter aortic valve replacement (TAVR).

Mitral valve leaflet response to ischaemic mitral regurgitation: from gene expression to tissue remodelling

Ischaemic mitral regurgitation (IMR), a frequent complication following myocardial infarction (MI), leads to higher mortality and poor clinical prognosis if untreated. Accumulating evidence suggests that mitral valve (MV) leaflets actively remodel post MI, and this remodelling increases both the severity of IMR and the occurrence of MV repair failures. However, the mechanisms of extracellular matrix maintenance and modulation by MV interstitial cells (MVICs) and their impact on MV leaflet tissue integrity and repair failure remain largely unknown. Herein, we sought to elucidate the multiscale behaviour of IMR-induced MV remodelling using an established ovine model. Leaflet tissue at eight weeks post MI exhibited significant permanent plastic radial deformation, eliminating mechanical anisotropy, accompanied by altered leaflet composition. Interestingly, no changes in effective collagen fibre modulus were observed, with MVICs slightly rounder, at eight weeks post MI. RNA sequencing indicated that YAP-induced genes were elevated at four weeks post MI, indicating elevated mechanotransduction. Genes related to extracellular matrix organization were downregulated at four weeks post MI when IMR occurred. Transcriptomic changes returned to baseline by eight weeks post MI. This multiscale study suggests that IMR induces plastic deformation of the MV with no functional damage to the collagen fibres, providing crucial information for computational simulations of the MV in IMR.

Geometric description for the anatomy of the mitral valve: A review

The mitral valve is a complex anatomical structure whose physiological functioning relies on the biomechanical properties and structural integrity of its components. Their compromise can lead to mitral valve dysfunction, associated with morbidity and mortality. Therefore, a review on the morphometry of the mitral valve is crucial, more specifically on the importance of valve dimensions and shape for its function. This review initially provides a brief background on the anatomy and physiology of the mitral valve, followed by an analysis of the morphological information available. A characterisation of mathematical descriptions of several parts of the valve is performed and the impact of different dimensions and shape changes in disease is then outlined. Finally, a section regarding future directions and recommendations for the use of morphometric information in clinical analysis of the mitral valve is presented.

Modeling the mitral valve

This work is concerned with modeling and simulation of the mitral valve, one of the four valves in the human heart. The valve is composed of leaflets, the free edges of which are supported by a system of chordae, which themselves are anchored to the papillary muscles inside the left ventricle. First, we examine valve anatomy and present the results of original dissections. These display the gross anatomy and information on fiber structure of the mitral valve. Next, we build a model valve following a design-based methodology, meaning that we derive the model geometry and the forces that are needed to support a given load and construct the model accordingly. We incorporate information from the dissections to specify the fiber topology of this model. We assume the valve achieves mechanical equilibrium while supporting a static pressure load. The solution to the resulting differential equations determines the pressurized configuration of the valve model. To complete the model, we then specify a constitutive law based on a stress-strain relation consistent with experimental data that achieves the necessary forces computed in previous steps. Finally, using the immersed boundary method, we simulate the model valve in fluid in a computer test chamber. The model opens easily and closes without leak when driven by physiological pressures over multiple beats. Further, its closure is robust to driving pressures that lack atrial systole or are much lower or higher than normal.

Analyzing valve interstitial cell mechanics and geometry with spatial statistics

Understanding cell geometric and mechanical properties is crucial to understanding how cells sense and respond to their local environment. Moreover, changes to cell mechanical properties under varied micro-environmental conditions can both influence and indicate fundamental changes to cell behavior. Atomic Force Microscopy (AFM) is a well established, powerful tool to capture geometric and mechanical properties of cells. We have previously demonstrated substantial functional and behavioral differences between aortic and pulmonary valve interstitial cells (VIC) using AFM and subsequent models of VIC mechanical response. In the present work, we extend these studies by demonstrating that to best interpret the spatially distributed AFM data, the use of spatial statistics is required. Spatial statistics includes formal techniques to analyze spatially distributed data, and has been used successfully in the analysis of geographic data. Thus, spatially mapped AFM studies of cell geometry and mechanics are analogous to more traditional forms of geospatial data. We are able to compare the spatial autocorrelation of stiffness in aortic and pulmonary valve interstitial cells, and more accurately capture cell geometry from height recordings. Specifically, we showed that pulmonary valve interstitial cells display higher levels of spatial autocorrelation of stiffness than aortic valve interstitial cells. This suggests that aortic VICs form different stress fiber structures than their pulmonary counterparts, in addition to being more highly expressed and stiffer on average. Thus, the addition of spatial statistics can contribute to our fundamental understanding of the differences between cell types. Moving forward, we anticipate that this work will be meaningful to enhance direct analysis of experimental data and for constructing high fidelity computational of VICs and other cell models.

On the simulation of mitral valve function in health, disease, and treatment

The mitral valve (MV) is the heart valve that regulates blood ?ow between the left atrium and left ventricle (LV). In situations where the MV fails to fully cover the left atrioventricular ori?ce during systole, the resulting regurgitation causes pulmonary congestion, leading to heart failure and/or stroke. The causes of MV insuf?ciency can be either primary (e.g. myxomatous degeneration) where the valvular tissue is organically diseased, or secondary (typically inducded by ischemic cardiomyopathy) termed ischemic mitral regurgitation (IMR), is brought on by adverse LV remodeling. IMR is present in up to 40% of patients and more than doubles the probability of cardiovascular morbidity after 3.5 years. There is now agreement that adjunctive procedures are required to treat IMR caused by lea?et tethering. However, there is no consensus regarding the best procedure. Multicenter registries and randomized trials would be necessary to prove which procedure is superior. Given the number of proposed procedures and the complexity and duration of such studies, it is highly unlikely that IMR procedure optimization will be achieved by prospective clinical trials. There is thus an urgent need for cell and tissue physiologically based quantitative assessments of MV function to better design surgical solutions and associated therapies. Novel computational approaches directed towards optimized surgical repair procedures can substantially reduce the need for such trial-and-error approaches. We present the details of our MV modeling techniques, with an emphasis on what is known and investigated at various length scales.

Development of a Functionally Equivalent Model of the Mitral Valve Chordae Tendineae Through Topology Optimization

Ischemic mitral regurgitation (IMR) is a currently prevalent disease in the US that is projected to become increasingly common as the aging population grows. In recent years, image-based simulations of mitral valve (MV) function have improved significantly, providing new tools to refine IMR treatment. However, clinical implementation of MV simulations has long been hindered as the in vivo MV chordae tendineae (MVCT) geometry cannot be captured with sufficient fidelity for computational modeling. In the current study, we addressed this challenge by developing a method to produce functionally equivalent MVCT models that can be built from the image-based MV leaflet geometry alone. We began our analysis using extant micron-resolution 3D imaging datasets to first build anatomically accurate MV models. We then systematically simplified the native MVCT structure to generate a series of synthetic models by consecutively removing key anatomic features, such as the thickness variations, branching patterns, and chordal origin distributions. In addition, through topology optimization, we identified the minimal structural complexity required to capture the native MVCT behavior. To assess the performance and predictive power of each synthetic model, we analyzed their performance by comparing the mismatch in simulated MV closed shape, as well as the strain and stress tensors, to ground-truth MV models. Interestingly, our results revealed a substantial redundancy in the anatomic structure of native chordal anatomy. We showed that the closing behavior of complete MV apparatus under normal, diseased, and surgically repaired scenarios can be faithfully replicated by a functionally equivalent MVCT model comprised of two representative papillary muscle heads, single strand chords, and a uniform insertion distribution with a density of 15 insertions/cm2. Hence, even though the complete sub-valvular structure is mostly missing in in vivo MV images, we believe our approach will allow for the development of patient-specific complete MV models for surgical repair planning.

The Three-Dimensional Microenvironment of the Mitral Valve: Insights into the Effects of Physiological Loads

INTRODUCTION

In the mitral valve (MV), numerous pathological factors, especially those resulting from changes in external loading, have been shown to affect MV structure and composition. Such changes are driven by the MV interstitial cell (MVIC) population via protein synthesis and enzymatic degradation of extracellular matrix (ECM) components.

METHODS

While cell phenotype, ECM composition and regulation, and tissue level changes in MVIC shape under stress have been studied, a detailed understanding of the three-dimensional (3D) microstructural mechanisms are lacking. As a first step in addressing this challenge, we applied focused ion beam scanning electron microscopy (FIB-SEM) to reveal novel details of the MV microenvironment in 3D.

RESULTS

We demonstrated that collagen is organized into large fibers consisting of an average of 605 ± 113 fibrils, with a mean diameter of 61.2 ± 9.8 nm. In contrast, elastin was organized into two distinct structural subtypes: (1) sheet-like lamellar elastin, and (2) circumferentially oriented elastin struts, based on both the aspect ratio and transmural tilt. MVICs were observed to have a large cytoplasmic volume, as evidenced by the large mean surface area to volume ratio 3.68 ± 0.35, which increased under physiological loading conditions to 4.98 ± 1.17.

CONCLUSIONS

Our findings suggest that each MVIC mechanically interacted only with the nearest 3-4 collagen fibers. This key observation suggests that in developing multiscale MV models, each MVIC can be considered a mechanically integral part of the local fiber ensemble and is unlikely to be influenced by more distant structures.