Models of the structure and gating mechanisms of the pore domain of the NaChBac ion channel

The NaChBac prokaryotic sodium channel appears to be a descendent of an evolutionary link between voltage-gated K_V and Ca_V channels. Like K_V channels, four identical six-transmembrane subunits comprise the NaChBac channel, but its selectivity filter possesses a signature sequence of eukaryotic Ca_V channels. We developed structural models of the NaChBac channel in closed and open conformations, using K⁺-channel crystal structures as initial templates. Our models were also consistent with numerous experimental results and modeling criteria. This study concerns the pore domain. The major differences between our models and K⁺ crystal structures involve the latter portion of the selectivity filter and the bend region in S6 of the open conformation. These NaChBac models may serve as a stepping stone between K⁺ channels of known structure and Na_V, Ca_V, and TRP channels of unknown structure.     

Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA,PIIIA, and KIIIA in a NavAb-based model of Nav1.4

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.     

Molecular basis of ion permeability in a voltage‐gated sodium channel

Voltage‐gated sodium channels are essential for electrical signalling across cell membranes. They exhibit strong selectivities for sodium ions over other cations, enabling the finely tuned cascade of events associated with action potentials. This paper describes the ion permeability characteristics and the crystal structure of a prokaryotic sodium channel, showing for the first time the detailed locations of sodium ions in the selectivity filter of a sodium channel. Electrostatic calculations based on the structure are consistent with the relative cation permeability ratios (Na⁺ ≈ Li⁺ ≫ K⁺, Ca²⁺, Mg²⁺) measured for these channels. In an E178D selectivity filter mutant constructed to have altered ion selectivities, the sodium ion binding site nearest the extracellular side is missing. Unlike potassium ions in potassium channels, the sodium ions in these channels appear to be hydrated and are associated with side chains of the selectivity filter residues, rather than polypeptide backbones.     

The Mechanism of Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels Based on Molecular Dynamics Simulation

Voltage-gated sodium (Na_v) channels and their Na⁺/K⁺ selectivity are of great importance in the mammalian neuronal signaling. According to mutational analysis, the Na⁺/K⁺ selectivity in mammalian Na_v channels is mainly determined by the Lys and Asp/Glu residues located at the constriction site within the selectivity filter. Despite successful molecular dynamics simulations conducted on the prokaryotic Na_v channels, the lack of Lys at the constriction site of prokaryotic Na_v channels limits how much can be learned about the Na⁺/K⁺ selectivity in mammalian Na_v channels. In this work, we modeled the mammalian Na_v channel by mutating the key residues at the constriction site in a prokaryotic Na_v channel (Na_vRh) to its mammalian counterpart. By simulating the mutant structure, we found that the Na⁺ preference in mammalian Na_v channels is collaboratively achieved by the deselection from Lys and the selection from Asp/Glu within the constriction site.     

Mutant bacterial sodium channels as models for local anesthetic block of eukaryotic proteins

Voltage gated sodium channels are the target of a range of local anesthetic, anti-epileptic and anti-arrhythmic compounds. But, gaining a molecular level understanding of their mode of action is difficult as we only have atomic resolution structures of bacterial sodium channels not their eukaryotic counterparts. In this study we used molecular dynamics simulations to demonstrate that the binding sites of both the local anesthetic benzocaine and the anti-epileptic phenytoin to the bacterial sodium channel NavAb can be altered significantly by the introduction of point mutations. Free energy techniques were applied to show that increased aromaticity in the pore of the channel, used to emulate the aromatic residues observed in eukaryotic Nav1.2, led to changes in the location of binding and dissociation constants of each drug relative to wild type NavAb. Further, binding locations and dissociation constants obtained for both benzocaine (660 μM) and phenytoin (1 μ M) in the mutant channels were within the range expected from experimental values obtained from drug binding to eukaryotic sodium channels, indicating that these mutant NavAb may be a better model for drug binding to eukaryotic channels than the wild type.     

