Analysis of the selectivity filter of the voltage-gated sodium channel Na(v)Rh

The chemistry of electrical signaling in sodium channels from bacteria and beyond

Electrical signaling is essential for all fast processes in biology, but its molecular mechanisms have been uncertain. This review article focuses on studies of bacterial sodium channels in order to home in on the essential molecular and chemical mechanisms underlying transmembrane ion conductance and voltage-dependent gating without the overlay of complex protein interactions and regulatory mechanisms in mammalian sodium channels. This minimalist approach has yielded a nearly complete picture of sodium channel function at the atomic level that are mostly conserved in mammalian sodium channels, including sodium selectivity and conductance, voltage sensing and activation, electromechanical coupling to pore opening and closing, slow inactivation, and pathogenic dysfunction in a debilitating channelopathy. Future studies of nature's simplest sodium channels may continue to yield key insights into the fundamental molecular and chemical principles of their function and further elucidate the chemical basis of electrical signaling.

Pore residues of transient receptor potential channels canonical 1 and 4 heteromer determine channel properties

Transient receptor potential channels canonical 1 and 4 (TRPC1 and TRPC4) are proteins belonging to the same TRPC channel family, and the two are known to form a heterotetrameric channel. TRPC4 can form a homotetrameric, nonselective cation channel by itself, but the involvement of the TRPC1 subunit changes several major characteristics of the channel. In this study, we focused on the pore region (selectivity filter, pore helix, and S6 helix) of TRPC1 and TRPC4 as a determinant of the identity and characteristics of a heteromeric TRPC1/4 channel: decreased calcium permeability of the channel and outward-rectifying current-voltage (I-V) curve. Mutants and chimeras of the pore residues were created, and their currents were recorded using whole cell patch clamp. The lower gate mutants of TRPC4 exhibited diminished calcium permeability as measured by GCaMP6 fluorescence. Also, chimeric channels substituting the pore region of TRPC1 to TRPC4 were made to locate the pore region that is critical in the production of an outward-rectifying I-V curve characteristic of TRPC1/4 heteromeric channels.NEW & NOTEWORTHY Heteromer research has been a challenging field due to lack of structural studies. Using chimeras and single mutants, we present evidence that the pore region of TRPC1/4 heteromer contributes to determining the channel's characteristics such as calcium permeability, I-V curve, and conductance.

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.

Abdominal TRPV1 channel desensitization enhances stress-induced hyperthermia during social stress in rats

AIMS

In rats, stress-induced hyperthermia caused by social interaction depends on brown adipose tissue (BAT) thermogenesis and peripheral vasoconstriction. However, the peripheral mechanisms responsible for regulating the level of hyperthermia during social stress are still unknown. The transient receptor potential vanilloid 1 (TRPV1) subfamily, expressed in sensory and visceral neurons, can serve as a thermoreceptor. Here, we tested the hypothesis that the abdominal TRPV1 is essential in regulating stress-induced hyperthermia during social stress.

MAIN METHODS

Male Wistar rats received an intraperitoneal injection of Resiniferatoxin (RTX) - an ultra-potent capsaicin analog, (i.e., to desensitize the TRPV1 channels) or vehicle. Seven days later, we evaluated the effects of abdominal TRPV1 channels desensitization on core body temperature (CBT), brown adipose tissue (BAT) temperature, tail skin temperature, and heart rate (HR) of rats subjected to a social stress protocol.

KEY FINDINGS

We found abdominal TRPV1 desensitization increased CBT and BAT temperature but did not change tail skin temperature and HR during rest. However, under social stress, we found that abdominal TRPV1 desensitization heightened the increase in CBT and BAT caused by stress. Also, it abolished the increase in tail skin temperature that occurs during and after social stress. TRPV1 desensitization also delayed the HR recovery after the exposure to the social stress.

SIGNIFICANCE

These results show that abdominal TRPV1 channels desensitization heightens stress-induced hyperthermia, causing heat dissipation during and after social stress, enabling optimal thermal control during social encounters.

Mechanisms and significance of Ca2+ entry through TRPC channels

