PNP equations with steric effects: a model of ion flow through channels

The flow of current through an ionic channel is studied using the energetic variational approach of Liu applied to the primitive (implicit solvent) model of ionic solutions. This approach allows the derivation of self-consistent (Euler-Lagrange) equations to describe the flow of spheres through channels. The partial differential equations derived involve the global interactions of the spheres and are replaced here with a local approximation that we call steric PNP (Poisson-Nernst-Planck) (Lin, T. C.; Eisenberg, B. To be submitted for publication, 2012). Kong combining rules are used and a range of values of steric interaction parameters are studied. These parameters change the energetics of steric interaction but have no effect on diffusion coefficients in models and simulations. Calculations are made for the calcium (EEEE, EEEA) and sodium channels (DEKA) previously studied in Monte Carlo simulations with comparable results. The biological function is quite sensitive to the steric interaction parameters, and we speculate that a wide range of the function of channels and transporters, even enzymes, might depend on such terms. We point out that classical theories of channels, transporters, and enzymes depend on ideal representations of ionic solutions in which nothing interacts with nothing, even in the enormous concentrations found near and in these proteins or near electrodes in electrochemical cells for that matter. We suggest that a theory designed to handle interactions might be more appropriate. We show that one such theory is feasible and computable: steric PNP allows a direct comparison with experiments measuring flows as well as equilibrium properties. Steric PNP combines atomic and macroscales in a computable formulation that allows the calculation of the macroscopic effects of changes in atomic scale structures (size ~/= 10(-10) meters) studied very extensively in channology and molecular biology.