A data-driven model of a modal gated ion channel: The inositol 1,4,5-trisphosphate receptor in insect Sf9 cells

The inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) channel is crucial for the generation and modulation of intracellular Ca(2+) signals in animal cells. To gain insight into the complicated ligand regulation of this ubiquitous channel, we constructed a simple quantitative continuous-time Markov-chain model from the data. Our model accounts for most experimentally observed gating behaviors of single native IP(3)R channels from insect Sf9 cells. Ligand (Ca(2+) and IP(3)) dependencies of channel activity established six main ligand-bound channel complexes, where a complex consists of one or more states with the same ligand stoichiometry and open or closed conformation. Channel gating in three distinct modes added one complex and indicated that three complexes gate in multiple modes. This also restricted the connectivity between channel complexes. Finally, latencies of channel responses to abrupt ligand concentration changes defined a model with specific network topology between 9 closed and 3 open states. The model with 28 parameters can closely reproduce the equilibrium gating statistics for all three gating modes over a broad range of ligand concentrations. It also captures the major features of channel response latency distributions. The model can generate falsifiable predictions of IP(3)R channel gating behaviors and provide insights to both guide future experiment development and improve IP(3)R channel gating analysis. Maximum likelihood estimates of the model parameters and of the parameters in the De Young-Keizer model yield strong statistical evidence in favor of our model. Our method is simple and easily applicable to the dynamics of other ion channels and molecules.     

Statistical properties of ion channel records. Part II: estimation from the macroscopic current

Macroscopic ion channel current is the summation of the stochastic records of individual channel currents and therefore relates to their statistical properties. As a consequence of this relationship, it may be possible to derive certain statistical properties of single channel records or even generate some estimates of the records themselves from the macroscopic current when the direct measurement of single channel currents is not applicable. We present a procedure for generating the single channel records of an ion channel from its macroscopic current when the stochastic process of channel gating has the following two properties: (I) the open duration is independent of the time of opening event and has a single exponential probability density function (pdf), (II) all the channels have the same probability to open at time t. The application of this procedure is considered for cases where direct measurement of single channel records is difficult or impossible. First, the probability density function (pdf) of opening events, a statistical property of single channel records, is derived from the normalized macroscopic current and mean channel open duration. Second, it is shown that under the conditions (I) and (II), a non-stationary Markov model can represent the stochastic process of channel gating. Third, the non-stationary Markov model is calibrated using the results of the first step. The non-stationary formulation increases the model ability to generate a variety of different single channel records compared to common stationary Markov models. The model is then used to generate single channel records and to obtain other statistical properties of the records. Experimental single channel records of inactivating BK potassium channels are used to evaluate how accurately this procedure reconstructs measured single channel sweeps.     

A novel hypothesis for the binding mode of HERG channel blockers

We present a new docking model for HERG channel blockade. Our new model suggests three key interactions such that (1) a protonated nitrogen of the channel blocker forms a hydrogen bond with the carbonyl oxygen of HERG residue T623; (2) an aromatic moiety of the channel blocker makes a pi-pi interaction with the aromatic ring of HERG residue Y652; and (3) a hydrophobic group of the channel blocker forms a hydrophobic interaction with the benzene ring of HERG residue F656. The previous model assumes two interactions such that (1) a protonated nitrogen of the channel blocker forms a cation-pi interaction with the aromatic ring of HERG residue Y652; and (2) a hydrophobic group of the channel blocker forms a hydrophobic interaction with the benzene ring of HERG residue F656. To test these models, we classified 69 known HERG channel blockers into eight binding types based on their plausible binding modes, and further categorized them into two groups based on the number of interactions our model would predict with the HERG channel (two or three). We then compared the pIC(50) value distributions between these two groups. If the old hypothesis is correct, the distributions should not differ between the two groups (i.e., both groups show only two binding interactions). If our novel hypothesis is correct, the distributions should differ between Groups 1 and 2. Consistent with our hypothesis, the two groups differed with regard to pIC(50), and the group having more predicted interactions with the HERG channel had a higher mean pIC(50) value. Although additional work will be required to further validate our hypothesis, this improved understanding of the HERG channel blocker binding mode may help promote the development of in silico predictions methods for identifying potential HERG channel blockers.     