The transient receptor potential (TRP) superfamily of plasma membrane cation channels has been recognized as a signaling hub in highly diverse settings of human physiopathology. In the past three decades of TRP research, attention was focused mainly on the channels Ca2+ signaling function, albeit additional cellular functions, aside of providing a Ca2+ entry pathway, have been identified. Our understanding of Ca2+ signaling by TRP proteins has recently been advanced by a gain in high-resolution structure information on these pore complexes, and by the development of novel tools to investigate their role in spatiotemporal Ca2+ handling. This review summarizes recent discoveries as well as remaining, unresolved aspects of the canonical subfamily of transient receptor potential channels (TRPC) research. We aim at a concise overview on current mechanistic concepts of Ca2+ entry through- and Ca2+ signaling by TRPC channels.

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.

Extremely Potent Block of Bacterial Voltage-Gated Sodium Channels by µ-Conotoxin PIIIA

µ-Conotoxin PIIIA, in the sub-picomolar, range inhibits the archetypal bacterial sodium channel NaChBac (NavBh) in a voltage- and use-dependent manner. Peptide µ-conotoxins were first recognized as potent components of the venoms of fish-hunting cone snails that selectively inhibit voltage-gated skeletal muscle sodium channels, thus preventing muscle contraction. Intriguingly, computer simulations predicted that PIIIA binds to prokaryotic channel NavAb with much higher affinity than to fish (and other vertebrates) skeletal muscle sodium channel (Nav 1.4). Here, using whole-cell voltage clamp, we demonstrate that PIIIA inhibits NavBac mediated currents even more potently than predicted. From concentration-response data, with [PIIIA] varying more than 6 orders of magnitude (10⁻¹² to 10⁻⁵ M), we estimated an IC₅₀ = ~5 pM, maximal block of 0.95 and a Hill coefficient of 0.81 for the inhibition of peak currents. Inhibition was stronger at depolarized holding potentials and was modulated by the frequency and duration of the stimulation pulses. An important feature of the PIIIA action was acceleration of macroscopic inactivation. Docking of PIIIA in a NaChBac (NavBh) model revealed two interconvertible binding modes. In one mode, PIIIA sterically and electrostatically blocks the permeation pathway. In a second mode, apparent stabilization of the inactivated state was achieved by PIIIA binding between P2 helices and trans-membrane S5s from adjacent channel subunits, partially occluding the outer pore. Together, our experimental and computational results suggest that, besides blocking the channel-mediated currents by directly occluding the conducting pathway, PIIIA may also change the relative populations of conducting (activated) and non-conducting (inactivated) states.

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.

Cross-kingdom auxiliary subunit modulation of a voltage-gated sodium channel

Voltage-gated, sodium ion-selective channels (NaV) generate electrical signals contributing to the upstroke of the action potential in animals. NaVs are also found in bacteria and are members of a larger family of tetrameric voltage-gated channels that includes CaVs, KVs, and NaVs. Prokaryotic NaVs likely emerged from a homotetrameric Ca2+-selective voltage-gated progenerator, and later developed Na+ selectivity independently. The NaV signaling complex in eukaryotes contains auxiliary proteins, termed beta (β) subunits, which are potent modulators of the expression profiles and voltage-gated properties of the NaV pore, but it is unknown whether they can functionally interact with prokaryotic NaV channels. Herein, we report that the eukaryotic NaVβ1-subunit isoform interacts with and enhances the surface expression as well as the voltage-dependent gating properties of the bacterial NaV, NaChBac in Xenopus oocytes. A phylogenetic analysis of the β-subunit gene family proteins confirms that these proteins appeared roughly 420 million years ago and that they have no clear homologues in bacterial phyla. However, a comparison between eukaryotic and bacterial NaV structures highlighted the presence of a conserved fold, which could support interactions with the β-subunit. Our electrophysiological, biochemical, structural, and bioinformatics results suggests that the prerequisites for β-subunit regulation are an evolutionarily stable and intrinsic property of some voltage-gated channels.

Modeling the human Nav1.5 sodium channel: structural and mechanistic insights of ion permeation and drug blockade

