Am J Physiol Heart Circ Physiol 286: H1726-H1735, 2004.
First published December 23, 2003; doi:10.1152/ajpheart.00943.2003
0363-6135/04 $5.00
Action potential modulation of connexin40 gap junctional conductance
Xianming Lin and
Richard D. Veenstra
Department of Pharmacology, State University of New York Upstate Medical University, Syracuse, New York 13210
Submitted 2 October 2003
; accepted in final form 19 December 2003
 |
ABSTRACT
|
|---|
Connexin40 (Cx40) is abundantly expressed in the atrial myocardium, ventricular conduction system, and vascular endothelial and smooth muscle cells of the mammalian cardiovascular system. Rapid conduction through cardiac tissues depends on electrotonic transfer of the action potential between neighboring cells. To determine whether transjunctional voltages (Vj) elicited by an action potential can modulate conductance of Cx40 gap junctions, simulated myocardial action potentials were applied as voltage-clamp waveforms to Cx40 gap junctions expressed in mouse neuro2A (N2A) cells. Junctional currents resembled the action potential morphology but declined by >50% from peak to near-constant plateau values. Kinetics of Cx40 voltage gating were examined at peak voltages
100 mV, and decay time constants changed e-fold per 17.6 mV for Vj > ±40 mV. Junctional conductance recovered during phase 3 repolarization and early diastole to initial values. These phasic changes in junctional conductance were due to rapid decay kinetics, increasing to tens of milliseconds at peak Vj of 130 mV, and the increase in the steady-state conductance curve as Vj returned toward 0 mV. Time-dependent conductance curves for Cx40 were modeled with one inactivation and two recovery Vj-dependent components. There was a temporal correlation between development of conduction delay or block and the inactivation phase of junctional conductance. Likewise, recovery of junctional conductance was coincident with recovery from refractoriness, suggesting that gap junctions may play a role in the genesis and propagation of cardiac arrhythmias.
gap junction; ion channel; electrophysiology; transjuctional voltage
MYOCARDIAL GAP JUNCTIONS COMPOSED of connexin40 (Cx40) and connexin43 (Cx43) transmit electrical impulses via transient intercellular voltage gradients and are essential to action potential propagation. Cx40 is most abundantly expressed in the atrium and Purkinje fibers, along with Cx43, but is absent from the adult ventricular myocardium, where Cx43 exclusively predominates (6, 7, 34). Targeted gene knockout studies of Cx40 in mice produced prolonged P-R intervals and bundle branch block (19, 29). Most recently, familial atrial standstill was found to be associated with the coinheritance of genetic mutations in the cardiac sodium channel gene SCN5A and the Cx40 gap junction gene, which results in a functional decrease for both channel types (11). These observations attest to the importance of Cx40 to conduction through the atria and His-Purkinje system. Cx40 and Cx43 are also coexpressed in vascular endothelium and, to a lesser extent, in smooth muscle cells, where they are vital to the regulation of vasomotor tone via pharmacomechanical and electromechanical coupling signaling pathways (9, 20, 22).
Previously, we reported that the transjunctional voltage (Vj) gating properties of Cx43 are sufficient to produce phasic changes in junctional resistance (Rj) between two cells experiencing voltage gradients equivalent to the amplitude of the ventricular cardiac action potential (21). Because increasing Rj between cardiac myocytes can support slower conduction velocities (
) than reduced excitability, we wanted to determine whether the junctional conductance (gj) of Cx40 was similarly modulated by time-dependent Vj gradients (15, 27, 28). Intercellular conduction times can increase from 200 µs to >20 ms before complete conduction block develops (18, 36, 37, 41). Under these conditions, the safety factor for propagation depends increasingly on the L-type calcium current (25, 28, 31). Cx40 gap junctions have a half-inactivation voltage (V1/2) of
50 mV, 10 mV less than values reported for Cx43 (1, 33, 38, 39). Therefore, Cx40 gap junctions should experience similar, if not greater, changes in Rj during the depolarization and repolarization phases of the cardiac action potential relative to Cx43. In this study, we report our findings on the modulation of Cx40 gj by action potential-generated Vj gradients and develop a model for calculating the time-dependent gj curves at six different steady-state cycle lengths (CL) of stimulation.
 |
MATERIALS AND METHODS
|
|---|
Stable rat Cx40 transfectants of mouse neuro2A neuroblastoma (N2A) cells were cloned as previously described (1). Double whole cell patch-clamp recordings were achieved using the procedures defined for correcting for junctional series resistance errors (35). Electrode resistance was 631 M
(16.0 ± 6.4 M
) after membrane patch disruption, with whole cell input resistances of 15 G
and cell input capacitances of 1.53.0 pF. All experimental results included in the final analysis were limited to <5% error in the applied Vj for the entire duration of the protocol. All reported gj measurements represent corrected values, and junctional current (Ij) was defined as
I2 (35). All whole cell current recordings were obtained using previously described methods with dual Axopatch 200B (Axon Instruments, Union City, CA) or BioLogic RK-400 (Molecular Kinetics/Bio-Logic, Pullman, WA) patch-clamp amplifiers (21, 35). All records were low-pass filtered at 500 Hz (LPF-202A 4-pole Bessel filter, Warner Instruments, Sarasota, FL) and digitized at 4 kHz using a Digidata 1320 analog-to-digital board and pCLAMP 8.2 software (Axon Instruments). Data analysis was performed using the CLAMPfit program, and Ij and Vj calculations were performed offline. Graphs were constructed using Origin version 6.1 or 7.0 software (OriginLab, Northampton, MA).
The Luo-Rudy cardiac action potential voltage-clamp waveforms were the same as those previously described (21, 23), and stimulus CL of 250, 500, 750, 1,000, 1,500, and 2,000 ms were applied until a steady state was achieved (
200 simulated beats). The output voltage action potential waveform was applied to cell 1 (V1), whereas the command voltage for cell 2 (V2) was set constant to the simulated diastolic resting potential for each CL (88.1 mV
V2
90.2 mV).
 |