Biophysics of mechanoreception

Several types of cells' skeletal, muscle, nerve, epithelia, and heart have been shown to contain ion channels which are sensitive to membrane tension. In chick skeletal muscle, the transduction persists in excised patches and involves no chemical messengers. Quantitative analysis of single channel records reveals that the sensitivity to stretch can be described by a linear four state model with three closed (C) and one open (O) state: (Formula: see text). Only the rate constant k12 is sensitive to tension (and membrane potential) following the law: k12 = kO12 exp/(theta T2 + alpha V) where theta is a constant describing the sensitivity to tension, T, and alpha is a constant describing the sensitivity to voltage, V, and kO12 is a constant. The form of the tension sensitivity can be accounted for by a model in which strain energy is used to gate the channel. Analysis of strain sensitivity, theta, indicates that the channel must concentrate energy from a large (ca. 500-nm diameter) area of membrane which suggests that the channel is in series with a component of the cytoskeleton. Treatment with cytochalasins suggests that actin is mechanically in parallel with the channel. When a channel with the above properties is incorporated into a simple model of mechanical transduction in hair cells, the resulting model is capable of explaining the kinetic features and the sensitivity found in the cochlear-vestibular system. The proposed gating mechanism of mechanical transduction appears to be general and can account for existing data on a variety of systems.     

Mechanisms of atrial-selective block of Na⁺ channels by ranolazine: II. Insights from a mathematical model

Block of Na(+) channel conductance by ranolazine displays marked atrial selectivity that is an order of magnitude higher that of other class I antiarrhythmic drugs. Here, we present a Markovian model of the Na(+) channel gating, which includes activation-inactivation coupling, aimed at elucidating the mechanisms underlying this potent atrial selectivity of ranolazine. The model incorporates experimentally observed differences between atrial and ventricular Na(+) channel gating, including a more negative position of the steady-state inactivation curve in atrial versus ventricular cells. The model assumes that ranolazine requires a hydrophilic access pathway to the channel binding site, which is modulated by both activation and inactivation gates of the channel. Kinetic rate constants were obtained using guarded receptor analysis of the use-dependent block of the fast Na(+) current (I(Na)). The model successfully reproduces all experimentally observed phenomena, including the shift of channel availability, the sensitivity of block to holding or diastolic potential, and the preferential block of slow versus fast I(Na.) Using atrial and ventricular action potential-shaped voltage pulses, the model confirms significantly greater use-dependent block of peak I(Na) in atrial versus ventricular cells. The model highlights the importance of action potential prolongation and of a steeper voltage dependence of the time constant of unbinding of ranolazine from the atrial Na(+) channel in the development of use-dependent I(Na) block. Our model predictions indicate that differences in channel gating properties as well as action potential morphology between atrial and ventricular cells contribute equally to the atrial selectivity of ranolazine. The model indicates that the steep voltage dependence of ranolazine interaction with the Na(+) channel at negative potentials underlies the mechanism of the predominant block of I(Na) in atrial cells by ranolazine.     

A three-barrier model for the hemocyanin channel

Keyhole limpet hemocyanin forms ion-conducting channels in planar lipid bilayer membranes. Ionic current through the open hemocyanin channel presents the following characteristics: (a) it is carried mainly by cations; (b) it is a nonlinear function of membrane potential; (c) channel conductance is a saturating function of ion activity; (d) it shows ionic competition. A model for the open hemocyanin channel is developed from absolute reaction rate theory. The model calls for three energy barriers in the channel. Two energy barriers represent the entrance and exit of the ion into and out of the channel. The third barrier separates two energy minima that represent two binding sites. Furthermore, only one ion is allowed inside the channel at a given time. This model is able to recreate all the hemocyanin characteristics found experimentally in negatively charged and neutral membranes.     

Cardiac sodium channel Markov model with temperature dependence and recovery from inactivation

