Minimal model for human ventricular action potentials in tissue

Modeling the dynamics of wave propagation in human ventricular tissue and studying wave stability require models that reproduce realistic characteristics in tissue. We present a minimal ventricular (MV) human model that is designed to reproduce important tissue-level characteristics of epicardial, endocardial and midmyocardial cells, including action potential (AP) amplitudes and morphologies, upstroke velocities, steady-state action potential duration (APD) and conduction velocity (CV) restitution curves, minimum APD, and minimum diastolic interval. The model is then compared with three previously published human ventricular cell models, the Priebe and Beuckelmann (PB), the Ten Tusscher–Noble–Noble–Panfilov (TNNP), and the Iyer–Mazhari–Winslow (IMW). For the first time, the stability of reentrant waves for all four models is analyzed, and quantitative comparisons are made among the models in single cells and in tissue. The PB, TNNP, and IMW models exhibit quantitative differences in APD and CV rate adaptation, as well as completely different reentrant wave dynamics of quasi-breakup, stability, and breakup, respectively. All the models have dominant frequencies comparable to clinical values except for the IMW model, which has a large range of frequencies extending beyond the clinical range for both ventricular tachycardia (VT) and ventricular fibrillation (VF). The TNNP and IMW models possess a large degree of short-term memory and we show for the first time the existence of memory in CV restitution. The MV model also can be fitted to reproduce the dynamics of other models and is computationally more efficient: the times required to simulate the MV, TNNP, PB and IMW models follow the ratio 1:31:50:8084.     

A computationally efficient electrophysiological model of human ventricular cells

Recent experimental and theoretical results have stressed the importance of modeling studies of reentrant arrhythmias in cardiac tissue and at the whole heart level. We introduce a six-variable model obtained by a reformulation of the Priebe-Beuckelmann model of a single human ventricular cell. The reformulated model is 4.9 times faster for numerical computations and it is more stable than the original model. It retains the action potential shape at various frequencies, restitution of action potential duration, and restitution of conduction velocity. We were able to reproduce the main properties of epicardial, endocardial, and M cells by modifying selected ionic currents. We performed a simulation study of spiral wave behavior in a two-dimensional sheet of human ventricular tissue and showed that spiral waves have a frequency of 3.3 Hz and a linear core of approximately 50-mm diameter that rotates with an average frequency of 0.62 rad/s. Simulation results agreed with experimental data. In conclusion, the proposed model is suitable for efficient and accurate studies of reentrant phenomena in human ventricular tissue.     

A model for human ventricular tissue

The experimental and clinical possibilities for studying cardiac arrhythmias in human ventricular myocardium are very limited. Therefore, the use of alternative methods such as computer simulations is of great importance. In this article we introduce a mathematical model of the action potential of human ventricular cells that, while including a high level of electrophysiological detail, is computationally cost-effective enough to be applied in large-scale spatial simulations for the study of reentrant arrhythmias. The model is based on recent experimental data on most of the major ionic currents: the fast sodium, L-type calcium, transient outward, rapid and slow delayed rectifier, and inward rectifier currents. The model includes a basic calcium dynamics, allowing for the realistic modeling of calcium transients, calcium current inactivation, and the contraction staircase. We are able to reproduce human epicardial, endocardial, and M cell action potentials and show that differences can be explained by differences in the transient outward and slow delayed rectifier currents. Our model reproduces the experimentally observed data on action potential duration restitution, which is an important characteristic for reentrant arrhythmias. The conduction velocity restitution of our model is broader than in other models and agrees better with available data. Finally, we model the dynamics of spiral wave rotation in a two-dimensional sheet of human ventricular tissue and show that the spiral wave follows a complex meandering pattern and has a period of 265 ms. We conclude that the proposed model reproduces a variety of electrophysiological behaviors and provides a basis for studies of reentrant arrhythmias in human ventricular tissue.     

