The drift-diffusion (Poisson-Nernst-Planck) model is applied to the potassium channel in a biological membrane plus surrounding solution baths. Two-dimensional cylindrically symmetric simulations of the K channel in KCl solutions are presented which show significant boundary layers at the ends of the channel and display the spreading of charge into the bath regions. The computed current-voltage...

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....

We test the validity of the mean-field approximation in Poisson-Nernst-Planck theory by contrasting its predictions with those of Brownian dynamics simulations in schematic cylindrical channels and in a realistic potassium channel. Equivalence of the two theories in bulk situations is demonstrated in a control study. In simple cylindrical channels, considerable differences are found between the...

The concept of electrodiffusion based on the Nernst-Planck equations for ionic fluxes coupled with the Poisson equation expressing relation between gradient of the electric field and the charge density is widely used in many areas of natural sciences and engineering. In contrast to the steady-state solutions of the Nernst-Planck-Poisson (abbreviated as NPP or PNP) equations, little is known abo...

Electrical conductance of an aqueous pore in the lipid bilayer has an important role in the process of membrane electroporation, i.e., formation of pores induced by electric pulses. In our present study we compare the pore conductance as predicted by a theoretical model based on the continuum Poisson-Nernst-Planck theory to the pore conductance obtained with molecular dynamics simulations (Casc...

A novel algorithmic scheme for numerical solution of the 3D Poisson-Nernst-Planck model is proposed. The algorithmic improvements are universal and independent of the detailed physical model. They include three major steps: an adjustable gradient-based step value, an adjustable relaxation coefficient, and an optimized segmentation of the modeled space. The enhanced algorithm significantly accel...

Modeling of biomolecular systems plays an essential role in understanding biological processes, such as ionic flow across channels, protein modification or interaction, and cell signaling. The continuum model described by the Poisson-Boltzmann (PB)/Poisson-Nernst-Planck (PNP) equations has made great contributions towards simulation of these processes. However, the model has shortcomings in its...

Incorporation of van der Waals interactions enables the continuum model of electrodiffusion in biomolecular system to avoid the artifacts of introducing a molecular surface and the painful task of the surface mesh generation. Calculation examples show that the electrostatics, diffusion-reaction kinetics, and molecular surface defined as an isosurface of a certain density distribution can be ext...