Voltage-gated Na + channels play an essential role in electrical signaling in the nervous system and are key pharmacological targets for a range of disorders. The recent solution of X-ray structures for the bacterial channel Na v Ab has provided an opportunity to study functional mechanisms at the atomic level. This channel's selectivity filter exhibits an EEEE ring sequence, characteristic of mammalian Ca 2+ , not Na + , channels. This raises the fundamentally important question: just what makes a Na + channel conduct Na + ions? Here we explore ion permeation on multimicrosecond timescales using the purpose-built Anton supercomputer. We isolate the likely protonation states of the EEEE ring and observe a striking flexibility of the filter that demonstrates the necessity for extended simulations to study conduction in this channel. We construct free energy maps to reveal complex multi-ion conduction via knock-on and "pass-by" mechanisms, involving concerted ion and glutamate side chain movements. Simulations in mixed ionic solutions reveal relative energetics for Na + , K + , and Ca 2+ within the pore that are consistent with the modest selectivity seen experimentally. We have observed conformational changes in the pore domain leading to asymmetrical collapses of the activation gate, similar to proposed inactivated structures of Na v Ab, with helix bending involving conserved residues that are critical for slow inactivation. These structural changes are shown to regulate access to fenestrations suggested to be pathways for lipophilic drugs and provide deeper insight into the molecular mechanisms connecting drug activity and slow inactivation.

Voltage-gated sodium (Na v ) channels are important targets in the treatment of a range of pathologies. Bacterial channels, for which crystal structures have been solved, exhibit modulation by local anesthetic and anti-epileptic agents, allowing molecular-level investigations into sodium channel-drug interactions. These structures reveal no basis for the "hinged lid"-based fast inactivation, seen in eukaryotic Na v channels. Thus, they enable examination of potential mechanisms of use-or state-dependent drug action based on activation gating, or slower pore-based inactivation processes. Multimicrosecond simulations of Na v Ab reveal high-affinity binding of benzocaine to F203 that is a surrogate for FS6, conserved in helix S6 of Domain IV of mammalian sodium channels, as well as low-affinity sites suggested to stabilize different states of the channel. Phenytoin exhibits a different binding distribution owing to preferential interactions at the membrane and water-protein interfaces. Two drug-access pathways into the pore are observed: via lateral fenestrations connecting to the membrane lipid phase, as well as via an aqueous pathway through the intracellular activation gate, despite being closed. These observations provide insight into drug modulation that will guide further developments of Na v inhibitors.bacterial sodium channel | drug binding

Increased extracellular proton concentrations during neurotransmission are converted to excitatory sodium influx by acid-sensing ion channels (ASICs). 10-fold sodium/potassium selectivity in ASICs has long been attributed to a central constriction in the channel pore, but experimental verification is lacking due to the sensitivity of this structure to conventional manipulations. Here, we explored the basis for ion selectivity by incorporating unnatural amino acids into the channel, engineering channel stoichiometry and performing free energy simulations. We observed no preference for sodium at the “GAS belt” in the central constriction. Instead, we identified a band of glutamate and aspartate side chains at the lower end of the pore that enables preferential sodium conduction.DOI: http://dx.doi.org/10.7554/eLife.24630.001

Membrane ion channels are the fundamental electrical components in the nervous system. Recent developments in X-ray crystallography and cryo-EM microscopy have revealed what these proteins look like in atomic detail but do not tell us how they function. Molecular dynamics simulations have progressed to the point that we can now simulate realistic molecular assemblies to produce quantitative calculations of the thermodynamic and kinetic quantities that control function. In this review, we summarize the state of atomistic simulation methods for ion channels to understand their conduction, activation, and drug modulation mechanisms. We are at a crossroads in atomistic simulation, where long time scale observation can provide unbiased exploration of mechanisms, supplemented by biased free energy methodologies. We illustrate the use of these approaches to describe ion conduction and selectivity in voltage-gated sodium and acid-sensing ion channels. Studies of channel gating present a significant challenge, as activation occurs on longer time scales. Enhanced sampling approaches can ensure convergence on minimum free energy pathways for activation, as illustrated here for pentameric ligand-gated ion channels that are principal to nervous system function and the actions of general anesthetics. We also examine recent studies of local anesthetic and antiepileptic drug binding to a sodium channel, revealing sites and pathways that may offer new targets for drug development. Modern simulations thus offer a range of molecular-level insights into ion channel function and modulation as a learning platform for mechanistic discovery and drug development.

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.