Two-dimensional components and hidden dependencies provide insight into ion channel gating mechanisms

Correlations between the durations of adjacent open and shut intervals recorded from ion channels contain information about the underlying gating mechanism. This study presents an additional approach to extracting the correlation information. Detailed correlation information is obtained directly from single-channel data and quantified in a manner that can provide insight into the connections among the states underlying the gating. The information is obtained independently of any specific kinetic scheme, except for the general assumption of Markov gating. The durations of adjacent open and shut intervals are binned into two-dimensional (2-D) dwell-time distributions. The 2-D (joint) distributions are fitted with sums of 2-D exponential components to determine the number of 2-D components, their volumes, and their open and closed time constants. The dependency of each 2-D component is calculated by comparing its observed volume to the volume that would be expected if open and shut intervals paired independently. The estimated component dependencies are then used to suggest gating mechanisms and to provide a powerful means of examining whether proposed gating mechanisms have the correct connections among states. The sensitivity of the 2-D method can identify hidden components and dependencies that can go undetected by previous correlation methods.     

Testing for microscopic reversibility in the gating of maxi K+ channels using two-dimensional dwell-time distributions.

An assumption usually made when developing kinetic models for the gating of ion channels is that the transitions among the various states involved in the gating obey microscopic reversibility. If this assumption is incorrect, then the models and estimated rate constants made with the assumption would be in error. This paper examines whether the gating of a large conductance Ca-activated K+ channel in skeletal muscle is consistent with microscopic reversibility. If microscopic reversibility is obeyed, then the number of forward and backward transitions per unit time for each individual reaction step will, on average, be identical and, consequently, the gating must show time reversibility. To look for time reversibility, two-dimensional dwell-time distributions of the durations of open and closed intervals were obtained from single-channel current records analyzed in the forward and in the backward directions. Two-dimensional dwell-time distributions of pairs of open intervals and of pairs of closed intervals were also analyzed to extend the resolution of the method to special circumstances in which intervals from different closed (or open) states might have similar durations. No significant differences were observed between the forward and backward analysis of the two-dimensional dwell-time distributions, suggesting time reversibility. Thus, we find no evidence to indicate that the gating of the maxi K+ channel violates microscopic reversibility.     

Inversion of Markov processes to determine rate constants from single-channel data.

The determination of rate constants from single-channel data can be very difficult, in part because the single-channel lifetime distributions commonly analyzed by experimenters often have a complicated mathematical relation to the channel gating mechanism. The standard treatment of channel gating as a Markov process leads to the prediction that lifetime distributions are exponential functions. As the number of states of a channel gating scheme increases, the number of exponential terms in the lifetime distribution increases, and the weights and decay constants of the lifetime distributions become progressively more complicated functions of the underlying rate constants. In the present study a mathematical strategy for inverting these functions is introduced in order to determine rate constants from single-channel lifetime distributions. This inversion is easy for channel gating schemes with two or fewer states of a given conductance, so the present study focuses on schemes with more states. The procedure is to derive explicit equations relating the parameters of the lifetime distribution to the rate constants of the scheme. Such equations can be derived using the equality between symmetric functions of eigenvalues of a matrix and sums over principle minors, as well as expressions for the moments, derivatives, and weights of a lifetime distribution. The rate constants are then obtained as roots to this system of equations. For a gating scheme with three sequential closed states and a single gateway state, exact analytical expressions were found for each rate constant in terms of the parameters of the three-exponential closed-time distribution. For several other gating schemes, systems of equations were found that could be solved numerically to obtain the rate constants. Lifetime distributions were shown to specify a unique set of real rate constants in sequential gating schemes with up to five closed or five open states. For kinetic schemes with multiple gating pathways, the analysis of simulated data revealed multiple solutions. These multiple solutions could be distinguished by examining two-dimensional probability density functions. The utility of the methods introduced here are demonstrated by analyzing published data on nicotinic acetylcholine receptors, GABA(A) receptors, and NMDA receptors.     

