Ventricular Cell Model Research Articles

Simulating cardiac ventricular action potentials in rat and mouse.

We have developed a computational model for the rat cardiac ventricular cell by using electrophysiology data. Our control model (Pandit et al., 2001). The model represents the bioelectric activity in the left ventricular cells in adult male rats. We have described the differences in the membrane properties within the left ventricle to simulate the action potential variations of the endocardial and epicardial cells. We also developed a right ventricular cell model from our control model (the left ventricular cell model) to study ionic mechanisms in diabetic rats (Pandit et al., 2003) Our right ventricular cell model was also the template for us to develop a mouse ventricular cell model by using experimental data. The simulation studies in our models of the rat and mouse ventricular cells conclude that (1) the action potential changes in the left and right ventricles of the normal, diabetic, aged or spontaneously hypertensive rats is mainly due to differences in the 4AP sensitive, Ca2+independent transient outward K+current (I<inf>t</inf>) current, and (2) the presence of the 4AP sensitive, slowly inactivating, delayed rectifier K+current (I<inf>Kslow</inf>) in mouse is one of the main factors contributing to the faster rate of repolarization seen in mouse compared to rats.

Gender Disparities in Cardiac Cellular Electrophysiology and Arrhythmia Susceptibility in Human Failing Ventricular Myocytes

Gender disparities in ECG variables and susceptibility to arrhythmia exist. The basis of these sex-related distinctions in cardiac electrophysiology has been extensively studied in various species, but is virtually unexplored in humans. The aim of this study was to clarify the cellular basis of electrophysiological gender disparities in human cardiac myocytes. Human midmyocardial left ventricular myocytes were isolated from explanted hearts of male and female patients in end-stage heart failure at the time of cardiac transplantation. The action potentials, sarcolemmal ion currents, and susceptibility to the generation of early afterdepolarizations were studied using whole-cell patch-clamp methodology. The functional effects of gender disparities in sarcolemmal ion currents were assessed by computer simulations using the Priebe-Beuckelmann or the ten Tusscher-Noble-Noble-Panfilov human ventricular cell models. Female myocytes had significantly longer action potentials and greater susceptibility to early afterdepolarizations than male myocytes. All other action potential parameters (resting membrane potential, amplitude, plateau level, upstroke velocity, maximal velocity of phase-1 and phase-3 repolarization) had similar values for both genders. In female myocytes, the transient outward potassium current (I(to1)) tended to be smaller, while the L-type calcium current (I(Ca,L)) and quasi-steady state current (I(QSS)) tended to be larger. Computer simulations showed that these subtle differences in sarcolemmal ion currents may conspire to cause the observed gender disparities in action potential properties. Female failing myocytes have longer action potentials and a greater susceptibility to early afterdepolarizations than male failing myocytes. These gender disparities may be due to slightly larger depolarizing I(Ca,L) in conjunction with slightly smaller repolarizing I(QSS) and I(to1) in female myocytes.

Rate Dependence and Regulation of Action Potential and Calcium Transient in a Canine Cardiac Ventricular Cell Model

Computational biology is a powerful tool for elucidating arrhythmogenic mechanisms at the cellular level, where complex interactions between ionic processes determine behavior. A novel theoretical model of the canine ventricular epicardial action potential and calcium cycling was developed and used to investigate ionic mechanisms underlying Ca2+ transient (CaT) and action potential duration (APD) rate dependence. The Ca2+/calmodulin-dependent protein kinase (CaMKII) regulatory pathway was integrated into the model, which included a novel Ca2+-release formulation, Ca2+ subspace, dynamic chloride handling, and formulations for major ion currents based on canine ventricular data. Decreasing pacing cycle length from 8000 to 300 ms shortened APD primarily because of I(Ca(L)) reduction, with additional contributions from I(to1), I(NaK), and late I(Na). CaT amplitude increased as cycle length decreased from 8000 to 500 ms. This positive rate-dependent property depended on CaMKII activity. CaMKII is an important determinant of the rate dependence of CaT but not of APD, which depends on ion-channel kinetics. The model of CaMKII regulation may serve as a paradigm for modeling effects of other regulatory pathways on cell function.

In silico study of action potential and QT interval shortening due to loss of inactivation of the cardiac rapid delayed rectifier potassium current

