Ion conduction and conformational flexibility of a bacterial voltage-gated sodium 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.

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

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

Simulating Current-Voltage Relationships for a Narrow Ion Channel Using the Weighted Ensemble Method

Ion channels are responsible for a myriad of fundamental biological processes via their role in controlling the flow of ions through water-filled membrane-spanning pores in response to environmental cues. Molecular simulation has played an important role in elucidating the mechanism of ion conduction, but connecting atomistically detailed structural models of the protein to electrophysiological measurements remains a broad challenge due to the computational cost of reaching the necessary time scales. Here, we introduce an enhanced sampling method for simulating the conduction properties of narrow ion channels using the Weighted ensemble (WE) sampling approach. We demonstrate the application of this method to calculate the current–voltage relationship as well as the nonequilibrium ion distribution at steady-state of a simple model ion channel. By direct comparisons with long brute force simulations, we show that the WE simulations rigorously reproduce the correct long-time scale kinetics of the system and are capable of determining these quantities using significantly less aggregate simulation time under conditions where permeation events are rare.

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 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.

K+ Block Is the Mechanism of Functional Asymmetry in Bacterial Na(v) Channels

Crystal structures of several bacterial Na_v channels have been recently published and molecular dynamics simulations of ion permeation through these channels are consistent with many electrophysiological properties of eukaryotic channels. Bacterial Na_v channels have been characterized as functionally asymmetric, and the mechanism of this asymmetry has not been clearly understood. To address this question, we combined non-equilibrium simulation data with two-dimensional equilibrium unperturbed landscapes generated by umbrella sampling and Weighted Histogram Analysis Methods for multiple ions traversing the selectivity filter of bacterial Na_vAb channel. This approach provided new insight into the mechanism of selective ion permeation in bacterial Na_v channels. The non-equilibrium simulations indicate that two or three extracellular K⁺ ions can block the entrance to the selectivity filter of Na_vAb in the presence of applied forces in the inward direction, but not in the outward direction. The block state occurs in an unstable local minimum of the equilibrium unperturbed free-energy landscape of two K⁺ ions that can be ‘locked’ in place by modest applied forces. In contrast to K⁺, three Na⁺ ions move favorably through the selectivity filter together as a unit in a loose “knock-on” mechanism of permeation in both inward and outward directions, and there is no similar local minimum in the two-dimensional free-energy landscape of two Na⁺ ions for a block state. The useful work predicted by the non-equilibrium simulations that is required to break the K⁺ block is equivalent to large applied potentials experimentally measured for two bacterial Na_v channels to induce inward currents of K⁺ ions. These results illustrate how inclusion of non-equilibrium factors in the simulations can provide detailed information about mechanisms of ion selectivity that is missing from mechanisms derived from either crystal structures or equilibrium unperturbed free-energy landscapes.

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

Refining the treatment of membrane proteins by coarse-grained models

Obtaining a quantitative description of the membrane proteins stability is crucial for understanding many biological processes. However the advance in this direction has remained a major challenge for both experimental studies and molecular modeling. One of the possible directions is the use of coarse-grained models but such models must be carefully calibrated and validated. Here we use a recent progress in benchmark studies on the energetics of amino acid residue and peptide membrane insertion and membrane protein stability in refining our previously developed coarse-grained model (Vicatos et al., Proteins 2014;82:1168). Our refined model parameters were fitted and/or tested to reproduce water/membrane partitioning energetics of amino acid side chains and a couple of model peptides. This new model provides a reasonable agreement with experiment for absolute folding free energies of several β-barrel membrane proteins as well as effects of point mutations on a relative stability for one of those proteins, OmpLA. The consideration and ranking of different rotameric states for a mutated residue was found to be essential to achieve satisfactory agreement with the reference data.

TMalphaDB and TMbetaDB: web servers to study the structural role of sequence motifs in α-helix and β-barrel domains of membrane proteins

Background

Membrane proteins represent over 25 % of human protein genes and account for more than 60 % of drug targets due to their accessibility from the extracellular environment. The increasing number of available crystal structures of these proteins in the Protein Data Bank permits an initial estimation of their structural properties.

Description

We have developed two web servers—TMalphaDB for α-helix bundles and TMbetaDB for β-barrels—to analyse the growing repertoire of available crystal structures of membrane proteins. TMalphaDB and TMbetaDB permit to search for these specific sequence motifs in a non-redundant structure database of transmembrane segments and quantify structural parameters such as ϕ and ψ backbone dihedral angles, χ₁ side chain torsion angle, unit bend and unit twist.

Conclusions

