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.     

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.     

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.     

Hydration valve controlled non-selective conduction of Na and K in the NaK channel

The Na⁺ and K⁺ channels are essential to neural signaling, but our current knowledge at the atomic level is mainly limited to the conducting mechanism of K⁺. Unlike a K⁺ channel having four equivalent K⁺-binding sites in its selectivity filter, a NaK channel has a vestibule in the middle part of its selectivity filter, and can conduct both Na⁺ and K⁺ ions. However, the underlying mechanism for non-selective ion conduction in NaK remains elusive. Here we find four small grottos connecting with the vestibule of the NaK selectivity filter, which form a vestibule-grotto complex perpendicular to the filter pore with a few water molecules within it. It is shown that two or more of the water molecules coming to the vestibule to coordinate the cation are necessary for conducting both Na⁺ and K⁺ ions, while only one water molecule in the vestibule will obstruct ion permeation. Thus, the complex with the aid of interior water movement forms a dynamic hydration valve which is flexible in conveying different cations through the vestibule. Similar exquisite hydration valve mechanisms are expected to be utilized by other non-selective cation channels, and the results should shed new light on the importance of water in neural signaling.     

Mechanism of tetrodotoxin block and resistance in sodium channels

Tetrodotoxin (TTX) has been used for many decades to characterize the structure and function of biological ion channels. Yet, the precise mechanism by which TTX blocks voltage-gated sodium (Na_V) channels is not fully understood. Here molecular dynamics simulations are used to elucidate how TTX blocks mammalian voltage-gated sodium (Nav) channels and why it fails to be effective for the bacterial sodium channel, Na_VAb. We find that, in Na_VAb, a sodium ion competes with TTX for the binding site at the extracellular end of the filter, thus reducing the blocking efficacy of TTX. Using a model of the skeletal muscle channel, Na_V1.4, we show that the conduction properties of the channel observed experimentally are faithfully reproduced. We find that TTX occludes the entrance of Na_V1.4 by forming a network of hydrogen-bonds at the outer lumen of the selectivity filter. The guanidine group of TTX adopts a lateral orientation, rather than pointing into the filter as proposed previously. The acidic residues just above the selectivity filter are important in stabilizing the hydrogen-bond network between TTX and Na_V1.4. The effect of two single mutations of a critical tyrosine residue in the filter of Na_V1.4 on TTX binding observed experimentally is reproduced using computational mutagenesis.     

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.     

K binding in the G-loop and water cavity facilitates Ba movement in the Kir2.1 channel

K⁺ are selectively coordinated in the selectivity filter and concerted K⁺ and water movements in this region ensure high conduction rates in K⁺ channels. In channels with long pores many K⁺ binding sites are located intracellular to the selectivity filter (inner vestibule), but their contribution to permeation has not been well studied. We investigated this phenomenon by slowing the ion permeation process via blocking inwardly rectifying Kir2.1 channels with Ba²⁺ in the selectivity filter and observing the effect of K⁺ in the inner vestibule on Ba²⁺ exit. The dose-response effect of the intracellular K⁺ concentration ([K⁺]_i) on Ba²⁺ exit was recorded with and without intracellular polyamines, which compete with K⁺ for binding sites. Ba²⁺ exit was facilitated by the cooperative binding of at least three K⁺. Site-directed mutagenesis studies suggest that K⁺ interacting with Ba²⁺ bound in the selectivity filter were located in the region between selectivity filter and cytoplasmic pore, i.e. the water cavity and G-loop. One of the K⁺ binding sites was located at residue D172 and another was possibly at M301. This study provides functional evidence for the three K⁺ binding sites in the inner vestibule previously identified by crystal structure study.     

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.     

Structure of a Prokaryotic Sodium Channel Pore Reveals Essential Gating Elements and an Outer Ion Binding Site Common to Eukaryotic Channels

