High-throughput functional curation of cellular electrophysiology models.

Effective reuse of a quantitative mathematical model requires not just access to curated versions of the model equations, but also an understanding of the functional capabilities of the model, and the advisable scope of its application. To enable this "functional curation" we have developed a simulation environment that provides high-throughput evaluation of a mathematical model's functional response to an arbitrary user-defined protocol, and optionally compares the results against experimental data. In this study we demonstrate the efficacy of this simulation environment on 31 cardiac electrophysiology cell models using two test cases. The S1-S2 response is evaluated to characterise the models' restitution curves, and their L-type calcium channel current-voltage curves are evaluated. The significant variation in the response of these models, even when the models represent the same species and temperature, demonstrates the importance of knowing the functional characteristics of a model prior to its reuse. We also discuss the wider implications for this approach, in improving the selection of models for reuse, enabling the identification of models that exhibit particular experimentally observed phenomena, and making the incremental development of models more robust.

[1]  C Antzelevitch,et al.  A subpopulation of cells with unique electrophysiological properties in the deep subepicardium of the canine ventricle. The M cell. , 1991, Circulation research.

[2]  Satoshi Matsuoka,et al.  Role of individual ionic current systems in ventricular cells hypothesized by a model study. , 2003, The Japanese journal of physiology.

[3]  H Zhang,et al.  Models of cardiac tissue electrophysiology: progress, challenges and open questions. , 2011, Progress in biophysics and molecular biology.

[4]  Joseph L Greenstein,et al.  Cardiac myocytes and local signaling in nano-domains. , 2011, Progress in biophysics and molecular biology.

[5]  Michal Pásek,et al.  A model of the guinea-pig ventricular cardiac myocyte incorporating a transverse-axial tubular system. , 2008, Progress in biophysics and molecular biology.

[6]  Yoram Rudy,et al.  Properties and ionic mechanisms of action potential adaptation, restitution, and accommodation in canine epicardium. , 2009, American journal of physiology. Heart and circulatory physiology.

[7]  Toni Giorgino,et al.  The Virtual Physiological Human ToolKit , 2010, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[8]  Donald M Bers,et al.  A mathematical treatment of integrated Ca dynamics within the ventricular myocyte. , 2004, Biophysical journal.

[9]  Yoram Rudy,et al.  Regulation of Ca2+ and electrical alternans in cardiac myocytes: role of CAMKII and repolarizing currents. , 2007, American journal of physiology. Heart and circulatory physiology.

[10]  T. Lincoln,et al.  Myosin phosphatase regulatory pathways: different functions or redundant functions? , 2007, Circulation research.

[11]  Frank Bergmann A Simple Nested Simulation for SED-ML , 2010 .

[12]  J. Nolasco,et al.  A graphic method for the study of alternation in cardiac action potentials. , 1968, Journal of applied physiology.

[13]  Rodrigo Weber dos Santos,et al.  CellML and associated tools and techniques , 2008, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[14]  David P Nickerson,et al.  A physiome standards-based model publication paradigm , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[15]  John Glauert,et al.  SISAL: streams and iteration in a single assignment language. Language reference manual, Version 1. 2. Revision 1 , 1985 .

[16]  James B. Bassingthwaighte,et al.  Strategies for the Physiome Project , 2000, Annals of Biomedical Engineering.

[17]  Sarah M. Keating,et al.  Evolving a lingua franca and associated software infrastructure for computational systems biology: the Systems Biology Markup Language (SBML) project. , 2004, Systems biology.

[18]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[19]  John W. Clark,et al.  Modeling CICR in rat ventricular myocytes: voltage clamp studies , 2010, Theoretical Biology and Medical Modelling.

[20]  Elizabeth M Cherry,et al.  A tale of two dogs: analyzing two models of canine ventricular electrophysiology. , 2007, American journal of physiology. Heart and circulatory physiology.

[21]  Denis Noble,et al.  Modulatory effect of calmodulin-dependent kinase II (CaMKII) on sarcoplasmic reticulum Ca2+ handling and interval–force relations: a modelling study , 2006, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  Y Rudy,et al.  Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study. , 2000, Biophysical journal.

[23]  Alexander V. Panfilov,et al.  human ventricular tissue , 2007 .

[24]  Peter J. Hunter,et al.  The CellML Model Repository , 2008, Bioinform..

[25]  Leon Espinosa L'échange Na+-Ca2+ dans l'hypertrophie ventriculaire d'altitude chez le rat : étude électrophysiologique et utilisation du modèle "oxsoft heart" , 1998 .

[26]  D. Noble,et al.  Excitation-contraction coupling and extracellular calcium transients in rabbit atrium: reconstruction of basic cellular mechanisms , 1987, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[27]  Alberto Corrias,et al.  Ionic mechanisms of electrophysiological properties and repolarization abnormalities in rabbit Purkinje fibers. , 2011, American journal of physiology. Heart and circulatory physiology.

[28]  R. Winslow,et al.  Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, II: model studies. , 1999, Circulation research.

[29]  Alexander G. Fletcher,et al.  Chaste: A test-driven approach to software development for biological modelling , 2009, Comput. Phys. Commun..

[30]  Melanie I. Stefan,et al.  BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models , 2010, BMC Systems Biology.