Comparison of electrophysiological models for human ventricular cells and tissues

In this paper we briefly review currently published models for human ventricular cells and tissues. We discuss the Priebe–Beuckelmann (PB) model and the reduced version of this model constructed by Bernus et al. (redPB), the Ten Tusscher–Noble–Noble–Panfilov (TNNP) model and the Iyer–Mazhari–Winslow (IMW) model. We compare several characteristics of these models such as: sources of experimental data the models are based on, action potential morphology, action potential duration (APD) and conduction velocity (CV) restitution and computational efficiency. Finally, we discuss the application of a subset of these models—the redPB and the TNNP model—to study simulated spiral wave dynamics in 2D tissue sheets and in the human ventricles. We discuss the suitability of the different models for particular research questions and their limitations.     

Effects of Mechano-Electric Feedback on Scroll Wave Stability in Human Ventricular Fibrillation

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, −60 mV and −10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of −60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of −10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of −10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.     

Electrical refractory period restitution and spiral wave reentry in simulated cardiac tissue

Theoretical and experimental studies have shown that restitution of the cardiac action potential (AP) duration (APD) plays a major role in predisposing ventricular tachycardia to degenerate to ventricular fibrillation, whereas its role in atrial fibrillation is unclear. We used the Courtemanche human atrial cell model and the Luo-Rudy guinea pig ventricular model to compare the roles of electrical restitution in destabilizing spiral wave reentry in simulated two-dimensional homogeneous atrial and ventricular tissue. Because atrial AP morphology is complex, we also validated the usefulness of effective refractory period (ERP) restitution. ERP restitution correlated best with APD restitution at transmembrane potentials greater than or equal to -62 mV, and its steepness was a reliable predictor of spiral wave phenotype (stable, meandering, hypermeandering, and breakup) in both atrial and ventricular tissue. Spiral breakup or single hypermeandering spirals occurred when the slope of ERP restitution exceeded 1 at short diastolic intervals. Thus ERP restitution, which is easier to measure clinically than APD restitution, is a reliable determinant of spiral wave stability in simulated atrial and ventricular tissue.     

A Discrete Electromechanical Model for Human Cardiac Tissue: Effects of Stretch-Activated Currents and Stretch Conditions on Restitution Properties and Spiral Wave Dynamics

We introduce an electromechanical model for human cardiac tissue which couples a biophysical model of cardiac excitation (Tusscher, Noble, Noble, Panfilov, 2006) and tension development (adjusted Niederer, Hunter, Smith, 2006 model) with a discrete elastic mass-lattice model. The equations for the excitation processes are solved with a finite difference approach, and the equations of the mass-lattice model are solved using Verlet integration. This allows the coupled problem to be solved with high numerical resolution. Passive mechanical properties of the mass-lattice model are described by a generalized Hooke''s law for finite deformations (Seth material). Active mechanical contraction is initiated by changes of the intracellular calcium concentration, which is a variable of the electrical model. Mechanical deformation feeds back on the electrophysiology via stretch-activated ion channels whose conductivity is controlled by the local stretch of the medium. We apply the model to study how stretch-activated currents affect the action potential shape, restitution properties, and dynamics of spiral waves, under constant stretch, and dynamic stretch caused by active mechanical contraction. We find that stretch conditions substantially affect these properties via stretch-activated currents. In constantly stretched medium, we observe a substantial decrease in conduction velocity, and an increase of action potential duration; whereas, with dynamic stretch, action potential duration is increased only slightly, and the conduction velocity restitution curve becomes biphasic. Moreover, in constantly stretched medium, we find an increase of the core size and period of a spiral wave, but no change in rotation dynamics; in contrast, in the dynamically stretching medium, we observe spiral drift. Our results may be important to understand how altered stretch conditions affect the heart''s functioning.     

Model Studies of the Role of Mechano-sensitive Currents in the Generation of Cardiac Arrhythmias