Voltage-gated sodium channels (Na_Vs) are central elements of cellular excitation. Notwithstanding advances from recent bacterial Na_V (BacNa_V) structures, key questions about gating and ion selectivity remain. Here, we present a closed conformation of Na_VAe1p, a pore-only BacNa_V derived from Na_VAe1, a BacNa_V from the arsenite oxidizer Alkalilimnicola ehrlichei found in Mono Lake, California, that provides insight into both fundamental properties. The structure reveals a pore domain in which the pore-lining S6 helix connects to a helical cytoplasmic tail. Electrophysiological studies of full-length BacNa_Vs show that two elements defined by the Na_VAe1p structure, an S6 activation gate position and the cytoplasmic tail “neck”, are central to BacNa_V gating. The structure also reveals the selectivity filter ion entry site, termed the “outer ion” site. Comparison with mammalian voltage-gated calcium channel (Ca_V) selectivity filters, together with functional studies, shows that this site forms a previously unknown determinant of Ca_V high-affinity calcium binding. Our findings underscore commonalities between BacNa_Vs and eukaryotic voltage-gated channels and provide a framework for understanding gating and ion permeation in this superfamily.     

Molecular Dynamics Study of Binding of μ-Conotoxin GIIIA to the Voltage-Gated Sodium Channel Nav1.4

Homology models of mammalian voltage-gated sodium (Na_V) channels based on the crystal structures of the bacterial counterparts are needed to interpret the functional data on sodium channels and understand how they operate. Such models would also be invaluable in structure-based design of therapeutics for diseases involving sodium channels such as chronic pain and heart diseases. Here we construct a homology model for the pore domain of the Na_V1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to Na_V1.4 to validate the model. The initial poses for the Na_V1.4–GIIIA complex are obtained using the HADDOCK protein docking program, which are then refined in molecular dynamics simulations. The binding mode for the final complex is shown to be in broad agreement with the available mutagenesis data. The standard binding free energy, determined from the potential of mean force calculations, is also in good agreement with the experimental value. Because the pore domains of Na_V1 channels are highly homologous, the model constructed for Na_V1.4 will provide an excellent template for other Na_V1 channels.     

Use-dependent block of the voltage-gated Na+ channel by tetrodotoxin and saxitoxin: Effect of pore mutations that change ionic selectivity

Voltage-gated Na⁺ channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca²⁺ permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca²⁺ or Na⁺ ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca²⁺ permeability, suggesting that ion–toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.     

Sensitivity of cloned muscle,heart and neuronal voltage-gated sodium channels to block by polyamines: A possible basis for modulation of excitability in vivo

Spermidine and spermine, are endogenous polyamines (PAs) that regulate cell growth and modulate the activity of numerous ion channel proteins. In particular, intracellular PAs are potent blockers of many different cation channels and are responsible for strong suppression of outward K⁺ current, a phenomenon known as inward rectification characteristic of a major class of K_IR K⁺ channels. We previously described block of heterologously expressed voltage-gated Na⁺ channels (Na_V) of rat muscle by intracellular PAs and PAs have recently been found to modulate excitability of brain neocortical neurons by blocking neuronal Na_V channels. In this study, we compared the sensitivity of four different cloned mammalian Na_V isoforms to PAs to investigate whether PA block is a common feature of Na_V channel pharmacology. We find that outward Na⁺ current of muscle (Na_V1.4), heart (Na_V1.5), and neuronal (Na_V1.2, Na_V1.7) Na_V isoforms is blocked by PAs, suggesting that PA metabolism may be linked to modulation of action potential firing in numerous excitable tissues. Interestingly, the cardiac Na_V1.5 channel is more sensitive to PA block than other isoforms. Our results also indicate that rapid binding of PAs to blocking sites in the Na_V1.4 channel is restricted to access from the cytoplasmic side of the channel, but plasma membrane transport pathways for PA uptake may contribute to long-term Na_V channel modulation. PAs may also play a role in drug interactions since spermine attenuates the use-dependent effect of the lidocaine, a typical local anesthetic and anti-arrhythmic drug.     

