Biological ion channels are protein nanotubes embedded in, and passing through, the bilipid membranes of cells. Physiologically, they are of crucial importance in that they allow ions to pass into and out of cells, fast and efficiently, though in a highly selective way. Here we show that the conduction and selectivity of calcium/sodium ion channels can be described in terms of ionic Coulomb blockade in a simplified electrostatic and Brownian dynamics model of the channel. The Coulomb blockade phenomenon arises from the discreteness of electrical charge, the strong electrostatic interaction, and an electrostatic exclusion principle. The model predicts a periodic pattern of Ca 2+ conduction versus the fixed charge Q f at the selectivity filter (conduction bands) with a period equal to the ionic charge. It thus provides provisional explanations of some observed and modelled conduction and valence selectivity phenomena, including the anomalous mole fraction effect and the calcium conduction bands. Ionic Coulomb blockade and resonant conduction are similar to electronic Coulomb blockade and resonant tunnelling in quantum dots. The same considerations may also be applicable to other kinds of channel, as well as to charged artificial nanopores.

Online at stacks.iop.org/PhysBio/10/026007 Arxiv at arxiv.org/abs/1209.2381We report self-consistent Brownian dynamics simulations of a simple electrostatic model of the selectivity filters (SF) of calcium ion channels. They reveal regular structure in the conductance and selectivity as functions of the fixed negative charge Q f at the SF. This structure comprises distinct regions of high conductance (conduction bands) M0, M1, M2 separated by regions of almost zero-conductance (stop-bands). Two of these conduction bands, M1 and M2, are related to the saturated calcium occupancies of P=1 and P=2, respectively and demonstrate self-sustained conductivity. Despite the model's limitations, its M1 and M2 bands show high calcium selectivity and prominent anomalous mole fraction effects and can be identified with the L-type and RyR calcium channels. The non-selective band M0 can be identified with a non-selective cation channel, or with OmpF porin.

We use Brownian dynamics (BD) simulations to study the ionic conduction and valence selectivity of a generic electrostatic model of a biological ion channel as functions of the fixed charge Q(f) at its selectivity filter. We are thus able to reconcile the discrete calcium conduction bands recently revealed in our BD simulations, M0 (Q(f)=1e), M1 (3e), M2 (5e), with a set of sodium conduction bands L0 (0.5e), L1 (1.5e), thereby obtaining a completed pattern of conduction and selectivity bands vs Q(f) for the sodium-calcium channels family. An increase of Q(f) leads to an increase of calcium selectivity: L0 (sodium-selective, nonblocking channel) → M0 (nonselective channel) → L1 (sodium-selective channel with divalent block) → M1 (calcium-selective channel exhibiting the anomalous mole fraction effect). We create a consistent identification scheme where the L0 band is putatively identified with the eukaryotic sodium channel The scheme created is able to account for the experimentally observed mutation-induced transformations between nonselective channels, sodium-selective channels, and calcium-selective channels, which we interpret as transitions between different rows of the identification table. By considering the potential energy changes during permeation, we show explicitly that the multi-ion conduction bands of calcium and sodium channels arise as the result of resonant barrierless conduction. The pattern of periodic conduction bands is explained on the basis of sequential neutralization taking account of self-energy, as Q(f)(z,i)=ze(1/2+i), where i is the order of the band and z is the valence of the ion. Our results confirm the crucial influence of electrostatic interactions on conduction and on the Ca(2+)/Na(+) valence selectivity of calcium and sodium ion channels. The model and results could be also applicable to biomimetic nanopores with charged walls.

A self-consistent analytic approach is introduced for the estimation of the access resistance and the current through an open ion channel for an arbitrary number of species. For an ion current flowing radially inward from infinity to the channel mouth, the Poisson-Boltzmann-Nernst-Planck equations are solved analytically in the bulk with spherical symmetry in three dimensions, by linearization. Within the channel, the Poisson-Nernst-Planck equation is solved analytically in a one-dimensional approximation. An iterative procedure is used to match the two solutions together at the channel mouth in a self-consistent way. It is shown that the current-voltage characteristics obtained are in good quantitative agreement with experimental measurements.

Electronic analog experiments on escape over a fluctuating potential barrier are performed for the case when the fluctuations are caused by Ornstein-Uhlenbeck noise (OUN). In its dependence on the relation between the two OUN parameters (the correlation time tau and noise strength Q) the nonmonotonic variation of the mean escape time T as a function of tau can exhibit either a minimum (resonant activation), or a maximum (inhibition of activation), or both these effects. The possible resonant nature of these features is discussed. We claim that T is not a good quantity to describe the resonancelike character of the problem. Independently of the specific relation between the OUN parameters, the resonance manifests itself as a maximal lowering of the potential barrier during the escape event, and it appears for tau of the order of the relaxation time toward the metastable state.