Molecular docking of the scorpion toxin Tc1 to the structural model of the voltage-gated potassium channel Kv1.1 from human Homo sapiens

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K+ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

Molecular Docking of the Scorpion Toxin Tc1 to the Structural Model of the Voltage-gated Potassium Channel Kv1.1 from Human Homo sapiens

Abstract

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K⁺ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

KcsA crystal structure as framework for a molecular model of the Na(+) channel pore

The crystal structure of the pore-forming part of the KcsA bacterial K(+)-selective channel suggests a possible motif for related voltage-gated channels. We examined the hypothesis that the spacial orientation of the KcsA M1 and M2 alpha-helices also predicts the backbone location of S5 and S6 helices of the voltage-gated Na(+) channel. That channel's P region structure is expected to be different because selectivity is determined by side-chain interactions rather than by main-chain carbonyls, and its outer vestibule accommodates relatively large toxin molecules, tetrodotoxin (TTX) and saxitoxin (STX), which interact with selectivity ring residues. The Na(+) channel P loop was well-modeled by the alpha-helix-turn-beta-strand motif, which preserves the relationships for toxin interaction with the Na(+) channel found experimentally. This outer vestibule was docked into the extracellular part of the inverted teepee structure formed by the S5 and S6 helices that were spacially located by coordinates of the KcsA M1 and M2 helix main chains [Doyle et al. (1998) Science 280, 69-74], but populated with side chains of the respective S5 and S6 structures. van der Waals contacts were optimized with minimal adjustment of the S5, S6, and P loop structures, forming a densely packed pore structure. Nonregular external S5-P and P-S6 segments were not modeled here, except the P-S6 segment of domain II. The resulting selectivity region structure is consistent with Na(+) channel permeation properties, offering suggestions for the molecular processes involved in selectivity. The ability to construct a Na(+) channel pore model consistent with most of the available biophysical and mutational information suggests that the KcsA structural framework may be conserved in voltage-gated channels.     

Probing the outer mouth structure of the HERG channel with peptide toxin footprinting and molecular modeling

Previous studies have shown that the unusually long S5-P linker lining human ether a-go-go related gene's (hERG's) outer vestibule is critical for its channel function: point mutations at high-impact positions here can interfere with the inactivation process and, in many cases, also reduce the pore's K+ selectivity. Because no data are available on the equivalent region in the available K channel crystal structures to allow for homology modeling, we used alternative approaches to model its three-dimensional structure. The first part of this article describes mutant cycle analysis used to identify residues on hERG's outer vestibule that interact with specific residues on the interaction surface of BeKm-1, a peptide toxin with known NMR structure and a high binding affinity to hERG. The second part describes molecular modeling of hERG's pore domain. The transmembrane region was modeled after the crystal structure of KvAP pore domain. The S5-P linker was docked to the transmembrane region based on data from previous NMR and mutagenesis experiments, as well as a set of modeling criteria. The models were further restrained by contact points between hERG's outer vestibule and the bound BeKm-1 toxin molecule deduced from the mutant cycle analysis. Based on these analyses, we propose a working model for the open conformation of the outer vestibule of the hERG channel, in which the S5-P linkers interact with the pore loops to influence ion flux through the pore.     

Structural similarities between glutamate receptor channels and K(+) channels examined by scanning mutagenesis

