Different inward and outward conduction mechanisms in NaVMs suggested by molecular dynamics simulations

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.

A cooperative knock-on mechanism underpins Ca2+-selective cation permeation in TRPV channels

The selective exchange of ions across cellular membranes is a vital biological process. Ca2+-mediated signaling is implicated in a broad array of physiological processes in cells, while elevated intracellular concentrations of Ca2+ are cytotoxic. Due to the significance of this cation, strict Ca2+ concentration gradients are maintained across the plasma and organelle membranes. Therefore, Ca2+ signaling relies on permeation through selective ion channels that control the flux of Ca2+ ions. A key family of Ca2+-permeable membrane channels is the polymodal signal-detecting transient receptor potential (TRP) ion channels. TRP channels are activated by a wide variety of cues including temperature, small molecules, transmembrane voltage, and mechanical stimuli. While most members of this family permeate a broad range of cations non-selectively, TRPV5 and TRPV6 are unique due to their strong Ca2+ selectivity. Here, we address the question of how some members of the TRPV subfamily show a high degree of Ca2+ selectivity while others conduct a wider spectrum of cations. We present results from all-atom molecular dynamics simulations of ion permeation through two Ca2+-selective and two non-selective TRPV channels. Using a new method to quantify permeation cooperativity based on mutual information, we show that Ca2+-selective TRPV channel permeation occurs by a three-binding site knock-on mechanism, whereas a two-binding site knock-on mechanism is observed in non-selective TRPV channels. Each of the ion binding sites involved displayed greater affinity for Ca2+ over Na+. As such, our results suggest that coupling to an extra binding site in the Ca2+-selective TRPV channels underpins their increased selectivity for Ca2+ over Na+ ions. Furthermore, analysis of all available TRPV channel structures shows that the selectivity filter entrance region is wider for the non-selective TRPV channels, slightly destabilizing ion binding at this site, which is likely to underlie mechanistic decoupling.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Determinants of conductance of a bacterial voltage-gated sodium channel

Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel Na_vMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na_v channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na⁺ in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.

The regulatory function of mixed lineage kinase 3 in tumor and host immunity

Protein kinases are the second most sought-after G-protein coupled receptors as drug targets because of their overexpression, mutations, and dysregulated catalytic activities in various pathological conditions. Till 2019, 48 protein kinase inhibitors have received FDA approval for the treatment of multiple illnesses, of which the majority of them are indicated for different malignancies. One of the attractive sub-group of protein kinases that has attracted attention for drug development is the family members of MAPKs that are recognized to play significant roles in different cancers. Several inhibitors have been developed against various MAPK members; however, none of them as monotherapy has shown sustainable efficacy. One of the MAPK members, called Mixed Lineage Kinase 3 (MLK3), has attracted considerable attention due to its role in inflammation and neurodegenerative diseases; however, its role in cancer is an emerging area that needs more investigation. Recent advances have shown that MLK3 plays a role in cancer cell survival, migration, drug resistance, cell death, and tumor immunity. This review describes how MLK3 regulates different MAPK pathways, cancer cell growth and survival, apoptosis, and host's immunity. We also discuss how MLK3 inhibitors can potentially be used along with immunotherapy for different malignancies.

IR Spectroscopy Can Reveal the Mechanism of K+ Transport in Ion Channels

Ion channels like KcsA enable ions to move across cell membranes at near diffusion-limited rates and with very high selectivity. Various mechanisms have been proposed to explain this phenomenon. Broadly, there is disagreement among the proposed mechanisms about whether ions occupy adjacent sites in the channel during the transport process. Here, using a mixed quantum-classical approach to calculate theoretical infrared spectra, we propose a set of infrared spectroscopy experiments that can discriminate between mechanisms with and without adjacent ions. These experiments differ from previous ones in that they independently probe specific ion binding sites within the selectivity filter. When ions occupy adjacent sites in the selectivity filter, the predicted spectra are significantly redshifted relative to when ions do not occupy adjacent sites. Comparisons between theoretical and experimental peak frequencies will therefore discriminate the mechanisms.