Gating kinetics of potassium channel in rat dorsal root ganglion neurons analyzed with fractal model

The kinetics of ion channels have been widely modeled as a Markov process. In these models it is assumed that the channel protein has a small number of discrete conformational states and kinetic rate constants connecting these states are constant. To study the gating kinetics of voltage-dependent K(+) channel in rat dorsal root ganglion neurons, K(+) channel current were recorded using cell-attached patch-clamp technique. The K(+) channel characteristic of kinetics were found to be statistically self-similar at different time scales as predicted by the fractal model. The fractal dimension D for the closed times and for the open times depend on the pipette potential. For the open and closed times of kinetic setpoint, it was found dependent on the applied pipette potential, which indicated that the ion channel gating kinetics had nonlinear kinetic properties. Thus, the open and closed durations, which had the voltage dependence of the gating of this ion channel, were well described by the fractal model.     

Single-channel currents from diethylpyrocarbonate-modified NMDA receptors in cultured rat brain cortical neurons

The role of histidine residues in the function of N-methyl-D-aspartate (NMDA)-activated channels was tested with the histidine-modifying reagent diethylpyrocarbonate (DEP) applied to cells and membrane patches from rat brain cortical neurons in culture. Channels in excised outside-out patches that were treated with 3 mM DEP for 15-30 s (pH 6.5) showed an average 3.4-fold potentiation in steady state open probability when exposed to NMDA and glycine. Analysis of the underlying alterations in channel gating revealed no changes in the numbers of kinetic states: distributions of open intervals were fitted with three exponential components, and four components described the shut intervals, in both control and DEP-modified channels. However, the distribution of shut intervals was obviously different after DEP treatment, consistent with the single-channel current record. After modification, the proportion of long shut states was decreased while the time constants were largely unaffected. Burst kinetics reflected these effects with an increase in the average number of openings/burst from 1.5 (control) to 2.2 (DEP), and a decrease in the average interburst interval from 54.1 to 38.2 ms. These effects were most likely due to histidine modification because other reagents (n- acetylimidazole and 2,4,6-trinitrobenzene 1-sulfonic acid) that are specific for residues other than histidine failed to reproduce the effects of DEP, whereas hydroxylamine could restore channel open probability to control levels. In contrast to these effects on channel gating, DEP had no effect on average single-channel conductance or reversal potential under bi-ionic (Na+:Cs+) conditions. Inhibition by zinc was also unaffected by DEP. We propose a channel gating model in which transitions between single- and multi-opening burst modes give rise to the channel activity observed under steady state conditions. When adjusted to account for the effects of DEP, this model suggests that one or more extracellular histidine residues involved in channel gating are associated with a single kinetic state.     

Gating of single non-Shaker A-type potassium channels in larval Drosophila neurons

The voltage-dependent gating of transient A2-type potassium channels from primary cultures of larval Drosophila central nervous system neurons was studied using whole-cell and single-channel voltage clamp. A2 channels are genetically distinct from the Shaker A1 channels observed in Drosophila muscle, and differ in single-channel conductance, voltage dependence, and gating kinetics. Single A2 channels were recorded and analyzed at -30, -10, +10, and +30 mV. The channels opened in bursts in response to depolarizing steps, with three to four openings per burst and two to three bursts per 480-ms pulse (2.8-ms burst criterion). Mean open durations were in a range of 2-4 ms and mean burst durations in a range of 9-17 ms. With the exception of the first latency distributions, none of the means of the distributions measured showed a consistent trend with voltage. Macroscopic inactivation of both whole-cell A currents and ensemble average currents of single A2 channels was well fitted by a sum of two exponentials. The fast time constants in different cells were in a range of 9-25 ms, and the slow time constants in a range of 60-140 ms. A six-state kinetic model (three closed, one open, two inactivated states) was tested at four command voltages by fitting frequency histograms of open durations, burst durations, burst closed durations, number of openings per burst, and number of bursts per trace. The model provided good fits to these data, as well as to the ensemble averages. With the exception of the rates leading to initial opening, the transitions in the model were largely independent of voltage.     