Abnormalities in the human Nav1.5 (hNav1.5) voltage-gated sodium ion channel (VGSC) are associated with a wide range of cardiac problems and diseases in humans. Current structural models of hNav1.5 are still far from complete and, consequently, their ability to study atomistic interactions of this channel is very limited. Here, we report a comprehensive atomistic model of the hNav1.5 ion channel, constructed using homology modeling technique and refined through long molecular dynamics simulations (680 ns) in the lipid membrane bilayer. Our model was comprehensively validated by using reported mutagenesis data, comparisons with previous models, and binding to a panel of known hNav1.5 blockers. The relatively long classical MD simulation was sufficient to observe a natural sodium permeation event across the channel's selectivity filters to reach the channel's central cavity, together with the identification of a unique role of the lysine residue. Electrostatic potential calculations revealed the existence of two potential binding sites for the sodium ion at the outer selectivity filters. To obtain further mechanistic insight into the permeation event from the central cavity to the intracellular region of the channel, we further employed "state-of-the-art" steered molecular dynamics (SMD) simulations. Our SMD simulations revealed two different pathways through which a sodium ion can be expelled from the channel. Further, the SMD simulations identified the key residues that are likely to control these processes. Finally, we discuss the potential binding modes of a panel of known hNav1.5 blockers to our structural model of hNav1.5. We believe that the data presented here will enhance our understanding of the structure-property relationships of the hNav1.5 ion channel and the underlying molecular mechanisms in sodium ion permeation and drug interactions. The results presented here could be useful for designing safer drugs that do not block the hNav1.5 channel.

Voltage-gated sodium channels viewed through a structural biology lens

Voltage-gated sodium (Nav) channels initiate and propagate action potentials in excitable cells, and are frequently dysregulated or mutated in human disease. Despite decades of intense physiological and biophysical research, eukaryotic Nav channels have so far eluded high-resolution structure determination because of their biochemical complexity. Recently, simpler bacterial voltage-gated sodium (BacNav) channels have provided templates to understand the structural basis of voltage-dependent activation, inactivation, ion selectivity, and drug block in eukaryotic Nav and related voltage-gated calcium (Cav) channels. Further breakthroughs employing BacNav channels have also enabled visualization of bound small molecule modulators that can guide the rational design of next generation therapeutics. This review will highlight the emerging structural biology of BacNav channels and its contribution to our understanding of the gating, ion selectivity, and pharmacological regulation of eukaryotic Nav (and Cav) 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.

Lysine and the Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_v) channels are critical in the generation and transmission of neuronal signals in mammals. The crystal structures of several prokaryotic Na_v channels determined in recent years inspire the mechanistic studies on their selection upon the permeable cations (especially between Na⁺ and K⁺ ions), a property that is proposed to be mainly determined by residues in the selectivity filter. However, the mechanism of cation selection in mammalian Na_v channels lacks direct explanation at atomic level due to the difference in amino acid sequences between mammalian and prokaryotic Na_v homologues, especially at the constriction site where the DEKA motif has been identified to determine the Na⁺/K⁺ selectivity in mammalian Na_v channels but is completely absent in the prokaryotic counterparts. Among the DEKA residues, Lys is of the most importance since its mutation to Arg abolishes the Na⁺/K⁺ selectivity. In this work, we modeled the pore domain of mammalian Na_v channels by mutating the four residues at the constriction site of a prokaryotic Na_v channel (Na_vRh) to DEKA, and then mechanistically investigated the contribution of Lys in cation selection using molecular dynamics simulations. The DERA mutant was generated as a comparison to understand the loss of ion selectivity caused by the K-to-R mutation. Simulations and free energy calculations on the mutants indicate that Lys facilitates Na⁺/K⁺ selection by electrostatically repelling the cation to a highly Na⁺-selective location sandwiched by the carboxylate groups of Asp and Glu at the constriction site. In contrast, the electrostatic repulsion is substantially weakened when Lys is mutated to Arg, because of two intrinsic properties of the Arg side chain: the planar geometric design and the sparse charge distribution of the guanidine group.

Understanding Sodium Channel Function and Modulation Using Atomistic Simulations of Bacterial Channel Structures