The pores of glutamate receptors and K(+) channels share sequence homology, suggesting a conserved secondary structure. Scanning mutagenesis with substitution of alanine and tryptophan in GluR6 channels was performed based on the structure of KcsA. Our assay used disruption of voltage-dependent polyamine block to test for changes in the packing of pore-forming regions. Alanine scanning from D567 to R603 revealed reduced rectification resulting from channel block in two regions. A periodic pattern from F575 to M589 aligned with the pore helix in KcsA, whereas a cluster of sensitive positions around Q590, a site regulated by RNA editing, mapped to the selectivity filter in KcsA. Tryptophan scanning from D567 to R603 revealed similar patterns, but with a complete disruption of spermine block for 7 out of the 37 positions and a pM dissociation constant for Q590W. Molecular modeling with KcsA coordinates showed that GluR6 pore helix mutants disrupting polyamine block pack against M1 and M2, and are not exposed in the ion channel pore. In the selectivity filter, tryptophan creates an aromatic cage consistent with the pM dissociation constant for Q590W. A scan with glutamate substitution was used to map the cytoplasmic entrance to the pore based on charge neutralization experiments, which established that E594 was uniquely required for high affinity polyamine block. In E594Q mutants, introduction of glutamate at positions S593-L600 restored polyamine block at positions corresponding to surface-exposed residues in KcsA. Our results reinforce proposals that the pore region of glutamate receptors contains a helix and pore loop analogous to that found in K(+) channels. At the cytoplasmic entrance of the channel, a negatively charged amino acid, located in an extended loop with solvent-exposed side chains, is required for high affinity polyamine block and probably attracts cations via a through space electrostatic mechanism.     

Light at the end of the Ca(2+)-release channel tunnel: structures and mechanisms involved in ion translocation in ryanodine receptor channels

RyR and InsP3R are Ca(2+)-release channels. When induced to open by the appropriate stimulus, these channels allow Ca2+ to leave intracellular storage organelles at an astonishing rate. Investigations of the ion-handling properties of isolated RyR channels have demonstrated that, at least in comparison to voltage-gated channels of surface membranes, these channels display limited powers of discrimination between physiologically relevant cations and this relative lack of selectivity is likely to contribute to the ability of Ca(2+)-release channels to maintain high rates of cation translocation without compromising function. A range of ion-handling properties in RyR are consistent with the proposal that this channel functions as a single-ion channel and theoretical considerations indicate that the high rates of ion translocation monitored for RyR would require the pore of such a structure to be short and possess a large capture radius. Measurements of the dimensions of regions of RyR involved in ion conduction and discrimination indicate that this is likely to be the case. In each monomer of RyR/InsP3R, residues making up the last two trans-membrane spanning domains and a luminal loop linking these two helices contribute to the formation of the channel pore. The luminal loops of both RyR and InsP3R contain amino acid sequences similar to those known to form the selectivity filter of K+ channels. In addition the luminal loops of both Ca(2+)-release channels contain sequences that are likely to form helices that may be analogous to the pore helix visualised in KcsA. The correlation in structural elements of the luminal loops of RyR/InsP3R and KcsA has prompted us to speculate on the tertiary arrangement for this region of the Ca(2+)-release channels using the established structure of KcsA as a framework.     

Probing the role of negatively charged amino acid residues in ion permeation of skeletal muscle ryanodine receptor

Sequence comparison suggests that the ryanodine receptors (RyRs) have pore architecture similar to that of the bacterial K+ channel KcsA. The lumenal loop linking the two most C-terminal transmembrane spanning segments in the RyRs has a predicted pore helix and an amino acid motif (GGGIG) similar to the selectivity filter (TVGYG) of KcsA identified by x-ray analysis. The RyRs have many negatively charged amino acid residues in the two regions linking the GGGIG motif and predicted pore helix with the two most C-terminal transmembrane spanning segments. We tested the role of these residues by generating single-site mutants, focusing on amino acid residues conserved among the mammalian RyRs. Replacement of two acidic residues immediately after the GGGIG motif in skeletal muscle ryanodine receptor (RyR1-D4899 and -E4900) with asparagine and glutamine profoundly affected ion permeation and selectivity. By comparison, mutagenesis of aspartate and glutamate residues in the putative linker regions showed a K+ conductance and selectivity for Ca2+ compared to K+ (P(Ca)/P(K)) close to wild-type. The results show that the negatively charged carboxyl oxygens of D4899 and E4900 side chains are major determinants of RyR ion conductance and selectivity.     

Models of ion pores in N-type voltage-gated calcium channels