The rapid delayed rectifier K+ current, IKr, plays a key role in repolarisation of cardiac ventricular action potentials (APs). In recent years, a novel clinical condition denoted the short QT syndrome (SQTS) has been identified and, very recently, gain in function mutations in the gene encoding the pore-forming sub-unit of the IKr channel have been proposed to underlie SQTS in some patients. Here, computer simulations were used to investigate the effects of the selective loss of voltage-dependent inactivation of IKr upon ventricular APs and on the QT interval of the electrocardiogram. IKr and inactivation-deficient IKr were incorporated into Luo–Rudy ventricular AP models. Inactivation-deficient IKr produced AP shortening that was heterogeneous between endocardial, mid-myocardial, and epicardial ventricular cell models, irrespective of whether heterogeneity between these sub-regions was incorporated of slow delayed rectifier K+ current (IKs) alone, or of IKs together with that of transient outward K+ current. The selective loss of rectification of IKr did not augment transmural dispersion of AP repolarisation, as AP shortening was greater in mid-myocardial than in endo- or epicardial cell models. Simulated conduction through a 1 D transmural ventricular strand was altered by incorporation of inactivation-deficient IKr and the reconstructed QT interval was shortened. Collectively, these results substantiate the notion that selective loss of IKr inactivation produces a gain in IKr function that causes QT interval shortening.

Controlling alternans in cardiac cells.

Electrical alternans, the alternation in action potential morphology, has been suggested as an important cause of potentially dangerous cardiac rhythm disorders. Previous studies have developed alternans control strategies based on the dynamics of the relationship between action potential duration and the previous diastolic interval. We demonstrate that alternans in a single cardiac cell can also be controlled by directly modifying the underlying ion channel dynamics. Surprisingly, we find that, for a detailed canine ventricular cell model, the best time to apply the control stimulus is not during the repolarization phase of the action potential, but rather during the early plateau phase, when the charge requirements are two orders of magnitude smaller. Computer simulations show that a single control stimulus applied during the early plateau can completely eliminate small-amplitude alternans, while a small number of stimuli can rapidly extinguish large-amplitude alternans. We have also developed an effective control algorithm that uses only the membrane potential as control input and requires no prior detailed knowledge of the cell dynamics. The study suggests that control strategies based on ion channel dynamics can provide new directions for the development of algorithms intended to control dangerous cardiac rhythm disorders.

Intercellular coupling through gap junctions masks M cells in the human heart.

M cells have been described in many mammalian species. They are thought to be relevant for the genesis of long QT intervals, afterdepolarizations and for dispersion in action potential duration and in repolarization time. Their role in the human heart is subject to debate. We simulated action potential propagation in a strand of transversally oriented myocytes running from endocardium towards epicardium through the left ventricular free wall. The characteristics of the myocytes were either based on the Priebe-Beuckelmann ventricular cell model or on the Luo-Rudy ventricular cell model. The former model is based on the latter and includes adaptations in order to mimic the human ventricular myocyte. The amount and location of M cells as well as the intercellular coupling through gap junctions were varied. Also, we assessed action potential duration in a Langendorff-perfused explanted human heart and in a wedge preparation obtained from such a heart. At low, but physiological intercellular coupling conductance, the inclusion of M cells leads to a much longer 'QT interval' in the simulations than in the in vivo or isolated human heart. Dispersion in repolarization time becomes unphysiologically large when M cells are included in the strand and is also substantially larger than in the in vivo or isolated human heart. At stronger intercellular coupling this effect disappears. The manifestation of M cells is absent in the human heart, probably by effective intercellular coupling, turning them functionally "invisible".

A SIMPLE MODEL FOR INTERACTION OF VOLTAGE AND CALCIUM DYNAMICS IN VIRTUAL VENTRICULAR TISSUE

A simple two-variable model is used to replace the formulation of calcium dynamics in the Luo–Rudy ventricular cell model. Virtual ventricular cell and tissue are developed and validated to reproduce restitution properties and calcium-dependent voltage patterns present in the original model. Basic interactions between the membrane potential and Ca 2+ dynamics in the virtual cell and a strand of the virtual tissue are studied. Intracellular calcium waves can be linked to both action potentials (APs) and delayed afterdepolarizations (DADs). An intracellular calcium wave propagating from the cell interior can induce an AP upon reaching the cell membrane. The voltage and the intracellular Ca 2+ patterns within the same cell can be highly desynchronized. In a one-dimensional strand of the virtual tissue calcium motion is driven by the AP propagation. However, calcium release can be induced upon certain conditions (e.g. Na + overload of the cells), which results in DADs propagating in the wake of AP. Such propagating DADs can reach the excitation threshold, generating a pair of extrasystolic APs. Collision of a propagating AP with a site of elevated intracellular Ca 2+ concentration does not affect the propagation under the normal conditions. Under Na + overload local elevation of the intracellular Ca 2+ leads to generation of an extrasystolic AP, which destroys the original propagating AP.

