Selective ion permeation involves complexation with carboxylates and lysine in a model human sodium channel

Sodium Channel Gene Variants in Fetuses with Abnormal Sonographic Findings: Expanding the Prenatal Phenotypic Spectrum of Sodium Channelopathies

Voltage-gated sodium channels (VGSCs) are responsible for the initiation and propagation of action potentials in the brain and muscle. Pathogenic variants in genes encoding VGSCs have been associated with severe disorders including epileptic encephalopathies and congenital myopathies. In this study, we identified pathogenic variants in genes encoding the α subunit of VGSCs in the fetuses of two unrelated families with the use of trio-based whole exome sequencing, as part of a larger cohort study. Sanger sequencing was performed for variant confirmation as well as parental phasing. The fetus of the first family carried a known de novo heterozygous missense variant in the SCN2A gene (NM_001040143.2:c.751G>A p.(Val251Ile)) and presented intrauterine growth retardation, hand clenching and ventriculomegaly. Neonatally, the proband also exhibited refractory epilepsy, spasms and MRI abnormalities. The fetus of the second family was a compound heterozygote for two parentally inherited novel missense variants in the SCN4A gene (NM_000334.4:c.4340T>C, p.(Phe1447Ser), NM_000334.4:c.3798G>C, p.(Glu1266Asp)) and presented a severe prenatal phenotype including talipes, fetal hypokinesia, hypoplastic lungs, polyhydramnios, ear abnormalities and others. Both probands died soon after birth. In a subsequent pregnancy of the latter family, the fetus was also a compound heterozygote for the same parentally inherited variants. This pregnancy was terminated due to multiple ultrasound abnormalities similar to the first pregnancy. Our results suggest a potentially crucial role of the VGSC gene family in fetal development and early lethality.

Calcium binding and permeation in TRPV channels: Insights from molecular dynamics simulations

Some calcium channels selectively permeate Ca2+, despite the high concentration of monovalent ions in the surrounding environment, which is essential for many physiological processes. Without atomistic and dynamical ion permeation details, the underlying mechanism of Ca2+ selectivity has long been an intensively studied, yet controversial, topic. This study takes advantage of the homologous Ca2+-selective TRPV6 and non-selective TRPV1 and utilizes the recently solved open-state structures and a newly developed multisite calcium model to investigate the ion binding and permeation features in TRPV channels by molecular dynamics simulations. Our results revealed that the open-state TRPV6 and TRPV1 show distinct ion binding patterns in the selectivity filter, which lead to different ion permeation features. Two Ca2+ ions simultaneously bind to the selectivity filter of TRPV6 compared with only one Ca2+ in the case of TRPV1. Multiple Ca2+ binding at the selectivity filter of TRPV6 permeated in a concerted manner, which could efficiently block the permeation of Na+. Cations of various valences differentiate between the binding sites at the entrance of the selectivity filter in TRPV6. Ca2+ preferentially binds to the central site with a higher probability of permeation, repelling Na+ to a peripheral site. Therefore, we believe that ion binding competition at the selectivity filter of calcium channels, including the binding strength and number of binding sites, determines Ca2+ selectivity under physiological conditions.

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac 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.

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel

Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K⁺ and Cl^- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K⁺, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl^- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Druggability of Voltage-Gated Sodium Channels-Exploring Old and New Drug Receptor Sites

Voltage-gated ion channels are important drug targets because they play crucial physiological roles in both excitable and non-excitable cells. About 15% of clinical drugs used for treating human diseases target ion channels. However, most of these drugs do not provide sufficient specificity to a single subtype of the channels and their off-target side effects can be serious and sometimes fatal. Recent advancements in imaging techniques have enabled us for the first time to visualize unique and hidden parts of voltage-gated sodium channels in different structural conformations, and to develop drugs that further target a selected functional state in each channel subtype with the potential for high precision and low toxicity. In this review we describe the druggability of voltage-gated sodium channels in distinct functional states, which could potentially be used to selectively target the channels. We review classical drug receptors in the channels that have recently been structurally characterized by cryo-electron microscopy with natural neurotoxins and clinical drugs. We further examine recent drug discoveries for voltage-gated sodium channels and discuss opportunities to use distinct, state-dependent receptor sites in the voltage sensors as unique drug targets. Finally, we explore potential new receptor sites that are currently unknown for sodium channels but may be valuable for future drug discovery. The advancement presented here will help pave the way for drug development that selectively targets voltage-gated sodium channels.

P-Loop Channels: Experimental Structures, and Physics-Based and Neural Networks-Based Models

The superfamily of P-loop channels includes potassium, sodium, and calcium channels, as well as TRP channels and ionotropic glutamate receptors. A rapidly increasing number of crystal and cryo-EM structures have revealed conserved and variable elements of the channel structures. Intriguing differences are seen in transmembrane helices of channels, which may include π-helical bulges. The bulges reorient residues in the helices and thus strongly affect their intersegment contacts and patterns of ligand-sensing residues. Comparison of the experimental structures suggests that some π-bulges are dynamic: they may appear and disappear upon channel gating and ligand binding. The AlphaFold2 models represent a recent breakthrough in the computational prediction of protein structures. We compared some crystal and cryo-EM structures of P-loop channels with respective AlphaFold2 models. Folding of the regions, which are resolved experimentally, is generally similar to that predicted in the AlphaFold2 models. The models also reproduce some subtle but significant differences between various P-loop channels. However, patterns of π-bulges do not necessarily coincide in the experimental and AlphaFold2 structures. Given the importance of dynamic π-bulges, further studies involving experimental and theoretical approaches are necessary to understand the cause of the discrepancy.