Mechano-electrical feedback is studied by incorporating linear, instantaneously activating mechano-sensitive conductances into single cardiac cell models, as well as one- and two-dimensional cardiac network models. The models qualitatively reproduce effects of maintained mechanical stretch on experimentally measured action potential characteristics such as amplitude, maximum diastolic potential, peak upstroke velocity, and conduction velocity. Models are also used to simulate stretch-induced depolarizations, action potentials, and arrhythmias produced by pulsatile volume changes in left ventricle of dog. The mechano-sensitive conductance threshold for a stretch-induced action potential is closely related to the magnitude of the time-independent K⁺current,I_K1, which offsets inward mechano-sensitive current. Activation of mechano-sensitive conductances in small, spatially localized region of cells can evoke graded depolarizations, propagating ectopic beats, and if timed appropriately, spiral reentrant waves. Mechano-sensitive conductance changes required to evoke these responses are well within the physiologically plausible range. Results therefore indicate that many mechano-electrical feedback effects can be modeled using linear, instantaneously activating mechano-sensitive conductances. As an example of how stretch can occur in real human hearts, magnetic resonance images with saturation tagging are used to reconstruct the three-dimensional left ventricular wall motion. In patients with infarcts or recent ischemic events, “paradoxical deformation” is observed in that regions of myocardium are stretched rather than contracted during systole. In contrast, normal hearts contract uniformly with no stretch during systole. Paradoxical deformations in ischemic hearts may therefore present one possible substrate for the mechanically induced arrhythmias modeled above.     

The Transfer Functions of Cardiac Tissue during Stochastic Pacing

The restitution properties of cardiac action potential duration (APD) and conduction velocity (CV) are important factors in arrhythmogenesis. They determine alternans, wavebreak, and the patterns of reentrant arrhythmias. We developed a novel approach to characterize restitution using transfer functions. Transfer functions relate an input and an output quantity in terms of gain and phase shift in the complex frequency domain. We derived an analytical expression for the transfer function of interbeat intervals (IBIs) during conduction from one site (input) to another site downstream (output). Transfer functions can be efficiently obtained using a stochastic pacing protocol. Using simulations of conduction and extracellular mapping of strands of neonatal rat ventricular myocytes, we show that transfer functions permit the quantification of APD and CV restitution slopes when it is difficult to measure APD directly. We find that the normally positive CV restitution slope attenuates IBI variations. In contrast, a negative CV restitution slope (induced by decreasing extracellular [K⁺]) amplifies IBI variations with a maximum at the frequency of alternans. Hence, it potentiates alternans and renders conduction unstable, even in the absence of APD restitution. Thus, stochastic pacing and transfer function analysis represent a powerful strategy to evaluate restitution and the stability of conduction.     

Steeper restitution slopes across right ventricular endocardium in patients with cardiomyopathy at high risk of ventricular arrhythmias

Steep action potential duration (APD) restitution slopes (>1) and spatial APD restitution heterogeneity provide the substrate for ventricular fibrillation in computational models and experimental studies. Their relationship to ventricular arrhythmia vulnerability in human cardiomyopathy has not been defined. Patients with cardiomyopathy [left ventricular (LV) ejection fraction <40%] and no history of ventricular arrhythmias underwent risk stratification with programmed electrical stimulation or T wave alternans (TWA). Low-risk patients (n = 10) had no inducible ventricular tachycardia (VT) or negative TWA, while high-risk patients (n = 8) had inducible VT or positive TWA. Activation recovery interval (ARI) restitution slopes were measured simultaneously from 10 right ventricular (RV) endocardial sites during an S1-S2 pacing protocol. ARI restitution slope heterogeneity was defined as the coefficient of variation of slopes. Mean ARI restitution slope was significantly steeper in the high-risk group compared with the low-risk group [1.16 (SD 0.31) vs. 0.59 (SD 0.19), P = 0.0002]. The proportion of endocardial recording sites with a slope >1 was significantly larger in the high-risk patients [47% (SD 35) vs. 13% (SD 21), P = 0.022]. Spatial heterogeneity of ARI restitution slopes was similar between the two groups [29% (SD 16) vs. 39% (SD 34), P = 0.48]. There was an inverse linear relationship between the ARI restitution slope and the minimum diastolic interval (P < 0.001). In cardiomyopathic patients at high risk of ventricular arrhythmias, ARI restitution slopes along the RV endocardium are steeper, but restitution slope heterogeneity is similar compared with those at low risk. Steeper ARI restitution slopes may increase the propensity for ventricular arrhythmias in patients with impaired left ventricular function.     