Abstracts

Abstracts

Transition from ventricular fibrillation to ventricular tachycardia: a simulation study on the role of Ca2+-channel blockers in human ventricular tissue

We study the effect of blocking the L-type Ca2+-channel on fibrillation in simulations in two-dimensional (2D) isotropic sheets of ventricular tissue and in a three-dimensional anisotropic anatomical model of human ventricles, using a previously developed model of human ventricular cells. Ventricular fibrillation (VF) was obtained as a result of spiral wave breakup and consisted of a varying number of chaotically wandering wavelets activating tissue at a frequency of about 6.0 Hz. We show that blocking the Ca2+-current by 75% can convert ventricular fibrillation into a periodic regime with a small number of stable spiral waves, ranging from six in 2D sheets of 25 × 25 cm to a single spiral in the anatomical model of human ventricles. The dominant frequency during this process changed to about 10.0 Hz in the 2D simulations, but to only 5.0 Hz in the whole heart simulations where a single spiral wave anchored around an anatomical obstacle. We show that the observed effects were due to a flattening of the electrical restitution curve, which prevented the generation of wave breaks and stabilized the activation patterns.

Novel cardioprotective effects of tetrahydrobiopterin after anoxia and reoxygenation: Identifying cellular targets for pharmacologic manipulation

Contemporary cardioprotective strategies to prevent perioperative ischemia-reperfusion injury have focused on the l-arginine nitric oxide pathway. Tetrahydrobiopterin is an absolute cofactor required for the enzyme nitric oxide synthase and is thus a critical determinant of nitric oxide production. We hypothesized that ischemia-reperfusion results in diminished levels of tetrahydrobiopterin, which might represent a key cellular defect underlying endothelial and myocyte dysfunction after ischemia-reperfusion. To this aim, we examined the effects of tetrahydrobiopterin supplementation in (1) an in vivo experimental model of global ischemia-reperfusion and (2) an in vitro human ventricular heart cell model of simulated ischemia-reperfusion. Measures of endothelial function, oxidant production, cell survival, and cardiac function were used to assess outcome. In study 1 Wistar rats were divided into one of 2 groups (n = 10 per group). One group received tetrahydrobiopterin (25 mg x kg(-1) x d(-1) for 7 days), and the other group served as the control group. Hearts were subjected to 30 minutes of ischemia followed by 30 minutes of reperfusion, and left ventricular developed pressure, left ventricular systolic pressure, and left ventricular end-diastolic pressure were determined by using the modified Langendorff technique. In study 2 we quantitated myocardial malondialdehyde, a marker of lipid peroxidation, in ventricular tissues from both groups of animals using butanol phase extraction and spectrophotometric analysis. In study 3 coronary vascular responses were determined in vascular segments of the left coronary artery in both groups of animals after ischemia-reperfusion. Endothelium-dependent and endothelium-independent vasodilatation to acetylcholine and sodium nitroprusside, respectively, were compared between groups. In study 4, using a human ventricular heart cell model of simulated ischemia-reperfusion, we studied the effects of tetrahydrobiopterin (20 micromol/L) on cellular injury (as assessed by means of trypan blue uptake). After ischemia-reperfusion, myocardial dysfunction was evidenced by a decrease in left ventricular developed pressure and an increase in left ventricular end-diastolic pressure (P =.01 compared with baseline). Hearts from tetrahydrobiopterin-treated rats exhibited protection against ischemia-reperfusion injury (left ventricular developed pressure: 74 +/- 4 vs control 42 +/- 8 mm Hg, P =.01; left ventricular end-diastolic pressure: 12 +/- 3 vs 34 +/- 7 mm Hg, P =.01). Furthermore, tetrahydrobiopterin treatment attenuated the rise in malondialdehyde levels after ischemia-reperfusion (P =.01). After reperfusion, coronary endothelial function to acetylcholine was attenuated (P =.003 vs sham-treated mice), whereas responses to sodium nitroprusside remained unchanged. Tetrahydrobiopterin-treated rats exhibited an improvement in acetylcholine-mediated vasorelaxation (P =.01 vs ischemia-reperfusion group). Cellular injury, as assessed by means of trypan blue uptake, was higher in human ventricular heart cells subjected to simulated ischemia-reperfusion; this effect was prevented with tetrahydrobiopterin treatment (P =.001). Supplemental tetrahydrobiopterin provides a novel cardioprotective effect on left ventricular function, endothelial-vascular reactivity, oxidative damage, and cardiomyocyte injury after ischemia-reperfusion injury and might represent an important cellular target for future operative myocardial protection strategies.

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.