Two computer models of the outer vestibule of the pore of the N-type voltage-gated Ca²⁺ channel are predicted. The models are constructed from β-hairpin peptide segments in the S5–S6 loops of each of the four domains that produce the channel. These hairpins together are modeled to form a short eight-stranded β barrel. The models contain a ring of glutamates at the base of the barrel, which have been shown by mutagenesis experiments to function as a selectivity filter. These filters are suggested by the models to be of the correct dimensions to allow the permeation of a hydrated calcium ion, where the filter glutamates may substitute for molecules of water from the hydration shell of the ion. The models also suggest that a ring of threonines and an aspartate might be present between the mouth of the pore and the filter, and hence the models may prove useful in suggesting future mutagenesis experiments.     

A model of sodium channels

We have explored the permeation and blockage of ions in sodium channels, relating the channel structure to function using electrostatic profiles and Brownian dynamics simulations. The model used resembles the KcsA potassium channel with an added external vestibule and a shorter selectivity filter. The electrostatic energy landscape seen by permeating ions is determined by solving Poisson's equation. The two charged amino acid rings of Glu-Glu-Asp-Asp (EEDD) and Asp-Glu-Lys-Ala (DEKA) around the selectivity filter region are seen to play a crucial role in making the channel sodium selective, and strongly binding calcium ions such that they block the channel. Our model closely reproduces a range of experimental data including the current-voltage curves, current-concentration curves and blockage of monovalent ions by divalent ions.     

Potassium and sodium binding to the outer mouth of the K+ channel.

Molecular dynamics simulations of the K+ channel from Streptomyces lividans (KcsA channel) were performed in a membrane-mimetic environment with Na+ and K+ in different initial locations. The structure of the channel remained stable and well preserved for simulations lasting up to 1.5 ns. Salt bridges between Asp80 and Arg89 of neighboring subunits, not detected in the X-ray structure, enhanced the stability of the tetrameric structure. Na+ or K+ ions located in the channel vestibule lost part of their hydration shell and diffused into the channel inner pore in less than a few hundred picoseconds. This powerful catalytic action was caused by strong electrostatic interactions with Asp80 and Glu71. The hydration state of the metal ions turned out to depend significantly on the conformational flexibility of the channel. Furthermore, Na+ entered the channel inner pore bound to more water molecules than K+. The different hydration state of the two ions may be a determinant factor in the ion selectivity of the channel.     

Influence of pore residues on permeation properties in the Kv2.1 potassium channel. Evidence for a selective functional interaction of K+ with the outer vestibule

The Kv2.1 potassium channel contains a lysine in the outer vestibule (position 356) that markedly reduces open channel sensitivity to changes in external [K(+)]. To investigate the mechanism underlying this effect, we examined the influence of this outer vestibule lysine on three measures of K(+) and Na(+) permeation. Permeability ratio measurements, measurements of the lowest [K(+)] required for interaction with the selectivity filter, and measurements of macroscopic K(+) and Na(+) conductance, were all consistent with the same conclusion: that the outer vestibule lysine in Kv2.1 interferes with the ability of K(+) to enter or exit the extracellular side of the selectivity filter. In contrast to its influence on K(+) permeation properties, Lys 356 appeared to be without effect on Na(+) permeation. This suggests that Lys 356 limited K(+) flux by interfering with a selective K(+) binding site. Combined with permeation studies, results from additional mutagenesis near the external entrance to the selectivity filter indicated that this site was located external to, and independent from, the selectivity filter. Protonation of a naturally occurring histidine in the same outer vestibule location in the Kv1.5 potassium channel produced similar effects on K(+) permeation properties. Together, these results indicate that a selective, functional K(+) binding site (e.g., local energy minimum) exists in the outer vestibule of voltage-gated K(+) channels. We suggest that this site is the location of K(+) hydration/dehydration postulated to exist based on the structural studies of KcsA. Finally, neutralization of position 356 enhanced outward K(+) current magnitude, but did not influence the ability of internal K(+) to enter the pore. These data indicate that in Kv2.1, exit of K(+) from the selectivity filter, rather than entry of internal K(+) into the channel, limits outward current magnitude. We discuss the implications of these findings in relation to the structural basis of channel conductance in different K(+) channels.     

