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

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.

Purification and Characterization of JZTx-14, a Potent Antagonist of Mammalian and Prokaryotic Voltage-Gated Sodium Channels

Exploring the interaction of ligands with voltage-gated sodium channels (Na_Vs) has advanced our understanding of their pharmacology. Herein, we report the purification and characterization of a novel non-selective mammalian and bacterial Na_Vs toxin, JZTx-14, from the venom of the spider Chilobrachys jingzhao. This toxin potently inhibited the peak currents of mammalian Na_V1.2–1.8 channels and the bacterial NaChBac channel with low IC₅₀ values (<1 µM), and it mainly inhibited the fast inactivation of the Na_V1.9 channel. Analysis of Na_V1.5/Na_V1.9 chimeric channel showed that the Na_V1.5 domain II S3–4 loop is involved in toxin association. Kinetics data obtained from studying toxin–Na_V1.2 channel interaction showed that JZTx-14 was a gating modifier that possibly trapped the channel in resting state; however, it differed from site 4 toxin HNTx-III by irreversibly blocking Na_V currents and showing state-independent binding with the channel. JZTx-14 might stably bind to a conserved toxin pocket deep within the Na_V1.2–1.8 domain II voltage sensor regardless of channel conformation change, and its effect on Na_Vs requires the toxin to trap the S3–4 loop in its resting state. For the NaChBac channel, JZTx-14 positively shifted its conductance-voltage (G–V) and steady-state inactivation relationships. An alanine scan analysis of the NaChBac S3–4 loop revealed that the 108th phenylalanine (F108) was the key residue determining the JZTx-14–NaChBac interaction. In summary, this study provided JZTx-14 with potent but promiscuous inhibitory activity on both the ancestor bacterial Na_Vs and the highly evolved descendant mammalian Na_Vs, and it is a useful probe to understand the pharmacology of Na_Vs.

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.

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.