Search

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.

Background. The selectivity of biological cation channels is defined by a short, narrow selectivity filter, having a negative net fixed charge Qf . Voltage gated bacterial channels (NaChBac and some others) are frequently used in biophysics as simplified models of mammalian calcium and sodium channels. We report an experimental, analytic and numerical study of the effects of Qf and bulk ionic concentrations of Ca2+ and Na+ on conduction and selectivity of NaChBac channels, wild type and Qf -varied mutants. Methods. Site-directed mutagenesis and voltage clamp recordings were used to investigate the Na+ /Ca2+ selectivity, divalent blockade and anomalous mole fraction effect (AMFE) for different NaChBac wild type/mutants channels and the properties dependence on Qf . Experimental results were compared with Brownian dynamics simulations and with analytic predictions of the ionic Coulomb blockade (ICB) model, which was extended to encompass bulk concentration effects. Results. It was shown that changing of Qf from –4e (for LESWAS wild type) to –8e (for LEDWAS mutant) leads to strong divalent blockade of the Na+ current by micromolar amounts of Ca2+ ions, similar to the effects seen in mammalian calcium channels. The BD simulations revealed a concentration-related logarithmic shift of the conduction bands. These results were shown to be consistent with ICB model predictions. Conclusions. The extended ICB model explains the experimental (divalent blockade and AMFE) and simulated (multi-ion bands and their concentration-related shifts) selectivity phenomena of NaChBac channel and its charge-varied mutants. These results extend the understanding of ion channel selectivity and may also be applicable to biomimetic nanopores with charged walls.

The highly selective permeation of ions through biological ion channels is an unsolved problem of noise and fluctuations. In this paper, we motivate and introduce a non-equilibrium and self-consistent multi-species kinetic model, with the express aims of comparing with experimental recordings of current versus voltage and concentration and extracting important permeation parameters. For self-consistency, the behavior of the model at the two-state, i.e., selective limit in linear response, must agree with recent results derived from an equilibrium statistical theory. The kinetic model provides a good fit to data, including the key result of an anomalous mole fraction effect.