Dmitrii G. Luchinsky

Multi-ion conduction bands in a simple model of calcium ion channels

We 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 Qf at the SF. With increasing Qf, there are 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 models 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.

Energetics of discrete selectivity bands and mutation-induced transitions in the calcium-sodium ion channels family

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.

NONEQUILIBRIUM RATE THEORY FOR CONDUCTION IN OPEN ION CHANNELS

We present a nonequilibrium reaction rate model of the ionic transition through an open ion channel, taking account of the interaction between an ion at the entrance of the channel and an ion at the binding site in a self-consistent way. The electrostatic potential is calculated by solution of the Poisson equation for a channel modeled as a cylindrical tube. The transition rate, and the binding site occupancy as a function of the left bulk concentration are compared to 1D Brownian dynamics simulations. The analysis is performed for a single binding site of high-affinity, with the exit rate influenced by barrier fluctuations at the channel exit. The results are compared with experimental data for the permeation of the Na+ ion through the Gramicidin A channel, with which they are shown to be in good agreement.

Resonant multi-ion conduction in a simple model of calcium channels

The ionic permeation of a biological ion channel is a multi-particle, non-equilibrium, stochastic process. Brownian dynamics simulations for a simple electrostatic model of the calcium channel reveal regular structure in the conductance and selectivity as functions of the negative fixed charge Qf on the protein wall at the selectivity filter. This structure consists of distinct high conductance regions (conduction bands) separated by regions of near non-conductance (stop-bands). We report self-consistent electrostatic calculations of single-file, double-ion, stochastic optimal trajectories, and of the energy profiles along these trajectories, for different Qf. We show that the energy difference ΔE along the optimal path exhibits a pronounced minimum near Qf = 3e corresponding to an almost barrier-less (ΔE ~ kBT) resonance-like form of conduction. We demonstrate explicitly that the sharply-defined conduction/selectivity peak of the L-type calcium channel is attributable to the barrier-less knock-on motion of a pair of calcium ions that can occur when their mutual electrostatic repulsion balances their electrostatic attraction to the charge at the selectivity filter. The electrostatics calculations agree well with the results of Brownian dynamics simulations. These results clarify the longstanding puzzle of how the L-type calcium channel exhibits, simultaneously, both high calcium selectivity and conduction at almost the rate of free diffusion.

Self-organized enhancement of conductivity in biological ion channels

We discuss an example of self-organization in a biological system. It arises from long-range ion–ion interactions, and it leads us to propose a new kind of enhanced conduction in ion channels. The underlying mechanism involves charge fluctuations near the channel mouth, amplified by the mismatch between the relative permittivities of water and the protein of the channel walls. We use Brownian dynamics simulations to show that, as in conventional knock on permeation, these interactions can strongly enhance the channel current; but unlike the conventional mechanism, the enhancement occurs without the instigating bath ion entering the channel. The transition between these two mechanisms is clearly demonstrated, emphasizing their distinction. A simple model accurately reproduces the observed phenomena. We point out that electrolyte plus protein of low relative permittivity are universal in living systems, so that long-range ion–ion correlations of the kind considered must be common.

On selectivity and gating of ionic channels

A novel conceptual model is introduced in which ion permeation is coupled to the protein wall vibration and the later in turn modulates exponentially strongly the permeation via radial oscillations of the potential of mean force. In the framework of this model of ion-wall-water interaction we discuss problems of selectivity between alike ions and coupling of ion permeation to gating.

Stochastic dynamics of remote knock-on permeation in biological ion channels

Brownian dynamics simulations provide evidence for a remote knock-on mechanism facilitating the permeation of a biological ion channel by an ion that is initially trapped at the selectivity filter (SF). Unlike the case of conventional direct knock-on, the second ion that instigates permeation does not need to enter the channel. Nor does it necessarily take the place of the permeating ion at the SF, and it can even be of a different ionic species. The study is based on the simultaneous, self-consistent, solution of the coupled Poisson and Langevin equations for a simple generic model, taking account of all the charges present. The new permeation mechanism involves electrostatic amplification attributable to the permittivity mismatch between water and protein: the arrival of the instigating ion at the channel entrance reduces the exit barrier for the ion trapped at the SF, facilitating escape.

Self-consistent analytic solution for the current and access resistance in open ionic channels.

Ionic motion in the bulk solution away from the mouth of a biological ion channel, and inside the channel, is analyzed using Poisson-Nernst-Planck (PNP) equation. The one-dimensional method allows us to connect in a self-consistent way ion dynamics in the bulk solution and inside the channel by taking into account access resistance to the channel. In order to glue the PNP solution in the bulk to that inside the channel, a continuity condition is used for the concentration and the current near the channel mouth at the surface of the hemisphere. The resulting one dimensional (1D) current-voltage characteristics are compared with the Kurnikova16 results which are in good agreement with experimental measurement on the channel, by using a filling factor as the only fitting parameter. The filling factor compensates the fact that the radial charge distribution is non-uniform in a real channel as compared to the cylindrically symmetrical channel used in the 1D approximation.

Polarization switches in vertical-cavity surface-emitting lasers

The polarization dynamics of a vertical cavity surface emitting laser is investigated as a nonlinear stochastic dynamical system. The polarization switches in the device are considered as activation processes in a two dimensional system with a saddle cycle; the optimal way of switching is determined as the solution of a boundary value problem. The theoretical results are in good agreement with the numerical simulations.