A Markov model of the cardiac sodium channel is presented. The model is similar to the CA1 hippocampal neuron sodium channel model developed by Kuo and Bean (1994. Neuron. 12:819-829) with the following modifications: 1) an additional open state is added; 2) open-inactivated transitions are made voltage-dependent; and 3) channel rate constants are exponential functions of enthalpy, entropy, and voltage and have explicit temperature dependence. Model parameters are determined using a simulated annealing algorithm to minimize the error between model responses and various experimental data sets. The model reproduces a wide range of experimental data including ionic currents, gating currents, tail currents, steady-state inactivation, recovery from inactivation, and open time distributions over a temperature range of 10 degrees C to 25 degrees C. The model also predicts measures of single channel activity such as first latency, probability of a null sweep, and probability of reopening.     

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Estimation of Na+ dwell time in the gramicidin A channel. Na+ ions as blockers of H+ currents

The permeation of Na+ through gramicidin A channels shows a simple saturation with increasing Na+ concentration that can be described by two different models. The first model assumes that one Na+ binds to the channel with high affinity (approximately 30 M-1) and that conduction occurs by a 'knock-on' mechanism requiring double occupancy of the channel; the other model assumes that Na+ binding is of low affinity (less than 1 M-1), and that double occupancy of the channel is rare. NMR measurements have shown tight Na+ binding, favoring the first model, but measurements of flux ratios and water transport support the second model. We present here a relatively model-independent measurement of the dwell time of Na+ inside the channel, in which we characterize the fluctuations in H+ current through the channel induced by 'block' from the more slowly permeating Na+ ions. The mean Na+ dwell time inside the channel is estimated to be approximately 10 ns at a membrane potential of 200 mV. This result is inconsistent with tight Na+ binding, thus favoring the second model.     

A new double-chamber model of ion channels. Beyond the Hodgkin and Huxley model

This paper proposes a new double-chamber model (DCM) of ion channels. The model ion channel consists of a series of three pores alternating with two chambers. The chambers are net negatively charged. The chamber's electric charge originates from dissociated amino acid side chains and is pH dependent. The chamber's net negative charge is compensated by cations present inside the chamber and in a diffuse electric layer outside the chamber. The pore's permeability is constant independent of time. One pore of the sodium channel and one of the potassium channel is a voltage-sensing pore. Due to the channel's structure, ions flow through the pores and chambers in a time-dependent manner. The model reproduces experimental voltage clamp and action potential data. The current flowing through a single sodium channel is less then one femtoampere. The DCM is considerably simpler then the Hodgkin and Huxley model (HHM) used to describe the electrophysiological properties of an axon. Unlike the HHM, the DCM can explain refractoriness, anode break excitation, accommodation and the effect of pH and temperature on the channels without additional parameters. In the DCM, the axon membrane shows repetitive activity depending on the channel density, sodium to potassium channel ratio and external potassium concentration. In the DCM, the action potential starts from 'hot spot areas' of higher channel densities and a higher sodium to potassium channel ratio, and then propagates through the whole axon.     

A Model of sodium channel-inactivation based upon the modulated blocker

An attempt is made to model sodium channel inactivation based upon real physical processes. The principle involved, which is supported by calculation and by direct appeal to experimental results, is that the gating dipole reversal or gating charge transfer that occurs when the channel is activated, markedly modulates the electrical properties of charged groups at the channel ends. Four examples of possible mechanisms that lead to channel inactivation are described. The simple four-state model that results is able to predict: (a) the steep voltage dependence of the equilibrium inactivation characteristic without the presence of any appreciable displacement current associated with inactivation; (b) the negative shift in membrane voltage of the equilibrium inactivation characteristic relative to the activation characteristic; (c) the bell-shaped dependence of inactivation time constant on membrane voltage; (d) the similarity of the membrane voltage dependence of the time constant of recovery from inactivation, to that of inactivation itself. A brief discussion of a model for sodium channel activation based upon the same physical principle is included.     

Stochastic description of the three-state model of the Na channel

The three-state model of the Na channel is formulated stochastically and various observable quantities calculated from this model. This model assumes that the Na channel can occupy one of three states--resting, open and inactivated--and that each channel is independent. This is a simplified representation of the model proposed by Aldrich et al. (1983) mainly with respect to the neuroblastoma. When the system contains many channels, the probability distribution of the numbers of the channels that occupy each of these states is a time-dependent multinomial distribution and the distribution of the first passage time from resting state to open state becomes an exponential decay function with higher components. The average and variance of the channel current calculated by the distribution show the time course with a single peak and its ratio converges to the current of a single channel. Stochastic treatment of the transient process developed here gives a method of estimating all of the transition rates and the number of channels in the system.     

