Abstract
Hyperpolarizationactivated (pacemaker) channels are dually gated by negative voltage and intracellular cAMP. Kinetics of native cardiac fchannels are not compatible with HH gating, and require closed/open multistate models. We verified that members of the HCN channel family (mHCN1, hHCN2, hHCN4) also have properties not complying with HH gating, such as sigmoidal activation and deactivation, activation deviating from fixed power of an exponential, removal of activation “delay” by preconditioning hyperpolarization. Previous work on native channels has indicated that the shifting action of cAMP on the open probability (Po) curve can be accounted for by an allosteric model, whereby cAMP binds more favorably to open than closed channels. We therefore asked whether not only cAMPdependent, but also voltagedependent gating of hyperpolarizationactivated channels could be explained by an allosteric model. We hypothesized that HCN channels are tetramers and that each subunit comprises a voltage sensor moving between “reluctant” and “willing” states, whereas voltage sensors are independently gated by voltage, channel closed/open transitions occur allosterically. These hypotheses led to a multistate scheme comprising five open and five closed channel states. We estimated model rate constants by fitting first activation delay curves and single exponential time constant curves, and then individual activation/deactivation traces. By simply using different sets of rate constants, the model accounts for qualitative and quantitative aspects of voltage gating of all three HCN isoforms investigated, and allows an interpretation of the different kinetic properties of different isoforms. For example, faster kinetics of HCN1 relative to HCN2/HCN4 are attributable to higher HCN1 voltage sensors' rates and looser voltageindependent interactions between subunits in closed/open transitions. It also accounts for experimental evidence that reduction of sensors' positive charge leads to negative voltage shifts of Po curve, with little change of curve slope. HCN voltage gating thus involves two processes: voltage sensor gating and allosteric opening/closing.
Introduction
Hyperpolarizationactivated f/h channels, also referred to as “pacemaker” channels, have been described in a variety of cell types including cardiac and neuronal cells (DiFrancesco 1985, DiFrancesco 1993; Pape 1996). In cells where pacemaker channels are expressed, their properties serve different purposes, normally associated to the ability to generate a depolarization upon activation by previous hyperpolarization or by the second messenger cAMP.
Four isoforms of hyperpolarizationactivated, cyclic nucleotide–gated (HCN) channels have been recently cloned in several organisms (Santoro et al. 1997, Santoro et al. 1998; Gauss et al. 1998; Ludwig et al. 1998, Ludwig et al. 1999; Seifert et al. 1999; Vaccari et al. 1999). Expression experiments have demonstrated that HCN channels are the cloned equivalents of native f/h channels (Santoro et al. 1998; Ishii et al. 1999; Ludwig et al. 1999; Seifert et al. 1999; Moroni et al. 2000). HCN sequences are homologous to those of Kv and cyclic nucleotide–gated (CNG) channels (Clapham 1998), and belong to the same superfamily of six transmembrane (TM) domain channels.
The availability of different isoforms of HCN channels allows a detailed analysis of isoformspecific gating and modulatory properties, and a comparison of these properties with those of native channels as a way to corroborate expression data. Further, the existence of HCN isoforms with different gating properties prompts the question whether a comprehensive model can be developed to describe the kinetic features of the HCN family.
In this work we have analyzed detailed properties of voltage gating of three HCN clones (mHCN1, hHCN2, and hHCN4) based on the following considerations. First, it is known that kinetic features of native fchannel activation/deactivation cannot be described satisfactorily by the use of HodgkinHuxley (HH) description (DiFrancesco 1984). To account for these features, a complex kinetic model based on the coexistence of a “delaying” and a proper “gating” process and involving five voltagedependent gating variables of three different types had been previously proposed (DiFrancesco 1984). Second, voltage gating of several ion channels has been recently interpreted on the basis of allosteric models (Marks and Jones 1992; Rìos et al. 1993; McCormack et al. 1994; Horrigan and Aldrich 1999; Horrigan et al. 1999; Jones 1999). Third, the structural homology of HCN channels to 6TM domain K^{+} channels favors the hypothesis that HCN channels, too, are tetramers.
Our starting point was based on a recent investigation in cardiac SAN cells, showing that basic kinetic and modulatory properties of native pacemaker fchannels can be explained by a hybrid HHMWC model (DiFrancesco 1999). In that model, voltage gating was described by a standard HH scheme with two independent subunits, while cAMPinduced activation occurred according to a cyclic allosteric MonodWymanChangeaux (MWC) model (Monod et al. 1965) involving binding of one cAMP molecule per gating subunit.
On the basis of the above considerations, we explored the possibility that, as well as by cAMP, gating of pacemaker channels is modulated allosterically by voltage, too. In this attempt, we extended the hypothesis made in the gating scheme proposed earlier (DiFrancesco 1999) that pacemaker channels are composed of subunits, each independently gated by voltage and cAMP. We made the assumption that, like mSlo channels (Horrigan and Aldrich 1999; Horrigan et al. 1999), HCN channels are homotetramers, with each subunit carrying one independent voltage sensor, and that closed/open channel transitions involve concerted structural modifications of all four subunits and occur allosterically.
We found that by simply selecting an appropriate set of rate constants for each isoform, the allosteric assumption accounts for quantitative features of voltagedependent gating of different isoforms. The allosteric hypothesis provides an interpretation of the possible mechanisms responsible for the different gating properties of different isoforms, and may represent an unifying model for HCN channel gating.
Materials And Methods
Molecular Biology
Functional Expression of HCN Clones in Phoenix Cells
Mouse HCN1 cDNA (mHCN1) (construct pSD64TF, provided by Dr. B. Santoro, Columbia University, New York), was cut with EcoRV and NotI restriction enzymes, filled blunt and ligated into the eukaryotic expression vector pIRESEGFP (CLONTECH Laboratories, Inc.) predigested with EcoRV. Human HCN2 minigene (hHCN2) was subcloned into the eukaryotic expression vector pcDNA3.1 (Invitrogen) as previously described (Vaccari et al. 1999). Human HCN4 (hHCN4) was provided by Dr. U.B. Kaupp (Forschungszentrum, Jülich, Germany) in the eukaryotic expression vector pcDNA1 (CLONTECH Laboratories, Inc.). For transient functional expression channelspecific cDNA was transfected into modified HEK 293 (Phoenix) cells (Kinsella and Nolan 1996; Moroni et al. 2000). hHCN2 and hHCN4 constructs were cotransfected with a GFPcontaining plasmid. In brief, for each 35mm Petri dish we used 7 μg of channelspecific construct and 3.5 μg of GFPcontaining plasmid. In the case of mHCN1/pIRESEGFP, we used 10 μg of construct, which coexpresses the channel and the GFP protein on a polycystronic mRNA. Phoenix cells were cultured in Dulbecco's modified essential medium supplemented with 10% FCS (GIBCO BRL) and antibiotics (SigmaAldrich). In the following text, the species symbol will be dropped for simplicity and the mHCN1, hHCN2, hHCN4 isoforms will be termed HCN1, HCN2, and HCN4, respectively.
Electrophysiology
Experiments were performed on Phoenix cells incubated after transfection at 37°C in 5% CO_{2} for 1– 5 d. Before electrophysiological measurements, cells were dispersed by trypsinization and plated onto 35mm plastic Petri dishes at low density, where they were allowed to settle for 3–4 h. Petri dishes were transferred under the stage of an inverted microscope, and single cells were patchclamped in the wholecell configuration according to standard methods (Moroni et al. 2000). For the analysis, we visually selected strongly GFPreactive cells. Cells were perfused at room temperature (24–25°C) with a control Tyrode solution containing (in mM): 140 NaCl, 5.4 KCl, 1.8 CaCl_{2}, 1 MgCl_{2}, 5.5 dglucose, and 5 HEPESNaOH, pH 7.4. The solution was delivered via a fastflow perfusion pipet.
Wholecell pipets were filled with an intracellularlike solution containing (in mM): 10 NaCl, 130 KCl, 1.0 EGTA, 5 HEPESKOH, 0.5 MgCl_{2}, 2 ATP (sodium salt), 0.1 GTP (sodium salt), and 5 phosphocreatine, pH 7.2. We added 1 mM BaCl_{2}, 2 mM MnCl_{2}, 100 μM NiCl_{2}, and 20 μM nifedipine to perfusing solutions (in mM: 110 NaCl, 30 KCl, 1.8 CaCl_{2}, 0.5 MgCl_{2}, and 5 HEPESNaOH, pH 7.4) to reduce interference from overlapping Ca^{2+} and K^{+} components in the measurement of HCN currents.
Data Analysis
Activation curves were obtained by standard twostep voltageclamp protocols previously used (DiFrancesco et al. 1986; Moroni et al. 2000). The protocols consisted of hyperpolarizations applied in 10mV steps from a fixed holding potential of −35 mV, followed by steps to test voltages where the amplitude of deactivation current tails was measured. The duration of activation steps varied with the test voltage to allow attainment of steadystate.
Curve Fitting
For all HCN isoforms investigated, activation/deactivation traces were fitted by single exponential curves after an initial delay; delay, timeconstant, and current amplitude were free to vary in the bestfitting routine.
Numerical Computations
Analytical treatment of the allosteric scheme proposed (see Fig. 2) is given in the . Calculations were performed by integration of the vectorial differential Formula A2: where the probability vector p, the vector k, and the matrix A are defined in the . The integration was performed by a fourth order RungeKutta approximation method, with an integration step set to 1/1,500th of the total integration time. The channel open probability at any time t was calculated as Po(t) = p_{6}(t) + p_{7}(t) + p_{8}(t) + p_{9}(t) + p_{10}(t). In computations reproducing experimental current kinetics, the predicted current time course at a given voltage V was calculated as I(V,t) = I_{max}(V) · Po(V,t), where I_{max}(V) = g_{max} · (V − V_{rev}) was the fully activated current extrapolated from experimental data.
To reproduce individual experimental traces (Fig. 4Fig. 5Fig. 6) we used a semiempirical fitting routine. Computations were initially performed by assigning firstguess values to the selected parameters, chosen by previous identification of general criteria for restricting parameter values (see results) and then by trialanderror approximation, until satisfactory simulations were obtained.
Results
HCN Channels Have Kinetic Properties That Are Not Compatible with HH Models
Early investigation of the kinetics of native hyperpolarizationactivated channels (fchannels of cardiac Purkinje fibres) has identified features that are not compatible with HH gating models. These include a sigmoidal activation time course that cannot be described by a fixed power of an exponential at all voltages, sigmoidal deactivation and removal of activation “delay” by conditioning prehyperpolarizing steps (DiFrancesco 1984).
To probe whether HCN channels also have features not complying with HH kinetics, we analyzed the time course of current activation in cells expressing HCN1, HCN2, or HCN4, and checked if it could be described by a fixed power of an exponential. In Fig. 1 A, HCN1, HCN2, and HCN4 current activation traces recorded on hyperpolarization at voltages in the range −75 to −115 mV are shown. Plots of the function F(t) = ln(1− (I(t)/I_{∞})^{1/n}) should result in straight lines with zero yintercepts and slope = −1/τ if the current activates according to the nth power of an exponential. To verify this, we plotted in Fig. 1 B the function F(t) for each current record with n = 1, 2 and 3, and in Fig. 1 C, we plotted the yintercepts of fits to the linear part of F(t) plots.
Clearly, in none of the three examples presented can a fixed power of an exponential adequately describe the time course at all voltages. This is especially evident for HCN2. The voltage dependence is less marked for HCN4 and absent for HCN1, but a delay is apparent also for these two isoforms (see also Fig. 3). The analysis of Fig. 1 shows that, for all isoforms investigated, the activation time course can be described as a single exponential after a delay. The delay is small in HCN1, larger in HCN4 and most marked in HCN2 records.
Other gating properties of native fchannels in Purkinje fibres that cannot be reconciled with HH kinetics include sigmoidal deactivation and the removal of activation delay by conditioning prehyperpolarizing steps (DiFrancesco 1984). Both these features are also observed in HCN channels, and are dealt with later (Fig. 4Fig. 5Fig. 6 and Fig. 8). These observations indicate that, as in native fchannels, voltage gating of HCN channels is not compatible with HH models, and require more complex gating schemes involving multiple open and closed states.
Description of Kinetic Properties of HCN Channels by Allosteric Voltagedependent Gating
In Fig. 2, the reaction scheme used to describe allosteric voltagedependent gating of HCN channels is shown. Channels are thought to be composed of four subunits, each possessing one voltage sensor that can assume two configurations (one “willing” and one “reluctant” to open channels) and that is gated independently by voltage. With the allosteric hypothesis, all subunits undergo simultaneous, concerted transitions between closed/open states.
The analytical treatment, and the methods for numerical computations of the kinetics predicted by the model in Fig. 2 are described in . Under the assumptions made, the kinetic behavior of the allosteric model is fully characterized by knowledge of the parameter a and the eight parameters αo, βo, γo, δo, z_{α}, z_{β}, z_{γ}, and z_{δ} describing the voltage dependence of rate constants α, β, γ, δ according to Formula A9. A further simplification is obtained by setting z_{β} = −z_{α}, z_{γ} = −z_{δ} (see ). This reduces the overall free parameters to n = 7.
HCN Channels Activation and Deactivation Kinetics: Selection of Firstguess Fitting Parameters
Since the number of free parameters needed to simulate channel kinetics (n = 7) was too large to allow adequate fitting of data, we searched for criteria to restrict parameters values. The parameter “a” represents the multiplying factor of equilibrium constants of closed/open transitions upon movement of one voltage sensor from the reluctant to the willing state.
The model assumes that the transition of a voltage sensor to the willing state increases the opentoclosed channel ratio by a fixed factor, regardless of the channel subunit the voltage sensor belongs to. This implies that the parameter a is a constant in the range 0 < a < 1. Values of “a” close to 1 determine the condition L ≈ L1 ≈ L2 ≈ L3 ≈ L4, corresponding to practically ineffective voltage sensors; this would tend to abolish delays in activation and deactivation. Small values of “a” (a ≪ 1), on the other hand, determine the condition L ≫L1 ≫ L2 ≫ L3 ≫ L4. This would imply that most openings occur when all four voltage sensors are in the willing state (C4→O4), and that activation kinetics could be approximated by the linear scheme C→C1→C2→C3→C4→O4, with well separated “delaying” (C→C4) and “opening” (C4→O4) processes. Because of the symmetry of the cyclic allosteric scheme proposed, small values of “a” would also favor closing to occur mainly when all voltage sensors are in the reluctant state (O→C).
To test the predictive ability of the model, we tentatively assumed a = 0.2 and ran pilot computations. We found that, indeed, as expected for small values of a, the single exponential time course (after a delay) of computed activation traces approximated that of C4→O4 transitions (time constant τ4L; see Formula A10). We also noticed that plots of the delay displayed a voltage dependence similar to that of the time constant of voltage sensor transitions (τK(V); see Formula A10), reflecting the fact that, in the model of Fig. 2, the process responsible for delaying channel opening is the displacement of voltage sensors from reluctant to willing states (data not shown).
These preliminary observations provided a method for a rough estimate of some of the required parameters. We fitted experimental traces by single exponentials after an initial delay. Plots of time constant and delay data, thus, measured for HCN1, HCN2, and HCN4 channels are shown in Fig. 3. We bestfitted the time constant data with the equation for τ4L(V) (see Formula A10; Fig. 3 A); this yielded the functions α4(V) and β4(V), from which firstguess parameters αo, βo, and z_{β} (= −z_{α}) could be calculated as reported in the legend to Fig. 3.
Delay data (Fig. 3 B) were also bestfitted with the equation for τK(V) (see Formula A10); this generated the functions γ(V) and δ(V), from which the parameters γo, δo, and z_{δ} (= −z_{γ}) were calculated as reported in Fig. 3 legend. Table 1 lists the fitting parameters that were used as an initial set for predicting the kinetic behavior of the three HCN isoforms investigated.
To verify if the parameters in Table 1 could represent a starting set for subsequent fitting of kinetic data, we selected representative experimental current traces recorded during activation/deactivation voltageclamp protocols, and used the parameter values in Table 1 to run model computations with the same protocols. We found that despite the elementary criteria used to select rate parameters, the computations appeared to succeed in reproducing qualitatively the main kinetic features of the three HCN isoforms, including the voltage dependence of rates of activation and deactivation and the different degrees of activation and deactivation delay for different isoforms, although, not surprisingly, they failed to reproduce accurately the time course of individual traces. To verify how the parameter a affected computed channel kinetics, prediction tests were also performed using a = 0.1, 0.3, and 0.5, and were less satisfactory (data not shown). Therefore, we selected a = 0.2.
Fitting of Activation and Deactivation Time Course for HCN1, HCN2, and HCN4 Channels
To verify if the model was able to more accurately mimic the kinetics of current activation and deactivation of HCN channels in specific cases, we used semiempirical fitting routines on individual sets of traces as shown in Fig. 4Fig. 5Fig. 6 for HCN1, HCN2, and HCN4 channels, respectively.
For each trace, we initially used Table 1 to assign starting values to selected parameters (y4, z, τ4L, and τK) and calculated the rate constants α, β, γ, and δ according to Formula A10 and Formula A6 of . These values (with a = 0.2) allowed calculation of the matrix A and of vector k (Formula A7b and Formula A8), and thus of the current time course. The process was repeated by trialanderror approximation, until a satisfactory simulation was obtained.
For all isoforms, the model was able to reproduce the detailed kinetic features of activation and deactivation in the whole voltage range investigated. The data points selected with the above procedure for y4, z, τ4L, and τK deviated from the firstguess values assumed on the basis of the preliminary analysis of Fig. 3 (compare closed circles with broken lines in Fig. 4Fig. 5Fig. 6, panels C–F). They were bestfitted using the equations for y4, z, τ4L, and τK (see Formula A9, and Formula A10) and the resulting rate parameters are reported in Table 2. These can be directly compared with those of Table 1.
Steadystate Open Probability
Steadystate solution of the allosteric model in Fig. 2 yields the open probabilityvoltage relation Po(V) expressed by Formula A13, which clearly does not represent a Boltzmann type of voltage dependence. Nonetheless, to verify if it can describe the steadystate open probability curve of HCN channels, we used it to fit in Fig. 7 the mean activation curves for HCN1, HCN2, and HCN4.
For all three isoforms, the allosteric Po relation was able to fit adequately the steadystate activation curves. Po values obtained independently by fitting activation/deactivation kinetics in Fig. 4Fig. 5Fig. 6 (open circles) fell near the Po curves (full lines) obtained with Formula A13, indicating that the two procedures led to comparable predictions.
Removal of the Delay in Current Activation by Conditioning Prehyperpolarization
In cardiac Purkinje fibres, activation delay can be removed by conditioning prehyperpolarizations (DiFrancesco 1984). To verify whether HCN channels have the same property, we applied twostep protocols to cells expressing HCN4 channels. As apparent in the example of Fig. 8 A, short prehyperpolarizations to −150 mV accelerated activation at −110 mV and progressively removed the delay until monoexponential time course was achieved with a 22ms prepulse (see also semilog plots in Fig. 8 B) monoexponential time course was achieved.
To verify the ability of the allosteric model to mimic this behavior, we fitted records at −110 mV and at −150 mV with the procedure used above (Fig. 6), and then ran twostep model computations reproducing experimental protocols (Fig. 8 C). The model accurately predicted a progressive current acceleration and removal of delay with increasing prepulse durations. With sufficiently long prepulses, the computed time course became a single exponential, as apparent from the semilog plot of Fig. 8 D. This is expected if the action of prehyperpolarizing steps is to “prime” the channel opening process by pushing the majority of voltage sensors to their willing state, such that at −110 mV, the current development reflects mostly C4→O4 transitions. The computed traces at −110 mV, however, did not perfectly overlap experimental ones particularly after the longest prepulses. This may be due to the fact that, especially with long prepulses, the experimental current at −110 mV underwent an early acceleration compared with a single exponential (also observed in native fchannel; DiFrancesco and Ferroni 1983; DiFrancesco 1984), as evident in the semilog plots of Fig. 8 D (46ms prepulse record). The model constraints do not allow to reproduce this behavior.
Discussion
Early investigation of native pacemaker fchannels has shown that description of detailed kinetic properties requires a complex multistate, nonHH model (DiFrancesco 1984). Experimental evidence presented here indicates that some features of the voltagedependent gating of various HCN channel isoforms (mHCN1, hHCN2, and hHCN4) also do not comply with HH kinetics. These features include the following: (1) sigmoidal activation that cannot be approximated by the nth power of an exponential; (2) sigmoidal deactivation; and (3) removal of the activation delay by prehyperpolarizing steps. This indicates that a multistate nonHH kinetic model may be required for HCN channels.
HCN channels belong to the superfamily of Kv and CNG channels. The amino acid sequence comprises six transmembrane (TM) domains with a pore region between S5 and S6 and a positively charged domain S4, which acts as the primary voltage sensor (Vaca et al. 2000; Chen et al. 2000). The structural similarities to 6TM domain channels favor the hypothesis that HCN channel subunits coassemble as tetramers, and that their voltagedependent closed/open transitions are governed by the voltage sensitivity of S4 movement.
A gating model taking into account the molecular information available and the structural similarities with Kv and CNG channels should, therefore, regard HCN channels as tetramers whose closed/open transitions are regulated by voltage sensors. Two additional observations help clarify the features of this model: first, previous work on native fchannels has shown that the cAMPinduced modulation can be described by an allosteric model (DiFrancesco 1999); second, allosteric voltagedependent gating has been proposed for other channels, including Ca^{2+} channels (Marks and Jones 1992), Shaker K^{+} channel (McCormack et al. 1994), and large conductance, Ca^{2+}activated (mSlo) K^{+} channels (Horrigan and Aldrich 1999; Horrigan et al. 1999).
Allosteric Model of HCN Channels Voltage Gating and Comparison with Previous Models
Based on the above considerations, we asked whether the voltage gating of HCN channels could be described as an allosteric process, which led to a 10state cyclic reaction scheme (Fig. 2). A major problem that we faced in using the model was to restrict the choice of free model parameters in a way that would allow the selection of parameter values by data fitting, since clearly too many free parameters would imply the need of arbitrary selection of some of the values.
One assumption of the model at one time intuitive and useful to the effect of limiting free parameters, was to consider that the action on closed/open channel ratio of any voltage sensor's displacement is always the same, regardless of the distribution of voltage sensors (i.e., the number of sensors in the willing state versus those in the reluctant state), or to which channel subunit the voltage sensor belongs. This meant that the parameter a could be considered as a constant (in the range 0–1). This assumption, together with the assumption of balanced rate constant changes (i.e., when the equilibrium constant is multiplied by a, then the opening rate constant is divided by √a, and the closing rate constant is multiplied by √a), was very effective in that it reduced the set of 10 originally free “vertical” rate constants in Fig. 2 (α, β, α1, β1, α2, β2, α3, β3, α4, and β4) to just three free rate constants (the parameter a and α, β or, equivalently, α4 and β4).
On this basis, after a preliminary selection of model parameters (Table 1), we used individual curve fitting to estimate rate constants (Fig. 4Fig. 5Fig. 6 and Table 2). In all three HCN isoforms investigated, the allosteric scheme of Fig. 2 could reproduce the observed kinetic behavior, including sigmoidal activation and deactivation (Fig. 4Fig. 5Fig. 6) and removal of activation delay by prehyperpolarizations (Fig. 8). Also, the model accounted for the position and shape of the open probability curve for the three isoforms (Fig. 7).
The allosteric scheme examined here is analogous to that proposed by Marks and Jones 1992 for dihydropyridinesensitive Ca^{2+} channels, by Rìos et al. 1993 for the Ca^{2+}release channel of sarcoplasmic reticulum of skeletal muscle, by McCormack et al. 1994 for Shaker K^{+} channels, and by Aldrich and collaborators for large conductance Ca^{2+}activated K^{+} (mSlo) channels (Horrigan and Aldrich 1999; Horrigan et al. 1999). A distinct property of mSlo channels is the presence of a region of shallow slope in the time constant curve at high negative voltages; according to the allosteric hypothesis, this can be interpreted as due to a weak voltage dependence of closed/open transitions, which are slower than voltage sensor transitions and ratelimiting at large negative voltages (Horrigan et al. 1999).
In our model of HCN channels, closed/open transitions are indeed severalfold slower than sensor transitions for all isoforms (see Fig. 10); however, fitting of the time constant curves in Fig. 3 A and curve reconstruction in Fig. 4Fig. 5Fig. 6 indicate that closed/open transitions have a substantial voltage dependence (z_{β} = −z_{α} = 0.94, 1.07, and 0.80 equivalent electronic charges for HCN1, HCN2, and HCN4 isoform, respectively, Table 2).
It is important to stress that the values of z_{α} and z_{β} reported in Table were not assumed, but resulted from data fitting, and, therefore, cannot be reduced arbitrarily. In other words, the assumption that each voltage sensor contributes equally to favor channel opening, together with the other model assumptions, leads to experimental data fitting indicating the presence of a substantial voltage dependence for closed/open transitions for all isoforms.
Investigating a model where closed/open transitions are weakly voltagedependent would imply a different set of assumptions, since most of the voltage dependence of kinetics would derive from voltage sensors' state changes. The operating framework of the model would be likely to change significantly. We would need to release the assumptions concerning the constancy of a and/or the balanced rate constant changes, and let most of the 10 vertical rate constants α, β, α1, β1, α2, β2, α3, β3, α4, and β4 free to vary. In fact, the situation could be even more complex since releasing the condition a = constant would also affect the horizontal rate constants.
We did not adopt this solution because it would have involved too much arbitrariness in the selection of parameter values. More importantly, we did not have any real indication that closed/open transitions are weakly voltagedependent or voltageindependent for HCN channels, since we did not have evidence for regions of shallow slope at negative voltages in time constant curves of any of the HCN isoforms investigated (Fig. 3 A). However, the model scheme of Fig. 2 is not incompatible with a situation where closed/open transitions are nearly voltageindependent (i.e., with a low z_{α}) and most of the voltage dependence of channel kinetics is attributable to sensor movements. Although we lack experimental evidence, we cannot rule out this possibility entirely, since, due to experimental constraints, the voltage range where recordings could be made was limited (positive to −135 mV; Fig. 3 A).
Implications of Allosteric Design
Recent data have shown that replacement of basic residues in the S4 segment of HCN2 with neutral ones leads to negative shifts of the activation curve voltage dependence (Chen et al. 2000; Vaca et al. 2000). However, the slope of the activation curve did not decrease as expected if it were correlated to the total number of equivalent charges moving through the electric field during activation/deactivation. The allosteric model may provide a simple explanation for these results, as illustrated in the example of Fig. 9 for HCN2.
In Fig. 9 A, the sensor distribution parameter z (reluctant/willing configurations) was assumed to undergo changes due to reduction of the equivalent charge number (z_{δ}) by 1/6, 2/6, 3/6, and 4/6 of the control value (labeled 1–4). As expected, reductions of z_{δ} led to proportional reductions in the slope of the z(V) curve. However, in Fig. 9 B, the Po(V) relations calculated with the different z(V) curves in Fig. 9 A were shifted to more negative voltages (by up to 46.9 mV with a foursixths reduction) with only a minor change in slope. Thus, a possible explanation for the lack of modification of the slope of the activation curve in experiments where charged residues are mutated is that the activation process can be separated, as in the allosteric model proposed here, into two distinct steps (i.e., redistribution of voltage sensors and closed/open transitions). A direct consequence is that the slope of the Po curve does not necessarily reflect the slope of gating charge voltage dependence. Since the model provides a molecular interpretation of gating properties of HCN channels, it will be interesting to use it to simulate expectations of experiments in which structural elements of HCN channels are modified (i.e., mutations and/or coexpression of different isoforms).
Model Interpretation of the Different Kinetic Features of HCN Isoforms
The kinetics and cAMP dependence of different HCN isoforms differ quantitatively. For example, HCN1 subunits form channels with fast voltagedependent gating and reduced cAMP sensitivity, whereas HCN2 and HCN4 channels respond more efficiently to cAMP and have slower kinetics (Santoro et al. 1998; Seifert et al. 1999; Ishii et al. 1999; Ludwig et al. 1999; this work). The allosteric scheme proposed accounts for quantitative differences in gating properties among HCN isoforms within the framework of a common reaction model. To help understand how the different kinetic properties of HCN isoforms can be interpreted according to the allosteric scheme, we superimposed in Fig. 10 plots of the rate parameters obtained for the three isoforms investigated by the fitting procedures in Fig. 4Fig. 5Fig. 6 ([closed circles] HCN1; [open circles] HCN2; [closed triangles] HCN4) along with their bestfitting curves. These plots are representative of both kinetic and steadystate channel properties according to the allosteric model (see superimposition of Po values in Fig. 7).
The comparison led to the following observations: (1) the y4(V) curves of all three isoforms are almost superimposable; (2) the voltage dependence of τ4L(V) curves is similar for all isoforms, the main difference being a scalingup factor when going from HCN1 to HCN2 and HCN4 curves; (3) the z(V) curves of HCN1 and HCN4 are little different, whereas that of HCN2 is some 40 mV more negative; and (4) a similarly large voltage shift is present for the τK(V) curve of HCN2 relative to HCN1 and HCN4 curves; the peak τK increases when going from HCN1 to HCN2 and to HCN4.
The above observations allow a few tentative speculations on the mechanisms responsible for the kinetics of HCN channels and their differences. First, a scalingup of time constant curve that leaves unaltered the distribution curve implies in a Boltzmann mechanism that on and offrates are divided by the same factor (>1), equivalent to adding a fixed positive energy term to the energy changes associated to on and offreactions. The data in Fig. 10 (A and B), thus, would be compatible with the assumption that closed/open transitions of HCN channels involve two mechanisms: (1) a voltagedependent one, similar for all isoforms and involving movement of an equivalent electric charge of ∼1; and (2) a voltageindependent “mechanical” one perhaps associated to intersubunit interactions that are loose in HCN1, tighter in HCN2, and tighter still in HCN4 channels. Assuming simple scalingup of time constant curves, the τ4L curve of HCN1 is ∼5 and 15 times smaller than those of HCN2 and HCN4, respectively. Compared with HCN1, the extra energies required by HCN2 and HCN4 closed/open transitions, thus, would be RT ln(5) = 0.956 kcal/mol and RT ln(15) = 1.609 kcal/mol (see Formula A9 and Formula A10).
Second, according to the data in Fig. 10, faster activation of HCN1 relative to HCN2 and HCN4 is attributable to both faster voltage sensor activation and channel opening, whereas the more negative position of the HCN2 activation curve relative to that of HCN1 or HCN4 (Fig. 7) essentially depends on the voltage dependence of voltage sensors. Third, transition rates of voltage sensors are several times faster than closed/open transition rates, implying that for all isoforms, the activation process involves early sensor movements (from reluctant to willing states) as a preliminary requirement for efficient channel opening. Finally, lack of major changes in the equivalent charge associated to movement of voltage sensors and closed/open transitions suggests common structural rearrangements responsible for voltagedependent state changes in the isoforms studied.
Acknowledgments
We thank Antonio Malgaroli, Richard B. Robinson and Arnaldo Ferroni for discussion, and Monica Beltrame for assistance in development of the molecular biology of HCN clones. We also wish to thank B. Santoro and S. Siegelbaum (Columbia University, New York, NY) for supplying mHCN1 cDNA, and U.B. Kaupp of the University of Jülich for supplying hHCN4 cDNA.
This work was supported by the MURST (Cofin 1999 to D. DiFrancesco) and by Telethon (Grant no 971 to D. DiFrancesco).
Appendix
Kinetics of Allosteric Model
A general set of equations governing the kinetics of the scheme in 2 is FD1 where p_{i} is the probability of state i, k_{i,j} is the rate constant for the transition i to j, and n is the number of possible states (DiFrancesco 1984). The system can be reduced by one degree by use of the relation ∑p_{i} = 1 and expressed in matrix notation as A2 where the vectors p and k and the matrix A are defined as: A3 A4 A5
To apply the general Eq. 5 to the scheme of Fig. 2, it is useful to first define the rate constants k_{i,j}. As illustrated in the legend to Fig. 2, the following relations hold for equilibrium constants and rate constants in the reaction scheme (note that a = 1/f^{2}): A6
Using the numbering of states 1–10 as labeled in Fig. 2 and the relations A6, the rate constants k_{i,j} can be tabulated as in Table A1.
We can now explicitly write the elements of Eq. 2 as follows: A7a A7b
Although integration of Eq. A2 is possible (DiFrancesco 1984), the complexity of the system makes an analytical solution too difficult and requires numerical integration. The system is fully described once the parameter f (or a = 1/f^{2}) and the rate constants α, β, γ and δ are known. The voltage dependence of rate constants can be described by Boltzmann functions: A8 A9 where αo, βo, γo, and δo are rate constants at zero voltage, z_{α}, z_{β}, z_{γ}, and z_{δ} are the equivalent units of electronic charge moving in the electric field during corresponding state transitions, and r = RT/F = 25.85 mV (T = 300°K). It is assumed that z_{β} = −z_{α} and z_{γ} = −z_{δ}. The former assumption is based on the experimental observation that activation and deactivation rate constants of native pacemaker channels have a nearly symmetrical dependence on voltage (Noble et al. 1989; DiFrancesco 1999); the latter considers that the charge moved across the electric field upon transition of a voltage sensor from a reluctant to a willing state is equal and opposite to that moved in the reverse reaction. Numerical description of the kinetics of Fig. 2 scheme, thus, depends on knowledge of the six parameters αo, βo, γo, δo, z_{β}, and z_{δ} and of the parameter a.
An equivalent description, useful to a more intuitive interpretation of the meaning of parameters, can be obtained by using fractional distribution parameters and time constants. One such set of parameters equivalent to Formula A9 is for example: A10
(note that α4 = α/a^{2} and β4 = β a^{2} according to A6). Here y4 and z represent the fractional distribution between channel states O4 and C4 and between sensor willing and reluctant states, and τ4L and τK represent the corresponding time constants of equilibrium changes.
An explicit solution of A2 describing the probability of open state can be obtained by the MWC formulation of cyclic allosteric reactions (Monod et al. 1965). At equilibrium: A11
Summing over open and closed states and rearranging yields total open ([Ot]) and total closed ([Ct]) concentrations: A12 from which the probability of open state is derived as A13
The voltage dependence of Po results from the voltage dependence of rate constants according to Formula A9 and the definition of equilibrium constants L and K in the scheme of Fig. 2: A14 where A15 are the midpoint voltages of the fractional distribution parameters y = 1/(1+L) (distribution between O/C states) and z = 1/(1+K) (distribution between willing/reluctant sensor states). As mentioned above (Eq. A10) an equivalent description can be obtained by referring to C4↔O4, rather than C↔O transitions. In this case, the variables used are (see Eq. A6): A16 A17 and from A13 A18
Footnotes

Abbreviations used in this paper: CNG, cyclic nucleotide–gated; HCN, hyperpolarizationactivated, cyclic nucleotide–gated; HH, HodgkinHuxley; MWC, MonodWhymanChangeaux; Po, open probability; SAN, sinoatrial node; TM, transmembrane.
 Submitted: 19 December 2000
 Revision requested 9 April 2001
 Accepted: 13 April 2001