1
|
Trayanova NA, Lyon A, Shade J, Heijman J. Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation. Physiol Rev 2024; 104:1265-1333. [PMID: 38153307 DOI: 10.1152/physrev.00017.2023] [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/05/2023] [Revised: 12/19/2023] [Accepted: 12/21/2023] [Indexed: 12/29/2023] Open
Abstract
The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland, United States
| | - Aurore Lyon
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
- Division of Heart and Lungs, Department of Medical Physiology, University Medical Centre Utrecht, Utrecht, The Netherlands
| | - Julie Shade
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland, United States
| | - Jordi Heijman
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
2
|
Wülfers EM, Moss R, Lehrmann H, Arentz T, Westermann D, Seemann G, Odening KE, Steinfurt J. Whole-heart computational modelling provides further mechanistic insights into ST-elevation in Brugada syndrome. INTERNATIONAL JOURNAL OF CARDIOLOGY. HEART & VASCULATURE 2024; 51:101373. [PMID: 38464963 PMCID: PMC10924145 DOI: 10.1016/j.ijcha.2024.101373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 02/08/2024] [Accepted: 02/21/2024] [Indexed: 03/12/2024]
Abstract
Background Brugada syndrome (BrS) is characterized by dynamic ST-elevations in right precordial leads and increased risk of ventricular fibrillation and sudden cardiac death. As the mechanism underlying ST-elevation and malignant arrhythmias is controversial computational modeling can aid in exploring the disease mechanism. Thus we aim to test the main competing hypotheses ('delayed depolarization' vs. 'early repolarization') of BrS in a whole-heart computational model. Methods In a 3D whole-heart computational model, delayed epicardial RVOT activation with local conduction delay was simulated by reducing conductivity in the epicardial RVOT. Early repolarization was simulated by instead increasing the transient outward potassium current (Ito) in the same region. Additionally, a reduction in the fast sodium current (INa) was incorporated in both models. Results Delayed depolarization with local conduction delay in the computational model resulted in coved-type ST-elevation with negative T-waves in the precordial surface ECG leads. 'Saddleback'-shaped ST-elevation was obtained with reduced substrate extent or thickness. Increased Ito simulations showed early repolarization in the RVOT with a descending but not coved-type ST-elevation. Reduced INa did not show a significant effect on ECG morphology. Conclusions In this whole-heart BrS computational model of both major hypotheses, realistic coved-type ECG resulted only from delayed epicardial RVOT depolarization with local conduction delay but not early repolarizing ion channel modifications. These simulations provide further support for the depolarization hypothesis as electrophysiological mechanism underlying BrS.
Collapse
Affiliation(s)
- Eike M Wülfers
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Department of Physics and Astronomy, Faculty of Sciences, Ghent University, Ghent, Belgium
| | - Robin Moss
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Heiko Lehrmann
- Department of Cardiology and Angiology, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Thomas Arentz
- Department of Cardiology and Angiology, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Dirk Westermann
- Department of Cardiology and Angiology, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Gunnar Seemann
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Katja E Odening
- Translational Cardiology, Department of Cardiology and Institute of Physiology, University Hospital Bern, University of Bern, Switzerland
| | - Johannes Steinfurt
- Department of Cardiology and Angiology, University Heart Center Freiburg - Bad Krozingen, and Faculty of Medicine, University of Freiburg, Freiburg, Germany
| |
Collapse
|
3
|
Gillette K, Gsell MAF, Nagel C, Bender J, Winkler B, Williams SE, Bär M, Schäffter T, Dössel O, Plank G, Loewe A. MedalCare-XL: 16,900 healthy and pathological synthetic 12 lead ECGs from electrophysiological simulations. Sci Data 2023; 10:531. [PMID: 37553349 PMCID: PMC10409805 DOI: 10.1038/s41597-023-02416-4] [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: 03/29/2023] [Accepted: 07/25/2023] [Indexed: 08/10/2023] Open
Abstract
Mechanistic cardiac electrophysiology models allow for personalized simulations of the electrical activity in the heart and the ensuing electrocardiogram (ECG) on the body surface. As such, synthetic signals possess known ground truth labels of the underlying disease and can be employed for validation of machine learning ECG analysis tools in addition to clinical signals. Recently, synthetic ECGs were used to enrich sparse clinical data or even replace them completely during training leading to improved performance on real-world clinical test data. We thus generated a novel synthetic database comprising a total of 16,900 12 lead ECGs based on electrophysiological simulations equally distributed into healthy control and 7 pathology classes. The pathological case of myocardial infraction had 6 sub-classes. A comparison of extracted features between the virtual cohort and a publicly available clinical ECG database demonstrated that the synthetic signals represent clinical ECGs for healthy and pathological subpopulations with high fidelity. The ECG database is split into training, validation, and test folds for development and objective assessment of novel machine learning algorithms.
Collapse
Affiliation(s)
- Karli Gillette
- Gottfried Schatz Research Center: Division of Medical Physics and Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Matthias A F Gsell
- Gottfried Schatz Research Center: Division of Medical Physics and Biophysics, Medical University of Graz, Graz, Austria
| | - Claudia Nagel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Jule Bender
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Benjamin Winkler
- Physikalisch-Technische Bundesanstalt, National Metrology Institute, Berlin, Germany
| | - Steven E Williams
- King's College London, London, United Kingdom
- University of Edinburgh, Edinburgh, United Kingdom
| | - Markus Bär
- Physikalisch-Technische Bundesanstalt, National Metrology Institute, Berlin, Germany
| | - Tobias Schäffter
- Physikalisch-Technische Bundesanstalt, National Metrology Institute, Berlin, Germany
- King's College London, London, United Kingdom
- Biomedical Engineering, Technische Universität Berlin, Einstein Centre Digital Future, Berlin, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Gernot Plank
- Gottfried Schatz Research Center: Division of Medical Physics and Biophysics, Medical University of Graz, Graz, Austria.
- BioTechMed-Graz, Graz, Austria.
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.
| |
Collapse
|
4
|
Rinné S, Oertli A, Nagel C, Tomsits P, Jenewein T, Kääb S, Kauferstein S, Loewe A, Beckmann BM, Decher N. Functional Characterization of a Spectrum of Novel Romano-Ward Syndrome KCNQ1 Variants. Int J Mol Sci 2023; 24:ijms24021350. [PMID: 36674868 PMCID: PMC9865342 DOI: 10.3390/ijms24021350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Revised: 12/20/2022] [Accepted: 12/23/2022] [Indexed: 01/13/2023] Open
Abstract
The KCNQ1 gene encodes the α-subunit of the cardiac voltage-gated potassium (Kv) channel KCNQ1, also denoted as Kv7.1 or KvLQT1. The channel assembles with the ß-subunit KCNE1, also known as minK, to generate the slowly activating cardiac delayed rectifier current IKs, a key regulator of the heart rate dependent adaptation of the cardiac action potential duration (APD). Loss-of-function variants in KCNQ1 cause the congenital Long QT1 (LQT1) syndrome, characterized by delayed cardiac repolarization and a QT interval prolongation in the surface electrocardiogram (ECG). Autosomal dominant loss-of-function variants in KCNQ1 result in the LQT syndrome called Romano-Ward syndrome (RWS), while autosomal recessive variants affecting function, lead to Jervell and Lange-Nielsen syndrome (JLNS), associated with deafness. The aim of this study was the characterization of novel KCNQ1 variants identified in patients with RWS to widen the spectrum of known LQT1 variants, and improve the interpretation of the clinical relevance of variants in the KCNQ1 gene. We functionally characterized nine human KCNQ1 variants using the voltage-clamp technique in Xenopus laevis oocytes, from which we report seven novel variants. The functional data was taken as input to model surface ECGs, to subsequently compare the functional changes with the clinically observed QTc times, allowing a further interpretation of the severity of the different LQTS variants. We found that the electrophysiological properties of the variants correlate with the severity of the clinically diagnosed phenotype in most cases, however, not in all. Electrophysiological studies combined with in silico modelling approaches are valuable components for the interpretation of the pathogenicity of KCNQ1 variants, but assessing the clinical severity demands the consideration of other factors that are included, for example in the Schwartz score.
Collapse
Affiliation(s)
- Susanne Rinné
- Institute of Physiology and Pathophysiology, Vegetative Physiology, University of Marburg, 35037 Marburg, Germany
| | - Annemarie Oertli
- Institute of Physiology and Pathophysiology, Vegetative Physiology, University of Marburg, 35037 Marburg, Germany
| | - Claudia Nagel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Philipp Tomsits
- Department of Medicine I, University Hospital, LMU Munich, 80802 Munich, Germany
- Deutsches Zentrum für Herz-Kreislauferkrankungen (DZHK), Partner Site Munich, 80636 Munich, Germany
- Member of the European Reference Network for Rare, Low Prevalance and Complex Diseases of the Heart (ERN GUARD-Heart), 81377 Munich, Germany
- Institute of Surgical Research at the Walter-Brendel-Centre of Experimental Medicine, University Hospital, LMU Munich, Marchioninistrasse 27, 81377 Munich, Germany
| | - Tina Jenewein
- Institute of Legal Medicine, Goethe University, University Hospital Frankfurt, 60590 Frankfurt, Germany
- Institute for Transfusion Medicine and Immunohematology, German Red Cross Blood Service Baden-Württemberg-Hessen, Goethe University Frankfurt, 60528 Frankfurt, Germany
| | - Stefan Kääb
- Department of Medicine I, University Hospital, LMU Munich, 80802 Munich, Germany
- Deutsches Zentrum für Herz-Kreislauferkrankungen (DZHK), Partner Site Munich, 80636 Munich, Germany
- Member of the European Reference Network for Rare, Low Prevalance and Complex Diseases of the Heart (ERN GUARD-Heart), 81377 Munich, Germany
| | - Silke Kauferstein
- Institute of Legal Medicine, Goethe University, University Hospital Frankfurt, 60590 Frankfurt, Germany
- Deutsches Zentrum für Herz-Kreislauferkrankungen (DZHK), Partner Site Frankfurt, 60596 Frankfurt, Germany
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Britt Maria Beckmann
- Department of Medicine I, University Hospital, LMU Munich, 80802 Munich, Germany
- Institute of Legal Medicine, Goethe University, University Hospital Frankfurt, 60590 Frankfurt, Germany
| | - Niels Decher
- Institute of Physiology and Pathophysiology, Vegetative Physiology, University of Marburg, 35037 Marburg, Germany
- Correspondence: ; Tel.: +49-(0)6421-28-62148
| |
Collapse
|
5
|
Moss R, Wülfers EM, Lewetag R, Hornyik T, Perez-Feliz S, Strohbach T, Menza M, Krafft A, Odening KE, Seemann G. A computational model of rabbit geometry and ECG: Optimizing ventricular activation sequence and APD distribution. PLoS One 2022; 17:e0270559. [PMID: 35771854 PMCID: PMC9246225 DOI: 10.1371/journal.pone.0270559] [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: 07/12/2021] [Accepted: 06/13/2022] [Indexed: 11/19/2022] Open
Abstract
Computational modeling of electrophysiological properties of the rabbit heart is a commonly used way to enhance and/or complement findings from classic lab work on single cell or tissue levels. Yet, thus far, there was no possibility to extend the scope to include the resulting body surface potentials as a way of validation or to investigate the effect of certain pathologies. Based on CT imaging, we developed the first openly available computational geometrical model not only of the whole heart but also the complete torso of the rabbit. Additionally, we fabricated a 32-lead ECG-vest to record body surface potential signals of the aforementioned rabbit. Based on the developed geometrical model and the measured signals, we then optimized the activation sequence of the ventricles, recreating the functionality of the Purkinje network, and we investigated different apico-basal and transmural gradients in action potential duration. Optimization of the activation sequence resulted in an average root mean square error between measured and simulated signal of 0.074 mV/ms for all leads. The best-fit T-Wave, compared to measured data (0.038 mV/ms), resulted from incorporating an action potential duration gradient from base to apex with a respective shortening of 20 ms and a transmural gradient with a shortening of 15 ms from endocardium to epicardium. By making our model and measured data openly available, we hope to give other researchers the opportunity to verify their research, as well as to create the possibility to investigate the impact of electrophysiological alterations on body surface signals for translational research.
Collapse
Affiliation(s)
- Robin Moss
- 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
- * E-mail:
| | - Eike M. Wülfers
- 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
| | - Raphaela Lewetag
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg ⋅ Bad Krozingen, Medical Center—University of Freiburg, Freiburg, Germany
- Department of Cardiology and Angiology I, Heart Center University of Freiburg, Medical Faculty, Freiburg, Germany
| | - Tibor Hornyik
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg ⋅ Bad Krozingen, Medical Center—University of Freiburg, Freiburg, Germany
- Translational Cardiology, Department of Cardiology and Department of Physiology, University Hospital Bern, Bern, Switzerland
| | - Stefanie Perez-Feliz
- 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
- Department of Cardiology and Angiology I, Heart Center University of Freiburg, Medical Faculty, Freiburg, Germany
| | - Tim Strohbach
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Freiburg, Germany
| | - Marius Menza
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Freiburg, Germany
| | - Axel Krafft
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Freiburg, Germany
| | - Katja E. Odening
- Department of Cardiology and Angiology I, Heart Center University of Freiburg, Medical Faculty, Freiburg, Germany
- Translational Cardiology, Department of Cardiology and Department of Physiology, University Hospital Bern, Bern, Switzerland
| | - 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
| |
Collapse
|
6
|
Moss R, Wülfers EM, Schuler S, Loewe A, Seemann G. A Fully-Coupled Electro-Mechanical Whole-Heart Computational Model: Influence of Cardiac Contraction on the ECG. Front Physiol 2022; 12:778872. [PMID: 34975532 PMCID: PMC8716847 DOI: 10.3389/fphys.2021.778872] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Accepted: 11/17/2021] [Indexed: 01/12/2023] Open
Abstract
The ECG is one of the most commonly used non-invasive tools to gain insights into the electrical functioning of the heart. It has been crucial as a foundation in the creation and validation of in silico models describing the underlying electrophysiological processes. However, so far, the contraction of the heart and its influences on the ECG have mainly been overlooked in in silico models. As the heart contracts and moves, so do the electrical sources within the heart responsible for the signal on the body surface, thus potentially altering the ECG. To illuminate these aspects, we developed a human 4-chamber electro-mechanically coupled whole heart in silico model and embedded it within a torso model. Our model faithfully reproduces measured 12-lead ECG traces, circulatory characteristics, as well as physiological ventricular rotation and atrioventricular valve plane displacement. We compare our dynamic model to three non-deforming ones in terms of standard clinically used ECG leads (Einthoven and Wilson) and body surface potential maps (BSPM). The non-deforming models consider the heart at its ventricular end-diastatic, end-diastolic and end-systolic states. The standard leads show negligible differences during P-Wave and QRS-Complex, yet during T-Wave the leads closest to the heart show prominent differences in amplitude. When looking at the BSPM, there are no notable differences during the P-Wave, but effects of cardiac motion can be observed already during the QRS-Complex, increasing further during the T-Wave. We conclude that for the modeling of activation (P-Wave/QRS-Complex), the associated effort of simulating a complete electro-mechanical approach is not worth the computational cost. But when looking at ventricular repolarization (T-Wave) in standard leads as well as BSPM, there are areas where the signal can be influenced by cardiac motion of the heart to an extent that should not be ignored.
Collapse
Affiliation(s)
- Robin Moss
- 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
| | - Eike Moritz Wülfers
- 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
| | - Steffen Schuler
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Axel Loewe
- 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
| |
Collapse
|
7
|
Sung E, Prakosa A, Trayanova NA. Analyzing the Role of Repolarization Gradients in Post-infarct Ventricular Tachycardia Dynamics Using Patient-Specific Computational Heart Models. Front Physiol 2021; 12:740389. [PMID: 34658925 PMCID: PMC8514757 DOI: 10.3389/fphys.2021.740389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Accepted: 09/06/2021] [Indexed: 11/20/2022] Open
Abstract
Aims: Disease-induced repolarization heterogeneity in infarcted myocardium contributes to VT arrhythmogenesis but how apicobasal and transmural (AB-TM) repolarization gradients additionally affect post-infarct VT dynamics is unknown. The goal of this study is to assess how AB-TM repolarization gradients impact post-infarct VT dynamics using patient-specific heart models. Method: 3D late gadolinium-enhanced cardiac magnetic resonance images were acquired from seven post-infarct patients. Models representing the patient-specific scar and infarct border zone distributions were reconstructed without (baseline) and with repolarization gradients along both the AB-TM axes. AB only and TM only models were created to assess the effects of each ventricular gradient on VT dynamics. VTs were induced in all models via rapid pacing. Results: Ten baseline VTs were induced. VT inducibility in AB-TM models was not significantly different from baseline (p>0.05). Reentry pathways in AB-TM models were different than baseline pathways due to alterations in the location of conduction block (p<0.05). VT exit sites in AB-TM models were different than baseline VT exit sites (p<0.05). VT inducibility of AB only and TM only models were not significantly different than that of baseline (p>0.05) or AB-TM models (p>0.05). Reentry pathways and VT exit sites in AB only and TM only models were different than in baseline (p<0.05). Lastly, repolarization gradients uncovered multiple VT morphologies with different reentrant pathways and exit sites within the same structural, conducting channels. Conclusion: VT inducibility was not impacted by the addition of AB-TM repolarization gradients, but the VT reentrant pathway and exit sites were greatly affected due to modulation of conduction block. Thus, during ablation procedures, physiological and pharmacological factors that impact the ventricular repolarization gradient might need to be considered when targeting the VTs.
Collapse
Affiliation(s)
- Eric Sung
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States.,Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, MD, United States
| | - Adityo Prakosa
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States.,Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, MD, United States
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States.,Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, MD, United States
| |
Collapse
|
8
|
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: 61] [Impact Index Per Article: 20.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
|
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
|
Zhang XD, Thai PN, Lieu DK, Chiamvimonvat N. Model Systems for Addressing Mechanism of Arrhythmogenesis in Cardiac Repair. Curr Cardiol Rep 2021; 23:72. [PMID: 34050853 PMCID: PMC8164614 DOI: 10.1007/s11886-021-01498-z] [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] [Accepted: 03/09/2021] [Indexed: 11/09/2022]
Abstract
PURPOSE OF REVIEW Cardiac cell-based therapy represents a promising approach for cardiac repair. However, one of the main challenges is cardiac arrhythmias associated with stem cell transplantation. The current review summarizes the recent progress in model systems for addressing mechanisms of arrhythmogenesis in cardiac repair. RECENT FINDINGS Animal models have been extensively developed for mechanistic studies of cardiac arrhythmogenesis. Advances in human induced pluripotent stem cells (hiPSCs), patient-specific disease models, tissue engineering, and gene editing have greatly enhanced our ability to probe the mechanistic bases of cardiac arrhythmias. Additionally, recent development in multiscale computational studies and machine learning provides yet another powerful tool to quantitatively decipher the mechanisms of cardiac arrhythmias. Advancing efforts towards the integrations of experimental and computational studies are critical to gain insights into novel mitigation strategies for cardiac arrhythmias in cell-based therapy.
Collapse
Affiliation(s)
- Xiao-Dong Zhang
- Division of Cardiovascular Medicine, Department of Internal Medicine, School of Medicine, University of California, Davis, Davis, CA 95616 USA
- Department of Veterans Affairs, Veterans Affairs Northern California Health Care System, Mather, CA 95655 USA
| | - Phung N. Thai
- Division of Cardiovascular Medicine, Department of Internal Medicine, School of Medicine, University of California, Davis, Davis, CA 95616 USA
- Department of Veterans Affairs, Veterans Affairs Northern California Health Care System, Mather, CA 95655 USA
| | - Deborah K. Lieu
- Division of Cardiovascular Medicine, Department of Internal Medicine, School of Medicine, University of California, Davis, Davis, CA 95616 USA
| | - Nipavan Chiamvimonvat
- Division of Cardiovascular Medicine, Department of Internal Medicine, School of Medicine, University of California, Davis, Davis, CA 95616 USA
- Department of Veterans Affairs, Veterans Affairs Northern California Health Care System, Mather, CA 95655 USA
- Department of Pharmacology, School of Medicine, University of California, Davis, Davis, CA 95616 USA
| |
Collapse
|
11
|
Ushenin K, Kalinin V, Gitinova S, Sopov O, Solovyova O. Parameter variations in personalized electrophysiological models of human heart ventricles. PLoS One 2021; 16:e0249062. [PMID: 33909606 PMCID: PMC8081243 DOI: 10.1371/journal.pone.0249062] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2020] [Accepted: 03/10/2021] [Indexed: 11/18/2022] Open
Abstract
The objectives of this study were to evaluate the accuracy of personalized numerical simulations of the electrical activity in human ventricles by comparing simulated electrocardiograms (ECGs) with real patients’ ECGs and analyzing the sensitivity of the model output to variations in the model parameters. We used standard 12-lead ECGs and up to 224 unipolar body-surface ECGs to record three patients with cardiac resynchronization therapy devices and three patients with focal ventricular tachycardia. Patient-tailored geometrical models of the ventricles, atria, large vessels, liver, and spine were created using computed tomography data. Ten cases of focal ventricular activation were simulated using the bidomain model and the TNNP 2006 cellular model. The population-based values of electrical conductivities and other model parameters were used for accuracy analysis, and their variations were used for sensitivity analysis. The mean correlation coefficient between the simulated and real ECGs varied significantly (from r = 0.29 to r = 0.86) among the simulated cases. A strong mean correlation (r > 0.7) was found in eight of the ten model cases. The accuracy of the ECG simulation varied widely in the same patient depending on the localization of the excitation origin. The sensitivity analysis revealed that variations in the anisotropy ratio, blood conductivity, and cellular apicobasal heterogeneity had the strongest influence on transmembrane potential, while variation in lung conductivity had the greatest influence on body-surface ECGs. Futhermore, the anisotropy ratio predominantly affected the latest activation time and repolarization time dispersion, while the cellular apicobasal heterogeneity mainly affected the dispersion of action potential duration, and variation in lung conductivity mainly led to changes in the amplitudes of ECGs and cardiac electrograms. We also found that the effects of certain parameter variations had specific regional patterns on the cardiac and body surfaces. These observations are useful for further developing personalized cardiac models.
Collapse
Affiliation(s)
- Konstantin Ushenin
- Institute of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg, Russia
- Institute of Immunology and Physiology of the Ural Branch of the RAS, Ekaterinburg, Russia
- * E-mail:
| | | | - Sukaynat Gitinova
- Department of Surgical Treatment of Tachyarrhythmias, A.N. Bakulev National Medical Research Center of Cardiovascular Surgery, Moscow, Russia
| | - Oleg Sopov
- Department of Surgical Treatment of Tachyarrhythmias, A.N. Bakulev National Medical Research Center of Cardiovascular Surgery, Moscow, Russia
| | - Olga Solovyova
- Institute of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg, Russia
- Institute of Immunology and Physiology of the Ural Branch of the RAS, Ekaterinburg, Russia
| |
Collapse
|
12
|
A Framework for the generation of digital twins of cardiac electrophysiology from clinical 12-leads ECGs. Med Image Anal 2021; 71:102080. [PMID: 33975097 DOI: 10.1016/j.media.2021.102080] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 02/15/2021] [Accepted: 04/06/2021] [Indexed: 12/21/2022]
Abstract
Cardiac digital twins (Cardiac Digital Twin (CDT)s) of human electrophysiology (Electrophysiology (EP)) are digital replicas of patient hearts derived from clinical data that match like-for-like all available clinical observations. Due to their inherent predictive potential, CDTs show high promise as a complementary modality aiding in clinical decision making and also in the cost-effective, safe and ethical testing of novel EP device therapies. However, current workflows for both the anatomical and functional twinning phases within CDT generation, referring to the inference of model anatomy and parameters from clinical data, are not sufficiently efficient, robust and accurate for advanced clinical and industrial applications. Our study addresses three primary limitations impeding the routine generation of high-fidelity CDTs by introducing; a comprehensive parameter vector encapsulating all factors relating to the ventricular EP; an abstract reference frame within the model allowing the unattended manipulation of model parameter fields; a novel fast-forward electrocardiogram (Electrocardiogram (ECG)) model for efficient and bio-physically-detailed simulation required for parameter inference. A novel workflow for the generation of CDTs is then introduced as an initial proof of concept. Anatomical twinning was performed within a reasonable time compatible with clinical workflows (<4h) for 12 subjects from clinically-attained magnetic resonance images. After assessment of the underlying fast forward ECG model against a gold standard bidomain ECG model, functional twinning of optimal parameters according to a clinically-attained 12 lead ECG was then performed using a forward Saltelli sampling approach for a single subject. The achieved results in terms of efficiency and fidelity demonstrate that our workflow is well-suited and viable for generating biophysically-detailed CDTs at scale.
Collapse
|
13
|
Pilia N, Severi S, Raimann JG, Genovesi S, Dössel O, Kotanko P, Corsi C, Loewe A. Quantification and classification of potassium and calcium disorders with the electrocardiogram: What do clinical studies, modeling, and reconstruction tell us? APL Bioeng 2020; 4:041501. [PMID: 33062908 PMCID: PMC7532940 DOI: 10.1063/5.0018504] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 09/13/2020] [Indexed: 11/14/2022] Open
Abstract
Diseases caused by alterations of ionic concentrations are frequently observed challenges and play an important role in clinical practice. The clinically established method for the diagnosis of electrolyte concentration imbalance is blood tests. A rapid and non-invasive point-of-care method is yet needed. The electrocardiogram (ECG) could meet this need and becomes an established diagnostic tool allowing home monitoring of the electrolyte concentration also by wearable devices. In this review, we present the current state of potassium and calcium concentration monitoring using the ECG and summarize results from previous work. Selected clinical studies are presented, supporting or questioning the use of the ECG for the monitoring of electrolyte concentration imbalances. Differences in the findings from automatic monitoring studies are discussed, and current studies utilizing machine learning are presented demonstrating the potential of the deep learning approach. Furthermore, we demonstrate the potential of computational modeling approaches to gain insight into the mechanisms of relevant clinical findings and as a tool to obtain synthetic data for methodical improvements in monitoring approaches.
Collapse
Affiliation(s)
- N Pilia
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - S Severi
- Department of Electrical, Electronic, and Information Engineering "Guglielmo Marconi," University of Bologna, 47522 Cesena, Italy
| | - J G Raimann
- Renal Research Institute, New York, New York 10065, USA
| | - S Genovesi
- Department of Medicine and Surgery, University of Milan-Bicocca, 20100 Milan, Italy
| | - O Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | | | - C Corsi
- Department of Electrical, Electronic, and Information Engineering "Guglielmo Marconi," University of Bologna, 47522 Cesena, Italy
| | - A Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
14
|
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
|
15
|
Kalinin A, Potyagaylo D, Kalinin V. Solving the Inverse Problem of Electrocardiography on the Endocardium Using a Single Layer Source. Front Physiol 2019; 10:58. [PMID: 30804802 PMCID: PMC6370732 DOI: 10.3389/fphys.2019.00058] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Accepted: 01/18/2019] [Indexed: 12/04/2022] Open
Abstract
The inverse problem of electrocardiography consists in reconstructing cardiac electrical activity from given body surface electrocardiographic measurements. Despite tremendous progress in the field over the last decades, the solution of this problem in terms of electrical potentials on both epi- and the endocardial heart surfaces with acceptable accuracy remains challenging. This paper presents a novel numerical approach aimed at improving the solution quality on the endocardium. Our method exploits the solution representation in the form of electrical single layer densities on the myocardial surface. We demonstrate that this representation brings twofold benefits: first, the inverse problem can be solved for the physiologically meaningful single layer densities. Secondly, a conventional transfer matrix for electrical potentials can be split into two parts, one of which turned out to posess regularizing properties leading to improved endocardial reconstructions. The method was tested in-silico for ventricular pacings utilizing realistic CT-based heart and torso geometries. The proposed approach provided more accurate solution on the ventricular endocardium compared to the conventional potential-based solutions with Tikhonov regularization of the 0th, 1st, and 2nd orders. Furthermore, we show a uniform spatio-temporal behavior of the single layer densities over the heart surface, which could be conveniently employed in the regularization procedure.
Collapse
|
16
|
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
|
17
|
Loewe A, Wülfers EM, Seemann G. Cardiac ischemia-insights from computational models. Herzschrittmacherther Elektrophysiol 2018; 29:48-56. [PMID: 29305703 DOI: 10.1007/s00399-017-0539-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Accepted: 10/26/2017] [Indexed: 06/07/2023]
Abstract
BACKGROUND Complementary to clinical and experimental studies, computational cardiac modeling serves to obtain a comprehensive understanding of the cardiovascular system in order to analyze dysfunction, evaluate existing, and develop novel treatment strategies. OBJECTIVES We describe the basics of multiscale computational modeling of cardiac electrophysiology from the molecular ion channel to the whole body scale. By modeling cardiac ischemia, we illustrate how in silico experiments can contribute to our understanding of how the pathophysiological mechanisms translate into changes observed in diagnostic tools such as the electrocardiogram (ECG). MATERIALS AND METHODS Quantitative in silico modeling spans a wide range of scales from ion channel biophysics to ECG signals. For each of the scales, a set of mathematical equations describes electrophysiology in relation to the other scales. Integration of ischemia-induced changes is performed on the ion channel, single-cell, and tissue level. This approach allows us to study how effects simulated at molecular scales translate to changes in the ECG. RESULTS Ischemia induces action potential shortening and conduction slowing. Hence, ischemic myocardium has distinct and significant effects on propagation and repolarization of excitation, depending on the intramural extent of the ischemic region. For transmural and subendocardial ischemic regions, ST segment elevation and depression, respectively, were observed, whereas intermediate ischemic regions were found to be electrically silent (NSTEMI). CONCLUSIONS In silico modeling contributes quantitative and mechanistic insight into fundamental ischemia-related arrhythmogenic mechanisms. In addition, computational modeling can help to translate experimental findings at the (sub-)cellular level to the organ and body context (e. g., ECG), thereby providing a thorough understanding of this routinely used diagnostic tool that may translate into optimized applications.
Collapse
Affiliation(s)
- Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Eike Moritz Wülfers
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg, Bad Krozingen, Medical Center, Computational Modeling Group, Albert-Ludwigs University of Freiburg, Elsässerstr. 2q, 79110, Freiburg, Germany
- Faculty of Medicine, University of Freiburg, Freiburg, 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, Medical Center, Computational Modeling Group, Albert-Ludwigs University of Freiburg, Elsässerstr. 2q, 79110, Freiburg, Germany.
- Faculty of Medicine, University of Freiburg, Freiburg, Germany.
| |
Collapse
|
18
|
Bayer J, Prassl AJ, Pashaei A, Gomez JF, Frontera A, Neic A, Plank G, Vigmond EJ. Universal ventricular coordinates: A generic framework for describing position within the heart and transferring data. Med Image Anal 2018; 45:83-93. [PMID: 29414438 DOI: 10.1016/j.media.2018.01.005] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 01/16/2018] [Accepted: 01/22/2018] [Indexed: 12/20/2022]
Abstract
Being able to map a particular set of cardiac ventricles to a generic topologically equivalent representation has many applications, including facilitating comparison of different hearts, as well as mapping quantities and structures of interest between them. In this paper we describe Universal Ventricular Coordinates (UVC), which can be used to describe position within any biventricular heart. UVC comprise four unique coordinates that we have chosen to be intuitive, well defined, and relevant for physiological descriptions. We describe how to determine these coordinates for any volumetric mesh by illustrating how to properly assign boundary conditions and utilize solutions to Laplace's equation. Using UVC, we transferred scalar, vector, and tensor data between four unstructured ventricular meshes from three different species. Performing the mappings was very fast, on the order of a few minutes, since mesh nodes were searched in a KD tree. Distance errors in mapping mesh nodes back and forth between meshes were less than the size of an element. Analytically derived fiber directions were also mapped across meshes and compared, showing < 5° difference over most of the ventricles. The ability to transfer gradients was also demonstrated. Topologically variable structures, like papillary muscles, required further definition outside of the UVC framework. In conclusion, UVC can aid in transferring many types of data between different biventricular geometries.
Collapse
Affiliation(s)
- Jason Bayer
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, avenue du Haut-Lévèque, Pessac 33600, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France.
| | - Anton J Prassl
- Gottfried Schatz Research Center, Biophysics, Medical University of Graz, Neue Stiftingtalstrasse 6, 8010 Graz, Austria.
| | - Ali Pashaei
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, avenue du Haut-Lévèque, Pessac 33600, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France.
| | - Juan F Gomez
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, avenue du Haut-Lévèque, Pessac 33600, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France.
| | - Antonio Frontera
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, avenue du Haut-Lévèque, Pessac 33600, France; Department of Electrophysiology, Hôpital Haut Lévèque, 1 avenue Magellan, Pessac 33100 France.
| | - Aurel Neic
- Gottfried Schatz Research Center, Biophysics, Medical University of Graz, Neue Stiftingtalstrasse 6, 8010 Graz, Austria.
| | - Gernot Plank
- Gottfried Schatz Research Center, Biophysics, Medical University of Graz, Neue Stiftingtalstrasse 6, 8010 Graz, Austria.
| | - Edward J Vigmond
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, avenue du Haut-Lévèque, Pessac 33600, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France.
| |
Collapse
|
19
|
Van Nieuwenhuyse E, Seemann G, Panfilov AV, Vandersickel N. Effects of early afterdepolarizations on excitation patterns in an accurate model of the human ventricles. PLoS One 2017; 12:e0188867. [PMID: 29216239 PMCID: PMC5720514 DOI: 10.1371/journal.pone.0188867] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Accepted: 11/14/2017] [Indexed: 12/17/2022] Open
Abstract
Early Afterdepolarizations, EADs, are defined as the reversal of the action potential before completion of the repolarization phase, which can result in ectopic beats. However, the series of mechanisms of EADs leading to these ectopic beats and related cardiac arrhythmias are not well understood. Therefore, we aimed to investigate the influence of this single cell behavior on the whole heart level. For this study we used a modified version of the Ten Tusscher-Panfilov model of human ventricular cells (TP06) which we implemented in a 3D ventricle model including realistic fiber orientations. To increase the likelihood of EAD formation at the single cell level, we reduced the repolarization reserve (RR) by reducing the rapid delayed rectifier Potassium current and raising the L-type Calcium current. Varying these parameters defined a 2D parametric space where different excitation patterns could be classified. Depending on the initial conditions, by either exciting the ventricles with a spiral formation or burst pacing protocol, we found multiple different spatio-temporal excitation patterns. The spiral formation protocol resulted in the categorization of a stable spiral (S), a meandering spiral (MS), a spiral break-up regime (SB), spiral fibrillation type B (B), spiral fibrillation type A (A) and an oscillatory excitation type (O). The last three patterns are a 3D generalization of previously found patterns in 2D. First, the spiral fibrillation type B showed waves determined by a chaotic bi-excitable regime, i.e. mediated by both Sodium and Calcium waves at the same time and in same tissue settings. In the parameter region governed by the B pattern, single cells were able to repolarize completely and different (spiral) waves chaotically burst into each other without finishing a 360 degree rotation. Second, spiral fibrillation type A patterns consisted of multiple small rotating spirals. Single cells failed to repolarize to the resting membrane potential hence prohibiting the Sodium channel gates to recover. Accordingly, we found that Calcium waves mediated these patterns. Third, a further reduction of the RR resulted in a more exotic parameter regime whereby the individual cells behaved independently as oscillators. The patterns arose due to a phase-shift of different oscillators as disconnection of the cells resulted in continuation of the patterns. For all patterns, we computed realistic 9 lead ECGs by including a torso model. The B and A type pattern exposed the behavior of Ventricular Tachycardia (VT). We conclude that EADs at the single cell level can result in different types of cardiac fibrillation at the tissue and 3D ventricle level.
Collapse
Affiliation(s)
| | - Gunnar Seemann
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg, Bad Krozingen, Germany
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | | | - Nele Vandersickel
- Department of Physics and Astronomy, Ghent University, Ghent, Belgium
| |
Collapse
|
20
|
Neic A, Campos FO, Prassl AJ, Niederer SA, Bishop MJ, Vigmond EJ, Plank G. Efficient computation of electrograms and ECGs in human whole heart simulations using a reaction-eikonal model. JOURNAL OF COMPUTATIONAL PHYSICS 2017; 346:191-211. [PMID: 28819329 PMCID: PMC5555399 DOI: 10.1016/j.jcp.2017.06.020] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Anatomically accurate and biophysically detailed bidomain models of the human heart have proven a powerful tool for gaining quantitative insight into the links between electrical sources in the myocardium and the concomitant current flow in the surrounding medium as they represent their relationship mechanistically based on first principles. Such models are increasingly considered as a clinical research tool with the perspective of being used, ultimately, as a complementary diagnostic modality. An important prerequisite in many clinical modeling applications is the ability of models to faithfully replicate potential maps and electrograms recorded from a given patient. However, while the personalization of electrophysiology models based on the gold standard bidomain formulation is in principle feasible, the associated computational expenses are significant, rendering their use incompatible with clinical time frames. In this study we report on the development of a novel computationally efficient reaction-eikonal (R-E) model for modeling extracellular potential maps and electrograms. Using a biventricular human electrophysiology model, which incorporates a topologically realistic His-Purkinje system (HPS), we demonstrate by comparing against a high-resolution reaction-diffusion (R-D) bidomain model that the R-E model predicts extracellular potential fields, electrograms as well as ECGs at the body surface with high fidelity and offers vast computational savings greater than three orders of magnitude. Due to their efficiency R-E models are ideally suitable for forward simulations in clinical modeling studies which attempt to personalize electrophysiological model features.
Collapse
Affiliation(s)
- Aurel Neic
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Fernando O. Campos
- Institute of Biophysics, Medical University of Graz, Graz, Austria
- Dept. of Congenital Heart Diseases and Pediatric Cardiology, German Heart Institute Berlin, Berlin, Germany
| | - Anton J. Prassl
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Steven A. Niederer
- Dept. Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College of London, London, United Kingdom
| | - Martin J. Bishop
- Dept. Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College of London, London, United Kingdom
| | | | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Austria
- Corresponding author. (G. Plank)
| |
Collapse
|
21
|
Walmsley J, van Everdingen W, Cramer MJ, Prinzen FW, Delhaas T, Lumens J. Combining computer modelling and cardiac imaging to understand right ventricular pump function. Cardiovasc Res 2017; 113:1486-1498. [DOI: 10.1093/cvr/cvx154] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/01/2017] [Accepted: 08/08/2017] [Indexed: 11/13/2022] Open
|
22
|
Augustin CM, Crozier A, Neic A, Prassl AJ, Karabelas E, Ferreira da Silva T, Fernandes JF, Campos F, Kuehne T, Plank G. Patient-specific modeling of left ventricular electromechanics as a driver for haemodynamic analysis. Europace 2017; 18:iv121-iv129. [PMID: 28011839 PMCID: PMC5386137 DOI: 10.1093/europace/euw369] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2016] [Accepted: 08/26/2016] [Indexed: 01/30/2023] Open
Abstract
Aims Models of blood flow in the left ventricle (LV) and aorta are an important tool for analysing the interplay between LV deformation and flow patterns. Typically, image-based kinematic models describing endocardial motion are used as an input to blood flow simulations. While such models are suitable for analysing the hemodynamic status quo, they are limited in predicting the response to interventions that alter afterload conditions. Mechano-fluidic models using biophysically detailed electromechanical (EM) models have the potential to overcome this limitation, but are more costly to build and compute. We report our recent advancements in developing an automated workflow for the creation of such CFD ready kinematic models to serve as drivers of blood flow simulations. Methods and results EM models of the LV and aortic root were created for four pediatric patients treated for either aortic coarctation or aortic valve disease. Using MRI, ECG and invasive pressure recordings, anatomy as well as electrophysiological, mechanical and circulatory model components were personalized. Results The implemented modeling pipeline was highly automated and allowed model construction and execution of simulations of a patient’s heartbeat within 1 day. All models reproduced clinical data with acceptable accuracy. Conclusion Using the developed modeling workflow, the use of EM LV models as driver of fluid flow simulations is becoming feasible. While EM models are costly to construct, they constitute an important and nontrivial step towards fully coupled electro-mechano-fluidic (EMF) models and show promise as a tool for predicting the response to interventions which affect afterload conditions.
Collapse
Affiliation(s)
- Christoph M Augustin
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria.,Department of Mechanical Engineering, University of California, 5126 Etcheverry Hall, Berkeley, CA 94720, USA
| | - Andrew Crozier
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria
| | - Aurel Neic
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria
| | - Anton J Prassl
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria
| | - Elias Karabelas
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria
| | - Tiago Ferreira da Silva
- Department of Congenital Heart Disease/Pediatric Cardiology, German Heart Institute Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Joao F Fernandes
- Department of Congenital Heart Disease/Pediatric Cardiology, German Heart Institute Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Fernando Campos
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria.,Department of Congenital Heart Disease/Pediatric Cardiology, German Heart Institute Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Titus Kuehne
- Department of Congenital Heart Disease/Pediatric Cardiology, German Heart Institute Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010 Graz, Austria
| |
Collapse
|
23
|
Heine T, Lenis G, Reichensperger P, Beran T, Doessel O, Deml B. Electrocardiographic features for the measurement of drivers' mental workload. APPLIED ERGONOMICS 2017; 61:31-43. [PMID: 28237018 DOI: 10.1016/j.apergo.2016.12.015] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Revised: 11/26/2016] [Accepted: 12/21/2016] [Indexed: 06/06/2023]
Abstract
This study examines the effect of mental workload on the electrocardiogram (ECG) of participants driving the Lane Change Task (LCT). Different levels of mental workload were induced by a secondary task (n-back task) with three levels of difficulty. Subjective data showed a significant increase of the experienced workload over all three levels. An exploratory approach was chosen to extract a large number of rhythmical and morphological features from the ECG signal thereby identifying those which differentiated best between the levels of mental workload. No single rhythmical or morphological feature was able to differentiate between all three levels. A group of parameters were extracted which were at least able to discriminate between two levels. For future research, a combination of features is recommended to achieve best diagnosticity for different levels of mental workload.
Collapse
Affiliation(s)
- Tobias Heine
- Institute of Human and Industrial Engineering (ifab), Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - Gustavo Lenis
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Patrick Reichensperger
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Tobias Beran
- Institute of Human and Industrial Engineering (ifab), Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Olaf Doessel
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Barbara Deml
- Institute of Human and Industrial Engineering (ifab), Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
24
|
Lenis G, Pilia N, Loewe A, Schulze WHW, Dössel O. Comparison of Baseline Wander Removal Techniques considering the Preservation of ST Changes in the Ischemic ECG: A Simulation Study. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2017; 2017:9295029. [PMID: 28373893 PMCID: PMC5361052 DOI: 10.1155/2017/9295029] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/12/2016] [Revised: 02/10/2017] [Accepted: 02/19/2017] [Indexed: 11/25/2022]
Abstract
The most important ECG marker for the diagnosis of ischemia or infarction is a change in the ST segment. Baseline wander is a typical artifact that corrupts the recorded ECG and can hinder the correct diagnosis of such diseases. For the purpose of finding the best suited filter for the removal of baseline wander, the ground truth about the ST change prior to the corrupting artifact and the subsequent filtering process is needed. In order to create the desired reference, we used a large simulation study that allowed us to represent the ischemic heart at a multiscale level from the cardiac myocyte to the surface ECG. We also created a realistic model of baseline wander to evaluate five filtering techniques commonly used in literature. In the simulation study, we included a total of 5.5 million signals coming from 765 electrophysiological setups. We found that the best performing method was the wavelet-based baseline cancellation. However, for medical applications, the Butterworth high-pass filter is the better choice because it is computationally cheap and almost as accurate. Even though all methods modify the ST segment up to some extent, they were all proved to be better than leaving baseline wander unfiltered.
Collapse
Affiliation(s)
- Gustavo Lenis
- Karlsruhe Institute of Technology (KIT), Institute of Biomedical Engineering (IBT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Nicolas Pilia
- Karlsruhe Institute of Technology (KIT), Institute of Biomedical Engineering (IBT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Axel Loewe
- Karlsruhe Institute of Technology (KIT), Institute of Biomedical Engineering (IBT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | | | - Olaf Dössel
- Karlsruhe Institute of Technology (KIT), Institute of Biomedical Engineering (IBT), Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| |
Collapse
|
25
|
Lines GT, de Oliveira BL, Skavhaug O, Maleckar MM. Simple T-Wave Metrics May Better Predict Early Ischemia as Compared to ST Segment. IEEE Trans Biomed Eng 2016; 64:1305-1309. [PMID: 27576235 DOI: 10.1109/tbme.2016.2600198] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
There is pressing clinical need to identify developing heart attack (infarction) in patients as early as possible. However, current state-of-the-art tools in clinical practice, underpinned by the evaluation of elevation of the ST segment of the 12-lead electrocardiogram (ECG), do not identify all patients suffering from lack of blood flow to the heart muscle (cardiac ischemia), worsening the risk for further adverse events and patient outcome overall. In this study, we aimed to explore and compare the portions of cardiac repolarization in the ECG that best capture the electrophysiological changes associated with ischemia. We developed three-dimensional electrophysiological models of the human ventricles and torso, incorporating biophysically-based membrane kinetics and realistic activation sequence, to compute simulated ECGs and their alteration with the application of simulated ischemia of differing severity in diverse regions of the heart. Results suggest that metrics based on the T-wave in addition to the ST segment may be more sensitive to detecting ischemia than those using the ST segment alone. Further research into how such simulation-aided risk assessment methods may aid workflows in extant clinical practice, with the ultimate goal of multimodality clinical support, is warranted.
Collapse
|
26
|
Alday EAP, Ni H, Zhang C, Colman MA, Gan Z, Zhang H. Comparison of Electric- and Magnetic-Cardiograms Produced by Myocardial Ischemia in Models of the Human Ventricle and Torso. PLoS One 2016; 11:e0160999. [PMID: 27556808 PMCID: PMC4996509 DOI: 10.1371/journal.pone.0160999] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Accepted: 07/28/2016] [Indexed: 11/18/2022] Open
Abstract
Myocardial ventricular ischemia arises from a lack of blood supply to the heart, which may cause abnormal repolarization and excitation wave conduction patterns in the tissue, leading to cardiac arrhythmias and even sudden death. Current diagnosis of cardiac ischemia by the 12-lead electrocardiogram (ECG) has limitations as they are insensitive in many cases and may show unnoticeable differences to normal patterns. As the magnetic field provides extra information on cardiac excitation and is more sensitive to tangential currents to the surface of the chest, whereas the electric field is more sensitive to flux currents, it has been hypothesized that the magnetocardiogram (MCG) may provide a complementary method to the ECG in ischemic diagnosis. However, it is unclear yet about the differences in sensitivity regions of body surface ECG and MCG signals to ischemic conditions. The aim of this study was to investigate such differences by using 12-, 36- ECG and 36-MCG computed from multi-scale biophysically detailed computational models of the human ventricles and torso in both control and ischemic conditions. It was shown that ischemia produced changes in the ECG and MCG signals in the QRS complex, T-wave and ST-segment, with greater relative differences seen in the 36-lead ECG and MCG as compared to the 12-leads ECG (34% and 37% vs 26%, respectively). The 36-lead ECG showed more averaged sensitivity than the MCG in the change of T-wave due to ischemia (37% vs 32%, respectively), whereas the MCG showed greater sensitivity than the ECG in the change of the ST-segment (50% vs 40%, respectively). In addition, both MCG and ECG showed regional-dependent changes to ischemia, but with MCG showing a stronger correlation between ischemic region in the heart. In conclusion, MCG shows more sensitivity than ECG in response to ischemia, which may provide an alternative method for the diagnosis of ischemia.
Collapse
Affiliation(s)
- Erick A. Perez Alday
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
| | - Haibo Ni
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
| | - Chen Zhang
- Applied superconductivity Research Center, School of Physics, Peking University, Beijing, China
| | - Michael A. Colman
- Theoretical Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
| | - Zizhao Gan
- Applied superconductivity Research Center, School of Physics, Peking University, Beijing, China
| | - Henggui Zhang
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
- * E-mail:
| |
Collapse
|
27
|
Colli Franzone P, Pavarino LF, Scacchi S. Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study. Math Biosci 2016; 280:71-86. [PMID: 27545966 DOI: 10.1016/j.mbs.2016.08.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Revised: 05/25/2016] [Accepted: 08/04/2016] [Indexed: 10/21/2022]
Abstract
The aim of this work is to investigate, by means of numerical simulations, the influence of myocardial deformation due to muscle contraction and relaxation on the cardiac repolarization process in presence of transmural intrinsic action potential duration (APD) heterogeneities. The three-dimensional electromechanical model considered consists of the following four coupled components: the quasi-static transversely isotropic finite elasticity equations for the deformation of the cardiac tissue; the active tension model for the intracellular calcium dynamics and cross-bridge binding; the anisotropic Bidomain model for the electrical current flow through the deforming cardiac tissue; the membrane model of ventricular myocytes, including stretch-activated channels. The numerical simulations are based on our finite element parallel solver, which employs Multilevel Additive Schwarz preconditioners for the solution of the discretized Bidomain equations and Newton-Krylov methods for the solution of the discretized non-linear finite elasticity equations. Our findings show that: (i) the presence of intrinsic transmural cellular APD heterogeneities is not fully masked by electrotonic current flow or by the presence of the mechanical deformation; (ii) despite the presence of transmural APD heterogeneities, the recovery process follows the activation sequence and there is no significant transmural repolarization gradient; (iii) with or without transmural APD heterogeneities, epicardial electrograms always display the same wave shape and discordance between the polarity of QRS complex and T-wave; (iv) the main effects of the mechanical deformation are an increase of the dispersion of repolarization time and APD, when computed over the total cardiac domain and over the endo- and epicardial surfaces, while there is a slight decrease along the transmural direction.
Collapse
Affiliation(s)
- P Colli Franzone
- Dipartimento di Matematica, Università di Pavia, Via Ferrata 1, Pavia 27100, Italy.
| | - L F Pavarino
- Dipartimento di Matematica, Università di Milano, Via Saldini 50, Milano 20133, Italy.
| | - S Scacchi
- Dipartimento di Matematica, Università di Milano, Via Saldini 50, Milano 20133, Italy.
| |
Collapse
|
28
|
Schenone E, Collin A, Gerbeau JF. Numerical simulation of electrocardiograms for full cardiac cycles in healthy and pathological conditions. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2016; 32:e02744. [PMID: 26249327 DOI: 10.1002/cnm.2744] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2014] [Revised: 07/29/2015] [Accepted: 08/03/2015] [Indexed: 06/04/2023]
Abstract
This work is dedicated to the simulation of full cycles of the electrical activity of the heart and the corresponding body surface potential. The model is based on a realistic torso and heart anatomy, including ventricles and atria. One of the specificities of our approach is to model the atria as a surface, which is the kind of data typically provided by medical imaging for thin volumes. The bidomain equations are considered in their usual formulation in the ventricles, and in a surface formulation on the atria. Two ionic models are used: the Courtemanche-Ramirez-Nattel model on the atria and the 'minimal model for human ventricular action potentials' by Bueno-Orovio, Cherry, and Fenton in the ventricles. The heart is weakly coupled to the torso by a Robin boundary condition based on a resistor-capacitor transmission condition. Various electrocardiograms (ECGs) are simulated in healthy and pathological conditions (left and right bundle branch blocks, Bachmann's bundle block, and Wolff-Parkinson-White syndrome). To assess the numerical ECGs, we use several qualitative and quantitative criteria found in the medical literature. Our simulator can also be used to generate the signals measured by a vest of electrodes. This capability is illustrated at the end of the article. Copyright © 2015 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Elisa Schenone
- Sorbonne Universités UPMC, Paris, France
- Inria Paris-Rocquencourt, Paris, France
| | | | | |
Collapse
|
29
|
Transmural, interventricular, apicobasal and anteroposterior action potential duration gradients are all essential to the genesis of the concordant and realistic T wave: A whole-heart model study. J Electrocardiol 2016; 49:569-78. [PMID: 27034121 DOI: 10.1016/j.jelectrocard.2016.03.010] [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: 10/12/2015] [Indexed: 12/12/2022]
Abstract
BACKGROUND It has been reported that ventricular repolarization dispersion resulting from transmural, apicobasal and interventricular action potential duration (APD) gradients makes the T wave concordant with the QRS complex. METHOD AND RESULTS A whole-heart model integrating transmural, apicobasal, interventricular and anteroposterior APD gradients was used, and the corresponding electrocardiograms were simulated to study the influence of these APD gradients on the T-wave amplitudes. The simulation results showed that changing a single APD gradient (e.g., interventricular APD gradient alone) only made substantial changes to the T-wave amplitudes in a limited number of leads and was not able to generate T waves with amplitudes comparable with clinical findings in all leads. A combination of transmural, apicobasal and interventricular APD gradients could simulate T waves with amplitudes similar to clinical values in the limb leads only. Adding the anteroposterior APD gradient into the model greatly improved the consistency between the simulated T-wave amplitudes and the clinical values. CONCLUSION The simulation results support that the transmural, apicobasal, interventricular and the anteroposterior APD gradient are all essential to the genesis of the clinical T wave.
Collapse
|
30
|
Augustin CM, Neic A, Liebmann M, Prassl AJ, Niederer SA, Haase G, Plank G. Anatomically accurate high resolution modeling of human whole heart electromechanics: A strongly scalable algebraic multigrid solver method for nonlinear deformation. JOURNAL OF COMPUTATIONAL PHYSICS 2016; 305:622-646. [PMID: 26819483 PMCID: PMC4724941 DOI: 10.1016/j.jcp.2015.10.045] [Citation(s) in RCA: 76] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Electromechanical (EM) models of the heart have been used successfully to study fundamental mechanisms underlying a heart beat in health and disease. However, in all modeling studies reported so far numerous simplifications were made in terms of representing biophysical details of cellular function and its heterogeneity, gross anatomy and tissue microstructure, as well as the bidirectional coupling between electrophysiology (EP) and tissue distension. One limiting factor is the employed spatial discretization methods which are not sufficiently flexible to accommodate complex geometries or resolve heterogeneities, but, even more importantly, the limited efficiency of the prevailing solver techniques which are not sufficiently scalable to deal with the incurring increase in degrees of freedom (DOF) when modeling cardiac electromechanics at high spatio-temporal resolution. This study reports on the development of a novel methodology for solving the nonlinear equation of finite elasticity using human whole organ models of cardiac electromechanics, discretized at a high para-cellular resolution. Three patient-specific, anatomically accurate, whole heart EM models were reconstructed from magnetic resonance (MR) scans at resolutions of 220 μm, 440 μm and 880 μm, yielding meshes of approximately 184.6, 24.4 and 3.7 million tetrahedral elements and 95.9, 13.2 and 2.1 million displacement DOF, respectively. The same mesh was used for discretizing the governing equations of both electrophysiology (EP) and nonlinear elasticity. A novel algebraic multigrid (AMG) preconditioner for an iterative Krylov solver was developed to deal with the resulting computational load. The AMG preconditioner was designed under the primary objective of achieving favorable strong scaling characteristics for both setup and solution runtimes, as this is key for exploiting current high performance computing hardware. Benchmark results using the 220 μm, 440 μm and 880 μm meshes demonstrate efficient scaling up to 1024, 4096 and 8192 compute cores which allowed the simulation of a single heart beat in 44.3, 87.8 and 235.3 minutes, respectively. The efficiency of the method allows fast simulation cycles without compromising anatomical or biophysical detail.
Collapse
Affiliation(s)
| | - Aurel Neic
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Manfred Liebmann
- Institute for Mathematics and Scientific Computing, Karl-Franzens-University Graz, Graz, Austria
| | - Anton J. Prassl
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Steven A. Niederer
- Dept. Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College of London, London, United Kingdom
| | - Gundolf Haase
- Institute for Mathematics and Scientific Computing, Karl-Franzens-University Graz, Graz, Austria
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Austria
- Corresponding author (Gernot Plank)
| |
Collapse
|
31
|
Crozier A, Augustin CM, Neic A, Prassl AJ, Holler M, Fastl TE, Hennemuth A, Bredies K, Kuehne T, Bishop MJ, Niederer SA, Plank G. Image-Based Personalization of Cardiac Anatomy for Coupled Electromechanical Modeling. Ann Biomed Eng 2016. [PMID: 26424476 DOI: 10.1007/sl0439-015-1474-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Computational models of cardiac electromechanics (EM) are increasingly being applied to clinical problems, with patient-specific models being generated from high fidelity imaging and used to simulate patient physiology, pathophysiology and response to treatment. Current structured meshes are limited in their ability to fully represent the detailed anatomical data available from clinical images and capture complex and varied anatomy with limited geometric accuracy. In this paper, we review the state of the art in image-based personalization of cardiac anatomy for biophysically detailed, strongly coupled EM modeling, and present our own tools for the automatic building of anatomically and structurally accurate patient-specific models. Our method relies on using high resolution unstructured meshes for discretizing both physics, electrophysiology and mechanics, in combination with efficient, strongly scalable solvers necessary to deal with the computational load imposed by the large number of degrees of freedom of these meshes. These tools permit automated anatomical model generation and strongly coupled EM simulations at an unprecedented level of anatomical and biophysical detail.
Collapse
Affiliation(s)
- A Crozier
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - C M Augustin
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - A Neic
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - A J Prassl
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - M Holler
- Institute for Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - T E Fastl
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - A Hennemuth
- Modeling and Simulation Group, Fraunhofer MEVIS, Bremen, Germany
| | - K Bredies
- Institute for Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - T Kuehne
- Non-Invasive Cardiac Imaging in Congenital Heart Disease Unit, Charité-Universitätsmedizin, Berlin, Germany
- German Heart Institute, Berlin, Germany
| | - M J Bishop
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - S A Niederer
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - G Plank
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria.
| |
Collapse
|
32
|
Johnstone RH, Chang ETY, Bardenet R, de Boer TP, Gavaghan DJ, Pathmanathan P, Clayton RH, Mirams GR. Uncertainty and variability in models of the cardiac action potential: Can we build trustworthy models? J Mol Cell Cardiol 2015; 96:49-62. [PMID: 26611884 PMCID: PMC4915860 DOI: 10.1016/j.yjmcc.2015.11.018] [Citation(s) in RCA: 84] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/31/2015] [Revised: 10/13/2015] [Accepted: 11/17/2015] [Indexed: 01/07/2023]
Abstract
Cardiac electrophysiology models have been developed for over 50 years, and now include detailed descriptions of individual ion currents and sub-cellular calcium handling. It is commonly accepted that there are many uncertainties in these systems, with quantities such as ion channel kinetics or expression levels being difficult to measure or variable between samples. Until recently, the original approach of describing model parameters using single values has been retained, and consequently the majority of mathematical models in use today provide point predictions, with no associated uncertainty. In recent years, statistical techniques have been developed and applied in many scientific areas to capture uncertainties in the quantities that determine model behaviour, and to provide a distribution of predictions which accounts for this uncertainty. In this paper we discuss this concept, which is termed uncertainty quantification, and consider how it might be applied to cardiac electrophysiology models. We present two case studies in which probability distributions, instead of individual numbers, are inferred from data to describe quantities such as maximal current densities. Then we show how these probabilistic representations of model parameters enable probabilities to be placed on predicted behaviours. We demonstrate how changes in these probability distributions across data sets offer insight into which currents cause beat-to-beat variability in canine APs. We conclude with a discussion of the challenges that this approach entails, and how it provides opportunities to improve our understanding of electrophysiology. Uncertainty and variability in action potential models can be quantified. A probabilistic method for inferring maximal current densities is developed and applied. We use this to infer the currents responsible for canine beat-to-beat variability. Emulation of mathematical models provides rich information at low computational cost. The importance of considering uncertainty and variability in future is discussed.
Collapse
Affiliation(s)
- Ross H Johnstone
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - Eugene T Y Chang
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield S1 4DP, UK
| | - Rémi Bardenet
- CNRS & CRIStAL, Université de Lille, 59651 Villeneuve d'Ascq, France
| | - Teun P de Boer
- Division of Heart & Lungs, Department of Medical Physiology, University Medical Center Utrecht, Utrecht, The Netherlands
| | - David J Gavaghan
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - Pras Pathmanathan
- U.S. Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20993, USA.
| | - Richard H Clayton
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield S1 4DP, UK.
| | - Gary R Mirams
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK.
| |
Collapse
|
33
|
Crozier A, Augustin CM, Neic A, Prassl AJ, Holler M, Fastl TE, Hennemuth A, Bredies K, Kuehne T, Bishop MJ, Niederer SA, Plank G. Image-Based Personalization of Cardiac Anatomy for Coupled Electromechanical Modeling. Ann Biomed Eng 2015; 44:58-70. [PMID: 26424476 PMCID: PMC4690840 DOI: 10.1007/s10439-015-1474-5] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Accepted: 09/24/2015] [Indexed: 11/26/2022]
Abstract
Computational models of cardiac electromechanics (EM) are increasingly being applied to clinical problems, with patient-specific models being generated from high fidelity imaging and used to simulate patient physiology, pathophysiology and response to treatment. Current structured meshes are limited in their ability to fully represent the detailed anatomical data available from clinical images and capture complex and varied anatomy with limited geometric accuracy. In this paper, we review the state of the art in image-based personalization of cardiac anatomy for biophysically detailed, strongly coupled EM modeling, and present our own tools for the automatic building of anatomically and structurally accurate patient-specific models. Our method relies on using high resolution unstructured meshes for discretizing both physics, electrophysiology and mechanics, in combination with efficient, strongly scalable solvers necessary to deal with the computational load imposed by the large number of degrees of freedom of these meshes. These tools permit automated anatomical model generation and strongly coupled EM simulations at an unprecedented level of anatomical and biophysical detail.
Collapse
Affiliation(s)
- A Crozier
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - C M Augustin
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - A Neic
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - A J Prassl
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria
| | - M Holler
- Institute for Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - T E Fastl
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - A Hennemuth
- Modeling and Simulation Group, Fraunhofer MEVIS, Bremen, Germany
| | - K Bredies
- Institute for Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - T Kuehne
- Non-Invasive Cardiac Imaging in Congenital Heart Disease Unit, Charité-Universitätsmedizin, Berlin, Germany
- German Heart Institute, Berlin, Germany
| | - M J Bishop
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - S A Niederer
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - G Plank
- Institute of Biophysics, Medical University of Graz, Harrachgasse 21/IV, 8010, Graz, Austria.
| |
Collapse
|
34
|
Perotti LE, Krishnamoorthi S, Borgstrom NP, Ennis DB, Klug WS. Regional segmentation of ventricular models to achieve repolarization dispersion in cardiac electrophysiology modeling. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2015; 31:10.1002/cnm.2718. [PMID: 25845576 PMCID: PMC4519348 DOI: 10.1002/cnm.2718] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2014] [Revised: 03/01/2015] [Accepted: 03/31/2015] [Indexed: 05/08/2023]
Abstract
The electrocardiogram (ECG) is one of the most significant outputs of a computational model of cardiac electrophysiology because it relates the numerical results to clinical data and is a universal tool for diagnosing heart diseases. One key features of the ECG is the T-wave, which is caused by longitudinal and transmural heterogeneity of the action potential duration (APD). Thus, in order to model a correct wave of repolarization, different cell properties resulting in different APDs must be assigned across the ventricular wall and longitudinally from apex to base. To achieve this requirement, a regional parametrization of the heart is necessary. We propose a robust approach to obtain the transmural and longitudinal segmentation in a general heart geometry without relying on ad hoc procedures. Our approach is based on auxiliary harmonic lifting analyses, already used in the literature to generate myocardial fiber orientations. Specifically, the solution of a sequence of Laplace boundary value problems allows parametrically controlled segmentation of both heart ventricles. The flexibility and simplicity of the proposed method is demonstrated through several representative examples, varying the locations and extents of the epicardial, midwall, and endocardial layers. Effects of the control parameters on the T-wave morphology are illustrated via computed ECGs.
Collapse
Affiliation(s)
- L. E. Perotti
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, California, United States of America
- Department of Bioengineering, University of California, Los Angeles, California, United States of America
- Department of Radiological Sciences, University of California, Los Angeles, California, United States of America
| | - S. Krishnamoorthi
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, California, United States of America
| | - N. P. Borgstrom
- Department of Bioengineering, University of California, Los Angeles, California, United States of America
| | - D. B. Ennis
- Department of Bioengineering, University of California, Los Angeles, California, United States of America
- Department of Radiological Sciences, University of California, Los Angeles, California, United States of America
| | - W. S. Klug
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, California, United States of America
- Correspondence to: Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, 420 Westwood Plaza, Los Angeles, CA 90095, United States of America.
| |
Collapse
|
35
|
Arteyeva NV, Azarov JE, Vityazev VA, Shmakov DN. Action potential duration gradients in the heart ventricles and the cardiac electric field during ventricular repolarization (a model study). J Electrocardiol 2015; 48:678-85. [PMID: 25818745 DOI: 10.1016/j.jelectrocard.2015.03.010] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2014] [Indexed: 12/14/2022]
Abstract
BACKGROUND We simulated contributions of transmural, apicobasal, anteroposterior and interventricular action potential duration (APD) gradients to the body surface potential distribution (BSPD) with constant or varied magnitudes of the transmural and apicobasal gradients. METHODS Simulations were done in the framework of the discrete computer model of the rabbit heart ventricles on the basis of realistic activation sequence and APDs. The APD gradients were set constant at 20 ms or varied in the range of ±80 ms. RESULTS The apicobasal, transmural and interventricular APD gradients of 20 ms produced similar BSPDs, whereas the BSPD inversion was caused by the inverted apicobasal or transmural 80 ms gradients. The transmural APD gradient produced transversal and mainly apicobasal T-wave vectors due to wall curvature and cancellation effects. The "normal" transversal and apicobasal repolarization gradients were decreased and increased by activation sequence, respectively. CONCLUSION The different APD gradients contributed consistently to the development of BSPD.
Collapse
Affiliation(s)
- Natalia V Arteyeva
- Laboratory of Cardiac Physiology, Institute of Physiology, Komi Science Center, Ural Branch, Russian Academy of Sciences, 50, Pervomayskaya St., Syktyvkar, Russia
| | - Jan E Azarov
- Laboratory of Cardiac Physiology, Institute of Physiology, Komi Science Center, Ural Branch, Russian Academy of Sciences, 50, Pervomayskaya St., Syktyvkar, Russia; Department of Physiology, Medical Institute of Syktyvkar State University, 11, Babushkin St., Syktyvkar, Russia.
| | - Vladimir A Vityazev
- Laboratory of Cardiac Physiology, Institute of Physiology, Komi Science Center, Ural Branch, Russian Academy of Sciences, 50, Pervomayskaya St., Syktyvkar, Russia; Department of Physiology, Medical Institute of Syktyvkar State University, 11, Babushkin St., Syktyvkar, Russia
| | - Dmitry N Shmakov
- Laboratory of Cardiac Physiology, Institute of Physiology, Komi Science Center, Ural Branch, Russian Academy of Sciences, 50, Pervomayskaya St., Syktyvkar, Russia
| |
Collapse
|
36
|
Krishnamoorthi S, Perotti LE, Borgstrom NP, Ajijola OA, Frid A, Ponnaluri AV, Weiss JN, Qu Z, Klug WS, Ennis DB, Garfinkel A. Simulation Methods and Validation Criteria for Modeling Cardiac Ventricular Electrophysiology. PLoS One 2014; 9:e114494. [PMID: 25493967 PMCID: PMC4262432 DOI: 10.1371/journal.pone.0114494] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Accepted: 11/07/2014] [Indexed: 01/24/2023] Open
Abstract
We describe a sequence of methods to produce a partial differential equation model of the electrical activation of the ventricles. In our framework, we incorporate the anatomy and cardiac microstructure obtained from magnetic resonance imaging and diffusion tensor imaging of a New Zealand White rabbit, the Purkinje structure and the Purkinje-muscle junctions, and an electrophysiologically accurate model of the ventricular myocytes and tissue, which includes transmural and apex-to-base gradients of action potential characteristics. We solve the electrophysiology governing equations using the finite element method and compute both a 6-lead precordial electrocardiogram (ECG) and the activation wavefronts over time. We are particularly concerned with the validation of the various methods used in our model and, in this regard, propose a series of validation criteria that we consider essential. These include producing a physiologically accurate ECG, a correct ventricular activation sequence, and the inducibility of ventricular fibrillation. Among other components, we conclude that a Purkinje geometry with a high density of Purkinje muscle junctions covering the right and left ventricular endocardial surfaces as well as transmural and apex-to-base gradients in action potential characteristics are necessary to produce ECGs and time activation plots that agree with physiological observations.
Collapse
Affiliation(s)
- Shankarjee Krishnamoorthi
- Department of Mechanical and Aerospace Engineering, University of California Los Angeles, Los Angeles, California, United States of America
| | - Luigi E. Perotti
- Department of Mechanical and Aerospace Engineering, University of California Los Angeles, Los Angeles, California, United States of America
| | - Nils P. Borgstrom
- Department of Bioengineering, University of California Los Angeles, Los Angeles, California, United States of America
| | - Olujimi A. Ajijola
- Department of Medicine (Cardiology), University of California Los Angeles, Los Angeles, California, United States of America
| | - Anna Frid
- Department of Integrative Biology and Physiology, University of California Los Angeles, Los Angeles, California, United States of America
| | - Aditya V. Ponnaluri
- Department of Mechanical and Aerospace Engineering, University of California Los Angeles, Los Angeles, California, United States of America
| | - James N. Weiss
- Department of Medicine (Cardiology), University of California Los Angeles, Los Angeles, California, United States of America
| | - Zhilin Qu
- Department of Medicine (Cardiology), University of California Los Angeles, Los Angeles, California, United States of America
| | - William S. Klug
- Department of Mechanical and Aerospace Engineering, University of California Los Angeles, Los Angeles, California, United States of America
| | - Daniel B. Ennis
- Department of Bioengineering, University of California Los Angeles, Los Angeles, California, United States of America
- Department of Radiological Sciences, University of California Los Angeles, Los Angeles, California, United States of America
| | - Alan Garfinkel
- Department of Medicine (Cardiology), University of California Los Angeles, Los Angeles, California, United States of America
- Department of Integrative Biology and Physiology, University of California Los Angeles, Los Angeles, California, United States of America
- * E-mail:
| |
Collapse
|
37
|
Rodriguez B. In Silico Organ Modelling in Predicting Efficacy and Safety of New Medicines. HUMAN-BASED SYSTEMS FOR TRANSLATIONAL RESEARCH 2014. [DOI: 10.1039/9781782620136-00219] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The development of new medicines faces important challenges due to difficulties in the assessment of their efficacy and their safety in the targeted human population. In silico approaches through the use of mathematical modelling and computer simulations are increasingly being used to overcome some of the limitations of current experimental methods used in the development of new medicines. This chapter describes state-of-the-art in silico approaches for the evaluation of the safety and efficacy of medicines targeting important causes of mortality such as cardiovascular disease. Firstly, we describe the in silico multi-scale mathematical models and simulation techniques required to describe drug-induced effects on physiological systems such as the heart from the subcellular to the whole organ level. Then we illustrate the power of in silico approaches used to augment experimental and clinical investigations, by providing the framework to unravel multi-scale mechanisms underlying variability in the response to medicines and to focus on effects in human rather than animal models. We devote the last part of the chapter to discussing the process of validation of in silico models and simulations, which is key in building up their credibility.
Collapse
Affiliation(s)
- Blanca Rodriguez
- Department of Computer Science, University of Oxford Parks Road Oxford OX1 3QD UK
| |
Collapse
|
38
|
Palamara S, Vergara C, Catanzariti D, Faggiano E, Pangrazzi C, Centonze M, Nobile F, Maines M, Quarteroni A. Computational generation of the Purkinje network driven by clinical measurements: the case of pathological propagations. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2014; 30:1558-77. [PMID: 25319252 DOI: 10.1002/cnm.2689] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2014] [Revised: 09/25/2014] [Accepted: 09/25/2014] [Indexed: 05/16/2023]
Abstract
To properly describe the electrical activity of the left ventricle, it is necessary to model the Purkinje fibers, responsible for the fast and coordinate ventricular activation, and their interaction with the muscular propagation. The aim of this work is to propose a methodology for the generation of a patient-specific Purkinje network driven by clinical measurements of the activation times related to pathological propagations. In this case, one needs to consider a strongly coupled problem between the network and the muscle, where the feedback from the latter to the former cannot be neglected as in a normal propagation. We apply the proposed strategy to data acquired on three subjects, one of them suffering from muscular conduction problems owing to a scar and the other two with a muscular pre-excitation syndrome (Wolff-Parkinson-White). To assess the accuracy of the proposed method, we compare the results obtained by using the patient-specific Purkinje network generated by our strategy with the ones obtained by using a non-patient-specific network. The results show that the mean absolute errors in the activation time is reduced for all the cases, highlighting the importance of including a patient-specific Purkinje network in computational models.
Collapse
Affiliation(s)
- Simone Palamara
- Modellistica e Calcolo Scientifico (MOX), Dipartimento di Matematica, Politecnico di Milano, Milan, Italy
| | | | | | | | | | | | | | | | | |
Collapse
|
39
|
Walmsley J, Mirams GR, Pitt-Francis J, Rodriguez B, Burrage K. Application of stochastic phenomenological modelling to cell-to-cell and beat-to-beat electrophysiological variability in cardiac tissue. J Theor Biol 2014; 365:325-36. [PMID: 25451525 PMCID: PMC4271765 DOI: 10.1016/j.jtbi.2014.10.029] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2014] [Revised: 10/21/2014] [Accepted: 10/23/2014] [Indexed: 01/08/2023]
Abstract
Variability in the action potential of isolated myocytes and tissue samples is observed in experimental studies. Variability is manifested as both differences in the action potential (AP) morphology between cells (extrinsic variability), and also ‘intrinsic’ or beat-to-beat variability of repolarization (BVR) in the AP duration of each cell. We studied the relative contributions of experimentally recorded intrinsic and extrinsic variability to dispersion of repolarization in tissue. We developed four cell-specific parameterizations of a phenomenological stochastic differential equation AP model exhibiting intrinsic variability using APs recorded from isolated guinea pig ventricular myocytes exhibiting BVR. We performed simulations in tissue using the four different model parameterizations in the presence and the absence of both intrinsic and extrinsic variability. We altered the coupling of the tissue to determine how inter-cellular coupling affected the dispersion of the AP duration in tissue. Both intrinsic and extrinsic variability were gradually revealed by reduction of tissue coupling. However, the recorded extrinsic variability between individual myocytes produced a greater degree of dispersion in repolarization in tissue than the intrinsic variability of each myocyte. We modelled inter-cell and beat-to-beat repolarization variability in cardiomyocytes. We coupled the cells together into cardiac tissue. Reducing tissue coupling increased repolarization dispersion in tissue. Inter-cell variability had a greater effect on repolarization dispersion.
Collapse
Affiliation(s)
- John Walmsley
- Department of Computer Science, University of Oxford, Oxford, United Kingdom.
| | - Gary R Mirams
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Joe Pitt-Francis
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Kevin Burrage
- Department of Computer Science, University of Oxford, Oxford, United Kingdom; School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
| |
Collapse
|
40
|
Zemzemi N, Bernabeu MO, Saiz J, Cooper J, Pathmanathan P, Mirams GR, Pitt-Francis J, Rodriguez B. Computational assessment of drug-induced effects on the electrocardiogram: from ion channel to body surface potentials. Br J Pharmacol 2013; 168:718-33. [PMID: 22946617 PMCID: PMC3579290 DOI: 10.1111/j.1476-5381.2012.02200.x] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2011] [Revised: 08/06/2012] [Accepted: 08/14/2012] [Indexed: 12/20/2022] Open
Abstract
Background and Purpose Understanding drug effects on the heart is key to safety pharmacology assessment and anti-arrhythmic therapy development. Here our goal is to demonstrate the ability of computational models to simulate the effect of drug action on the electrical activity of the heart, at the level of the ion-channel, cell, heart and ECG body surface potential. Experimental Approach We use the state-of-the-art mathematical models governing the electrical activity of the heart. A drug model is introduced using an ion channel conductance block for the hERG and fast sodium channels, depending on the IC50 value and the drug dose. We simulate the ECG measurements at the body surface and compare biomarkers under different drug actions. Key Results Introducing a 50% hERG-channel current block results in 8% prolongation of the APD90 and 6% QT interval prolongation, hERG block does not affect the QRS interval. Introducing 50% fast sodium current block prolongs the QRS and the QT intervals by 12% and 5% respectively, and delays activation times, whereas APD90 is not affected. Conclusions and Implications Both potassium and sodium blocks prolong the QT interval, but the underlying mechanism is different: for potassium it is due to APD prolongation; while for sodium it is due to a reduction of electrical wave velocity. This study shows the applicability of in silico models for the investigation of drug effects on the heart, from the ion channel to the ECG-based biomarkers.
Collapse
Affiliation(s)
- Nejib Zemzemi
- Department of Computer Science, University of Oxford, Oxford, UK.
| | | | | | | | | | | | | | | |
Collapse
|
41
|
Bishop MJ, Vigmond EJ, Plank G. The functional role of electrophysiological heterogeneity in the rabbit ventricle during rapid pacing and arrhythmias. Am J Physiol Heart Circ Physiol 2013; 304:H1240-52. [PMID: 23436328 PMCID: PMC3652087 DOI: 10.1152/ajpheart.00894.2012] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Electrophysiological heterogeneity in action potential recordings from healthy intact hearts remains highly variable and, where present, is almost entirely abolished at fast pacing rates. Consequently, the functional importance of intrinsic action potential duration (APD) heterogeneity in healthy ventricles, and particularly its role during rapidly activating reentrant arrhythmias, remain poorly understood. By incorporating both transmural and apicobasal APD heterogeneity within a biventricular rabbit computational model and comparing with an equivalent homogeneous model, we directly investigated the functional importance of intrinsic APD heterogeneity under fast pacing and arrhythmogenic protocols. Although differences in APD were significantly modulated at the tissue level during pacing and further reduced as pacing frequency increased, small differences were still noticeable. Such differences were further marginally accentuated/attenuated via electrotonic effects relative to wavefront propagation directions. The remaining small levels of APD heterogeneity under the fastest pacing frequencies resulted in arrhythmia initiation via heterogeneous conduction block, in contrast to complete block in the homogeneous model. Such induction mechanisms were more evident during premature stimuli at slower paced rhythms where intrinsic heterogeneity remained to a greater degree. During sustained arrhythmias, however, intrinsic heterogeneity made little difference to overall reentrant behavior, either visually, or in terms of duration, metrics quantifying filament/phase singularity dynamics, and global electrocardiogram characteristics. These findings suggest that, despite being important during arrhythmia initiation, intrinsic electrophysiological heterogeneity plays little functional role during rapid pacing and sustained arrhythmia dynamics in the healthy ventricle and thus questions the need to incorporate such detail in computational models when simulating rapid arrhythmias.
Collapse
Affiliation(s)
- Martin J Bishop
- Biomedical Engineering Department, Division of Imaging Sciences, King's College London, London, United Kingdom.
| | | | | |
Collapse
|
42
|
Wilhelms M, Rombach C, Scholz EP, Dossel O, Seemann G. Impact of amiodarone and cisapride on simulated human ventricular electrophysiology and electrocardiograms. Europace 2012; 14 Suppl 5:v90-v96. [DOI: 10.1093/europace/eus281] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
43
|
Potse M, Krause D, Bacharova L, Krause R, Prinzen FW, Auricchio A. Similarities and differences between electrocardiogram signs of left bundle-branch block and left-ventricular uncoupling. Europace 2012; 14 Suppl 5:v33-v39. [DOI: 10.1093/europace/eus272] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
44
|
Bayer JD, Blake RC, Plank G, Trayanova NA. A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models. Ann Biomed Eng 2012; 40:2243-54. [PMID: 22648575 DOI: 10.1007/s10439-012-0593-5] [Citation(s) in RCA: 265] [Impact Index Per Article: 22.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2012] [Accepted: 05/10/2012] [Indexed: 12/15/2022]
Abstract
Electrical waves traveling throughout the myocardium elicit muscle contractions responsible for pumping blood throughout the body. The shape and direction of these waves depend on the spatial arrangement of ventricular myocytes, termed fiber orientation. In computational studies simulating electrical wave propagation or mechanical contraction in the heart, accurately representing fiber orientation is critical so that model predictions corroborate with experimental data. Typically, fiber orientation is assigned to heart models based on Diffusion Tensor Imaging (DTI) data, yet few alternative methodologies exist if DTI data is noisy or absent. Here we present a novel Laplace-Dirichlet Rule-Based (LDRB) algorithm to perform this task with speed, precision, and high usability. We demonstrate the application of the LDRB algorithm in an image-based computational model of the canine ventricles. Simulations of electrical activation in this model are compared to those in the same geometrical model but with DTI-derived fiber orientation. The results demonstrate that activation patterns from simulations with LDRB and DTI-derived fiber orientations are nearly indistinguishable, with relative differences ≤6%, absolute mean differences in activation times ≤3.15 ms, and positive correlations ≥0.99. These results convincingly show that the LDRB algorithm is a robust alternative to DTI for assigning fiber orientation to computational heart models.
Collapse
Affiliation(s)
- J D Bayer
- Department of Biomedical Engineering, The Johns Hopkins University, Baltimore, MD, USA.
| | | | | | | |
Collapse
|
45
|
Carusi A, Burrage K, Rodríguez B. Bridging experiments, models and simulations: an integrative approach to validation in computational cardiac electrophysiology. Am J Physiol Heart Circ Physiol 2012; 303:H144-55. [PMID: 22582088 DOI: 10.1152/ajpheart.01151.2011] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Computational models in physiology often integrate functional and structural information from a large range of spatiotemporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and skepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace, and refine animal experiments. A fundamental requirement to fulfill these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations among experiments, models, and simulations in cardiac electrophysiology. We describe the processes, data, and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. We argue that validation is part of the whole MSE system and is contingent upon 1) understanding and coping with sources of biovariability; 2) testing and developing robust techniques and tools as a prerequisite to conducting physiological investigations; 3) defining and adopting standards to facilitate the interoperability of experiments, models, and simulations; 4) and understanding physiological validation as an iterative process that contributes to defining the specific aspects of cardiac electrophysiology the MSE system targets, rather than being only an external test, and that this is driven by advances in experimental and computational methods and the combination of both.
Collapse
|