A naturally occurring amino acid substitution in the voltage-dependent sodium channel selectivity filter affects channel gating

The pore of sodium channels contains a selectivity filter made of 4 amino acids, D/E/K/A. In voltage sensitive sodium channel (Nav) channels from jellyfish to human the fourth amino acid is Ala. This Ala, when mutated to Asp, promotes slow inactivation. In some Nav channels of pufferfishes, the Ala is replaced with Gly. We studied the biophysical properties of an Ala-to-Gly substitution (A1529G) in rat Nav1.4 channel expressed in Xenopus oocytes alone or with a β1 subunit. The Ala-to-Gly substitution does not affect monovalent cation selectivity and positively shifts the voltage-dependent inactivation curve, although co-expression with a β1 subunit eliminates the difference between A1529G and WT. There is almost no difference in channel fast inactivation, but the β1 subunit accelerates WT current inactivation significantly more than it does the A1529G channels. The Ala-to-Gly substitution mainly influences the rate of recovery from slow inactivation. Again, the β1 subunit is less effective on speeding recovery of A1529G than the WT. We searched Nav channels in numerous databases and noted at least four other independent Ala-to-Gly substitutions in Nav channels in teleost fishes. Thus, the Ala-to-Gly substitution occurs more frequently than previously realized, possibly under selection for alterations of channel gating.     

Analysis of the selectivity filter of the voltage-gated sodium channel NavRh

NaChBac is a bacterial voltage-gated sodium (Na_v) channel that shows sequence similarity to voltage-gated calcium channels. To understand the ion-permeation mechanism of Na_v channels, we combined molecular dynamics simulation, structural biology and electrophysiological approaches to investigate the recently determined structure of Na_vRh, a marine bacterial NaChBac ortholog. Two Na⁺ binding sites are identified in the selectivity filter (SF) in our simulations: The extracellular Na⁺ ion first approaches site 1 constituted by the side groups of Ser181 and Glu183, and then spontaneously arrives at the energetically more favorable site 2 formed by the carbonyl oxygens of Leu179 and Thr178. In contrast, Ca²⁺ ions are prone to being trapped by Glu183 at site 1, which then blocks the entrance of both Na⁺ and Ca²⁺ to the vestibule of the SF. In addition, Na⁺ permeates through the selective filter in an asymmetrical manner, a feature that resembles that of the mammalian Na_v orthologs. The study reported here provides insights into the mechanism of ion selectivity on Na⁺ over Ca²⁺ in mammalian Na_v channels.     

Conduction of Na+ and K+ through the NaK channel: molecular and Brownian dynamics studies

Conduction of ions through the NaK channel, with M0 helix removed, was studied using both Brownian dynamics and molecular dynamics. Brownian dynamics simulations predict that the truncated NaK has approximately a third of the conductance of the related KcsA K⁺ channel, is outwardly rectifying, and has a Michaelis-Menten current-concentration relationship. Current magnitude increases when the glutamine residue located near the intracellular gate is replaced with a glutamate residue. The channel is blocked by extracellular Ca²⁺. Molecular dynamics simulations show that, under the influence of a strong applied potential, both Na⁺ and K⁺ move across the selectivity filter, although conduction rates for Na⁺ ions are somewhat lower. The mechanism of conduction of Na⁺ differs significantly from that of K⁺ in that Na⁺ is preferentially coordinated by single planes of pore-lining carbonyl oxygens, instead of two planes as in the usual K⁺ binding sites. The water-containing filter pocket resulting from a single change in the selectivity filter sequence (compared to potassium channels) disrupts several of the planes of carbonyl oxygens, and thus reduces the filter's ability to discriminate against sodium.     

Molecular and functional determinants of local anesthetic inhibition of NaChBac