The structural information offered by TMalphaDB and TMbetaDB permits to quantify structural distortions induced by specific sequence motifs, and to elucidate their role in the 3D structure. This specific structural information has direct implications in homology modeling of the growing sequences of membrane proteins lacking experimental structure. TMalphaDB and TMbetaDB are freely available at http://lmc.uab.cat/TMalphaDB and http://lmc.uab.cat/TMbetaDB.

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.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

A structural, functional, and computational analysis suggests pore flexibility as the base for the poor selectivity of CNG channels

Cyclic nucleotide-gated (CNG) ion channels, despite a significant homology with the highly selective K(+) channels, do not discriminate among monovalent alkali cations and are permeable also to several organic cations. We combined electrophysiology, molecular dynamics (MD) simulations, and X-ray crystallography to demonstrate that the pore of CNG channels is highly flexible. When a CNG mimic is crystallized in the presence of a variety of monovalent cations, including Na(+), Cs(+), and dimethylammonium (DMA(+)), the side chain of Glu66 in the selectivity filter shows multiple conformations and the diameter of the pore changes significantly. MD simulations indicate that Glu66 and the prolines in the outer vestibule undergo large fluctuations, which are modulated by the ionic species and the voltage. This flexibility underlies the coupling between gating and permeation and the poor ionic selectivity of CNG channels.

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.

Improving the description of interactions between Ca2+ and protein carboxylate groups, including γ-carboxyglutamic acid: revised CHARMM22* parameters

We show that the CHARMM22* force-field over-binds the interaction between aqueous carboxylates and Ca²⁺, and introduce a modification that can recover experimentally-determined binding free energies for these systems.

Effects of the protonation state of the EEEE motif of a bacterial Na(+)-channel on conduction and pore structure

A distinctive feature of prokaryotic Na(+)-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na(+) ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na(+) ions. These higher-energy barriers preclude Na(+) ions to permeate the selectivity filter of prokaryotic Na(+)-channels when one or more Glu177 residues are protonated.

Ion selectivity strategies of sodium channel selectivity filters

CONSPECTUS: Sodium ion channels selectively transport Na(+) cations across the cell membrane. These integral parts of the cell machinery are implicated in regulating the cardiac, skeletal and smooth muscle contraction, nerve impulses, salt and water homeostasis, as well as pain and taste perception. Their malfunction often results in various channelopathies of the heart, brain, skeletal muscles, and lung; thus, sodium channels are key drug targets for various disorders including cardiac arrhythmias, heart attack, stroke, migraine, epilepsy, pain, cancer, and autoimmune disorders. The ability of sodium channels to discriminate the native Na(+) among other competing ions in the surrounding fluids is crucial for proper cellular functions. The selectivity filter (SF), the narrowest part of the channel's open pore, lined with amino acid residues that specifically interact with the permeating ion, plays a major role in determining Na(+) selectivity. Different sodium channels have different SFs, which vary in the symmetry, number, charge, arrangement, and chemical type of the metal-ligating groups and pore size: epithelial/degenerin/acid-sensing ion channels have generally trimeric SFs lined with three conserved neutral serines and/or backbone carbonyls; eukaryotic sodium channels have EKEE, EEKE, DKEA, and DEKA SFs with an invariant positively charged lysine from the second or third domain; and bacterial voltage-gated sodium (Nav) channels exhibit symmetrical EEEE SFs, reminiscent of eukaryotic voltage-gated calcium channels. How do these different sodium channel SFs achieve high selectivity for Na(+) over its key rivals, K(+) and Ca(2+)? What factors govern the metal competition in these SFs and which of these factors are exploited to achieve Na(+) selectivity in the different sodium channel SFs? The free energies for replacing K(+) or Ca(2+) bound inside different model SFs with Na(+), evaluated by a combination of density functional theory and continuum dielectric calculations, have shed light on these questions. The SFs of epithelial and eukaryotic Nav channels select Na(+) by providing an optimal number and ligating strength of metal ligands as well as a rigid pore whose size fits the cognate Na(+) ideally. On the other hand, the SFs of bacterial Nav channels select Na(+), as the protein matrix attenuates ion-protein interactions relative to ion-solvent interactions by enlarging the pore and allowing water to enter, so the ion interacts indirectly with the conserved glutamates via bridging water molecules. This shows how these various SFs have adapted to the specific physicochemical properties of the native ion, using different strategies to select Na(+) among its contenders.

Permeation and dynamics of an open-activated TRPV1 channel