Modulation of gramicidin A open channel lifetime by ion occupancy.

The hypothesis that the gramicidin A channel stability depends on the level of ion occupancy of the channel was used to derive a mathematical model relating channel lifetime to channel occupancy. Eyring barrier permeation models were examined for their ability to fit the zero-voltage conductance, current-voltage, as well as lifetime data. The simplest permeation model required to explain the major features of the experimental data consists of three barriers and four sites (3B4S) with a maximum of two ions occupying the channel. The average lifetime of the channel was calculated from the barrier model by assuming the closing rate constant to be proportional to the probability of the internal channel sites being empty. The link between permeation and lifetime has as its single parameter the experimentally determined averaged lifetime of gramicidin A channels in the limit of infinitely dilute solutions and has therefore no adjustable parameters. This simple assumption that one or more ions inside the channel completely stabilize the dimer conformation is successful in explaining the experimental data considering the fact that this model for stabilization is independent of ion species and configurational occupancy. The model is used to examine, by comparison with experimental data, the asymmetrical voltage dependence of the lifetime in asymmetrical solutions, the effects of blockers, and the effects of elevated osmotic pressure.     

Gramicidin channel selectivity. Molecular mechanics calculations for formamidinium, guanidinium, and acetamidinium.

Empirical energy function calculations were used to evaluate the effects of minimization on the structure of a gramicidin A channel and to analyze the energies of interaction between three cations (guanidinium, acetamidinium, formamidinium) and the channel as a function of position along the channel axis. The energy minimized model of the gramicidin channel, which was based on the results of Venkatachalam and Urry (1983), has a constriction at the channel entrance. If the channel is not allowed to relax in the presence of the ions (rigid model), there is a large potential energy barrier for all three cations. The barrier varies with cation size and is due to high van der Waals and ion deformation energies. If the channel is minimized in the presence of the ions, the potential energy barrier to formamidinium entry is almost eliminated, but a residual barrier remains for guanidinium and acetamidinium. The residual barrier is primarily due, not to the expansion of the helix, but, to the disruption of hydrogen bonds between the terminal ethanoloamine and the next turn of the helix which occurs when the carbonyls of the outer turn of the helix librate inward toward the ion as it enters the channel. The residual potential energy barriers could be a possible explanation for the measured selectivity of gramicidin for formamidinium over guanidinium. The results of this full-atomic model address the applicability of the size-exclusion concept for the selectivity of the gramicidin channel.     

Subunit interactions in the activation of cyclic nucleotide-gated ion channels.

Cyclic nucleotide-gated (CNG) ion channels of retinal photoreceptors and olfactory neurons are multimeric proteins of unknown stoichiometry. To investigate the subunit interactions that occur during CNG channel activation, we have used tandem cDNA constructs of the rod CNG channel to generate heteromultimeric channels composed of wild-type and mutant subunits. We introduced point mutations that affect channel activation: 1) D604M, which alters the relative ability of agonists to promote the allosteric conformational change(s) associated with channel opening, and 2) T560A, which primarily affects the initial binding affinity for cGMP, and to a lesser extent, the allosteric transition. At saturating concentrations of agonist, heteromultimeric channels were intermediate between wild-type and mutant homomultimers in agonist efficacy and apparent affinity for cGMP, cIMP, and cAMP, consistent with a model for the allosteric transition involving a concerted conformational change in all of the channel subunits. Results were also consistent with a model involving independent transitions in two or three, but not one or four, of the channel subunits. The behavior of the heterodimers implies that the channel stoichiometry is some multiple of 2 and is consistent with a tetrameric quaternary structure for the functional channel complex. Steady-state dose-response relations for homomultimeric and heteromultimeric channels were well fit by a Monod, Wyman, and Changeux model with a concerted allosteric opening transition stabilized by binding of agonist.     

