Probability fluxes and transition paths in a Markovian model describing complex subunit cooperativity in HCN2 channels.

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are voltage-gated tetrameric cation channels that generate electrical rhythmicity in neurons and cardiomyocytes. Activation can be enhanced by the binding of adenosine-3',5'-cyclic monophosphate (cAMP) to an intracellular cyclic nucleotide binding domain. Based on previously determined rate constants for a complex Markovian model describing the gating of homotetrameric HCN2 channels, we analyzed probability fluxes within this model, including unidirectional probability fluxes and the probability flux along transition paths. The time-dependent probability fluxes quantify the contributions of all 13 transitions of the model to channel activation. The binding of the first, third and fourth ligand evoked robust channel opening whereas the binding of the second ligand obstructed channel opening similar to the empty channel. Analysis of the net probability fluxes in terms of the transition path theory revealed pronounced hysteresis for channel activation and deactivation. These results provide quantitative insight into the complex interaction of the four structurally equal subunits, leading to non-equality in their function.

Introduction

In a protein the time scale of conformational changes ranges from picoseconds for the thermal vibration of the atoms to seconds, or even more, for structural arrangements [1]. Since proteins consist typically of hundreds to several thousand amino acids and each amino acid consists of seven or more atoms, the complete space of energetically possible states is huge, and certainly too huge for the experimentalist to study. It turned out, however, that often a set of a great many of these states can be viewed as metastable. This means that a protein typically fluctuates within a set of multiple different structures for a very long time if compared to the picosecond time scale of atom vibrations. Only extremely rarely the thermal energy suffices to leave a metastable set of states to reach another metastable set of states. A striking proof for this concept has been the discovery of the digital behavior of ion channels: Their pore is either closed or open, with dwell times often in the millisecond or even second time range [2]–[4].

Over the past decade there has been considerable progress in studying protein function theoretically by molecular dynamics (MD) simulations [5]. However, the maximum time range of presently possible MD simulations is still too short to become useful for the description of ion channel function. Kinetic analyses of ion channel gating are therefore presently conducted in terms of Markovian models, thereby treating a metastable set of states as one state and specifying transitions between two of these states according to Eyring's transition state theory [6]. Often, however, the accuracy of the experimental data does not suffice to differentiate between plausible models of sufficient complexity, e.g. to describe the action of three or more subunits. As a consequence the models used are often implausibly small or many simplifying assumptions for the equilibrium or rate constants have to be adopted.

Markovian models frequently used to describe the action of ion channels are either of the sequential or the cyclic type. Among the cyclic models the Monod-Wyman-Changeux (MWC) model [7] is of outstanding relevance because of its simple elegance: It assumes that a protein exists in two global conformations, taut and relaxed, and that each binding step systematically shifts the taut-relaxed isomerization and increases the binding affinity. This model, which has been most widely used to describe the sigmoidal oxygen-binding curve to hemoglobin (for review see [8]), has also been applied to describe the gating of ion channels by ligands, as e.g. the nicotinic acetylcholine receptor [9], cyclic nucleotide-gated (CNG) channels [10] or hyperpolarization activated cyclic nucleotide-modulated (HCN) channels [11]. However, as elegant as this model is, it does neither describe all aspects of the oxygen binding to hemoglobin [12], [13] nor is there any evidence that it fully describes the action of any ion channel [14]. For nicotinic acetylcholine and glycine receptors, the analysis of single-channel recordings has allowed the investigators to determine models with substantially more free parameters than typically employed by the MWC model. As a consequence, also the type of cooperativity differed substantially [15]–[18]. However, in all these approaches the ligand binding affinity could only be determined indirectly by the global fit of the channel activity, i.e. by the pore action, but not by direct measurement.