Bases of Bacterial Sodium Channel Selectivity Among Organic Cations

Hille's (1971) seminal study of organic cation selectivity of eukaryotic voltage-gated sodium channels showed a sharp size cut-off for ion permeation, such that no ion possessing a methyl group was permeant. Using the prokaryotic channel, NaChBac, we found some similarity and two peculiar differences in the selectivity profiles for small polyatomic cations. First, we identified a diverse group of minimally permeant cations for wildtype NaChBac, ranging in sizes from ammonium to guanidinium and tetramethylammonium; and second, for both ammonium and hydrazinium, the charge-conserving selectivity filter mutation (E191D) yielded substantial increases in relative permeability (PX/PNa). The relative permeabilities varied inversely with relative Kd calculated from 1D Potential of Mean Force profiles (PMFs) for the single cations traversing the channel. Several of the cations bound more strongly than Na+, and hence appear to act as blockers, as well as charge carriers. Consistent with experimental observations, the E191D mutation had little impact on Na+ binding to the selectivity filter, but disrupted the binding of ammonium and hydrazinium, consequently facilitating ion permeation across the NaChBac-like filter. We concluded that for prokaryotic sodium channels, a fine balance among filter size, binding affinity, occupancy, and flexibility seems to contribute to observed functional differences.

Molecular Dynamics of Ion Conduction through the Selectivity Filter of the NaVAb Sodium Channel

The determination of the atomic structures of voltage-gated bacterial sodium channels using X-ray crystallography has provided a first view of this family of membrane proteins. Molecular dynamics simulations offer one approach to clarify the underlying mechanism of permeation and selectivity in these channels. However, it appears that the intracellular gate of the pore domain is either closed or only open partially in the available X-ray structures. The lack of structure with a fully open intracellular gate poses a special challenge to computational studies aimed at simulating ion conduction. To circumvent this problem, we simulated a model of the Na_VAb channel in which the transmembrane S5 and S6 helices of the pore domain have been truncated to provide direct open access to the intracellular entryway to the pore. Molecular dynamics simulations were carried out over a range of membrane potential and ion concentration of sodium and potassium. The simulations show that the Na_VAb selectivity filter is essentially a cationic pore supporting the conduction of ions at a rate comparable to aqueous diffusion with no significant selectivity for sodium. Conductance and selectivity vary as a function of ion concentration for both cations. Permeation occurs primarily via a knock-on mechanism for both sodium and potassium, although the ion ordering in single file along the pore is not strictly maintained. The character of the outward current appears quite different from the inward current, with a buildup on ions in the selectivity filter prior to escape toward the extracellular side, indicating the presence of a rectification effect that is overcome by nonphysiological applied voltages.

Role of the Interaction Motif in Maintaining the Open Gate of an Open Sodium Channel

Voltage-gated sodium channels undergo transitions between open, closed, and inactivated states, enabling regulation of the translocation of sodium ions across membranes. A recently published crystal structure of the full-length prokaryotic NavMs crystal structure in the activated open conformation has revealed the presence of a novel motif consisting of an extensive network of salt bridges involving residues in the voltage sensor, S4-S5 linker, pore, and C-terminal domains. This motif has been proposed to be responsible for maintaining an open conformation that enables ion translocation through the channel. In this study, we have used long-time molecular dynamics calculations without artificial restraints to demonstrate that the interaction network of full-length NavMs indeed prevents a rapid collapse and closure of the gate, in marked difference to earlier studies of the pore-only construct in which the gate had to be restrained to remain open. Interestingly, a frequently discussed "hydrophobic gating" mechanism at nanoscopic level is also observed in our simulations, in which the discontinuous water wire close to the gate region leads to an energetic barrier for ion conduction. In addition, we demonstrate the effects of in silico mutations of several of the key residues in the motif on the open channel's stability and functioning, correlating them with existing functional studies on this channel and homologous disease-associated mutations in human sodium channels; we also examine the effects of truncating/removing the voltage sensor and C-terminal domains in maintaining an open gate.

