High-throughput cardiac science on the Grid

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.

[1]  Godfrey L. Smith,et al.  Heterogeneous changes in action potential and intracellular Ca2+ in left ventricular myocyte sub-types from rabbits with heart failure. , 2000, Cardiovascular research.

[2]  D Gavaghan,et al.  Grid computing simulations of ion channel block effects on the ECG using 3D anatomically-based models , 2009, 2009 36th Annual Computers in Cardiology Conference (CinC).

[3]  Steve McKeever,et al.  On the application of partial evaluation to the optimisation of cardiac electrophysiological simulations , 2006, PEPM '06.

[4]  David Abramson,et al.  Incorporating Local Ca2+ Dynamics into Single Cell Ventricular Models , 2008, ICCS.

[5]  Ian Foster,et al.  The Grid 2 - Blueprint for a New Computing Infrastructure, Second Edition , 1998, The Grid 2, 2nd Edition.

[6]  D. Noble.  A modification of the Hodgkin—Huxley equations applicable to Purkinje fibre action and pacemaker potentials , 1962, The Journal of physiology.

[7]  D. Noble,et al.  Rectifying Properties of Heart Muscle , 1960, Nature.

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

[9]  David Abramson,et al.  Robust workflows for science and engineering , 2009, MTAGS '09.

[10]  David Abramson,et al.  High performance parametric modeling with Nimrod/G: killer application for the global grid? , 2000, Proceedings 14th International Parallel and Distributed Processing Symposium. IPDPS 2000.

[11]  Julian Cummings,et al.  Comparison of C++ and Fortran 90 for object-oriented scientific programming , 1997 .

[12]  Bennet Vance,et al.  On the application of partial evaluation to database optimization , 1993 .

[13]  Andrew Lewis,et al.  NIMROD/O: A TOOL FOR AUTOMATIC DESIGN OPTIMISATION USING PARALLEL AND DISTRIBUTED SYSTEMS , 2000 .

[14]  M.O. Bernabeu,et al.  High performance computer simulations for the study of biological function in 3D heart models incorporating fibre orientation and realistic geometry at para-cellular resolution , 2008, 2008 Computers in Cardiology.

[15]  D. Noble.  Cardiac Action and Pacemaker Potentials based on the Hodgkin-Huxley Equations , 1960, Nature.

[16]  D. Noble,et al.  Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. , 1998, The Canadian journal of cardiology.

[17]  David Abramson,et al.  Nimrod: a tool for performing parametrised simulations using distributed workstations , 1995, Proceedings of the Fourth IEEE International Symposium on High Performance Distributed Computing.

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

[19]  Bertram Ludäscher,et al.  A Framework for the Design and Reuse of Grid Workflows , 2004, SAG.

[20]  James P. Keener,et al.  Mathematical physiology , 1998 .

[21]  A. Garny,et al.  Chaste: incorporating a novel multi-scale spatial and temporal algorithm into a large-scale open source library , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  A Fitton,et al.  Sotalol , 1993, Drugs.

[23]  David Abramson,et al.  Application deployment over heterogeneous grids using distributed ant , 2005, First International Conference on e-Science and Grid Computing (e-Science'05).

[24]  M. Sansom,et al.  Ion channels: structural bioinformatics and modelling. , 2002, Human molecular genetics.

[25]  A. Tanskanen,et al.  A simplified local control model of calcium-induced calcium release in cardiac ventricular myocytes. , 2004, Biophysical journal.

[26]  Yoram Rudy,et al.  Computational biology in the study of cardiac ion channels and cell electrophysiology , 2006, Quarterly Reviews of Biophysics.

[27]  David Abramson,et al.  Nimrod/K: Towards massively parallel dynamic Grid workflows , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[28]  G Plank,et al.  Evaluating intramural virtual electrodes in the myocardial wedge preparation: simulations of experimental conditions. , 2008, Biophysical journal.

[29]  David Abramson,et al.  Fractional factorial design for parameter sweep experiments using Nimrod/E , 2008, Sci. Program..

[30]  Ami Marowka,et al.  The GRID: Blueprint for a New Computing Infrastructure , 2000, Parallel Distributed Comput. Pract..

[31]  D. Abramson,et al.  An Automatic Design Optimization Tool and its Application to Computational Fluid Dynamics , 2001, ACM/IEEE SC 2001 Conference (SC'01).

[32]  Roy C. P. Kerckhoffs,et al.  Computational Methods for Cardiac Electromechanics , 2006, Proceedings of the IEEE.

[33]  David Abramson,et al.  Fractional factorial design for parameter sweep experiments using Nimrod/E , 2008 .

[34]  Joseph L Greenstein,et al.  Mechanisms of excitation-contraction coupling in an integrative model of the cardiac ventricular myocyte. , 2006, Biophysical journal.