Improved 3D continuum calculations of ion flux through membrane channels

A continuum model, based on the Poisson–Nernst–Planck (PNP) theory, is applied to simulate steady-state ion flux through protein channels. The PNP equations are modified to explicitly account (1) for the desolvation of mobile ions in the membrane pore and (2) for effects related to ion sizes. The proposed algorithm for a three-dimensional self-consistent solution of PNP equations, in which final results are refined by a focusing technique, is shown to be suitable for arbitrary channel geometry and arbitrary protein charge distribution. The role of the pore shape and protein charge distribution in formation of basic electrodiffusion properties, such as channel conductivity and selectivity, as well as concentration distributions of mobile ions in the pore region, are illustrated by simulations on model channels. The influence of the ionic strength in the bulk solution and of the externally applied electric field on channel properties are also discussed.