Selective ion permeation involves complexation with carboxylates and lysine in a model human sodium channel

Bacterial and human voltage-gated sodium channels (Navs) exhibit similar cation selectivity, despite their distinct EEEE and DEKA selectivity filter signature sequences. Recent high-resolution structures for bacterial Navs have allowed us to learn about ion conduction mechanisms in these simpler homo-tetrameric channels, but our understanding of the function of their mammalian counterparts remains limited. To probe these conduction mechanisms, a model of the human Nav1.2 channel has been constructed by grafting residues of its selectivity filter and external vestibular region onto the bacterial NavRh channel with atomic-resolution structure. Multi-μs fully atomistic simulations capture long time-scale ion and protein movements associated with the permeation of Na+ and K+ ions, and their differences. We observe a Na+ ion knock-on conduction mechanism facilitated by low energy multi-carboxylate/multi-Na+ complexes, akin to the bacterial channels. These complexes involve both the DEKA and vestibular EEDD rings, acting to draw multiple Na+ into the selectivity filter and promote permeation. When the DEKA ring lysine is protonated, we observe that its ammonium group is actively participating in Na+ permeation, presuming the role of another ion. It participates in the formation of a stable complex involving carboxylates that collectively bind both Na+ and the Lys ammonium group in a high-field strength site, permitting pass-by translocation of Na+. In contrast, multiple K+ ion complexes with the DEKA and EEDD rings are disfavored by up to 8.3 kcal/mol, with the K+-lysine-carboxylate complex non-existent. As a result, lysine acts as an electrostatic plug that partially blocks the flow of K+ ions, which must instead wait for isomerization of lysine downward to clear the path for K+ passage. These distinct mechanisms give us insight into the nature of ion conduction and selectivity in human Nav channels, while uncovering high field strength carboxylate binding complexes that define the more general phenomenon of Na+-selective ion transport in nature.

Comparison of permeation mechanisms in sodium-selective ion channels

Voltage-gated sodium channels are the molecular components of electrical signaling in the body, yet the molecular origins of Na⁺-selective transport remain obscured by diverse protein chemistries within this family of ion channels. In particular, bacterial and mammalian sodium channels are known to exhibit similar relative ion permeabilities for Na⁺ over K⁺ ions, despite their distinct signature EEEE and DEKA sequences. Atomic-level molecular dynamics simulations using high-resolution bacterial channel structures and mammalian channel models have begun to describe how these sequences lead to analogous high field strength ion binding sites that drive Na⁺ conduction. Similar complexes have also been identified in unrelated acid sensing ion channels involving glutamate and aspartate side chains that control their selectivity. These studies suggest the possibility of a common origin for Na⁺ selective binding and transport.

Faster Simulations with a 5 fs Time Step for Lipids in the CHARMM Force Field

The performance of all-atom molecular dynamics simulations is limited by an integration time step of 2 fs, which is needed to resolve the fastest degrees of freedom in the system, namely, the vibration of bonds and angles involving hydrogen atoms. The virtual interaction sites (VIS) method replaces hydrogen atoms by massless virtual interaction sites to eliminate these degrees of freedom while keeping intact nonbonded interactions and the explicit treatment of hydrogen atoms. We have modified the existing VIS algorithm for most lipids in the popular CHARMM36 force field by increasing the hydrogen atom masses at regular intervals in the lipid acyl chains and obtained lipid properties and pore formation free energies in very good agreement with those calculated in simulations without VIS. Our modified VIS scheme enables a 5 fs time step resulting in a significant performance gain for all-atom simulations of membranes. The method has the potential to make longer time and length scales accessible in all-atom simulations of membrane-protein complexes.

Ion-triggered selectivity in bacterial sodium channels