In our recent publication, we describe the local anesthetic (LA) inhibition of the prokaryotic voltage gated sodium channel NaChBac. Despite the numerous functional and putative structural differences with the mammalian sodium channels, the data show that LA compounds effectively and reversibly inhibit NaChBac channels in a concentration range similar to resting blockade on eukaryotic Navs. In addition to current reduction, LA application accelerated channel inactivation kinetics of NaChBac which could be accounted for in a simple state-model whereby local anesthetics increase the probability of entering the inactivated state. We have further explored what state (or states) local anesthetic blockade of NaChBac could pertain to eukaryotic sodium channels, and what molecular similarities exist between these disparate channel families. Here we show that the rate of recovery from inactivation remains unaffected in the presence of local anesthetics. Further, we show that two sites that support use-dependent inhibition in eukaryotic channels, do not affect block to the same extent when mutated in NaChBac channels. The data indicate that the molecular determinants and the inherent mechanisms for LA block are likely to be divergent between bacterial and eukaryotic Navs, but future experiments will help define possible similarities.     

Molecular and functional determinants of local anesthetic inhibition of NaChBac

In our recent publication, we describe the local anesthetic (LA) inhibition of the prokaryotic voltage gated sodium channel NaChBac. Despite the numerous functional and putative structural differences with the mammalian sodium channels, the data show that LA compounds effectively and reversibly inhibit NaChBac channels in a concentration range similar to resting blockade on eukaryotic Navs. In addition to current reduction, LA application accelerated channel inactivation kinetics of NaChBac which could be accounted for in a simple state-model whereby local anesthetics increase the probability of entering the inactivated state. We have further explored what state (or states) local anesthetic blockade of NaChBac could pertain to eukaryotic sodium channels, and what molecular similarities exist between these disparate channel families. Here we show that the rate of recovery from inactivation remains unaffected in the presence of local anesthetics. Further, we show that two sites that support use-dependent inhibition in eukaryotic channels, do not affect block to the same extent when mutated in NaChBac channels. The data indicate that the molecular determinants and the inherent mechanisms for LA block are likely to be divergent between bacterial and eukaryotic Navs, but future experiments will help define possible similarities.     

Local anesthetic inhibition of a bacterial sodium channel

Recent structural breakthroughs with the voltage-gated sodium channel from Arcobacter butzleri suggest that such bacterial channels may provide a structural platform to advance the understanding of eukaryotic sodium channel gating and pharmacology. We therefore set out to determine whether compounds known to interact with eukaryotic Na(V)s could also inhibit the bacterial channel from Bacillus halodurans and NaChBac and whether they did so through similar mechanisms as in their eukaryotic homologues. The data show that the archetypal local anesthetic (LA) lidocaine inhibits resting NaChBac channels with a dissociation constant (K(d)) of 260 μM, and channels displayed a left-shifted steady-state inactivation gating relationship in the presence of the drug. Extracellular application of QX-314 to expressed NaChBac channels had no effect on sodium current, whereas internal exposure via injection of a bolus of the quaternary derivative rapidly reduced sodium conductance, consistent with a hydrophilic cytoplasmic access pathway to an internal binding site. However, the neutral derivative benzocaine applied externally inhibited NaChBac channels, suggesting that hydrophobic pathways can also provide drug access to inhibit channels. Alternatively, ranolazine, a putative preopen state blocker of eukaryotic Na(V)s, displayed a K(d) of 60 μM and left-shifted the NaChBac activation-voltage relationship. In each case, block enhanced entry into the inactivated state of the channel, an effect that is well described by a simple kinetic scheme. The data suggest that although significant differences exist, LA block of eukaryotic Na(V)s also occurs in bacterial sodium channels and that NaChBac shares pharmacological homology to the resting state of vertebrate Na(V) homologues.     

Differential Lipid Dependence of the Function of Bacterial Sodium Channels