Sodium channels are chief proteins involved in electrical signaling in the nervous system, enabling critical functions like heartbeat and brain activity. New high-resolution X-ray structures for bacterial sodium channels have created an opportunity to see how these proteins operate at the molecular level. An important challenge to overcome is establishing relationships between the structures and functions of mammalian and bacterial channels. Bacterial sodium channels are known to exhibit the main structural features of their mammalian counterparts, as well as several key functional characteristics, including selective ion conduction, voltage-dependent gating, pore-based inactivation and modulation by local anesthetic, antiarrhythmic and antiepileptic drugs. Simulations have begun to shed light on each of these features in the past few years. Despite deviations in selectivity signatures for bacterial and mammalian channels, simulations have uncovered the nature of the multiion conduction mechanism associated with Na(+) binding to a high-field strength site established by charged glutamate side chains. Simulations demonstrated a surprising level of flexibility of the protein, showing that these side chains are active participants in the permeation process. They have also uncovered changes in protein structure, leading to asymmetrical collapses of the activation gate that have been proposed to correspond to inactivated structures. These observations offer the potential to examine the mechanisms of state-dependent drug activity, focusing on pore-blocking and pore-based slow inactivation in bacterial channels, without the complexities of inactivation on multiple timescales seen in eukaryotic channels. Simulations have provided molecular views of the interactions of drugs, consistent with sites predicted in mammalian channels, as well as a wealth of other sites as potential new drug targets. In this chapter, we survey the new insights into sodium channel function that have emerged from studies of simpler bacterial channels, which provide an excellent learning platform, and promising avenues for mechanistic discovery and pharmacological development.

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.

Molecular Modeling of the Structural and Dynamical Changes in Calcium Channel TRPV5 Induced by the African-Specific A563T Variation

Transient receptor potential cation channels, vanilloid subfamily, member 5 (TRPV5) plays a key role in active Ca(2+) reabsorption in the kidney. Variations in TRPV5 occur at high frequency in African populations and may contribute to their higher efficiency of Ca(2+) reabsorption. One of the African specific variations, A563T, exhibits increased Ca(2+) transport ability. However, it is unclear how this variation influences the channel pore. On the basis of the structure of TRPV1, a TRPV5 model was generated to simulate the structural and dynamical changes induced by the A563T variation. On the basis of this model, amino acid residue 563 interacts with V540, which is one residue away from the key residue, D542, involved in Ca(2+) selectivity and Mg(2+) blockade. The A563T variation increases secondary structure stability and reduces dynamical motion of D542. In addition, the A563T variation alters the electrostatic potential of the outer surface of the pore. Differences in contact between selective filter residues and residue 563 and in electrostatic potential between the two TRPV5 variants were also observed in another model derived from an alternative alignment in the selective filters between TRPV5 and TRPV1. These findings indicate that the A563T variation induces structural, dynamical, and electrostatic changes in the TRPV5 pore, providing structural insight into the functional alterations associated with the A563T variation.

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.

Mechanism of Ion Permeation in Mammalian Voltage-Gated Sodium Channels

Recent determination of the crystal structures of bacterial voltage-gated sodium (NaV) channels have raised hopes that modeling of the mammalian counterparts could soon be achieved. However, there are substantial differences between the pore domains of the bacterial and mammalian NaV channels, which necessitates careful validation of mammalian homology models constructed from the bacterial NaV structures. Such a validated homology model for the NaV1.4 channel was constructed recently using the extensive mutagenesis data available for binding of μ-conotoxins. Here we use this NaV1.4 model to study the ion permeation mechanism in mammalian NaV channels. Linking of the DEKA residues in the selectivity filter with residues in the neighboring domains is found to be important for keeping the permeation pathway open. Molecular dynamics simulations and potential of mean force calculations reveal that there is a binding site for a Na+ ion just inside the DEKA locus, and 1-2 Na+ ions can occupy the vestibule near the EEDD ring. These sites are separated by a low free energy barrier, suggesting that inward conduction occurs when a Na+ ion in the vestibule goes over the free energy barrier and pushes the Na+ ion in the filter to the intracellular cavity, consistent with the classical knock-on mechanism. The NaV1.4 model also provides a good description of the observed Na+/K+ selectivity.

