Bayesian optimisation for efficient parameter inference in a cardiac mechanics model of the left ventricle

Advanced statistical inference of myocardial stiffness: A time series Gaussian process approach of emulating cardiac mechanics for real-time clinical decision support

Cardiac mechanics modelling promises to revolutionize personalized health care; however, inferring patient-specific biophysical parameters, which are critical for understanding myocardial functions and performance, poses substantial methodological challenges. Our work is primarily motivated to determine the passive stiffness of the myocardium from the measurement of the left ventricle (LV) volume at various time points, which is crucial for diagnosing cardiac physiological conditions. Although there have been significant advancements in cardiac mechanics modelling, the tasks of inference and uncertainty quantification of myocardial stiffness remain challenging, with high computational costs preventing real-time decision support. We adapt Gaussian processes to construct a statistical surrogate model for emulating LV cavity volume during diastolic filling to overcome this challenge. As the LV volumes, obtained at different time points in diastole, constitute a time series, we apply the Kronecker product trick to decompose the complex covariance matrix of the whole system into two separate covariance matrices, one for time and the other for biophysical parameters. To proceed towards personalized health care, we further integrate patient-specific LV geometries into the Gaussian process emulator using principal component analysis (PCA). Utilizing a deep learning neural network for extracting time-series left ventricle volumes from magnetic resonance images, Bayesian inference is applied to determine the posterior probability distribution of critical cardiac mechanics parameters. Tests on real-patient data illustrate the potential for real-time estimation of myocardial properties for clinical decision-making. These advancements constitute a crucial step towards clinical impact, offering valuable insights into posterior uncertainty quantification for complex cardiac mechanics models.

Personalized evaluation of the passive myocardium in ischemic cardiomyopathy via computational modeling using Bayesian optimization

Biomechanics-based patient-specific modeling is a promising approach that has proved invaluable for its clinical potential to assess the adversities caused by ischemic heart disease (IHD). In the present study, we propose a framework to find the passive material properties of the myocardium and the unloaded shape of cardiac ventricles simultaneously in patients diagnosed with ischemic cardiomyopathy (ICM). This was achieved by minimizing the difference between the simulated and the target end-diastolic pressure-volume relationships (EDPVRs) using black-box Bayesian optimization, based on the finite element analysis (FEA). End-diastolic (ED) biventricular geometry and the location of the ischemia were determined from cardiac magnetic resonance (CMR) imaging. We employed our pipeline to model the cardiac ventricles of three patients aged between 57 and 66 years, with and without the inclusion of valves. An excellent agreement between the simulated and the target EDPVRs has been reached. Our results revealed that the incorporation of valvular springs typically leads to lower hyperelastic parameters for both healthy and ischemic myocardium, as well as a higher fiber Green strain in the viable regions compared to models without valvular stiffness. Furthermore, the addition of valve-related effects did not result in significant changes in myofiber stress after optimization. We concluded that more accurate results could be obtained when cardiac valves were considered in modeling ventricles. The present novel and practical methodology paves the way for developing digital twins of ischemic cardiac ventricles, providing a non-invasive assessment for designing optimal personalized therapies in precision medicine.

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.

Image-based estimation of the left ventricular cavity volume using deep learning and Gaussian process with cardio-mechanical applications

In this investigation, an image-based method has been developed to estimate the volume of the left ventricular cavity using cardiac magnetic resonance (CMR) imaging data. Deep learning and Gaussian processes have been applied to bring the estimations closer to the cavity volumes manually extracted. CMR data from 339 patients and healthy volunteers have been used to train a stepwise regression model that can estimate the volume of the left ventricular cavity at the beginning and end of diastole. We have decreased the root mean square error (RMSE) of cavity volume estimation approximately from 13 to 8 ml compared to the common practice in the literature. Considering the RMSE of manual measurements is approximately 4 ml on the same dataset, 8 ml of error is notable for a fully automated estimation method, which needs no supervision or user-hours once it has been trained. Additionally, to demonstrate a clinically relevant application of automatically estimated volumes, we inferred the passive material properties of the myocardium given the volume estimates using a well-validated cardiac model. These material properties can be further used for patient treatment planning and diagnosis.