Properties of two human atrial cell models in tissue: restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.     

Vulnerable window for conduction block in a one-dimensional cable of cardiac cells, 2: multiple extrasystoles

Unidirectional conduction block of premature extrasystoles can lead to initiation of cardiac reentry, causing lethal arrhythmias including ventricular fibrillation. Multiple extrasystoles are often more effective at inducing unidirectional conduction block and reentry than a single extrasystole. Since the substrate for conduction block is spatial dispersion of refractoriness, in this study we investigate how the first extrasystole modulates this dispersion to influence the "vulnerable window" for conduction block by subsequent extrasystoles, particularly in relation to action potential duration restitution and conduction velocity restitution properties. Using a kinematic model to represent wavefront-waveback interactions and simulations with the Luo-Rudy model in a one-dimensional cable of cardiac cells, we show that in homogeneous tissue, a premature extrasystole can create a large dispersion of refractoriness leading to conduction block of a subsequent extrasystole. In heterogeneous tissue, however, a premature extrasystole can either reduce or enhance the dispersion of refractoriness depending on its propagation direction with respect to the previous beat. With multiple extrasystoles at random coupling intervals, vulnerability to conduction block is proportional to their number. In general, steep action potential duration restitution and broad conduction velocity restitution promote dispersion of refractoriness in response to multiple extrasystoles, and thus enhance vulnerability to conduction block. These restitution properties also promote spatially discordant alternans, a setting which is particularly prone to conduction block. The equivalent dispersion of refractoriness created dynamically in homogeneous tissue by spatially discordant alternans is more likely to cause conduction block than a comparable degree of preexisting dispersion in heterogeneous tissue.     

Intracellular Ca(2+) dynamics and the stability of ventricular tachycardia

Ventricular fibrillation (VF), the major cause of sudden cardiac death, is typically preceded by ventricular tachycardia (VT), but the mechanisms underlying the transition from VT to VF are poorly understood. Intracellular Ca(2+) overload occurs during rapid heart rates typical of VT and is also known to promote arrhythmias. We therefore studied the role of intracellular Ca(2+) dynamics in the transition from VT to VF, using a combined experimental and mathematical modeling approach. Our results show that 1) rapid pacing of rabbit ventricular myocytes at 35 degrees C led to increased intracellular Ca(2+) levels and complex patterns of action potential (AP) configuration and the intracellular Ca(2+) transients; 2) the complex patterns of the Ca(2+) transient arose directly from the dynamics of intracellular Ca(2+) cycling, and were not merely passive responses to beat-to-beat alterations in AP; 3) the complex Ca(2+) dynamics were simulated in a modified version of the Luo-Rudy (LR) ventricular action potential with improved intracellular Ca(2+) dynamics, and showed good agreement with the experimental findings in isolated myocytes; and 4) when incorporated into simulated two-dimensional cardiac tissue, this action potential model produced a form of spiral wave breakup from VT to a VF-like state in which intracellular Ca(2+) dynamics played a key role through its influence on Ca(2+)-sensitive membrane currents such as I(Ca), I(NaCa), and I(ns(Ca)). To the extent that spiral wave breakup is useful as a model for the transition from VT to VF, these findings suggest that intracellular Ca(2+) dynamics may play an important role in the destabilization of VT and its degeneration into VF.     