A model of sodium channels

We have explored the permeation and blockage of ions in sodium channels, relating the channel structure to function using electrostatic profiles and Brownian dynamics simulations. The model used resembles the KcsA potassium channel with an added external vestibule and a shorter selectivity filter. The electrostatic energy landscape seen by permeating ions is determined by solving Poisson's equation. The two charged amino acid rings of Glu-Glu-Asp-Asp (EEDD) and Asp-Glu-Lys-Ala (DEKA) around the selectivity filter region are seen to play a crucial role in making the channel sodium selective, and strongly binding calcium ions such that they block the channel. Our model closely reproduces a range of experimental data including the current-voltage curves, current-concentration curves and blockage of monovalent ions by divalent ions.     

A model of the putative pore region of the cardiac ryanodine receptor channel

Using the bacterial K+ channel KcsA as a template, we constructed models of the pore region of the cardiac ryanodine receptor channel (RyR2) monomer and tetramer. Physicochemical characteristics of the RyR2 model monomer were compared with the template, including homology, predicted secondary structure, surface area, hydrophobicity, and electrostatic potential. Values were comparable with those of KcsA. Monomers of the RyR2 model were minimized and assembled into a tetramer that was, in turn, minimized. The assembled tetramer adopts a structure equivalent to that of KcsA with a central pore. Characteristics of the RyR2 model tetramer were compared with the KcsA template, including average empirical energy, strain energy, solvation free energy, solvent accessibility, and hydrophobic, polar, acid, and base moments. Again, values for the model and template were comparable. The pores of KcsA and RyR2 have a common motif with a hydrophobic channel that becomes polar at both entrances. Quantitative comparisons indicate that the assembled structure provides a plausible model for the pore of RyR2. Movement of Ca2+, K+, and tetraethylammonium (TEA+) through the model RyR2 pore were simulated with explicit solvation. These simulations suggest that the model RyR2 pore is permeable to Ca2+ and K+ with rates of translocation greater for K+. In contrast, simulations indicate that tetraethylammonium blocks movement of metal cations.     

Mechanism of Cd(2+) Coordination during Slow Inactivation in Potassium Channels

In K(+) channels, rearrangements of the pore outer vestibule have been associated with C-type inactivation gating. Paradoxically, the crystal structure of Open/C-type inactivated KcsA suggests these movements to be modest in magnitude. In this study, we show that under physiological conditions, the KcsA outer vestibule undergoes relatively large dynamic rearrangements upon inactivation. External Cd(2+) enhances the rate of C-type inactivation in an cysteine mutant (Y82C) via metal-bridge formation. This effect is not present in a non-inactivating mutant (E71A/Y82C). Tandem dimer and tandem tetramer constructs of equivalent cysteine mutants in KcsA and Shaker K(+) channels demonstrate that these Cd(2+) metal bridges are formed only between adjacent subunits. This is well supported by molecular dynamics simulations. Based on the crystal structure of Cd(2+)-bound Y82C-KcsA in the closed state, together with electron paramagnetic resonance distance measurements in the KcsA outer vestibule, we suggest that subunits must dynamically come in close proximity as the channels undergo inactivation.     

Exploring the open pore of the potassium channel from Streptomyces lividans

The tetrameric potassium channel from Streptomyces lividans (KcsA) embedded in planar bilayers exhibits the following electrophysiological characteristics: (i) K+ ions can cross the pore in a highly hydrated state (nH2O > or = 6), (ii) the selectivity for K+ exceeds that for Na+ ions by 11 times, and both Ca2+ and Mg2+ are permeant, (iii) the internal side is blocked by Ba2+ ions in a voltage-dependent manner, (iv) intrinsic rectification is due to gating, depending on the direction of the electric field, (v) the internal side is pH-sensitive, and (vi) the open pore has a diameter of approximately 5.8 A. In conclusion, our results show that ion conduction and selectivity of KcsA cannot easily be reconciled with the properties deduced from the rigid crystal structure [Doyle et al., Science 280 (1998) 69-77], which must be concluded to have the pore trapped in its closed state.     