[31]  K. T. ten Tusscher,et al.  Alternans and spiral breakup in a human ventricular tissue model. , 2006, American journal of physiology. Heart and circulatory physiology.

[32]  D. Noble,et al.  A model of the single atrial cell: relation between calcium current and calcium release , 1990, Proceedings of the Royal Society of London. B. Biological Sciences.

[33]  D. Bers,et al.  A novel computational model of the human ventricular action potential and Ca transient. , 2010, Journal of Molecular and Cellular Cardiology.

[34]  A V Holden,et al.  Towards a computational reconstruction of the electrodynamics of premature and full term human labour. , 2011, Progress in biophysics and molecular biology.

[35]  D Noble,et al.  A meta‐analysis of cardiac electrophysiology computational models , 2009, Experimental physiology.

[36]  E. Rowland,et al.  Dispersion of monophasic action potential duration: demonstrable in humans after premature ventricular extrastimulation but not in steady state. , 1992, Journal of the American College of Cardiology.

[37]  Steve McKeever,et al.  A model‐driven approach to automatic conversion of physical units , 2008, Softw. Pract. Exp..

[38]  J. Restrepo,et al.  A rabbit ventricular action potential model replicating cardiac dynamics at rapid heart rates. , 2007, Biophysical journal.

[39]  W. Giles,et al.  A mathematical model of action potential heterogeneity in adult rat left ventricular myocytes. , 2001, Biophysical journal.

[40]  D. Noble,et al.  The Role of Sodium ‐ Calcium Exchange during the Cardiac Action Potential a , 1991, Annals of the New York Academy of Sciences.

[41]  Henggui Zhang,et al.  Optimal velocity and safety of discontinuous conduction through the heterogeneous Purkinje-ventricular junction. , 2009, Biophysical journal.

[42]  Denis Noble,et al.  Cellular Open Resource (COR): a Public CellML Based Environment for Modeling Biological Function , 2003, Int. J. Bifurc. Chaos.

[43]  Antonis A Armoundas,et al.  Mechanisms of Abnormal Calcium Homeostasis in Mutations Responsible for Catecholaminergic Polymorphic Ventricular Tachycardia , 2007, Circulation research.

[44]  Kenneth E. Iverson,et al.  A programming language , 1899, AIEE-IRE '62 (Spring).

[45]  D. Noble,et al.  A model for human ventricular tissue. , 2004, American journal of physiology. Heart and circulatory physiology.

[46]  Nicolas Le Novère,et al.  SED-ML - An XML Format for the Implementation of the MIASE Guidelines , 2008, CMSB.

[47]  Edmund J. Crampin,et al.  Computational biology of cardiac myocytes: proposed standards for the physiome , 2007, Journal of Experimental Biology.

[48]  R. Winslow,et al.  A computational model of the human left-ventricular epicardial myocyte. , 2004, Biophysical journal.

[49]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[50]  Peter J. Hunter,et al.  CellML metadata standards, associated tools and repositories , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[51]  C. Luo,et al.  A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. , 1991, Circulation research.

[52]  Nicolas Le Novère,et al.  A simple scheme for annotating SBML with references to controlled vocabularies and database entries , 2022 .

[53]  Denis Noble,et al.  Contributions of HERG K+ current to repolarization of the human ventricular action potential. , 2008, Progress in biophysics and molecular biology.

[54]  Norman W. Paton,et al.  SBRML: a markup language for associating systems biology data with models , 2010, Bioinform..

[55]  Sven-Bodo Scholz,et al.  Single Assignment C: efficient support for high-level array operations in a functional setting , 2003, Journal of Functional Programming.

[56]  Stefano Severi,et al.  Theoretical investigation of action potential duration dependence on extracellular Ca2+ in human cardiomyocytes. , 2009, Journal of molecular and cellular cardiology.

[57]  J. Clark,et al.  A model of the L‐type Ca2+ channel in rat ventricular myocytes: ion selectivity and inactivation mechanisms , 2000, The Journal of physiology.

[58]  G. W. Beeler,et al.  Reconstruction of the action potential of ventricular myocardial fibres , 1977, The Journal of physiology.

[59]  Robert F Gilmour,et al.  Ionic mechanism of electrical alternans. , 2002, American journal of physiology. Heart and circulatory physiology.

[60]  G. Bett,et al.  Computer model of action potential of mouse ventricular myocytes. , 2004, American journal of physiology. Heart and circulatory physiology.

[61]  Catherine M Lloyd,et al.  CellML: its future, present and past. , 2004, Progress in biophysics and molecular biology.

[62]  David Gavaghan,et al.  Considerations for the use of cellular electrophysiology models within cardiac tissue simulations. , 2011, Progress in biophysics and molecular biology.

[63]  D. Noble,et al.  Distribution of a Persistent Sodium Current Across the Ventricular Wall in Guinea Pigs , 2000, Circulation research.

[64]  J. Clark,et al.  A model of the action potential and underlying membrane currents in a rabbit atrial cell. , 1996, The American journal of physiology.

[65]  M. Courtemanche,et al.  Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. , 1998, The American journal of physiology.

[66]  Yoram Rudy,et al.  Rate Dependence and Regulation of Action Potential and Calcium Transient in a Canine Cardiac Ventricular Cell Model , 2004, Circulation.

[67]  D. Beuckelmann,et al.  Simulation study of cellular electric properties in heart failure. , 1998, Circulation research.