Sensitivity of cloned muscle,heart and neuronal voltage-gated sodium channels to block by polyamines

Spermidine and spermine, are endogenous polyamines (PAs) that regulate cell growth and modulate the activity of numerous ion channel proteins. In particular, intracellular PAs are potent blockers of many different cation channels and are responsible for strong suppression of outward K⁺ current, a phenomenon known as inward rectification characteristic of a major class of K_IR K⁺ channels. We previously described block of heterologously expressed voltage-gated Na⁺ channels (Na_V) of rat muscle by intracellular PAs and PAs have recently been found to modulate excitability of brain neocortical neurons by blocking neuronal Na_V channels. In this study, we compared the sensitivity of four different cloned mammalian Na_V isoforms to PAs to investigate whether PA block is a common feature of Na_V channel pharmacology. We find that outward Na⁺ current of muscle (Na_V1.4), heart (Na_V1.5), and neuronal (Na_V1.2, Na_V1.7) Na_V isoforms is blocked by PAs, suggesting that PA metabolism may be linked to modulation of action potential firing in numerous excitable tissues. Interestingly, the cardiac Na_V1.5 channel is more sensitive to PA block than other isoforms. Our results also indicate that rapid binding of PAs to blocking sites in the Na_V1.4 channel is restricted to access from the cytoplasmic side of the channel, but plasma membrane transport pathways for PA uptake may contribute to long-term Na_V channel modulation. PAs may also play a role in drug interactions since spermine attenuates the use-dependent effect of the lidocaine, a typical local anesthetic and anti-arrhythmic drug.     

Energy distribution and ion selectivity of the bacterial potassium channel

The ion selectivity of the bacterial potassium channel K_CSA is explained upon comparing the energy characteristics of the interaction of cations (Li⁺, Na⁺, K⁺) with atoms of the selectivity filter of the protein pore. Quantum-chemical calculations reveal a deeper potential well for potassium ions, which accounts for preferred K⁺ permeation. It is shown that the conventional methods with AMBER, CHARMM, OPLS force fields in standard parametrization as well as partial re-parametrization give incorrect estimates of ion energy distribution in the channel.     

Ancient Origin of Four-Domain Voltage-gated Na<Superscript>+</Superscript> Channels Predates the Divergence of Animals and Fungi

The four-domain voltage-gated Na⁺ channels are believed to have arisen in multicellular animals, possibly during the evolution of the nervous system. Recent genomic studies reveal that many ion channels, including Na⁺ channels and Ca²⁺ channels previously thought to be restricted to animals, can be traced back to one of the unicellular ancestors of animals, Monosiga brevicollis. The eukaryotic supergroup Opisthokonta contains animals, fungi, and a diverse group of their unicellular relatives including M. brevicollis. Here, we demonstrate the presence of a putative voltage-gated Na⁺ channel homolog (TtrNa_V) in the apusozoan protist Thecamonas trahens, which belongs to the unicellular sister group to Opisthokonta. TtrNa_V displays a unique selectivity motif distinct from most animal voltage-gated Na⁺ channels. The identification of TtrNa_V suggests that voltage-gated Na⁺ channels might have evolved before the divergence of animals and fungi. Furthermore, our analyses reveal that Na_V channels have been lost independently in the amoeboid holozoan Capsaspora owczarzaki of the animal lineage and in several basal fungi. These findings provide novel insights into the evolution of four-domain voltage-gated ion channels, ion selectivity, and membrane excitability in the Opisthokonta lineage.     

Contribution of residues in second transmembrane domain of ASIC1a protein to ion selectivity