Deciphering voltage-gated Na(+) and Ca(2+) channels by studying prokaryotic ancestors

Voltage-gated sodium channels (NaVs) and calcium channels (CaVs) are involved in electrical signaling, contraction, secretion, synaptic transmission, and other physiological processes activated in response to depolarization. Despite their physiological importance, the structures of these closely related proteins have remained elusive because of their size and complexity. Bacterial NaVs have structures analogous to a single domain of eukaryotic NaVs and CaVs and are their likely evolutionary ancestor. Here we review recent work that has led to new understanding of NaVs and CaVs through high-resolution structural studies of their prokaryotic ancestors. New insights into their voltage-dependent activation and inactivation, ion conductance, and ion selectivity provide realistic structural models for the function of these complex membrane proteins at the atomic level.

Voltage-Gated Sodium Channels: Structure, Function, Pharmacology, and Clinical Indications

The tremendous therapeutic potential of voltage-gated sodium channels (Na(v)s) has been the subject of many studies in the past and is of intense interest today. Na(v)1.7 channels in particular have received much attention recently because of strong genetic validation of their involvement in nociception. Here we summarize the current status of research in the Na(v) field and present the most relevant recent developments with respect to the molecular structure, general physiology, and pharmacology of distinct Na(v) channel subtypes. We discuss Na(v) channel ligands such as small molecules, toxins isolated from animal venoms, and the recently identified Na(v)1.7-selective antibody. Furthermore, we review eight characterized ligand binding sites on the Na(v) channel α subunit. Finally, we examine possible therapeutic applications of Na(v) ligands and provide an update on current clinical studies.

Theoretical and simulation studies on voltage-gated sodium channels

Voltage-gated sodium (Na_v) channels are indispensable membrane elements for the generation and propagation of electric signals in excitable cells. The successes in the crystallographic studies on prokaryotic Na_v channels in recent years greatly promote the mechanistic investigation of these proteins and their eukaryotic counterparts. In this paper, we mainly review the progress in computational studies, especially the simulation studies, on these proteins in the past years.

Bacterial voltage-gated sodium channels (BacNa(V)s) from the soil, sea, and salt lakes enlighten molecular mechanisms of electrical signaling and pharmacology in the brain and heart

Voltage-gated sodium channels (Na(V)s) provide the initial electrical signal that drives action potential generation in many excitable cells of the brain, heart, and nervous system. For more than 60years, functional studies of Na(V)s have occupied a central place in physiological and biophysical investigation of the molecular basis of excitability. Recently, structural studies of members of a large family of bacterial voltage-gated sodium channels (BacNa(V)s) prevalent in soil, marine, and salt lake environments that bear many of the core features of eukaryotic Na(V)s have reframed ideas for voltage-gated channel function, ion selectivity, and pharmacology. Here, we analyze the recent advances, unanswered questions, and potential of BacNa(V)s as templates for drug development efforts.

Molecular dynamics study of binding of µ-conotoxin GIIIA to the voltage-gated sodium channel Na(v)1.4

Homology models of mammalian voltage-gated sodium (NaV) 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 NaV1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to NaV1.4 to validate the model. The initial poses for the NaV1.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 NaV1 channels are highly homologous, the model constructed for NaV1.4 will provide an excellent template for other NaV1 channels.

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

