1
|
Zaman MS, Dhamala J, Bajracharya P, Sapp JL, Horácek BM, Wu KC, Trayanova NA, Wang L. Fast Posterior Estimation of Cardiac Electrophysiological Model Parameters via Bayesian Active Learning. Front Physiol 2021; 12:740306. [PMID: 34759835 PMCID: PMC8573318 DOI: 10.3389/fphys.2021.740306] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 09/24/2021] [Indexed: 11/13/2022] Open
Abstract
Probabilistic estimation of cardiac electrophysiological model parameters serves an important step toward model personalization and uncertain quantification. The expensive computation associated with these model simulations, however, makes direct Markov Chain Monte Carlo (MCMC) sampling of the posterior probability density function (pdf) of model parameters computationally intensive. Approximated posterior pdfs resulting from replacing the simulation model with a computationally efficient surrogate, on the other hand, have seen limited accuracy. In this study, we present a Bayesian active learning method to directly approximate the posterior pdf function of cardiac model parameters, in which we intelligently select training points to query the simulation model in order to learn the posterior pdf using a small number of samples. We integrate a generative model into Bayesian active learning to allow approximating posterior pdf of high-dimensional model parameters at the resolution of the cardiac mesh. We further introduce new acquisition functions to focus the selection of training points on better approximating the shape rather than the modes of the posterior pdf of interest. We evaluated the presented method in estimating tissue excitability in a 3D cardiac electrophysiological model in a range of synthetic and real-data experiments. We demonstrated its improved accuracy in approximating the posterior pdf compared to Bayesian active learning using regular acquisition functions, and substantially reduced computational cost in comparison to existing standard or accelerated MCMC sampling.
Collapse
Affiliation(s)
- Md Shakil Zaman
- Rochester Institute of Technology, Rochester, NY, United States
| | - Jwala Dhamala
- Rochester Institute of Technology, Rochester, NY, United States
| | | | - John L Sapp
- Department of Medicine, Dalhousie University, Halifax, NS, Canada
| | - B Milan Horácek
- Department of Electrical and Computer Engineering, Halifax, NS, Canada
| | - Katherine C Wu
- Department of Medicine, Johns Hopkins University, Baltimore, MD, United States
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States
| | - Linwei Wang
- Rochester Institute of Technology, Rochester, NY, United States
| |
Collapse
|
2
|
Dabiri Y, Van der Velden A, Sack KL, Choy JS, Kassab GS, Guccione JM. Prediction of Left Ventricular Mechanics Using Machine Learning. FRONTIERS IN PHYSICS 2019; 7:117. [PMID: 31903394 PMCID: PMC6941671 DOI: 10.3389/fphy.2019.00117] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The goal of this paper was to provide a real-time left ventricular (LV) mechanics simulator using machine learning (ML). Finite element (FE) simulations were conducted for the LV with different material properties to obtain a training set. A hyperelastic fiber-reinforced material model was used to describe the passive behavior of the myocardium during diastole. The active behavior of the heart resulting from myofiber contractions was added to the passive tissue during systole. The active and passive properties govern the LV constitutive equation. These mechanical properties were altered using optimal Latin hypercube design of experiments to obtain training FE models with varied active properties (volume and pressure predictions) and varied passive properties (stress predictions). For prediction of LV pressures, we used eXtreme Gradient Boosting (XGboost) and Cubist, and XGBoost was used for predictions of LV pressures, volumes as well as LV stresses. The LV pressure and volume results obtained from ML were similar to FE computations. The ML results could capture the shape of LV pressure as well as LV pressure-volume loops. The results predicted by Cubist were smoother than those from XGBoost. The mean absolute errors were as follows: XGBoost volume: 1.734 ± 0.584 ml, XGBoost pressure: 1.544 ± 0.298 mmHg, Cubist volume: 1.495 ± 0.260 ml, Cubist pressure: 1.623 ± 0.191 mmHg, myofiber stress: 0.334 ± 0.228 kPa, cross myofiber stress: 0.075 ± 0.024 kPa, and shear stress: 0.050 ± 0.032 kPa. The simulation results show ML can predict LV mechanics much faster than the FE method. The ML model can be used as a tool to predict LV behavior. Training of our ML model based on a large group of subjects can improve its predictability for real world applications.
Collapse
Affiliation(s)
- Yaghoub Dabiri
- California Medical Innovations Institute, San Diego, CA, United States
- Department of Surgery, University of California, San Francisco, San Francisco, CA, United States
| | | | - Kevin L. Sack
- Department of Surgery, University of California, San Francisco, San Francisco, CA, United States
- Division of Biomedical Engineering, Department of Human Biology, University of Cape Town, Cape Town, South Africa
| | - Jenny S. Choy
- California Medical Innovations Institute, San Diego, CA, United States
| | - Ghassan S. Kassab
- California Medical Innovations Institute, San Diego, CA, United States
| | - Julius M. Guccione
- Department of Surgery, University of California, San Francisco, San Francisco, CA, United States
| |
Collapse
|
3
|
Molléro R, Pennec X, Delingette H, Ayache N, Sermesant M. Population-based priors in cardiac model personalisation for consistent parameter estimation in heterogeneous databases. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3158. [PMID: 30239175 DOI: 10.1002/cnm.3158] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Revised: 09/10/2018] [Accepted: 09/16/2018] [Indexed: 06/08/2023]
Abstract
Personalised cardiac models are a virtual representation of the patient heart, with parameter values for which the simulation fits the available clinical measurements. Models usually have a large number of parameters while the available data for a given patient are typically limited to a small set of measurements; thus, the parameters cannot be estimated uniquely. This is a practical obstacle for clinical applications, where accurate parameter values can be important. Here, we explore an original approach based on an algorithm called Iteratively Updated Priors (IUP), in which we perform successive personalisations of a full database through maximum a posteriori (MAP) estimation, where the prior probability at an iteration is set from the distribution of personalised parameters in the database at the previous iteration. At the convergence of the algorithm, estimated parameters of the population lie on a linear subspace of reduced (and possibly sufficient) dimension in which for each case of the database, there is a (possibly unique) parameter value for which the simulation fits the measurements. We first show how this property can help the modeller select a relevant parameter subspace for personalisation. In addition, since the resulting priors in this subspace represent the population statistics in this subspace, they can be used to perform consistent parameter estimation for cases where measurements are possibly different or missing in the database, which we illustrate with the personalisation of a heterogeneous database of 811 cases.
Collapse
Affiliation(s)
- Roch Molléro
- Inria, Epione Research Project, Sophia Antipolis, France
| | - Xavier Pennec
- Inria, Epione Research Project, Sophia Antipolis, France
| | | | | | | |
Collapse
|
4
|
Gray RA, Pathmanathan P. Patient-Specific Cardiovascular Computational Modeling: Diversity of Personalization and Challenges. J Cardiovasc Transl Res 2018; 11:80-88. [PMID: 29512059 PMCID: PMC5908828 DOI: 10.1007/s12265-018-9792-2] [Citation(s) in RCA: 73] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Accepted: 02/02/2018] [Indexed: 02/07/2023]
Abstract
Patient-specific computer models have been developed representing a variety of aspects of the cardiovascular system spanning the disciplines of electrophysiology, electromechanics, solid mechanics, and fluid dynamics. These physiological mechanistic models predict macroscopic phenomena such as electrical impulse propagation and contraction throughout the entire heart as well as flow and pressure dynamics occurring in the ventricular chambers, aorta, and coronary arteries during each heartbeat. Such models have been used to study a variety of clinical scenarios including aortic aneurysms, coronary stenosis, cardiac valvular disease, left ventricular assist devices, cardiac resynchronization therapy, ablation therapy, and risk stratification. After decades of research, these models are beginning to be incorporated into clinical practice directly via marketed devices and indirectly by improving our understanding of the underlying mechanisms of health and disease within a clinical context.
Collapse
Affiliation(s)
- Richard A Gray
- Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, Food and Drug Administration, Silver Spring, MD, 20993, USA.
- , Silver Spring, USA.
| | - Pras Pathmanathan
- Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, Food and Drug Administration, Silver Spring, MD, 20993, USA
| |
Collapse
|
5
|
Balaban G, Finsberg H, Odland HH, Rognes ME, Ross S, Sundnes J, Wall S. High-resolution data assimilation of cardiac mechanics applied to a dyssynchronous ventricle. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2017; 33:e2863. [PMID: 28039961 DOI: 10.1002/cnm.2863] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Revised: 10/31/2016] [Accepted: 12/28/2016] [Indexed: 06/06/2023]
Abstract
Computational models of cardiac mechanics, personalized to a patient, offer access to mechanical information above and beyond direct medical imaging. Additionally, such models can be used to optimize and plan therapies in-silico, thereby reducing risks and improving patient outcome. Model personalization has traditionally been achieved by data assimilation, which is the tuning or optimization of model parameters to match patient observations. Current data assimilation procedures for cardiac mechanics are limited in their ability to efficiently handle high-dimensional parameters. This restricts parameter spatial resolution, and thereby the ability of a personalized model to account for heterogeneities that are often present in a diseased or injured heart. In this paper, we address this limitation by proposing an adjoint gradient-based data assimilation method that can efficiently handle high-dimensional parameters. We test this procedure on a synthetic data set and provide a clinical example with a dyssynchronous left ventricle with highly irregular motion. Our results show that the method efficiently handles a high-dimensional optimization parameter and produces an excellent agreement for personalized models to both synthetic and clinical data.
Collapse
Affiliation(s)
- Gabriel Balaban
- Simula Research Laboratory, P.O. Box 134 1325 Lysaker, Norway
- Department of Informatics, University of Oslo, P.O. Box 1080, Blindern 0316 Oslo, Norway
- Center for Cardiological Innovation, Songsvannsveien 9, 0372 Oslo, Norway
| | - Henrik Finsberg
- Simula Research Laboratory, P.O. Box 134 1325 Lysaker, Norway
- Department of Informatics, University of Oslo, P.O. Box 1080, Blindern 0316 Oslo, Norway
- Center for Cardiological Innovation, Songsvannsveien 9, 0372 Oslo, Norway
| | - Hans Henrik Odland
- Faculty of Medicine, University of Oslo, P.O. Box 1078 Blindern, 0316 Oslo, Norway
- Department of Pediatrics, Oslo University Hospital, PO Nydalen, Oslo, Norway
| | - Marie E Rognes
- Simula Research Laboratory, P.O. Box 134 1325 Lysaker, Norway
- Department of Mathematics, University of Oslo, P.O Box 1053, Blindern 0316 Oslo, Norway
| | - Stian Ross
- Faculty of Medicine, University of Oslo, P.O. Box 1078 Blindern, 0316 Oslo, Norway
- Center for Cardiological Innovation, Songsvannsveien 9, 0372 Oslo, Norway
| | - Joakim Sundnes
- Simula Research Laboratory, P.O. Box 134 1325 Lysaker, Norway
- Department of Informatics, University of Oslo, P.O. Box 1080, Blindern 0316 Oslo, Norway
- Center for Cardiological Innovation, Songsvannsveien 9, 0372 Oslo, Norway
| | - Samuel Wall
- Simula Research Laboratory, P.O. Box 134 1325 Lysaker, Norway
- Center for Cardiological Innovation, Songsvannsveien 9, 0372 Oslo, Norway
- Department of Mathematical Science and Technology, Norwegian University of Life Sciences, Universitetstunet 3 1430 Ås, Norway
| |
Collapse
|
6
|
Duchateau N, De Craene M, Allain P, Saloux E, Sermesant M. Infarct Localization From Myocardial Deformation: Prediction and Uncertainty Quantification by Regression From a Low-Dimensional Space. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:2340-2352. [PMID: 27164583 DOI: 10.1109/tmi.2016.2562181] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Diagnosing and localizing myocardial infarct is crucial for early patient management and therapy planning. We propose a new method for predicting the location of myocardial infarct from local wall deformation, which has value for risk stratification from routine examinations such as (3D) echocardiography. The pipeline combines non-linear dimensionality reduction of deformation patterns and two multi-scale kernel regressions. Confidence in the diagnosis is assessed by a map of local uncertainties, which integrates plausible infarct locations generated from the space of reduced dimensionality. These concepts were tested on 500 synthetic cases generated from a realistic cardiac electromechanical model, and 108 pairs of 3D echocardiographic sequences and delayed-enhancement magnetic resonance images from real cases. Infarct prediction is made at a spatial resolution around 4 mm, more than 10 times smaller than the current diagnosis, made regionally. Our method is accurate, and significantly outperforms the clinically-used thresholding of the deformation patterns (on real data: sensitivity/specificity of 0.828/0.804, area under the curve: 0.909 versus 0.742 for the most predictive strain component). Uncertainty adds value to refine the diagnosis and eventually re-examine suspicious cases.
Collapse
|
7
|
Chabiniok R, Wang VY, Hadjicharalambous M, Asner L, Lee J, Sermesant M, Kuhl E, Young AA, Moireau P, Nash MP, Chapelle D, Nordsletten DA. Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics. Interface Focus 2016; 6:20150083. [PMID: 27051509 PMCID: PMC4759748 DOI: 10.1098/rsfs.2015.0083] [Citation(s) in RCA: 139] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.
Collapse
Affiliation(s)
- Radomir Chabiniok
- Division of Imaging Sciences and Biomedical Engineering, King's College London, St Thomas’ Hospital, London SE1 7EH, UK
- Inria and Paris-Saclay University, Bâtiment Alan Turing, 1 rue Honoré d'Estienne d'Orves, Campus de l'Ecole Polytechnique, Palaiseau 91120, France
| | - Vicky Y. Wang
- Auckland Bioengineering Institute, University of Auckland, 70 Symonds Street, Auckland, New Zealand
| | - Myrianthi Hadjicharalambous
- Division of Imaging Sciences and Biomedical Engineering, King's College London, St Thomas’ Hospital, London SE1 7EH, UK
| | - Liya Asner
- Division of Imaging Sciences and Biomedical Engineering, King's College London, St Thomas’ Hospital, London SE1 7EH, UK
| | - Jack Lee
- Division of Imaging Sciences and Biomedical Engineering, King's College London, St Thomas’ Hospital, London SE1 7EH, UK
| | - Maxime Sermesant
- Inria, Asclepios team, 2004 route des Lucioles BP 93, Sophia Antipolis Cedex 06902, France
| | - Ellen Kuhl
- Departments of Mechanical Engineering, Bioengineering, and Cardiothoracic Surgery, Stanford University, 496 Lomita Mall, Durand 217, Stanford, CA 94306, USA
| | - Alistair A. Young
- Auckland Bioengineering Institute, University of Auckland, 70 Symonds Street, Auckland, New Zealand
| | - Philippe Moireau
- Inria and Paris-Saclay University, Bâtiment Alan Turing, 1 rue Honoré d'Estienne d'Orves, Campus de l'Ecole Polytechnique, Palaiseau 91120, France
| | - Martyn P. Nash
- Auckland Bioengineering Institute, University of Auckland, 70 Symonds Street, Auckland, New Zealand
- Department of Engineering Science, University of Auckland, 70 Symonds Street, Auckland, New Zealand
| | - Dominique Chapelle
- Inria and Paris-Saclay University, Bâtiment Alan Turing, 1 rue Honoré d'Estienne d'Orves, Campus de l'Ecole Polytechnique, Palaiseau 91120, France
| | - David A. Nordsletten
- Division of Imaging Sciences and Biomedical Engineering, King's College London, St Thomas’ Hospital, London SE1 7EH, UK
| |
Collapse
|
8
|
Kayvanpour E, Mansi T, Sedaghat-Hamedani F, Amr A, Neumann D, Georgescu B, Seegerer P, Kamen A, Haas J, Frese KS, Irawati M, Wirsz E, King V, Buss S, Mereles D, Zitron E, Keller A, Katus HA, Comaniciu D, Meder B. Towards Personalized Cardiology: Multi-Scale Modeling of the Failing Heart. PLoS One 2015; 10:e0134869. [PMID: 26230546 PMCID: PMC4521877 DOI: 10.1371/journal.pone.0134869] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2015] [Accepted: 07/14/2015] [Indexed: 01/14/2023] Open
Abstract
BACKGROUND Despite modern pharmacotherapy and advanced implantable cardiac devices, overall prognosis and quality of life of HF patients remain poor. This is in part due to insufficient patient stratification and lack of individualized therapy planning, resulting in less effective treatments and a significant number of non-responders. METHODS AND RESULTS State-of-the-art clinical phenotyping was acquired, including magnetic resonance imaging (MRI) and biomarker assessment. An individualized, multi-scale model of heart function covering cardiac anatomy, electrophysiology, biomechanics and hemodynamics was estimated using a robust framework. The model was computed on n=46 HF patients, showing for the first time that advanced multi-scale models can be fitted consistently on large cohorts. Novel multi-scale parameters derived from the model of all cases were analyzed and compared against clinical parameters, cardiac imaging, lab tests and survival scores to evaluate the explicative power of the model and its potential for better patient stratification. Model validation was pursued by comparing clinical parameters that were not used in the fitting process against model parameters. CONCLUSION This paper illustrates how advanced multi-scale models can complement cardiovascular imaging and how they could be applied in patient care. Based on obtained results, it becomes conceivable that, after thorough validation, such heart failure models could be applied for patient management and therapy planning in the future, as we illustrate in one patient of our cohort who received CRT-D implantation.
Collapse
Affiliation(s)
- Elham Kayvanpour
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Tommaso Mansi
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Farbod Sedaghat-Hamedani
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Ali Amr
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Dominik Neumann
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Bogdan Georgescu
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Philipp Seegerer
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Ali Kamen
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Jan Haas
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Karen S. Frese
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
| | - Maria Irawati
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
| | - Emil Wirsz
- Siemens AG, Corporate Technology, Erlangen, Germany
| | - Vanessa King
- Siemens Corporation, Corporate Technology, Sensor Technologies, Princeton, New Jersey, United States of America
| | - Sebastian Buss
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
| | - Derliz Mereles
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
| | - Edgar Zitron
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
| | - Andreas Keller
- Biomarker Discovery Center Heidelberg, Heidelberg, Germany
- Department of Human Genetics, Saarland University, Homburg, Germany
| | - Hugo A. Katus
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
- Klaus Tschira Institute for Computational Cardiology, Heidelberg, Germany
| | - Dorin Comaniciu
- Siemens Corporation, Corporate Technology, Imaging and Computer Vision, Princeton, New Jersey, United States of America
| | - Benjamin Meder
- Department of Medicine III, University of Heidelberg, Heidelberg, Germany
- DZHK (German Centre for Cardiovascular Research), Heidelberg, Germany
- Klaus Tschira Institute for Computational Cardiology, Heidelberg, Germany
| |
Collapse
|