The lipid bilayer is important for maintaining the integrity of cellular compartments and plays a vital role in providing the hydrophobic and charged interactions necessary for membrane protein structure, conformational flexibility and function. To directly assess the lipid dependence of activity for voltage-gated sodium channels, we compared the activity of three bacterial sodium channel homologues (NaChBac, NavMs, and NavSp) by cumulative ²²Na⁺ uptake into proteoliposomes containing a 3∶1 ratio of 1-palmitoyl 2-oleoyl phosphatidylethanolamine and different “guest” glycerophospholipids. We observed a unique lipid profile for each channel tested. NavMs and NavSp showed strong preference for different negatively-charged lipids (phosphatidylinositol and phosphatidylglycerol, respectively), whilst NaChBac exhibited a more modest variation with lipid type. To investigate the molecular bases of these differences we used synchrotron radiation circular dichroism spectroscopy to compare structures in liposomes of different composition, and molecular modeling and electrostatics calculations to rationalize the functional differences seen. We then examined pore-only constructs (with voltage sensor subdomains removed) and found that in these channels the lipid specificity was drastically reduced, suggesting that the specific lipid influences on voltage-gated sodium channels arise primarily from their abilities to interact with the voltage-sensing subdomains.     

Exploring the Structure of the Voltage-gated Na+ Channel by an Engineered Drug Access Pathway to the Receptor Site for Local Anesthetics

Despite the availability of several crystal structures of bacterial voltage-gated Na⁺ channels, the structure of eukaryotic Na⁺ channels is still undefined. We used predictions from available homology models and crystal structures to modulate an external access pathway for the membrane-impermeant local anesthetic derivative QX-222 into the internal vestibule of the mammalian rNa_V1.4 channel. Potassium channel-based homology models predict amino acid Ile-1575 in domain IV segment 6 to be in close proximity to Lys-1237 of the domain III pore-loop selectivity filter. The mutation K1237E has been shown previously to increase the diameter of the selectivity filter. We found that an access pathway for external QX-222 created by mutations of Ile-1575 was abolished by the additional mutation K1237E, supporting the notion of a close spatial relationship between sites 1237 and 1575. Crystal structures of bacterial voltage-gated Na⁺ channels predict that the side chain of rNa_V1.4 Trp-1531 of the domain IV pore-loop projects into the space between domain IV segment 6 and domain III pore-loop and, therefore, should obstruct the putative external access pathway. Indeed, mutations W1531A and W1531G allowed for exceptionally rapid access of QX-222. In addition, W1531G created a second non-selective ion-conducting pore, bypassing the outer vestibule but probably merging into the internal vestibule, allowing for control by the activation gate. These data suggest a strong structural similarity between bacterial and eukaryotic voltage-gated Na⁺ channels.     

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

The Fibroblast Growth Factor 14·Voltage-gated Sodium Channel Complex Is a New Target of Glycogen Synthase Kinase 3 (GSK3)

The FGF14 protein controls biophysical properties and subcellular distribution of neuronal voltage-gated Na⁺ (Nav) channels through direct binding to the channel C terminus. To gain insights into the dynamic regulation of this protein/protein interaction complex, we employed the split luciferase complementation assay to screen a small molecule library of kinase inhibitors against the FGF14·Nav1.6 channel complex and identified inhibitors of GSK3 as hits. Through a combination of a luminescence-based counter-screening, co-immunoprecipitation, patch clamp electrophysiology, and quantitative confocal immunofluorescence, we demonstrate that inhibition of GSK3 reduces the assembly of the FGF14·Nav channel complex, modifies FGF14-dependent regulation of Na⁺ currents, and induces dissociation and subcellular redistribution of the native FGF14·Nav channel complex in hippocampal neurons. These results further emphasize the role of FGF14 as a critical component of the Nav channel macromolecular complex, providing evidence for a novel GSK3-dependent signaling pathway that might control excitability through specific protein/protein interactions.     

Non-Equilibrium Dynamics Contribute to Ion Selectivity in the KcsA Channel