RESULTS
|
|---|
Time-dependent changes in Ij and gj during an action potential. Cx40 cell pairs were stimulated at a constant CL of 1,000 ms for 200 s as a control. The initial Ij response to the first applied action potential and the steady-state average of the last 100 action potentials are shown in Fig. 1A. There was a 12-ms delay before the onset of the action potential, and the action potential delay at 95% repolarization (APD95) was 153 ms at CL = 1,000 ms. The peak gj of this cell pair was 1.56 nS. The initial and average steady-state Ij traces qualitatively resemble the shape of the action potential but decline from the action potential peak to an almost constant steady-state value achieved near the midpoint of the action potential plateau. To determine the decline in gj during constant pacing at CL = 1,000 ms, Ij was divided by the corrected Vj value for the initial and average steady-state action potentials, and the time-dependent changes in gj during the train of 200 action potentials were plotted. Figure 1B illustrates that, for this example, gj declined by >50% from the peak of the action potential to a minimum quasi-steady-state plateau value. From this minimum plateau value, gj gradually began to recover during the final phase of repolarization. To average the results from all eight experiments at CL = 1,000 ms, gj was normalized (Gj) to the peak gj of the first action potential in each experiment, and the initial and average steady-state Gj curves were plotted as described above (Fig. 1C).

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 1. Connexin40 (Cx40) junctional current (Ij) and conductance (gj) responses to a cardiac action potential waveform. Ij signals (A) and gj calculations (B) from the 1st action potential (AP#1) and the ensemble average of the last 100 action potentials (AP#SS) from a train of 200 action potentials are illustrated for the action potential with cycle length (CL) = 1,000 ms. Action potential duration at 95% repolarization (APD95) was 153 ms, and there was a 12-ms delay before the onset of the action potential. Peak gj was 1.56 nS initially, declined to a quasi-steady-state value of 0.51 nS, and recovered toward initial values during action potential repolarization. C: average normalized junctional conductance (Gj) from 8 experiments. All gj values were normalized to the peak gj of the 1st action potential. Average steady-state Gj declined by 55% during the action potential before beginning its recovery phase.
|
|
At a frequency of 1 Hz, the average behavior of a Cx40 gap junction in response to constant pacing at 1 beat/s closely resembled the behavior of the single experiment shown in Fig. 1, A and B. During the first 25 ms of the action potential, Gj declined by 37% from the normalized peak Gj. Gj declined by another 18% during the plateau phase of the action potential to a nearly constant minimum of 0.45. From this minimum of 0.45, Gj rose gradually to 0.65 at APD95 for CL = 1,000 ms (153 mS; Fig. 1C). Gj recovered to its full initial value over the next 30 ms, 182 ms after the onset of the action potential.
The behavior of individual Cx40 gap junction channels was assessed from low-gj experiments, which permitted the resolution of quantal channel current fluctuations. The Ij recordings from one low-gj experiment (<0.6 nS) are illustrated in Fig. 2. Ten individual Ij traces, the ensemble average of those 10 traces, and the ensemble average of all 200 traces are displayed in Fig. 2A. The ensemble average of as few as 10 Ij traces, inclusive of all channel conductance states, reproduced the behavior of the higher-gj macroscopic recordings, where distinct channel current fluctuations were not observed because of the activity of
10 gap junction channels. A maximum of four 150-pS Cx40 gap junction channels were estimated to be present in this experiment (Fig. 2B). On average, approximately two 150-pS channels were open initially, with only one equivalent channel remaining open during the plateau phase of the action potential in this experiment. The activity of the inactivated channel began to recover during phase 3 repolarization. These data provide direct evidence for the Vj-dependent gating of individual Cx40 gap junction channels during the action potential and indicate that this is the primary mechanism for the modulation of Gj.

View larger version (47K):
[in this window]
[in a new window]
|
Fig. 2. Ensemble Cx40 single-channel activity. A: 10 Ij traces from a low-gj ( 0.6 nS) experiment, ensemble average of the same 10 traces, and ensemble average from all 200 Ij traces demonstrating high incidence of open channels at the beginning of the action potential with CL = 1,000 ms. B: gj calculations for traces in A. A maximum of four and an average of two 150-pS channels were initially open with an average of one channel closing during the action potential and remaining closed until final repolarization began during the cardiac action potential.
|
|
The same procedures were followed for the five other CL, and the results are summarized in Fig. 3A. APD95 for each CL was used to plot the steady-state Gj curve, because the Gj calculations became more variable as
, thus excluding the last 67 mV of repolarization from the final Gj analysis. Gj declined to a minimum of 0.450.50 relative to the peak Gj, except at CL = 250 ms, at which a minimum time-dependent Gj of only 0.60 was achieved, owing to the shorter action potential duration. In all cases, Gj increased toward initial peak values as Vj decreased from 85 mV toward 0 mV. Gj returned to the peak value at Vj
10 mV during the return phase of the CL-dependent Gj curve. The mean gj was 1.75 ± 1.33 (SD) for all experiments (n = 28).

View larger version (31K):
[in this window]
[in a new window]
|
Fig. 3. Time-dependent and steady-state Gj curves for Cx40. A: ensemble-averaged Gj over the last 100 action potentials from 58 experiments calculated relative to APD95 for each CL. Gj declined to a minimum plateau value of 0.450.60 and recovered to 0.701.00 of its initial value at ADP95 for each CL. B: steady-state Gj-junctional voltage (Vj) curve for Cx40 obtained from continuous Vj ramps from 0 to ±100 mV in 200 ms/mV increments. Solid curved line is best fit of a Boltzmann distribution (Eq. 1) to data from 4 experiments. Ensemble-averaged Gj curves from 5 Vj ramps per experiment were pooled to calculate mean Gj data points (symbols).
|
|
To determine the time-independent relation between Gj and Vj for Cx40, a steady-state Gj-Vj curve was obtained by pooling the ensemble average of five 200 ms/mV Vj ramps from 0 to ±100 mV from each of four different experiments (Fig. 3B). The average Gj-Vj curve was fitted with a Boltzmann equation as follows
 | (1) |
where Gmaxss is 1(Gj = gj/gj,max), i.e., the resting normalized slope conductance for each experiment; Gminss is the minimum value of gj/gj,max; A is the slope factor for the Boltzmann curve (zF/RT where F is Faraday's constant, R is the gas constant, and T is temperature); and V1/2 is the half-inactivation voltage. The slope factor is proportional to the gating charge movement (z) of the state transition (14, 30). For the curve shown in Fig. 3B that best describes the averaged data, Gminss = 0.24 and 0.22, V1/2 =52 and +48 mV, and z =2.9 and +3.7 elementary charges for negative and positive Vj values, respectively. These results are in close agreement with previous findings using this Vj protocol (1, 35). The minimum Gj achieved during the action potential plateau does not achieve the Gminss of the Cx40 steady-state Gj-Vj curve. The recovery phase of the CL-dependent Gj curves coincided with the increase in steady-state Gj as Vj decreased from +85 mV.
Voltage-dependent changes in Ij decay constants. The basis for the decline in gj during the early phases of the cardiac action potential to the quasi-steady-state minimum Gj was examined with Vj pulses between 40 and 140 mV to determine the decay constants at fixed Vj values. V1 was made negative, because there is a slow time- and voltage-dependent membrane conductance that activates above 20 mV in N2A cell membranes that can become quite large by the end of a 2.5-s Vj pulse (7.5 s for Vj
60 mV). A train of 5 or 10 Vj pulses of equal amplitude was applied at a rate of 1 pulse every 30 s, and the ensemble-averaged current was fitted with an exponentially decaying function. The decay time constants of selected Vj pulses from a single experiment are shown in Fig. 4A. Only the initial 400 ms of the 2.5-s pulses are displayed to better illustrate the rapid decay phase of the Ij signal recorded from cell 2. At some Vj values, there was a second slower decay phase with time constants on the order of 1 s that were not analyzed, because they amounted to <20% of the total decay in Ij and were insignificant within the time frame of the action potential. The decay time constants became progressively faster with increasing Vj, and the mathematical fits closely match the data in amplitude and time. The results from 411 experiments at the absolute value for each Vj (|Vj|) were averaged, and the reciprocals of the decay time constants are plotted in Fig. 4B. The >1,000-ms decay constants reported for Vj values near the V1/2 (i.e., 4060 mV) of the steady-state Gj-Vj curves of cardiac connexin gap junctions declined to
10 ms at potentials equivalent to the peak of the cardiac action potential. The curved line indicates the best fit of the data as follows
 | (2) |
The decay constants decline e-fold per 17.6 mV (=Vc, the voltage constant) from an initial amplitude (A0) of 0.0004 ms1 [time constant (
) = 2,500 ms] at Vj = 40 mV. This Vj-dependent decline in the Ij decay constants was determined only by examining Vj > ±100 mV. The increasingly rapid first-order decay kinetics indicate that the Vj-gating properties of Cx40 are fast enough to modulate gj during the time course of a cardiac action potential.
Time-dependent recovery of Ij and Gj. The novel observation that the CL-dependent Gj recovers during the late phases of the cardiac action potential was examined further by applying full-amplitude premature action potentials with incrementally increasing 10-ms delays to Cx40 gap junctions at CL = 1,000 ms. Figure 5A demonstrates that the Ij amplitude clearly increases with each additional delay. Each Ij trace from this single experiment represents the ensemble average of 20 sweeps for each stimulus delay. The premature stimulus delay was increased from 120 to 190 ms from the onset of the normal CL = 1,000 ms action potential. The average Gj from 5 or 10 experiments is indicative of a consistent increase in Gj from the minimum Gj of 0.450.50 obtained during the action potential to 85% of the initial Gj value during the first 30 ms of the diastolic interval (Fig. 5B). These observations confirm the increase in Gj observed in Figs. 1C and 3A.
To assess the time and Vj dependence of the Gj recovery phase, a pulse protocol that stepped from the plateau action potential Vj of approximately +80 to +70 mV
Vj
+10 mV was utilized. Results from this protocol revealed that Ij increased in a time-dependent manner with time constants of 1040 ms for Vj
+60 mV (Fig. 5C, Table 1). Small time-independent increases in Gj were also observed at Vj
+40 mV. There was no increase in Gj at +70 mV Vj. The time constants were poorly correlated with the Vj of the recovery pulse (r = 0.60 for linear or exponential fits). The time-dependent increase in Gj approached the values of the steady-state Cx40 Gj-Vj curve for each Vj (Fig. 5D). These data were consistent with the hypothesis that the rise in Gj resulted from a proportional Vj-dependent increase in steady-state open probability as observed in the steady state Gj-Vj curve (Eq. 1 and Fig. 3B). This occurred despite the convergence of the minimum Gj attained during the action potential to twice the Gmin of the steady-state Gj-Vj curve. The recovery time constants were
40 ms in the Vj
60-mV range, at least two orders of magnitude faster than the decay time constants at these same voltages. The rapid inactivation and recovery time constants at high and low Vj values reinforce the concept that the gating kinetics of gap junctions can achieve values in the range of 10 ms.
Modeling phasic Gj alterations. Figure 6 illustrates the temporal correlations between the inactivation of Gj during the large Vj gradients associated with the onset of the action potential and the recovery phases of Gj with the two phases of action potential repolarization. Projection of the average Gj curve at CL = 1,000 ms into the Vj plane clearly demonstrates that inactivation predominantly occurred at >100 mV and that Gj rose at Vj < 90 mV. If the changes in Gj during the cardiac action potential at each CL are due to the time-dependent kinetics and steady-state Gj-Vj relation of Cx40 gap junctions (Figs. 3B and 4B), it should be possible to model these phasic alterations in Gj under the same conditions. The time integral of Gj was calculated using the two-component model developed for Cx43 (21). The inactivation components for Cx40 were defined as follows
 | (3) |
where the two components described for Cx43 were set equal to each other for Cx40, because the decay of Ij was adequately described by a single-exponential decay (Fig. 4A). The initial conditions for each inactivation component were defined as follows
 | (4) |
so that
 | (5) |

View larger version (39K):
[in this window]
[in a new window]
|
Fig. 6. Functional correlation of Gj in time and voltage. A: temporal correlation between Gj and Vj for action potential at CL = 1,000 ms. Vertically hatched bar indicates the first 25 ms of an action potential during which conduction block develops; Cx40 Gj inactivates by 37% during this critical time period. Horizontally hatched bar indicates initial phase of Gj recovery during phase 3 repolarization, which produces a 15% increase in Gj from a minimum of 0.45. The diagonally hatched bar indicates the final phase of Gj recovery during final repolarization. Phases of Gj recovery also coincide with relative refractory and supernormal excitability periods of the cardiac action potential. B: autocorrelation plot of Gj and Vj best indicates voltages where inactivation, initial recovery, and final recovery phases of Gj occur. Inactivation is driven by high Vj values above +100 mV, and complete recovery requires repolarization to Vj < 10 mV.
|
|
The equivalent time-dependent inactivation components (Gt1 and Gt2) were each given a Gmin equal to one-half the Gmin of the steady-state Gj Boltzmann curve (Eq. 1 and Fig. 3B), and the Vj-dependent time constants were computed on the basis of the exponential decay time constants described in Fig. 4B (Eq. 2). Equation 2 provides an accurate description of the decline in Gj from initial resting values at each CL, as illustrated in Fig. 7. The dashed line in Fig. 7 represents the numerical solution to Eq. 3 for each CL. The exact Gmax, Gmin, and A0 values used to fit the decay phases of the Gj curve at each CL are listed in Table 2. The value of A0 was decreased from the estimated value of 0.0004 from Eq. 2 (Fig. 4B) for most CL owing to the different response times of the voltage-clamp recordings. The time lag between the model and Gj curves in Fig. 7, AF, reflects the difference between the model calculated from the real-time Vj values and the finite response time of the recorded Ij values.

View larger version (53K):
[in this window]
[in a new window]
|
Fig. 7. Model CL-dependent Gj curves. Average Gj curve (thin solid line) and calculated time-dependent changes in Gj according to Eq. 10 (model, thick solid line) for each CL are illustrated. Integration time step ( t) was 250 µs. Computed output of the 2 identical Vj-dependent inactivation components (Eq. 3; G1 & G2, thick dashed line), initial recovery (Eq. 6; R1, long-dotted line), and final recovery (Eq. 7, R2, short-dotted line) phases of Gj are also indicated. Mathematical model provides an accurate description of Cx40 data averaged from 58 experiments at each CL. Parameter values for fitted curves are provided in Tables 2 and 3.
|
|
Figure 6A further indicates that the two phases of Gj recovery coincided with the two rates of repolarization of the action potential during phase 3 and early diastole. The rates of repolarization for these two temporal phases are approximately 2.0 and 0.2 mV/ms, yet the majority of the increase in Gj occurs at Vj
10 mV. There is a time dependence to the increase in Gj at Vj
+60 mV, but the time constants varied only between 10 and 40 ms (Table 1 and Fig. 5C). Gj exceeded the steady-state Gj values at Vj
50 mV, because the minimum Gj achieved during the action potential was approximately twice the steady-state Gmin, and the recovery is due to the increasing open probability of only the inactivated channels. As shown previously for Cx43 (21), a dual-exponential function of voltage, related to the two phases of repolarization, adequately described the increase in Gj relative to Vj (Fig. 6B). The two components of the recovery phases of Gj were defined as follows
 | (6) |
and
 | (7) |
where
 | (8) |
and
 | (9) |
are the final values of Rt1 and Rt2. The actual values for Rmax1, Rmax2, AR1, and AR2 used to provide the fitted recovery curves in Fig. 7 are provided in Table 3. At CL = 250 ms, the contribution of Rt2 was minimal, whereas the amplitudes of Rt1 and Rt2 were relatively constant at the other CL. The voltage constant in Eq. 6 was changed from 24.5 to 32.0 mV for CL = 250 and 1,000 ms to provide a better fit of the data in Fig. 7, A and D. The numerical solutions to Eqs. 6 and 7 are illustrated by the short- and long-dotted lines in Fig. 7.
The Vj-dependent solutions to Eqs. 3, 4, 6, and 7 provide an accurate description of the phasic changes in Cx40 Gj observed during the cardiac action potential over frequencies of 30240 beats/min (Fig. 7). The final expression for Gj is
 | (10) |
which can be solved for continuously in time and readily lends itself to incorporation into existing models of cardiac action potential propagation. These results demonstrate the flexibility of the dual-component model to fit continuous time-dependent Gj curves as a function of Vj for a variety of CL and cardiac gap junctions comprising different homotypic connexin combinations.
 |
DISCUSSION
|
|---|
Our data demonstrate that GJ is in a non-steady state during the first 2040 ms of a cardiac action potential and declines with a time constant beginning in the 10-ms range until a quasi-steady-state GJ is achieved during the slowly changing plateau voltages. At this juncture, GJ increases in proportion to the steady-state GJ-VJ curve for Cx40 gap junctions. This results in a slow and steady increase in GJ as the VJ gradient decreases from 85 to 10 mV. At VJ < 10 mV, GJ returned to peak values. According to the original model for voltage gating of GJ, each gap junction channel has two gates in series that close with one specific voltage polarity (14, 30). As VJ increases from 0 mV with one polarity, the open probability will decrease toward a minimum residual conductance state (Gmin) on one side while remaining fully open on the other (Fig. 3B). As VJ declines, the closed gates reopen with a VJ-dependent rate constant (
) that is different from that for inactivation (
) (14, 30, 39) (Figs. 4 and 5 and Table 1). The modulation of GJ by the cardiac action potential occurs by similar mechanisms in Cx40- and Cx43-containing gap junctions (21). Although they are qualitatively similar, there are some quantitative differences in the rates of inactivation and recovery reported for these two cardiovascular connexins. The most significant difference is the monoexponential inactivation of Cx40 gap junctions relative to the biexponential decay observed for Cx43. There was also a difference in GJ recovery at +70 mV, where Cx40 exhibited no increase in GJ while Cx43 did. These differences are likely due to the modestly different VJ dependencies of Cx40 and Cx43 on VJ.
In the atrium, depolarization of an isolated atrial myocyte from 90 mV to the threshold for the sodium current of 60 mV will require
1.4 pC of charge. This calculation is based on the average atrial myocyte diameter of 15 µm, length of 100 µm, membrane surface area of 4,700 µm2, membrane resistivity of 2.63.0 k
·cm2, and specific membrane capacitance of 1 µF/cm2, all of which yield an estimated input resistance of 55 M
and membrane capacitance of 47 pF (2, 3, 16). With the assumption of a peak Vj of 130 mV, an Ij of 14 nA would deliver enough charge (Qj) to activate the cell in 100 µs. This requires a gj of 108 nS or an Rj of <9.3 M
. From Rj measurements between isolated pairs of adult ventricular myocytes, an Rj of 1.7 M
(gj = 590 nS) can be estimated (40, 41). There are no similar measurements for adult atrial myocyte pairs, but assuming a twofold reduction in cell surface area and a proportional increase in Rj provides a reasonable estimate. Hence, Rj = 4M
or gj = 250 nS, only a 2.5-fold excess of gap junction channels.
Values of
of 0.5 m/s in the atrial myocardium and 1.0 m/s in a Purkinje fiber translate into a conduction time of
100200 µs/cell. At these
values, up to 50% of the 200-µs conduction time per cell can be attributed to the junctional delay in activation (10). An action potential amplitude of 130 mV and maximum upstroke velocity (d
m/dt) of 100 V/s in the atrium indicates that the action potential upstroke has a duration of 1.3 ms, or approximately six cell lengths. Hence, the maximum Vj under normal conditions is only
20 mV, so the intercellular conduction delay and Vj gradient are not sufficient to induce Vj-dependent inactivation in well-coupled myocardium (4). However, this changes when
is reduced to
0.1 m/s. At
< 0.1 m/s, the intercellular conduction delay approaches the full duration of the action potential upstroke. This indicates that the slow conduction often associated with reentrant arrhythmias is associated with Vj gradients >100 mV and intercellular conduction delays >1 ms. It is under these conditions of decremental conduction that this newly described Vj-dependent gating mechanism can influence action potential propagation on a beat-to-beat basis. Furthermore, if, under normal conditions, each cell experiences only a fraction of the action potential amplitude during propagation (4), then the number of open gap junction channels required to deliver the 14 nA of excitatory current within 100 µs increases reciprocally. Hence, the number of open gap junction channels plays a more important role than previously suggested under normal conditions as well. Furthermore, recent evidence suggests that cardiac sodium channels are locally distributed near the Z lines and intercalated disks of the sarcolemma (43). This would serve to maximize the depolarization of the t tubule and intercalated disk membranes necessary to sustain excitation-contraction coupling and intercellular propagation.
Models of linear longitudinal propagation demonstrate that d
m/dt increases and then decreases as gj is reduced and that
varies linearly with the number of gap junction channels from >10,000 to <1,000 channels (5, 28). These models vary quantitatively in terms of the value of gj that must be achieved before
= 0.1 mm/ms, but there is general agreement that this occurs when <1,000 channels are present. This represents a >95% reduction in the estimated number of open gap junctions between ventricular myocytes (21) but only an estimated 40% reduction in the number of open gap junctions between atrial myocytes. Below
= 0.1 mm/ms,
begins to vary nonlinearly with gj (5, 28). Given that d
m/dt is not yet decreasing at the corresponding value of gj, this implies that the junctional delay has increased >10-fold, to >90% of the observed conduction time, for this 5-fold reduction in
. It is estimated that reductions in gj can support 200-fold reductions in
, whereas decreased excitability can produce only a 3-fold reduction in
before complete block develops (28). Hence, in regions of slowed conduction between 0.1 and 0.5 m/s, conduction delays of
1 ms occur between myocytes, Vj-dependent inactivation would be minimal, and the reduction in gj associated with the 10-fold increase in junctional conduction delay must occur by other regulatory mechanisms (13). However, further reductions in gj support extremely slow
and allow sufficient time for development of Vj-dependent inactivation. A doubling of the conduction time from this condition leads to a decreased safety factor for propagation and eventual conduction block (28).
Conduction delays of a few milliseconds normally occur in nodal tissues and between the Purkinje and ventricular cell layers of the ventricular subendocardial surface. These conduction delays were postulated to occur as a result of higher Rj values in these tissues (8, 17, 32). Conduction delays between ventricular cardiomyocytes can also increase to >20 ms before conduction fails completely (18, 36, 37, 41). Hence, conditions can exist where Vj-dependent changes in gj become prominent within regions of the mammalian heart, particularly in transitional border zones between tissues where physiological or pathophysiological differences in gj and membrane excitability occur. It is the purpose of this model to become incorporated into models of action potential propagation, where the effects of Vj-dependent activation on
can be examined at different levels of resting gj and various cellular properties.
The rapid decline in gj during early repolarization is especially relevant to slowed conduction, because it was demonstrated that phase 1 repolarization and the slow inward calcium current play critical roles in the propagation of the cardiac action potential under conditions resulting from high Rj values (25, 28, 31). The transient increase in gj that occurs during phase 3 repolarization is also significant to cardiac electrophysiology, because this phenomenon coincides with the vulnerable period of cardiac excitability when premature extrasystoles can arise (26, 27, 32). During this window in time and voltage, a single extrasystole can induce unidirectional conduction block, which elicits electrical disorganization and fibrillation. Therefore, in poorly coupled tissues with heterogeneous depolarization and repolarization patterns, changes in Rj may facilitate the occurrence and spread of a premature impulse in the time and space domains. This increases the probability of unidirectional conduction block and fibrillation in the ventricular conduction system and working myocardium.
Another possible explanation for the differences in the dynamic model for the gating of Cx40 and Cx43 gap junctions is the slightly higher whole cell access resistances of the Cx40 Ij recordings. The channel access resistance determines how the current is divided among the individual gap junction channels in a plaque. One account for the observed loss of Vj sensitivity with increasing gj was the postulated increase in access resistance as gap junction plaques become larger (42). This means that gap junction channels in the center of a gap junction plaque would contribute less current to the overall Ij. It is estimated that the normal-sized myocardial gap junction plaque would possess three times more access resistance if it existed as a uniformly packed 220-µm circular patch than if it took the form of a 20 nm x 2 µm strip (12). Nonuniform packing of gap junction channels within large-diameter macular gap junction plaques and reduction of plaque width to <0.6 µm serve to minimize the effect of access resistance and maximize conductance efficiency (12). These models indicate that the spatial distribution of gap junction channels within a plaque will influence a channel's access resistance (12, 42).
The loss of Vj-dependent gating in the dual whole cell configuration is unavoidably influenced by patch electrode resistance (35). To illustrate how patch electrode series resistance (Rel) affects the kinetics and steady-state Vj gating parameters, whole cell currents (I) were modeled as follows
 | (11) |
 | (12) |
where the input resistance (Rin) of each cell is equal to the parallel combination of the membrane (Rm) and gigaseal (Rs) resistances: Rin = (Rin·Rs)/(Rin + Rs). For this demonstration, Rin was kept constant at 1 G
by making Rs and Rm = 2 G
each. Actual experimentally obtained Rin values were slightly higher: 15 G
. The resting potential (Erest) of each cell was defined as the resting diastolic potential (89.77 mV) during the action potential at CL = 1,000 ms. Figure 8 illustrates the results obtained for Ij =
I2 under various experimental dual whole cell recording conditions during application of the action potential at CL = 1,000 ms. The dynamic gap junction model for Cx40 was added to the whole cell current model according to the parameters listed in Tables 2 and 3 for Eqs. 3, 6, and 7. Slowing of the inactivation kinetics is evident at gj = 1 nS as the resistance of each whole cell patch electrode increases from 0 (ideal) to 25 M
(upper limit of experimental recording conditions). Diminished inactivation is more pronounced under normal recording conditions (Rel1 = Rel2 = 10 M
) as gj increases from 1 to 40 nS (Rj = 1 G
25 M
). The operative factor for the reduced inactivation is the diminution of the applied Vj gradient in proportion to [Rj/(Rj + Rel1 + Rel2)]·(V1 V2), ignoring the contribution of the nonjunctional membrane currents to the total whole cell current (35). The decay time constants for Cx40 gap junctions will decrease e-fold for every 17.6-mV drop in applied Vj (Eq. 2). Our experimental results include data from low-conductance pairs consisting of a few Cx40 gap junction channels to cell pairs with gj
8 nS with minimal quantitative differences in the relative steady state to peak gj values (Figs. 2 and 3) (21, 35). The results with Cx40 also demonstrate the flexibility of the basic four-component Gj model to fit continuous time-dependent Gj curves as a function of Vj for a variety of CL and distinct connexin gap junctions by simple modification of the component amplitudes and voltage constants.
The importance of Cx40 to myocardial conduction is manifested in the conduction disturbances observed in the targeted Cx40 gene knockout studies in transgenic mice and in the recently reported incidence of familial atrial standstill (12, 19, 29). Additionally, the absence of Cx40 and/or Cx43 from vascular endothelium and smooth muscle impairs the regulation and conduction of vasomotor tone via coupled electrical and chemical signaling pathways (9, 20). Although tonic contractions are more prevalent in the vasculature than phasic contractions associated with action potentials, the possibility exists that long-duration intercellular gradients affect junctional communication by the same Vj-dependent mechanisms described here. It is also probable that long-duration Vj gradients >20 mV will produce significant reductions of gj by spermine block of Cx40 gap junction channels (24). In conclusion, these data and the dynamic model for the regulation of Cx40 gj by the ventricular cardiac action potential suggest that phasic changes in gj may contribute to the development of conduction block and fibrillation between poorly coupled excitable and inexcitable cells in a tissue. Implementation of these mathematical models for homotypic Cx40 and Cx43 gap junctions will allow for the future development of models related to various cardiac action potential waveforms, connexin composition of cardiac gap junctions, and propagation through cardiac tissues possessing different geometries and ionic currents.
 |
ACKNOWLEDGMENTS
|
|---|
We gratefully acknowledge the technical assistance of Mark Crye for maintenance of the Cx40-transfected N2a cell cultures. The simulated action potential waveforms used for voltage clamping were provided by Dr. Yoram Rudy's laboratory at Case Western Reserve University.
GRANTS
This study was supported by National Heart, Lung, and Blood Institute Grants HL-42220 and HL-45466 (to R. D. Veenstra).
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: R. D. Veenstra, Dept. of Pharmacology, SUNY Upstate Medical Univ., 750 E. Adams St., Syracuse, NY 13210 (E-mail: veenstrr{at}upstate.edu).
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
 |
REFERENCES
|
|---|
- Beblo DA, Wang HZ, Beyer EC, Westphale EM, and Veenstra RD. Unique conductance, gating, and selective permeability properties of gap junction channels formed by connexin40. Circ Res 77: 813822, 1995.[Abstract/Free Full Text]
- Bonke FIM. Passive electrical properties of atrial fibers of the rabbit heart. Pflügers Arch 339: 115, 1973.[CrossRef][ISI][Medline]
- Bonke FIM. Electrotonic spread in the sinoatrial node of the rabbit heart. Pflügers Arch 339: 1723, 1973.[CrossRef][ISI][Medline]
- Brink PR, Cronin K, and Ramanan SV. Gap junctions in excitable cells. J Bioenerg Biomembr 28: 351358, 1996.[CrossRef][ISI][Medline]
- Cole WC, Picone JB, and Sperelakis N. Gap junction uncoupling and discontinuous propagation in the heart. A comparison of experimental data with computer simulations. Biophys J 53: 809818, 1988.[Abstract/Free Full Text]
- Davis LM, Kanter L, Beyer EC, and Saffitz JE. Distinct gap junction phenotypes in cardiac tissues with disparate conduction properties. J Am Coll Cardiol 24: 11241132, 1994.[Abstract]
- Davis LM, Rodefeld ME, Green K, Beyer EC, and Saffitz JE. Gap junction protein phenotypes of the human heart and conduction system. J Cardiovasc Electrophysiol 6: 813822, 1995.[ISI][Medline]
- DeMello WC. Passive electrical properties of the atrio-ventricular node. Pflügers Arch 371: 135139, 1977.[CrossRef][ISI][Medline]
- De Wit C, Roos F, Bolz SS, Kirchhoff S, Kruger O, Willecke K, and Pohl U. Impaired conduction of vasodilation along arterioles in connexin40-deficient mice. Circ Res 86: 649655, 2000.[Abstract/Free Full Text]
- Fast VG and Kléber AG. Microscopic conduction in cultured strands of neonatal rat heart cells measured with voltage-sensitive dyes. Circ Res 73: 914925, 1993.[Abstract/Free Full Text]
- Groenewegen WA, Firouzi M, Bezzina CR, Vliex S, van Langen IM, Sandkuijl L, Smits JPP, Hulsbeek M, Rook MB, Jongsma HJ, and Wilde AAM. A cardiac sodium channel cosegregates with a rare connexin40 genotype in familial atrial standstill. Circ Res 92: 1422, 2003.[Abstract/Free Full Text]
- Hall JE and Gourdie RG. Spatial organization of cardiac gap junctions can affect access resistance. Microsc Res Tech 31: 446451, 1995.[CrossRef][ISI][Medline]
- Harris AL. Emerging issues of connexin channels: biophysics fills the gap. Q Rev Biophys 34: 325472, 2001.[ISI][Medline]
- Harris AL, Spray DC, and Bennett MVL. Kinetic properties of a voltage-dependent junctional conductance. J Gen Physiol 77: 95117, 1981.[Abstract/Free Full Text]
- Henriquez AP, Vogel R, Muller-Borer BJ, Henriquez CS, Weingart R, and Cascio W. Influence of dynamic gap junction resistance on impulse propagation in ventricular myocardium: a computer simulation study. Biophys J 81: 21122121, 2001.[Abstract/Free Full Text]
- Irisawa H. Electrophysiology of single cardiac cells. Jpn J Physiol 34: 375388, 1984.[ISI][Medline]
- Joyner RW, Overholt ED, Ramsa B, and Veenstra RD. Propagation through electrically coupled cells: two inhomogeneously coupled cardiac tissue layers. Am J Physiol Heart Circ Physiol 247: H596H609, 1984.[Abstract/Free Full Text]
- Joyner RW, Suguira H, and Tan RC. Unidirectional block between isolated rabbit ventricular cells coupled by a variable resistor. Biophys J 60: 10381045, 1991.[Abstract/Free Full Text]
- Kirchhoff S, Nelles E, Hagendorff A, Krüger O, Traub O, and Willecke K. Reduced cardiac conduction velocity and predisposition to arrhythmias in connexin40-deficient mice. Curr Biol 8: 299302, 1998.[CrossRef][ISI][Medline]
- Liao Y, Day KH, Damon DN, and Duling BR. Endothelial cell-specific knockout of connexin 43 causes hypotension and bradycardia in mice. Proc Natl Acad Sci USA 98: 99899994, 2001.[Abstract/Free Full Text]
- Lin X, Crye M, and Veenstra RD. Regulation of connexin43 gap junctional conductance by cardiac action potentials. Circ Res 93: E63E73, 2003.
- Little TL, Beyer EC, and Duling BR. Connexin-43 and connexin-40 gap junctional proteins are present in arteriolar smooth muscle and endothelium in vivo. Am J Physiol Heart Circ Physiol 268: H729H739, 1995.[Abstract/Free Full Text]
- Luo CH and Rudy Y. A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. Circ Res 74: 10711096, 1994.[Abstract/Free Full Text]
- Musa H and Veenstra RD. Voltage-dependent blockade of connexin40 gap junctions by spermine. Biophys J 84: 205219, 2003.[Abstract/Free Full Text]
- Rohr S and Kucera JP. Involvement of the calcium inward current in cardiac impulse propagation: induction of unidirectional conduction block by nifedipine and reversal by Bay K 8644. Biophys J 72: 754766, 1997.[ISI][Medline]
- Rudy Y. Reentry: insights from theoretical simulations in a fixed pathway. J Cardiovasc Electrophysiol 6: 294312, 1995.[ISI][Medline]
- Shaw RM and Rudy Y. The vulnerable window for unidirectional block in cardiac tissue: characterization and dependence on membrane excitability and intercellular coupling. J Cardiovasc Electrophysiol 6: 115131, 1995.[ISI][Medline]
- Shaw RM and Rudy Y. Ionic mechanisms of propagation in cardiac tissue. Roles of the sodium and L-type calcium currents during reduced excitability and decreased gap junction coupling. Circ Res 81: 727741, 1997.[Abstract/Free Full Text]
- Simon AM, Goodenough DA, and Paul DL. Mice lacking connexin40 have cardiac conduction abnormalities characteristic of atrioventricular block and bundle branch block. Curr Biol 8: 295298, 1998.[CrossRef][ISI][Medline]
- Spray DC, Harris AL, and Bennett MVL. Equilibrium properties of a voltage-dependent junctional conductance. J Gen Physiol 77: 7793, 1981.[Abstract/Free Full Text]
- Sugiura H and Joyner RW. Action potential conduction between guinea pig ventricular cells can be modulated by calcium current. Am J Physiol Heart Circ Physiol 263: H1591H1604, 1992.[Abstract/Free Full Text]
- Quan W and Rudy Y. Unidirectional block and reentry of cardiac excitation: a model study. Circ Res 66: 367382, 1990.[Abstract/Free Full Text]
- Valiunas V, Weingart R, and Brink PR. Formation of heterotypic gap junction channels by connexins 40 and 43. Circ Res 86: E42E49, 2000.[Abstract/Free Full Text]
- Van Kempen MJA, Ten Velde I, Wessels A, Oosthoek PW, Gros D, Jongsma HJ, Moorman AFM, and Lamers WH. Differential connexin distribution accommodates cardiac function in different species. Microsc Res Tech 31: 420436, 1995.[CrossRef][ISI][Medline]
- Veenstra RD. Voltage clamp limitations of dual whole-cell gap junction current and voltage recordings. I. Conductance measurements. Biophys J 80: 22312247, 2001.[Abstract/Free Full Text]
- Veenstra RD and DeHaan RL. Electrotonic interactions between aggregates of chick embryo cardiac pacemaker cells. Am J Physiol Heart Circ Physiol 250: H453H463, 1986.[Abstract/Free Full Text]
- Veenstra RD, Joyner RW, and Rawling DA. Purkinje and ventricular activation sequences of canine papillary muscle: effects of quinidine and calcium on the Purkinje-ventricular conduction delay. Circ Res 54: 500515, 1984.[Abstract/Free Full Text]
- Veenstra RD, Wang HZ, Westphale EM, and Beyer EC. Multiple connexins confer distinct regulatory and conductance properties of gap junctions in developing heart. Circ Res 71: 12771283, 1992.[Abstract/Free Full Text]
- Wang HZ, Li J, Lemanski LF, and Veenstra RD. Gating of mammalian cardiac gap junction channels by transjunctional voltage. Biophys J 63: 139151, 1992.[Abstract/Free Full Text]
- Weingart R. Electrical properties of the nexal membrane studied in rat ventricular cell pairs. J Physiol 370: 267284, 1986.[Abstract/Free Full Text]
- Weingart R and Maurer P. Action potential transfer in cell pairs isolated from adult rat and guinea pig ventricles. Circ Res 63: 7280, 1988.[Abstract/Free Full Text]
- Wilders R and Jongsma HJ. Limitations of the dual whole cell voltage clamp method in assaying conductance and kinetics of gap junction channels. Biophys J 63: 942953, 1992.[Abstract/Free Full Text]
- Zimmer T, Biskup C, Dugarmaa S, Vogel F, Steinbis M, Böhle T, Wu YS, Dumaine R, and Benndorf K. Functional expression of GFP-linked human heart sodium channel (hH1) and subcellular localization of the a subunit in HEK293 cells and dog cardiac myocytes. J Membr Biol 186: 112, 2002.[CrossRef][ISI][Medline]
This article has been cited by other articles:

|
 |

|
 |
 
X. Lin, J. Gemel, E. C. Beyer, and R. D. Veenstra
Dynamic model for ventricular junctional conductance during the cardiac action potential
Am J Physiol Heart Circ Physiol,
March 1, 2005;
288(3):
H1113 - H1123.
[Abstract]
[Full Text]
[PDF]
|
 |
|