For CNGA2 and HCN2 channels we recently developed a strategy to simultaneously measure ligand binding and activation gating by combining confocal microscopy with patch-clamp fluorometry [19], thereby employing a fluorescently labeled ligand [20], [21]. In a subsequent study on HCN2 channels, pre-activated by a voltage pulse to −130 mV, we combined this approach with the method of concentration jumps. By globally fitting multiple time courses of ligand binding and channel activation, this approach allowed us to determine the equilibrium and rate constants in a Markovian model with 4 binding steps in both the closed and the open channel and 5 closed-open isomerizations (Fig. 1) [22]. The analysis revealed pronounced cooperativity with respect to the microscopic binding affinity in the surprising sequence ‘positive – negative – positive’ for the binding of the second, third, and fourth ligand, respectively. Moreover, we considered the population of all closed and open states as function of time when jumping the ligand concentration. As a result, the total open probability is dominated by open states with either zero, two or four ligands bound whereas states with one or three ligands bound are only transiently populated (Fig. 1) [22].

Figure 1

The C4L-O4L model.

(A) The C4L-O4L model is composed of five closed (C x) and five open states (O x; x = 0…4). Ligand (L) binding is possible in both the closed and the open channel. The closed-open isomerization can proceed at each degree of liganding. The rate constants and their errors were determined previously [22]. k O3C3 was set equal to k O4C4. Since for these rate constants only a lower border could be estimated, they were set to 3 s−1 herein if not otherwise noted. This results in k C3O3 = k C4O4 = 2.0×102 s.1. (B) Time-dependent probability to stay in a closed state, P C,x; (x = 0…4) for a ligand jump (green bar) from zero to 7.5 µM and back to zero. P c denotes the sum of all P C,x at each time. (C) Time-dependent probability to stay in an open state, P O,x; (x = 0…4) for a ligand jump (green bar) from zero to 7.5 µM and back to zero. P o, the sum of all P O,x at each time, indicates the measured time course of the total open probability. (D) Scheme of the C4L-O4L model with means of rate constants provided by Table S1.

Herein we analyze the transition pathways in the C4L-O4L model [23], [24] when disturbing an equilibrium by a sudden change of a parameter. In our case this sudden change is either the application or the removal of a saturating ligand concentration. As a result the probability fluxes as function of time and the net probability fluxes are quantified, thereby identifying relevant and irrelevant transition pathways in channel activation and deactivation.

Results

Net probability flux densities in the C4L-O4L model

To specify the probability flux for any transition, we first consider a simple model with the states A and B and the two rate constants, k 1 and k −1, specifying the transition rates. A and B can also be interpreted as probabilities (with A+B = 1) which change under non-equilibrium conditions with time. At any time, a net probability flux density, f, can be defined by the sum of two opposed unidirectional fluxes. For example, the net probability flux density from A to B is given by f AB = k 1 A−k −1 B. If applying this to the individual transitions in the C4L-O4L model (for rate constants see Fig. 1, Table S1), the net probability flux density of all transitions can be calculated at each time. In case of the binding reactions the rate constants have to be multiplied by the actual ligand concentration (either 0 or 7.5 µM). The arithmetic signs were chosen such that a probability flux is positive for the binding and opening reactions and, accordingly, negative for the unbinding and closing reactions. Because the probability of the states is time-dependent, the net probability flux density of each transition is also time-dependent.

First, the net probability flux density is considered when stepping from zero to 7.5 µM fcAMP. For the four binding steps in the closed channel, the net probability flux density is given by L is the ligand concentration. As expected, the net probability flux density moves like a wave from C 0 to C 4 (Fig. 2, left): f C0C1 ceases after less than 100 ms and f C1C2 after about 500 ms. f C2C3 and f C3C4 are slower and not finished after 1s yet. The net probability flux density from the closed states to the respective open states contributes to the reduction of the net probability flux density between the closed states.

Figure 2

Time courses of net probability flux density in the C4L-O4L model.