Rapid and selective ion transport is essential for the generation and regulation of electrical signaling pathways in living organisms. Here, we use molecular dynamics (MD) simulations with an applied membrane potential to investigate the ion flux of bacterial sodium channel NaVMs. 5.9 µs simulations with 500 mM NaCl suggest different mechanisms for inward and outward flux. The predicted inward conductance rate of ∼27±3 pS, agrees with experiment. The estimated outward conductance rate is 15±3 pS, which is considerably lower. Comparing inward and outward flux, the mean ion dwell time in the selectivity filter (SF) is prolonged from 13.5±0.6 ns to 20.1±1.1 ns. Analysis of the Na+ distribution revealed distinct patterns for influx and efflux events. In 32.0±5.9% of the simulation time, the E53 side chains adopted a flipped conformation during outward conduction, whereas this conformational change was rarely observed (2.7±0.5%) during influx. Further, simulations with dihedral restraints revealed that influx is less affected by the E53 conformational flexibility. In contrast, during outward conduction, our simulations indicate that the flipped E53 conformation provides direct coordination for Na+. The free energy profile (potential of mean force calculations) indicates that this conformational change lowers the putative barriers between sites SCEN and SHFS during outward conduction. We hypothesize that during an action potential, the increased Na+ outward transition propensities at depolarizing potentials might increase the probability of E53 conformational changes in the SF. Subsequently, this might be a first step towards initiating slow inactivation.

The twins K+ and Na+ in plants

In the earth's crust and in seawater, K(+) and Na(+) are by far the most available monovalent inorganic cations. Physico-chemically, K(+) and Na(+) are very similar, but K(+) is widely used by plants whereas Na(+) can easily reach toxic levels. Indeed, salinity is one of the major and growing threats to agricultural production. In this article, we outline the fundamental bases for the differences between Na(+) and K(+). We present the foundation of transporter selectivity and summarize findings on transporters of the HKT type, which are reported to transport Na(+) and/or Na(+) and K(+), and may play a central role in Na(+) utilization and detoxification in plants. Based on the structural differences in the hydration shells of K(+) and Na(+), and by comparison with sodium channels, we present an ad hoc mechanistic model that can account for ion permeation through HKTs.

K(+) and Na(+) conduction in selective and nonselective ion channels via molecular dynamics simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.

Bacterial sodium channels: models for eukaryotic sodium and calcium channels

Eukaryotic sodium and calcium channels are made up of four linked homologous but different transmembrane domains. Bacteria express sodium channels comprised of four identical subunits, each being analogous to a single homologous domain of their eukaryotic counterparts. Key elements of primary structure are conserved between bacterial and eukaryotic sodium and calcium channels. The simple protein structure of the bacterial channels has allowed extensive structure-function probes of key regions as well as allowing determination of several X-ray crystallographic structures of these channels. The structures have revealed novel features of sodium and calcium channel pores and elucidated the structural importance of many of the conserved features of primary sequence. The structural information has also formed the basis for computational studies probing the basis for sodium and calcium selectivity and gating.

Direct evidence that scorpion α-toxins (site-3) modulate sodium channel inactivation by hindrance of voltage-sensor movements

The position of the voltage-sensing transmembrane segment, S4, in voltage-gated ion channels as a function of voltage remains incompletely elucidated. Site-3 toxins bind primarily to the extracellular loops connecting transmembrane helical segments S1-S2 and S3-S4 in Domain 4 (D4) and S5-S6 in Domain 1 (D1) and slow fast-inactivation of voltage-gated sodium channels. As S4 of the human skeletal muscle voltage-gated sodium channel, hNav1.4, moves in response to depolarization from the resting to the inactivated state, two D4S4 reporters (R2C and R3C, Arg1451Cys and Arg1454Cys, respectively) move from internal to external positions as deduced by reactivity to internally or externally applied sulfhydryl group reagents, methane thiosulfonates (MTS). The changes in reporter reactivity, when cycling rapidly between hyperpolarized and depolarized voltages, enabled determination of the positions of the D4 voltage-sensor and of its rate of movement. Scorpion α-toxin binding impedes D4S4 segment movement during inactivation since the modification rates of R3C in hNav1.4 with methanethiosulfonate (CH3SO2SCH2CH2R, where R = -N(CH3)3 (+) trimethylammonium, MTSET) and benzophenone-4-carboxamidocysteine methanethiosulfonate (BPMTS) were slowed ~10-fold in toxin-modified channels. Based upon the different size, hydrophobicity and charge of the two reagents it is unlikely that the change in reactivity is due to direct or indirect blockage of access of this site to reagent in the presence of toxin (Tx), but rather is the result of inability of this segment to move outward to the normal extent and at the normal rate in the toxin-modified channel. Measurements of availability of R3C to internally applied reagent show decreased access (slower rates of thiol reaction) providing further evidence for encumbered D4S4 movement in the presence of toxins consistent with the assignment of at least part of the toxin binding site to the region of D4S4 region of the voltage-sensor module.

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 (NaVs) are central elements of cellular excitation. Notwithstanding advances from recent bacterial NaV (BacNaV) structures, key questions about gating and ion selectivity remain. Here, we present a closed conformation of NaVAe1p, a pore-only BacNaV derived from NaVAe1, a BacNaV 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 BacNaVs show that two elements defined by the NaVAe1p structure, an S6 activation gate position and the cytoplasmic tail "neck", are central to BacNaV gating. The structure also reveals the selectivity filter ion entry site, termed the "outer ion" site. Comparison with mammalian voltage-gated calcium channel (CaV) selectivity filters, together with functional studies, shows that this site forms a previously unknown determinant of CaV high-affinity calcium binding. Our findings underscore commonalities between BacNaVs and eukaryotic voltage-gated channels and provide a framework for understanding gating and ion permeation in this superfamily.

