Application of feed forward and recurrent neural networks in simulation of left ventricular mechanics

Material characterization of human middle ear using machine-learning-based surrogate models

This study aims to introduce a novel non-invasive method for rapid material characterization of middle-ear structures, taking into consideration the invaluable insights provided by the mechanical properties of ear tissues. Valuable insights into various ear pathologies can be gleaned from the mechanical properties of ear tissues, yet conventional techniques for assessing these properties often entail invasive procedures that preclude their use on living patients. In this study, in the first step, we developed machine-learning models of the middle ear to predict its responses with a significantly lower computational cost in comparison to finite-element models. Leveraging findings from prior research, we focused on the most influential model parameters: the Young's modulus and thickness of the tympanic membrane and the Young's modulus of the stapedial annular ligament. The eXtreme Gradient Boosting (XGBoost) method was implemented for creating the machine-learning models. Subsequently, we combined the created machine-learning models with Bayesian optimization (BoTorch) for fast and efficient estimation of the Young's moduli of the tympanic membrane and the stapedial annular ligament. We demonstrate that the resultant surrogate models can fairly represent the vibrational responses of the umbo, stapes footplate, and vibration patterns of the tympanic membrane at most frequencies. Also, our proposed material characterization approach successfully estimated the Young's moduli of the tympanic membrane and stapedial annular ligament (separately and simultaneously) with values of mean absolute percentage error of less than 7%. The remarkable accuracy achieved through the proposed material characterization method underscores its potential for eventual clinical applications of estimating mechanical properties of the middle-ear structures for diagnostic purposes.

Machine learning used for simulation of MitraClip intervention: A proof-of-concept study

Introduction: Severe mitral regurgitation (MR) is a mitral valve disease that can lead to lifethreatening complications. MitraClip (MC) therapy is a percutaneous solution for patients who cannot tolerate surgical solutions. In MC therapy, a clip is implanted in the heart to reduce MR. To achieve optimal MC therapy, the cardiologist needs to foresee the outcomes of different scenarios for MC implantation, including the location of the MC. Although finite element (FE) modeling can simulate the outcomes of different MC scenarios, it is not suitable for clinical usage because it requires several hours to complete.

Methods: In this paper, we used machine learning (ML) to predict the outcomes of MC therapy in less than 1 s. Two ML algorithms were used: XGBoost, which is a decision tree model, and a feed-forward deep learning (DL) model. The MC location, the geometrical attributes of the models and baseline stress and MR were the features of the ML models, and the predictions were performed for MR and maximum von Mises stress in the leaflets. The parameters of the ML models were determined to achieve the minimum errors obtained by applying the ML models on the validation set.

Results: The results for the test set (not used during training) showed relative agreement between ML predictions and ground truth FE predictions. The accuracy of the XGBoost models were better than DL models. Mean absolute percentage error (MAPE) for the XGBoost predictions were 0.115 and 0.231, and the MAPE for DL predictions were 0.154 and 0.310, for MR and stress, respectively.

Discussion: The ML models reduced the FE runtime from 6 hours (on average) to less than 1 s. The accuracy of ML models can be increased by increasing the dataset size. The results of this study have important implications for improving the outcomes of MC therapy by providing information about the outcomes of MC implantation in real-time.

A machine learning model to estimate myocardial stiffness from EDPVR