Since the availability of the first crystal structure of a bacterial Na⁺ channel in 2011, understanding selectivity across this family of membrane proteins has been the subject of intense research efforts. Initially, free energy calculations based on molecular dynamics simulations revealed that although sodium ions can easily permeate the channel with their first hydration shell almost intact, the selectivity filter is too narrow for efficient conduction of hydrated potassium ions. This steric view of selectivity was subsequently questioned by microsecond atomic trajectories, which proved that the selectivity filter appears to the permeating ions as a highly degenerate, liquid-like environment. Although this liquid-like environment looks optimal for rapid conduction of Na⁺, it seems incompatible with efficient discrimination between similar ion species, such as Na⁺ and K⁺, through steric effects. Here extensive molecular dynamics simulations, combined with Markov state model analyses, reveal that at positive membrane potentials, potassium ions trigger a conformational change of the selectivity toward a nonconductive metastable state. It is this transition of the selectivity filter, and not steric effects, that prevents the outward flux of K⁺ at positive membrane potentials. This description of selectivity, triggered by the nature of the permeating ions, might have implications on the current understanding of how ion channels, and in particular bacterial Na⁺ channels, operate at the atomic scale.

The hitchhiker's guide to the voltage-gated sodium channel galaxy

Eukaryotic voltage-gated sodium (Na_v) channels contribute to the rising phase of action potentials and served as an early muse for biophysicists laying the foundation for our current understanding of electrical signaling. Given their central role in electrical excitability, it is not surprising that (a) inherited mutations in genes encoding for Na_v channels and their accessory subunits have been linked to excitability disorders in brain, muscle, and heart; and (b) Na_v channels are targeted by various drugs and naturally occurring toxins. Although the overall architecture and behavior of these channels are likely to be similar to the more well-studied voltage-gated potassium channels, eukaryotic Na_v channels lack structural and functional symmetry, a notable difference that has implications for gating and selectivity. Activation of voltage-sensing modules of the first three domains in Na_v channels is sufficient to open the channel pore, whereas movement of the domain IV voltage sensor is correlated with inactivation. Also, structure–function studies of eukaryotic Na_v channels show that a set of amino acids in the selectivity filter, referred to as DEKA locus, is essential for Na⁺ selectivity. Structures of prokaryotic Na_v channels have also shed new light on mechanisms of drug block. These structures exhibit lateral fenestrations that are large enough to allow drugs or lipophilic molecules to gain access into the inner vestibule, suggesting that this might be the passage for drug entry into a closed channel. In this Review, we will synthesize our current understanding of Na_v channel gating mechanisms, ion selectivity and permeation, and modulation by therapeutics and toxins in light of the new structures of the prokaryotic Na_v channels that, for the time being, serve as structural models of their eukaryotic counterparts.

Bias-Exchange Metadynamics Simulations: An Efficient Strategy for the Analysis of Conduction and Selectivity in Ion Channels

Conduction through ion channels possesses two interesting features: (i) different ionic species are selected with high-selectivity and (ii) ions travel across the channel with rates approaching free-diffusion. Molecular dynamics simulations have the potential to reveal how these processes take place at the atomic level. However, analysis of conduction and selectivity at atomistic detail is still hampered by the short time scales accessible by computer simulations. Several algorithms have been developed to "accelerate" sampling along the slow degrees of freedom of the process under study and thus to probe longer time scales. In these algorithms, the slow degrees of freedom need to be defined in advance, which is a well-known shortcoming. In the particular case of ion conduction, preliminary assumptions about the number and type of ions participating in the permeation process need to be made. In this study, a novel approach for the analysis of conduction and selectivity based on bias-exchange metadynamics simulations was tested. This approach was compared with umbrella sampling simulations, using a model of a Na(+)-selective channel. Analogous conclusions resulted from both techniques, but the computational cost of bias-exchange simulations was lower. In addition, with bias-exchange metadynamics it was possible to calculate free energy profiles in the presence of a variable number and type of permeating ions. This approach might facilitate the definition of the set of collective variables required to analyze conduction and selectivity in ion channels.