Acid-sensing ion channels (ASICs) are proton-gated cation-selective channels expressed in the peripheral and central nervous systems. The ion permeation pathway of ASIC1a is defined by residues 426–450 in the second transmembrane (TM2) segment. The gate, formed by the intersection of the TM2 segments, localizes near the extracellular boundary of the plasma membrane. We explored the contribution to ion permeation and selectivity of residues in the TM2 segment of ASIC1a. Studies of accessibility with positively charged methanethiosulfonate reagents suggest that the permeation pathway in the open state constricts below the gate, restricting the passage to large ions. Substitution of residues in the intracellular vestibule at positions 437, 438, 443, or 446 significantly increased the permeability to K⁺ versus Na⁺. ASIC1a shows a selectivity sequence for alkali metals of Na⁺>Li⁺>K⁺≫Rb⁺>Cs⁺. Alanine and cysteine substitutions at position 438 increased, to different extents, the relative permeability to Li⁺, K⁺, Rb⁺, and Cs⁺. For these mutants, ion permeation was not a function of the diameter of the nonhydrated ion, suggesting that Gly-438 encompasses an ion coordination site that is essential for ion selectivity. M437C and A443C mutants showed slightly increased permeability to K⁺, Rb⁺, and Cs⁺, suggesting that substitutions at these positions influence ion discrimination by altering molecular sieving. Our results indicate that ion selectivity is accomplished by the contribution of multiple sites in the pore of ASIC1a.     

Ion behavior in the selectivity filter of HCN1 channels

Hyperpolarization-activated cyclic-nucleotide gated channels (HCNs) are responsible for the generation of pacemaker currents (I_f or I_h) in cardiac and neuronal cells. Despite the overall structural similarity to voltage-gated potassium (Kv) channels, HCNs show much lower selectivity for K⁺ over Na⁺ ions. This increased permeability to Na⁺ is critical to their role in membrane depolarization. HCNs can also select between Na⁺ and Li⁺ ions. Here, we investigate the unique ion selectivity properties of HCNs using molecular-dynamics simulations. Our simulations suggest that the HCN1 pore is flexible and dilated compared with Kv channels with only one stable ion binding site within the selectivity filter. We also observe that ion coordination and hydration differ within the HCN1 selectivity filter compared with those in Kv and cyclic-nucleotide gated channels. Additionally, the C358T mutation further stabilizes the symmetry of the binding site and provides a more fit space for ion coordination, particularly for Li⁺.     

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.     

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na⁺ ions, and helps to explain the block of KcsA by internal Na⁺ ions.     

Of cilia,titin, and neurosteroids

Missense mutations at arginine residues in the S4 voltage-sensor domains of Na_V1.4 are an established cause of hypokalemic periodic paralysis, an inherited disorder of skeletal muscle involving recurrent episodes of weakness in conjunction with low serum K⁺. Expression studies in oocytes have revealed anomalous, hyperpolarization-activated gating pore currents in mutant channels. This aberrant gating pore conductance creates a small inward current at the resting potential that is thought to contribute to susceptibility to depolarization in low K⁺ during attacks of weakness. A critical component of this hypothesis is the magnitude of the gating pore conductance relative to other conductances that are active at the resting potential in mammalian muscle: large enough to favor episodes of paradoxical depolarization in low K⁺, yet not so large as to permanently depolarize the fiber. To improve the estimate of the specific conductance for the gating pore in affected muscle, we sequentially measured Na⁺ current through the channel pore, gating pore current, and gating charge displacement in oocytes expressing R669H, R672G, or wild-type Na_V1.4 channels. The relative conductance of the gating pore to that of the pore domain pathway for Na⁺ was 0.03%, which implies a specific conductance in muscle from heterozygous patients of ∼10 µS/cm² or 1% of the total resting conductance.Unexpectedly, our data also revealed a substantial decoupling between gating charge displacement and peak Na⁺ current for both R669H and R672G mutant channels. This decoupling predicts a reduced Na⁺ current density in affected muscle, consistent with the observations that the maximal dV/dt and peak amplitude of the action potential are reduced in fibers from patients with R672G and in a knock-in mouse model of R669H. The defective coupling between gating charge displacement and channel activation identifies a previously unappreciated mechanism that contributes to the reduced excitability of affected fibers seen with these mutations and possibly with other R/X mutations of S4 of Na_V, Ca_V, and K_V channels associated with human disease.