|
||
Article |
| Discovery of Single-File Pores |
|---|
|
|
|---|
|
-helices appear poised to rigidly hold the selectivity filter structure so that it perfectly fits a K+ ion (
3 Å), but cannot collapse to fit a Na+ ion (
2 Å). This supports the classic "perfect fit" hypothesis for potassium channel selectivity (Mullins, 1959
Four decades apart, Hodgkin and Keynes (1955)
and Doyle et al. (1998)
form fundamental pillars to the field of ion channel permeation. They unequivocally demonstrate multi-ion, single-file permeation. Theoretical work between these two landmarks has explored the functional consequences of this mechanism.
| Modeling Single-File Pores |
|---|
|
|
|---|
Fig. 2 illustrates the paradigm for calculating selective flux through a single-file pore using elementary chemical kinetics. Fig. 2 A represents the hills and valleys of potential energy experienced by one kind of ion as it traverses the pore. Fig. 2 B indicates the various ways that the two energy minima of the model might be occupied (or unoccupied, 0) by two different ions, A and B. The lines show the allowed transitions between occupancy states. For example, A0 (state 2) has ion A in the left site and nothing in the right; the transition to 0A (state 4) involves the ion moving from left to right. Two rules determine which transitions can occur: (a) an ion cannot enter a site that is already occupied by another, and (b) an ion cannot pass over a site without binding to it. Transitions between sites are considered to be instantaneous.
|
The rule that transitions are instantaneous is analogous to "collision theory," the most simple-minded description of chemical reactions. Accordingly, flux of an ion over an energy barrier diminishes exponentially as the barrier height increases and is described by the rate constant: r =
(exp[–G/kT]), where [G/kT] is the height of the energy barrier (G) relative to the atom's thermal energy (kT). The frequency factor,
, is the maximum possible value of a rate constant, the value when there is no energy barrier. In single-file models, the frequency generally is set to kT/h (
6 x 1012 s–1 at biological temperatures; k and h are the Boltzmann and Planck constants, respectively, and T is absolute temperature). The rationale for this and the difficulties with it are discussed below (Limitations of Single-File Models).
The touchstone of single-file modeling studies is Hille and Schwarz (1978)
. The paper successfully mimicked a variety of unintuitive potassium channel permeation properties, including the failure of unidirectional fluxes to obey the Ussing ratio, a minimum in conductance when the relative fractions of two permeant ions are varied (the "anomalous mole fraction effect"), and rectification due to blocking ions. Hille and Schwarz (1978)
found that repulsion between intrapore ions had profound effects on the calculated flux and concluded that two properties of the model were necessary to reproduce physiological phenomena. First, external energy barriers had to be high compared with internal ones. This allows ions to equilibrate readily between intrapore sites—to rattle about the pore more easily than they can exit it. Second, the pore had to have at least three sites and contain at least two ions most of the time.
Now, Doyle et al. (1998)
find two ions in the potassium channel selectivity filter sitting at opposite extremes of the filter's length. This image is consistent with repulsion pushing the two ions to the extremes of a pore from which exit is more difficult than is intrapore mobility. Thus, the image is consistent with the physical principles suggested by Hille and Schwarz (1978)
, as well as their number of ions. This is the value of single-file calculations. Although the models must be grossly simplistic representations of permeation through ion channels, they still provide a means to detect and test the basic principles involved. Other important single-file papers have described permeation and selectivity in gramicidin channels (Andersen and Procopio, 1980
), sodium channels (Begenisich, 1987
), and calcium channels (below).
| Calcium Channel Models |
|---|
|
|
|---|
1 µM. When the high affinity site is unoccupied by Ca2+, Na+ is freely permeant; when Ca2+ occupies the site, as it must at physiological Ca2+ concentrations, Na+ current through the channel is blocked. Thus, Ca2+ confers selectivity on calcium channels because it blocks permeation of competing ions. Similarly, Na+ current through potassium channels is blocked by K+ (Korn and Ikeda, 1995
Discovery of high affinity binding in the calcium channel raised a fundamental issue. Since Kd = koff/kon, and since the maximum kon is the diffusion limit of 109/M–1s–1 (Almers and McCleskey, 1984
), the maximum off rate from a 1-µM binding site should be only 1,000/s. However, picoampere currents (millions of ions per second) pass through calcium channels. A goal of permeation models has been to suggest mechanisms for solving this 1,000-fold discrepancy between expected and observed flux.
This "sticky pore problem" is diagrammed in Fig. 3 A. Binding affinity is determined by the difference in energy between the outside environment and the well (arrow 1); flux is limited by the height of the greatest barrier to exit from the pore (arrow 2). In such a single site pore, tighter binding (a longer arrow 1) necessarily causes slower flux (a longer arrow 2); this makes selectivity by selective binding impossible. Models of multi-ion, single-file calcium channels have deduced two separate solutions to this problem, although the two mechanisms might well work together in a real pore. The repulsion model (Fig. 3 B) imagined two high affinity sites within the pore and allowed ions to repel each other when they occupied the pore simultaneously (Almers and McCleskey, 1984
; Hess and Tsien, 1984
). Foreign ions are blocked at micromolar Ca2+ concentrations when a single Ca2+ ion occupies the pore (left). When two Ca2+ ions occupy the pore, they repel each other, and this effectively diminishes the pore's affinity for Ca2+ (apparent Kd changes from
1 µM to 10 mM, right). Once the pore is in this low affinity state, Ca2+ flux occurs readily. Recent mutation studies (see next section) cleanly disprove one aspect of this model because they demonstrate only a single high affinity site in the calcium channel. Nevertheless, the more fundamental point—that flux occurs because multiple occupancy lowers the affinity of the pore for Ca2+—remains reasonable.
|
The PNP model of calcium channels (Nonner and Eisenberg, 1998
) provides another major, and unexpected, insight. The model's pore is an electronegative tube that is 10-Å long and 6-Å wide. Despite the absence of an explicit high affinity binding site, the model exhibits high affinity Ca2+ block of Na+ current as well as high Ca2+ flux. This discovery is very timely: it shows that a rigid, diffusely negative tube (much like the crystal structure in Doyle et al., 1998
) is sufficient to generate selectivity through high affinity binding. The pore binds Ca2+ without ligands that wrap around the ion, as is generally assumed necessary for high affinity binding; surely, this has great significance for ion channel permeation.
In an important sense, the differences between the single-file and PNP models have been overstated. The PNP/calcium channel demonstrates an unexpected mechanism for intrapore, high affinity binding. Single-file models do not, and cannot, address the mechanism or structures of binding (see Limitations of Single-File Modeling). Rather, they focus on the forces that allow rapid exit from the high affinity site. Two forces were described and here the models agree: both repulsion and potential energy steps contribute to rapid efflux from the PNP/calcium channel pore according to the authors. It is the mechanism of flux that seems so different in the two models, but this may largely be a matter of mathematical appearances. As discussed above, both types of models envision ion movement to be similar to that through crystals, but one describes this with discrete steps and the other with continuous functions. If the models are ever complete (and neither pretends to be), the two approaches must give identical results just as diffusion is described equally well through random walk or continuum equations (Einstein, 1908
). For the time being, it is encouraging that such different approaches agree in the description of forces that underlay high flux out of high affinity pores.
Like every other model, the PNP/calcium channel comes with ad hoc assumptions. A critical one necessary to allow Ca2+ to block Na+ current is that the pore has an excess chemical potential that is negative for Ca2+ and positive for Na+; in other words, some undefined chemistry of the model pore is attractive for Ca2+ but repulsive for Na+. What chemistry could attract one positive ion while repelling another of the same size is unclear. Another issue, discussed below, is an explicit disagreement with observations made in mutation studies. Such failings confirm only that the PNP model is an incomplete approximation of the calcium channel. But, by definition, all models are incomplete. In my view, such specific inconsistencies should not detract from the more general discovery inherent in the PNP/ calcium channel: that selectivity by high affinity binding can occur in a rigid pore.
| Calcium Channel Facts |
|---|
|
|
|---|
500 amino acids of sequence, these glutamates are in the putative pore-lining region of the four internally conserved repeats. In the analogous sodium channel sequence, two of the negative charges are present, but a neutral and a positive residue substitute at the other two sites. Mutation of these two residues to glutamates made the sodium channel behave like a calcium channel in that Ca2+ could block Na+ current at low concentrations. The four glutamate residues in calcium channels are called the EEEE locus.
Yang et al. (1993)
mutated the EEEE locus in L-type calcium channels, and thereby demonstrated the role of these glutamates in high affinity Ca2+ binding (Fig. 4 A). Neutralizing any of the glutamates to a glutamine shifted the block of monovalent current to higher Ca2+ concentrations. The extent of the shift differed for glutamates in the different repeats, and neutralizing multiple residues shifted the blocking concentration by as much as 1,000-fold (Ellinor et al., 1995
). Because a single point mutation can shift the dissociation constant for Ca2+, there must not be two independent high affinity sites in the pore. Moreover, pairwise mutations of the glutamates rigorously refute the two sites postulated by the repulsion model (Ellinor et al., 1995
).
|
These mutation studies potently argue against the PNP model for the calcium channel. The PNP model pore contains four carboxyl groups that are meant to mimic the EEEE locus. However, there is no difference in the concentration at which Ca2+ blocks Na+ current in the PNP/calcium channel when the number of intrapore carboxyl residues is set to 1, 2, 3, or 4 (see Figure 6 A in Nonner and Eisenberg, 1998
). Thus, the PNP model starkly fails to reproduce the data in Fig. 4 A. Whatever may be the actual structure of the EEEE locus, the four carboxyls in the PNP model do not mimic it.
| Limitations of Single-File Modeling |
|---|
|
|
|---|
The Appendix in Nonner and Eisenberg (1998)
takes a long step further in criticizing the single-file approach. It uses the framework of the PNP model to recalculate flux through the step model and finds that: (a) Ca2+ flux through this step-PNP model is many hundred-fold less than experimental observations, and (b) the pore fails to select Ca2+ over Na+. The results are not shown for other models, but it is stated that similar PNP recalculations of the repulsion models are "even less consistent" with experimental observations, and that "the reasons for failure are generic for rate theory models in the tradition of Hille and Schwarz (1978)
." The problem, it says, is that single-file models use a frequency factor for calculating rate constants that is 10,000-fold too high. If so, all single-file models are incorrect by many orders of magnitude, implying that their insights are based on false physics and should be abandoned.
The primary problem with the Appendix in Nonner and Eisenberg (1998)
is that it fails to recognize that frequency factors and barrier heights are not treated as separable entities in single-file models. For example, our models set the external barriers (the values that limit entry and exit of the pore and are therefore the most critical) so that they agree with experimental observation. Lansman et al. (1986)
demonstrated that Ca2+ enters the pore at a rate of
109 M–1s–1, the value expected from diffusion-limited access. Assuming a prefactor of 6 x 1012 s–1, the external barriers are set no lower than 8.7 kT, thereby creating an entry rate of 109 M–1s–1 (Almers and McCleskey, 1984
). If we had used a lower prefactor, we would have set the external barrier energies lower, always keeping the rate constant (the product of the prefactor and the exponential) consistent with experimental observations. If one diminishes the prefactor by 10,000 while keeping the barriers the same, the entry rate will be 10,000-fold less than experimental observation. In this trivial way, calculated flux would be unphysiologically low. This is the essence of the error in the Appendix in Nonner and Eisenberg (1998)
: it uses precise energies from one calculation scheme while applying assumptions from another.
Although the appendix's calculation is artificial (Nonner and Eisenberg, 1998
), the point about frequency factors is important. To understand the flaw in setting the frequency factor,
, equal to kT/h, one should understand its basis. kT/h is an atom's theoretical vibration frequency, obtained by assuming that all thermal energy (kT) becomes kinetic, or vibrational, energy (h
). The assumption confers a simple physical interpretation on the equation for a rate constant: the atom can attempt to surpass the reaction barrier each time that it vibrates, but the probability of a successful attempt is very low. The attempt rate is kT/h, about six times per picosecond, and the probability of success is exp[–G/kT]. The problem with this is nicely explained by Andersen and Koeppe (1992)
, and by Andersen's introduction to these Perspectives. In short, the frequency of vibration is only legitimate to use if the distance to surpass the barrier is similar to the length of a vibration, and this is just a small fraction of an angstrom. To traverse a selectivity filter that is many angstroms long, many such individual steps are necessary.
Assuming that movement within the filter is diffusive, these individual steps would have low energy barriers and transit would be fast compared with exit from the pore. Therefore, most single-file models simply lump these many intrapore steps into one or two individual steps that are not rate limiting. The models thereby focus on the high barriers that limit exit from the pore; these external barriers, because their size is determined by experimental measurements (see above), are handled appropriately by the models. However, there is no denying that lumping the internal steps together represents the pore in a physically inaccurate manner.
A perspective on the insights from single-file models arises from considering the ultimate aims of the permeation field. A major one is to construct an accurate map of the chemical potential energy of an ion as it traverses a pore of precisely known structure. This is analogous to converting photographs of a hillside into a topographic map, which describes your gravitational potential energy as you hike the hill. Single-file models have skipped directly to the topographic map without having the photographs. Unlinked to real structures, these energy diagrams are nothing more than quantitative cartoons, but the quantitation has been valuable. The models do not pretend to tell us about protein structure, but, by finding which forces allow calculated fluxes to mimic real experiments, they suggest general principles relevant to permeation. These principles include: the importance of having multiple ions moving in a coordinated fashion through a pore, that ions likely equilibrate within the pore far more easily than they exit it, that channels can select their ions through binding sites and still have high flux, and that exit from a tightly binding pore is promoted through either ion– ion interactions in a multiply occupied pore or by a series of lower affinity sites.
| What Are the Real Questions? |
|---|
|
|
|---|
(1) How can Ca2+ have diffusion-limited access to the pore? Three decades ago, Eigen and Winkler (1971)
noted that ions enter carriers at the diffusion limit; Lansman et al. (1986)
demonstrated the same for calcium channels. The high rate is unexpected: ions must completely dehydrate to enter pores and dehydration requires so much energy that it should occur only rarely. In addition to replacing an ion's hydrating waters with solvating intrapore ligands, channels and carriers must also rapidly catalyze the dehydration process. How is this done? Eigen and Winkler (1971)
argue that rapid dehydration requires systematic, stepwise replacement of water molecules by substituting ligands, rather than a loss of all waters at once. Andersen and Koeppe (1992)
note that this requires a pore entrance that is not rigid. What structures at the mouth of the pore are responsible for stepwise dehydration, and over how long a distance does this process occur?
(2) Is the strong binding of the selectivity filter of the calcium channel due to a flexible binding site analogous to Ca2+ chelators? The calcium channel mutation studies of Yang et al. (1993)
and Ellinor et al. (1995)
are most easily explained this way. However, the potassium channel selectivity filter described by Doyle et al. (1998)
is a rigid pore and the PNP model for calcium channels (Nonner and Eisenberg, 1998
) shows that high affinity binding can occur in such a structure. Is the selectivity filter of calcium channels rigid or flexible? Is it lined by backbone carbonyl oxygens or by amino acid side chains?
(3) What forces allow ions to exit this high affinity binding area so much faster than expected? Two have been suggested in various models: ion–ion interactions and stepwise increases in energy during exit. What might cause ion–ion interactions in the real pore— competition for a binding site and/or electrostatic repulsion? What might cause stairsteps of energy in the real pore—specific binding structures of successively lower affinity and/or stepwise rehydration? What is the relative importance of these two forces, and are there others that have not previously been considered?
| references |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
X. Gong and P. Linsdell Coupled movement of permeant and blocking ions in the CFTR chloride channel pore J. Physiol., June 1, 2003; 549(2): 375 - 385. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gong and P. Linsdell Molecular determinants and role of an anion binding site in the external mouth of the CFTR chloride channel pore J. Physiol., June 1, 2003; 549(2): 387 - 397. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Talavera, A. Janssens, N. Klugbauer, G. Droogmans, and B. Nilius Pore Structure Influences Gating Properties of the T-type Ca2+ Channel {alpha}1G J. Gen. Physiol., May 27, 2003; 121(6): 529 - 540. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. S. L. Yu, A. H. Enck, W. I. Lencer, and E. E. Schneeberger Claudin-8 Expression in Madin-Darby Canine Kidney Cells Augments the Paracellular Barrier to Cation Permeation J. Biol. Chem., May 2, 2003; 278(19): 17350 - 17359. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Barreiro, C. R. W. Guimaraes, and R. B. de Alencastro Potential of mean force calculations on an L-type calcium channel model Protein Eng. Des. Sel., March 1, 2003; 16(3): 209 - 215. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Perez-Reyes Molecular Physiology of Low-Voltage-Activated T-type Calcium Channels Physiol Rev, January 1, 2003; 83(1): 117 - 161. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Watanabe, C. Beck, T. Kuner, L. S. Premkumar, and L. P. Wollmuth DRPEER: A Motif in the Extracellular Vestibule Conferring High Ca2+ Flux Rates in NMDA Receptor Channels J. Neurosci., December 1, 2002; 22(23): 10209 - 10216. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Yue, B. Navarro, D. Ren, A. Ramos, and D. E. Clapham The Cation Selectivity Filter of the Bacterial Sodium Channel, NaChBac J. Gen. Physiol., November 25, 2002; 120(6): 845 - 853. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Cataldi, E. Perez-Reyes, and R. W. Tsien Differences in Apparent Pore Sizes of Low and High Voltage-activated Ca2+ Channels J. Biol. Chem., November 22, 2002; 277(48): 45969 - 45976. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Barreiro, C. R. W. Guimaraes, and R. B. de Alencastro A molecular dynamics study of an L-type calcium channel model Protein Eng. Des. Sel., February 1, 2002; 15(2): 109 - 122. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-S. Fan, Y. Yuan, and P. Palade FPL-64176 modifies pore properties of L-type Ca2+ channels Am J Physiol Cell Physiol, March 1, 2001; 280(3): C565 - C572. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. K. Cloues, S. M. Cibulsky, and W. A. Sather Ion Interactions in the High-Affinity Binding Locus of a Voltage-Gated Ca2+ Channel J. Gen. Physiol., October 1, 2000; 116(4): 569 - 586. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Cibulsky and W. A. Sather The Eeee Locus Is the Sole High-Affinity Ca2+ Binding Structure in the Pore of a Voltage-Gated Ca2+ Channel: Block by Ca2+ Entering from the Intracellular Pore Entrance J. Gen. Physiol., September 1, 2000; 116(3): 349 - 362. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. McCarty Permeation through the CFTR chloride channel J. Exp. Biol., January 7, 2000; 203(13): 1947 - 1962. [Abstract] [PDF] |
||||
![]() |
Forces Determining Ion Permeation J. Gen. Physiol., October 1, 1999; 114(4): 593 - 596. [Full Text] [PDF] |
||||
![]() |
Ion Permeation and Chemical Kinetics J. Gen. Physiol., October 1, 1999; 114(4): 601 - 604. [Full Text] [PDF] |
||||
![]() |
O. Mignen and T. J. Shuttleworth Permeation of Monovalent Cations through the Non-capacitative Arachidonate-regulated Ca2+ Channels in HEK293 Cells. COMPARISON WITH ENDOGENOUS STORE-OPERATED CHANNELS J. Biol. Chem., June 8, 2001; 276(24): 21365 - 21374. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|