Advanced Search
Article Contents
Article Contents

Ion size effects on individual fluxes via Poisson-Nernst-Planck systems with Bikerman's local hard-sphere potential: Analysis without electroneutrality boundary conditions

Abstract Full Text(HTML) Figure(6) Related Papers Cited by
  • A quasi-one-dimensional steady-state Poisson-Nernst-Planck model with Bikerman's local hard-sphere potential for ionic flows of two oppositely charged ion species through a membrane channel is analyzed. Of particular interest is the qualitative properties of ionic flows in terms of individual fluxes without the assumption of electroneutrality conditions, which is more realistic to study ionic flow properties of interest. This is the novelty of this work. Our result shows that ⅰ) boundary concentrations and relative size of ion species play critical roles in characterizing ion size effects on individual fluxes; ⅱ) the first order approximation $\mathcal{J}_{k1} = D_kJ_{k1}$ in ion volume of individual fluxes $\mathcal{ J}_k = D_kJ_k$ is linear in boundary potential, furthermore, the signs of $\partial_V \mathcal{ J}_{k1}$ and $\partial^2_{Vλ} \mathcal{J}_{k1}$ , which play key roles in characterizing ion size effects on ionic flows can be both negative depending further on boundary concentrations while they are always positive and independent of boundary concentrations under electroneutrality conditions (see Corollaries 3.2-3.3, Theorems 3.4-3.5 and Proposition 3.7). Numerical simulations are performed to identify some critical potentials defined in (2). We believe our results will provide useful insights for numerical and even experimental studies of ionic flows through membrane channels.

    Mathematics Subject Classification: Primary: 34A26, 34B16, 34D15; Secondary: 37D10, 92C35.


    \begin{equation} \\ \end{equation}
  • 加载中
  • Figure 1.  Numerical detection of critical values for $\lambda$ (left graph) and $m$ (right one) in Theorem 3.2.

    Figure 2.  Numerical detection of critical values for $\sigma$, which corresponds to statement (Ⅰ) in Theorem 3.2. The left graph is for $\lambda<\lambda_1^* = 0.072$, and the right one is for $\lambda>\lambda_2^* = 13.93.$

    Figure 3.  Numerical detection of critical values $\sigma$ corresponding to statement (Ⅱ) in Theorem 3.2 with $\lambda_1^*<\lambda<\lambda_2^*$. The left graph is for $0<m<m^* = 0.4934$, and the right one is for $m^*<m<\frac{1}{2}.$

    Figure 4.  Numerical identification of six critical potentials in (15) with $z_1 = -z_2 = 1$. In the left column, the vertical axis actually represents, from top to bottom, ${\mathcal I}(V;\nu, \lambda)-{\mathcal I}_0(V), \ {\mathcal J}_1(V;\nu, \lambda)-{\mathcal J}_{10}(V)$ and ${\mathcal J}_2(V;\nu, \lambda)-{\mathcal J}_{20}(V)$, respectively. In particular, the x-axis for all figures actually represents $\frac{e}{k_BT}V$.

    Figure 5.  Numerical identification of critical potentials $V_{1c}$ and $ V_1^c$ for individual flux ${\mathcal J}_1$ with $z_1 = -z_2 = 1$ and nonzero permanent charge. The x-axis for all figures actually represents $\frac{e}{k_BT}V$.

    Figure 6.  Numerical approximations of critical potentials $V_1^c$ for individual flux ${\mathcal J}_1$ with $z_1 = -z_2 = 1$ and nonzero permanent charge as illustrated in Proposition 3.14. The x-axis for all figures actually represents $\frac{e}{k_BT}V$.

  • [1] N. AbaidR. S. Eisenberg and W. Liu, Asymptotic expansions of Ⅰ-Ⅴ relations via a Poisson-Nernst-Planck system, SIAM J. Appl. Dyn. Syst., 7 (2008), 1507-1526.  doi: 10.1137/070691322.
    [2] V. Barcilon, Ion flow through narrow membrane channels: Part Ⅰ, SIAM J. Appl. Math., 52 (1992), 1391-1404.  doi: 10.1137/0152080.
    [3] D. BodaD. BusathB. EisenbergD. Henderson and W. Nonner, Monte Carlo simulations of ion selectivity in a biological Na+ channel: Charge-space competition, Phys. Chem. Chem. Phys., 4 (2002), 5154-5160.  doi: 10.1039/B203686J.
    [4] V. BarcilonD.-P. Chen and R. S. Eisenberg, Ion flow through narrow membrane channels: Part Ⅱ, SIAM J. Appl. Math., 52 (1992), 1405-1425.  doi: 10.1137/0152081.
    [5] V. BarcilonD.-P. ChenR. S. Eisenberg and J. W. Jerome, Qualitative properties of steady-state Poisson-Nernst-Planck systems: Perturbation and simulation study, SIAM J. Appl. Math., 57 (1997), 631-648.  doi: 10.1137/S0036139995312149.
    [6] M. BurgerR. S. Eisenberg and H. W. Engl, Inverse problems related to ion channel selectivity, SIAM J. Appl. Math., 67 (2007), 960-989.  doi: 10.1137/060664689.
    [7] D. BodaD. GillespieW. NonnerD. Henderson and B. Eisenberg, Computing induced charges in inhomogeneous dielectric media: application in a Monte Carlo simulation of complex ionic systems, Phys. Rev. E, 69 (2004), 046702 (1-10).  doi: 10.1103/PhysRevE.69.046702.
    [8] J. J. Bikerman, Structure and capacity of the electrical double layer, Philos. Mag., 33 (1942), 384-397.  doi: 10.1080/14786444208520813.
    [9] P. W. BatesY. JiaG. LinH. Lu and M. Zhang, Individual flux study via steady-state Poisson-Nernst-Planck systems: Effects from boundary conditions, SIAM J. Appl. Dyn. Syst., 16 (2017), 410-430.  doi: 10.1137/16M1071523.
    [10] P. W. BatesW. LiuH. Lu and M. Zhang, Ion size and valence effects on ionic flows via Poisson-Nernst-Planck models, Commu. Math. Sci., 15 (2017), 881-901.  doi: 10.4310/CMS.2017.v15.n4.a1.
    [11] J. H. ChaudhryS. D. Bond and L. N. Olson, Finite Element Approximation to a Finite-Size Modified Poisson-Boltzmann Equation, J. Sci. Comput, 47 (2011), 347-364.  doi: 10.1007/s10915-010-9441-7.
    [12] A. E. CardenasR. D. Coalson and M. G. Kurnikova, Three-dimensional poisson-nernst-planck theory studies: Influence of membrane electrostatics on gramicidin a channel conductance, Biophys. J., 79 (2000), 80-93.  doi: 10.1016/S0006-3495(00)76275-8.
    [13] D. P. Chen and R. S. Eisenberg, Charges, currents and potentials in ionic channels of one conformation, Biophys. J., 64 (1993), 1405-1421.  doi: 10.1016/S0006-3495(93)81507-8.
    [14] S. Chung and S. Kuyucak, Predicting channel function from channel structure using Brownian dynamics simulations, Clin. Exp. Pharmacol Physiol., 28 (2001), 89-94. 
    [15] R. Coalson and M. Kurnikova, Poisson-Nernst-Planck theory approach to the calculation of current through biological ion channels, IEEE Transaction on NanoBioscience, 4 (2005), 81-93.  doi: 10.1109/TNB.2004.842495.
    [16] B. Eisenberg, Ion channels as devices, J. Comp. Electro., 2 (2003), 245-249.  doi: 10.1023/B:JCEL.0000011432.03832.22.
    [17] B. Eisenberg, Proteins, channels, and crowded ions, Biophys. Chem., 100 (2003), 507-517.  doi: 10.1016/S0301-4622(02)00302-22.
    [18] R. S. Eisenberg, Channels as enzymes, J. Memb. Biol., 115 (1990), 1-12.  doi: 10.1007/BF01869101.
    [19] R. S. Eisenberg and A. Biology, Electrostatics and Ionic Channels, In New Developments and Theoretical Studies of Proteins, R. Elber, Editor, 269-357, World Scientific, Philadelphia, (1996). 
    [20] R. S. Eisenberg, From structure to function in open ionic channels, J. Memb. Biol., 171 (1999), 1-24.  doi: 10.1007/s002329900554.
    [21] B. Eisenberg and W. Liu, Poisson-Nernst-Planck systems for ion channels with permanent charges, SIAM J. Math. Anal., 38 (2007), 1932-1966.  doi: 10.1137/060657480.
    [22] B. EisenbergW. Liu and H. Xu, Reversal permanent charge and reversal potential: Case studies via classical Poisson-Nernst-Planck models, Nonlinearity, 28 (2015), 103-127. 
    [23] A. ErnR. Joubaud and T. Leliévre, Mathematical study of non-ideal electrostatic correlations in equilibrium electrolytes, Nonlinearity, 25 (2012), 1635-1652.  doi: 10.1088/0951-7715/25/6/1635.
    [24] J. Fischer and U. Heinbuch, Relationship between free energy density functional, Born-Green-Yvon, and potential distribution approaches for inhomogeneous fluids, J. Chem. Phys., 88 (1988), 1909-1913.  doi: 10.1063/1.454114.
    [25] D. Gillespie and R. S. Eisenberg, Physical descriptions of experimental selectivity measurements in ion channels, European Biophys. J., 31 (2002), 454-466.  doi: 10.1007/s00249-002-0239-x.
    [26] D. Gillespie, A Singular Perturbation Analysis of the Poisson-Nernst-Planck System: Applications to Ionic Channels, Ph. D Dissertation, Rush University at Chicago, 1999.
    [27] D. GillespieW. Nonner and R. S. Eisenberg, Coupling Poisson-Nernst-Planck and density functional theory to calculate ion flux, J. Phys.: Condens. Matter, 14 (2002), 12129-12145.  doi: 10.1088/0953-8984/14/46/317.
    [28] D. Gillespie, W. Nonner and R. S. Eisenberg, Density functional theory of charged, hard-sphere fluids Phys. Rev. E, 68 (2003), 0313503 (1-10). doi: 10.1103/PhysRevE.68.031503.
    [29] D. GillespieW. Nonner and R. S. Eisenberg, Crowded charge in biological ion channels, Nanotech, 3 (2003), 435-438. 
    [30] E. GongadzeU. van RienenV. Kralj-lglič and A. lglič, Spatial variation of permittivity of an electrolyte solution in contact with a charged metal surface: A mini review, Computer Method in Biomechanics and Biomedical Engineering, 16 (2013), 463-480. 
    [31] D. GillespieL. XuY. Wang and G. Meissner, (De)constructing the Ryanodine Receptor: Modeling Ion Permeation and Selectivity of the Calcium Release Channel, J. Phys. Chem. B, 109 (2005), 15598-15610.  doi: 10.1021/jp052471j.
    [32] U. HollerbachD.-P. Chen and R. S. Eisenberg, Two-and three-dimensional Poisson-Nernst-Planck simulations of current flow through gramicidin-A, J. Comp. Science, 16 (2002), 373-409. 
    [33] U. Hollerbach, D. Chen, W. Nonner and B. Eisenberg, Three-dimensional Poisson-Nernst-Planck Theory of Open Channels, Biophys. J. , 76 (1999), p. A205.
    [34] Y. HyonB. Eisenberg and C. Liu, A mathematical model for the hard sphere repulsion in ionic solutions, Commun. Math. Sci., 9 (2011), 459-475.  doi: 10.4310/CMS.2011.v9.n2.a5.
    [35] Y. HyonJ. FonsecaB. Eisenberg and C. Liu, Energy variational approach to study charge inversion (layering) near charged walls, Discrete Contin. Dyn. Syst. Ser. B, 17 (2012), 2725-2743.  doi: 10.3934/dcdsb.2012.17.2725.
    [36] Y. HyonC. Liu and B. Eisenberg, PNP equations with steric effects: A model of ion flow through channels, J. Phys. Chem. B, 116 (2012), 11422-11441. 
    [37] W. ImD. Beglov and B. Roux, Continuum solvation model: Electrostatic forces from numerical solutions to the Poisson-Bolztmann equation, Comp. Phys. Comm., 111 (1998), 59-75. 
    [38] W. Im and B. Roux, Ion permeation and selectivity of OmpF porin: A theoretical study based on molecular dynamics, Brownian dynamics, and continuum electrodiffusion theory, J. Mol. Biol., 322 (2002), 851-869.  doi: 10.1016/S0022-2836(02)00778-7.
    [39] S. Ji and W. Liu, Poisson-Nernst-Planck systems for ion flow with density functional theory for hard-sphere potential: Ⅰ-Ⅴ relations and critical potentials. Part Ⅰ: Analysis, J. Dyn. Diff. Equat., 24 (2012), 955-983.  doi: 10.1007/s10884-012-9277-y.
    [40] S. JiW. Liu and M. Zhang, Effects of (small) permanent charges and channel geometry on ionic flows via classical Poisson-Nernst-Planck models, SIAM J. on Appl. Math., 75 (2015), 114-135.  doi: 10.1137/140992527.
    [41] Y. JiaW. Liu and M. Zhang, Qualitative properties of ionic flows via Poisson-Nernst-Planck systems with Bikerman's local hard-sphere potential: Ion size effects, Discrete and Continuous Dynamical Systems-B, 21 (2016), 1775-1802.  doi: 10.3934/dcdsb.2016022.
    [42] M. S. Kilic, M. Z. Bazant and A. Ajdari, Steric effects in the dynamics of electrolytes at large applied voltages. Ⅱ. Modified Poisson-Nernst-Planck equations Phys. Rev. E, 75 (2007), 021503 (11 pages). doi: 10.1103/PhysRevE.75.021503.
    [43] A. S. Khair and T. M. Squires, Ion steric effects on electrophoresis of a colloidal particle, J. Fluid Mech., 640 (2009), 343-356.  doi: 10.1017/S0022112009991728.
    [44] C.-C. LeeH. LeeY. HyonT.-C. Lin and C. Liu, New Poisson-Boltzmann type equations: One-dimensional solutions, Nonlinearity, 24 (2011), 431-458.  doi: 10.1088/0951-7715/24/2/004.
    [45] B. Li, Continuum electrostatics for ionic solutions with non-uniform ionic sizes, Nonlinearity, 22 (2009), 811-833.  doi: 10.1088/0951-7715/22/4/007.
    [46] W. Liu, Geometric singular perturbation approach to steady-state Poisson-Nernst-Planck systems, SIAM J. Appl. Math., 65 (2005), 754-766.  doi: 10.1137/S0036139903420931.
    [47] W. Liu, One-dimensional steady-state Poisson-Nernst-Planck systems for ion channels with multiple ion species, J. Differential Equations, 246 (2009), 428-451.  doi: 10.1016/j.jde.2008.09.010.
    [48] G. LinW. LiuY. Yi and M. Zhang, Poisson-Nernst-Planck systems for ion flow with density functional theory for local hard-sphere potential, SIAM J. Appl. Dyn. Syst., 12 (2013), 1613-1648.  doi: 10.1137/120904056.
    [49] W. LiuX. Tu and M. Zhang, Poisson-Nernst-Planck Systems for Ion Flow with Density Functional Theory for Hard-Sphere Potential: Ⅰ-Ⅴ relations and Critical Potentials. Part Ⅱ: Numerics, J. Dyn. Diff. Equat., 24 (2012), 985-1004.  doi: 10.1007/s10884-012-9278-x.
    [50] W. Liu and B. Wang, Poisson-Nernst-Planck systems for narrow tubular-like membrane channels, J. Dyn. Diff. Equat., 22 (2010), 413-437.  doi: 10.1007/s10884-010-9186-x.
    [51] W. Liu and H. Xu, A complete analysis of a classical Poisson-Nernst-Planck model for ionic flow, J. Differential Equations, 258 (2015), 1192-1228. 
    [52] W. Nonner and R. S. Eisenberg, Ion permeation and glutamate residues linked by Poisson-Nernst-Planck theory in L-type Calcium channels, Biophys. J., 75 (1998), 1287-1305.  doi: 10.1016/S0006-3495(98)74048-2.
    [53] S. Y. NoskovW. Im and B. Roux, Ion Permeation through the $z_1$-Hemolysin Channel: Theoretical Studies Based on Brownian Dynamics and Poisson-Nernst-Planck Electrodiffusion Theory, Biophys. J., 87 (2004), 2299-2309. 
    [54] S. Y. Noskov and B. Roux, Ion selectivity in potassium channels, Biophys. Chem., 124 (2006), 279-291.  doi: 10.1016/j.bpc.2006.05.033.
    [55] B. NadlerZ. SchussA. Singer and B. Eisenberg, Diffusion through protein channels: From molecular description to continuum equations, Nanotech., 3 (2003), 439-442. 
    [56] J. K. Percus, Equilibrium state of a classical fluid of hard rods in an external field, J. Stat. Phys., 15 (1976), 505-511.  doi: 10.1007/BF01020803.
    [57] J. K. Percus, Model grand potential for a nonuniform classical fluid, J. Chem. Phys., 75 (1981), 1316-1319.  doi: 10.1063/1.442136.
    [58] J.-K. Park and J. W. Jerome, Qualitative properties of steady-state Poisson-Nernst-Planck systems: Mathematical study, SIAM J. Appl. Math., 57 (1997), 609-630.  doi: 10.1137/S0036139995279809.
    [59] B. RouxT. W. AllenS. Berneche and W. Im, Theoretical and computational models of biological ion channels, Quat. Rev. Biophys., 37 (2004), 15-103.  doi: 10.1017/S0033583504003968.
    [60] Y. Rosenfeld, Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing, Phys. Rev. Lett., 63 (1989), 980-983.  doi: 10.1103/PhysRevLett.63.980.
    [61] Y. Rosenfeld, Free energy model for the inhomogeneous fluid mixtures: Yukawa-charged hard spheres, general interactions, and plasmas, J. Chem. Phys., 98 (1993), 8126-8148.  doi: 10.1063/1.464569.
    [62] R. Roth, Fundamental measure theory for hard-sphere mixtures: A review, J. Phys.: Condens. Matter, 22 (2010), 063102 (1-18).  doi: 10.1088/0953-8984/22/6/063102.
    [63] A. SingerD. GillespieJ. Norbury and R. S. Eisenberg, Singular perturbation analysis of the steady-state Poisson-Nernst-Planck system: applications to ion channels, European J. Appl. Math., 19 (2008), 541-560. 
    [64] A. Singer and J. Norbury, A Poisson-Nernst-Planck model for biological ion channels-an asymptotic analysis in a three-dimensional narrow funnel, SIAM J. Appl. Math., 70 (2009), 949-968.  doi: 10.1137/070687037.
    [65] Z. Schuss, B. Nadler and R. S. Eisenberg, Derivation of Poisson and Nernst-Planck equations in a bath and channel from a molecular model Phys. Rev. E, 64 (2001), 036116. doi: 10.1103/PhysRevE.64.036116.
    [66] X.-S. WangD. HeJ. Wylie and H. Huang, Singular perturbation solutions of steady-state Poisson-Nernst-Planck systems, Phys. Rev. E, 89 (2014), 022722 (1-14).  doi: 10.1103/PhysRevE.89.022722.
    [67] G. W. WeiQ. ZhengZ. Chen and K. Xia, Variational multiscale models for charge transport, SIAM Review, 54 (2012), 699-754.  doi: 10.1137/110845690.
    [68] J. ZhangD. Acheampong and M. Zhang, Geometric singular approach to Poisson-Nernst-Planck models with excess chemical potentials: Ion size effects on individual fluxes, Mol. Based Math. Biol., 5 (2017), 58-77. 
    [69] M. Zhang, Asymptotic expansions and numerical simulations of Ⅰ-Ⅴ relations via a steady-state Poisson-Nernst-Planck system, Rocky Mountain J. Math., 45 (2015), 1681-1708.  doi: 10.1216/RMJ-2015-45-5-1681.
    [70] Q. Zheng and W. Wei, Poisson-Boltzmann-Nernst-Planck model, J. Chem. Phys., 134 (2011), 194101 (1-17).  doi: 10.1063/1.3581031.
  • 加载中



Article Metrics

HTML views(1651) PDF downloads(291) Cited by(0)

Access History



    DownLoad:  Full-Size Img  PowerPoint