Structural conservation of the pores of calcium-activated and voltage-gated potassium channels determined by a sea anemone toxin.

The structurally defined sea anemone peptide toxins ShK and BgK potently block the intermediate conductance, Ca(2+)-activated potassium channel IKCa1, a well recognized therapeutic target present in erythrocytes, human T-lymphocytes, and the colon. The well characterized voltage-gated Kv1.3 channel in human T-lymphocytes is also blocked by both peptides, although ShK has a approximately 1,000-fold greater affinity for Kv1.3 than IKCa1. To gain insight into the architecture of the toxin receptor in IKCa1, we used alanine-scanning in combination with mutant cycle analyses to map the ShK-IKCa1 interface, and compared it with the ShK-Kv1.3 interaction surface. ShK uses the same five core residues, all clustered around the critical Lys(22), to interact with IKCa1 and Kv1.3, although it relies on a larger number of contacts to stabilize its weaker interactions with IKCa1 than with Kv1.3. The toxin binds to IKCa1 in a region corresponding to the external vestibule of Kv1.3, and the turret and outer pore of the structurally defined bacterial potassium channel, KcsA. Based on the NMR structure of ShK, we deduce the toxin receptor in IKCa1 to have x-y dimensions of approximately 22 A, a diameter of approximately 31 A, and a depth of approximately 8 A; we estimate that the ion selectivity lies approximately 13 A below the outer lip of the toxin receptor. These dimensions are in good agreement with those of the KcsA channel determined from its crystal structure, and the inferred structure of Kv1.3 based on mapping with scorpion toxins. Thus, these distantly related channels exhibit architectural similarities in the outer pore region. This information could facilitate development of specific and potent modulators of the therapeutically important IKCa1 channel.     

Molecular architecture of the voltage-dependent Na channel: functional evidence for alpha helices in the pore

The permeation pathway of the Na channel is formed by asymmetric loops (P segments) contributed by each of the four domains of the protein. In contrast to the analogous region of K channels, previously we (Yamagishi, T., M. Janecki, E. Marban, and G. Tomaselli. 1997. Biophys. J. 73:195-204) have shown that the P segments do not span the selectivity region, that is, they are accessible only from the extracellular surface. The portion of the P-segment NH(2)-terminal to the selectivity region is referred to as SS1. To explore further the topology and functional role of the SS1 region, 40 amino acids NH(2)-terminal to the selectivity ring (10 in each of the P segments) of the rat skeletal muscle Na channel were substituted by cysteine and expressed in tsA-201 cells. Selected mutants in each domain could be blocked with high affinity by externally applied Cd(2)+ and were resistant to tetrodotoxin as compared with the wild-type channel. None of the externally applied sulfhydryl-specific methanethiosulfonate reagents modified the current through any of the mutant channels. Both R395C and R750C altered ionic selectivity, producing significant increases in K(+) and NH(4)(+) currents. The pattern of side chain accessibility is consistent with a pore helix like that observed in the crystal structure of the bacterial K channel, KcsA. Structure prediction of the Na channel using the program PHDhtm suggests an alpha helix in the SS1 region of each domain channel. We conclude that each of the P segments undergoes a hairpin turn in the permeation pathway, such that amino acids on both sides of the putative selectivity filter line the outer mouth of the pore. Evolutionary conservation of the pore helix motif from bacterial K channels to mammalian Na channels identifies this structure as a critical feature in the architecture of ion selective pores.     

Sequence-function analysis of the K+-selective family of ion channels using a comprehensive alignment and the KcsA channel structure