Hyperglycemia exaggerates ischemia-reperfusion-induced cardiomyocyte injury: reversal with endothelin antagonism.

We have previously demonstrated an importance of endothelin-1 in diabetic patients undergoing bypass surgery. Recent evidence suggests that cardiomyocytes might also produce endothelin-1, which might directly impair myocyte contractility by increasing intracellular calcium levels. Because hyperglycemia is a potent stimulus of endothelin-1 production, we hypothesized that increased production, action, or both of endothelin-1 might be a mediator of direct cardiomyocyte injury in diabetes. Therefore we studied the effects of endothelin receptor blockers (BQ-123 and bosentan) on hyperglycemia-induced endothelin-1 production and cellular injury after ischemia-reperfusion. Using a human ventricular heart cell model of simulated ischemia-reperfusion, we studied the effects of normoglycemia (5 mmol/L, 48 hours) and hyperglycemia (25 mmol/L, 48 hours) on cellular injury and endothelin-1 production. Furthermore, the effects of selective endothelin-A and mixed endothelin-A/B receptor antagonism (with BQ-123 and bosentan, respectively) were evaluated. Cellular injury, as assessed by means of trypan blue uptake, was higher in human ventricular heart cells subjected to hyperglycemia and simulated ischemia-reperfusion injury (P =.01); this effect was prevented with both BQ-123 and bosentan (P =.01). In addition, heart cells from the hyperglycemic group elaborated more endothelin-1 after ischemia-reperfusion (P =.02). Endothelin-1 production and cellular injury were greater in human ventricular heart cells subjected to hyperglycemic conditions and simulated ischemia-reperfusion. These effects are mediated by endothelin-A receptors because both BQ-123 and bosentan exerted similar degrees of protection. Endothelin receptor blockade is a novel strategy to improve the resistance of the diabetic heart to cardioplegic arrest and reperfusion.

L-arginine protects human heart cells from low-volume anoxia and reoxygenation.

Protective effects of L-arginine were evaluated in a human ventricular heart cell model of low-volume anoxia and reoxygenation independent of alternate cell types. Cell cultures were subjected to 90 min of low-volume anoxia and 30 min of reoxygenation. L-Arginine (0-5.0 mM) was administered during the preanoxic period or the reoxygenation phase. Nitric oxide (NO) production, NO synthase (NOS) activity, cGMP levels, and cellular injury were assessed. To evaluate the effects of the L-arginine on cell signaling, the effects of the NOS antagonist N(G)-nitro-L-arginine methyl ester, NO donor S-nitroso-N-acetyl-penicillamine, guanylate cyclase inhibitor methylene blue, cGMP analog 8-bromo-cGMP, and ATP-sensitive K+ channel antagonist glibenclamide were examined. Our data indicate that low-volume anoxia and reoxygenation increased NOS activity and facilitated the conversion of L-arginine to NO, which provided protection against cellular injury in a dose-dependent fashion. In addition, L-arginine cardioprotection was achieved by the activation of guanylate cyclase, leading to increased cGMP levels in human heart cells. This action involves a glibenclamide-sensitive, NO-cGMP-dependent pathway.

Influence of transmural repolarization gradients on the electrophysiology and pharmacology of ventricular myocardium. Cellular basis for the Brugada and long–QT syndromes