Electrophysiological Properties from Computations at a Single Voltage: Testing Theory with Stochastic Simulations

We use stochastic simulations to investigate the performance of two recently developed methods for calculating the free energy profiles of ion channels and their electrophysiological properties, such as current–voltage dependence and reversal potential, from molecular dynamics simulations at a single applied voltage. These methods require neither knowledge of the diffusivity nor simulations at multiple voltages, which greatly reduces the computational effort required to probe the electrophysiological properties of ion channels. They can be used to determine the free energy profiles from either forward or backward one-sided properties of ions in the channel, such as ion fluxes, density profiles, committor probabilities, or from their two-sided combination. By generating large sets of stochastic trajectories, which are individually designed to mimic the molecular dynamics crossing statistics of models of channels of trichotoxin, p7 from hepatitis C and a bacterial homolog of the pentameric ligand-gated ion channel, GLIC, we find that the free energy profiles obtained from stochastic simulations corresponding to molecular dynamics simulations of even a modest length are burdened with statistical errors of only 0.3 kcal/mol. Even with many crossing events, applying two-sided formulas substantially reduces statistical errors compared to one-sided formulas. With a properly chosen reference voltage, the current–voltage curves can be reproduced with good accuracy from simulations at a single voltage in a range extending for over 200 mV. If possible, the reference voltages should be chosen not simply to drive a large current in one direction, but to observe crossing events in both directions.

Structural Pharmacology of Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_V) channels initiate and propagate action potentials in excitable tissues to mediate key physiological processes including heart contraction and nervous system function. Accordingly, Na_V channels are major targets for drugs, toxins and disease-causing mutations. Recent breakthroughs in cryo-electron microscopy have led to the visualization of human Na_V1.1, Na_V1.2, Na_V1.4, Na_V1.5 and Na_V1.7 channel subtypes at high-resolution. These landmark studies have greatly advanced our structural understanding of channel architecture, ion selectivity, voltage-sensing, electromechanical coupling, fast inactivation, and the molecular basis underlying Na_V channelopathies. Na_V channel structures have also been increasingly determined in complex with toxin and small molecule modulators that target either the pore module or voltage sensor domains. These structural studies have provided new insights into the mechanisms of pharmacological action and opportunities for subtype-selective Na_V channel drug design. This review will highlight the structural pharmacology of human Na_V channels as well as the potential use of engineered and chimeric channels in future drug discovery efforts.

Changes in Ion Selectivity Following the Asymmetrical Addition of Charge to the Selectivity Filter of Bacterial Sodium Channels

Voltage-gated sodium channels (NaVs) play fundamental roles in eukaryotes, but their exceptional size hinders their structural resolution. Bacterial NaVs are simplified homologues of their eukaryotic counterparts, but their use as models of eukaryotic Na⁺ channels is limited by their homotetrameric structure at odds with the asymmetric Selectivity Filter (SF) of eukaryotic NaVs. This work aims at mimicking the SF of eukaryotic NaVs by engineering radial asymmetry into the SF of bacterial channels. This goal was pursued with two approaches: the co-expression of different monomers of the NaChBac bacterial channel to induce the random assembly of heterotetramers, and the concatenation of four bacterial monomers to form a concatemer that can be targeted by site-specific mutagenesis. Patch-clamp measurements and Molecular Dynamics simulations showed that an additional gating charge in the SF leads to a significant increase in Na⁺ and a modest increase in the Ca²⁺ conductance in the NavMs concatemer in agreement with the behavior of the population of random heterotetramers with the highest proportion of channels with charge -5e. We thus showed that charge, despite being important, is not the only determinant of conduction and selectivity, and we created new tools extending the use of bacterial channels as models of eukaryotic counterparts.

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 (P_X/P_Na). The relative permeabilities varied inversely with relative K_d 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.

Covalent linkage of bacterial voltage-gated sodium channels

Background

Bacterial sodium channels are important models for understanding ion permeation and selectivity. However, their homotetrameric structure limits their use as models for understanding the more complex eukaryotic voltage-gated sodium channels (which have a pseudo-heterotetrameric structure formed from an oligomer composed of four domains). To bridge this gap we attempted to synthesise oligomers made from four covalently linked bacterial sodium channel monomers and thus resembling their eukaryotic counterparts.

Results

Western blot analyses revealed NaChBac oligomers to be inherently unstable whereas intact expression of NavMs oligomers was possible. Immunodectection using confocal microscopy and electrophysiological characterisation of NavMs tetramers confirmed plasma membrane localisation and equivalent functionality with wild type NavMs channels when expressed in human embryonic kidney cells.

Conclusion

This study has generated new tools for the investigation of eukaryotic channels. The successful covalent linkage of four bacterial Nav channel monomers should permit the introduction of radial asymmetry into the structure of bacterial Nav channels and enable the known structures of these channels to be used to gain unique insights into structure-function relationships of their eukaryotic counterparts.

Electronic supplementary material

The online version of this article (10.1186/s13628-019-0049-5) contains supplementary material, which is available to authorized users.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

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.