The ability of biological ion channels to conduct selected ions across cell membranes is critical for the survival of both animal and bacterial cells. Numerous investigations of ion selectivity have been conducted over more than 50 years, yet the mechanisms whereby the channels select certain ions and reject others are not well understood. Here we report a new application of Jarzynski’s Equality to investigate the mechanism of ion selectivity using non-equilibrium molecular dynamics simulations of Na⁺ and K⁺ ions moving through the KcsA channel. The simulations show that the selectivity filter of KcsA adapts and responds to the presence of the ions with structural rearrangements that are different for Na⁺ and K⁺. These structural rearrangements facilitate entry of K⁺ ions into the selectivity filter and permeation through the channel, and rejection of Na⁺ ions. A mechanistic model of ion selectivity by this channel based on the results of the simulations relates the structural rearrangement of the selectivity filter to the differential dehydration of ions and multiple-ion occupancy and describes a mechanism to efficiently select and conduct K⁺. Estimates of the K⁺/Na⁺ selectivity ratio and steady state ion conductance for KcsA from the simulations are in good quantitative agreement with experimental measurements. This model also accurately describes experimental observations of channel block by cytoplasmic Na⁺ ions, the “punch through” relief of channel block by cytoplasmic positive voltages, and is consistent with the knock-on mechanism of ion permeation.     

Ion Selectivity of the KcsA Channel: A Perspective from Multi-Ion Free Energy Landscapes

Potassium (K⁺) channels are specialized membrane proteins that are able to facilitate and regulate the conduction of K⁺ through cell membranes. Comprising five specific cation binding sites (S₀-S₄) formed by the backbone carbonyl groups of conserved residues common to all K⁺ channels, the narrow selectivity filter allows fast conduction of K⁺ while being highly selective for K⁺ over Na⁺. To extend our knowledge of the microscopic mechanism underlying selectivity in K⁺ channels, we characterize the free energy landscapes governing the entry and translocation of a Na⁺ or a K⁺ from the extracellular side into the selectivity filter of KcsA. The entry process of an extracellular ion is examined in the presence of two additional K⁺ in the pore, and the three-ion potential of mean force is computed using extensive all-atom umbrella sampling molecular dynamics simulations. A comparison of the potentials of mean force yields a number of important results. First, the free energy minima corresponding to configurations with extracellular K⁺ or Na⁺ in binding site S₀ or S₁ are similar in depth, suggesting that the thermodynamic selectivity governed by the free energy minima for those two binding sites is insignificant. Second, the free energy barriers between stable multi-ion configurations are generally higher for Na⁺ than for K⁺, implying that the kinetics of ion conduction is slower when a Na⁺ enters the pore. Third, the region corresponding to binding site S₂ near the center of the narrow pore emerges as the most selective for K⁺ over Na⁺. In particular, while there is a stable minimum for K⁺ in site S₂, Na⁺ faces a steep free energy increase with no local free energy well in this region. Lastly, analysis shows that selectivity is not correlated with the overall coordination number of the ion entering the pore, but is predominantly affected by changes in the type of coordinating ligands (carbonyls versus water molecules). These results further highlight the importance of the central region near binding site S₂ in the selectivity filter of K⁺ channels.     

K+ Block Is the Mechanism of Functional Asymmetry in Bacterial Nav Channels

Crystal structures of several bacterial Na_v channels have been recently published and molecular dynamics simulations of ion permeation through these channels are consistent with many electrophysiological properties of eukaryotic channels. Bacterial Na_v channels have been characterized as functionally asymmetric, and the mechanism of this asymmetry has not been clearly understood. To address this question, we combined non-equilibrium simulation data with two-dimensional equilibrium unperturbed landscapes generated by umbrella sampling and Weighted Histogram Analysis Methods for multiple ions traversing the selectivity filter of bacterial Na_vAb channel. This approach provided new insight into the mechanism of selective ion permeation in bacterial Na_v channels. The non-equilibrium simulations indicate that two or three extracellular K⁺ ions can block the entrance to the selectivity filter of Na_vAb in the presence of applied forces in the inward direction, but not in the outward direction. The block state occurs in an unstable local minimum of the equilibrium unperturbed free-energy landscape of two K⁺ ions that can be ‘locked’ in place by modest applied forces. In contrast to K⁺, three Na⁺ ions move favorably through the selectivity filter together as a unit in a loose “knock-on” mechanism of permeation in both inward and outward directions, and there is no similar local minimum in the two-dimensional free-energy landscape of two Na⁺ ions for a block state. The useful work predicted by the non-equilibrium simulations that is required to break the K⁺ block is equivalent to large applied potentials experimentally measured for two bacterial Na_v channels to induce inward currents of K⁺ ions. These results illustrate how inclusion of non-equilibrium factors in the simulations can provide detailed information about mechanisms of ion selectivity that is missing from mechanisms derived from either crystal structures or equilibrium unperturbed free-energy landscapes.     

