Abstract
We present a new size-modified Poisson-Boltzmann ion channel (SMPBIC) model and use it to calculate the electrostatic potential, ionic concentrations, and electrostatic solvation free energy for a voltage-dependent anion channel (VDAC) on a biological membrane in a solution mixture of multiple ionic species. In particular, the new SMPBIC model adopts a membrane surface charge density and a natural Neumann boundary condition to reflect the charge effect of the membrane on the electrostatics of VDAC. To avoid the singularity difficulties caused by the atomic charges of VDAC, the new SMPBIC model is split into three submodels such that the solution of one of the submodels is obtained analytically and contains all the singularity points of the SMPBIC model. The other two submodels are then solved numerically much more efficiently than the original SMPBIC model. As an application of this SMPBIC submodel partitioning scheme, we derive a new formula for computing the electrostatic solvation free energy. Numerical results for a human VDAC isoform 1 (hVDAC1) in three different salt solutions, each with up to five different ionic species, confirm the significant effects of membrane surface charges on both the electrostatics and ionic concentrations. The results also show that the new SMPBIC model can describe well the anion selectivity property of hVDAC1, and that the new electrostatic solvation free energy formula can significantly improve the accuracy of the currently used formula. © 2019 Wiley Periodicals, Inc.