Mechanisms of destabilization and early termination of spiral wave reentry in the ventricle by a class III antiarrhythmic agent, nifekalant

Nifekalant (NF) is a novel class III antiarrhythmic agent that is effective in preventing life-threatening ventricular tachycardia/fibrillation (VT/VF). We investigated mechanisms of destabilization and early termination of spiral-type reentrant VT by NF in a two-dimensional subepicardial myocardial layer of Langendorff-perfused rabbit hearts (n = 21) using a high-resolution optical action potential mapping system. During basic stimulation, NF (0.1 microM) caused uniform prolongation of action potential duration (APD) without affecting conduction velocity and an increase of APD restitution slope. VTs induced by direct current stimulation in the presence of NF were of shorter duration (VTs > 30 s: 2/54 NF vs. 19/93 control). During VTs in control (with visible rotors), the wave front chased its own tail with a certain distance (repolarized zone), and they seldom met each other. The average number of phase singularity (PS) points was 1.31 +/- 0.14 per 665 ms (n = 7). In the presence of NF, the wave front frequently encountered its own tail, causing a transient breakup of the spiral wave or sudden movement of the rotation center (spatial jump of PS). The average number of PS was increased to 1.63 +/- 0.22 per 665 ms (n = 7, P < 0.05) after NF. The mode of spontaneous termination of rotors in control was in most cases (9/10, 90.0%) the result of mutual annihilation of counterrotating wave fronts. With NF, rotors frequently terminated by wave front collision with the atrioventricular groove (12/19, 63.2%) or by trapping the spiral tip in a refractory zone (7/19, 36.8%). Destabilization and early termination of spiral wave reentry induced by NF are the result of a limited proportion of excitable tissue after modulation of repolarization.     

Sildenafil-nitric oxide donor combination promotes ventricular tachyarrhythmias in the swine right ventricle

We tested the hypothesis that sildenafil, singly or in combination with nitric oxide (NO) donors, promotes ventricular tachycardia (VT) and ventricular fibrillation (VF). Vulnerability to VT/VF was tested by rapid pacing in eight isolated normal swine right ventricles (RV). The endocardial activation was optically mapped, and the dynamic action potential duration (APD) restitution curves were constructed with metal microelectrodes. At baseline, no VT/VF could be induced. Sildenafil (0.2 microg/ml) or NO donor singly or in combination did not alter VT/VF vulnerability. However, when 2 microg/ml sildenafil was combined with NO donors, the incidence of VT and VF rose significantly (P < 0.01). VT with a single periodic wavefront was induced in five of eight RVs, and VF with multiple wavefronts was induced in all eight RVs. The sildenafil-NO donor pro-VT/VF combination significantly increased the maximum slope of the APD restitution curve and the amplitude of the APD alternans. The pro-VT/VF effects of sildenafil were reversible after drug-free Tyrode solution perfusion. We conclude that a sildenafil (2 microg/ml) and NO donor combination increases VT/VF vulnerability in the normal RV by a mechanism compatible with the restitution hypothesis.     

The pinwheel experiment revisited: effects of cellular electrophysiological properties on vulnerability to cardiac reentry

In normal heart, ventricular fibrillation can be induced by a single properly timed strong electrical or mechanical stimulus. A mechanism first proposed by Winfree and coined the "pinwheel experiment" emphasizes the timing and strength of the stimulus in inducing figure-of-eight reentry. However, the effects of cellular electrophysiological properties on vulnerability to reentry in the pinwheel scenario have not been investigated. In this study, we extend Winfree's pinwheel experiment to show how the vulnerability to reentry is affected by the graded action potential responses induced by a strong premature stimulus, action potential duration (APD), and APD restitution in simulated monodomain homogeneous two-dimensional tissue. We find that a larger graded response, longer APD, or steeper APD restitution slope reduces the vulnerable window of reentry. Strong graded responses and long APD promote tip-tip interactions at long coupling intervals, causing the two initiated spiral wave tips to annihilate. Steep APD restitution promotes wave front-wave back interaction, causing conduction block in the central common pathway of figure-of-eight reentry. We derive an analytical treatment that shows good agreement with numerical simulation results.     