The time courses for the net probability flux density, f XY, in the C4L-O4L model were computed for the event of applying the ligand fcAMP (7.5 µM; left diagrams) and removing it again (right diagrams). The inset diagrams show the original diagram with an either higher amplitude or time resolution. (A) Net probability flux densities between closed states following ligand application (left) and removal (right) obtained by equation (1). (B) Net probability flux densities between open states following ligand application (left) and removal (right) obtained by equation (2). (C) Net probability flux density of the closed-open isomerizations following ligand application (left) and removal (right) obtained by equation (3).

Accordingly, the net probability flux density for the four binding steps of the open channel, that has been pre-activated by voltage, is given by

The respective time courses of f Ox−1Ox are basically similar to those of the closed states but larger in amplitude (Fig. 2, left). For the closed-open isomerizations the net probability flux density is given by

These time courses have a robust amplitude for f C1O1, f C3O3, and f C4O4 and are negligible for f C0O0 and f C2O2 (Fig. 2, left). The obstruction of the pathway C 2→O 2 tells that the double liganded channel has a similarly taut structure as the non-liganded channel. The reduced degree of determinateness for f C3O3, and f C4O4 will be considered below.

Second, the net probability flux density of stepping from 7.5 µM fcAMP back to zero is considered. The respective time courses for the unbinding steps in the closed and open channel are given by equation (1) and (2), respectively, by setting L = 0 (Fig. 2, right). The time courses for the open-closed isomerizations are given by equation (3) (Fig. 2, right). The result is that the net probability flux density is much bigger in the unbinding steps of the open than of the closed channel and that f O3O4 and f O2O3 are much faster than f O1O2 and f O0O1, i.e. the state O 2 is a severe obstacle in the net probability flux on the way to lower liganded states. Concerning the open-closed isomerizations, f O4C4 is negligible and f O3C3 is rapidly finished after about 1 s. In contrast, the other net probability fluxes are much slower. However, they all notably contribute to the closing process despite the rate constants for closing, k O2C2, k O1C1, and k O1C1, differ substantially. This shows how important it is to consider the model as a whole, including the rate constants and the time-dependent population of the states.

Total net probability fluxes and transition paths in the C4L-O4L model

From the net probability flux densities, f XY, one can easily compute the total net probability fluxes, F XY, as the time integral over the time interval from the concentration jump (t = 0 s) to an end time, t end, t end is either 5 s for the fcAMP pulse or 20 s after removal of fcAMP. Adapted to our C4L-O4L model, F XY indicates how much of the total net probability flux moves along a transition X↔Y. The main results are: (1) The total net probability flux is bigger in the binding steps of the open than the closed channel in both the presence of fcAMP and after its removal (Fig. 3). (2) The total net probability flux between the closed states is larger for the binding than for the unbinding process (Fig. 3). (3) F C0O0 is negligible in the binding-induced relaxation but significantly present in the unbinding-induced relaxation whereas F C4O4 is negligible in the unbinding-induced relaxation but significantly present in the binding-induced relaxation (Fig. 3). Together, these results suggest different pathways for activation and deactivation.

Figure 3

Total net probability fluxes in the individual transitions of the C4L-O4L model.

The total net probability fluxes per transition, F XY (given as fraction of unity), were obtained by integrating (equation (4)) the time courses of the net probability flux densities obtained by equation (1) to (3). The data after stepping to 7.5 µM fcAMP and after its removal are shown on the left and right side, respectively. The dotted bars indicate the effect of fitting with k O3C3 to 30, instead of 3, on F C3C4, F O3O4, F C3O3, and F C4C4. The sum of F C3O3 and F C4C4 is practically unchanged. The other total net probability fluxes were unchanged. (A) Total net probability fluxes for the transitions between neighbored closed states. (B) Total net probability fluxes for the transitions between neighbored open states. (C) Total net probability fluxes for the five closed-open transitions.