Fractal model of ion-channel kinetics

Markov models with discrete states, such as closed in equilibrium with closed in equilibrium with open have been widely used to model the kinetics of ion channels in the cell membrane. In these models the transition probabilities per unit time (the kinetic rate constants) are independent of the time scale on which they are measured. However, in many physical systems, a property, L, depends on the scale, epsilon, at which it is measured such that L(epsilon) alpha epsilon 1-D where D is the fractal dimension. Such systems are said to be 'fractal'. Based on the assumption that the kinetic rates are given by k(t) alpha t1-D we derive a fractal model of ion-channel kinetics. This fractal model has fewer adjustable parameters, is more consistent with the dynamics of protein conformations, and fits the single-channel recordings from the corneal endothelium better than the discrete-state Markov model.     

Fractal and Markov behavior in ion channel kinetics

Kinetic analysis of ion channel recordings attempts to distinguish the number and lifetimes of channel molecular states. Most kinetic analysis assumes that the lifetime of each state is independent of previous channel history, so that open and closed durations are Markov processes whose probability densities are sums of exponential decays. An alternative approach assumes that channel molecules have many configurtions with widely varying lifetimes. Rates of opening and closing then vary with the time scale of observation, leading to fractal kinetics. We have examined kinetic behavior in two types of channels from human and avian fibroblasts, using a maximum likehood method to test the dependence of rates on observational time scale. For both channels, openings showed mixed fractal and Markov behavior, while closings gave mainly fractal kinetics.     

Theory of the kinetic analysis of patch-clamp data.

This paper describes a theory of the kinetic analysis of patch-clamp data. We assume that channel gating is a Markov process that can be described by a model consisting of n kinetic states and n(n - 1) rate constants at each voltage, and that patch-clamp data describe the occupancy of x different conductance levels over time. In general, all the kinetic information in a set of patch-clamp data is found in either two-dimensional dwell time histograms describing the frequency of observation of sequential dwell times of durations tau 1 and tau 2 (Fredkin, D. R., M. Montal, and J. A. Rice, 1985, Proceedings of the Berkeley Conference in Honor of Jerzy Neyman and Jack Kiefer, vol. 1, 269-289) or in three-point joint probability functions describing the probability that a channel is in a given conductance at time t, and at time t + tau 1, and at time t + tau 1 + tau 2. For the special case of channels with a single open state plus multiple closed states, one-dimensional analyses provide all of the kinetic information. Stationary patch-clamp data have information that can be used to determine H rate constants, where H = n(n - 1) - G and G is the number of intraconductance rate constants. Thus, to calculate H rate constants, G rate constants must be fixed. In general there are multiple sets of G rate constants that can be fixed to allow the calculation of H rate constants although not every set of G rate constants will work. Arbitrary assignment of the G intraconductance rate constants equal to zero always provides a solution and the calculation of H rate constants. Nonstationary patch-clamp data have information for the determination of H rate constants at a reference voltage plus n(n - 1) rate constants at all test voltages. Thus, nonstationary data have extra information about the voltage dependencies of rate constants that can be used to rule out kinetic models that cannot be disqualified on the basis of stationary data.     

Linking exponential components to kinetic states in Markov models for single-channel gating

Discrete state Markov models have proven useful for describing the gating of single ion channels. Such models predict that the dwell-time distributions of open and closed interval durations are described by mixtures of exponential components, with the number of exponential components equal to the number of states in the kinetic gating mechanism. Although the exponential components are readily calculated (Colquhoun and Hawkes, 1982, Phil. Trans. R. Soc. Lond. B. 300:1-59), there is little practical understanding of the relationship between components and states, as every rate constant in the gating mechanism contributes to each exponential component. We now resolve this problem for simple models. As a tutorial we first illustrate how the dwell-time distribution of all closed intervals arises from the sum of constituent distributions, each arising from a specific gating sequence. The contribution of constituent distributions to the exponential components is then determined, giving the relationship between components and states. Finally, the relationship between components and states is quantified by defining and calculating the linkage of components to states. The relationship between components and states is found to be both intuitive and paradoxical, depending on the ratios of the state lifetimes. Nevertheless, both the intuitive and paradoxical observations can be described within a consistent framework. The approach used here allows the exponential components to be interpreted in terms of underlying states for all possible values of the rate constants, something not previously possible.     