Analysis of damped oscillations during reentry: a new approach to evaluate cardiac restitution

Reentry is a mechanism underlying numerous cardiac arrhythmias. During reentry, head-tail interactions of the action potential can cause cycle length (CL) oscillations and affect the stability of reentry. We developed a method based on a difference-delay equation to determine the slopes of the action potential duration and conduction velocity restitution functions, known to be major determinants of reentrant arrhythmogenesis, from the spatial period P and the decay length D of damped CL oscillations. Using this approach, we analyzed CL oscillations after the induction of reentry and the resetting of reentry with electrical stimuli in rings of cultured neonatal rat ventricular myocytes grown on microelectrode arrays and in corresponding simulations with the Luo-Rudy model. In the experiments, P was larger and D was smaller after resetting impulses compared to the induction of reentry, indicating that reentry became more stable. Both restitution slopes were smaller. Consistent with the experimental findings, resetting of simulated reentry caused oscillations with gradually increasing P, decreasing D, and decreasing restitution slopes. However, these parameters remained constant when ion concentrations were clamped, revealing that intracellular ion accumulation stabilizes reentry. Thus, the analysis of CL oscillations during reentry opens new perspectives to gain quantitative insight into action potential restitution.     

Propagation of action potentials in squid giant axons. Repetitive firing at regions of membrane inhomogeneities

Effects of reduction in potassium conductance on impulse conduction were studied in squid giant axons. Internal perfusion of axons with tetraethylammonium (TEA) ions reduces G K and causes the duration of action potential to be increased up to 300 ms. This prolongation of action potentials does not change their conduction velocity. The shape of these propagating action potentials is similar to membrane action potentials in TEA. Axons with regions of differing membrane potassium conductances are obtained by perfusing the axon trunk and one of its two main branches with TEA after the second branch has been filled with normal perfusing solution. Although the latter is initially free of TEA, this ion diffuses in slowly. Up until a large amount of TEA has diffused into the second branch, action potentials in the two branches have very different durations. During this period, membrane regions with prolonged action potentials are a source of depolarizing current for the other, and repetitive activity may be initiated at transitional regions. After a single stimulus in either axon region, interactions between action potentials of different durations usually led to rebound, or a short burst, of action potentials. Complex interactions between two axon regions whose action potentials have different durations resembles electric activity recorded during some cardiac arrhythmias.     

The electrophysiological effects of disopyramide phosphate on canine ventricular muscle and Purkinje fibers in normal and low potassium

We studied the effect of lowering the extracellular potassium concentration ([K+]o) on the electrophysiological actions of disopyramide phosphate, a new antiarrhythmic drug. At low [K+]o, therapeutic concentrations of disopyramide phosphate caused significantly less depression of action potential amplitude and maximum upstroke velocity of both Purkinje fiber and ventricular muscle action potentials. The drug shifted the membrane responsiveness curve along the voltage axis to more negative membrane potentials regardless of [K+]o. However, a greater shift occurred when [K+]o was normal. Disopyramide phosphate prolonged both action potential duration and effective refractory period in all fibers but there was consistently greater prolongation of these parameters at low [K+]o. More importantly, disopyramide phosphate altered repolarization time course of action potentials in such a way that action potentials with dissimilar durations throughout the ventricular conducting system became more equal. The drug was less effective in decreasing this disparity in action potential durations throughout the ventricles in the presence of low [K+]o. These modifications of the electrophysiological actions of disopyramide by low [K+]o suggest that a therapeutic concentration of disopyramide might have less of an antiarrhythmic effect in the presence of hypokalemia.