In-vivo estimation of mechanical properties of the myocardium is essential for patient-specific diagnosis and prognosis of cardiac disease involving myocardial remodeling, including myocardial infarction and heart failure with preserved ejection fraction. Current approaches use time-consuming finite-element (FE) inverse methods that involve reconstructing and meshing the heart geometry, imposing measured loading, and conducting computationally expensive iterative FE simulations. In this paper, we propose a machine learning (ML) model that feasibly and accurately predicts passive myocardial properties directly from select geometric, architectural, and hemodynamic measures, thus bypassing exhaustive steps commonly required in cardiac FE inverse problems. Geometric and fiber-orientation features were chosen to be readily obtainable from standard cardiac imaging protocols. The end-diastolic pressure-volume relationship (EDPVR), which can be obtained using a single-point pressure-volume measurement, was used as a hemodynamic (loading) feature. A comprehensive ML training dataset in the geometry-architecture-loading space was generated, including a wide variety of partially synthesized rodent heart geometry and myofiber helicity possibilities, and a broad range of EDPVRs obtained using forward FE simulations. Latin hypercube sampling was used to create 2500 examples for training, validation, and testing. A multi-layer feed-forward neural network (MFNN) was used as a deep learning agent to train the ML model. The model showed excellent performance in predicting stiffness parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_f$$\end{document}af and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_f$$\end{document}bf associated with fiber direction (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^2_{a_f}=99.471\%$$\end{document}Raf2=99.471% and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^2_{b_f}=92.837\%$$\end{document}Rbf2=92.837%). After conducting permutation feature importance analysis, the ML performance further improved for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_f$$\end{document}bf (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^2_{b_f}=96.240\%$$\end{document}Rbf2=96.240%), and the left ventricular volume and endocardial area were found to be the most critical geometric features for accurate predictions. The ML model predictions were evaluated further in two cases: (i) rat-specific stiffness data measured using ex-vivo mechanical testing, and (ii) patient-specific estimation using FE inverse modeling. Excellent agreements with ML predictions were found for both cases. The trained ML model offers a feasible technology to estimate patient-specific myocardial properties, thus, bridging the gap between EDPVR, as a confounded organ-level metric for tissue stiffness, and intrinsic tissue-level properties. These properties provide incremental information relative to traditional organ-level indices for cardiac function, improving the clinical assessment and prognosis of cardiac diseases.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.

Mitral Valve Atlas for Artificial Intelligence Predictions of MitraClip Intervention Outcomes

Severe mitral regurgitation (MR) is a cardiac disease that can lead to fatal consequences. MitraClip (MC) intervention is a percutaneous procedure whereby the mitral valve (MV) leaflets are connected along the edge using MCs. The outcomes of the MC intervention are not known in advance, i.e., the outcomes are quite variable. Artificial intelligence (AI) can be used to guide the cardiologist in selecting optimal MC scenarios. In this study, we describe an atlas of shapes as well as different scenarios for MC implantation for such an AI analysis. We generated the MV geometrical data from three different sources. First, the patients' 3-dimensional echo images were used. The pixel data from six key points were obtained from three views of the echo images. Using PyGem, an open-source morphing library in Python, these coordinates were used to create the geometry by morphing a template geometry. Second, the dimensions of the MV, from the literature were used to create data. Third, we used machine learning methods, principal component analysis, and generative adversarial networks to generate more shapes. We used the finite element (FE) software ABAQUS to simulate smoothed particle hydrodynamics in different scenarios for MC intervention. The MR and stresses in the leaflets were post-processed. Our physics-based FE models simulated the outcomes of MC intervention for different scenarios. The MR and stresses in the leaflets were computed by the FE models for a single clip at different locations as well as two and three clips. Results from FE simulations showed that the location and number of MCs affect subsequent residual MR, and that leaflet stresses do not follow a simple pattern. Furthermore, FE models need several hours to provide the results, and they are not applicable for clinical usage where the predicted outcomes of MC therapy are needed in real-time. In this study, we generated the required dataset for the AI models which can provide the results in a matter of seconds.

Multiscale simulations of left ventricular growth and remodeling

Cardiomyocytes can adapt their size, shape, and orientation in response to altered biomechanical or biochemical stimuli. The process by which the heart undergoes structural changes-affecting both geometry and material properties-in response to altered ventricular loading, altered hormonal levels, or mutant sarcomeric proteins is broadly known as cardiac growth and remodeling (G&R). Although it is likely that cardiac G&R initially occurs as an adaptive response of the heart to the underlying stimuli, prolonged pathological changes can lead to increased risk of atrial fibrillation, heart failure, and sudden death. During the past few decades, computational models have been extensively used to investigate the mechanisms of cardiac G&R, as a complement to experimental measurements. These models have provided an opportunity to quantitatively study the relationships between the underlying stimuli (primarily mechanical) and the adverse outcomes of cardiac G&R, i.e., alterations in ventricular size and function. State-of-the-art computational models have shown promise in predicting the progression of cardiac G&R. However, there are still limitations that need to be addressed in future works to advance the field. In this review, we first outline the current state of computational models of cardiac growth and myofiber remodeling. Then, we discuss the potential limitations of current models of cardiac G&R that need to be addressed before they can be utilized in clinical care. Finally, we briefly discuss the next feasible steps and future directions that could advance the field of cardiac G&R.