Simulation of QRST integral maps with a membrane-based computer heart model employing parallel processing

The simulation of the propagation of electrical activity in a membrane-based realistic-geometry computer model of the ventricles of the human heart, using the governing monodomain reaction-diffusion equation, is described. Each model point is represented by the phase 1 Luo-Rudy membrane model, modified to represent human action potentials. A separate longer duration action potential was used for the M cells found in the ventricular midwall. Cardiac fiber rotation across the ventricular wall was implemented via an analytic equation, resulting in a spatially varying anisotropic conductivity tensor and, consequently, anisotropic propagation. Since the model comprises approximately 12.5 million points, parallel processing on a multiprocessor computer was used to cut down on simulation time. The simulation of normal activation as well as that of ectopic beats is described. The hypothesis that in situ electrotonic coupling in the myocardium can diminish the gradients of action-potential duration across the ventricular wall was also verified in the model simulations. Finally, the sensitivity of QRST integral maps to local alterations in action-potential duration was investigated.

[1]  S. Abboud,et al.  Simulation of cardiac activity and the ECG using a heart model with a reaction-diffusion action potential. , 1996, Medical engineering & physics.

[2]  K.T. Ng,et al.  A new three-dimensional finite-difference bidomain formulation for inhomogeneous anisotropic cardiac tissues , 1998, IEEE Transactions on Biomedical Engineering.

[3]  R. Beyar,et al.  A Computer Study of the Left Ventricular Performance Based on Fiber Structure, Sarcomere Dynamics, and Transmural Electrical Propagation Velocity , 1984, Circulation research.

[4]  S. Rush,et al.  Qualitative effects of thoracic resistivity variations on the interpretation of electrocardiograms: the low resistance surface layer. , 1968, American Heart Journal.

[5]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. I. Model construction and simulation of normal activation. , 1993, Journal of electrocardiology.

[6]  M. Rosen,et al.  The Controversial M Cell , 1999, Journal of cardiovascular electrophysiology.

[7]  L. Clerc.  Directional differences of impulse spread in trabecular muscle from mammalian heart. , 1976, The Journal of physiology.

[8]  G. Huiskamp,et al.  Simulation of depolarization in a membrane-equations-based model of the anisotropic ventricle , 1998, IEEE Transactions on Biomedical Engineering.

[9]  Y. Rudy,et al.  Imaging Dispersion of Myocardial Repolarization, I: Comparison of Body-Surface and Epicardial Measures , 2001, Circulation.

[10]  A. M. Scher,et al.  Influence of Cardiac Fiber Orientation on Wavefront Voltage, Conduction Velocity, and Tissue Resistivity in the Dog , 1979, Circulation research.

[11]  R. Barr,et al.  Determining surface potentials from current dipoles, with application to electrocardiography. , 1966, IEEE transactions on bio-medical engineering.

[12]  B. Roth.  Electrical conductivity values used with the bidomain model of cardiac tissue , 1997, IEEE Transactions on Biomedical Engineering.

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

[14]  R W Joyner,et al.  Modulation of repolarization by electrotonic interactions. , 1986, Japanese heart journal.

[15]  C. Henriquez,et al.  Modeling impulse propagation and extracellular potential distributions in anisotropic cardiac tissue using a finite volume element discretization , 2002 .

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

[17]  Martin Buist,et al.  Torso Coupling Techniques for the Forward Problem of Electrocardiography , 2002, Annals of Biomedical Engineering.

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

[19]  Craig S. Henriquez,et al.  Using computer models to understand the roles of tissue structure and membrane dynamics in arrhythmogenesis , 1996, Proc. IEEE.

[20]  C. Henriquez,et al.  Bipolar stimulation of a three-dimensional bidomain incorporating rotational anisotropy , 1998, IEEE Transactions on Biomedical Engineering.

[21]  Natalia A. Trayanova,et al.  Computational techniques for solving the bidomain equations in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[22]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[23]  D P Zipes,et al.  Effects of Spatial Segmentation in the Continuous Model of Excitation Propagation in Cardiac Muscle , 1999, Journal of cardiovascular electrophysiology.

[24]  R L Lux,et al.  The Unidentified Information Content of the Electrocardiogram , 1977, Circulation research.

[25]  C S Henriquez,et al.  A space-time adaptive method for simulating complex cardiac dynamics. , 2000, Physical review letters.

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

[27]  J Jalife,et al.  Purkinje-muscle reentry as a mechanism of polymorphic ventricular arrhythmias in a 3-dimensional model of the ventricles. , 1998, Circulation research.

[28]  Michael D. Lesh,et al.  Cellular Uncoupling Can Unmask Dispersion of Action Potential Duration in Ventricular Myocardium A Computer Modeling Study , 1989, Circulation research.

[29]  C. Antzelevitch,et al.  Transient Outward Current Prominent in Canine Ventricular Epicardium but Not Endocardium , 1988, Circulation research.

[30]  R. Gulrajani,et al.  Simulation of propagation in a realistic-geometry computer heart model with parallel processing , 2001, 2001 Conference Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[31]  A. M. Scher,et al.  Effect of Tissue Anisotropy on Extracellular Potential Fields in Canine Myocardium in Situ , 1982, Circulation research.

[32]  B M Horácek,et al.  Computer model of excitation and recovery in the anisotropic myocardium. II. Excitation in the simplified left ventricle. , 1991, Journal of electrocardiology.

[33]  B M Horácek,et al.  Computer model of excitation and recovery in the anisotropic myocardium. I. Rectangular and cubic arrays of excitable elements. , 1991, Journal of electrocardiology.

[34]  E. R. Smith,et al.  Vulnerability to ventricular arrhythmia: assessment by mapping of body surface potential. , 1986, Circulation.

[35]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[36]  C. Antzelevitch,et al.  Characteristics and distribution of M cells in arterially perfused canine left ventricular wedge preparations. , 1998, Circulation.

[37]  R. Penland,et al.  A flexible method for simulating cardiac conduction in three-dimensional complex geometries. , 2000, Journal of electrocardiology.

[38]  D. Geselowitz,et al.  Simulation Studies of the Electrocardiogram: I. The Normal Heart , 1978, Circulation research.

[39]  A. J. Pullan,et al.  Mathematical models and numerical methods for the forward problem in cardiac electrophysiology , 2002 .

[40]  Ramesh M. Gulrajani,et al.  A Membrane-Based Computer Heart Model Employing Parallel Processing , 2001 .

[41]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. II. Simulations of conduction block. , 1993, Journal of electrocardiology.

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

[43]  N. Hayat,et al.  Asymptomatic severe aortic regurgitation in the young. Hemodynamic and angiographic data in 37 patients with echocardiographic correlation. , 1986, Japanese heart journal.