Dynamic properties of the colicin E1 ion channel

Abstract The mechanism of channel formation and action of channel-forming colicins is a paradigm for the study of dynamic aspects of membrane-protein interactions. The following experimental results concerning interaction of the colicin E1 channel domain with target membranes, in vitro and in vivo, are discussed: (1) the nature of the translocation-competent state of the channel-forming domain; (2) unfolding of the colicin channel peptide during in vitro binding and anchoring of the channel to liposome membranes at acidic pH; (3) reversal of channel peptide binding to liposomes by an alkaline-directed pH shift; (4) voltage-driven translocation and gating of the ion channel, discussed in the context of a four-helix model for a monomeric channel; (5) rescue of colicin-treated cells by high levels of external K⁺; (6) trypsin rescue of cells depolarized by the colicin ion channel; and (7) interaction of the channel domain with its immunity protein.     

Electromechanical coupling model of gating the large mechanosensitive ion channel (MscL) of Escherichia coli by mechanical force.

We have developed a theoretical electromechanical coupling (EMC) model of gating of the large-conductance mechanosensitive ion channel (MscL). The model presents the first attempt to explain the pressure-dependent transitions between the closed and open channel conformations on a molecular level by assuming 1) a homohexameric structural model of the channel, 2) electrostatic interactions between various domains of the homohexamer, 3) structural flexibility of the N-terminal portion of the monomer, and 4) mechanically and electrostatically induced displacement of the N-terminal domain relative to other structural domains of the protein. In the EMC model, 12 membrane-spanning alpha-helices (six each of the M1 and M2 transmembrane domains of the MscL monomer), are envisaged to line the channel pore with a diameter of 40 A, whereas the N- and C-termini are oriented toward each other inside the pore when the channel is closed. The model proposes that stretching the membrane bilayer by mechanical force causes the monomers to be pulled away from and slightly tilted toward each other. This relative movement of alpha-helices could serve as a trigger to initiate a "swing-like" motion of the N-terminus around the glycine residue G14 that may act as a pivot. The analysis of the attractive and repulsive coulomb forces between all domains of the channel homohexamer suggested that an inclination angle of approximately 3.0 degrees - 4.1 degrees between the oppositely oriented channel monomers should suffice for the N-terminus to turn away from other domains causing the channel to open. According to the EMC model the minimal free energy change, deltaG, that could initiate the opening of the channel was 2 kT. Also, the model predicted that the negative pressure required for channel open probability, Po = 0.5, should be between 50 and 80 mmHg. These values were in a good agreement with the experimentally estimated pressures of 60-70 mmHg obtained with the MscL reconstituted in liposomes. Furthermore, consistent with a notion that the N-terminus may present a mechanosensitive structural element providing a mechanism to open the MscL by mechanical force, the model provides a simple explanation for the variations in pressure sensitivity observed with several MscL mutants having either deletions or substitutions in N- or C-terminus, or site-directed mutations in the S2-S3 loop.     

The peptaibol antiamoebin as a model ion channel: similarities to bacterial potassium channels.

Antiamoebin (AAM) is a polypeptide antibiotic that is capable of forming ion channels in phospholipid membranes: planar bilayer studies have suggested the channels are octamers. The crystal structure of a monomeric form of AAM has provided the basis for molecular modelling of an octameric helical bundle channel. The channel model is funnel-shaped due to a substantial bend in the middle of the polypeptide chain caused by the presence of several imino acids. Inter-monomer hydrogen bonds orientate a ring of glutamine side chains to form a constriction in the pore lumen. The channel lumen is lined both by side chains of Gln11 and by polypeptide backbone carbonyl groups. Electrostatic calculations on the model are compatible with a channel that transports cations across membranes. The AAM channel model is compared with the crystal structures of two bacterial (KcsA andMthK) potassium channels. AAM and the potassium channels exhibit common functional features, such as cation-selectivity and similar single channel conductances. Common structural features include being multimers, each formed from a bundle of eight transmembrane helices, with lengths roughly comparable to the thickness of lipid bilayers. In addition, they all have aromatic amino acids that lie at the bilayer interfaces and which may aid in the stabilization of the transmembrane helices, as well as narrower constrictions that define the ion binding sites or selectivity filters in the pore lumen. The commonality of structural and functional features in these channels thus suggests that antiamoebin is a good, simple model for more complex bacterial and eukaryotic ion channels, capable of providing insight into details of the mechanisms of ion transport and multimeric channel stability.     

