Toward real-time modeling of human heart ventricles at cellular resolution: Simulation of drug-induced arrhythmias

We have developed a highly efficient and scalable cardiac electrophysiology simulation capability that supports groundbreaking resolution and detail to elucidate the mechanisms of sudden cardiac death from arrhythmia. We can simulate thousands of heartbeats at a resolution of 0.1 mm, comparable to the size of cardiac cells, thereby enabling scientific inquiry not previously possible. Based on scaling results from the partially deployed Sequoia IBM Blue Gene/Q machine at Lawrence Livermore National Laboratory and planned optimizations, we estimate that by SC12 we will simulate 8 -- 10 heartbeats per minute -- a time-to-solution 400 -- 500 times faster than the state-of-the-art. Performance between 8 and 11 PFlop/s on the full 1,572,864 cores is anticipated, representing 40 -- 55 percent of peak. The power of the model is demonstrated by illuminating the subtle arrhythmogenic mechanisms of anti-arrhythmic drugs that paradoxically increase arrhythmias in some patient populations.

[1]  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.

[2]  N. Trayanova,et al.  A Computational Model to Predict the Effects of Class I Anti-Arrhythmic Drugs on Ventricular Rhythms , 2011, Science Translational Medicine.

[3]  F. Charpentier,et al.  Electrophysiologic characteristics of cells spanning the left ventricular wall of human heart: evidence for presence of M cells. , 1995, Journal of the American College of Cardiology.

[4]  Eric Kerfoot,et al.  Verification of cardiac tissue electrophysiology simulators using an N-version benchmark , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[5]  Lawrence Mitchell,et al.  Simulating Human Cardiac Electrophysiology on Clinical Time-Scales , 2011, Front. Physio..

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

[7]  I. Efimov,et al.  Transmural Dispersion of Repolarization in Failing and Nonfailing Human Ventricle , 2010, Circulation research.

[8]  CHARLES ANTZELEVITCH,et al.  The M Cell: , 1999, Journal of cardiovascular electrophysiology.

[9]  R. Virmani,et al.  Sudden cardiac death. , 1987, Human pathology.

[10]  C Antzelevitch,et al.  Larger late sodium conductance in M cells contributes to electrical heterogeneity in canine ventricle. , 2001, American journal of physiology. Heart and circulatory physiology.

[11]  C. Antzelevitch,et al.  Unique Topographical Distribution of M Cells Underlies Reentrant Mechanism of Torsade de Pointes in the Long-QT Syndrome , 2002, Circulation.

[12]  Matthias Reumann,et al.  Orthogonal recursive bisection as data decomposition strategy for massively parallel cardiac simulations , 2011, Biomedizinische Technik. Biomedical engineering.

[13]  N. Estes,et al.  Predicting and preventing sudden cardiac death. , 2011, Circulation.

[14]  D. Durrer,et al.  Total Excitation of the Isolated Human Heart , 1970, Circulation.

[15]  D. Zipes,et al.  Antiarrhythmic and Nonantiarrhythmic Drugs for Sudden Cardiac Death Prevention , 2010, Journal of cardiovascular pharmacology.

[16]  Olaf Dössel,et al.  Influence of ${I_{Ks}}$ Heterogeneities on the Genesis of the T-wave: A Computational Evaluation , 2012, IEEE Transactions on Biomedical Engineering.

[17]  C. Antzelevitch,et al.  Overview of Basic Mechanisms of Cardiac Arrhythmia. , 2011, Cardiac electrophysiology clinics.

[18]  Mark Potse,et al.  Poster: hybrid parallelization of a realistic heart model , 2011, SC '11 Companion.

[19]  Matthias Reumann,et al.  Performance of Hybrid Programming Models for Multiscale Cardiac Simulations: Preparing for Petascale Computation , 2011, IEEE Transactions on Biomedical Engineering.

[20]  Borje Darpo,et al.  Electrocardiographic Identification of Drug‐Induced QT Prolongation: Assessment by Different Recording and Measurement Methods , 2004, Annals of noninvasive electrocardiology : the official journal of the International Society for Holter and Noninvasive Electrocardiology, Inc.

[21]  Michael Gschwind,et al.  The IBM Blue Gene/Q Compute Chip , 2012, IEEE Micro.

[22]  Mark Potse,et al.  A Comparison of Monodomain and Bidomain Reaction-Diffusion Models for Action Potential Propagation in the Human Heart , 2006, IEEE Transactions on Biomedical Engineering.

[23]  S. Rush,et al.  A Practical Algorithm for Solving Dynamic Membrane Equations , 1978, IEEE Transactions on Biomedical Engineering.