To gain a more thorough insight into the transition pathways in our C4L-O4L model, they were analyzed according to the principles of the transition path theory [23], [25]. In this type of analysis one is interested in trajectories between two selected states within the model. This implies that all repeated jumps between two neighbored states in both directions are not counted; i.e. considered are only transitions in one direction. This information is hidden in the total net probability flux for the individual transitions, F XY, which we determined in the previous section. In any model the amount of net flux produced by one set of states must equal the amount of net flux collected by the other set of states. When applying the ligand, the only flux producers are the states C 0 and O 0. Because the applied ligand concentration was saturating, the main flux collector state is O 4 while O 3 and O 2 also contribute but to an only small extent (Fig. 1). After removing the ligand, the main flux producer state is O 4, and to a small extent O 3 and O 2. The main flux collector states are C 0 and O 0. The states C 2 and O 2 are, to a minor extent, also flux collector states because after 20 s the steady state was not reached (Fig. 1).

In case of a fully reversible Markovian model, which holds for the C4L-O4L model, the total probability flux can easily be decomposed into the pathway net probability fluxes, and no cycles remain because the condition of the detailed balance is fulfilled [22]. The procedure is to choose first a pathway from state X 0 to state X k along the states X i. Then for this chosen pathway, X 0→X k, the net probability flux F p,X0Xk is computed according to F p,X0X is then removed from the flux along all edges of this pathway and the procedure is repeated for the remaining pathway net probability fluxes until the total possible probability flux between these states is zero [25]. It should be noted that the sequence of such a decomposition is not unique because different paths can be chosen. It has become useful to start with the strongest pathway [25].

Let us first consider the pathway net probability flux after applying the ligand at 7.5 µM fcAMP (Fig. 4). The pathway net probability flux from C 0, one of the two flux producers, to O 4, the main flux collector, results in five defined pathway fluxes, F p,, which are, according to equation (5), given by the respective minimum values of F CxOx (x = 0…4). The pathway net probability flux from O 0, the other flux producer, to O 4 has only one available path because the fluxes of all closed-open transitions are directed to the respective open states. Illustration of the weights of the pathway net probability fluxes in the C4L-O4L model after applying the ligand leads to the result that ligand-induced channel opening proceeds mainly from C 3 and C 4 (with the uncertainty of the exact attribution; Fig. 3, left), and additionally from C 1, but not relevantly from C 0 and C 2. The same type of analysis was then performed for the pathway net probability flux after removing the ligand from O 4, the main flux producer, to C 0 and O 0, the main flux collectors. Closing of the channels proceeds predominantly from O 1 and O 0, and to a minor extent from O 2 and O 3, but not relevantly from O 4. Together these results show that there is pronounced hysteresis for ligand-induced activation and deactivation. The pathway flux from O 4 to O 0, the other flux collector, has again only one available path because the fluxes of all closed-open transitions are directed to the respective closed states.

Figure 4

Transition pathways in the C4L-O4L model at saturating fcAMP.

The weight of the pathway net probability fluxes, obtained by equation (5), is given by the numbers besides the colored graphs (as fraction of unity) and, for values >0.030, also encoded by the thickness of the arrows. Fluxes from 0.002 to 0.030 are illustrated by equally thick dotted arrows. All other fluxes are only very small (<4.5×10−4) and therefore not represented by arrows. For further explanation see text. (A) Pathway net probability fluxes along the pathways C 0 to O 4 (red) and O 0 to O 4 (blue) after switching to 7.5 µM fcAMP. (B) Pathway net probability fluxes along the pathways O 4 to C 0 (red) and O 4 to O 0 (blue) after removing fcAMP.

