Towards real-time simulation of cardiac electrophysiology in a human heart at high resolution

We have developed the capability to rapidly simulate cardiac electrophysiological phenomena in a human heart discretised at a resolution comparable with the length of a cardiac myocyte. Previous scientific investigation has generally invoked simplified geometries or coarse-resolution hearts, with simulation duration limited to 10s of heartbeats. Using state-of-the-art high-performance computing techniques coupled with one of the most powerful computers available (the 20 PFlop/s IBM BlueGene/Q at Lawrence Livermore National Laboratory), high-resolution simulation of the human heart can now be carried out over 1200 times faster compared with published results in the field. We demonstrate the utility of this capability by simulating, for the first time, the formation of transmural re-entrant waves in a 3D human heart. Such wave patterns are thought to underlie Torsades de Pointes, an arrhythmia that indicates a high risk of sudden cardiac death. Our new simulation capability has the potential to impact a multitude of applications in medicine, pharmaceuticals and implantable devices.

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

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

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

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

[5]  Lin Wu,et al.  Late Sodium Current Contributes to the Reverse Rate-Dependent Effect of IKr Inhibition on Ventricular Repolarization , 2011, Circulation.

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

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

[8]  G. Plank,et al.  A Novel Rule-Based Algorithm for Assigning Myocardial Fiber Orientation to Computational Heart Models , 2012, Annals of Biomedical Engineering.

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

[10]  Mark Potse,et al.  Cardiac anisotropy in boundary-element models for the electrocardiogram , 2009, Medical & Biological Engineering & Computing.

[11]  Olaf Dössel,et al.  Ranking the Influence of Tissue Conductivities on Forward-Calculated ECGs , 2010, IEEE Transactions on Biomedical Engineering.

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

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

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

[15]  Viatcheslav Gurev,et al.  Toward real-time modeling of human heart ventricles at cellular resolution: Simulation of drug-induced arrhythmias , 2012, 2012 International Conference for High Performance Computing, Networking, Storage and Analysis.

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

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