On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation

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.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Application of a Statistical and Linear Response Theory to Multi-Ion Na+ Conduction in NaChBac

Biological ion channels are fundamental to maintaining life. In this manuscript we apply our recently developed statistical and linear response theory to investigate Na+ conduction through the prokaryotic Na+ channel NaChBac. This work is extended theoretically by the derivation of ionic conductivity and current in an electrochemical gradient, thus enabling us to compare to a range of whole-cell data sets performed on this channel. Furthermore, we also compare the magnitudes of the currents and populations at each binding site to previously published single-channel recordings and molecular dynamics simulations respectively. In doing so, we find excellent agreement between theory and data, with predicted energy barriers at each of the four binding sites of ∼4,2.9,3.6, and 4kT.

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.

Ionic Coulomb blockade and the determinants of selectivity in the NaChBac bacterial sodium channel

Mutation-induced transformations of conductivity and selectivity in NaChBac bacterial channels are studied experimentally and interpreted within the framework of ionic Coulomb blockade (ICB), while also taking account of resonant quantised dehydration (QD) and site protonation. Site-directed mutagenesis and whole-cell patch-clamp experiments are used to investigate how the fixed charge Q_f at the selectivity filter (SF) affects both valence selectivity and same-charge selectivity. The new ICB/QD model predicts that increasing ∣Q_f∣ should lead to a shift in selectivity sequences toward larger ion sizes, in agreement with the present experiments and with earlier work. Comparison of the model with experimental data leads to the introduction of an effective charge Q_f^∗ at the SF, which was found to differ between Aspartate and Glutamate charged rings, and also to depend on position within the SF. It is suggested that protonation of the residues within the restricted space of the SF is important in significantly reducing the effective charge of the EEEE ring. Values of Q_f^∗ derived from experiments on divalent blockade agree well with expectations based on the ICB/QD model and have led to the first demonstration of ICB oscillations in Ca²⁺ conduction as a function of the fixed charge. Preliminary studies of the dependence of Ca²⁺ conduction on pH are qualitatively consistent with the predictions of the model.

Different roles for aspartates and glutamates for cation permeation in bacterial sodium channels

A key driving force for ion channel selectivity is represented by the negative charge of the Selectivity Filter carried by aspartate (D) and glutamate (E) residues. However, the structural effects and specific properties of D and E residues have not been extensively studied. In order to investigate this issue we studied the mutants of NaChBac channel with all possible combinations of D and E in the charged rings in position 191 and 192. Electrophysiological measurements showed significant Ca²⁺ currents only when position 191 was occupied by E. Equilibrium Molecular Dynamics simulations revealed the existence of two binding sites, corresponding to the charged rings and another one, more internal, at the level of L190. The simulations showed that the ion in the innermost site can interact with the residue in position 191 only when this is glutamate. Based on the MD simulations, we suggest that a D in position 191 leads to a high affinity Ca²⁺ block site resulting from a significant drop in the free energy of binding for an ion moving between the binding sites; in contrast, the free energy change is more gradual when an E residue occupies position 191, resulting in Ca²⁺ permeability. This scenario is consistent with the model of ion channel selectivity through stepwise changes in binding affinity proposed by Dang and McCleskey. Our study also highlights the importance of the structure of the selectivity filter which should contribute to the development of more detailed physical models for ion channel selectivity.