Bimodal voltage dependence of TRPA1: mutations of a key pore helix residue reveal strong intrinsic voltage-dependent inactivation

Transient receptor potential A1 (TRPA1) is implicated in somatosensory processing and pathological pain sensation. Although not strictly voltage-gated, ionic currents of TRPA1 typically rectify outwardly, indicating channel activation at depolarized membrane potentials. However, some reports also showed TRPA1 inactivation at high positive potentials, implicating voltage-dependent inactivation. Here we report a conserved leucine residue, L906, in the putative pore helix, which strongly impacts the voltage dependency of TRPA1. Mutation of the leucine to cysteine (L906C) converted the channel from outward to inward rectification independent of divalent cations and irrespective to stimulation by allyl isothiocyanate. The mutant, but not the wild-type channel, displayed exclusively voltage-dependent inactivation at positive potentials. The L906C mutation also exhibited reduced sensitivity to inhibition by TRPA1 blockers, HC030031 and ruthenium red. Further mutagenesis of the leucine to all natural amino acids individually revealed that most substitutions at L906 (15/19) resulted in inward rectification, with exceptions of three amino acids that dramatically reduced channel activity and one, methionine, which mimicked the wild-type channel. Our data are plausibly explained by a bimodal gating model involving both voltage-dependent activation and inactivation of TRPA1. We propose that the key pore helix residue, L906, plays an essential role in responding to the voltage-dependent gating.

The mechanism of Na⁺/K⁺ selectivity in mammalian voltage-gated sodium channels based on molecular dynamics simulation

Voltage-gated sodium (Nav) 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 Nav 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 Nav channels, the lack of Lys at the constriction site of prokaryotic Nav channels limits how much can be learned about the Na⁺/K⁺ selectivity in mammalian Nav channels. In this work, we modeled the mammalian Nav channel by mutating the key residues at the constriction site in a prokaryotic Nav channel (NavRh) to its mammalian counterpart. By simulating the mutant structure, we found that the Na⁺ preference in mammalian Nav channels is collaboratively achieved by the deselection from Lys and the selection from Asp/Glu within the constriction site.

Conduction in a biological sodium selective channel

The crystal structure of NavAb, a bacterial voltage gated Na(+) channel, exhibits a selectivity filter (SF) wider than that of K(+) channels. This new structure provides the opportunity to explore the mechanism of conduction and help rationalize its selectivity for sodium. Recent molecular dynamics (MD) simulations of single- and two-ion permeation processes have revealed that a partially hydrated Na(+) permeates the channel by exploring three SF binding sites while being loosely coupled to other ions and/or water molecules; a finding that differs significantly from the behavior of K(+) selective channels. Herein, we present results derived from a combination of metadynamics and voltage-biased MD simulations that throws more light on the nature of the Na(+) conduction mechanism. Conduction under 0 mV bias explores several distinct pathways involving the binding of two ions to three possible SF sites. While these pathways are very similar to those observed in the presence of a negative potential (inward conduction), a completely different mechanism operates for outward conduction at positive potentials.