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
|
Trayanova NA, Boyle PM, Nikolov PP. Personalized Imaging and Modeling Strategies for Arrhythmia Prevention and Therapy. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2018; 5:21-28. [PMID: 29546250 PMCID: PMC5847279 DOI: 10.1016/j.cobme.2017.11.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The goal of this article is to review advances in computational modeling of the heart, with a focus on recent non-invasive clinical imaging- and simulation-based strategies aimed at improving the diagnosis and treatment of patients with arrhythmias and structural heart disease. Following a brief overview of the field of computational cardiology, we present recent applications of the personalized virtual-heart approach in predicting the optimal targets for infarct-related ventricular tachycardia and atrial fibrillation ablation, and in determining risk of sudden cardiac death in myocardial infarction patients. The hope is that with such models at the patient bedside, therapies could be improved, invasiveness of diagnostic procedures minimized, and health-care costs reduced.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| | - Patrick M Boyle
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| | - Plamen P Nikolov
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| |
Collapse
|
3
|
Karathanos TV, Boyle PM, Trayanova NA. Light-based Approaches to Cardiac Arrhythmia Research: From Basic Science to Translational Applications. CLINICAL MEDICINE INSIGHTS-CARDIOLOGY 2016; 10:47-60. [PMID: 27840581 PMCID: PMC5094582 DOI: 10.4137/cmc.s39711] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Revised: 09/27/2016] [Accepted: 10/09/2016] [Indexed: 02/06/2023]
Abstract
Light has long been used to image the heart, but now it can be used to modulate its electrophysiological function. Imaging modalities and techniques have long constituted an indispensable part of arrhythmia research and treatment. Recently, advances in the fields of optogenetics and photodynamic therapy have provided scientists with more effective approaches for probing, studying and potentially devising new treatments for cardiac arrhythmias. This article is a review of research toward the application of these techniques. It contains (a) an overview of advancements in technology and research that have contributed to light-based cardiac applications and (b) a summary of current and potential future applications of light-based control of cardiac cells, including modulation of heart rhythm, manipulation of cardiac action potential morphology, quantitative analysis of arrhythmias, defibrillation and cardiac ablation.
Collapse
Affiliation(s)
- Thomas V. Karathanos
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Patrick M. Boyle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Natalia A. Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
- Department of Medicine, Johns Hopkins University, Baltimore, MD, USA
| |
Collapse
|
4
|
Arevalo HJ, Boyle PM, Trayanova NA. Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2016; 121:185-94. [PMID: 27334789 DOI: 10.1016/j.pbiomolbio.2016.06.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2016] [Accepted: 06/13/2016] [Indexed: 12/29/2022]
Abstract
Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.
Collapse
Affiliation(s)
- Hermenegild J Arevalo
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA; Simula Research Laboratory, Oslo, Norway
| | - Patrick M Boyle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
5
|
Trayanova NA, Chang KC. How computer simulations of the human heart can improve anti-arrhythmia therapy. J Physiol 2016; 594:2483-502. [PMID: 26621489 DOI: 10.1113/jp270532] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2015] [Accepted: 11/25/2015] [Indexed: 01/26/2023] Open
Abstract
Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Johns Hopkins School of Medicine, Baltimore, MD, 21205, USA.,Johns Hopkins University, Baltimore, MD, 21218, USA
| | - Kelly C Chang
- Johns Hopkins School of Medicine, Baltimore, MD, 21205, USA
| |
Collapse
|
6
|
Trayanova NA, Boyle PM, Arevalo HJ, Zahid S. Exploring susceptibility to atrial and ventricular arrhythmias resulting from remodeling of the passive electrical properties in the heart: a simulation approach. Front Physiol 2014; 5:435. [PMID: 25429272 PMCID: PMC4228852 DOI: 10.3389/fphys.2014.00435] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Accepted: 10/24/2014] [Indexed: 12/19/2022] Open
Abstract
Under diseased conditions, remodeling of the cardiac tissue properties (“passive properties”) takes place; these are aspects of electrophysiological behavior that are not associated with active ion transport across cell membranes. Remodeling of the passive electrophysiological properties most often results from structural remodeling, such as gap junction down-regulation and lateralization, fibrotic growth infiltrating the myocardium, or the development of an infarct scar. Such structural remodeling renders atrial or ventricular tissue as a major substrate for arrhythmias. The current review focuses on these aspects of cardiac arrhythmogenesis. Due to the inherent complexity of cardiac arrhythmias, computer simulations have provided means to elucidate interactions pertinent to this spatial scale. Here we review the current state-of-the-art in modeling atrial and ventricular arrhythmogenesis as arising from the disease-induced changes in the passive tissue properties, as well as the contributions these modeling studies have made to our understanding of the mechanisms of arrhythmias in the heart. Because of the rapid advance of structural imaging methodologies in cardiac electrophysiology, we chose to present studies that have used such imaging methodologies to construct geometrically realistic models of cardiac tissue, or the organ itself, where the regional remodeling properties of the myocardium can be represented in a realistic way. We emphasize how the acquired knowledge can be used to pave the way for clinical applications of cardiac organ modeling under the conditions of structural remodeling.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University Baltimore, MD, USA
| | - Patrick M Boyle
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University Baltimore, MD, USA
| | - Hermenegild J Arevalo
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University Baltimore, MD, USA
| | - Sohail Zahid
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University Baltimore, MD, USA
| |
Collapse
|
7
|
Trayanova NA, Rantner LJ. New insights into defibrillation of the heart from realistic simulation studies. Europace 2014; 16:705-13. [PMID: 24798960 PMCID: PMC4010179 DOI: 10.1093/europace/eut330] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2013] [Accepted: 09/17/2013] [Indexed: 11/12/2022] Open
Abstract
Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices, constitutes the most important means of combating sudden cardiac death. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart, particularly under diseased conditions, is a promising approach to achieve an optimal therapy. The aim of this article is to assess the current state-of-the-art in whole-heart defibrillation modelling, focusing on major insights that have been obtained using defibrillation models, primarily those of realistic heart geometry and disease remodelling. The article showcases the contributions that modelling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modelling of the heart (i.e. cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level, and has the potential to contribute to the betterment of the clinical practice of defibrillation.
Collapse
Affiliation(s)
- Natalia A. Trayanova
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, 3400 N Charles Street, 216 Hackerman Hall, Baltimore, MD 21218, USA
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, Baltimore, MD 21218, USA
| | - Lukas J. Rantner
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, 3400 N Charles Street, 216 Hackerman Hall, Baltimore, MD 21218, USA
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, Baltimore, MD 21218, USA
| |
Collapse
|
8
|
Weinberg SH, Chang KC, Zhu R, Tandri H, Berger RD, Trayanova NA, Tung L. Defibrillation success with high frequency electric fields is related to degree and location of conduction block. Heart Rhythm 2013; 10:740-8. [PMID: 23354078 DOI: 10.1016/j.hrthm.2013.01.016] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/18/2012] [Indexed: 11/27/2022]
Abstract
BACKGROUND We recently demonstrated that high frequency alternating current (HFAC) electric fields can reversibly block propagation in the heart by inducing an oscillating, elevated transmembrane potential (Vm) that maintains myocytes in a refractory state for the field duration and can terminate arrhythmias, including ventricular fibrillation (VF). OBJECTIVES To quantify and characterize conduction block (CB) induced by HFAC fields and to determine whether the degree of CB can be used to predict defibrillation success. METHODS Optical mapping was performed in adult guinea pig hearts (n = 14), and simulations were performed in an anatomically accurate rabbit ventricular model. HFAC fields (50-500 Hz) were applied to the ventricles. A novel power spectrum metric of CB-the loss of spectral power in the 1-30 Hz range, termed loss of conduction power (LCP)-was assessed during the HFAC field and compared with defibrillation success and VF vulnerability. RESULTS LCP increased with field strength and decreased with frequency. Optical mapping experiments conducted on the epicardial surface showed that LCP and the size of CB regions were significantly correlated with VF initiation and termination. In simulations, subsurface myocardial LCP and CB sizes were more closely correlated with VF termination than surface values. Multilinear regression analysis of simulation results revealed that while CB on both the surface and the subsurface myocardium was predictive, subsurface myocardial CB was the better predictor of defibrillation success. CONCLUSIONS HFAC fields induce a field-dependent state of CB, and defibrillation success is related to the degree and location of the CB.
Collapse
Affiliation(s)
- Seth H Weinberg
- Department of Applied Science, College of William and Mary, Williamsburg, Virginia, USA
| | | | | | | | | | | | | |
Collapse
|
9
|
Trayanova N, Constantino J, Ashihara T, Plank G. Modeling defibrillation of the heart: approaches and insights. IEEE Rev Biomed Eng 2012; 4:89-102. [PMID: 22273793 DOI: 10.1109/rbme.2011.2173761] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices (ICDs), constitutes the most important means of combating sudden cardiac death. While ICD therapy has proved to be efficient and reliable, defibrillation is a traumatic experience. Thus, research on defibrillation mechanisms, particularly aimed at lowering defibrillation voltage, remains an important topic. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart is the most promising approach to achieve this goal. The aim of this paper is to assess the current state-of-the-art in ventricular defibrillation modeling, focusing on both numerical modeling approaches and major insights that have been obtained using defibrillation models, primarily those of realistic ventricular geometry. The paper showcases the contributions that modeling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modeling of the heart (i.e., cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level and has the potential to contribute to the betterment of the clinical practice of defibrillation.
Collapse
Affiliation(s)
- Natalia Trayanova
- Department of Biomedical Engineering and Institute for Computational Medicine, The Johns Hopkins University, Baltimore, MD 20218, USA.
| | | | | | | |
Collapse
|
10
|
SUZUKI TOHRU, SATO SHUNSUKE, OHE TOHRU, SUZUKI RYOJI, KAJIYA FUMIHIKO. ANALYSIS OF THE VIRTUAL ELECTRODE PHENOMENA USING BIDOMAIN MODEL: BASIC CHARACTERISTICS FOR PASSIVE MEMBRANE. J MECH MED BIOL 2011. [DOI: 10.1142/s0219519406002023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The virtual electrode (VE) has been recognized as an important factor for success or failure of cardiac defibrillation. Many researches have been performed to study characteristics of the VE. However, there are some questions which remain unanswered. In this study, we developed a simulator to solve a three-dimensional bidomain model and performed several simulations to elucidate the basic characteristics of VE in a simplified cardiac tissue with passive membrane when a constant unipolar cathodal stimulus was applied. The results showed that for smaller electrodes, VE has a typical dog-bone shaped virtual cathode (VC) and two egg-shaped virtual anodes (VAs). The distributions both in intra- and extracellular potentials have concentric ellipsoidal isosurfaces, but their ellipticities are subtly different, producing VE. For larger electrodes, VC becomes larger and has a flat-dish shape rather than dog-bone, and VA becomes smaller and also flattens and collapses. The peak values of VE are larger for smaller electrodes, but their time courses show similar tendency among the different sized electrodes. The change of stimulus strength and polarity only affects the magnitude of VE in a linear manner and the distribution pattern is unchanged. These results provide us fundamental knowledge about VE.
Collapse
Affiliation(s)
- TOHRU SUZUKI
- Department of Cardiovascular Physiology, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama 700-8558, Japan
| | - SHUNSUKE SATO
- Department of Physical Therapy, Aino University, Ibaraki, Osaka 567-0012, Japan
| | - TOHRU OHE
- Department of Cardiovascular Medicine, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama 700-8558, Japan
| | - RYOJI SUZUKI
- Human Information System Laboratory, Kanazawa Institute of Technology, Hakusan, Ishikawa 924-0838, Japan
| | - FUMIHIKO KAJIYA
- Department of Cardiovascular Physiology, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama 700-8558, Japan
| |
Collapse
|
11
|
Trayanova NA, Rice JJ. Cardiac electromechanical models: from cell to organ. Front Physiol 2011; 2:43. [PMID: 21886622 PMCID: PMC3154390 DOI: 10.3389/fphys.2011.00043] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2011] [Accepted: 07/12/2011] [Indexed: 11/13/2022] Open
Abstract
The heart is a multiphysics and multiscale system that has driven the development of the most sophisticated mathematical models at the frontiers of computational physiology and medicine. This review focuses on electromechanical (EM) models of the heart from the molecular level of myofilaments to anatomical models of the organ. Because of the coupling in terms of function and emergent behaviors at each level of biological hierarchy, separation of behaviors at a given scale is difficult. Here, a separation is drawn at the cell level so that the first half addresses subcellular/single-cell models and the second half addresses organ models. At the subcellular level, myofilament models represent actin–myosin interaction and Ca-based activation. The discussion of specific models emphasizes the roles of cooperative mechanisms and sarcomere length dependence of contraction force, considered to be the cellular basis of the Frank–Starling law. A model of electrophysiology and Ca handling can be coupled to a myofilament model to produce an EM cell model, and representative examples are summarized to provide an overview of the progression of the field. The second half of the review covers organ-level models that require solution of the electrical component as a reaction–diffusion system and the mechanical component, in which active tension generated by the myocytes produces deformation of the organ as described by the equations of continuum mechanics. As outlined in the review, different organ-level models have chosen to use different ionic and myofilament models depending on the specific application; this choice has been largely dictated by compromises between model complexity and computational tractability. The review also addresses application areas of EM models such as cardiac resynchronization therapy and the role of mechano-electric coupling in arrhythmias and defibrillation.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Department of Biomedical Engineering and Institute for Computational Medicine, Johns Hopkins University Baltimore, MD, USA
| | | |
Collapse
|
12
|
Bernabeu MO, Wallman M, Rodriguez B. Shock-induced arrhythmogenesis in the human heart: A computational modelling study. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2011; 2010:760-3. [PMID: 21095904 DOI: 10.1109/iembs.2010.5626338] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Electrical defibrillation by application of a strong shock to the heart is the only effective treatment against lethal cardiac arrhythmias such as ventricular fibrillation. A large body of experimental and computational research has been devoted to understanding shock-induced effects on the heart in an attempt to improve defibrillation efficacy. However, most of the research has been performed in small animal hearts, and in particular rabbits. The difference in size between rabbits and humans might limit the extrapolation of the results to the clinical setting. In this paper, we present, for the first time, computer simulations of shock-induced effects on a human ventricular model with realistic ion channel dynamics and fibre architecture. Bidomain simulations using the human ventricular model were performed using the Chaste open source simulation package. The parallel performance of the software package was highly improved in order to meet the computational requirements of these kind of studies.
Collapse
Affiliation(s)
- Miguel O Bernabeu
- Oxford University Computing Laboratory, Wolfson Building, Parks Rd, OX1 3QD, United Kingdom.
| | | | | |
Collapse
|
13
|
How hyperpolarization and the recovery of excitability affect propagation through a virtual anode in the heart. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2011; 2011:375059. [PMID: 21331264 PMCID: PMC3038668 DOI: 10.1155/2011/375059] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2010] [Accepted: 12/16/2010] [Indexed: 11/17/2022]
Abstract
Researchers have suggested that the fate of a shock-induced wave front at the edge of a "virtual anode" (a region hyperpolarized by the shock) is a key factor determining success or failure during defibrillation of the heart. In this paper, we use a simple one-dimensional computer model to examine propagation speed through a hyperpolarized region. Our goal is to test the hypothesis that rapid propagation through a virtual anode can cause failure of propagation at the edge of the virtual anode. The calculations support this hypothesis and suggest that the time constant of the sodium inactivation gate is an important parameter. These results may be significant in understanding the mechanism of the upper limit of vulnerability.
Collapse
|
14
|
Abstract
Recent developments in cardiac simulation have rendered the heart the most highly integrated example of a virtual organ. We are on the brink of a revolution in cardiac research, one in which computational modeling of proteins, cells, tissues, and the organ permit linking genomic and proteomic information to the integrated organ behavior, in the quest for a quantitative understanding of the functioning of the heart in health and disease. The goal of this review is to assess the existing state-of-the-art in whole-heart modeling and the plethora of its applications in cardiac research. General whole-heart modeling approaches are presented, and the applications of whole-heart models in cardiac electrophysiology and electromechanics research are reviewed. The article showcases the contributions that whole-heart modeling and simulation have made to our understanding of the functioning of the heart. A summary of the future developments envisioned for the field of cardiac simulation and modeling is also presented. Biophysically based computational modeling of the heart, applied to human heart physiology and the diagnosis and treatment of cardiac disease, has the potential to dramatically change 21st century cardiac research and the field of cardiology.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
15
|
Kanlop N, Shinlapawittayatorn K, Sungnoon R, Weerateerangkul P, Chattipakorn S, Chattipakorn N. Cilostazol attenuates ventricular arrhythmia induction and improves defibrillation efficacy in swine. Can J Physiol Pharmacol 2010; 88:422-8. [PMID: 20555410 DOI: 10.1139/y09-127] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Previous reports demonstrated that cilostazol, a phosphodiesterase 3 inhibitor, affected cellular electrophysiology and reduced episodes of ventricular fibrillation (VF) in patients with Brugada syndrome. However, its effects on VF induction and defibrillation efficacy have never been investigated. We tested the hypothesis that cilostazol increases the VF threshold (VFT) and decreases the upper limit of vulnerability (ULV) and the defibrillation threshold (DFT). A total of 48 pigs were randomly assigned to defibrillation and VF induction studies. The diastolic pacing threshold (DPT), VFT, ULV, DFT, and effective refractory period were determined before and after the infusion of cilostazol at 6 mg/kg, 3 mg/kg, or vehicle. The DPT was significantly increased after administration of 3 and 6 mg/kg cilostazol. The ULV and DFT were significantly decreased after administration of 6 mg/kg cilostazol only. The ULV in the 6 mg/kg group had 12% lower peak voltage and 25% lower total energy, and the DFT had 13% lower peak voltage and 25% lower total energy. The VFT was not altered in any experimental group. This study shows that cilostazol administration significantly increased the DPT, which was associated with significantly reduced DFT and ULV.
Collapse
Affiliation(s)
- Natnicha Kanlop
- Cardiac Electrophysiology Unit, Department of Physiology, and Cardiac Electrophysiology Research and Training Center, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand
| | | | | | | | | | | |
Collapse
|
16
|
Trayanova NA, Tice BM. Integrative computational models of cardiac arrhythmias -- simulating the structurally realistic heart. ACTA ACUST UNITED AC 2009; 6:85-91. [PMID: 20628585 DOI: 10.1016/j.ddmod.2009.08.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Simulation of cardiac electrical function, and specifically, simulation aimed at understanding the mechanisms of cardiac rhythm disorders, represents an example of a successful integrative multiscale modeling approach, uncovering emergent behavior at the successive scales in the hierarchy of structural complexity. The goal of this article is to present a review of the integrative multiscale models of realistic ventricular structure used in the quest to understand and treat ventricular arrhythmias. It concludes with the new advances in image-based modeling of the heart and the promise it holds for the development of individualized models of ventricular function in health and disease.
Collapse
|
17
|
Vigmond E, Vadakkumpadan F, Gurev V, Arevalo H, Deo M, Plank G, Trayanova N. Towards predictive modelling of the electrophysiology of the heart. Exp Physiol 2009; 94:563-77. [PMID: 19270037 DOI: 10.1113/expphysiol.2008.044073] [Citation(s) in RCA: 67] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The simulation of cardiac electrical function is an example of a successful integrative multiscale modelling approach that is directly relevant to human disease. Today we stand at the threshold of a new era, in which anatomically detailed, tomographically reconstructed models are being developed that integrate from the ion channel to the electromechanical interactions in the intact heart. Such models hold high promise for interpretation of clinical and physiological measurements, for improving the basic understanding of the mechanisms of dysfunction in disease, such as arrhythmias, myocardial ischaemia and heart failure, and for the development and performance optimization of medical devices. The goal of this article is to present an overview of current state-of-art advances towards predictive computational modelling of the heart as developed recently by the authors of this article. We first outline the methodology for constructing electrophysiological models of the heart. We then provide three examples that demonstrate the use of these models, focusing specifically on the mechanisms for arrhythmogenesis and defibrillation in the heart. These include: (1) uncovering the role of ventricular structure in defibrillation; (2) examining the contribution of Purkinje fibres to the failure of the shock; and (3) using magnetic resonance imaging reconstructed heart models to investigate the re-entrant circuits formed in the presence of an infarct scar.
Collapse
Affiliation(s)
- Edward Vigmond
- Department of Electrical and Computer Engineering, University of Calgary, Calgary, Alberta, Canada
| | | | | | | | | | | | | |
Collapse
|
18
|
Abstract
BACKGROUND The strongest shock that induces reentry in the heart is the upper limit of vulnerability (ULV). In order to understand defibrillation, one must know what causes the ULV. OBJECTIVE The goal of this study was to examine the mechanism of the upper limit of vulnerability. METHODS Numerical simulations of cardiac tissue were performed using the bidomain model. An S2 shock was applied during the refractory period of the S1 action potential, and results using a smooth curving fiber geometry were compared with results using a smooth plus random fiber geometry. RESULTS When using a smooth fiber geometry only, no ULV was observed. However, when a random fiber geometry was included, the ULV was present. The difference arises from the fate of the shock-induced break wave front when it reaches the edge of the tissue hyperpolarized by the shock (the virtual anode). CONCLUSION Our numerical simulations suggest that local heterogeneities throughout the tissue may be crucial for determining the fate of the shock-induced wave front at the edge of the virtual anode, and therefore play an important role in the mechanism underlying the ULV.
Collapse
Affiliation(s)
- Nachaat Mazeh
- Department of Physics, Oakland University, Rochester, Michigan 48309, USA
| | | |
Collapse
|
19
|
Brennan T, Fink M, Rodriguez B. Multiscale modelling of drug-induced effects on cardiac electrophysiological activity. Eur J Pharm Sci 2008; 36:62-77. [PMID: 19061955 DOI: 10.1016/j.ejps.2008.09.013] [Citation(s) in RCA: 75] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2008] [Accepted: 09/08/2008] [Indexed: 01/09/2023]
Abstract
Many drugs fail in the clinical trials and therefore do not reach the market due to adverse effects on cardiac electrical function. This represents a growing concern for both regulatory and pharmaceutical agencies as it translates into important socio-economic costs. Drugs affecting cardiac activity come from diverse pharmacological groups and their interaction with cardiac electrophysiology can result in increased risk of potentially life threatening arrhythmias, such as Torsade de Pointes. The mechanisms of drug interaction with the heart are very complex and the effects span from the ion channel to the whole organ level. This makes their investigation using solely experimental in vitro and in vivo techniques very difficult. Computational modelling of cardiac electrophysiological behaviour has provided insight into the mechanisms of cardiac arrhythmogenesis, with high spatio-temporal resolution, from the ion channel to the whole organ level. It therefore represents a powerful tool in investigating mechanisms of drug-induced changes in cardiac behaviour and in their pro-arrhythmic potential. This article presents a comprehensive review of the recent advances in detailed models of drug action on cardiac electrophysiological activity.
Collapse
Affiliation(s)
- T Brennan
- Department of Engineering Science, University of Oxford, Oxford, United Kingdom
| | | | | |
Collapse
|
20
|
Tice BM, Rodríguez B, Eason J, Trayanova N. Mechanistic investigation into the arrhythmogenic role of transmural heterogeneities in regional ischaemia phase 1A. Europace 2008; 9 Suppl 6:vi46-58. [PMID: 17959693 DOI: 10.1093/europace/eum204] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
AIMS Studies of arrhythmogenesis during ischemia have focused primarily on reentrant mechanisms manifested on the epicardial surface. The goal of this study was to use a physiologically-accurate model of acute regional ischemia phase 1A to determine the contribution of ischaemia-induced transmural electrophysiological heterogeneities to arrhythmogenesis following left anterior descending artery occlusion. METHODS AND RESULTS A slice through a geometrical model of the rabbit ventricles was extracted and a model of regional ischaemia developed. The model included a central ischaemic zone incorporating transmural gradients of I(K(ATP)) activation and [K+]o, surrounded by ischaemic border zones (BZs), with the degree of ischaemic effects varied to represent progression of ischaemia 2-10 min post-occlusion. Premature stimulation was applied over a range of coupling intervals to induce re-entry. The presence of ischaemic BZs and a transmural gradient in I(K(ATP)) activation provided the substrate for re-entrant arrhythmias. Increased dispersion of refractoriness and conduction velocity in the BZs with time post-occlusion led to a progressive increase in arrhythmogenesis. In the absence of a transmural gradient of I(K(ATP)) activation, re-entry was rarely sustained. CONCLUSION Knowledge of the mechanism by which specific electrophysiological heterogeneities underlie arrhythmogenesis during acute ischaemia could be useful in developing preventative treatments for patients at risk of coronary vascular disease.
Collapse
Affiliation(s)
- Brock M Tice
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, CSEB 216, Baltimore, MD 21218, USA
| | | | | | | |
Collapse
|
21
|
Bishop MJ, Rodriguez B, Trayanova N, Gavaghan DJ. Modulation of shock-end virtual electrode polarisation as a direct result of 3D fluorescent photon scattering. CONFERENCE PROCEEDINGS : ... ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL CONFERENCE 2008; 2006:1556-9. [PMID: 17946049 DOI: 10.1109/iembs.2006.259243] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Due to the large transmural variation in transmembrane potential following the application of strong electric shocks, it is thought that fluorescent photon scattering from depth plays a significant role in optical signal modulation at shock-end. For the first time, a model of photon scattering is used to accurately synthesize fluorescent signals over the irregular geometry of the rabbit ventricles following the application of such strong shocks. A bidomain representation of electrical activity is combined with finite element solutions to the photon diffusion equation, simulating both the excitation and emission processes, over an anatomically-based model of rabbit ventricular geometry and fiber orientation. Photon scattering from within a 3D volume beneath the epicardial optical recording site is shown to transduce differences in transmembrane potential within this volume through the myocardial wall. This leads directly to a significantly modulated optical signal response with respect to that predicted by the bidomain simulations, distorting epicardial virtual electrode polarization produced at shock-end. Furthermore, we show that this degree of distortion is very sensitive to the optical properties of the tissue, an important variable to consider during experimental mapping set-ups. These findings provide an essential first-step in aiding the interpretation of experimental optical mapping recordings following strong defibrillation shocks.
Collapse
Affiliation(s)
- M J Bishop
- Computational Biology Group, University of Oxford, Oxford, UK
| | | | | | | |
Collapse
|
22
|
Li W, Gurev V, McCulloch AD, Trayanova NA. The role of mechanoelectric feedback in vulnerability to electric shock. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2008; 97:461-78. [PMID: 18374394 DOI: 10.1016/j.pbiomolbio.2008.02.020] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
Experimental and clinical studies have shown that ventricular dilatation is associated with increased arrhythmogenesis and elevated defibrillation threshold; however, the underlying mechanisms remain poorly understood. The goal of the present study was to test the hypothesis that (1) stretch-activated channel (SAC) recruitment and (2) geometrical deformations in organ shape and fiber architecture lead to increased arrhythmogenesis by electric shocks following acute ventricular dilatation. To elucidate the contribution of these two factors, the study employed, for the first time, a combined electro-mechanical simulation approach. Acute dilatation was simulated in a model of rabbit ventricular mechanics by raising the LV end-diastolic pressure from 0.6 (control) to 4.2 kPa (dilated). The output of the mechanics model was used in the electrophysiological model. Vulnerability to shocks was examined in the control, the dilated ventricles, and in the dilated ventricles that also incorporated currents through SAC as a function of local strain, by constructing vulnerability grids. Results showed that dilatation-induced deformation alone decreased upper limit of vulnerability (ULV) slightly and did not result in increased vulnerability. With SAC recruitment in the dilated ventricles, the number of shock-induced arrhythmia episodes increased by 37% (from 41 to 56) and the lower limit of vulnerability (LLV) decreased from 9 to 7 V/cm, while ULV did not change. The heterogeneous activation of SAC caused by the heterogeneous fiber strain in the ventricular walls was the main reason for increased vulnerability to electric shocks since it caused dispersion of electrophysiological properties in the tissue, resulting in postshock unidirectional block and establishment of reentry.
Collapse
Affiliation(s)
- Weihui Li
- Department of Pathology, Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, USA
| | | | | | | |
Collapse
|
23
|
Yang MJ, Tran DX, Weiss JN, Garfinkel A, Qu Z. The pinwheel experiment revisited: effects of cellular electrophysiological properties on vulnerability to cardiac reentry. Am J Physiol Heart Circ Physiol 2007; 293:H1781-90. [PMID: 17586622 DOI: 10.1152/ajpheart.00014.2007] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
In normal heart, ventricular fibrillation can be induced by a single properly timed strong electrical or mechanical stimulus. A mechanism first proposed by Winfree and coined the “pinwheel experiment” emphasizes the timing and strength of the stimulus in inducing figure-of-eight reentry. However, the effects of cellular electrophysiological properties on vulnerability to reentry in the pinwheel scenario have not been investigated. In this study, we extend Winfree's pinwheel experiment to show how the vulnerability to reentry is affected by the graded action potential responses induced by a strong premature stimulus, action potential duration (APD), and APD restitution in simulated monodomain homogeneous two-dimensional tissue. We find that a larger graded response, longer APD, or steeper APD restitution slope reduces the vulnerable window of reentry. Strong graded responses and long APD promote tip-tip interactions at long coupling intervals, causing the two initiated spiral wave tips to annihilate. Steep APD restitution promotes wave front-wave back interaction, causing conduction block in the central common pathway of figure-of-eight reentry. We derive an analytical treatment that shows good agreement with numerical simulation results.
Collapse
Affiliation(s)
- Ming-Jim Yang
- Cardiovascular Research Laboratory, David Geffen School of Medicine, University of California, Los Angeles 90095, USA
| | | | | | | | | |
Collapse
|
24
|
Deo M, Bauer S, Plank G, Vigmond E. Reduced-order preconditioning for bidomain simulations. IEEE Trans Biomed Eng 2007; 54:938-42. [PMID: 17518292 DOI: 10.1109/tbme.2006.889203] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Simulations of the bidomain equations involve solving large, sparse, linear systems of the form Ax = b. Being an initial value problems, it is solved at every time step. Therefore, efficient solvers are essential to keep simulations tractable. Iterative solvers, especially the preconditioned conjugate gradient (PCG) method, are attractive since memory demands are minimized compared to direct methods, albeit at the cost of solution speed. However, a proper preconditioner can drastically speed up the solution process by reducing the number of iterations. In this paper, a novel preconditioner for the PCG method based on system order reduction using the Arnoldi method (A-PCG) is proposed. Large order systems, generated during cardiac bidomain simulations employing a finite element method formulation, are solved with the A-PCG method. Its performance is compared with incomplete LU (ILU) preconditioning. Results indicate that the A-PCG estimates an approximate solution considerably faster than the ILU, often within a single iteration. To reduce the computational demands in terms of memory and run time, the use of a cascaded preconditioner was suggested. The A-PCG was applied to quickly obtain an approximate solution, and subsequently a cheap iterative method such as successive overrelaxation (SOR) is applied to further refine the solution to arrive at a desired accuracy. The memory requirements are less than those of direct LU but more than ILU method. The proposed scheme is shown to yield significant speedups when solving time evolving systems.
Collapse
Affiliation(s)
- Makarand Deo
- Department of Electrical and Computer Engineering, University of Calgary, 2500 University Drive NW, Calgary, AB T2N 1N4, Canada.
| | | | | | | |
Collapse
|
25
|
Deo M, Vigmond E. Arnoldi preconditioning for solving large linear biomedical systems. CONFERENCE PROCEEDINGS : ... ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL CONFERENCE 2007; 2005:2914-7. [PMID: 17282853 DOI: 10.1109/iembs.2005.1617084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Simulations of biomedical systems often involve solving large, sparse, linear systems of the form Ax = b. In initial value problems, this system is solved at every time step, so a quick solution is essential for tractability. Iterative solvers, especially preconditioned conjugate gradient, are attractive since memory demands are minimized compared to direct methods, albeit at a cost of solution speed. A proper preconditioner can drastically reduce computation and remains an area of active research. In this paper, we propose a novel preconditioner based on system order reduction using the Arnoldi method. Systems of orders up to a million, generated from a finite element method formulation of the elliptic portion of the bidomain equations, are solved with the new preconditioner and performance is compared with that of other preconditioners. Results indicate that the new method converges considerably faster, often within a single iteration. It also uses less memory than an incomplete LU decomposition (ILU). For solving a system repeatedly, the Arnoldi transformation must be continually recomputed, unlike ILU, but this can be done quickly. In conclusion, for solving a system once, the Arnoldi preconditioner offers a greatly reduced solution time, and for repeated solves, will still be faster than an ILU preconditioner.
Collapse
Affiliation(s)
- Makarand Deo
- Dept. of Electr. & Comput. Eng., Calgary Univ., Alta
| | | |
Collapse
|
26
|
Abstract
Myocardial ischemia is one of the main causes of sudden cardiac death, with 80% of victims suffering from coronary heart disease. In acute myocardial ischemia, the obstruction of coronary flow leads to the interruption of oxygen flow, glucose, and washout in the affected tissue. Cellular metabolism is impaired and severe electrophysiological changes in ionic currents and concentrations ensue, which favor the development of lethal cardiac arrhythmias such as ventricular fibrillation. Due to the burden imposed by ischemia in our societies, a large body of research has attempted to unravel the mechanisms of initiation, sustenance, and termination of cardiac arrhythmias in acute ischemia, but the rapidity and complexity of ischemia-induced changes as well as the limitations in current experimental techniques have hampered evaluation of ischemia-induced alterations in cardiac electrical activity and understanding of the underlying mechanisms. Over the last decade, computer simulations have demonstrated the ability to provide insight, with high spatiotemporal resolution, into ischemic abnormalities in cardiac electrophysiological behavior from the ionic channel to the whole organ. This article aims to review and summarize the results of these studies and to emphasize the role of computer simulations in improving the understanding of ischemia-related arrhythmias and how to efficiently terminate them.
Collapse
Affiliation(s)
- Blanca Rodríguez
- Oxford University Computing Laboratory, Wolfson Building, Parks Road, Oxford, UK.
| | | | | |
Collapse
|
27
|
Li W, Kohl P, Trayanova N. Myocardial ischemia lowers precordial thump efficacy: an inquiry into mechanisms using three-dimensional simulations. Heart Rhythm 2006; 3:179-86. [PMID: 16443533 DOI: 10.1016/j.hrthm.2005.10.033] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/19/2005] [Accepted: 10/31/2005] [Indexed: 11/20/2022]
Abstract
BACKGROUND Precordial thump is the first International Liaison Committee on Resuscitation-prescribed procedure for advanced life support in witnessed cardiac arrest. Success rates vary and, according to clinical evidence, are significantly reduced under ischemic conditions. OBJECTIVES The purpose of this study was to elucidate the mechanisms involved in termination of ventricular tachycardia (VT) by precordial thump and its decreased rate of success in ischemia using a three-dimensional realistic model of electrical behavior in the rabbit ventricles. METHODS The electrophysiologic effect of precordial thump was represented by recruitment of mechanosensitive channels in the regions affected by precordial thump. In normoxia, precordial thump opened cation nonselective stretch-activated channels (SAC-NS, reversal potential -20 mV). In ischemia, precordial thump was assumed to additionally activate ATP-sensitive K+ (K(ATP)) channels (reversal potential -95 mV). Ten randomly selected cases of VT were used, and for each case the effect of precordial thump on VT was examined for normoxia and under ischemic conditions of varying severity. RESULTS Precordial thump was found to have a 60% success rate in normoxia and 30% in ischemia. Results demonstrate that precordial thump-induced SAC-NS opening in normoxia reduced heterogeneity in transmembrane potential by partially repolarizing excited tissue and depolarizing resting myocardium, potentially causing foci of excitation that eradicate the excitable gap, thus facilitating VT termination. Decreased precordial thump efficacy in ischemia was caused by recruitment of K(ATP), which diminished the depolarizing effect of SAC-NS on resting tissue and caused pronounced action potential shortening, thus facilitating establishment of reentry. CONCLUSION This study provides mechanistic insight into precordial thump mechanisms and its reduced clinical utility in patients with myocardial ischemia.
Collapse
Affiliation(s)
- Weihui Li
- Biomedical Engineering, Tulane University, New Orleans, Louisiana 70118, USA.
| | | | | |
Collapse
|
28
|
Beaudoin DL, Roth BJ. The effect of the fiber curvature gradient on break excitation in cardiac tissue. PACING AND CLINICAL ELECTROPHYSIOLOGY: PACE 2006; 29:496-501. [PMID: 16689845 DOI: 10.1111/j.1540-8159.2006.00382.x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
BACKGROUND Break excitation has been hypothesized as a mechanism for the initiation of reentry in cardiac tissue. One way break excitation can occur is by virtual electrodes formed due to a curving fiber geometry. In this article, we are concerned with the relationship between the peak gradient of fiber curvature and the threshold for break stimulation and the initiation of reentry. METHODS We calculate the maximum gradient of fiber curvature for different scales of fiber geometry in a constant tissue size (20x20 mm), and also examine the mechanisms by which reentry initiation fails. RESULTS For small peak gradients, reentry fails because break excitation does not occur. For larger peak gradients, reentry fails because break excitation fails to develop into full-scale reentry. For strong stimuli above the upper limit of vulnerability, reentry fails because the break excitation propagates through the hyperpolarized region and then encounters refractory tissue, causing the wave front to die.
Collapse
|
29
|
Abstract
Despite its critical role in restoring cardiac rhythm and thus in saving human life, cardiac defibrillation remains poorly understood. Further mechanistic inquiry is hampered by the inability of presently available experimental techniques to resolve, with sufficient accuracy, electrical behaviour confined to the depth of the ventricles. The objective of this review article is to demonstrate that realistic 3-D simulations of the ventricular defibrillation process in close conjunction with experimental observations are capable of bringing a new level of understanding of the electrical events that ensue from the interaction between fibrillating myocardium and applied shock. The article does this by reviewing the results of two studies, one on vulnerability to electric shocks and another on defibrillation. An overview of the modelling tools used in these studies is also provided.
Collapse
Affiliation(s)
- Natalia Trayanova
- Department of Biomedical Engineering, 500 Lindy Boggs Center, Suite 500, Tulane University, New Orleans, LA 70118, USA.
| |
Collapse
|
30
|
Deo M, Bauer S, Plank G, Vigmond E. Accelerating large cardiac bidomain simulations by arnoldi preconditioning. CONFERENCE PROCEEDINGS : ... ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL CONFERENCE 2006; 2006:3923-3926. [PMID: 17946209 DOI: 10.1109/iembs.2006.259271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Bidomain simulations of cardiac systems often in volve solving large, sparse, linear systems of the form Ax=b. These simulations are computationally very expensive in terms of run time and memory requirements. Therefore, efficient solvers are essential to keep simulations tractable. In this paper, an efficient preconditioner for the conjugate gradient (CG) method based on system order reduction using the Arnoldi method (A-PCG) is explained. Large order systems generated during cardiac bidomain simulations using a finite element method formulation, are solved using the A-PCG method. Its performance is compared with incomplete LU (ILU) preconditioning. Results indicate that the A-PCG estimates an approximate solution considerably faster than the ILU, often within a single iteration. To reduce the computational demands in terms of memory and run time, the use of a cascaded preconditioner is suggested. The A-PCG can be applied to quickly obtain an approximate solution, subsequently a cheap iterative method such as successive overrelaxation (SOR) is applied to further refine the solution to arrive at a desired accuracy. The memory requirements are less than direct LU but more than ILU method. The proposed scheme is shown to yield significant speedups when solving time evolving systems.
Collapse
|
31
|
Beaudoin DL, Roth BJ. How the spatial frequency of polarization influences the induction of reentry in cardiac tissue. J Cardiovasc Electrophysiol 2005; 16:748-52. [PMID: 16050833 DOI: 10.1111/j.1540-8167.2005.40651.x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
UNLABELLED Influences of spatial frequency of polarization. INTRODUCTION The mechanism by which an electric field induces a rotor during cross-field stimulation of cardiac tissue is not entirely known. Different heterogeneous aspects of cardiac tissue have been offered as possible theories, a prominent one being fiber curvature. The polarization produced when an electric field is applied to a sheet of tissue is varied over many spatial frequencies, depending upon the fiber angle. This article compares the effect of high and low spatial frequencies of polarization on reentry induction. METHODS AND RESULTS We incorporate a randomized fiber angle geometry into a two-dimensional active cardiac tissue model with unequal anisotropy ratios already exhibiting smooth, curving fibers. We simulate cross-field stimulation to initiate reentry in the tissue model, and compare the electric field thresholds at different S1-S2 intervals for tissue with randomized fiber angles, tissue with a smooth fiber geometry, and tissue with randomized fiber angles plus smooth, curving fibers. The tissue with both small, random fiber angles and curving fibers has a significantly lower threshold for reentry at certain intervals on the strength-interval curve than for the two cases individually. CONCLUSION Cardiac tissue exhibiting a random fiber geometry in conjunction with a smooth fiber geometry includes high and low spatial frequencies of polarization that may have an effect on the mechanism for reentry at certain S1-S2 intervals. Low spatial frequency regions of hyperpolarization carve out excitable pathways, and high spatial frequency regions provide the large gradient of transmembrane potential required to initiate break excitation.
Collapse
|
32
|
McCulloch AD. Functionally and structurally integrated computational modeling of ventricular physiology. ACTA ACUST UNITED AC 2005; 54:531-9. [PMID: 15760485 DOI: 10.2170/jjphysiol.54.531] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Computational biology is integrative in several ways. Functionally, computational models are valuable for integrating the many interacting processes within biochemical networks and the many interacting physiological subsystems within the cell. Structurally detailed models provide a way of integrating across scales of biological organization from molecule to organism. Data integration across diverse laboratory and clinical measurements is another unique strength of computational biology. We describe examples of all three categories of integration by using recent advances in modeling cardiac excitation-contraction coupling and whole-heart electromechanics in health and disease.
Collapse
Affiliation(s)
- Andrew D McCulloch
- Department of Bioengineering, Whitaker Institute for Biomedical Engineering, University of California San Diego, La Jolla, CA 92093-0412, USA.
| |
Collapse
|
33
|
Rodríguez B, Tice BM, Eason JC, Aguel F, Trayanova N. Cardiac vulnerability to electric shocks during phase 1A of acute global ischemia. Heart Rhythm 2005; 1:695-703. [PMID: 15851241 DOI: 10.1016/j.hrthm.2004.08.018] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/02/2004] [Accepted: 08/24/2004] [Indexed: 11/16/2022]
Abstract
OBJECTIVES The purpose of this study is to characterize the changes in vulnerability to electric shocks during phase 1A of global ischemia in the rabbit ventricles and to determine the mechanisms responsible for these changes. BACKGROUND Mechanisms responsible for the changes in cardiac vulnerability over the course of ischemia phase 1A remain poorly understood. The lack of understanding results from the rapid ischemic change in cardiac electrophysiologic properties, which renders experimental evaluation of vulnerability difficult. METHODS To examine dynamic changes in vulnerability to electric shocks over the course of acute global ischemia phase 1A, this study used a three-dimensional anatomically accurate bidomain model of ischemic rabbit ventricles. Monophasic shocks are applied at various coupling intervals to construct vulnerability grids in normoxia and at various stages of ischemia phase 1A. RESULTS Our simulations demonstrate that 2 to 3 minutes after the onset of ischemia, the upper limit of vulnerability remains at its normoxic value (12.75 V/cm); however, arrhythmias are induced at shorter coupling intervals. As ischemia progresses, the upper limit of vulnerability decreases, reaching 6.4 V/cm in the advanced stage of ischemia phase 1A, and the vulnerable window shifts towards longer coupling intervals. CONCLUSIONS Changes in the upper limit of vulnerability result from an increase in the spatial extent of the shock-end excitation wavefronts and the slower recovery from shock-induced positive polarization. Shifts in the vulnerable window stem from decreases in local repolarization times and the occurrence of postshock conduction failure caused by prolonged postrepolarization refractoriness.
Collapse
|
34
|
Rodríguez B, Li L, Eason JC, Efimov IR, Trayanova NA. Differences between left and right ventricular chamber geometry affect cardiac vulnerability to electric shocks. Circ Res 2005; 97:168-75. [PMID: 15976315 PMCID: PMC2925187 DOI: 10.1161/01.res.0000174429.00987.17] [Citation(s) in RCA: 119] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Although effects of shock strength and waveform on cardiac vulnerability to electric shocks have been extensively documented, the contribution of ventricular anatomy to shock-induced polarization and postshock propagation and thus, to shock outcome, has never been quantified; this is caused by lack of experimental methodology capable of mapping 3-D electrical activity. The goal of this study was to use optical imaging experiments and 3-D bidomain simulations to investigate the role of structural differences between left and right ventricles in vulnerability to electric shocks in rabbit hearts. The ventricles were paced apically, and uniform-field, truncated-exponential, monophasic shocks of reversed polarity were applied over a range of coupling intervals (CIs) in experiment and model. Experiments and simulations revealed that reversing the direction of externally-applied field (RV- or LV- shocks) alters the shape of the vulnerability area (VA), the 2-D grid encompassing episodes of arrhythmia induction. For RV- shocks, VA was nearly rectangular indicating little dependence of postshock arrhythmogenesis on CI. For LV- shocks, the probability of arrhythmia induction was higher for longer than for shorter CIs. The 3-D simulations demonstrated that these effects stem from the fact that reversal of field direction results in relocation of the main postshock excitable area from LV wall (RV- shocks) to septum (LV- shocks). Furthermore, the effect of septal (but not LV) excitable area in postshock propagation was found to strongly depend on preshock state. Knowledge regarding the location of the main postshock excitable area within the 3-D ventricular volume could be important for improving defibrillation efficacy.
Collapse
|
35
|
Li W, Kohl P, Trayanova N. Induction of ventricular arrhythmias following mechanical impact: a simulation study in 3D. J Mol Histol 2005; 35:679-86. [PMID: 15614623 DOI: 10.1007/s10735-004-2666-8] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Commotio cordis, mechanical induction of heart rhythm disturbances, including sudden cardiac death, in the absence of corresponding structural damage, has been reported with increasing frequency in young individuals participating in sporting activities. Recently, the electrophysiological changes during c. cordis have been attributed to mechano-electric feedback, and particularly, to the recruitment of stretch-activated ion channels. The underlying mechanisms, however, by which a mechanical impact results in ventricular fibrillation, remain unknown. This study employs a 3D realistic model of rabbit ventricular geometry and fiber orientation to elucidate the electrophysiological mechanisms involved in arrhythmia induction following acute mechanical stimulation of the heart. Impact effects are modeled through stretch-activated ion channel activation in a 3D region of the ventricles representing the impact profile. Both cation-nonselective and potassium-selective stretch-activated ion channels are recruited upon mechanical impact. The impact is administered at various coupling intervals following pacing at the apex. To aid in the interpretation of results, the effect of mechanical stimulation on single cell action potentials is also examined. The results demonstrate that the region of impact is characterized by different types of cellular responses, including generation of a new action potential, shortening, or lengthening of action potential duration. The impact induces sustained reentry only when (1) a new activation is elicited by mechanical stimulation (caused by activation of cation-nonselective stretch-activated ion channels), and (2) upon return to the original region of impact, this activation does not encounter an extension of action potential duration (prevented by activation of potassium-selective stretch-activated ion channels).
Collapse
Affiliation(s)
- Weihui Li
- Biomedical Engineering, Tulane University, New Orleans, LA 70118, USA
| | | | | |
Collapse
|
36
|
Kuijpers NHL, Keldermann RH, Arts T, Hilbers PAJ. Computer simulations of successful defibrillation in decoupled and non-uniform cardiac tissue. ACTA ACUST UNITED AC 2005; 7 Suppl 2:166-77. [PMID: 16102514 DOI: 10.1016/j.eupc.2005.03.021] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2004] [Revised: 02/03/2005] [Accepted: 05/03/2005] [Indexed: 11/19/2022]
Abstract
Abstract
Aim
The aim of the present study is to investigate the origin and effect of virtual electrode polarization in uniform, decoupled and non-uniform cardiac tissue during field stimulation.
Methods
A discrete bidomain model with active membrane behaviour was used to simulate normal cardiac tissue as well as cardiac tissue that is decoupled due to fibrosis and gap junction remodelling. Various uniform and non-uniform electric fields were applied to the external domain of uniform, decoupled and non-uniform resting cardiac tissue as well as cardiac tissue in which spiral waves were induced.
Results
Field stimulation applied on non-uniform tissue results in more virtual electrodes compared with uniform tissue. The spiral waves were terminated in decoupled tissue, but not in uniform, homogeneous tissue. By gradually increasing local differences in intracellular conductivities, the amount and spread of virtual electrodes increased and the spiral waves were terminated.
Conclusion
Fast depolarization of the tissue after field stimulation may be explained by intracellular decoupling and spatial heterogeneity present in normal and pathological cardiac tissue. We demonstrated that termination of spiral waves by means of field stimulation can be achieved when the tissue is modelled as a non-uniform, anisotropic bidomain with active membrane behaviour.
Collapse
Affiliation(s)
- N H L Kuijpers
- Department of Biomedical Engineering, Technische Universiteit Eindhoven, The Netherlands.
| | | | | | | |
Collapse
|