Biophysical Adaptations of Prokaryotic Voltage-Gated Sodium Channels

This chapter describes the adaptive features found in voltage-gated sodium channels (NaVs) of prokaryotes and eukaryotes. These two families are distinct, having diverged early in evolutionary history but maintain a surprising degree of convergence in function. While prokaryotic NaVs are required for growth and motility, eukaryotic NaVs selectively conduct fast electrical currents for short- and long-range signaling across cell membranes in mammalian organs. Current interest in prokaryotic NaVs is stoked by their resolved high-resolution structures and functional features which are reminiscent of eukaryotic NaVs. In this chapter, comparisons between eukaryotic and prokaryotic NaVs are made to highlight the shared and unique aspects of ion selectivity, voltage sensitivity, and pharmacology. Examples of prokaryotic and eukaryotic NaV convergent evolution will be discussed within the context of their structural features.

Computational Study of Binding of μ-Conotoxin GIIIA to Bacterial Sodium Channels NaVAb and NaVRh

Structures of several voltage-gated sodium (NaV) channels from bacteria have been determined recently, but the same feat might not be achieved for the mammalian counterparts in the near future. Thus, at present, computational studies of the mammalian NaV channels have to be performed using homology models based on the bacterial crystal structures. A successful homology model for the mammalian NaV1.4 channel was recently constructed using the extensive mutation data for binding of μ-conotoxin GIIIA to NaV1.4, which was further validated through studies of binding of other μ-conotoxins and ion permeation. Understanding the similarities and differences between the bacterial and mammalian NaV channels is an important issue, and the NaV1.4-GIIIA system provides a good opportunity for such a comparison. To this end, we study the binding of GIIIA to the bacterial channels NaVAb and NaVRh. The complex structures are obtained using docking and molecular dynamics simulations, and the dissociation of GIIIA is studied through umbrella sampling simulations. The results are compared to those obtained from the NaV1.4-GIIIA system, and the differences in the binding modes arising from the changes in the selectivity filters are highlighted.

Voltage-Gated Sodium Channels: Mechanistic Insights From Atomistic Molecular Dynamics Simulations

The permeation of ions and other molecules across biological membranes is an inherent requirement of all cellular organisms. Ion channels, in particular, are responsible for the conduction of charged species, hence modulating the propagation of electrical signals. Despite the universal physiological implications of this property, the molecular functioning of ion channels remains ambiguous. The combination of atomistic structural data with computational methodologies, such as molecular dynamics (MD) simulations, is now considered routine to investigate structure-function relationships in biological systems. A fuller understanding of conduction, selectivity, and gating, therefore, is steadily emerging due to the applicability of these techniques to ion channels. However, because their structure is known at atomic resolution, studies have consistently been biased toward K(+) channels, thus the molecular determinants of ionic selectivity, activation, and drug blockage in Na(+) channels are often overlooked. The recent increase of available crystallographic data has eminently encouraged the investigation of voltage-gated sodium (NaV) channels via computational methods. Here, we present an overview of simulation studies that have contributed to our understanding of key principles that underlie ionic conduction and selectivity in Na(+) channels, in comparison to the K(+) channel analogs.

K+ Block Is the Mechanism of Functional Asymmetry in Bacterial Na(v) 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.

In this paper we show that bacterial voltage-gated sodium channels (Na_v) conduct ions differently in the inward direction and in the outward direction because extracellular K⁺ blocks these channels. The mechanism of block by K⁺ is revealed in the simulations only when effects of applied forces in the inward direction are considered appropriately. The block of bacterial Na_v channels by K⁺ explains the experimental findings that inward K⁺ current in bacterial Na_v channels is very small under conditions where most mammalian Na_v channels conduct significant K⁺ current. The block by K⁺ ions can occur in the bacterial Na_v channels but is very unlikely to occur in mammalian channels due to differences in the amino acid sequences of the selectivity filters of the different channels.