A Cell-Based Framework for Numerical Modeling of Electrical Conduction in Cardiac Tissue
Open Access
- 10 October 2017
- journal article
- research article
- Published by Frontiers Media SA in Frontiers in Physics
Abstract
In this paper, we study a mathematical model of cardiac tissue based on explicit representation of individual cells. In this EMI model, the extracellular (E) space, the cell membrane (M) and the intracellular (I) space are represented as separate geometrical domains. This representation introduces modelling flexibility needed for detailed representation of the properties of cardiac cells including their membrane. In particular, we will show that the model allows ion channels to be non-uniformly distributed along the membrane of the cell. Such features are difficult to include in classical homogenized models like the monodomain and bidomain models frequently used in computational analyses of cardiac electrophysiology. The EMI model is solved using a finite difference method (FDM) and two variants of the finite element method (FEM). We compare the three schemes numerically, reporting on CPU-efforts and convergence rates. Finally, we illustrate the distinctive capabilities of the EMI model compared to classical models by simulating monolayers of cardiac cells with heterogeneous distributions of ionic channels along the cell membrane. Because of the detailed representation of every cell, the computational problems that result from using the EMI model are much larger than for the classical homogenized models, and thus represent a computational challenge. However, our numerical simulations indicate that the FDM scheme is optimal in the sense that the computational complexity increases proportionally to the number of cardiac cells in the model. Moreover, we present simulations, based on systems of equations involving ~ 117 million unknowns, representing up to ~ 16000 cells. We conclude that collections of cardiac cells can be simulated using the EMI model, and that the EMI model enable greater modeling flexibility than the classical monodomain and bidomain models.This publication has 78 references indexed in Scilit:
- Instabilities of the resting state in a mathematical model of calcium handling in cardiac myocytesMathematical Biosciences, 2012
- Note on a Possible Proarrhythmic Property of Antiarrhythmic Drugs Aimed at Improving Gap-Junction CouplingBiophysical Journal, 2012
- A Computational Model to Predict the Effects of Class I Anti-Arrhythmic Drugs on Ventricular RhythmsScience Translational Medicine, 2011
- So Little Source, So Much Sink: Requirements for Afterdepolarizations to Propagate in TissueBiophysical Journal, 2010
- A novel computational model of the human ventricular action potential and Ca transientJournal of Molecular and Cellular Cardiology, 2010
- Incorporating Histology into a 3D Microscopic Computer Model of Myocardium to Study Propagation at a Cellular LevelAnnals of Biomedical Engineering, 2010
- K+current changes account for the rate dependence of the action potential in the human atrial myocyteAmerican Journal of Physiology-Heart and Circulatory Physiology, 2009
- Loss of T-tubules and other changes to surface topography in ventricular myocytes from failing human and rat heartProceedings of the National Academy of Sciences of the United States of America, 2009
- Towards predictive modelling of the electrophysiology of the heartExperimental Physiology, 2009
- Computational biology in the study of cardiac ion channels and cell electrophysiologyQuarterly Reviews of Biophysics, 2006