To demonstrate the function of our model at subsaturating fcAMP concentrations we considered the fluxes to the two main collectors at 0.75 µM fcAMP (O 2 and O 4; Fig. 5 A,B) and 0.075 µM fcAMP (O 1 and O 2; Fig. 5 C,D) and the respective reverse fluxes when removing fcAMP. Notably, these selected fluxes are only the dominating fluxes. The results show that at the intermediate concentration of 0.75 µM fcAMP activation proceeds along C 1→O 1 predominantly to O 4 (Fig. 5 A) and to a minor extent to O 2 (Fig. 5 B). In addition there is a remarkably big net probability flux in the open channel along O 0→O 2 and in the reverse direction along O 2→O 0 which is absent along O 0→O 4 and O 4→O 0, respectively. This indicates that O 2 is a metastable state. Moreover, this result corresponds to the high energy barrier for the transition O 2→O 3 [22]. At 0.075 µM fcAMP the predominant activation proceeds along C 1→O 1 to O 2 but not anymore to O 4. O 1 is only passed in the activation pathway. However, it is of importance as a collector for the probability flux in the open channel (Fig. 5 D).

Figure 5

Transition pathways in the C4L-O4L model at subsaturating fcAMP.

The illustration is analogous to Fig. 4. The two main net probability fluxes at each concentration are shown. Pathway net probability fluxes associated with activation (C 0→O x, O 0→O x) and deactivation (C x→O 0, O x→O 0) are shown in red and blue color, respectively. (A) 0.75 µM fcAMP, C 0↔O 4. (B) 0.75 µM fcAMP, net probability flux along C 0↔O 2 in addition to that included in C 0↔O 4. (C) 0.075 µM fcAMP, C 0↔O 2. (D) 0.075 µM fcAMP, net probability flux along C 0↔O 1 in addition to that included in C 0↔O 2.

Unidirectional fluxes in the closed-open isomerization as function of time

In our considerations on probability fluxes so far the focus was set on net probability fluxes. If a flux is reversible, as in case of ligand binding and closed-open isomerizations, a net flux consists of two opposed unidirectional fluxes. Generally, very different opposed unidirectional fluxes can cause the same net flux. Large opposed unidirectional fluxes indicate conformational flexibility between two states, whereas small opposed unidirectional fluxes indicate conformational tautness between two states.

For the principal conformational change of the closed-open isomerization we considered the time courses of the unidirectional probability flux densities, f U,CxOx>0 and f U,OxCx<0 (x = 0…4), at the saturating fcAMP concentration of 7.5 µM fcAMP and related them to the net probability flux density, f CxOx (Fig. 6; c.f. Fig. 2), according toFor the time after application of fcAMP, these time courses show that the transitions C 0↔O 0 and C 2↔O 2 are functionally irrelevant. In contrast, the net probability flux densities in the transitions C 1↔O 1 and C 3↔O 3 are robust and the unidirectional probability flux densities exceed the respective net probability flux density substantially. An extreme surplus of the unidirectional probability flux densities with respect to the net probability flux density is inherent in the transition C 4↔O 4, suggesting pronounced conformational flexibility when the channel is fully liganded. It is also notable that the empty activated channel (Fig. 6, top left) is much less flexible in this sense compared to the fully liganded channel (Fig. 6, bottom left).

Figure 6

Unidirectional probability flux densities in the closed-open isomerizations.

The time courses of the unidirectional probability flux densities, f U,CxOx and f U,OxCx (x = 0…4) for the closed-open isomerizations are plotted together with the respective net probability flux densities, f CxOx, when applying 7.5 µM fcAMP (left) and removing it (right). The arrows indicate the time point of concentration change. For explanation see text.

After removal of fcAMP (Fig. 6, right), the fully liganded state is left rapidly (bottom) and the triple liganded state is rapidly passed. For the remaining transitions C 2↔O 2, C 1↔O 1, and C 0↔O 0 both the maximum net and the maximum unidirectional probability flux densities are much smaller than for C 4↔O 4 and C 3↔O 3. The finally large total net probability flux in the transitions C 2↔O 2, C 1↔O 1, and C 0↔O 0 (c.f. Fig. 3, right) is only reached after many seconds. Remarkably, in the transition C 2↔O 2 the unidirectional probability flux density in the opening direction approximates zero, which is due to the low occupancy of C 2, resulting in a close similarity of the unidirectional closing and the net probability flux density. This contrasts to C 1↔O 1 and C 0↔O 0 which both show larger unidirectional probability flux densities than the respective net probability flux density. Both the net and the unidirectional flux densities in the transition C 0↔O 0 reach finally the respective initial values before applying fcAMP (Fig. 6, top left). The results in Fig. 6 also suggest that the binding of two ligands reduces the conformational flexibility of the closed-open isomerization maximally and, more generally, how different the effects of the four binding steps on the closed-open isomerization are.