Ventricular myocardium comprises at least three electrophysiologically distinct cell types: epicardial, endocardial and M cells. Epicardial and M cells, but not endocardial cells, display action potentials with a notched or spike–and–dome morphology: the result of a prominent, transient, outward current–mediated phase 1. M cells are distinguished from endocardial and epicardial cells by the ability of their action potential to disproportionately prolong in response to a slowing down of rate and/or in response to agents with class III actions. This intrinsic electrical heterogeneity contributes to the inscription of the electrocardiogram (ECG) as well as to the development of a variety of cardiac arrhythmias. Heterogeneous response of the three cell types to pharmacological agents and/or pathophysiological states results in amplification of intrinsic electrical heterogeneities, thus providing a substrate as well as a trigger for the development of re–entrant arrhythmias, including Torsade de Pointes, commonly associated with the long–QT syndrome (LQTS), and the polymorphic ventricular tachycardia/ventricular fibrillation (VT/VF) encountered in the Brugada syndrome. Despite an abundance of experimental data describing the heterogeneity of cellular electrophysiology that exists across the ventricular wall, relatively few computer models have been developed to investigate the physiological and pathophysiological consequences of such electrical heterogeneity. As computer power increases and numerical algorithms improve, three–dimensional computer models of ventricular conduction that combine physiological membrane kinetics with realistic descriptions of myocardial structure and geometry will become more feasible. With sufficient detail and accuracy, these models should illuminate the complex mechanisms underlying the initiation and maintenance of Torsade de Pointes and other arrhythmias.

Models of cardiac ventricular action potentials: iterative interaction between experiment and simulation

The development of cardiac ventricular cell models since 1960 is reviewed with the focus on the interaction between simulation and experimental work. At each stage, the reasons for new models are e...

Modeling short-term interval-force relations in cardiac muscle.

This study employs two modeling approaches to investigate short-term interval-force relations. The first approach is to develop a low-order, discrete-time model of excitation-contraction coupling to determine which parameter combinations produce the degree of postextrasystolic potentiation seen experimentally. Potentiation is found to increase 1) for low recirculation fraction, 2) for high releasable fraction, i.e., the maximum fraction of Ca(2+) released from the sarcoplasmic reticulum (SR) given full restitution, and 3) for strong negative feedback of the SR release on sarcolemmal Ca(2+) influx. The second modeling approach is to develop a more detailed single ventricular cell model that simulates action potentials, Ca(2+)-handling mechanisms, and isometric force generation by the myofilaments. A slow transition from the adapted state of the ryanodine receptor produces a gradual recovery of the SR release and restitution behavior. For potentiation, a small extrasystolic release leaves more Ca(2+) in the SR but also increases the SR loading by two mechanisms: 1) less Ca(2+)-induced inactivation of L-type channels and 2) reduction of action potential height by residual activation of the time-dependent delayed rectifier K(+) current, which increases Ca(2+) influx. The cooperativity of the myofilaments amplifies the relatively small changes in the Ca(2+) transient amplitude to produce larger changes in isometric force. These findings suggest that short-term interval-force relations result mainly from the interplay of the ryanodine receptor adaptation and the SR Ca(2+) loading, with additional contributions from membrane currents and myofilament activation.

Protein kinase C-dependent activation of KATP channel enhances adenosine-induced cardioprotection.

Prior activation of protein kinase C (PKC) can precondition the cardiac cell against injury during subsequent ischaemia. By using cultured chick ventricular cell model for simulated ischaemia and preconditioning, the present study investigated the biochemical mechanism underlying the PKC-mediated preconditioning. A 5 min exposure to PMA enhanced the ability of pinacidil to mediate cardioprotection during a subsequent 90 min period of ischaemia, which is consistent with a sustained activation of the KATP channel initiated by PKC. The brief prior exposure to PMA was also associated with an enhanced ability of the adenosine A1 or A3 receptor agonist 2-chloro-N6-cyclopentyladenosine or N6-(3-iodobenzyl)adenosine-5'-N-methyluronamide to elicit a cardioprotective response during the subsequent ischaemia. In myocytes pretreated with PMA, the cardioprotection mediated by receptor agonist was blocked by the concomitant presence of KATP-channel antagonists glibenclamide or 5-hydroxydecanoic acid during the ischaemia. Thus the KATP channel acts downstream of the adenosine A1 and A3 receptors in mediating the protective effect due to prior PMA exposure. KATP channel activation is responsible for the adenosine receptor-mediated effect. PMA treatment had no effect on other A1 or A3 receptor-mediated effects such as the inhibition of adenylate cyclase, ruling out a direct stimulation of the receptor or G-protein by PMA. The present results indicate that prior stimulation of PKC causes a sustained KATP channel activation, which in turn renders the myocyte more responsive to the protective action of adenosine A1 and A3 receptor agonists during the subsequent ischaemia.

Cardiac Ca 2+ Dynamics: The Roles of Ryanodine Receptor Adaptation and Sarcoplasmic Reticulum Load

