1
|
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
|
2
|
Leenknegt L, Panfilov AV, Dierckx H. Impact of electrode orientation, myocardial wall thickness, and myofiber direction on intracardiac electrograms: numerical modeling and analytical solutions. Front Physiol 2023; 14:1213218. [PMID: 37492643 PMCID: PMC10364610 DOI: 10.3389/fphys.2023.1213218] [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: 04/27/2023] [Accepted: 06/14/2023] [Indexed: 07/27/2023] Open
Abstract
Intracardiac electrograms (iEGMs) are time traces of the electrical potential recorded close to the heart muscle. We calculate unipolar and bipolar iEGMs analytically for a myocardial slab with parallel myofibers and validate them against numerical bidomain simulations. The analytical solution obtained via the method of mirrors is an infinite series of arctangents. It goes beyond the solid angle theory and is in good agreement with the simulations, even though bath loading effects were not accounted for in the analytical calculation. At a large distance from the myocardium, iEGMs decay as 1/R (unipolar), 1/R 2 (bipolar and parallel), and 1/R 3 (bipolar and perpendicular to the endocardium). At the endocardial surface, there is a mathematical branch cut. Here, we show how a thicker myocardium generates iEGMs with larger amplitudes and how anisotropy affects the iEGM width and amplitude. If only the leading-order term of our expansion is retained, it can be determined how the conductivities of the bath, torso, myocardium, and myofiber direction together determine the iEGM amplitude. Our results will be useful in the quantitative interpretation of iEGMs, the selection of thresholds to characterize viable tissues, and for future inferences of tissue parameters.
Collapse
Affiliation(s)
- Lore Leenknegt
- Department of Mathematics, KU Leuven Campus KULAK, KU Leuven, Kortrijk, Belgium
- iSi Health–KU Leuven Institute of Physics-based Modeling for In Silico Health, KU Leuven, Leuven, Belgium
| | | | - Hans Dierckx
- Department of Mathematics, KU Leuven Campus KULAK, KU Leuven, Kortrijk, Belgium
- iSi Health–KU Leuven Institute of Physics-based Modeling for In Silico Health, KU Leuven, Leuven, Belgium
| |
Collapse
|
3
|
Verzaal NJ, van Deursen CJM, Pezzuto S, Wecke L, van Everdingen WM, Vernooy K, Delhaas T, Auricchio A, Prinzen FW. Synchronization of repolarization after cardiac resynchronization therapy: a combined clinical and modeling study. J Cardiovasc Electrophysiol 2022; 33:1837-1846. [PMID: 35662306 PMCID: PMC9539692 DOI: 10.1111/jce.15581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Revised: 05/27/2022] [Accepted: 05/30/2022] [Indexed: 11/29/2022]
Abstract
Introduction The changes in ventricular repolarization after cardiac resynchronization therapy (CRT) are poorly understood. This knowledge gap is addressed using a multimodality approach including electrocardiographic and echocardiographic measurements in patients and using patient‐specific computational modeling. Methods In 33 patients electrocardiographic and echocardiographic measurements were performed before and at various intervals after CRT, both during CRT‐ON and temporary CRT‐OFF. T‐wave area was calculated from vectorcardiograms, and reconstructed from the 12‐lead electrocardiography (ECG). Computer simulations were performed using a patient‐specific eikonal model of cardiac activation with spatially varying action potential duration (APD) and repolarization rate, fit to a patient's ECG. Results During CRT‐ON T‐wave area diminished within a day and remained stable thereafter, whereas QT‐interval did not change significantly. During CRT‐OFF T‐wave area doubled within 5 days of CRT, while QT‐interval and peak‐to‐end T‐wave interval hardly changed. Left ventricular (LV) ejection fraction only increased significantly increased after 1 month of CRT. Computer simulations indicated that the increase in T‐wave area during CRT‐OFF can be explained by changes in APD following chronic CRT that are opposite to the change in CRT‐induced activation time. These APD changes were associated with a reduction in LV dispersion in repolarization during chronic CRT. Conclusion T‐wave area during CRT‐OFF is a sensitive marker for adaptations in ventricular repolarization during chronic CRT that may include a reduction in LV dispersion of repolarization.
Collapse
Affiliation(s)
- Nienke J Verzaal
- Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, P.O. Box 616, 6200 MD, Maastricht, the Netherlands
| | - Caroline J M van Deursen
- Department of Cardiology, Maastricht University Medical Center, P.O. Box 5800, 6202 AZ, Maastricht, The Netherlands
| | - Simone Pezzuto
- Center for Computational Medicine in Cardiology, Euler Institute, Università della Svizzera Italiana, via via Buffi 13, 6900, Lugano, Switzerland
| | - Liliane Wecke
- Heart Clinic, Capio St. Göran's Hospital, Sankt Göransplan 1, 11281, Stockholm, Sweden
| | | | - Kevin Vernooy
- Department of Cardiology, Maastricht University Medical Center, P.O. Box 5800, 6202 AZ, Maastricht, The Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, Maastricht University, P.O. Box 616, 6200 MD, Maastricht, the Netherlands
| | - Angelo Auricchio
- Center for Computational Medicine in Cardiology, Euler Institute, Università della Svizzera Italiana, via via Buffi 13, 6900, Lugano, Switzerland.,Department of Cardiology, Istituto Cardiocentro Ticino, Lugano, Switzerland
| | - Frits W Prinzen
- Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, P.O. Box 616, 6200 MD, Maastricht, the Netherlands
| |
Collapse
|
4
|
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
|
5
|
Plank G, Loewe A, Neic A, Augustin C, Huang YL, Gsell MAF, Karabelas E, Nothstein M, Prassl AJ, Sánchez J, Seemann G, Vigmond EJ. The openCARP simulation environment for cardiac electrophysiology. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2021; 208:106223. [PMID: 34171774 DOI: 10.1016/j.cmpb.2021.106223] [Citation(s) in RCA: 67] [Impact Index Per Article: 22.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Accepted: 05/28/2021] [Indexed: 06/13/2023]
Abstract
BACKGROUND AND OBJECTIVE Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community. METHODS AND RESULTS openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential. CONCLUSION As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.
Collapse
Affiliation(s)
- Gernot Plank
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria.
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | | | - Christoph Augustin
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Yung-Lin Huang
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg. Bad Krozingen, Medical Center - University of Freiburg, Freiburg, Germany; Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Matthias A F Gsell
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Elias Karabelas
- Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - Mark Nothstein
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Anton J Prassl
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Jorge Sánchez
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Gunnar Seemann
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg. Bad Krozingen, Medical Center - University of Freiburg, Freiburg, Germany; Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Edward J Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France; Université Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| |
Collapse
|
6
|
Pezzuto S, Prinzen FW, Potse M, Maffessanti F, Regoli F, Caputo ML, Conte G, Krause R, Auricchio A. Reconstruction of three-dimensional biventricular activation based on the 12-lead electrocardiogram via patient-specific modelling. Europace 2021; 23:640-647. [PMID: 33241411 PMCID: PMC8025079 DOI: 10.1093/europace/euaa330] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Accepted: 10/01/2020] [Indexed: 12/27/2022] Open
Abstract
Aims Non-invasive imaging of electrical activation requires high-density body surface potential mapping. The nine electrodes of the 12-lead electrocardiogram (ECG) are insufficient for a reliable reconstruction with standard inverse methods. Patient-specific modelling may offer an alternative route to physiologically constraint the reconstruction. The aim of the study was to assess the feasibility of reconstructing the fully 3D electrical activation map of the ventricles from the 12-lead ECG and cardiovascular magnetic resonance (CMR). Methods and results Ventricular activation was estimated by iteratively optimizing the parameters (conduction velocity and sites of earliest activation) of a patient-specific model to fit the simulated to the recorded ECG. Chest and cardiac anatomy of 11 patients (QRS duration 126–180 ms, documented scar in two) were segmented from CMR images. Scar presence was assessed by magnetic resonance (MR) contrast enhancement. Activation sequences were modelled with a physiologically based propagation model and ECGs with lead field theory. Validation was performed by comparing reconstructed activation maps with those acquired by invasive electroanatomical mapping of coronary sinus/veins (CS) and right ventricular (RV) and left ventricular (LV) endocardium. The QRS complex was correctly reproduced by the model (Pearson’s correlation r = 0.923). Reconstructions accurately located the earliest and latest activated LV regions (median barycentre distance 8.2 mm, IQR 8.8 mm). Correlation of simulated with recorded activation time was very good at LV endocardium (r = 0.83) and good at CS (r = 0.68) and RV endocardium (r = 0.58). Conclusion Non-invasive assessment of biventricular 3D activation using the 12-lead ECG and MR imaging is feasible. Potential applications include patient-specific modelling and pre-/per-procedural evaluation of ventricular activation.
Collapse
Affiliation(s)
- Simone Pezzuto
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - Frits W Prinzen
- Department of Physiology, CARIM, Maastricht University, Maastricht, The Netherlands
| | - Mark Potse
- University of Bordeaux, IMB, UMR 5251, Talence, France.,CARMEN Research Team, Inria Bordeaux - Sud-Ouest, Talence, France.,IHU Liryc, Fondation Bordeaux Université, Pessac, France
| | - Francesco Maffessanti
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - François Regoli
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland.,Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Maria Luce Caputo
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland.,Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Giulio Conte
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland.,Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Rolf Krause
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - Angelo Auricchio
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland.,Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| |
Collapse
|
7
|
Vincent KP, Forsch N, Govil S, Joblon JM, Omens JH, Perry JC, McCulloch AD. Atlas-based methods for efficient characterization of patient-specific ventricular activation patterns. Europace 2021; 23:i88-i95. [PMID: 33751079 DOI: 10.1093/europace/euaa397] [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: 11/23/2020] [Accepted: 12/03/2020] [Indexed: 11/15/2022] Open
Abstract
AIMS Ventricular activation patterns can aid clinical decision-making directly by providing spatial information on cardiac electrical activation or indirectly through derived clinical indices. The aim of this work was to derive an atlas of the major modes of variation of ventricular activation from model-predicted 3D bi-ventricular activation time distributions and to relate these modes to corresponding vectorcardiograms (VCGs). We investigated how the resulting dimensionality reduction can improve and accelerate the estimation of activation patterns from surface electrogram measurements. METHODS AND RESULTS Atlases of activation time (AT) and VCGs were derived using principal component analysis on a dataset of simulated electrophysiology simulations computed on eight patient-specific bi-ventricular geometries. The atlases provided significant dimensionality reduction, and the modes of variation in the two atlases described similar features. Utility of the atlases was assessed by resolving clinical waveforms against them and the VCG atlas was able to accurately reconstruct the patient VCGs with fewer than 10 modes. A sensitivity analysis between the two atlases was performed by calculating a compact Jacobian. Finally, VCGs generated by varying AT atlas modes were compared with clinical VCGs to estimate patient-specific activation maps, and the resulting errors between the clinical and atlas-based VCGs were less than those from more computationally expensive method. CONCLUSION Atlases of activation and VCGs represent a new method of identifying and relating the features of these high-dimensional signals that capture the major sources of variation between patients and may aid in identifying novel clinical indices of arrhythmia risk or therapeutic outcome.
Collapse
Affiliation(s)
- Kevin P Vincent
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA
| | - Nickolas Forsch
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA
| | - Sachin Govil
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA
| | - Jake M Joblon
- Department of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Jeffrey H Omens
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA.,Department of Medicine, University of California San Diego, La Jolla, CA, USA
| | - James C Perry
- Department of Pediatrics, University of California San Diego, La Jolla, CA, USA.,Rady Children's Hospital, San Diego, CA, USA
| | - Andrew D McCulloch
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA.,Department of Medicine, University of California San Diego, La Jolla, CA, USA
| |
Collapse
|
8
|
Grandits T, Effland A, Pock T, Krause R, Plank G, Pezzuto S. GEASI: Geodesic-based earliest activation sites identification in cardiac models. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3505. [PMID: 34170082 PMCID: PMC8459297 DOI: 10.1002/cnm.3505] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 06/19/2021] [Accepted: 06/22/2021] [Indexed: 05/18/2023]
Abstract
The identification of the initial ventricular activation sequence is a critical step for the correct personalization of patient-specific cardiac models. In healthy conditions, the Purkinje network is the main source of the electrical activation, but under pathological conditions the so-called earliest activation sites (EASs) are possibly sparser and more localized. Yet, their number, location and timing may not be easily inferred from remote recordings, such as the epicardial activation or the 12-lead electrocardiogram (ECG), due to the underlying complexity of the model. In this work, we introduce GEASI (Geodesic-based Earliest Activation Sites Identification) as a novel approach to simultaneously identify all EASs. To this end, we start from the anisotropic eikonal equation modeling cardiac electrical activation and exploit its Hamilton-Jacobi formulation to minimize a given objective function, for example, the quadratic mismatch to given activation measurements. This versatile approach can be extended to estimate the number of activation sites by means of the topological gradient, or fitting a given ECG. We conducted various experiments in 2D and 3D for in-silico models and an in-vivo intracardiac recording collected from a patient undergoing cardiac resynchronization therapy. The results demonstrate the clinical applicability of GEASI for potential future personalized models and clinical intervention.
Collapse
Affiliation(s)
- Thomas Grandits
- Institute of Computer Graphics and VisionTU GrazGrazAustria
- BioTechMed‐GrazGrazAustria
| | - Alexander Effland
- Institute of Computer Graphics and VisionTU GrazGrazAustria
- Silicon Austria Labs (TU Graz SAL DES Lab)GrazAustria
- Institute for Applied MathematicsUniversity of BonnBonnGermany
| | - Thomas Pock
- Institute of Computer Graphics and VisionTU GrazGrazAustria
- BioTechMed‐GrazGrazAustria
| | - Rolf Krause
- Center for Computational Medicine in Cardiology, Euler InstituteUniversità della Svizzera ItalianaLuganoSwitzerland
| | - Gernot Plank
- BioTechMed‐GrazGrazAustria
- Gottfried Schatz Research Center—Division of BiophysicsMedical University of GrazGrazAustria
| | - Simone Pezzuto
- Center for Computational Medicine in Cardiology, Euler InstituteUniversità della Svizzera ItalianaLuganoSwitzerland
| |
Collapse
|
9
|
Abstract
Computer modeling of the electrophysiology of the heart has undergone significant progress. A healthy heart can be modeled starting from the ion channels via the spread of a depolarization wave on a realistic geometry of the human heart up to the potentials on the body surface and the ECG. Research is advancing regarding modeling diseases of the heart. This article reviews progress in calculating and analyzing the corresponding electrocardiogram (ECG) from simulated depolarization and repolarization waves. First, we describe modeling of the P-wave, the QRS complex and the T-wave of a healthy heart. Then, both the modeling and the corresponding ECGs of several important diseases and arrhythmias are delineated: ischemia and infarction, ectopic beats and extrasystoles, ventricular tachycardia, bundle branch blocks, atrial tachycardia, flutter and fibrillation, genetic diseases and channelopathies, imbalance of electrolytes and drug-induced changes. Finally, we outline the potential impact of computer modeling on ECG interpretation. Computer modeling can contribute to a better comprehension of the relation between features in the ECG and the underlying cardiac condition and disease. It can pave the way for a quantitative analysis of the ECG and can support the cardiologist in identifying events or non-invasively localizing diseased areas. Finally, it can deliver very large databases of reliably labeled ECGs as training data for machine learning.
Collapse
|
10
|
Kunisch K, Trautmann P. An Inverse Problem Involving a Viscous Eikonal Equation with Applications in Electrophysiology. VIETNAM JOURNAL OF MATHEMATICS 2021; 50:301-317. [PMID: 34901261 PMCID: PMC8632858 DOI: 10.1007/s10013-021-00509-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Accepted: 02/24/2021] [Indexed: 06/14/2023]
Abstract
In this work we discuss the reconstruction of cardiac activation instants based on a viscous Eikonal equation from boundary observations. The problem is formulated as a least squares problem and solved by a projected version of the Levenberg-Marquardt method. Moreover, we analyze the well-posedness of the state equation and derive the gradient of the least squares functional with respect to the activation instants. In the numerical examples we also conduct an experiment in which the location of the activation sites and the activation instants are reconstructed jointly based on an adapted version of the shape gradient method from (J. Math. Biol. 79, 2033-2068, 2019). We are able to reconstruct the activation instants as well as the locations of the activations with high accuracy relative to the noise level.
Collapse
Affiliation(s)
- Karl Kunisch
- Institute for Mathematics and Scientific Computing, University of Graz, Heinrichstrasse 36, A-8010 Graz, Austria
- Johann Radon Institute for Computational and Applied Mathematics, Austrian Academy of Science, Linz, Austria
| | - Philip Trautmann
- Institute for Mathematics and Scientific Computing, University of Graz, Heinrichstrasse 36, A-8010 Graz, Austria
| |
Collapse
|
11
|
Pagani S, Dede’ L, Manzoni A, Quarteroni A. Data integration for the numerical simulation of cardiac electrophysiology. Pacing Clin Electrophysiol 2021; 44:726-736. [PMID: 33594761 PMCID: PMC8252775 DOI: 10.1111/pace.14198] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Revised: 01/26/2021] [Accepted: 02/07/2021] [Indexed: 12/20/2022]
Abstract
The increasing availability of extensive and accurate clinical data is rapidly shaping cardiovascular care by improving the understanding of physiological and pathological mechanisms of the cardiovascular system and opening new frontiers in designing therapies and interventions. In this direction, mathematical and numerical models provide a complementary relevant tool, able not only to reproduce patient-specific clinical indicators but also to predict and explore unseen scenarios. With this goal, clinical data are processed and provided as inputs to the mathematical model, which quantitatively describes the physical processes that occur in the cardiac tissue. In this paper, the process of integration of clinical data and mathematical models is discussed. Some challenges and contributions in the field of cardiac electrophysiology are reported.
Collapse
Affiliation(s)
- Stefano Pagani
- MOX‐Department of MathematicsPolitecnico di MilanoMilanItaly
| | - Luca Dede’
- MOX‐Department of MathematicsPolitecnico di MilanoMilanItaly
| | - Andrea Manzoni
- MOX‐Department of MathematicsPolitecnico di MilanoMilanItaly
| | - Alfio Quarteroni
- MOX‐Department of MathematicsPolitecnico di MilanoMilanItaly
- Institute of MathematicsEPFLLausanneSwitzerland
| |
Collapse
|
12
|
Auricchio A, Özkartal T, Salghetti F, Neumann L, Pezzuto S, Gharaviri A, Demarchi A, Caputo ML, Regoli F, De Asmundis C, Chierchia GB, Brugada P, Klersy C, Moccetti T, Schotten U, Conte G. Short P-Wave Duration is a Marker of Higher Rate of Atrial Fibrillation Recurrences after Pulmonary Vein Isolation: New Insights into the Pathophysiological Mechanisms Through Computer Simulations. J Am Heart Assoc 2021; 10:e018572. [PMID: 33410337 PMCID: PMC7955300 DOI: 10.1161/jaha.120.018572] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Background Short ECG P-wave duration has recently been demonstrated to be associated with higher risk of atrial fibrillation (AF). The aim of this study was to assess the rate of AF recurrence after pulmonary vein isolation in patients with a short P wave, and to mechanistically elucidate the observation by computer modeling. Methods and Results A total of 282 consecutive patients undergoing a first single-pulmonary vein isolation procedure for paroxysmal or persistent AF were included. Computational models studied the effect of adenosine and sodium conductance on action potential duration and P-wave duration (PWD). About 16% of the patients had a PWD of 110 ms or shorter (median PWD 126 ms, interquartile range, 115 ms-138 ms; range, 71 ms-180 ms). At Cox regression, PWD was significantly associated with AF recurrence (P=0.012). Patients with a PWD <110 ms (hazard ratio [HR], 2.20; 95% CI, 1.24-3.88; P=0.007) and patients with a PWD ≥140 (HR, 1.87, 95% CI, 1.06-3.30; P=0.031) had a nearly 2-fold increase in risk with respect to the other group. In the computational model, adenosine yielded a significant reduction of action potential duration 90 (52%) and PWD (7%). An increased sodium conductance (up to 200%) was robustly accompanied by an increase in conduction velocity (26%), a reduction in action potential duration 90 (28%), and PWD (22%). Conclusions One out of 5 patients referred for pulmonary vein isolation has a short PWD which was associated with a higher rate of AF after the index procedure. Computer simulations suggest that shortening of atrial action potential duration leading to a faster atrial conduction may be the cause of this clinical observation.
Collapse
Affiliation(s)
- Angelo Auricchio
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland.,Center for Computational Medicine in Cardiology Università della Svizzera Italiana Lugano Switzerland
| | - Tardu Özkartal
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland.,Department of Cardiology Ospedale San Giovanni Bellinzona Switzerland
| | - Francesca Salghetti
- Heart Rhythm Management Centre Postgraduate course in Cardiac Electrophysiology and PacingVrije Universiteit BrusselUniversitair Ziekenhuis Brussel Brussels Belgium
| | - Laura Neumann
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland
| | - Simone Pezzuto
- Center for Computational Medicine in Cardiology Università della Svizzera Italiana Lugano Switzerland
| | - Ali Gharaviri
- Center for Computational Medicine in Cardiology Università della Svizzera Italiana Lugano Switzerland
| | - Andrea Demarchi
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland
| | | | - François Regoli
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland
| | - Carlo De Asmundis
- Heart Rhythm Management Centre Postgraduate course in Cardiac Electrophysiology and PacingVrije Universiteit BrusselUniversitair Ziekenhuis Brussel Brussels Belgium
| | - Gian-Battista Chierchia
- Heart Rhythm Management Centre Postgraduate course in Cardiac Electrophysiology and PacingVrije Universiteit BrusselUniversitair Ziekenhuis Brussel Brussels Belgium
| | - Pedro Brugada
- Heart Rhythm Management Centre Postgraduate course in Cardiac Electrophysiology and PacingVrije Universiteit BrusselUniversitair Ziekenhuis Brussel Brussels Belgium
| | - Catherine Klersy
- Service of Clinical Epidemiology & Biometry Fondazione IRCCS Policlinico San Matteo Pavia Italy
| | | | - Ulrich Schotten
- Department of Physiology Cardiovascular Research Institute MaastrichtMaastricht University Maastricht The Netherlands
| | - Giulio Conte
- Division of Cardiology Cardiocentro Ticino Lugano Switzerland
| |
Collapse
|
13
|
Grandits T, Gillette K, Neic A, Bayer J, Vigmond E, Pock T, Plank G. An Inverse Eikonal Method for Identifying Ventricular Activation Sequences from Epicardial Activation Maps. JOURNAL OF COMPUTATIONAL PHYSICS 2020; 419:109700. [PMID: 32952215 PMCID: PMC7116090 DOI: 10.1016/j.jcp.2020.109700] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
A key mechanism controlling cardiac function is the electrical activation sequence of the heart's main pumping chambers termed the ventricles. As such, personalization of the ventricular activation sequences is of pivotal importance for the clinical utility of computational models of cardiac electrophysiology. However, a direct observation of the activation sequence throughout the ventricular volume is virtually impossible. In this study, we report on a novel method for identification of activation sequences from activation maps measured at the outer surface of the heart termed the epicardium. Conceptually, the method attempts to identify the key factors governing the ventricular activation sequence - the timing of earliest activation sites (EAS) and the velocity tensor field within the ventricular walls - from sparse and noisy activation maps sampled from the epicardial surface and fits an Eikonal model to the observations. Regularization methods are first investigated to overcome the severe ill-posedness of the inverse problem in a simplified 2D example. These methods are then employed in an anatomically accurate biventricular model with two realistic activation models of varying complexity - a simplified trifascicular model (3F) and a topologically realistic model of the His-Purkinje system (HPS). Using epicardial activation maps at full resolution, we first demonstrate that reconstructing the volumetric activation sequence is, in principle, feasible under the assumption of known location of EAS and later evaluate robustness of the method against noise and reduced spatial resolution of observations. Our results suggest that the FIMIN algorithm is able to robustly recover the full 3D activation sequence using epicardial activation maps at a spatial resolution achievable with current mapping systems and in the presence of noise. Comparing the accuracy achieved in the reconstructed activation maps with clinical data uncertainties suggests that the FIMIN method may be suitable for the patient- specific parameterization of activation models.
Collapse
Affiliation(s)
- Thomas Grandits
- Institute of Computer Graphics and Vision, Graz University of Technology
- BioTechMed-Graz, Austria
| | - Karli Gillette
- Institute of Biophysics, Medical University of Graz
- BioTechMed-Graz, Austria
| | - Aurel Neic
- Institute of Biophysics, Medical University of Graz
| | - Jason Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux
| | - Edward Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux
| | - Thomas Pock
- Institute of Computer Graphics and Vision, Graz University of Technology
- BioTechMed-Graz, Austria
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz
- BioTechMed-Graz, Austria
| |
Collapse
|
14
|
Ganellari D, Haase G, Zumbusch G, Lotz J, Peltzer P, Leppkes K, Naumann U. Patient-Specific Cardiac Parametrization from Eikonal Simulations. LECTURE NOTES IN COMPUTER SCIENCE 2020. [PMCID: PMC7302314 DOI: 10.1007/978-3-030-50371-0_21] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
Simulations in cardiac electrophysiology use the bidomain equations to describe the electrical potential in the heart. If only the electrical activation sequence in the heart is needed, then the full bidomain equations can be substituted by the Eikonal equation which allows much faster responses w.r.t. the changed material parameters in the equation. We use our Eikonal solver optimized for memory usage and parallelization. Patient-specific simulations in cardiac electrophysiology require patient-specific conductivity parameters which are not accurately available in vivo. One chance to improve the given conductivity parameters consists in comparing the computed activation sequence on the heart surface with the measured ECG on the torso mapped onto this surface. By minimizing the squared distance between the measured solution and the Eikonal computed solution we are able to determine the material parameters more accurately. To reduce the number of optimization parameters in this process, we group the material parameters and introduce a specific scaling parameter \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\gamma _k$$\end{document} for each group. The minimization takes place w.r.t. the scaling \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\gamma }$$\end{document}. We solve the minimization problem by the BFGS method and adaptive step size control. The required gradient \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla _\gamma f({\gamma })$$\end{document} is computed either via finite differences or algorithmic differentiation using
in tangent as well as in adjoint mode. We present convergence behavior as well as runtime and scaling results.
Collapse
|
15
|
Loewe A, Poremba E, Oesterlein T, Luik A, Schmitt C, Seemann G, Dössel O. Patient-Specific Identification of Atrial Flutter Vulnerability-A Computational Approach to Reveal Latent Reentry Pathways. Front Physiol 2019; 9:1910. [PMID: 30692934 PMCID: PMC6339942 DOI: 10.3389/fphys.2018.01910] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2018] [Accepted: 12/18/2018] [Indexed: 11/23/2022] Open
Abstract
Atypical atrial flutter (AFlut) is a reentrant arrhythmia which patients frequently develop after ablation for atrial fibrillation (AF). Indeed, substrate modifications during AF ablation can increase the likelihood to develop AFlut and it is clinically not feasible to reliably and sensitively test if a patient is vulnerable to AFlut. Here, we present a novel method based on personalized computational models to identify pathways along which AFlut can be sustained in an individual patient. We build a personalized model of atrial excitation propagation considering the anatomy as well as the spatial distribution of anisotropic conduction velocity and repolarization characteristics based on a combination of a priori knowledge on the population level and information derived from measurements performed in the individual patient. The fast marching scheme is employed to compute activation times for stimuli from all parts of the atria. Potential flutter pathways are then identified by tracing loops from wave front collision sites and constricting them using a geometric snake approach under consideration of the heterogeneous wavelength condition. In this way, all pathways along which AFlut can be sustained are identified. Flutter pathways can be instantiated by using an eikonal-diffusion phase extrapolation approach and a dynamic multifront fast marching simulation. In these dynamic simulations, the initial pattern eventually turns into the one driven by the dominant pathway, which is the only pathway that can be observed clinically. We assessed the sensitivity of the flutter pathway maps with respect to conduction velocity and its anisotropy. Moreover, we demonstrate the application of tailored models considering disease-specific repolarization properties (healthy, AF-remodeled, potassium channel mutations) as well as applicabiltiy on a clinical dataset. Finally, we tested how AFlut vulnerability of these substrates is modulated by exemplary antiarrhythmic drugs (amiodarone, dronedarone). Our novel method allows to assess the vulnerability of an individual patient to develop AFlut based on the personal anatomical, electrophysiological, and pharmacological characteristics. In contrast to clinical electrophysiological studies, our computational approach provides the means to identify all possible AFlut pathways and not just the currently dominant one. This allows to consider all relevant AFlut pathways when tailoring clinical ablation therapy in order to reduce the development and recurrence of AFlut.
Collapse
Affiliation(s)
- Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Emanuel Poremba
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Tobias Oesterlein
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Claus Schmitt
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Gunnar Seemann
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg Bad Krozingen, Freiburg, Germany
- Faculty of Medicine, Albert-Ludwigs University, Freiburg, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
16
|
Pezzuto S, Gharaviri A, Schotten U, Potse M, Conte G, Caputo ML, Regoli F, Krause R, Auricchio A. Beat-to-beat P-wave morphological variability in patients with paroxysmal atrial fibrillation: anin silicostudy. Europace 2018; 20:iii26-iii35. [DOI: 10.1093/europace/euy227] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2018] [Accepted: 09/19/2018] [Indexed: 11/12/2022] Open
Affiliation(s)
- Simone Pezzuto
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - Ali Gharaviri
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - Ulrich Schotten
- Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, the Netherlands
| | - Mark Potse
- CARMEN Research Team, INRIA, Talence, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Foundation Bordeaux Université, Pessac, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Giulio Conte
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Maria Luce Caputo
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Francois Regoli
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Rolf Krause
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Institute of Computational Science, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
| | - Angelo Auricchio
- Center for Computational Medicine in Cardiology, Università della Svizzera italiana, Lugano, Switzerland
- Fondazione Cardiocentro Ticino, Lugano, Switzerland
| |
Collapse
|
17
|
Siregar P, Julen N, Hufnagl P, Mutter G. A general framework dedicated to computational morphogenesis Part I - Constitutive equations. Biosystems 2018; 173:298-313. [PMID: 30005999 DOI: 10.1016/j.biosystems.2018.07.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2018] [Revised: 05/30/2018] [Accepted: 07/05/2018] [Indexed: 01/14/2023]
Abstract
In order to understand living organisms, considerable experimental efforts and resources have been devoted to correlate genes and their expressions with cell, tissue, organ and whole organisms' phenotypes. This data driven approach to knowledge discovery has led to many breakthrough in our understanding of healthy and diseased states, and is paving the way to improve the diagnosis and treatment of diseases. Complementary to this data-driven approach, computational models of biological systems based on first principles have been developed in order to deepen our understanding of the multi-scale dynamics that drives normal and pathological biological functions. In this paper we describe the biological, physical and mathematical concepts that led to the design of a Computational Morphogenesis (CM) platform baptized Generic Modeling and Simulating Platform (GMSP). Its role is to generate realistic 3D multi-scale biological tissues from virtual stem cells and the intended target applications include in virtuo studies of normal and abnormal tissue (re)generation as well as the development of complex diseases such as carcinogenesis. At all space-scales of interest, biological agents interact with each other via biochemical, bioelectrical, and mechanical fields that operate in concert during embryogenesis, growth and adult life. The spatio-temporal dependencies of these fields can be modeled by physics-based constitutive equations that we propose to examine in relation to the landmark biological events that occur during embryogenesis.
Collapse
Affiliation(s)
| | | | - Peter Hufnagl
- Department of Digital Pathology and IT, Institute of Pathology, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - George Mutter
- Department of Pathology, Harvard Medical School and Brigham and Women's Hospital, Boston, MA, USA
| |
Collapse
|
18
|
Quaglino A, Pezzuto S, Koutsourelakis PS, Auricchio A, Krause R. Fast uncertainty quantification of activation sequences in patient-specific cardiac electrophysiology meeting clinical time constraints. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e2985. [PMID: 29577657 DOI: 10.1002/cnm.2985] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2017] [Revised: 01/16/2018] [Accepted: 03/15/2018] [Indexed: 06/08/2023]
Abstract
We present a fast, patient-specific methodology for uncertainty quantification in electrophysiology, aimed at meeting the time constraints of clinical practitioners. We focus on computing the statistics of the activation map, given the uncertainties associated with the conductivity tensor modeling the fiber orientation in the heart. We use a fast parallel solution method implemented on a graphics processing unit for the eikonal approximation, in order to compute the activation map and to sample the random fiber field with correlation on the basis of geodesic distances. While this enables to perform uncertainty quantification studies with a manageable computational effort, the required time frame still exceeds clinically suitable time expectations. In order to reduce it further by 2 orders of magnitude, we rely on Bayesian multifidelity methods. In particular, we propose a low-fidelity model that is patient-specific and free from the additional training cost associated with reduced models. This is achieved by a sound physics-based simplification of the full eikonal model. The low-fidelity output is then corrected by the standard multifidelity framework. In practice, the complete procedure only requires approximately 100 new runs of our eikonal graphics processing unit solver for producing the sought estimates and their associated credible intervals, enabling a full online analysis in less than 5 minutes.
Collapse
Affiliation(s)
- A Quaglino
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | - S Pezzuto
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | | | - A Auricchio
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
- Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - R Krause
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| |
Collapse
|
19
|
Potse M. Scalable and Accurate ECG Simulation for Reaction-Diffusion Models of the Human Heart. Front Physiol 2018; 9:370. [PMID: 29731720 PMCID: PMC5920200 DOI: 10.3389/fphys.2018.00370] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2018] [Accepted: 03/27/2018] [Indexed: 11/13/2022] Open
Abstract
Realistic electrocardiogram (ECG) simulation with numerical models is important for research linking cellular and molecular physiology to clinically observable signals, and crucial for patient tailoring of numerical heart models. However, ECG simulation with a realistic torso model is computationally much harder than simulation of cardiac activity itself, so that many studies with sophisticated heart models have resorted to crude approximations of the ECG. This paper shows how the classical concept of electrocardiographic lead fields can be used for an ECG simulation method that matches the realism of modern heart models. The accuracy and resource requirements were compared to those of a full-torso solution for the potential and scaling was tested up to 14,336 cores with a heart model consisting of 11 million nodes. Reference ECGs were computed on a 3.3 billion-node heart-torso mesh at 0.2 mm resolution. The results show that the lead-field method is more efficient than a full-torso solution when the number of simulated samples is larger than the number of computed ECG leads. While the initial computation of the lead fields remains a hard and poorly scalable problem, the ECG computation itself scales almost perfectly and, even for several hundreds of ECG leads, takes much less time than the underlying simulation of cardiac activity.
Collapse
Affiliation(s)
- Mark Potse
- CARMEN Research Team, Inria Bordeaux Sud-Ouest, Talence, France.,Institut de Mathématiques de Bordeaux, UMR 5251, Université de Bordeaux, Talence, France.,IHU Liryc, Electrophysiology and Heart Modeling Institute, Foundation Bordeaux Université, Pessac-Bordeaux, France
| |
Collapse
|