Discussion

The C4L-O4L model

We present a detailed analysis of probability fluxes within a Markovian model describing ligand-induced activation of HCN2 channels. The analysis is based on the combined recording of ligand binding and channel activation [22]. The main results are that significant activation of the channel proceeds with one, three, and four ligands bound whereas significant deactivation proceeds relevantly with two, one and zero ligands bound. The consequence of this result is a pronounced hysteresis for channel activation and deactivation. In addition to this we show that the channel is in a flexible conformation with one, three, and four ligands bound whereas it is in a taught conformation with zero and two ligands bound. Notably, we herein considered at each degree of liganding probability fluxes for the whole channel, i.e. we did not use specific stoichiometric factors. This allowed us to avoid any further assumptions concerning equivalence or non-equivalence of the available binding sites.

The high degree of determinateness of our approach was certainly not only caused by analyzing data of ligand binding and gating simultaneously, but also by the specific nature of the effect induced by the concentration jumps: Applying the ligand abruptly transformed the simple C 0↔O 0 model into the complex C4L-O4L model. Conversely, removing the ligand emptied the C4L-O4L model and led to the initial C 0↔O 0 model. Only few parameters were not determined. The limited time resolution of our method did not allow us to distinguish some of the rates with three and four ligands bound. It should be emphasized, however, that none of the present conclusions depends on specific assumptions in the fit.

Transition path analysis in ion channel gating

For analyzing the action of proteins great progress has been achieved over the past years by molecular dynamics (MD) simulations and building energy landscapes. In these landscapes often a small set of energy basins can be identified that are separated by energy barriers (for review see [26]). However, present MD simulations are typically limited by the available computer technologies to one microsecond [5] or, very recently, even to a millisecond [27]–[29], thereby already touching the time range of the observables in ion channel gating. In contrast, functional measurements are real and they provide a different kind of information about the action of ion channels compared to MD simulations: The pore opening can be monitored by the ion current, distinct conformational changes by changes in the fluorescence intensity of introduced labels [30]–[32] and, in case of voltage-gated channels, the movement of the gating apparatus by gating currents [33]–[35]. Markovian state models are often used to interpret functional data.

Over the past decade there has been great progress in applying the transition path theory in combination with Markovian state models to learn more about the action of proteins. However, this was done only in theory by MD simulations and in the mentioned short time range [23], [25], [36]. Also graphical visualization of the paths has been performed [37]. So far, the transition path theory has not been applied to the gating of ion channels as performed herein.

Our approach quantifies at which degree of liganding the closed-open isomerization proceeds and demonstrates a complex type of hysteresis in channel gating (Fig. 4). Though this information is, of course, determined by the rate constants, the complexity of the C4L-O4L model precludes an immediate identification. Such a hysteresis might have physiological consequences for the regulation of the channels. For example, molecules regulating the activity of these channels might have different affinities at different conformations. Also, knowledge of the transition pathways might become very helpful for future strategies to develop drugs stimulating or blocking these channels. For example, if one intends to reduce the cAMP effect on the open probability of HCN2 channels in a living cell, it might be a good idea to selectively affect the C 1↔O 1 transition because it is this transition which mediates the main effect in the lower concentration range.

Net and unidirectional probability flux densities

In addition to the sum of the probability fluxes determined by the transition path analysis, our approach provided us dynamic information: the net and unidirectional probability flux densities. The net probability flux density shows the flux for the different transitions as function of time, i.e. how often per time interval a net transition appears (Fig. 2). There are two important aspects of this information: First, these time courses show the maximum intensities of the transitions. Second, these time courses show how differently rapid these transitions are. This might be also of interest for the development of highly specific drugs modulating the channel activities. Moreover, this knowledge might be of relevance for learning how the subunits interact, how the identified complex cooperativity [22] is generated by subunits that are equal in their amino acid sequence.