Fractal analysis of a voltage-dependent potassium channel from cultured mouse hippocampal neurons.

The kinetics of ion channels have been widely modeled as a Markov process. In these models it is assumed that the channel protein has a small number of discrete conformational states and the kinetic rate constants connecting these states are constant. In the alternative fractal model the spontaneous fluctuations of the channel protein at many different time scales are represented by a kinetic rate constant k = At1-D, where A is the kinetic setpoint and D the fractal dimension. Single-channel currents were recorded at 146 mM external K+ from an inwardly rectifying, 120 pS, K+ selective, voltage-sensitive channel in cultured mouse hippocampal neurons. The kinetics of these channels were found to be statistically self-similar at different time scales as predicted by the fractal model. The fractal dimensions were approximately 2 for the closed times and approximately 1 for the open times and did not depend on voltage. For both the open and closed times the logarithm of the kinetic setpoint was found to be proportional to the applied voltage, which indicates that the gating of this channel involves the net inward movement of approximately one negative charge when this channel opens. Thus, the open and closed times and the voltage dependence of the gating of this channel are well described by the fractal model.     

Two-dimensional kinetic analysis suggests nonsequential gating of mechanosensitive channels in Xenopus oocytes

Xenopus oocytes express mechanosensitive (MS(XO)) channels that can be studied in excised patches of membrane with the patch-clamp technique. This study examines the steady-state kinetic gating properties of MS(XO) channels using detailed single-channel analysis. The open and closed one-dimensional dwell-time distributions were described by the sums of 2-3 open and 5-7 closed exponential components, respectively, indicating that the channels enter at least 2-3 open and 5-7 closed kinetic states during gating. Dependency plots revealed that the durations of adjacent open and closed intervals were correlated, indicating two or more gateway states in the gating mechanism for MS channels. Maximum likelihood fitting of two-dimensional dwell-time distributions to both generic and specific models was used to examine gating mechanism and rank models. A kinetic scheme with five closed and five open states, in which each closed state could make a direct transition to an open state (two-tiered model) could account for the major features of the single-channel data. Two-tiered models that allowed direct transitions to subconductance open states in addition to the fully open state were also consistent with multiple gateway states. Thus, the gating mechanism of MS(XO) channels differs from the sequential (linear) gating mechanisms considered for MS channels in bacteria, chick skeletal muscle, and Necturus proximal tubule.     

The Estimation of Rapid Rate Constants from Current-Amplitude Frequency Distributions of Single-Channel Recordings

A method is described for estimating rapid rate constants from the distributions of current amplitude observed in single-channel electrical recordings. It has the advantages over previous, similar approaches that it can accommodate both multistate kinetic models and adjustable filtering of the data using an 8-pole Bessel filter. The method is conceptually straightforward: the observed distributions of current amplitude are compared with theoretical distributions derived by combining several simplifying assumptions about the underlying stochastic process with a model of the filter and electrical noise. Parameters are estimated by approximate maximum likelihood. The method was used successfully to estimate rate constants for both a simple two-state kinetic model (the transitions between open and closed states during the rapid gating of an outward-rectifying K⁺-selective channel in the plasma membrane of Acetabularia) and a complex multistate kinetic model (the blockade of the maxi cation channel in the plasma membrane of rye roots by verapamil). For the two-state model, parameters were estimated well, provided that they were not too fast or too slow in relation to the sampling rate. In the three-state model the precision of estimates depended in a complex way on the values of all rate parameters in the model. Received: 4 October 1996/Revised: 2 September 1997     

Dependency plots suggest the kinetic structure of ion channels.