Fast and slow gating are inherent properties of the pore module of the K+ channel Kcv

Kcv from the chlorella virus PBCV-1 is a viral protein that forms a tetrameric, functional K⁺ channel in heterologous systems. Kcv can serve as a model system to study and manipulate basic properties of the K⁺ channel pore because its minimalistic structure (94 amino acids) produces basic features of ion channels, such as selectivity, gating, and sensitivity to blockers. We present a characterization of Kcv properties at the single-channel level. In symmetric 100 mM K⁺, single-channel conductance is 114 ± 11 pS. Two different voltage-dependent mechanisms are responsible for the gating of Kcv. “Fast” gating, analyzed by β distributions, is responsible for the negative slope conductance in the single-channel current–voltage curve at extreme potentials, like in MaxiK potassium channels, and can be explained by depletion-aggravated instability of the filter region. The presence of a “slow” gating is revealed by the very low (in the order of 1–4%) mean open probability that is voltage dependent and underlies the time-dependent component of the macroscopic current.     

Altered Ionic Selectivity of the Sodium Channel Revealed by Cysteine Mutations within the Pore

To explore the role of pore-lining amino acids in Na⁺ channel ion-selectivity, pore residues were  replaced serially with cysteine in cloned rat skeletal muscle Na⁺ channels. Ionic selectivity was determined by measuring permeability and ionic current ratios of whole-cell currents in Xenopus oocytes. The rSkM1 channels displayed an ionic selectivity sequence Na⁺>Li⁺>NH₄ ⁺>>K⁺>>Cs⁺ and were impermeable to divalent cations.  Replacement of residues in domain IV showed significantly enhanced current and permeability ratios of NH₄ ⁺ and K⁺, and negative shifts in the reversal potentials recorded in the presence of external Na⁺ solutions when compared to cysteine mutants in domains I, II, and III (except K1237C). Mutants in domain IV showed altered selectivity sequences: W1531C (NH₄ ⁺>K⁺>Na⁺≥Li⁺≈Cs⁺), D1532C, and G1533C (Na⁺>Li⁺≥NH₄ ⁺>K⁺>Cs⁺). Conservative replacement of the aromatic residue in domain IV (W1531) with phenylalanine or tyrosine retained Na⁺ selectivity of the channel while the alanine mutant (W1531A) reduced ion selectivity. A single mutation within the third pore forming region (K1237C) dramatically altered the selectivity sequence of the rSkM1 channel (NH₄ ⁺>K⁺>Na⁺≥Li⁺≈Cs⁺) and was permeable to divalent cations having the selectivity sequence Ca²⁺≥Sr²⁺>Mg²⁺>Ba²⁺. Sulfhydryl modification of K1237C, W1531C or D1532C with methanethiosulfonate derivatives that introduce a positively charged ammonium group, large trimethylammonium moiety, or a negatively charged sulfonate group within the pore was ineffective in restoring Na⁺ selectivity to these channels. Selectivity of D1532C mutants could be largely restored by increasing extracellular pH suggesting altering the ionized state at this position influences selectivity. These data suggest that K1237 in domain III and W1531, D1532, and G1533 in domain IV play a critical role in determining the ionic selectivity of the Na⁺ channel.