1
|
An Automata-Based Cardiac Electrophysiology Simulator to Assess Arrhythmia Inducibility. MATHEMATICS 2022. [DOI: 10.3390/math10081293] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Personalized cardiac electrophysiology simulations have demonstrated great potential to study cardiac arrhythmias and help in therapy planning of radio-frequency ablation. Its application to analyze vulnerability to ventricular tachycardia and sudden cardiac death in infarcted patients has been recently explored. However, the detailed multi-scale biophysical simulations used in these studies are very demanding in terms of memory and computational resources, which prevents their clinical translation. In this work, we present a fast phenomenological system based on cellular automata (CA) to simulate personalized cardiac electrophysiology. The system is trained on biophysical simulations to reproduce cellular and tissue dynamics in healthy and pathological conditions, including action potential restitution, conduction velocity restitution and cell safety factor. We show that a full ventricular simulation can be performed in the order of seconds, emulate the results of a biophysical simulation and reproduce a patient’s ventricular tachycardia in a model that includes a heterogeneous scar region. The system could be used to study the risk of arrhythmia in infarcted patients for a large number of scenarios.
Collapse
|
2
|
Bakir AA, Al Abed A, Lovell NH, Dokos S. Multiphysics computational modelling of the cardiac ventricles. IEEE Rev Biomed Eng 2021; 15:309-324. [PMID: 34185649 DOI: 10.1109/rbme.2021.3093042] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Development of cardiac multiphysics models has progressed significantly over the decades and simulations combining multiple physics interactions have become increasingly common. In this review, we summarise the progress in this field focusing on various approaches of integrating ventricular structures. electrophysiological properties, myocardial mechanics, as well as incorporating blood hemodynamics and the circulatory system. Common coupling approaches are discussed and compared, including the advantages and shortcomings of each. Currently used strategies for patient-specific implementations are highlighted and potential future improvements considered.
Collapse
|
3
|
Kojic M, Milosevic M, Simic V, Geroski V, Ziemys A, Filipovic N, Ferrari M. Smeared multiscale finite element model for electrophysiology and ionic transport in biological tissue. Comput Biol Med 2019; 108:288-304. [PMID: 31015049 DOI: 10.1016/j.compbiomed.2019.03.023] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2018] [Revised: 03/22/2019] [Accepted: 03/23/2019] [Indexed: 10/27/2022]
Abstract
Basic functions of living organisms are governed by the nervous system through bidirectional signals transmitted from the brain to neural networks. These signals are similar to electrical waves. In electrophysiology the goal is to study the electrical properties of biological cells and tissues, and the transmission of signals. From a physics perspective, there exists a field of electrical potential within the living body, the nervous system, extracellular space and cells. Electrophysiological problems can be investigated experimentally and also theoretically by developing appropriate mathematical or computational models. Due to the enormous complexity of biological systems, it would be almost impossible to establish a detailed computational model of the electrical field, even for only a single organ (e.g. heart), including the entirety of cells comprising the neural network. In order to make computational models feasible for practical applications, we here introduce the concept of smeared fields, which represents a generalization of the previously formulated multiscale smeared methodology for mass transport in blood vessels, lymph, and tissue. We demonstrate the accuracy of the smeared finite element computational models for the electric field in numerical examples. The electrical field is further coupled with ionic mass transport within tissue composed of interstitial spaces extracellularly and by cytoplasm and organelles intracellularly. The proposed methodology, which couples electrophysiology and molecular ionic transport, is applicable to a variety of biological systems.
Collapse
Affiliation(s)
- M Kojic
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7-117, Houston, TX, 77030, USA; Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia; Serbian Academy of Sciences and Arts, Knez Mihailova 35, 11000, Belgrade, Serbia.
| | - M Milosevic
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia; Belgrade Metropolitan University, Tadeuša Košćuška 63, 11000, Belgrade, Serbia
| | - V Simic
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
| | - V Geroski
- Bioengineering Research and Development Center BioIRC Kragujevac, Prvoslava Stojanovica 6, 3400, Kragujevac, Serbia
| | - A Ziemys
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7-117, Houston, TX, 77030, USA
| | - N Filipovic
- University of Kragujevac, Faculty for Engineering Sciences, Sestre Janic 6, 34000, Kragujevac, Serbia
| | - M Ferrari
- Houston Methodist Research Institute, The Department of Nanomedicine, 6670 Bertner Ave., R7-117, Houston, TX, 77030, USA
| |
Collapse
|
4
|
Yang PC, Purawat S, Ieong PU, Jeng MT, DeMarco KR, Vorobyov I, McCulloch AD, Altintas I, Amaro RE, Clancy CE. A demonstration of modularity, reuse, reproducibility, portability and scalability for modeling and simulation of cardiac electrophysiology using Kepler Workflows. PLoS Comput Biol 2019; 15:e1006856. [PMID: 30849072 PMCID: PMC6426265 DOI: 10.1371/journal.pcbi.1006856] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2018] [Revised: 03/20/2019] [Accepted: 02/08/2019] [Indexed: 01/18/2023] Open
Abstract
Multi-scale computational modeling is a major branch of computational biology as evidenced by the US federal interagency Multi-Scale Modeling Consortium and major international projects. It invariably involves specific and detailed sequences of data analysis and simulation, often with multiple tools and datasets, and the community recognizes improved modularity, reuse, reproducibility, portability and scalability as critical unmet needs in this area. Scientific workflows are a well-recognized strategy for addressing these needs in scientific computing. While there are good examples if the use of scientific workflows in bioinformatics, medical informatics, biomedical imaging and data analysis, there are fewer examples in multi-scale computational modeling in general and cardiac electrophysiology in particular. Cardiac electrophysiology simulation is a mature area of multi-scale computational biology that serves as an excellent use case for developing and testing new scientific workflows. In this article, we develop, describe and test a computational workflow that serves as a proof of concept of a platform for the robust integration and implementation of a reusable and reproducible multi-scale cardiac cell and tissue model that is expandable, modular and portable. The workflow described leverages Python and Kepler-Python actor for plotting and pre/post-processing. During all stages of the workflow design, we rely on freely available open-source tools, to make our workflow freely usable by scientists. We present a computational workflow as a proof of concept for integration and implementation of a reusable and reproducible cardiac multi-scale electrophysiology model that is expandable, modular and portable. This framework enables scientists to create intuitive, user-friendly and flexible end-to-end automated scientific workflows using a graphical user interface. Kepler is an advanced open-source platform that supports multiple models of computation. The underlying workflow engine handles scalability, provenance, reproducibility aspects of the code, performs orchestration of data flow, and automates execution on heterogeneous computing resources. One of the main advantages of workflow utilization is the integration of code written in multiple languages Standardization occurs at the interfaces of the workflow elements and allows for general applications and easy comparison and integration of code from different research groups or even multiple programmers coding in different languages for various purposes from the same group. A workflow driven problem-solving approach enables domain scientists to focus on resolving the core science questions, and delegates the computational and process management burden to the underlying Workflow. The workflow driven approach allows scaling the computational experiment with distributed data-parallel execution on multiple computing platforms, such as, HPC resources, GPU clusters, Cloud etc. The workflow framework tracks software version information along with hardware information to allow users an opportunity to trace any variation in workflow outcome to the system configurations.
Collapse
Affiliation(s)
- Pei-Chi Yang
- Department of Physiology and Membrane Biology, Department of Pharmacology, School of Medicine, University of California Davis, Davis, California, United States of America
| | - Shweta Purawat
- San Diego Supercomputer Center (SDSC), University of California, San Diego, La Jolla, California, United States of America
| | - Pek U. Ieong
- Department of Chemistry and Biochemistry, National Biomedical Computation Resource, Drug Design Data Resource (D3R), University of California San Diego, La Jolla, California, United States of America
| | - Mao-Tsuen Jeng
- Department of Physiology and Membrane Biology, Department of Pharmacology, School of Medicine, University of California Davis, Davis, California, United States of America
| | - Kevin R. DeMarco
- Department of Physiology and Membrane Biology, Department of Pharmacology, School of Medicine, University of California Davis, Davis, California, United States of America
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, Department of Pharmacology, School of Medicine, University of California Davis, Davis, California, United States of America
| | - Andrew D. McCulloch
- Departments of Bioengineering and Medicine, University of California, San Diego, La Jolla, California, United States of America
| | - Ilkay Altintas
- San Diego Supercomputer Center (SDSC), University of California, San Diego, La Jolla, California, United States of America
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, National Biomedical Computation Resource, Drug Design Data Resource (D3R), University of California San Diego, La Jolla, California, United States of America
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, Department of Pharmacology, School of Medicine, University of California Davis, Davis, California, United States of America
- * E-mail:
| |
Collapse
|
5
|
Clancy CE, Chen-Izu Y, Bers DM, Belardinelli L, Boyden PA, Csernoch L, Despa S, Fermini B, Hool LC, Izu L, Kass RS, Lederer WJ, Louch WE, Maack C, Matiazzi A, Qu Z, Rajamani S, Rippinger CM, Sejersted OM, O'Rourke B, Weiss JN, Varró A, Zaza A. Deranged sodium to sudden death. J Physiol 2015; 593:1331-45. [PMID: 25772289 DOI: 10.1113/jphysiol.2014.281204] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2014] [Accepted: 10/14/2014] [Indexed: 12/19/2022] Open
Abstract
In February 2014, a group of scientists convened as part of the University of California Davis Cardiovascular Symposium to bring together experimental and mathematical modelling perspectives and discuss points of consensus and controversy on the topic of sodium in the heart. This paper summarizes the topics of presentation and discussion from the symposium, with a focus on the role of aberrant sodium channels and abnormal sodium homeostasis in cardiac arrhythmias and pharmacotherapy from the subcellular scale to the whole heart. Two following papers focus on Na(+) channel structure, function and regulation, and Na(+)/Ca(2+) exchange and Na(+)/K(+) ATPase. The UC Davis Cardiovascular Symposium is a biannual event that aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The focus on Na(+) in the 2014 symposium stemmed from the multitude of recent studies that point to the importance of maintaining Na(+) homeostasis in the heart, as disruption of homeostatic processes are increasingly identified in cardiac disease states. Understanding how disruption in cardiac Na(+)-based processes leads to derangement in multiple cardiac components at the level of the cell and to then connect these perturbations to emergent behaviour in the heart to cause disease is a critical area of research. The ubiquity of disruption of Na(+) channels and Na(+) homeostasis in cardiac disorders of excitability and mechanics emphasizes the importance of a fundamental understanding of the associated mechanisms and disease processes to ultimately reveal new targets for human therapy.
Collapse
Affiliation(s)
- Colleen E Clancy
- Department of Pharmacology, University of California, Davis, Genome Building Rm 3503, Davis, CA, 95616-8636, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
6
|
Cooper J, Spiteri RJ, Mirams GR. Cellular cardiac electrophysiology modeling with Chaste and CellML. Front Physiol 2015; 5:511. [PMID: 25610400 PMCID: PMC4285015 DOI: 10.3389/fphys.2014.00511] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2014] [Accepted: 12/09/2014] [Indexed: 11/18/2022] Open
Abstract
Chaste is an open-source C++ library for computational biology that has well-developed cardiac electrophysiology tissue simulation support. In this paper, we introduce the features available for performing cardiac electrophysiology action potential simulations using a wide range of models from the Physiome repository. The mathematics of the models are described in CellML, with units for all quantities. The primary idea is that the model is defined in one place (the CellML file), and all model code is auto-generated at compile or run time; it never has to be manually edited. We use ontological annotation to identify model variables describing certain biological quantities (membrane voltage, capacitance, etc.) to allow us to import any relevant CellML models into the Chaste framework in consistent units and to interact with them via consistent interfaces. This approach provides a great deal of flexibility for analysing different models of the same system. Chaste provides a wide choice of numerical methods for solving the ordinary differential equations that describe the models. Fixed-timestep explicit and implicit solvers are provided, as discussed in previous work. Here we introduce the Rush–Larsen and Generalized Rush–Larsen integration techniques, made available via symbolic manipulation of the model equations, which are automatically rearranged into the forms required by these approaches. We have also integrated the CVODE solvers, a ‘gold standard’ for stiff systems, and we have developed support for symbolic computation of the Jacobian matrix, yielding further increases in the performance and accuracy of CVODE. We discuss some of the technical details of this work and compare the performance of the available numerical methods. Finally, we discuss how this is generalized in our functional curation framework, which uses a domain-specific language for defining complex experiments as a basis for comparison of model behavior.
Collapse
Affiliation(s)
- Jonathan Cooper
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| | - Raymond J Spiteri
- Numerical Simulation Research Lab, Department of Computer Science, University of Saskatchewan Saskatoon, SK, Canada
| | - Gary R Mirams
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| |
Collapse
|
7
|
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
|
8
|
Krishnamoorthi S, Sarkar M, Klug WS. Numerical quadrature and operator splitting in finite element methods for cardiac electrophysiology. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2013; 29:1243-66. [PMID: 23873868 PMCID: PMC4519349 DOI: 10.1002/cnm.2573] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2012] [Revised: 06/06/2013] [Accepted: 06/07/2013] [Indexed: 05/05/2023]
Abstract
We study the numerical accuracy and computational efficiency of alternative formulations of the finite element solution procedure for the monodomain equations of cardiac electrophysiology, focusing on the interaction of spatial quadrature implementations with operator splitting and examining both nodal and Gauss quadrature methods and implementations that mix nodal storage of state variables with Gauss quadrature. We evaluate the performance of all possible combinations of 'lumped' approximations of consistent capacitance and mass matrices. Most generally, we find that quadrature schemes and lumped approximations that produce decoupled nodal ionic equations allow for the greatest computational efficiency, this being afforded through the use of asynchronous adaptive time-stepping of the ionic state variable ODEs. We identify two lumped approximation schemes that exhibit superior accuracy, rivaling that of the most expensive variationally consistent implementations. Finally, we illustrate some of the physiological consequences of discretization error in electrophysiological simulation relevant to cardiac arrhythmia and fibrillation. These results suggest caution with the use of semi-automated free-form tetrahedral and hexahedral meshing algorithms available in most commercially available meshing software, which produce nonuniform meshes having a large distribution of element sizes.
Collapse
|
9
|
Roberts BN, Yang PC, Behrens SB, Moreno JD, Clancy CE. Computational approaches to understand cardiac electrophysiology and arrhythmias. Am J Physiol Heart Circ Physiol 2012; 303:H766-83. [PMID: 22886409 DOI: 10.1152/ajpheart.01081.2011] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.
Collapse
Affiliation(s)
- Byron N Roberts
- Tri-Institutional MD-PhD Program, Physiology, Biophysics and Systems Biology Graduate Program, Weill Cornell Medical College/The Rockefeller University/Sloan-Kettering Cancer Institute, Weill Medical College of Cornell University, New York, New York, USA
| | | | | | | | | |
Collapse
|
10
|
Burton BM, Tate JD, Erem B, Swenson DJ, Wang DF, Steffen M, Brooks DH, van Dam PM, Macleod RS. A toolkit for forward/inverse problems in electrocardiography within the SCIRun problem solving environment. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2012; 2011:267-70. [PMID: 22254301 DOI: 10.1109/iembs.2011.6090052] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Computational modeling in electrocardiography often requires the examination of cardiac forward and inverse problems in order to non-invasively analyze physiological events that are otherwise inaccessible or unethical to explore. The study of these models can be performed in the open-source SCIRun problem solving environment developed at the Center for Integrative Biomedical Computing (CIBC). A new toolkit within SCIRun provides researchers with essential frameworks for constructing and manipulating electrocardiographic forward and inverse models in a highly efficient and interactive way. The toolkit contains sample networks, tutorials and documentation which direct users through SCIRun-specific approaches in the assembly and execution of these specific problems.
Collapse
Affiliation(s)
- Brett M Burton
- Department of Bioengineering, University of Utah, Salt Lake City, UT 84112, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|
11
|
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
|
12
|
Corrias A, Pathmanathan P, Gavaghan DJ, Buist ML. Modelling tissue electrophysiology with multiple cell types: applications of the extended bidomain framework. Integr Biol (Camb) 2012; 4:192-201. [PMID: 22222297 DOI: 10.1039/c2ib00100d] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The bidomain framework has been extensively used to model tissue electrophysiology in a variety of applications. One limitation of the bidomain model is that it describes the activity of only one cell type interacting with the extracellular space. If more than one cell type contributes to the tissue electrophysiology, then the bidomain model is not sufficient. Recently, evidence has suggested that this is the case for at least two important applications: cardiac and gastrointestinal tissue electrophysiology. In the heart, fibroblasts ubiquitously interact with myocytes and are believed to play an important role in the organ electrophysiology. Along the GI tract, interstitial cells of Cajal (ICC) generate electrical waves that are passed on to surrounding smooth muscle cells (SMC), which are interconnected with the ICC and with each other. Because of the contribution of more than one cell type to the overall organ electrophysiology, investigators in different fields have independently proposed similar extensions of the bidomain model to incorporate multiple cell types and tested it on simplified geometries. In this paper, we provide a general derivation of such an extended bidomain framework applicable to any tissue and provide a generic and efficient implementation applicable to any geometry. Proof-of-concept results of tissue electrophysiology on realistic 3D organ geometries using the extended bidomain framework are presented for the heart and the stomach.
Collapse
Affiliation(s)
- Alberto Corrias
- National University of Singapore - Bioengineering, 9 Engineering Drive 1 Block EA #03-12, Singapore 117576, Singapore.
| | | | | | | |
Collapse
|
13
|
O'Hara T, Virág L, Varró A, Rudy Y. Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput Biol 2011; 7:e1002061. [PMID: 21637795 PMCID: PMC3102752 DOI: 10.1371/journal.pcbi.1002061] [Citation(s) in RCA: 701] [Impact Index Per Article: 53.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2010] [Accepted: 04/05/2011] [Indexed: 11/19/2022] Open
Abstract
Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca(2+) versus voltage dependent inactivation of L-type Ca(2+) current (I(CaL)); kinetics for the transient outward, rapid delayed rectifier (I(Kr)), Na(+)/Ca(2+) exchange (I(NaCa)), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca(2+) (including peak and decay) and intracellular sodium ([Na(+)](i)) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by I(Kr) block during slow pacing, and AP and Ca(2+) alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca(2+)/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca(2+) cycling. I(NaCa) linked Ca(2+) alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na(+)](i), via its modulation of the electrogenic Na(+)/K(+) ATPase current. At fast pacing rates, late Na(+) current and I(CaL) were also contributors. APD shortening during restitution was primarily dependent on reduced late Na(+) and I(CaL) currents due to inactivation at short diastolic intervals, with additional contribution from elevated I(Kr) due to incomplete deactivation.
Collapse
Affiliation(s)
- Thomas O'Hara
- Cardiac Bioelectricity and Arrhythmia Center, Department of Biomedical
Engineering, Washington University in St. Louis, St. Louis, Missouri, United
States of America
| | - László Virág
- Department of Pharmacology and Pharmacotherapy, University of Szeged,
Szeged, Hungary
| | - András Varró
- Department of Pharmacology and Pharmacotherapy, University of Szeged,
Szeged, Hungary
- Division of Cardiovascular Pharmacology, Hungarian Academy of Sciences,
Szeged, Hungary
| | - Yoram Rudy
- Cardiac Bioelectricity and Arrhythmia Center, Department of Biomedical
Engineering, Washington University in St. Louis, St. Louis, Missouri, United
States of America
- * E-mail:
| |
Collapse
|
14
|
Corrias A, Giles W, Rodriguez B. Ionic mechanisms of electrophysiological properties and repolarization abnormalities in rabbit Purkinje fibers. Am J Physiol Heart Circ Physiol 2011; 300:H1806-13. [PMID: 21335469 DOI: 10.1152/ajpheart.01170.2010] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Purkinje cells play an important role in drug-induced arrhythmogenesis and are widely used in preclinical drug safety assessments. Repolarization abnormalities such as action potential (AP) prolongation and early afterdeploarizations (EAD) are often observed in vitro upon pharmacological interventions. However, because drugs do not act on only one defined target, it is often difficult to fully explain the mechanisms of action and their potential arrhythmogenicity. Computational models, when appropriately detailed and validated, can be used to gain mechanistic insights into the mechanisms of action of certain drugs. Nevertheless, no model of Purkinje electrophysiology that is able to reproduce characteristic Purkinje responses to drug-induced changes in ionic current conductances such as AP prolongation and EAD generation currently exists. In this study, a novel biophysically detailed model of rabbit Purkinje electrophysiology was developed by integration of data from voltage-clamp and AP experimental recordings. Upon validation, we demonstrate that the model reproduces many key electrophysiological properties of rabbit Purkinje cells. These include: AP morphology and duration, both input resistance and rate dependence properties as well as response to hyperkalemia. Pharmacological interventions such as inward rectifier K(+) current and rapid delayed rectifier K(+) current block as well as late Na(+) current increase result in significant AP changes. However, enhanced L-type Ca(2+) current (i(CaL)) dominates in EAD genesis in Purkinje fibers. In addition, i(CaL) inactivation dynamics and intercellular coupling in tissue strongly modulate EAD formation. We conclude that EAD generation in Purkinje cells is mediated by an increase in i(CaL) and modulated by its inactivation kinetics.
Collapse
|
15
|
Dada JO, Mendes P. Multi-scale modelling and simulation in systems biology. Integr Biol (Camb) 2011; 3:86-96. [DOI: 10.1039/c0ib00075b] [Citation(s) in RCA: 129] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Abramson D, Bernabeu MO, Bethwaite B, Burrage K, Corrias A, Enticott C, Garic S, Gavaghan D, Peachey T, Pitt-Francis J, Pueyo E, Rodriguez B, Sher A, Tan J. High-throughput cardiac science on the Grid. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2010; 368:3907-3923. [PMID: 20643684 DOI: 10.1098/rsta.2010.0170] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Cardiac electrophysiology is a mature discipline, with the first model of a cardiac cell action potential having been developed in 1962. Current models range from single ion channels, through very complex models of individual cardiac cells, to geometrically and anatomically detailed models of the electrical activity in whole ventricles. A critical issue for model developers is how to choose parameters that allow the model to faithfully reproduce observed physiological effects without over-fitting. In this paper, we discuss the use of a parametric modelling toolkit, called Nimrod, that makes it possible both to explore model behaviour as parameters are changed and also to tune parameters by optimizing model output. Importantly, Nimrod leverages computers on the Grid, accelerating experiments by using available high-performance platforms. We illustrate the use of Nimrod with two case studies, one at the cardiac tissue level and one at the cellular level.
Collapse
Affiliation(s)
- David Abramson
- Faculty of Information Technology, Monash University, Caulfield East, Melbourne, Victoria 3145, Australia.
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
17
|
A numerical guide to the solution of the bi-domain equations of cardiac electrophysiology. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2010; 102:136-55. [PMID: 20553747 DOI: 10.1016/j.pbiomolbio.2010.05.006] [Citation(s) in RCA: 62] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2010] [Accepted: 05/19/2010] [Indexed: 11/22/2022]
Abstract
Simulation of cardiac electrical activity using the bi-domain equations can be a massively computationally demanding problem. This study provides a comprehensive guide to numerical bi-domain modelling. Each component of bi-domain simulations--discretization, ODE-solution, linear system solution, and parallelization--is discussed, and previously-used methods are reviewed, new methods are proposed, and issues which cause particular difficulty are highlighted. Particular attention is paid to the choice of stimulus currents, compatibility conditions for the equations, the solution of singular linear systems, and convergence of the numerical scheme.
Collapse
|
18
|
Heidenreich EA, Ferrero JM, Doblaré M, Rodríguez JF. Adaptive macro finite elements for the numerical solution of monodomain equations in cardiac electrophysiology. Ann Biomed Eng 2010; 38:2331-45. [PMID: 20238165 DOI: 10.1007/s10439-010-9997-2] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2009] [Accepted: 03/03/2010] [Indexed: 11/26/2022]
Abstract
Many problems in biology and engineering are governed by anisotropic reaction-diffusion equations with a very rapidly varying reaction term. This usually implies the use of very fine meshes and small time steps in order to accurately capture the propagating wave while avoiding the appearance of spurious oscillations in the wave front. This work develops a family of macro finite elements amenable for solving anisotropic reaction-diffusion equations with stiff reactive terms. The developed elements are incorporated on a semi-implicit algorithm based on operator splitting that includes adaptive time stepping for handling the stiff reactive term. A linear system is solved on each time step to update the transmembrane potential, whereas the remaining ordinary differential equations are solved uncoupled. The method allows solving the linear system on a coarser mesh thanks to the static condensation of the internal degrees of freedom (DOF) of the macroelements while maintaining the accuracy of the finer mesh. The method and algorithm have been implemented in parallel. The accuracy of the method has been tested on two- and three-dimensional examples demonstrating excellent behavior when compared to standard linear elements. The better performance and scalability of different macro finite elements against standard finite elements have been demonstrated in the simulation of a human heart and a heterogeneous two-dimensional problem with reentrant activity. Results have shown a reduction of up to four times in computational cost for the macro finite elements with respect to equivalent (same number of DOF) standard linear finite elements as well as good scalability properties.
Collapse
Affiliation(s)
- Elvio A Heidenreich
- Group of Structural Mechanics and Materials Modeling, Aragón Institute of Engineering Research (I3A), Universidad de Zaragoza, Zaragoza, Spain.
| | | | | | | |
Collapse
|
19
|
|
20
|
Gavaghan D, Coveney PV, Kohl P. The virtual physiological human: tools and applications I. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2009; 367:1817-1821. [PMID: 19380313 DOI: 10.1098/rsta.2009.0070] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
|