Transient receptor potential (TRP) ion channels constitute a large and diverse protein family, found in yeast and widespread in the animal kingdom. TRP channels work as sensors for a wide range of cellular and environmental signals. Understanding how these channels respond to physical and chemical stimuli has been hindered by the limited structural information available until now. The three-dimensional structure of the vanilloid receptor 1 (TRPV1) was recently determined by single particle electron cryo-microscopy, offering for the first time the opportunity to explore ionic conduction in TRP channels at atomic detail. In this study, we present molecular dynamics simulations of the open-activated pore domain of TRPV1 in the presence of three cationic species: Na(+), Ca(2+) and K(+). The dynamics of these ions while interacting with the channel pore allowed us to rationalize their permeation mechanism in terms of a pathway involving three binding sites at the intracellular cavity, as well as the extracellular and intracellular entrance of the selectivity filter. Furthermore, conformational analysis of the pore in the presence of these ions reveals specific ion-mediated structural changes in the selectivity filter, which influences the permeability properties of the TRPV1 channel.

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.

Local anesthetic and antiepileptic drug access and binding to a bacterial voltage-gated sodium channel

Voltage-gated sodium (Nav) 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 Nav 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 NavAb 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 Nav inhibitors.

Computational Models for Predictive Cardiac Ion Channel Pharmacology

A wealth of experimental data exists describing the elementary building blocks of complex physiological systems. However, it is increasingly apparent in the biomedical sciences that mechanisms of biological function cannot be observed or readily predicted via study of constituent elements alone. This is especially clear in the longstanding failures in prediction of effects of drug treatment for heart rhythm disturbances. These failures stem in part from classical assumptions that have been made in cardiac antiarrhythmic drug development - that a drug operates by one mechanism via one target receptor that arises from one gene.

Side-chain conformation at the selectivity filter shapes the permeation free-energy landscape of an ion channel

On the basis of single-channel currents recorded from the muscle nicotinic acetylcholine receptor (AChR), we have recently hypothesized that the conformation adopted by the glutamate side chains at the first turn of the pore-lining α-helices is a key determinant of the rate of ion permeation. In this paper, we set out to test these ideas within a framework of atomic detail and stereochemical rigor by conducting all-atom molecular dynamics and Brownian dynamics simulations on an extensively validated model of the open-channel muscle AChR. Our simulations provided ample support to the notion that the different rotamers of these glutamates partition into two classes that differ markedly in their ability to catalyze ion conduction, and that the conformations of the four wild-type glutamates are such that two of them "fall" in each rotamer class. Moreover, the simulations allowed us to identify the mm (χ(1) ≅ -60°; χ(2) ≅ -60°) and tp (χ(1) ≅ 180°; χ(2) ≅ +60°) rotamers as the likely conduction-catalyzing conformations of the AChR's selectivity-filter glutamates. More generally, our work shows an example of how experimental benchmarks can guide molecular simulations into providing a type of structural and mechanistic insight that seems otherwise unattainable.

Locating the route of entry and binding sites of benzocaine and phenytoin in a bacterial voltage gated sodium channel

Sodium channel blockers are used to control electrical excitability in cells as a treatment for epileptic seizures and cardiac arrhythmia, and to provide short term control of pain. Development of the next generation of drugs that can selectively target one of the nine types of voltage-gated sodium channel expressed in the body requires a much better understanding of how current channel blockers work. Here we make use of the recently determined crystal structure of the bacterial voltage gated sodium channel NavAb in molecular dynamics simulations to elucidate the position at which the sodium channel blocking drugs benzocaine and phenytoin bind to the protein as well as to understand how these drugs find their way into resting channels. We show that both drugs have two likely binding sites in the pore characterised by nonspecific, hydrophobic interactions: one just above the activation gate, and one at the entrance to the the lateral lipid filled fenestrations. Three independent methods find the same sites and all suggest that binding to the activation gate is slightly more favourable than at the fenestration. Both drugs are found to be able to pass through the fenestrations into the lipid with only small energy barriers, suggesting that this can represent the long posited hydrophobic entrance route for neutral drugs. Our simulations highlight the importance of a number of residues in directing drugs into and through the fenestration, and in forming the drug binding sites.

The treatment of cardiac arrhythmia, epilepsy and pain usually involves blocking the protein channels responsible for initiating electrical activity in nerves and muscles. Current drugs block all such channels, but improved medication requires compounds that can differentiate between the channels present in different parts of the body. Achieving this goal calls for a better understanding of the interactions of current drugs with the proteins. Here we use computer simulation to understand where a local anesthetic and an anti-epileptic bind to a bacterial sodium channel and how they find their way to this position, helping to uncover ways to selectively target different human 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.