The unidirectional probability flux density (Fig. 6) provides another type of information: It specifies all transitions from one state to a neighbor state, including also the repetitive transitions. The unidirectional probability flux density can be much bigger than the net probability flux density. Hence the unidirectional probability flux density provides information how taut or relaxed an HCN2 channel is at a given degree of liganding. The most challenging result is that the tautness is not a monotonous function of the degree of liganding, it is high when none or two ligands are bound and low when one, three and four ligands are bound (Fig. 6, left). In analogy to the T and R conformation in hemoglobin [8], this result suggests that the binding of the first, third and fourth ligand leads to a break of hydrogen bonds whereas the binding of the second ligand promotes the formation of hydrogen bonds.

One might be inclined to relate our results to entropy changes associated with the transition state upon channel activation under the assumption of reducing the complex activation to a single transition and employing Eyring's rate theory [6]. This has been performed for multiple channels, e.g. voltage-gated K+ channels [38]–[40] and voltage-gated Na+ channels [41]. For related CNGA2 channels we reported previously that the entropy of the open channel plus its environment is smaller than that of the closed channel plus its environment [42], similar to the results in the voltage-gated K+ and Na+ channels. For HCN channels respective data are missing.

The fact that the largest unidirectional probability flux density observed herein was observed for the fully liganded channel seems to suggest that the entropy of the channel is higher in the fully liganded state than in the incompletely liganded state. However, this measure of channel flexibility is fundamentally different from a thermodynamic entropy because the environment is not considered. Therefore, the elevated channel flexibility at full liganding indicates a selective property of closed-open isomerization, irrespective of its environment. For further analysis of the gating process in HCN2 channels it would be a good idea to repeat our experiments at different temperatures and study the temperature dependency of the individual rate constants.

Conclusion

We analyzed the gating of homotetrameric HCN2 channels by probability fluxes in a complex Markovian model, the C4L-O4L model. This analysis was not based on simulations but on a fit to functional experimental data. Studying time courses of the net probability flux density, the unidirectional probability flux density, the total net probability flux, and the transition paths within the C4L-O4L model provided us an unusual view on the gating of the channels. Most remarkably, there is considerable ligand-induced channel opening already after the first ligand has bound whereas there is practically no further opening after the second ligand has bound. Moreover, our analysis shows pronounced hysteresis associated with channel opening and closure. Our results should help to better understand the physiological function of HCN channels and, possibly, to develop strategies for a pharmacological modulation of special functional states.

Materials and Methods

All calculations are based on a kinetic model for the ligand-induced activation of HCN2 channels containing four ligand binding steps in both the closed and open channel and five closed-open isomerizations (C4L-O4L model termed herein; Fig. 1) [22]. The channels were activated by a voltage pulse from −30 mV to −130 mV and fcAMP was applied only after the voltage-induced activation was maximal. Fixing the voltage to −130 mV provided us the advantage to study the ligand-induced gating independent of any intervening effects of a changed voltage. The rate constants, determined previously by a global fit strategy of multiple time courses of ligand binding and unbinding as well as activation and deactivation gating, are listed together with their s.e.m. in Table S1. The time courses in the present study were computed by using the mean rate constants. The time-dependent occupancies of the states were computed with the Eigenvalue method using Matlab®. Initially, the channel was assumed to be in the equilibrium C 0↔O 0 because the ligand concentration, L, was zero (Fig. 1). Then L was set to the value of the applied fcAMP concentration (activation) and back to zero (deactivation). For the numerical computation of time integrals, Simpson coefficients were used.

Rate constants for the C4L-O4L model. The rate constants were computed by the global fit as described [22].

(DOC)

Click here for additional data file.