1
|
Koopsen T, van Osta N, van Loon T, Meiburg R, Huberts W, Beela AS, Kirkels FP, van Klarenbosch BR, Teske AJ, Cramer MJ, Bijvoet GP, van Stipdonk A, Vernooy K, Delhaas T, Lumens J. Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination. Biomed Eng Online 2024; 23:46. [PMID: 38741182 DOI: 10.1186/s12938-024-01232-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Accepted: 04/02/2024] [Indexed: 05/16/2024] Open
Abstract
BACKGROUND Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error (χ 2 ) of LV myocardial strain, strain rate, and cavity volume. RESULTS A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients (χ 2 < 1.6), but minimum parameter reproducibility was poor (ICC min = 0.01). Iterative reduction yielded a reproducible (ICC min = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs (χ 2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). CONCLUSIONS By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.
Collapse
Affiliation(s)
- Tijmen Koopsen
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands.
| | - Nick van Osta
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
| | - Tim van Loon
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
| | - Roel Meiburg
- Group SIMBIOTX, Institut de Recherche en Informatique et en Automatique (INRIA), Paris, France
| | - Wouter Huberts
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Ahmed S Beela
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
- Department of Cardiology, Suez Canal University, Ismailia, Egypt
| | - Feddo P Kirkels
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht (UMCU), Utrecht, The Netherlands
| | - Bas R van Klarenbosch
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht (UMCU), Utrecht, The Netherlands
| | - Arco J Teske
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht (UMCU), Utrecht, The Netherlands
| | - Maarten J Cramer
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht (UMCU), Utrecht, The Netherlands
| | - Geertruida P Bijvoet
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
- Department of Cardiology, Maastricht University Medical Center (MUMC), Maastricht, The Netherlands
| | - Antonius van Stipdonk
- Department of Cardiology, Maastricht University Medical Center (MUMC), Maastricht, The Netherlands
| | - Kevin Vernooy
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
- Department of Cardiology, Maastricht University Medical Center (MUMC), Maastricht, The Netherlands
- Department of Cardiology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands.
| |
Collapse
|
2
|
Salvador M, Strocchi M, Regazzoni F, Augustin CM, Dede' L, Niederer SA, Quarteroni A. Whole-heart electromechanical simulations using Latent Neural Ordinary Differential Equations. NPJ Digit Med 2024; 7:90. [PMID: 38605089 PMCID: PMC11009296 DOI: 10.1038/s41746-024-01084-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Accepted: 03/22/2024] [Indexed: 04/13/2024] Open
Abstract
Cardiac digital twins provide a physics and physiology informed framework to deliver personalized medicine. However, high-fidelity multi-scale cardiac models remain a barrier to adoption due to their extensive computational costs. Artificial Intelligence-based methods can make the creation of fast and accurate whole-heart digital twins feasible. We use Latent Neural Ordinary Differential Equations (LNODEs) to learn the pressure-volume dynamics of a heart failure patient. Our surrogate model is trained from 400 simulations while accounting for 43 parameters describing cell-to-organ cardiac electromechanics and cardiovascular hemodynamics. LNODEs provide a compact representation of the 3D-0D model in a latent space by means of an Artificial Neural Network that retains only 3 hidden layers with 13 neurons per layer and allows for numerical simulations of cardiac function on a single processor. We employ LNODEs to perform global sensitivity analysis and parameter estimation with uncertainty quantification in 3 hours of computations, still on a single processor.
Collapse
Affiliation(s)
- Matteo Salvador
- Institute for Computational and Mathematical Engineering, Stanford University, California, CA, USA.
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy.
| | - Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
| | | | - Christoph M Augustin
- Institute of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Luca Dede'
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
- The Alan Turing Institute, London, UK
| | - Alfio Quarteroni
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy
- École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
3
|
Rodero C, Baptiste TMG, Barrows RK, Lewalle A, Niederer SA, Strocchi M. Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways. FRONTIERS IN PHYSICS 2023; 11:1306210. [PMID: 38500690 PMCID: PMC7615748 DOI: 10.3389/fphy.2023.1306210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.
Collapse
Affiliation(s)
- Cristobal Rodero
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| | - Tiffany M. G. Baptiste
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King’s College London, London, United Kingdom
| | - Rosie K. Barrows
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King’s College London, London, United Kingdom
| | - Alexandre Lewalle
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| | - Steven A. Niederer
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Turing Research and Innovation Cluster in Digital Twins (TRIC: DT), The Alan Turing Institute, London, United Kingdom
| | - Marina Strocchi
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| |
Collapse
|
4
|
Tunedal K, Viola F, Garcia BC, Bolger A, Nyström FH, Östgren CJ, Engvall J, Lundberg P, Dyverfeldt P, Carlhäll CJ, Cedersund G, Ebbers T. Haemodynamic effects of hypertension and type 2 diabetes: Insights from a 4D flow MRI-based personalized cardiovascular mathematical model. J Physiol 2023; 601:3765-3787. [PMID: 37485733 DOI: 10.1113/jp284652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 06/29/2023] [Indexed: 07/25/2023] Open
Abstract
Type 2 diabetes (T2D) and hypertension increase the risk of cardiovascular diseases mediated by whole-body changes to metabolism, cardiovascular structure and haemodynamics. The haemodynamic changes related to hypertension and T2D are complex and subject-specific, however, and not fully understood. We aimed to investigate the haemodynamic mechanisms in T2D and hypertension by comparing the haemodynamics between healthy controls and subjects with T2D, hypertension, or both. For all subjects, we combined 4D flow magnetic resonance imaging data, brachial blood pressure and a cardiovascular mathematical model to create a comprehensive subject-specific analysis of central haemodynamics. When comparing the subject-specific haemodynamic parameters between the four groups, the predominant haemodynamic difference is impaired left ventricular relaxation in subjects with both T2D and hypertension compared to subjects with only T2D, only hypertension and controls. The impaired relaxation indicates that, in this cohort, the long-term changes in haemodynamic load of co-existing T2D and hypertension cause diastolic dysfunction demonstrable at rest, whereas either disease on its own does not. However, through subject-specific predictions of impaired relaxation, we show that altered relaxation alone is not enough to explain the subject-specific and group-related differences; instead, a combination of parameters is affected in T2D and hypertension. These results confirm previous studies that reported more adverse effects from the combination of T2D and hypertension compared to either disease on its own. Furthermore, this shows the potential of personalized cardiovascular models in providing haemodynamic mechanistic insights and subject-specific predictions that could aid in the understanding and treatment planning of patients with T2D and hypertension. KEY POINTS: The combination of 4D flow magnetic resonance imaging data and a cardiovascular mathematical model allows for a comprehensive analysis of subject-specific haemodynamic parameters that otherwise cannot be derived non-invasively. Using this combination, we show that diastolic dysfunction in subjects with both type 2 diabetes (T2D) and hypertension is the main group-level difference between controls, subjects with T2D, subjects with hypertension, and subjects with both T2D and hypertension. These results suggest that, in this relatively healthy population, the additional load of both hypertension and T2D affects the haemodynamic function of the left ventricle, whereas each disease on its own is not enough to cause significant effects under resting conditions. Finally, using the subject-specific model, we show that the haemodynamic effects of diastolic dysfunction alone are not sufficient to explain all the observed haemodynamic differences. Instead, additional subject-specific variations in cardiac and vascular function combine to explain the complex haemodynamics of subjects affected by hypertension and/or T2D.
Collapse
Affiliation(s)
- Kajsa Tunedal
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Federica Viola
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Belén Casas Garcia
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Ann Bolger
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Department of Medicine, University of California San Francisco, San Francisco, CA, USA
| | - Fredrik H Nyström
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Carl Johan Östgren
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Prevention, Rehabilitation and Community Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Jan Engvall
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Department of Clinical Physiology in Linköping, and Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Peter Lundberg
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Department of Radiation Physics, and Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Petter Dyverfeldt
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Carl-Johan Carlhäll
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Department of Clinical Physiology in Linköping, and Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Tino Ebbers
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- Division of Diagnostics and Specialist Medicine, Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| |
Collapse
|
5
|
Sander J, de Vos BD, Bruns S, Planken N, Viergever MA, Leiner T, Išgum I. Reconstruction and completion of high-resolution 3D cardiac shapes using anisotropic CMRI segmentations and continuous implicit neural representations. Comput Biol Med 2023; 164:107266. [PMID: 37494823 DOI: 10.1016/j.compbiomed.2023.107266] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Revised: 06/26/2023] [Accepted: 07/16/2023] [Indexed: 07/28/2023]
Abstract
Since the onset of computer-aided diagnosis in medical imaging, voxel-based segmentation has emerged as the primary methodology for automatic analysis of left ventricle (LV) function and morphology in cardiac magnetic resonance images (CMRI). In standard clinical practice, simultaneous multi-slice 2D cine short-axis MR imaging is performed under multiple breath-holds resulting in highly anisotropic 3D images. Furthermore, sparse-view CMRI often lacks whole heart coverage caused by large slice thickness and often suffers from inter-slice misalignment induced by respiratory motion. Therefore, these volumes only provide limited information about the true 3D cardiac anatomy which may hamper highly accurate assessment of functional and anatomical abnormalities. To address this, we propose a method that learns a continuous implicit function representing 3D LV shapes by training an auto-decoder. For training, high-resolution segmentations from cardiac CT angiography are used. The ability of our approach to reconstruct and complete high-resolution shapes from manually or automatically obtained sparse-view cardiac shape information is evaluated by using paired high- and low-resolution CMRI LV segmentations. The results show that the reconstructed LV shapes have an unconstrained subvoxel resolution and appear smooth and plausible in through-plane direction. Furthermore, Bland-Altman analysis reveals that reconstructed high-resolution ventricle volumes are closer to the corresponding reference volumes than reference low-resolution volumes with bias of [limits of agreement] -3.51 [-18.87, 11.85] mL, and 12.96 [-10.01, 35.92] mL respectively. Finally, the results demonstrate that the proposed approach allows recovering missing shape information and can indirectly correct for limited motion-induced artifacts.
Collapse
Affiliation(s)
- Jörg Sander
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Center location University of Amsterdam, The Netherlands; Amsterdam Cardiovascular Sciences, Amsterdam, The Netherlands; Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands.
| | - Bob D de Vos
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Center location University of Amsterdam, The Netherlands
| | - Steffen Bruns
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Center location University of Amsterdam, The Netherlands
| | - Nils Planken
- Department of Radiology and Nuclear Medicine,Amsterdam University Medical Center location University of Amsterdam, Amsterdam, The Netherlands
| | - Max A Viergever
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Tim Leiner
- Department of Radiology, Mayo Clinic, Rochester, MN, USA
| | - Ivana Išgum
- Department of Biomedical Engineering and Physics, Amsterdam University Medical Center location University of Amsterdam, The Netherlands; Amsterdam Cardiovascular Sciences, Amsterdam, The Netherlands; Department of Radiology and Nuclear Medicine,Amsterdam University Medical Center location University of Amsterdam, Amsterdam, The Netherlands; Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
6
|
Dokuchaev A, Chumarnaya T, Bazhutina A, Khamzin S, Lebedeva V, Lyubimtseva T, Zubarev S, Lebedev D, Solovyova O. Combination of personalized computational modeling and machine learning for optimization of left ventricular pacing site in cardiac resynchronization therapy. Front Physiol 2023; 14:1162520. [PMID: 37497440 PMCID: PMC10367108 DOI: 10.3389/fphys.2023.1162520] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Accepted: 06/26/2023] [Indexed: 07/28/2023] Open
Abstract
Introduction: The 30-50% non-response rate to cardiac resynchronization therapy (CRT) calls for improved patient selection and optimized pacing lead placement. The study aimed to develop a novel technique using patient-specific cardiac models and machine learning (ML) to predict an optimal left ventricular (LV) pacing site (ML-PS) that maximizes the likelihood of LV ejection fraction (LVEF) improvement in a given CRT candidate. To validate the approach, we evaluated whether the distance DPS between the clinical LV pacing site (ref-PS) and ML-PS is associated with improved response rate and magnitude. Materials and methods: We reviewed retrospective data for 57 CRT recipients. A positive response was defined as a more than 10% LVEF improvement. Personalized models of ventricular activation and ECG were created from MRI and CT images. The characteristics of ventricular activation during intrinsic rhythm and biventricular (BiV) pacing with ref-PS were derived from the models and used in combination with clinical data to train supervised ML classifiers. The best logistic regression model classified CRT responders with a high accuracy of 0.77 (ROC AUC = 0.84). The LR classifier, model simulations and Bayesian optimization with Gaussian process regression were combined to identify an optimal ML-PS that maximizes the ML-score of CRT response over the LV surface in each patient. Results: The optimal ML-PS improved the ML-score by 17 ± 14% over the ref-PS. Twenty percent of the non-responders were reclassified as positive at ML-PS. Selection of positive patients with a max ML-score >0.5 demonstrated an improved clinical response rate. The distance DPS was shorter in the responders. The max ML-score and DPS were found to be strong predictors of CRT response (ROC AUC = 0.85). In the group with max ML-score > 0.5 and DPS< 30 mm, the response rate was 83% compared to 14% in the rest of the cohort. LVEF improvement in this group was higher than in the other patients (16 ± 8% vs. 7 ± 8%). Conclusion: A new technique combining clinical data, personalized heart modelling and supervised ML demonstrates the potential for use in clinical practice to assist in optimizing patient selection and predicting optimal LV pacing lead position in HF candidates for CRT.
Collapse
Affiliation(s)
- Arsenii Dokuchaev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
| | - Tatiana Chumarnaya
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Anastasia Bazhutina
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Svyatoslav Khamzin
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
| | | | - Tamara Lyubimtseva
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Stepan Zubarev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Dmitry Lebedev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Olga Solovyova
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| |
Collapse
|
7
|
Meiburg R, Rijks JHJ, Beela AS, Bressi E, Grieco D, Delhaas T, Luermans JGLM, Prinzen FW, Vernooy K, Lumens J. Comparison of novel ventricular pacing strategies using an electro-mechanical simulation platform. Europace 2023; 25:euad144. [PMID: 37306315 PMCID: PMC10259067 DOI: 10.1093/europace/euad144] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 05/06/2023] [Indexed: 06/13/2023] Open
Abstract
AIMS Focus of pacemaker therapy is shifting from right ventricular (RV) apex pacing (RVAP) and biventricular pacing (BiVP) to conduction system pacing. Direct comparison between the different pacing modalities and their consequences to cardiac pump function is difficult, due to the practical implications and confounding variables. Computational modelling and simulation provide the opportunity to compare electrical, mechanical, and haemodynamic consequences in the same virtual heart. METHODS AND RESULTS Using the same single cardiac geometry, electrical activation maps following the different pacing strategies were calculated using an Eikonal model on a three-dimensional geometry, which were then used as input for a lumped mechanical and haemodynamic model (CircAdapt). We then compared simulated strain, regional myocardial work, and haemodynamic function for each pacing strategy. Selective His-bundle pacing (HBP) best replicated physiological electrical activation and led to the most homogeneous mechanical behaviour. Selective left bundle branch (LBB) pacing led to good left ventricular (LV) function but significantly increased RV load. RV activation times were reduced in non-selective LBB pacing (nsLBBP), reducing RV load but increasing heterogeneity in LV contraction. LV septal pacing led to a slower LV and more heterogeneous LV activation than nsLBBP, while RV activation was similar. BiVP led to a synchronous LV-RV, but resulted in a heterogeneous contraction. RVAP led to the slowest and most heterogeneous contraction. Haemodynamic differences were small compared to differences in local wall behaviour. CONCLUSION Using a computational modelling framework, we investigated the mechanical and haemodynamic outcome of the prevailing pacing strategies in hearts with normal electrical and mechanical function. For this class of patients, nsLBBP was the best compromise between LV and RV function if HBP is not possible.
Collapse
Affiliation(s)
- Roel Meiburg
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Universiteitssingel 40, 6200 MD, Maastricht, The Netherlands
| | - Jesse H J Rijks
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+ (MUMC+), Maastricht, The Netherlands
| | - Ahmed S Beela
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Universiteitssingel 40, 6200 MD, Maastricht, The Netherlands
- Department of Cardiovascular Diseases, Faculty of Medicine, Suez Canal University, Ismailia, Egypt
| | - Edoardo Bressi
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+ (MUMC+), Maastricht, The Netherlands
- Department of Cardiovascular Sciences, Policlinico Casilino of Rome, Rome, Italy
| | - Domenico Grieco
- Department of Cardiovascular Sciences, Policlinico Casilino of Rome, Rome, Italy
| | - Tammo Delhaas
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Universiteitssingel 40, 6200 MD, Maastricht, The Netherlands
| | - Justin G LM Luermans
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+ (MUMC+), Maastricht, The Netherlands
- Department of Cardiology, Radboud University Medical Centre (Radboudumc), Nijmegen, The Netherlands
| | - Frits W Prinzen
- Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
| | - Kevin Vernooy
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center+ (MUMC+), Maastricht, The Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Universiteitssingel 40, 6200 MD, Maastricht, The Netherlands
| |
Collapse
|
8
|
Viola F, Del Corso G, De Paulis R, Verzicco R. GPU accelerated digital twins of the human heart open new routes for cardiovascular research. Sci Rep 2023; 13:8230. [PMID: 37217483 DOI: 10.1038/s41598-023-34098-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Accepted: 04/24/2023] [Indexed: 05/24/2023] Open
Abstract
The recruitment of patients for rare or complex cardiovascular diseases is a bottleneck for clinical trials and digital twins of the human heart have recently been proposed as a viable alternative. In this paper we present an unprecedented cardiovascular computer model which, relying on the latest GPU-acceleration technologies, replicates the full multi-physics dynamics of the human heart within a few hours per heartbeat. This opens the way to extensive simulation campaigns to study the response of synthetic cohorts of patients to cardiovascular disorders, novel prosthetic devices or surgical procedures. As a proof-of-concept we show the results obtained for left bundle branch block disorder and the subsequent cardiac resynchronization obtained by pacemaker implantation. The in-silico results closely match those obtained in clinical practice, confirming the reliability of the method. This innovative approach makes possible a systematic use of digital twins in cardiovascular research, thus reducing the need of real patients with their economical and ethical implications. This study is a major step towards in-silico clinical trials in the era of digital medicine.
Collapse
Affiliation(s)
- Francesco Viola
- Gran Sasso Science Institute (GSSI), L'Aquila, Italy
- INFN-Laboratori Nazionali del Gran Sasso, Assergi (AQ), Italy
| | - Giulio Del Corso
- Gran Sasso Science Institute (GSSI), L'Aquila, Italy
- Institute of Information Science and Technologies A. Faedo, CNR, Pisa, Italy
| | - Ruggero De Paulis
- European Hospital, Rome, Italy
- UniCamillus International University of Health Sciences, Rome, Italy
| | - Roberto Verzicco
- Gran Sasso Science Institute (GSSI), L'Aquila, Italy.
- University of Rome Tor Vergata, Rome, Italy.
- POF Group, University of Twente, Enschede, The Netherlands.
| |
Collapse
|
9
|
Salvador M, Regazzoni F, Dede' L, Quarteroni A. Fast and robust parameter estimation with uncertainty quantification for the cardiac function. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 231:107402. [PMID: 36773593 DOI: 10.1016/j.cmpb.2023.107402] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Revised: 01/31/2023] [Accepted: 02/01/2023] [Indexed: 06/18/2023]
Abstract
BACKGROUND AND OBJECTIVES Parameter estimation and uncertainty quantification are crucial in computational cardiology, as they enable the construction of digital twins that faithfully replicate the behavior of physical patients. Many model parameters regarding cardiac electromechanics and cardiovascular hemodynamics need to be robustly fitted by starting from a few, possibly non-invasive, noisy observations. Moreover, short execution times and a small amount of computational resources are required for the effective clinical translation. METHODS In the framework of Bayesian statistics, we combine Maximum a Posteriori estimation and Hamiltonian Monte Carlo to find an approximation of model parameters and their posterior distributions. Fast simulations and minimal memory requirements are achieved by using an accurate and geometry-specific Artificial Neural Network surrogate model for the cardiac function, matrix-free methods, automatic differentiation and automatic vectorization. Furthermore, we account for the surrogate modeling error and measurement error. RESULTS We perform three different in silico test cases, ranging from the ventricular function to the entire cardiocirculatory system, involving whole-heart mechanics, arterial and venous hemodynamics. By employing a single central processing unit on a standard laptop, we attain highly accurate estimations for all model parameters in short computational times. Furthermore, we obtain posterior distributions that contain the true values inside the 90% credibility regions. CONCLUSIONS Many model parameters regarding the entire cardiovascular system can be fastly and robustly identified with minimal hardware requirements. This can be achieved when a small amount of non-invasive data is available and when high levels of signal-to-noise ratio are present in the quantities of interest. With these features, our approach meets the requirements for clinical exploitation, while being compliant with Green Computing practices.
Collapse
Affiliation(s)
- Matteo Salvador
- MOX-Dipartimento di Matematica, P.zza Leonardo da Vinci 32, Milan, 20133, Italy.
| | - Francesco Regazzoni
- MOX-Dipartimento di Matematica, P.zza Leonardo da Vinci 32, Milan, 20133, Italy
| | - Luca Dede'
- MOX-Dipartimento di Matematica, P.zza Leonardo da Vinci 32, Milan, 20133, Italy
| | - Alfio Quarteroni
- MOX-Dipartimento di Matematica, P.zza Leonardo da Vinci 32, Milan, 20133, Italy; Institute of Mathematics, École Polytechnique Fédérale de Lausanne, Av. Piccard, Lausanne, 1015, Switzerland
| |
Collapse
|
10
|
Galli E, Galand V, Le Rolle V, Taconne M, Wazzan AA, Hernandez A, Leclercq C, Donal E. The saga of dyssynchrony imaging: Are we getting to the point. Front Cardiovasc Med 2023; 10:1111538. [PMID: 37063957 PMCID: PMC10103462 DOI: 10.3389/fcvm.2023.1111538] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 02/27/2023] [Indexed: 04/03/2023] Open
Abstract
Cardiac resynchronisation therapy (CRT) has an established role in the management of patients with heart failure, reduced left ventricular ejection fraction (LVEF < 35%) and widened QRS (>130 msec). Despite the complex pathophysiology of left ventricular (LV) dyssynchrony and the increasing evidence supporting the identification of specific electromechanical substrates that are associated with a higher probability of CRT response, the assessment of LVEF is the only imaging-derived parameter used for the selection of CRT candidates.This review aims to (1) provide an overview of the evolution of cardiac imaging for the assessment of LV dyssynchrony and its role in the selection of patients undergoing CRT; (2) highlight the main pitfalls and advantages of the application of cardiac imaging for the assessment of LV dyssynchrony; (3) provide some perspectives for clinical application and future research in this field.Conclusionthe road for a more individualized approach to resynchronization therapy delivery is open and imaging might provide important input beyond the assessment of LVEF.
Collapse
|
11
|
Galappaththige S, Gray RA, Costa CM, Niederer S, Pathmanathan P. Credibility assessment of patient-specific computational modeling using patient-specific cardiac modeling as an exemplar. PLoS Comput Biol 2022; 18:e1010541. [PMID: 36215228 PMCID: PMC9550052 DOI: 10.1371/journal.pcbi.1010541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Accepted: 09/02/2022] [Indexed: 11/07/2022] Open
Abstract
Reliable and robust simulation of individual patients using patient-specific models (PSMs) is one of the next frontiers for modeling and simulation (M&S) in healthcare. PSMs, which form the basis of digital twins, can be employed as clinical tools to, for example, assess disease state, predict response to therapy, or optimize therapy. They may also be used to construct virtual cohorts of patients, for in silico evaluation of medical product safety and/or performance. Methods and frameworks have recently been proposed for evaluating the credibility of M&S in healthcare applications. However, such efforts have generally been motivated by models of medical devices or generic patient models; how best to evaluate the credibility of PSMs has largely been unexplored. The aim of this paper is to understand and demonstrate the credibility assessment process for PSMs using patient-specific cardiac electrophysiological (EP) modeling as an exemplar. We first review approaches used to generate cardiac PSMs and consider how verification, validation, and uncertainty quantification (VVUQ) apply to cardiac PSMs. Next, we execute two simulation studies using a publicly available virtual cohort of 24 patient-specific ventricular models, the first a multi-patient verification study, the second investigating the impact of uncertainty in personalized and non-personalized inputs in a virtual cohort. We then use the findings from our analyses to identify how important characteristics of PSMs can be considered when assessing credibility with the approach of the ASME V&V40 Standard, accounting for PSM concepts such as inter- and intra-user variability, multi-patient and “every-patient” error estimation, uncertainty quantification in personalized vs non-personalized inputs, clinical validation, and others. The results of this paper will be useful to developers of cardiac and other medical image based PSMs, when assessing PSM credibility. Patient-specific models are computational models that have been personalized using data from a patient. After decades of research, recent computational, data science and healthcare advances have opened the door to the fulfilment of the enormous potential of such models, from truly personalized medicine to efficient and cost-effective testing of new medical products. However, reliability (credibility) of patient-specific models is key to their success, and there are currently no general guidelines for evaluating credibility of patient-specific models. Here, we consider how frameworks and model evaluation activities that have been developed for generic (not patient-specific) computational models, can be extended to patient specific models. We achieve this through a detailed analysis of the activities required to evaluate cardiac electrophysiological models, chosen as an exemplar field due to its maturity and the complexity of such models. This is the first paper on the topic of reliability of patient-specific models and will help pave the way to reliable and trusted patient-specific modeling across healthcare applications.
Collapse
Affiliation(s)
- Suran Galappaththige
- Center for Devices and Radiological Health, US Food and Drug Administration, Silver Spring, Maryland, United States of America
| | - Richard A. Gray
- Center for Devices and Radiological Health, US Food and Drug Administration, Silver Spring, Maryland, United States of America
| | - Caroline Mendonca Costa
- School of Biomedical Engineering & Imaging Sciences, King’s College London, London, United Kingdom
| | - Steven Niederer
- School of Biomedical Engineering & Imaging Sciences, King’s College London, London, United Kingdom
| | - Pras Pathmanathan
- Center for Devices and Radiological Health, US Food and Drug Administration, Silver Spring, Maryland, United States of America
- * E-mail:
| |
Collapse
|
12
|
Marlevi D, Mariscal-Harana J, Burris NS, Sotelo J, Ruijsink B, Hadjicharalambous M, Asner L, Sammut E, Chabiniok R, Uribe S, Winter R, Lamata P, Alastruey J, Nordsletten D. Altered Aortic Hemodynamics and Relative Pressure in Patients with Dilated Cardiomyopathy. J Cardiovasc Transl Res 2022; 15:692-707. [PMID: 34882286 PMCID: PMC9622552 DOI: 10.1007/s12265-021-10181-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Accepted: 10/20/2021] [Indexed: 12/05/2022]
Abstract
Ventricular-vascular interaction is central in the adaptation to cardiovascular disease. However, cardiomyopathy patients are predominantly monitored using cardiac biomarkers. The aim of this study is therefore to explore aortic function in dilated cardiomyopathy (DCM). Fourteen idiopathic DCM patients and 16 controls underwent cardiac magnetic resonance imaging, with aortic relative pressure derived using physics-based image processing and a virtual cohort utilized to assess the impact of cardiovascular properties on aortic behaviour. Subjects with reduced left ventricular systolic function had significantly reduced aortic relative pressure, increased aortic stiffness, and significantly delayed time-to-pressure peak duration. From the virtual cohort, aortic stiffness and aortic volumetric size were identified as key determinants of aortic relative pressure. As such, this study shows how advanced flow imaging and aortic hemodynamic evaluation could provide novel insights into the manifestation of DCM, with signs of both altered aortic structure and function derived in DCM using our proposed imaging protocol.
Collapse
Affiliation(s)
- David Marlevi
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, USA
- Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, Huddinge, Sweden
- Department of Clinical Sciences, Karolinska Institutet, Danderyd, Sweden
| | - Jorge Mariscal-Harana
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | | | - Julio Sotelo
- School of Biomedical Engineering, Universidad de Valparaíso, Valparaíso, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus in Cardiovascular Magnetic Resonance, Santiago, Cardio MR, Chile
| | - Bram Ruijsink
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Myrianthi Hadjicharalambous
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| | - Liya Asner
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Eva Sammut
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Faculty of Health Science, Bristol Heart Institute and Translational Biomedical Research Centre, University of Bristol, Bristol, UK
| | - Radomir Chabiniok
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Inria, Palaiseau, France
- LMS, Ecole Polytechnique, CNRS, Institut Polytechnique de Paris, Paris, France
- Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, , Prague, Czech Republic
| | - Sergio Uribe
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus in Cardiovascular Magnetic Resonance, Santiago, Cardio MR, Chile
- Department of Radiology, School of Medicine, Pontifica Universidad Católica de Chile, Santiago, Chile
| | - Reidar Winter
- Department of Clinical Sciences, Karolinska Institutet, Danderyd, Sweden
| | - Pablo Lamata
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Jordi Alastruey
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- World-Class Research Center "Digital Biodesign and Personlized Healthcare", Sechenov University, Moscow, Russia
| | - David Nordsletten
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
- Department of Cardiac Surgery and Biomedical Engineering, University of Michigan, Plymouth Rd, Ann Arbor, MI, 48109, USA.
| |
Collapse
|
13
|
Marx L, Niestrawska JA, Gsell MA, Caforio F, Plank G, Augustin CM. Robust and efficient fixed-point algorithm for the inverse elastostatic problem to identify myocardial passive material parameters and the unloaded reference configuration. JOURNAL OF COMPUTATIONAL PHYSICS 2022; 463:111266. [PMID: 35662800 PMCID: PMC7612790 DOI: 10.1016/j.jcp.2022.111266] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Image-based computational models of the heart represent a powerful tool to shed new light on the mechanisms underlying physiological and pathological conditions in cardiac function and to improve diagnosis and therapy planning. However, in order to enable the clinical translation of such models, it is crucial to develop personalized models that are able to reproduce the physiological reality of a given patient. There have been numerous contributions in experimental and computational biomechanics to characterize the passive behavior of the myocardium. However, most of these studies suffer from severe limitations and are not applicable to high-resolution geometries. In this work, we present a novel methodology to perform an automated identification of in vivo properties of passive cardiac biomechanics. The highly-efficient algorithm fits material parameters against the shape of a patient-specific approximation of the end-diastolic pressure-volume relation (EDPVR). Simultaneously, an unloaded reference configuration is generated, where a novel line search strategy to improve convergence and robustness is implemented. Only clinical image data or previously generated meshes at one time point during diastole and one measured data point of the EDPVR are required as an input. The proposed method can be straightforwardly coupled to existing finite element (FE) software packages and is applicable to different constitutive laws and FE formulations. Sensitivity analysis demonstrates that the algorithm is robust with respect to initial input parameters.
Collapse
Affiliation(s)
- Laura Marx
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Justyna A. Niestrawska
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Matthias A.F. Gsell
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Federica Caforio
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
- Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - Gernot Plank
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Christoph M. Augustin
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
- Corresponding author at: Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Neue Stiftingtalstrasse 6/D04, 8010 Graz, Austria. (C.M.Augustin)
| |
Collapse
|
14
|
Mendonca Costa C, Gemmell P, Elliott MK, Whitaker J, Campos FO, Strocchi M, Neic A, Gillette K, Vigmond E, Plank G, Razavi R, O'Neill M, Rinaldi CA, Bishop MJ. Determining anatomical and electrophysiological detail requirements for computational ventricular models of porcine myocardial infarction. Comput Biol Med 2022; 141:105061. [PMID: 34915331 PMCID: PMC8819160 DOI: 10.1016/j.compbiomed.2021.105061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Revised: 11/04/2021] [Accepted: 11/20/2021] [Indexed: 12/01/2022]
Abstract
BACKGROUND Computational models of the heart built from cardiac MRI and electrophysiology (EP) data have shown promise for predicting the risk of and ablation targets for myocardial infarction (MI) related ventricular tachycardia (VT), as well as to predict paced activation sequences in heart failure patients. However, most recent studies have relied on low resolution imaging data and little or no EP personalisation, which may affect the accuracy of model-based predictions. OBJECTIVE To investigate the impact of model anatomy, MI scar morphology, and EP personalisation strategies on paced activation sequences and VT inducibility to determine the level of detail required to make accurate model-based predictions. METHODS Imaging and EP data were acquired from a cohort of six pigs with experimentally induced MI. Computational models of ventricular anatomy, incorporating MI scar, were constructed including bi-ventricular or left ventricular (LV) only anatomy, and MI scar morphology with varying detail. Tissue conductivities and action potential duration (APD) were fitted to 12-lead ECG data using the QRS duration and the QT interval, respectively, in addition to corresponding literature parameters. Paced activation sequences and VT induction were simulated. Simulated paced activation and VT inducibility were compared between models and against experimental data. RESULTS Simulations predict that the level of model anatomical detail has little effect on simulated paced activation, with all model predictions comparing closely with invasive EP measurements. However, detailed scar morphology from high-resolution images, bi-ventricular anatomy, and personalized tissue conductivities are required to predict experimental VT outcome. CONCLUSION This study provides clear guidance for model generation based on clinical data. While a representing high level of anatomical and scar detail will require high-resolution image acquisition, EP personalisation based on 12-lead ECG can be readily incorporated into modelling pipelines, as such data is widely available.
Collapse
Affiliation(s)
- Caroline Mendonca Costa
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK.
| | - Philip Gemmell
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Mark K Elliott
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - John Whitaker
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Fernando O Campos
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Marina Strocchi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | | | - Karli Gillette
- Gottfried Schatz Research Center, Biophysics, Medical University of Graz, Austria; Medical University of Graz, Austria and BioTechMed, Graz, Austria
| | - Edward Vigmond
- Institut de Rythmologie et de modélisation cardiaque (LIRYC), University of Bordeaux, France
| | - Gernot Plank
- Medical University of Graz, Austria and BioTechMed, Graz, Austria
| | - Reza Razavi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Mark O'Neill
- Department of Cardiology, Guy's and St Thomas' Hospital, London, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK; Department of Cardiology, Guy's and St Thomas' Hospital, London, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| |
Collapse
|
15
|
Fan L, Choy JS, Raissi F, Kassab GS, Lee LC. Optimization of cardiac resynchronization therapy based on a cardiac electromechanics-perfusion computational model. Comput Biol Med 2022; 141:105050. [PMID: 34823858 PMCID: PMC8810745 DOI: 10.1016/j.compbiomed.2021.105050] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 11/10/2021] [Accepted: 11/15/2021] [Indexed: 02/03/2023]
Abstract
Cardiac resynchronization therapy (CRT) is an established treatment for left bundle branch block (LBBB) resulting in mechanical dyssynchrony. Approximately 1/3 of patients with CRT, however, are non-responders. To understand factors affecting CRT response, an electromechanics-perfusion computational model based on animal-specific left ventricular (LV) geometry and coronary vascular networks located in the septum and LV free wall is developed. The model considers contractility-flow and preload-activation time relationships, and is calibrated to simultaneously match the experimental measurements in terms of the LV pressure, volume waveforms and total coronary flow in the left anterior descending and left circumflex territories from 2 swine models under right atrium and right ventricular pacing. The model is then applied to investigate the responses of CRT indexed by peak LV pressure and (dP/dt)max at multiple pacing sites with different degrees of perfusion in the LV free wall. Without the presence of ischemia, the model predicts that basal-lateral endocardial region is the optimal pacing site that can best improve (dP/dt)max by 20%, and is associated with the shortest activation time. In the presence of ischemia, a non-ischemic region becomes the optimal pacing site when coronary flow in the ischemic region fell below 30% of its original value. Pacing at the ischemic region produces little response at that perfusion level. The optimal pacing site is associated with one that optimizes the LV activation time. These findings suggest that CRT response is affected by both pacing site and coronary perfusion, which may have clinical implication in improving CRT responder rates.
Collapse
Affiliation(s)
- Lei Fan
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA.
| | - Jenny S Choy
- California Medical Innovations Institute, San Diego, CA, USA
| | - Farshad Raissi
- Department of Medicine, University of California, San Diego, La Jolla, CA, USA
| | | | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
16
|
Khamzin S, Dokuchaev A, Bazhutina A, Chumarnaya T, Zubarev S, Lyubimtseva T, Lebedeva V, Lebedev D, Gurev V, Solovyova O. Machine Learning Prediction of Cardiac Resynchronisation Therapy Response From Combination of Clinical and Model-Driven Data. Front Physiol 2022; 12:753282. [PMID: 34970154 PMCID: PMC8712879 DOI: 10.3389/fphys.2021.753282] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Accepted: 11/22/2021] [Indexed: 11/17/2022] Open
Abstract
Background: Up to 30–50% of chronic heart failure patients who underwent cardiac resynchronization therapy (CRT) do not respond to the treatment. Therefore, patient stratification for CRT and optimization of CRT device settings remain a challenge. Objective: The main goal of our study is to develop a predictive model of CRT outcome using a combination of clinical data recorded in patients before CRT and simulations of the response to biventricular (BiV) pacing in personalized computational models of the cardiac electrophysiology. Materials and Methods: Retrospective data from 57 patients who underwent CRT device implantation was utilized. Positive response to CRT was defined by a 10% increase in the left ventricular ejection fraction in a year after implantation. For each patient, an anatomical model of the heart and torso was reconstructed from MRI and CT images and tailored to ECG recorded in the participant. The models were used to compute ventricular activation time, ECG duration and electrical dyssynchrony indices during intrinsic rhythm and BiV pacing from the sites of implanted leads. For building a predictive model of CRT response, we used clinical data recorded before CRT device implantation together with model-derived biomarkers of ventricular excitation in the left bundle branch block mode of activation and under BiV stimulation. Several Machine Learning (ML) classifiers and feature selection algorithms were tested on the hybrid dataset, and the quality of predictors was assessed using the area under receiver operating curve (ROC AUC). The classifiers on the hybrid data were compared with ML models built on clinical data only. Results: The best ML classifier utilizing a hybrid set of clinical and model-driven data demonstrated ROC AUC of 0.82, an accuracy of 0.82, sensitivity of 0.85, and specificity of 0.78, improving quality over that of ML predictors built on clinical data from much larger datasets by more than 0.1. Distance from the LV pacing site to the post-infarction zone and ventricular activation characteristics under BiV pacing were shown as the most relevant model-driven features for CRT response classification. Conclusion: Our results suggest that combination of clinical and model-driven data increases the accuracy of classification models for CRT outcomes.
Collapse
Affiliation(s)
- Svyatoslav Khamzin
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Arsenii Dokuchaev
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Anastasia Bazhutina
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia.,Ural Federal University, Yekaterinburg, Russia
| | - Tatiana Chumarnaya
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Stepan Zubarev
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | | | | | - Dmitry Lebedev
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | | | - Olga Solovyova
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia.,Ural Federal University, Yekaterinburg, Russia
| |
Collapse
|
17
|
Prediction of Ventricular Mechanics After Pulmonary Valve Replacement in Tetralogy of Fallot by Biomechanical Modeling: A Step Towards Precision Healthcare. Ann Biomed Eng 2021; 49:3339-3348. [PMID: 34853921 DOI: 10.1007/s10439-021-02895-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 11/12/2021] [Indexed: 10/19/2022]
Abstract
Clinical indicators of heart function are often limited in their ability to accurately evaluate the current mechanical state of the myocardium. Biomechanical modeling has been shown to be a promising tool in addition to clinical indicators. By providing a patient-specific measure of myocardial active stress (contractility), biomechanical modeling can enhance the precision of the description of patient's pathophysiology at any given point in time. In this work we aim to explore the ability of biomechanical modeling to predict the response of ventricular mechanics to the progressively decreasing afterload in repaired tetralogy of Fallot (rTOF) patients undergoing pulmonary valve replacement (PVR) for significant residual right ventricular outflow tract obstruction (RVOTO). We used 19 patient-specific models of patients with rTOF prior to pulmonary valve replacement (PVR), denoted as PSMpre, and patient-specific models of the same patients created post-PVR (PSMpost)-both created in our previous published work. Using the PSMpre and assuming cessation of the pulmonary regurgitation and a progressive decrease of RVOT resistance, we built relationships between the contractility and RVOT resistance post-PVR. The predictive value of such in silico obtained relationships were tested against the PSMpost, i.e. the models created from the actual post-PVR datasets. Our results show a linear 1-dimensional relationship between the in silico predicted contractility post-PVR and the RVOT resistance. The predicted contractility was close to the contractility in the PSMpost model with a mean (± SD) difference of 6.5 (± 3.0)%. The relationships between the contractility predicted by in silico PVR vs. RVOT resistance have a potential to inform clinicians about hypothetical mechanical response of the ventricle based on the degree of pre-operative RVOTO.
Collapse
|
18
|
Augustin CM, Gsell MA, Karabelas E, Willemen E, Prinzen FW, Lumens J, Vigmond EJ, Plank G. A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2021; 386:114092. [PMID: 34630765 PMCID: PMC7611781 DOI: 10.1016/j.cma.2021.114092] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.
Collapse
Affiliation(s)
- Christoph M. Augustin
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Matthias A.F. Gsell
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Elias Karabelas
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Erik Willemen
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, Netherlands
| | - Frits W. Prinzen
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, Netherlands
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
| | - Gernot Plank
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
- Correspondence to: Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Neue Stiftingtalstrasse 6/IV, Graz 8010, Austria. (G. Plank)
| |
Collapse
|
19
|
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
|
20
|
Miller R, Kerfoot E, Mauger C, Ismail TF, Young AA, Nordsletten DA. An Implementation of Patient-Specific Biventricular Mechanics Simulations With a Deep Learning and Computational Pipeline. Front Physiol 2021; 12:716597. [PMID: 34603077 PMCID: PMC8481785 DOI: 10.3389/fphys.2021.716597] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 08/06/2021] [Indexed: 02/04/2023] Open
Abstract
Parameterised patient-specific models of the heart enable quantitative analysis of cardiac function as well as estimation of regional stress and intrinsic tissue stiffness. However, the development of personalised models and subsequent simulations have often required lengthy manual setup, from image labelling through to generating the finite element model and assigning boundary conditions. Recently, rapid patient-specific finite element modelling has been made possible through the use of machine learning techniques. In this paper, utilising multiple neural networks for image labelling and detection of valve landmarks, together with streamlined data integration, a pipeline for generating patient-specific biventricular models is applied to clinically-acquired data from a diverse cohort of individuals, including hypertrophic and dilated cardiomyopathy patients and healthy volunteers. Valve motion from tracked landmarks as well as cavity volumes measured from labelled images are used to drive realistic motion and estimate passive tissue stiffness values. The neural networks are shown to accurately label cardiac regions and features for these diverse morphologies. Furthermore, differences in global intrinsic parameters, such as tissue anisotropy and normalised active tension, between groups illustrate respective underlying changes in tissue composition and/or structure as a result of pathology. This study shows the successful application of a generic pipeline for biventricular modelling, incorporating artificial intelligence solutions, within a diverse cohort.
Collapse
Affiliation(s)
- Renee Miller
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Eric Kerfoot
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Charlène Mauger
- Auckland MR Research Group, University of Auckland, Auckland, New Zealand
| | - Tevfik F. Ismail
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Alistair A. Young
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- Auckland MR Research Group, University of Auckland, Auckland, New Zealand
| | - David A. Nordsletten
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- Department of Biomedical Engineering and Cardiac Surgery, University of Michigan, Ann Arbor, MI, United States
| |
Collapse
|
21
|
Sermesant M, Delingette H, Cochet H, Jaïs P, Ayache N. Applications of artificial intelligence in cardiovascular imaging. Nat Rev Cardiol 2021; 18:600-609. [PMID: 33712806 DOI: 10.1038/s41569-021-00527-2] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 02/08/2021] [Indexed: 01/31/2023]
Abstract
Research into artificial intelligence (AI) has made tremendous progress over the past decade. In particular, the AI-powered analysis of images and signals has reached human-level performance in many applications owing to the efficiency of modern machine learning methods, in particular deep learning using convolutional neural networks. Research into the application of AI to medical imaging is now very active, especially in the field of cardiovascular imaging because of the challenges associated with acquiring and analysing images of this dynamic organ. In this Review, we discuss the clinical questions in cardiovascular imaging that AI can be used to address and the principal methodological AI approaches that have been developed to solve the related image analysis problems. Some approaches are purely data-driven and rely mainly on statistical associations, whereas others integrate anatomical and physiological information through additional statistical, geometric and biophysical models of the human heart. In a structured manner, we provide representative examples of each of these approaches, with particular attention to the underlying computational imaging challenges. Finally, we discuss the remaining limitations of AI approaches in cardiovascular imaging (such as generalizability and explainability) and how they can be overcome.
Collapse
Affiliation(s)
| | | | - Hubert Cochet
- IHU Liryc, CHU Bordeaux, Université Bordeaux, Inserm 1045, Pessac, France
| | - Pierre Jaïs
- IHU Liryc, CHU Bordeaux, Université Bordeaux, Inserm 1045, Pessac, France
| | | |
Collapse
|
22
|
Dell'Era G, Gravellone M, Scacchi S, Franzone PC, Pavarino LF, Boggio E, Prenna E, De Vecchi F, Occhetta E, Devecchi C, Patti G. A clinical-in silico study on the effectiveness of multipoint bicathodic and cathodic-anodal pacing in cardiac resynchronization therapy. Comput Biol Med 2021; 136:104661. [PMID: 34332350 DOI: 10.1016/j.compbiomed.2021.104661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 07/16/2021] [Accepted: 07/16/2021] [Indexed: 11/18/2022]
Abstract
Up to one-third of patients undergoing cardiac resynchronization therapy (CRT) are nonresponders. Multipoint bicathodic and cathodic-anodal left ventricle (LV) stimulations could overcome this clinical challenge, but their effectiveness remains controversial. Here we evaluate the performance of such stimulations through both in vivo and in silico experiments, the latter based on computer electromechanical modeling. Seven patients, all candidates for CRT, received a quadripolar LV lead. Four stimulations were tested: right ventricular (RVS); conventional single point biventricular (S-BS); multipoint biventricular bicathodic (CC-BS) and multipoint biventricular cathodic-anodal (CA-BS). The following parameters were processed: QRS duration; maximal time derivative of arterial pressure (dPdtmax); systolic arterial pressure (Psys); and stroke volume (SV). Echocardiographic data of each patient were then obtained to create an LV geometric model. Numerical simulations were based on a strongly coupled Bidomain electromechanical coupling model. Considering the in vivo parameters, when comparing S-BS to RVS, there was no significant decrease in SV (from 45 ± 11 to 44 ± 20 ml) and 6% and 4% increases of dPdtmax and Psys, respectively. Focusing on in silico parameters, with respect to RVS, S-BS exhibited a significant increase of SV, dPdtmax and Psys. Neither the in vivo nor in silico results showed any significant hemodynamic and electrical difference among S-BS, CC-BS and CA-BS configurations. These results show that CC-BS and CA-BS yield a comparable CRT performance, but they do not always yield improvement in terms of hemodynamic parameters with respect to S-BS. The computational results confirmed the in vivo observations, thus providing theoretical support to the clinical experiments.
Collapse
Affiliation(s)
- G Dell'Era
- Cardiologia 1, Azienda Ospedaliera Universitaria "Maggiore Della Carità", Novara, Italy
| | - M Gravellone
- Divisione di Cardiologia, Ospedale Degli Infermi, Biella, Italy
| | - S Scacchi
- Dipartimento di Matematica, Università Degli Studi di Milano, Via Saldini 50, 20133, Milano, Italy.
| | - P Colli Franzone
- Dipartimento di Matematica, Università Degli Studi di Pavia, Via Ferrata 1, 27100, Pavia, Italy
| | - L F Pavarino
- Dipartimento di Matematica, Università Degli Studi di Pavia, Via Ferrata 1, 27100, Pavia, Italy
| | - E Boggio
- Divisione di Cardiologia, Ospedale Degli Infermi, Biella, Italy
| | - E Prenna
- Cardiologia 1, Azienda Ospedaliera Universitaria "Maggiore Della Carità", Novara, Italy
| | - F De Vecchi
- Divisione di Cardiologia, Ospedale Sant'Andrea, Vercelli, Italy
| | - E Occhetta
- Divisione di Cardiologia, Ospedale Sant'Andrea, Vercelli, Italy
| | - C Devecchi
- Divisione di Cardiologia, Ospedale Sant'Andrea, Vercelli, Italy
| | - G Patti
- Cardiologia 1, Azienda Ospedaliera Universitaria "Maggiore Della Carità", Novara, Italy; Dipartimento di Medicina Traslazionale, Università Del Piemonte Orientale, Novara, Italy
| |
Collapse
|
23
|
Coveney S, Corrado C, Oakley JE, Wilkinson RD, Niederer SA, Clayton RH. Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators. Front Physiol 2021; 12:693015. [PMID: 34366883 PMCID: PMC8339909 DOI: 10.3389/fphys.2021.693015] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Accepted: 06/28/2021] [Indexed: 11/13/2022] Open
Abstract
Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.
Collapse
Affiliation(s)
- Sam Coveney
- Insigneo Institute for In-Silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, United Kingdom
| | - Cesare Corrado
- Division of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Jeremy E. Oakley
- School of Mathematics and Statistics, University of Sheffield, Sheffield, United Kingdom
| | - Richard D. Wilkinson
- School of Mathematical Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Steven A. Niederer
- Division of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Richard H. Clayton
- Insigneo Institute for In-Silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, United Kingdom
| |
Collapse
|
24
|
Hadjicharalambous M, Stoeck CT, Weisskopf M, Cesarovic N, Ioannou E, Vavourakis V, Nordsletten DA. Investigating the reference domain influence in personalised models of cardiac mechanics : Effect of unloaded geometry on cardiac biomechanics. Biomech Model Mechanobiol 2021; 20:1579-1597. [PMID: 34047891 DOI: 10.1007/s10237-021-01464-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Accepted: 05/03/2021] [Indexed: 01/23/2023]
Abstract
A major concern in personalised models of heart mechanics is the unknown zero-pressure domain, a prerequisite for accurately predicting cardiac biomechanics. As the reference configuration cannot be captured by clinical data, studies often employ in-vivo frames which are unlikely to correspond to unloaded geometries. Alternatively, zero-pressure domain is approximated through inverse methodologies, which, however, entail assumptions pertaining to boundary conditions and material parameters. Both approaches are likely to introduce biases in estimated biomechanical properties; nevertheless, quantification of these effects is unattainable without ground-truth data. In this work, we assess the unloaded state influence on model-derived biomechanics, by employing an in-silico modelling framework relying on experimental data on porcine hearts. In-vivo images are used for model personalisation, while in-situ experiments provide a reliable approximation of the reference domain, creating a unique opportunity for a validation study. Personalised whole-cycle cardiac models are developed which employ different reference domains (image-derived, inversely estimated) and are compared against ground-truth model outcomes. Simulations are conducted with varying boundary conditions, to investigate the effect of data-derived constraints on model accuracy. Attention is given to modelling the influence of the ribcage on the epicardium, due to its close proximity to the heart in the porcine anatomy. Our results find merit in both approaches for dealing with the unknown reference domain, but also demonstrate differences in estimated biomechanical quantities such as material parameters, strains and stresses. Notably, they highlight the importance of a boundary condition accounting for the constraining influence of the ribcage, in forward and inverse biomechanical models.
Collapse
Affiliation(s)
| | - Christian T Stoeck
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Miriam Weisskopf
- Center for Surgical Research, University Hospital Zurich, University of Zurich, Zurich, Switzerland
| | - Nikola Cesarovic
- Center for Surgical Research, University Hospital Zurich, University of Zurich, Zurich, Switzerland.,Translational Cardiovascular Technologies, Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland.,Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany
| | - Eleftherios Ioannou
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| | - Vasileios Vavourakis
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus.,Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - David A Nordsletten
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.,Department of Biomedical Engineering and Cardiac Surgery, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
25
|
Building Models of Patient-Specific Anatomy and Scar Morphology from Clinical MRI Data. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11663-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|
26
|
Albatat M, Arevalo H, Bergsland J, Strøm V, Balasingham I, Odland HH. Optimal pacing sites in cardiac resynchronization by left ventricular activation front analysis. Comput Biol Med 2020; 128:104159. [PMID: 33301952 DOI: 10.1016/j.compbiomed.2020.104159] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Revised: 11/14/2020] [Accepted: 11/29/2020] [Indexed: 10/22/2022]
Abstract
Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, about one-third of patients who are implanted, derive no measurable benefit from CRT. Non-response may partly be due to suboptimal activation of the left ventricle (LV) caused by electrophysiological heterogeneities. The goal of this study is to investigate the performance of a newly developed method used to analyze electrical wavefront propagation in a heart model including myocardial scar and compare this to clinical benchmark studies. We used computational models to measure the maximum activation front (MAF) in the LV during different pacing scenarios. Different heart geometries and scars were created based on cardiac MR images of three patients. The right ventricle (RV) was paced from the apex and the LV was paced from 12 different sites, single site, dual-site and triple site. Our results showed that for single LV site pacing, the pacing site with the largest MAF corresponded with the latest activated regions of the LV demonstrated during RV pacing, which also agrees with previous markers used for predicting optimal single-site pacing location. We then demonstrated the utility of MAF in predicting optimal electrode placements in more complex scenarios including scar and multi-site LV pacing. This study demonstrates the potential value of computational simulations in understanding and planning CRT.
Collapse
Affiliation(s)
- Mohammad Albatat
- Intervention Centre, Oslo University Hospital, Oslo, Norway; Institute of Clinical Medicine, University of Oslo, Oslo, Norway.
| | - Hermenegild Arevalo
- Department of Computational Physiology, Simula Research Laboratory, Fornebu, Norway
| | | | - Vilde Strøm
- Department of Computational Physiology, Simula Research Laboratory, Fornebu, Norway
| | - Ilangko Balasingham
- Intervention Centre, Oslo University Hospital, Oslo, Norway; Department of Electronic Systems, Norwegian University of Science and Technology, Trondheim, Norway
| | | |
Collapse
|
27
|
Del Corso G, Verzicco R, Viola F. Sensitivity analysis of an electrophysiology model for the left ventricle. J R Soc Interface 2020; 17:20200532. [PMID: 33109017 DOI: 10.1098/rsif.2020.0532] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Modelling the cardiac electrophysiology entails dealing with the uncertainties related to the input parameters such as the heart geometry and the electrical conductivities of the tissues, thus calling for an uncertainty quantification (UQ) of the results. Since the chambers of the heart have different shapes and tissues, in order to make the problem affordable, here we focus on the left ventricle with the aim of identifying which of the uncertain inputs mostly affect its electrophysiology. In a first phase, the uncertainty of the input parameters is evaluated using data available from the literature and the output quantities of interest (QoIs) of the problem are defined. According to the polynomial chaos expansion, a training dataset is then created by sampling the parameter space using a quasi-Monte Carlo method whereas a smaller independent dataset is used for the validation of the resulting metamodel. The latter is exploited to run a global sensitivity analysis with nonlinear variance-based indices and thus reduce the input parameter space accordingly. Thereafter, the uncertainty probability distribution of the QoIs are evaluated using a direct UQ strategy on a larger dataset and the results discussed in the light of the medical knowledge.
Collapse
Affiliation(s)
| | - Roberto Verzicco
- Gran Sasso Science Institute (GSSI), L'Aquila, Italy.,University of Rome Tor Vergata, Rome, Italy.,POF Group, University of Twente, Enschede, The Netherlands
| | | |
Collapse
|
28
|
Integration of activation maps of epicardial veins in computational cardiac electrophysiology. Comput Biol Med 2020; 127:104047. [PMID: 33099220 DOI: 10.1016/j.compbiomed.2020.104047] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2020] [Revised: 10/06/2020] [Accepted: 10/06/2020] [Indexed: 12/16/2022]
Abstract
In this work we address the issue of validating the monodomain equation used in combination with the Bueno-Orovio ionic model for the prediction of the activation times in cardiac electro-physiology of the left ventricle. To this aim, we consider four patients who suffered from Left Bundle Branch Block (LBBB). We use activation maps performed at the septum as input data for the model and maps at the epicardial veins for the validation. In particular, a first set (half) of the latter are used to estimate the conductivities of the patient and a second set (the remaining half) to compute the errors of the numerical simulations. We find an excellent agreement between measures and numerical results. Our validated computational tool could be used to accurately predict activation times at the epicardial veins with a short mapping, i.e. by using only a part (the most proximal) of the standard acquisition points, thus reducing the invasive procedure and exposure to radiation.
Collapse
|
29
|
Fan L, Namani R, Choy JS, Kassab GS, Lee LC. Effects of Mechanical Dyssynchrony on Coronary Flow: Insights From a Computational Model of Coupled Coronary Perfusion With Systemic Circulation. Front Physiol 2020; 11:915. [PMID: 32922304 PMCID: PMC7457036 DOI: 10.3389/fphys.2020.00915] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Accepted: 07/08/2020] [Indexed: 01/01/2023] Open
Abstract
Mechanical dyssynchrony affects left ventricular (LV) mechanics and coronary perfusion. Due to the confounding effects of their bi-directional interactions, the mechanisms behind these changes are difficult to isolate from experimental and clinical studies alone. Here, we develop and calibrate a closed-loop computational model that couples the systemic circulation, LV mechanics, and coronary perfusion. The model is applied to simulate the impact of mechanical dyssynchrony on coronary flow in the left anterior descending artery (LAD) and left circumflex artery (LCX) territories caused by regional alterations in perfusion pressure and intramyocardial pressure (IMP). We also investigate the effects of regional coronary flow alterations on regional LV contractility in mechanical dyssynchrony based on prescribed contractility-flow relationships without considering autoregulation. The model predicts that LCX and LAD flows are reduced by 7.2%, and increased by 17.1%, respectively, in mechanical dyssynchrony with a systolic dyssynchrony index of 10% when the LAD's IMP is synchronous with the arterial pressure. The LAD flow is reduced by 11.6% only when its IMP is delayed with respect to the arterial pressure by 0.07 s. When contractility is sensitive to coronary flow, mechanical dyssynchrony can affect global LV mechanics, IMPs and contractility that in turn, further affect the coronary flow in a feedback loop that results in a substantial reduction of dPLV/dt, indicative of ischemia. Taken together, these findings imply that regional IMPs play a significant role in affecting regional coronary flows in mechanical dyssynchrony and the changes in regional coronary flow may produce ischemia when contractility is sensitive to the changes in coronary flow.
Collapse
Affiliation(s)
- Lei Fan
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Ravi Namani
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Jenny S Choy
- California Medical Innovation Institute, San Diego, CA, United States
| | - Ghassan S Kassab
- California Medical Innovation Institute, San Diego, CA, United States
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| |
Collapse
|
30
|
Isotani A, Yoneda K, Iwamura T, Watanabe M, Okada JI, Washio T, Sugiura S, Hisada T, Ando K. Patient-specific heart simulation can identify non-responders to cardiac resynchronization therapy. Heart Vessels 2020; 35:1135-1147. [PMID: 32166443 PMCID: PMC7332486 DOI: 10.1007/s00380-020-01577-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Accepted: 02/28/2020] [Indexed: 11/30/2022]
Abstract
To identify non-responders to cardiac resynchronization therapy (CRT), various biomarkers have been proposed, but these attempts have not been successful to date. We tested the clinical applicability of computer simulation of CRT for the identification of non-responders. We used the multi-scale heart simulator “UT-Heart,” which can reproduce the electrophysiology and mechanics of the heart based on a molecular model of the excitation–contraction mechanism. Patient-specific heart models were created for eight heart failure patients who were treated with CRT, based on the clinical data recorded before treatment. Using these heart models, bi-ventricular pacing simulations were performed at multiple pacing sites adopted in clinical practice. Improvement in pumping function measured by the relative change of maximum positive derivative of left ventricular pressure (%ΔdP/dtmax) was compared with the clinical outcome. The operators of the simulation were blinded to the clinical outcome. In six patients, the relative reduction in end-systolic volume exceeded 15% in the follow-up echocardiogram at 3 months (responders) and the remaining two patients were judged as non-responders. The simulated %ΔdP/dtmax at the best lead position could identify responders and non-responders successfully. With further refinement of the model, patient-specific simulation could be a useful tool for identifying non-responders to CRT.
Collapse
Affiliation(s)
- Akihiro Isotani
- Department of Cardiovascular Medicine, Kokura Memorial Hospital, Asano 3-2-1, Kokurakita-ku, Kitakyushu, Fukuoka, 802-8555, Japan
| | - Kazunori Yoneda
- Healthcare System Unit, Fujitsu Ltd, Ota-ku, Kamata, 144-8588, Japan
| | - Takashi Iwamura
- Healthcare System Unit, Fujitsu Ltd, Ota-ku, Kamata, 144-8588, Japan
| | - Masahiro Watanabe
- Healthcare System Unit, Fujitsu Ltd, Ota-ku, Kamata, 144-8588, Japan
| | - Jun-Ichi Okada
- Future Center Initiative, The University of Tokyo, Wakashiba 178-4-4, Kashiwa, Chiba, 277-0871, Japan
- UT-Heart Inc. Nozawa, 3-25-8, Setagaya, Tokyo, 154-0003, Japan
| | - Takumi Washio
- Future Center Initiative, The University of Tokyo, Wakashiba 178-4-4, Kashiwa, Chiba, 277-0871, Japan
- UT-Heart Inc. Nozawa, 3-25-8, Setagaya, Tokyo, 154-0003, Japan
| | - Seiryo Sugiura
- UT-Heart Inc. Nozawa, 3-25-8, Setagaya, Tokyo, 154-0003, Japan.
- Future Center #304, Wakashiba 178-4-4, Kashiwa, Chiba, 277-0871, Japan.
| | - Toshiaki Hisada
- UT-Heart Inc. Nozawa, 3-25-8, Setagaya, Tokyo, 154-0003, Japan
| | - Kenji Ando
- Department of Cardiovascular Medicine, Kokura Memorial Hospital, Asano 3-2-1, Kokurakita-ku, Kitakyushu, Fukuoka, 802-8555, Japan
| |
Collapse
|
31
|
Kimmig F, Caruel M. Hierarchical modeling of force generation in cardiac muscle. Biomech Model Mechanobiol 2020; 19:2567-2601. [PMID: 32681201 DOI: 10.1007/s10237-020-01357-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2019] [Accepted: 06/10/2020] [Indexed: 11/25/2022]
Abstract
Performing physiologically relevant simulations of the beating heart in clinical context requires to develop detailed models of the microscale force generation process. These models, however, may reveal difficult to implement in practice due to their high computational costs and complex calibration. We propose a hierarchy of three interconnected muscle contraction models-from the more refined to the more simplified-that are rigorously and systematically related to each other, offering a way to select, for a specific application, the model that yields a good trade-off between physiological fidelity, computational cost and calibration complexity. The three model families are compared to the same set of experimental data to systematically assess what physiological indicators can be reproduced or not and how these indicators constrain the model parameters. Finally, we discuss the applicability of these models for heart simulation.
Collapse
Affiliation(s)
- François Kimmig
- LMS, CNRS, École polytechnique, Institut Polytechnique de Paris, Paris, France.
- Inria, Inria Saclay-Ile-de-France, Palaiseau, France.
| | | |
Collapse
|
32
|
Lei CL, Ghosh S, Whittaker DG, Aboelkassem Y, Beattie KA, Cantwell CD, Delhaas T, Houston C, Novaes GM, Panfilov AV, Pathmanathan P, Riabiz M, dos Santos RW, Walmsley J, Worden K, Mirams GR, Wilkinson RD. Considering discrepancy when calibrating a mechanistic electrophysiology model. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190349. [PMID: 32448065 PMCID: PMC7287333 DOI: 10.1098/rsta.2019.0349] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 04/21/2020] [Indexed: 05/21/2023]
Abstract
Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions-that of uncertainty in the model structure or the equations themselves. The difference between imperfect models and reality is termed model discrepancy, and we are often uncertain as to the size and consequences of this discrepancy. Here, we provide two examples of the consequences of discrepancy when calibrating models at the ion channel and action potential scales. Furthermore, we attempt to account for this discrepancy when calibrating and validating an ion channel model using different methods, based on modelling the discrepancy using Gaussian processes and autoregressive-moving-average models, then highlight the advantages and shortcomings of each approach. Finally, suggestions and lines of enquiry for future work are provided. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- Chon Lok Lei
- Computational Biology and Health Informatics, Department of Computer Science, University of Oxford, Oxford, UK
| | - Sanmitra Ghosh
- MRC Biostatistics Unit, University of Cambridge, Cambridge, UK
| | - Dominic G. Whittaker
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Yasser Aboelkassem
- Department of Bioengineering, University of California San Diego, La Jolla, CA, USA
| | - Kylie A. Beattie
- Systems Modeling and Translational Biology, GlaxoSmithKline R&D, Stevenage, UK
| | - Chris D. Cantwell
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Tammo Delhaas
- CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Charles Houston
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Gustavo Montes Novaes
- Graduate Program in Computational Modeling, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
| | - Alexander V. Panfilov
- Department of Physics and Astronomy, Ghent University, Ghent, Belgium
- Laboratory of Computational Biology and Medicine, Ural Federal University, Ekaterinburg, Russia
| | - Pras Pathmanathan
- US Food and Drug Administration, Center for Devices and Radiological Health, Office of Science and Engineering Laboratories, Silver Spring, MD, USA
| | - Marina Riabiz
- Department of Biomedical Engineering King’s College London and Alan Turing Institute, London, UK
| | - Rodrigo Weber dos Santos
- Graduate Program in Computational Modeling, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
| | - John Walmsley
- James T. Willerson Center for Cardiovascular Modeling and Simulation, Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, USA
| | - Keith Worden
- Dynamics Research Group, Department of Mechanical Engineering, University of Sheffield, Sheffield, UK
| | - Gary R. Mirams
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | | |
Collapse
|
33
|
Clayton RH, Aboelkassem Y, Cantwell CD, Corrado C, Delhaas T, Huberts W, Lei CL, Ni H, Panfilov AV, Roney C, dos Santos RW. An audit of uncertainty in multi-scale cardiac electrophysiology models. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190335. [PMID: 32448070 PMCID: PMC7287340 DOI: 10.1098/rsta.2019.0335] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/16/2020] [Indexed: 05/21/2023]
Abstract
Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- Richard H. Clayton
- Insigneo institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, UK
- e-mail:
| | - Yasser Aboelkassem
- Department of Bioengineering, University of California, San Diego, CA, USA
| | | | - Cesare Corrado
- Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | - Tammo Delhaas
- School of Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Wouter Huberts
- School of Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Chon Lok Lei
- Computational Biology and Health Informatics, Department of Computer Science, University of Oxford, Oxford, UK
| | - Haibo Ni
- Department of Pharmacology, University of California, Davis, CA, USA
| | - Alexander V. Panfilov
- Department of Physics and Astronomy, University of Gent, Gent, Belgium
- Laboratory of Computational Biology and Medicine, Ural Federal University, Ekaterinburg, Russia
| | - Caroline Roney
- Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | | |
Collapse
|
34
|
Le Gall A, Vallée F, Pushparajah K, Hussain T, Mebazaa A, Chapelle D, Gayat É, Chabiniok R. Monitoring of cardiovascular physiology augmented by a patient-specific biomechanical model during general anesthesia. A proof of concept study. PLoS One 2020; 15:e0232830. [PMID: 32407353 PMCID: PMC7224549 DOI: 10.1371/journal.pone.0232830] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 04/22/2020] [Indexed: 12/29/2022] Open
Abstract
During general anesthesia (GA), direct analysis of arterial pressure or aortic flow waveforms may be inconclusive in complex situations. Patient-specific biomechanical models, based on data obtained during GA and capable to perform fast simulations of cardiac cycles, have the potential to augment hemodynamic monitoring. Such models allow to simulate Pressure-Volume (PV) loops and estimate functional indicators of cardiovascular (CV) system, e.g. ventricular-arterial coupling (Vva), cardiac efficiency (CE) or myocardial contractility, evolving throughout GA. In this prospective observational study, we created patient-specific biomechanical models of heart and vasculature of a reduced geometric complexity for n = 45 patients undergoing GA, while using transthoracic echocardiography and aortic pressure and flow signals acquired in the beginning of GA (baseline condition). If intraoperative hypotension (IOH) appeared, diluted norepinephrine (NOR) was administered and the model readjusted according to the measured aortic pressure and flow signals. Such patients were a posteriori assigned into a so-called hypotensive group. The accuracy of simulated mean aortic pressure (MAP) and stroke volume (SV) at baseline were in accordance with the guidelines for the validation of new devices or reference measurement methods in all patients. After NOR administration in the hypotensive group, the percentage of concordance with 10% exclusion zone between measurement and simulation was >95% for both MAP and SV. The modeling results showed a decreased Vva (0.64±0.37 vs 0.88±0.43; p = 0.039) and an increased CE (0.8±0.1 vs 0.73±0.11; p = 0.042) in hypotensive vs normotensive patients. Furthermore, Vva increased by 92±101%, CE decreased by 13±11% (p < 0.001 for both) and contractility increased by 14±11% (p = 0.002) in the hypotensive group post-NOR administration. In this work we demonstrated the application of fast-running patient-specific biophysical models to estimate PV loops and functional indicators of CV system using clinical data available during GA. The work paves the way for model-augmented hemodynamic monitoring at operating theatres or intensive care units to enhance the information on patient-specific physiology.
Collapse
Affiliation(s)
- Arthur Le Gall
- Inria, Paris, France
- LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Paris, France
- Anesthesiology and Intensive Care Department, Lariboisière - Saint Louis - Fernand Widal University Hospitals, Paris, France
- INSERM, Paris, France
| | - Fabrice Vallée
- Inria, Paris, France
- LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Paris, France
- Anesthesiology and Intensive Care Department, Lariboisière - Saint Louis - Fernand Widal University Hospitals, Paris, France
- INSERM, Paris, France
| | - Kuberan Pushparajah
- School of Biomedical Engineering & Imaging Sciences, St Thomas’ Hospital, King’s College London, London, United Kingdom
| | - Tarique Hussain
- Department of Pediatrics, Division of Pediatric Cardiology, UT Southwestern Medical Center, Dallas, TX, United States of America
| | - Alexandre Mebazaa
- Anesthesiology and Intensive Care Department, Lariboisière - Saint Louis - Fernand Widal University Hospitals, Paris, France
- INSERM, Paris, France
| | - Dominique Chapelle
- Inria, Paris, France
- LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Paris, France
| | - Étienne Gayat
- Anesthesiology and Intensive Care Department, Lariboisière - Saint Louis - Fernand Widal University Hospitals, Paris, France
- INSERM, Paris, France
| | - Radomír Chabiniok
- Inria, Paris, France
- LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Paris, France
- School of Biomedical Engineering & Imaging Sciences, St Thomas’ Hospital, King’s College London, London, United Kingdom
- Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Prague, Czech Republic
| |
Collapse
|
35
|
Giannakidis A, Gullberg GT. Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI. Front Physiol 2020; 11:265. [PMID: 32296341 PMCID: PMC7136532 DOI: 10.3389/fphys.2020.00265] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Accepted: 03/09/2020] [Indexed: 11/16/2022] Open
Abstract
The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λL), transversal diffusivity (λT), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.
Collapse
Affiliation(s)
- Archontis Giannakidis
- School of Science and Technology, Nottingham Trent University, Nottingham, United Kingdom.,Life Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States.,National Heart & Lung Institute, Imperial College London, London, United Kingdom
| | - Grant T Gullberg
- Life Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States.,Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
| |
Collapse
|
36
|
Dhamala J, Bajracharya P, Arevalo HJ, Sapp JL, Horácek BM, Wu KC, Trayanova NA, Wang L. Embedding high-dimensional Bayesian optimization via generative modeling: Parameter personalization of cardiac electrophysiological models. Med Image Anal 2020; 62:101670. [PMID: 32171168 DOI: 10.1016/j.media.2020.101670] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2019] [Revised: 12/16/2019] [Accepted: 02/24/2020] [Indexed: 11/28/2022]
Abstract
The estimation of patient-specific tissue properties in the form of model parameters is important for personalized physiological models. Because tissue properties are spatially varying across the underlying geometrical model, it presents a significant challenge of high-dimensional (HD) optimization at the presence of limited measurement data. A common solution to reduce the dimension of the parameter space is to explicitly partition the geometrical mesh. In this paper, we present a novel concept that uses a generative variational auto-encoder (VAE) to embed HD Bayesian optimization into a low-dimensional (LD) latent space that represents the generative code of HD parameters. We further utilize VAE-encoded knowledge about the generative code to guide the exploration of the search space. The presented method is applied to estimating tissue excitability in a cardiac electrophysiological model in a range of synthetic and real-data experiments, through which we demonstrate its improved accuracy and substantially reduced computational cost in comparison to existing methods that rely on geometry-based reduction of the HD parameter space.
Collapse
Affiliation(s)
- Jwala Dhamala
- Rochester Institute of Technology, Rochester, NY, USA. http://www.jwaladhamala.com
| | | | | | | | | | | | | | - Linwei Wang
- Rochester Institute of Technology, Rochester, NY, USA.
| |
Collapse
|
37
|
Dobutamine stress testing in patients with Fontan circulation augmented by biomechanical modeling. PLoS One 2020; 15:e0229015. [PMID: 32084180 PMCID: PMC7034893 DOI: 10.1371/journal.pone.0229015] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Accepted: 01/28/2020] [Indexed: 02/02/2023] Open
Abstract
Understanding (patho)physiological phenomena and mechanisms of failure in patients with Fontan circulation-a surgically established circulation for patients born with a functionally single ventricle-remains challenging due to the complex hemodynamics and high inter-patient variations in anatomy and function. In this work, we present a biomechanical model of the heart and circulation to augment the diagnostic evaluation of Fontan patients with early-stage heart failure. The proposed framework employs a reduced-order model of heart coupled with a simplified circulation including venous return, creating a closed-loop system. We deploy this framework to augment the information from data obtained during combined cardiac catheterization and magnetic resonance exams (XMR), performed at rest and during dobutamine stress in 9 children with Fontan circulation and 2 biventricular controls. We demonstrate that our modeling framework enables patient-specific investigation of myocardial stiffness, contractility at rest, contractile reserve during stress and changes in vascular resistance. Hereby, the model allows to identify key factors underlying the pathophysiological response to stress in these patients. In addition, the rapid personalization of the model to patient data and fast simulation of cardiac cycles make our framework directly applicable in a clinical workflow. We conclude that the proposed modeling framework is a valuable addition to the current clinical diagnostic XMR exam that helps to explain patient-specific stress hemodynamics and can identify potential mechanisms of failure in patients with Fontan circulation.
Collapse
|
38
|
Phung TKN, Waters CD, Holmes JW. Open-Source Routines for Building Personalized Left Ventricular Models From Cardiac Magnetic Resonance Imaging Data. J Biomech Eng 2020; 142:024504. [PMID: 31141592 PMCID: PMC7104752 DOI: 10.1115/1.4043876] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 05/20/2019] [Indexed: 11/08/2022]
Abstract
Creating patient-specific models of the heart is a promising approach for predicting outcomes in response to congenital malformations, injury, or disease, as well as an important tool for developing and customizing therapies. However, integrating multimodal imaging data to construct patient-specific models is a nontrivial task. Here, we propose an approach that employs a prolate spheroidal coordinate system to interpolate information from multiple imaging datasets and map those data onto a single geometric model of the left ventricle (LV). We demonstrate the mapping of the location and transmural extent of postinfarction scar segmented from late gadolinium enhancement (LGE) magnetic resonance imaging (MRI), as well as mechanical activation calculated from displacement encoding with stimulated echoes (DENSE) MRI. As a supplement to this paper, we provide MATLAB and Python versions of the routines employed here for download from SimTK.
Collapse
Affiliation(s)
- Thien-Khoi N. Phung
- Department of Biomedical Engineering, University of
Virginia, Charlottesville, VA 22908
| | - Christopher D. Waters
- Department of Biomedical Engineering, University of
Virginia, Charlottesville, VA 22908
| | - Jeffrey W. Holmes
- Department of Biomedical Engineering, University of
Virginia, Charlottesville, VA 22908;
Department of Medicine, University of Virginia,
Charlottesville, VA 22908; Robert M. Berne Cardiovascular
Center, University of Virginia, 415 Lane Road,
Charlottesville, VA 22908
| |
Collapse
|
39
|
Strocchi M, Gsell MAF, Augustin CM, Razeghi O, Roney CH, Prassl AJ, Vigmond EJ, Behar JM, Gould JS, Rinaldi CA, Bishop MJ, Plank G, Niederer SA. Simulating ventricular systolic motion in a four-chamber heart model with spatially varying robin boundary conditions to model the effect of the pericardium. J Biomech 2020; 101:109645. [PMID: 32014305 PMCID: PMC7677892 DOI: 10.1016/j.jbiomech.2020.109645] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2019] [Revised: 01/15/2020] [Accepted: 01/15/2020] [Indexed: 12/11/2022]
Abstract
The pericardium affects cardiac motion by limiting epicardial displacement normal to the surface. In computational studies, it is important for the model to replicate realistic motion, as this affects the physiological fidelity of the model. Previous computational studies showed that accounting for the effect of the pericardium allows for a more realistic motion simulation. In this study, we describe the mechanism through which the pericardium causes improved cardiac motion. We simulated electrical activation and contraction of the ventricles on a four-chamber heart in the presence and absence of the effect of the pericardium. We simulated the mechanical constraints imposed by the pericardium by applying normal Robin boundary conditions on the ventricular epicardium. We defined a regional scaling of normal springs stiffness based on image-derived motion from CT images. The presence of the pericardium reduced the error between simulated and image-derived end-systolic configurations from 12.8±4.1 mm to 5.7±2.5 mm. First, the pericardium prevents the ventricles from spherising during isovolumic contraction, reducing the outward motion of the free walls normal to the surface and the upwards motion of the apex. Second, by restricting the inward motion of the free and apical walls of the ventricles the pericardium increases atrioventricular plane displacement by four folds during ejection. Our results provide a mechanistic explanation of the importance of the pericardium in physiological simulations of electromechanical cardiac function.
Collapse
Affiliation(s)
- Marina Strocchi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
| | | | | | - Orod Razeghi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Caroline H Roney
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Anton J Prassl
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Edward J Vigmond
- University of Bordeaux, Talence, France; LIRYC Electrophysiology and Heart Modeling Institute, Campus Xavier Arnozan, Pessac, France
| | - Jonathan M Behar
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Justin S Gould
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Gernot Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Steven A Niederer
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| |
Collapse
|
40
|
Abstract
The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.
Collapse
|
41
|
Genet M. A relaxed growth modeling framework for controlling growth-induced residual stresses. Clin Biomech (Bristol, Avon) 2019; 70:270-277. [PMID: 31831206 DOI: 10.1016/j.clinbiomech.2019.08.015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Revised: 08/26/2019] [Accepted: 08/27/2019] [Indexed: 02/07/2023]
Abstract
BACKGROUND Constitutive models of the mechanical response of soft tissues have been established and are widely accepted, but models of soft tissues remodeling are more controversial. Specifically for growth, one important question arises pertaining to residual stresses: existing growth models inevitably introduce residual stresses, but it is not entirely clear if this is physiological or merely an artifact of the modeling framework. As a consequence, in simulating growth, some authors have chosen to keep growth-induced residual stresses, and others have chosen to remove them. METHODS In this paper, we introduce a novel "relaxed growth" framework allowing for a fine control of the amount of residual stresses generated during tissue growth. It is a direct extension of the classical framework of the multiplicative decomposition of the transformation gradient, to which an additional sub-transformation is introduced in order to let the original unloaded configuration evolve, hence relieving some residual stresses. We provide multiple illustrations of the framework mechanical response, on time-driven constrained growth as well as the strain-driven growth problem of the artery under internal pressure, including the opening angle experiment. FINDINGS The novel relaxed growth modeling framework introduced in this paper allows for a better control of growth-induced residual stresses compared to standard growth models based on the multiplicative decomposition of the transformation gradient. INTERPRETATION Growth-induced residual stresses should be better handled in soft tissues biomechanical models, especially in patient-specific models of diseased organs that are aimed at augmented diagnosis and treatment optimization.
Collapse
Affiliation(s)
- M Genet
- Laboratoire de Mécanique des Solides, École Polytechnique/Institut Polytechnique de Paris/CNRS, Palaiseau, France; M3DISIM Team, INRIA/Université Paris-Saclay, Palaiseau, France.
| |
Collapse
|
42
|
Technological and Clinical Challenges in Lead Placement for Cardiac Rhythm Management Devices. Ann Biomed Eng 2019; 48:26-46. [DOI: 10.1007/s10439-019-02376-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2019] [Accepted: 09/25/2019] [Indexed: 01/29/2023]
|
43
|
Lee AWC, Nguyen UC, Razeghi O, Gould J, Sidhu BS, Sieniewicz B, Behar J, Mafi-Rad M, Plank G, Prinzen FW, Rinaldi CA, Vernooy K, Niederer S. A rule-based method for predicting the electrical activation of the heart with cardiac resynchronization therapy from non-invasive clinical data. Med Image Anal 2019; 57:197-213. [PMID: 31326854 PMCID: PMC6746621 DOI: 10.1016/j.media.2019.06.017] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2018] [Revised: 06/20/2019] [Accepted: 06/27/2019] [Indexed: 12/13/2022]
Abstract
Background Cardiac Resynchronization Therapy (CRT) is one of the few effective treatments for heart failure patients with ventricular dyssynchrony. The pacing location of the left ventricle is indicated as a determinant of CRT outcome. Objective Patient specific computational models allow the activation pattern following CRT implant to be predicted and this may be used to optimize CRT lead placement. Methods In this study, the effects of heterogeneous cardiac substrate (scar, fast endocardial conduction, slow septal conduction, functional block) on accurately predicting the electrical activation of the LV epicardium were tested to determine the minimal detail required to create a rule based model of cardiac electrophysiology. Non-invasive clinical data (CT or CMR images and 12 lead ECG) from eighteen patients from two centers were used to investigate the models. Results Validation with invasive electro-anatomical mapping data identified that computer models with fast endocardial conduction were able to predict the electrical activation with a mean distance errors of 9.2 ± 0.5 mm (CMR data) or (CT data) 7.5 ± 0.7 mm. Conclusion This study identified a simple rule-based fast endocardial conduction model, built using non-invasive clinical data that can be used to rapidly and robustly predict the electrical activation of the heart. Pre-procedural prediction of the latest electrically activating region to identify the optimal LV pacing site could potentially be a useful clinical planning tool for CRT procedures.
Collapse
Affiliation(s)
- A W C Lee
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.
| | - U C Nguyen
- Department of Physiology, Maastricht University Medical Center (MUMC+), Cardiovascular Research Institute Maastricht (CARIM), Maastricht, the Netherlands; Department of Cardiology, Maastricht University Medical Center (MUMC+), Cardiovascular Research Institute Maastricht (CARIM), Maastricht, the Netherlands
| | - O Razeghi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - J Gould
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - B S Sidhu
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - B Sieniewicz
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - J Behar
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Bart's Heart Centre, St. Bartholomew's Hospital, London, United Kingdom
| | - M Mafi-Rad
- Department of Cardiology, Maastricht University Medical Center (MUMC+), Cardiovascular Research Institute Maastricht (CARIM), Maastricht, the Netherlands
| | - G Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - F W Prinzen
- Department of Physiology, Maastricht University Medical Center (MUMC+), Cardiovascular Research Institute Maastricht (CARIM), Maastricht, the Netherlands
| | - C A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - K Vernooy
- Department of Cardiology, Maastricht University Medical Center (MUMC+), Cardiovascular Research Institute Maastricht (CARIM), Maastricht, the Netherlands; Department of Cardiology, Radboud University Medical Center, Nijmegen, the Netherlands
| | - S Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| |
Collapse
|
44
|
Katbeh A, Van Camp G, Barbato E, Galderisi M, Trimarco B, Bartunek J, Vanderheyden M, Penicka M. Cardiac Resynchronization Therapy Optimization: A Comprehensive Approach. Cardiology 2019; 142:116-128. [PMID: 31117077 DOI: 10.1159/000499192] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 02/26/2019] [Indexed: 11/19/2022]
Abstract
Since the first report on biventricular pacing in 1994, cardiac resynchronization therapy (CRT) has become standard for patients with advanced heart failure (HF) and ventricular conduction delay. CRT improves myocardial function by resynchronizing myocardial contraction, which results in reverse left ventricular remodeling and improves symptoms and clinical outcomes. Despite the accelerated development of CRT device technology and its increased application in treating HF patients, almost one-third of these patients do not respond to the therapy or gain any clinical benefit from device implantation. Over the last decade, multiple cardiac imaging modalities have provided a deeper understanding of myocardial pathophysiology, thereby improving HF treatment management. However, the optimal strategy for improving the CRT response remains debatable. This article provides an updated overview of the electropathophysiology of myocardial dysfunction in ventricular conduction delay and the diagnostic approaches involving the use of multiple modalities.
Collapse
Affiliation(s)
- Asim Katbeh
- Cardiovascular Center Aalst, OLV Clinic, Aalst, Belgium.,Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Guy Van Camp
- Cardiovascular Center Aalst, OLV Clinic, Aalst, Belgium
| | - Emanuele Barbato
- Cardiovascular Center Aalst, OLV Clinic, Aalst, Belgium.,Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Maurizio Galderisi
- Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Bruno Trimarco
- Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | | | | | - Martin Penicka
- Cardiovascular Center Aalst, OLV Clinic, Aalst, Belgium,
| |
Collapse
|
45
|
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
|
46
|
Niederer SA, Campbell KS, Campbell SG. A short history of the development of mathematical models of cardiac mechanics. J Mol Cell Cardiol 2019; 127:11-19. [PMID: 30503754 PMCID: PMC6525149 DOI: 10.1016/j.yjmcc.2018.11.015] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/30/2018] [Revised: 11/02/2018] [Accepted: 11/21/2018] [Indexed: 11/15/2022]
Abstract
Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.
Collapse
Affiliation(s)
| | - Kenneth S Campbell
- Department of Physiology and Division of Cardiovascular Medicine, University of Kentucky, Lexington, USA
| | - Stuart G Campbell
- Departments of Biomedical Engineering and Cellular and Molecular Physiology, Yale University, New Haven, USA
| |
Collapse
|
47
|
The importance of the pericardium for cardiac biomechanics: from physiology to computational modeling. Biomech Model Mechanobiol 2018; 18:503-529. [DOI: 10.1007/s10237-018-1098-4] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Accepted: 11/18/2018] [Indexed: 10/27/2022]
|
48
|
Kirn B, Walmsley J, Lumens J. Uniqueness of local myocardial strain patterns with respect to activation time and contractility of the failing heart: a computational study. Biomed Eng Online 2018; 17:182. [PMID: 30518387 PMCID: PMC6280493 DOI: 10.1186/s12938-018-0614-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2018] [Accepted: 11/27/2018] [Indexed: 01/26/2023] Open
Abstract
Background Myocardial deformation measured by strain is used to detect electro-mechanical abnormalities in cardiac tissue. Estimation of myocardial properties from regional strain patterns when multiple pathologies are present is therefore a promising application of computer modelling. However, if different tissue properties lead to indistinguishable strain patterns (‘degeneracy’), the applicability of any such method will be limited. We investigated whether estimation of local activation time (AT) and contractility from myocardial strain patterns is theoretically possible. Methods For four different global cardiac pathologies local myocardial strain patterns for 1025 combinations of AT and contractility were simulated with a computational model (CircAdapt). For each strain pattern, a cohort of similar patterns was found within estimated measurement error using the sum of least-squared differences. Cohort members came from (1) the same pathology only, and (2) all four pathologies. Uncertainty was calculated as accuracy and precision of cohort members in parameter space. Connectedness within the cohorts was also studied. Results We found that cohorts drawn from one pathology had parameters with adjacent values although their distribution was neither constant nor symmetrical. In comparison cohorts drawn from four pathologies had disconnected components with drastically different parameter values and accuracy and precision values up to three times higher. Conclusions Global pathology must be known when extracting AT and contractility from strain patterns, otherwise degeneracy occurs causing unacceptable uncertainty in derived parameters. Electronic supplementary material The online version of this article (10.1186/s12938-018-0614-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Borut Kirn
- Department of Physiology, Medical Faculty, University of Ljubljana, Zaloska 4, 1000, Ljubljana, Slovenia.
| | - John Walmsley
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University Medical Center, Maastricht, The Netherlands
| |
Collapse
|
49
|
Santiago A, Aguado-Sierra J, Zavala-Aké M, Doste-Beltran R, Gómez S, Arís R, Cajas JC, Casoni E, Vázquez M. Fully coupled fluid-electro-mechanical model of the human heart for supercomputers. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e3140. [PMID: 30117302 DOI: 10.1002/cnm.3140] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Revised: 04/28/2018] [Accepted: 07/22/2018] [Indexed: 05/12/2023]
Abstract
In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.
Collapse
Affiliation(s)
- Alfonso Santiago
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Jazmín Aguado-Sierra
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Miguel Zavala-Aké
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | | | - Samuel Gómez
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Ruth Arís
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Juan C Cajas
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Eva Casoni
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Mariano Vázquez
- Department of Computer Applications in Science and Engineering, Barcelona Supercomputing Center (BSC), Barcelona, Spain
- Instituto de Investigación en Inteligencia Artificial (IIIA), Consejo Superior de Investigaciones Científicas (CSIC), Barcelona, Spain
| |
Collapse
|
50
|
Albatat M, King DR, Unger LA, Arevalo H, Wall S, Sundnes J, Bergsland J, Balasingham I. Electromechanical Model to Predict Cardiac Resynchronization Therapy. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2018; 2018:5446-5459. [PMID: 30441569 DOI: 10.1109/embc.2018.8513539] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, one-third of the CRT patients derive no measurable benefit from CRT, due to suboptimal placement of the left ventricular (LV) lead. We introduce a pipeline for improved CRT-therapy by creating an electromechanical model using patient-specific geometric parameters allowing individualization of therapy. The model successfully mimics expected changes when variables for tension, stiffness, and conduction are entered. Changing LV pacing site had a notable effect on maximum pressure gradient (dP/dtmax) in the presence of cardiac scarring, causing non-uniform excitation propagation through the LV. Tailoring CRT to the individual requires simulations with patient-specific biventricular meshes including cardiac geometry and conductivity properties.
Collapse
|