Ion channels are integral membrane proteins that regulate ionic flux through cell membranes by opening and closing (gating) their pores. The gating can be monitored by observing step changes in the current flowing through single channels, and analysis of the observed open and closed interval durations has provided a window to develop kinetic models for the gating process. One difficulty in developing such models has been to determine the connections (transition pathways) among the various kinetic states involved in the gating. To help overcome this difficulty we present a transform (dependency plot) of the single-channel data that can give immediate insight into the connections. A dependency plot is derived by calculating a contingency table from a two-dimensional (joint density) dwell-time distribution of adjacent open and closed intervals by assuming that the two classified criteria are the open and closed durations of each pair of adjacent intervals. A three-dimensional surface plot of the fractional difference between the numbers of observed interval pairs and the numbers expected if the durations of adjacent intervals are independent then gives the dependency plot. An excess of interval pairs in the dependency plot suggests that the open and closed states (or compound states) that give rise to the interval pairs in excess are directly connected. A deficit of interval pairs suggests that the open and closed states (or compound states) that give rise to the interval pairs in deficit are either not directly connected or that there are additional open-closed transition pathways arising from the directly connected states.     

Gating Models of the Anomalous Mole-Fraction Effect of Single-Channel Current in <Emphasis Type="Italic">Chara</Emphasis>

The dependence of single-channel current on the Tl+/K+ mole fraction exhibiting a minimum at [Tl+]/[K+] of about 1:15 is proportional to open probability in bursts. Five models are suggested to explain modulation of gating by the Tl+/K+ ratio. Three models start from a channel with 4 identical subunits, each with an allosteric binding site for K+ or Tl+. In the first model, ion binding is directly observable as a transition from one Markov state to another. This model can explain the dependence of the apparent single-channel current on Tl+ concentrations. However, the predicted linear dependence on ion concentrations of the apparent rate constants was not observed in measurements in 25 or 250 mM KNO3 and 250 mM Tl NO3. The second model can overcome this problem by introducing saturation kinetics for ion binding. In the third model, gating is caused by inherent vibrations of the protein, and the rate constants of the related transitions depend on the occupation of the allosteric sites. The fourth model is based on the foot-in-the-door approach with the essential feature that two K+ ions in the selectivity filter are necessary to keep the pore radius suitable for K+ ions. The fifth model is also a foot-in-the-door model, but non-Markovian because, similar to model 3, it is assumed that the conformation of the protein (and thus the rate constants of the Markov model of the time series) depends on the force exerted by the temporal average over the states of a Markov model of ion occupation. These ions may reside in the pore itself or outside.     

Gating kinetics of batrachotoxin-modified Na+ channels in the squid giant axon. Voltage and temperature effects.

The gating kinetics of batrachotoxin-modified Na+ channels were studied in outside-out patches of axolemma from the squid giant axon by means of the cut-open axon technique. Single channel kinetics were characterized at different membrane voltages and temperatures. The probability of channel opening (Po) as a function of voltage was well described by a Boltzmann distribution with an equivalent number of gating particles of 3.58. The voltage at which the channel was open 50% of the time was a function of [Na+] and temperature. A decrease in the internal [Na+] induced a shift to the right of the Po vs. V curve, suggesting the presence of an integral negative fixed charge near the activation gate. An increase in temperature decreased Po, indicating a stabilization of the closed configuration of the channel and also a decrease in entropy upon channel opening. Probability density analysis of dwell times in the closed and open states of the channel at 0 degrees C revealed the presence of three closed and three open states. The slowest open kinetic component constituted only a small fraction of the total number of transitions and became negligible at voltages greater than -65 mV. Adjacent interval analysis showed that there is no correlation in the duration of successive open and closed events. Consistent with this analysis, maximum likelihood estimation of the rate constants for nine different single-channel models produced a preferred model (model 1) having a linear sequence of closed states and two open states emerging from the last closed state. The effect of temperature on the rate constants of model 1 was studied. An increase in temperature increased all rate constants; the shift in Po would be the result of an increase in the closing rates predominant over the change in the opening rates. The temperature study also provided the basis for building an energy diagram for the transitions between channel states.