We construct a detailed mathematical model for Ca 2+ regulation in the ventricular myocyte that includes novel descriptions of subcellular mechanisms based on recent experimental findings: 1) the Keizer-Levine model for the ryanodine receptor (RyR), which displays adaptation at elevated Ca 2+; 2) a model for the L-type Ca 2+ channel that inactivates by mode switching; and 3) a restricted subspace into which the RyRs and L-type Ca 2+ channels empty and interact via Ca 2+. We add membrane currents from the Luo-Rudy Phase II ventricular cell model to our description of Ca 2+ handling to formulate a new model for ventricular action potentials and Ca 2+ regulation. The model can simulate Ca 2+ transients during an action potential similar to those seen experimentally. The subspace [Ca 2+] rises more rapidly and reaches a higher level (10–30 μM) than the bulk myoplasmic Ca 2+ (peak [Ca 2+] i ≈ 1 μM). Termination of sarcoplasmic reticulum (SR) Ca 2+ release is predominately due to emptying of the SR, but is influenced by RyR adaptation. Because force generation is roughly proportional to peak myoplasmic Ca 2+, we use [Ca 2+] i in the model to explore the effects of pacing rate on force generation. The model reproduces transitions seen in force generation due to changes in pacing that cannot be simulated by previous models. Simulation of such complex phenomena requires an interplay of both RyR adaptation and the degree of SR Ca 2+ loading. This model, therefore, shows improved behavior over existing models that lack detailed descriptions of subcellular Ca 2+ regulatory mechanisms.

Electrotonic modulation of electrical activity in rabbit atrioventricular node myocytes.

Electrotonic effects of electrically coupling atrioventricular (AV) nodal cells to each other and to real and passive models of atrial and ventricular cells were studied using a technique that does not require functional gap junctions. Membrane potential was measured in each cell using suction pipettes. Mutual entrainment of two spontaneously firing AV nodal cells was achieved with a junctional resistance (Rj) of 500 M omega, which corresponds to only 39 junctional channels, assuming a single-channel conductance of 50 pS. Coupling of AV nodal and atrial cells at Rj of 50 M omega caused hyperpolarization of the nodal cell, decreasing its action potential duration and either slowing or blocking diastolic depolarization in the AV node myocyte. Opposite changes occurred in the atrial action potential. When AV nodal and ventricular cells were coupled at Rj of 50 M omega, nodal diastolic potential was markedly hyperpolarized and diastolic depolarization was completely blocked with little change in ventricular diastolic potential. However, coupling did elicit marked changes in the action potential duration of both cells, with prolongation in the nodal cell and shortening in the ventricular cell. Nodal maximum upstroke velocity was increased by both atrial and ventricular coupling, as expected from the hyperpolarization that occurred. With an Rj of 50 M omega, spontaneous firing was blocked in all single AV nodal pacemaker cells during coupling to a real or passive model of an atrial or ventricular cell. These results demonstrate that action potential formation and waveform in a single AV nodal cell is significantly affected by electrical coupling to other myocytes.

Action potential conduction between a ventricular cell model and an isolated ventricular cell

We used the Luo and Rudy (LR) mathematical model of the guinea pig ventricular cell coupled to experimentally recorded guinea pig ventricular cells to investigate the effects of geometrical asymmetry on action potential propagation. The overall correspondence of the LR cell model with the recorded real cell action potentials was quite good, and the strength-duration curves for the real cells and for the LR model cell were in general correspondence. The experimental protocol allowed us to modify the effective size of either the simulation model or the real cell. 1) When we normalized real cell size to LR model cell size, required conductance for propagation between model cell and real cell was greater than that found for conduction between two LR model cells (5.4 nS), with a greater disparity when we stimulated the LR model cell (8.3 +/- 0.6 nS) than when we stimulated the real cell (7.0 +/- 0.2 nS). 2) Electrical loading of the action potential waveform was greater for real cell than for LR model cell even when real cell size was normalized to be equal to that of LR model cell. 3) When the size of the follower cell was doubled, required conductance for propagation was dramatically increased; but this increase was greatest for conduction from real cell to LR model cell, less for conduction from LR model cell to real cell, and least for conduction from LR model cell to LR model cell. The introduction of this "model clamp" technique allows testing of proposed membrane models of cardiac cells in terms of their source-sink behavior under conditions of extreme coupling by examining the symmetry of conduction of a cell pair composed of a model cell and a real cardiac cell. We have focused our experimental work with this technique on situations of extreme uncoupling that can lead to conduction block. In addition, the analysis of the geometrical factors that determine success or failure of conduction is important in the understanding of the process of discontinuous conduction, which occurs in myocardial infarction.