A scheme to account for the effects of Rb+ and K+ on inward rectifier K channels of bovine artery endothelial cells

An electrochemical gating model is presented to account for the effects described in the companion paper by M. R. Silver, M. S. Shapiro, and T. E. DeCoursey (1994. Journal of General Physiology, 103:519-548) of Rb+ and Rb+/K+ mixtures on the kinetics and voltage dependence of an inwardly rectifying (IR) K+ channel. The model proposes that both Rb+ and K+ act as allosteric modulators of an intrinsically voltage dependent isomerization between open and closed states. Occupancy of binding sites on the outside of the channel promotes channel opening and stabilizes the open state. Rb+ binds to separate sites within the pore and plugs IR channels. Occupancy of the pore by Rb+ can modify the rates of isomerization and the affinity of the allosteric sites for activator ions. The model also incorporates the proposed triple- barreled nature of the IR channel (Matsuda, H., 1988. Journal of Physiology. 397:237-258.) by proposing that plugging of the channel is a cooperative process involving a single site in each of the three bores, 80% of the way through the membrane field. Interaction between bores during plugging and permeation is consistent with correlated flux models of the properties of the IR channel. Parallel bores multiply the number allosteric sites associated with the macromolecular channel and allow for steep voltage dependence without compromising the parallel shift of the half-activation potential with reversal potential. Our model proposes at least six and possibly 12 such allosteric binding sites for activator ions. We derive algebraic relations that permit derivation of parameters that define simple versions of our model from the data of Silver et al. (1994). Numerical simulations based on those parameters closely reproduce that data. The model reproduces the RS+ induced slowing of IR kinetics and the negative shift of the relation between the half-activation voltage (V1/2) and reversal potential when channel plugging is associated with (a) a slowing of the isomerization rates; (b) an increase in the affinity of allosteric sites on closed channels that promote opening; and (c) a decrease in the affinity of sites on open channels that slow closing. Rb+ also slows closing at positive potentials where open channel blockade is unlikely. Allowing Rb+ to be 1.5 times more potent than K+ as an activator in the model can account for this effect and improves the match between the predicted and observed relation between the Rb+ to K+ mole fraction and the opening rate at V1/2.(ABSTRACT TRUNCATED AT 400 WORDS)     

Cysteine modification of a putative pore residue in ClC-0: implication for the pore stoichiometry of ClC chloride channels

The ClC channel family consists of chloride channels important for various physiological functions. Two members in this family, ClC-0 and ClC-1, share approximately 50-60% amino acid identity and show similar gating behaviors. Although they both contain two subunits, the number of pores present in the homodimeric channel is controversial. The double-barrel model proposed for ClC-0 was recently challenged by a one-pore model partly based on experiments with ClC-1 exploiting cysteine mutagenesis followed by modification with methanethiosulfonate (MTS) reagents. To investigate the pore stoichiometry of ClC-0 more rigorously, we applied a similar strategy of MTS modification in an inactivation-suppressed mutant (C212S) of ClC-0. Mutation of lysine 165 to cysteine (K165C) rendered the channel nonfunctional, but modification of the introduced cysteine by 2-aminoethyl MTS (MTSEA) recovered functional channels with altered properties of gating-permeation coupling. The fast gate of the MTSEA-modified K165C homodimer responded to external Cl(-) less effectively, so the P(o)-V curve was shifted to a more depolarized potential by approximately 45 mV. The K165C-K165 heterodimer showed double-barrel-like channel activity after MTSEA modification, with the fast-gating behaviors mimicking a combination of those of the mutant and the wild-type pore, as expected for the two-pore model. Without MTSEA modification, the heterodimer showed only one pore, and was easier to inactivate than the two-pore channel. These results showed that K165 is important for both the fast and slow gating of ClC-0. Therefore, the effects of MTS reagents on channel gating need to be carefully considered when interpreting the apparent modification rate.