Sequence-function analysis of K(+)-selective channels was carried out in the context of the 3.2 A crystal structure of a K(+) channel (KcsA) from Streptomyces lividans (Doyle et al., 1998). The first step was the construction of an alignment of a comprehensive set of K(+)-selective channel sequences forming the putative permeation path. This pathway consists of two transmembrane segments plus an extracellular linker. Included in the alignment are channels from the eight major classes of K(+)-selective channels from a wide variety of species, displaying varied rectification, gating, and activation properties. Segments of the alignment were assigned to structural motifs based on the KcsA structure. The alignment's accuracy was verified by two observations on these motifs: 1), the most variability is shown in the turret region, which functionally is strongly implicated in susceptibility to toxin binding; and 2), the selectivity filter and pore helix are the most highly conserved regions. This alignment combined with the KcsA structure was used to assess whether clusters of contiguous residues linked by hydrophobic or electrostatic interactions in KcsA are conserved in the K(+)-selective channel family. Analysis of sequence conservation patterns in the alignment suggests that a cluster of conserved residues is critical for determining the degree of K(+) selectivity. The alignment also supports the near-universality of the "glycine hinge" mechanism at the center of the inner helix for opening K channels. This mechanism has been suggested by the recent crystallization of a K channel in the open state. Further, the alignment reveals a second highly conserved glycine near the extracellular end of the inner helix, which may be important in minimizing deformation of the extracellular vestibule as the channel opens. These and other sequence-function relationships found in this analysis suggest that much of the permeation path architecture in KcsA is present in most K(+)-selective channels. Because of this finding, the alignment provides a robust starting point for homology modeling of the permeation paths of other K(+)-selective channel classes and elucidation of sequence-function relationships therein. To assay these applications, a homology model of the Shaker A channel permeation path was constructed using the alignment and KcsA as the template, and its structure evaluated in light of established structural criteria.     

Ion binding to KcsA: implications in ion selectivity and channel gating

Binding of K+ and Na+ to the potassium channel KcsA has been characterized from the stabilization observed in the heat-induced denaturation of the protein as the ion concentration is increased. KcsA thermal denaturation is known to include (i) dissociation of the homotetrameric channel into its constituent subunits and (ii) protein unfolding. The ion concentration-dependent changes in the thermal stability of the protein, evaluated as the Tm value for thermal-induced denaturation of the protein, may suggest the existence of both high- and low-affinity K+ binding sites of KcsA, which lend support to the tenet that channel gating may be governed by K+ concentration-dependent transitions between different affinity states of the channel selectivity filter. We also found that Na+ binds to KcsA with a KD similar to that estimated electrophysiologically from channel blockade. Therefore, our findings on ion binding to KcsA partly account for K+ over Na+ selectivity and Na+ blockade and argue against the strict “snug fit” hypothesis used initially to explain ion selectivity from the X-ray channel structure. Furthermore, the remarkable effects of increasing the ion concentration, K+ in particular, on the Tm of the denaturation process evidence that synergistic effects of the metal-mediated intersubunit interactions at the channel selectivity filter are a major contributor to the stability of the tetrameric protein. This observation substantiates the notion of a role for ions as structural “effectors” of ion channels.     

Chemical modification of the bacterial porin OmpF: gain of selectivity by volume reduction

OmpF is an essentially nonselective porin isolated from the outer membrane of Escherichia coli. Here we report on the manipulation of the ion selectivity of OmpF by chemical modification with MTS reagents (MTSET, MTSEA, and MTSES) and the (rather bulky) tripeptide glutathione, all cysteine specific. When recorded in a gradient of 0.1//1 M CaCl2 or 0.1//1 M NaCl, pH 7.4 solutions, measured reversal potentials of the most cation-selective modified mutants were (virtually) identical to the Nernst potential of Ca2+ or Na+. Compared to this full cation selectivity, the anion-selective modified mutants performed somewhat less but nevertheless showed high anion selectivity. We conclude that a low permanent charge in combination with a narrow pore can render the same selectivity as a highly charged but wider pore. These results favor the view that both the electrostatic potential arising form the fixed charge in the pore and the space available at the selectivity filter contribute to the charge selection (i.e., cation versus anion selectivity) of a biological ion channel.