Second-order Poisson-Nernst-Planck solver for ion transport

Journal of Computational Physics

Volumn 230, Issue 13, 2011, Pages 5239-5262

  Zheng, Qiong ^a   Chen, Duan ^a   Wei, Guo Wei ^a,b  

^a MICHIGAN STATE UNIVERSITY   (United States)

^b Michigan State University   (United States)

Author keywords

Dirichlet to Neumann mapping; Gramicidin A; Ion channels; Matched interface and boundary; Poisson Nernst Planck equations

Indexed keywords

COMPUTATION THEORY; CONTINUUM MECHANICS; GEOMETRY; IONS; ITERATIVE METHODS; NONLINEAR EQUATIONS; PROTEINS;

DIRICHLET TO NEUMANN MAPPING; DIRICHLET-TO-NEUMANN; GEOMETRIC SINGULARITY; GRAMICIDIN A; ION CHANNEL; MATCHED INTERFACE AND BOUNDARY; NUMERICAL ALGORITHMS; POISSON-NERNST-PLANCK; POISSON-NERNST-PLANCK EQUATIONS; SECOND ORDERS;

MAPPING;

EID: 79955476670     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2011.03.020     Document Type: Article

References (81)

1. Abaid N., Eisenberg R.S., Liu W.S. Asymptotic expansions of I-V relations via a Poisson-Nernst-Planck system. SIAM J. Appl. Dyn. Syst. 2008, 7(4):1507-1526.
2. Allen R., Hansen J.-P., Melchionna S. Electrostatic potential inside ionic solutions confined by dielectrics: a variational approach. Phys. Chem. Chem. Phys. 2001, 3:4177-4186.
3. Babuška I. The finite element method for elliptic equations with discontinuous coefficients. Computing 1970, 5:207-213.
4. Barcilon V. Ion flow through narrow membran channels: part I. SIAM J. Appl. Math. 1992, 52(5):1391-1404.
5. Barcilon V., Chen D., Eisenberg R. Ion flow through narrow membrane channels: Part ii. SIAM J. Appl. Math. 1992, 52:1405-1425.
6. Bardhan J.P., Eisenberg R.S., Gillespie D. Discretization of the induced-charge boundary integral equation. Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 2009, 80(1):11906-11910.
7. Bates P.W., Chen Z., Sun Y.H., Wei G.W., Zhao S. Geometric and potential driving formation and evolution of biomolecular surfaces. J. Math. Biol. 2009, 59:193-231.
8. Bates P.W., Wei G.W., Zhao S. Minimal molecular surfaces and their applications. J. Comput. Chem. 2008, 29(3):380-391.
9. Bazant M.Z., Thornton K., Ajdari A. Diffuse-charge dynamics in electrochemical systems. Phys. Rev. E 2004, 70:021506.
10. Boda D., Gillespie D., Nonner W., Henderson D., Eisenberg B. Computing induced charges in inhomogeneous dielectric media: application in a Monte Carlo simulation of complex ionic systems. Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 2004, 69:046702.
11. Boda D., Nonner W., Henderson D., Eisenberg B., Gillespie D. Volume exclusion in calcium selective channels. Biophys. J. 2008, 94:3486-3496.
12. Boda D., Nonner W., Valisko M., Henderson D., Eisenberg B., Gillespie D. Steric selectivity in Na channels arising from protein polarization and mobile side chains. Biophys. J. 2007, 93(6):1960-1980.
13. Boda D., Valisko M., Eisenberg B., Nonner W., Henderson D., Gillespie D. The combined effect of pore radius and protein dielectric coefficient on the selectivity of a calcium channel. Phys. Rev. Lett. 2007, 98:168102.
14. Boda D., Varga T., Henderson D., Busath D., Nonner W., Gillespie D., Eisenberg B. Monte carlo simulation study of a system with a dielectric boundary: application to calcium channel selectivity. Mol. Simul. 2004, 30:89-96.
15. Bramble J., King J. A finite element method for interface problems in domains with smooth boundaries and interfaces. Adv. Comput. Math. 1996, 6:109-138.
16. Cardenas A.E., Coalson R.D., Kurnikova M.G. Three-dimensional Poisson-Nernst-Planck theory studies: influence of membrane electrostatics on Gramicidin A channel conductance. Biophys. J. 2000, 79:80-93.
17. Chen D., Chen Z., Chen C., Geng W.H., Wei G.W. MIBPB: a software package for electrostatic analysis. J. Comput. Chem. 2011, 32:657-670.
18. Chen T., Strain J. Piecewise-polynomial discretization and Krylov-accelerated multigrid for elliptic interface problems. J. Comput. Phys. 2008, 16:7503-7542.
19. Chen Z., Baker N.A., Wei G.W. Differential geometry based solvation models I: Eulerian formulation. J. Comput. Phys. 2010, 229:8231-8258.
20. Cheng M.H., Coalson R.D. An accurate and efficient empirical approach for calculating the dielectric self-energy and ion-ion pair potential in continuum models of biological ion channels. J. Phys. Chem. B 2005, 109(1):488-498.
21. Cheng M.H., Coalson R.D., Tang P. Molecular dynamics and Brownian dynamics investigation of ion permeation and anesthetic halothane effects on a proton-gated ion channel. J. Am. Chem. Soc. 2010, 132(46):16442-16449.
22. Chern I.L., Liu J.-G., Weng W.-C. Accurate evaluation of electrostatics for macromolecules in solution. Methods Appl. Anal. 2003, 10(2):309-328.
23. Choudhary O.P., Ujwal R., Kowallis W., Coalson R., Abramson J., Grabe M. The electrostatics of VDAC: implications for selectivity and gating. J. Mol. Biol. 2010, 396(3):580-592.
24. Chung S.-H., Kuyucak S. Recent advances in ion channel research. Biochim. Biophys. Acta 2002, 1565:267-286.
25. Coalson R.D., Kurnikova M.G. Poisson-Nernst-Planck theory approach to the calculation of current through biological ion channels. IEEE Trans. Nanobiosci. 2005, 4(1):81-93.
26. Cole C.D., Frost A.S., Thompson N., Cotten M., Cross T.A., Busath D.D. Noncontact dipole effects on channel permeation. vi. 5f- and 6F-Trp gramicidin channel currents. Biophys. J. 2002, 83:1974-1986.
27. Dolinsky T.J., Czodrowski P., Li H., Nielsen J.E., Jensen J.H., Klebe G., Baker N.A. PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucl. Acids Res. 2007, 35(S2):W522-W525.
28. Duncan A., Sedgewick R.D., Coalson R.D. Improved local lattice approach for Coulombic simulations. Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 2005, 71:046702.
29. Eisenberg B., Hyon Y., Liu C. Energy variational analysis of ions in water and channels: field theory for primitive models of complex ionic fluids. J. Chem. Phys. 2010, 133(104104).
30. Eisenberg B., Liu W.S. Poisson-Nernst-Planck systems for ion channels with permanent charges. SIAM J. Appl. Dyn. Syst. 2006, 38(6):1932-1966.
31. Eisenberg R., Chen D. Poisson-Nernst-Planck (PNP) theory of an open ionic channel. Biophys. J. 1993, 64:A22.
32. Eisenberg R.S. Computing the field in proteins and channels. J. Membrane Biol. 1996, 150:1-25.
33. Fedkiw R.P., Aslam T., Merriman B., Osher S. A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method). J. Comput. Phys. 1999, 152:457-492.
34. Gabashvili I.S., Sokolowski B.H.A., Morton C.C., Giersch A.B.S. Ion channel gene expression in the inner ear. J. Assoc. Res. Otolaryngol. 2007, 8(3):305-328.
35. Geng W., Yu S., Wei G.W. Treatment of charge singularities in implicit solvent models. J. Chem. Phys. 2007, 127:114106.
36. Gillespie D. Energetics of divalent selectivity in a calcium channel: the ryanodine receptor case study. Biophys. J. 2008, 94(4):1169-1184.
37. Gillespie D., Valisko M., Boda D. Density functional theory of the electrical double layer: the RFD functional. J. Phys.: Condens. Matter 2005, 17:6609-6626.
38. Gummel H.K. A self-consistent iterative scheme for the one-dimensional steady-state transistor calculations. IEEE Trans. Electron Dev. 1964, 11:455-465.
39. Hacker M., Messer W., Bachmann K. Pharmacology: Principles and Practice 2009, Academic Press.
40. Hollerbach U., Chen D.P., Eisenberg R.S. Two- and three-dimensional Poisson-Nernst-Planck simulations of current flow through Gramicidin A. J. Sci. Comput. 2002, 16(4):373-409.
41. M. Holst, The Poisson-Boltzmann equation: analysis and multilevel numerical solution, PhD Thesis, California Institute of Technology, 1994.
42. Hwang H., Schatz G.C., Ratner M.A. Incorporation of inhomogeneous ion diffusion coefficients into kinetic lattice grand canonical Monte Carlo simulations and application to ion current calculations in a simple model ion channel. J. Phys. Chem. A 2007, 111(49):12506-12512.
43. Im W., Roux B. Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics, Brownian dynamics, and continuum electrodiffusion theory. J. Mol. Biol. 2002, 322:851-869.
44. Jacoboni C., Lugli P. The Monte Carlo Method for Semiconductor Device Simulation 1989, Springer-Verlag, New York.
45. Jerome J. Analysis of Charge Transport. Mathematical Theory and Approximation of Semiconductor Models 1995, Springer-Verlag, New York.
46. Jung Y.W., Lu B.Z., Mascagni M. A computational study of ion conductance in the KcsA K+ channel using a Nernst-Planck model with explicit resident ions. J. Chem. Phys. 2009, 131(215101).
47. Kurnikova M.G., Coalson R.D., Graf P., Nitzan A. A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the Gramicidin A channel. Biophys. J. 1999, 76:642-656.
48. Kuyucak S., Andersen O.S., Chung S.-H. Models of permeation in ion channels. Rep. Prog. Phys. 2001, 64:1427-1472.
49. LeVeque R.J., Li Z.L. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. Anal. 1994, 31:1019-1044.
50. Levitt D.G. Interpretation of biological ion channel flux data-reaction-rate versus continuum theory. Annu. Rev. Biophys. Biophys. Chem. 1986, 15:29-57.
51. Levitt D.G. Modeling of ion channels. J. Gen. Physiol. 1999, 113(6):789-794.
52. Li Z.L., Ito K. Maximum principle preserving schemes for interface problems with discontinuous coefficients. SIAM J. Sci. Comput. 2001, 23:339-361.
53. Liu W.K., Liu Y., Farrell D., Zhang L., Wang X., Fukui Y., Patankar N., Zhang Y., Bajaj C., Chen X., Hsu H. Immersed finite element method and its applications to biological systems. Comput. Methods Appl. Mech. Eng. 2006, 195:1722-1749.
54. Liu W.S. Geometric singular perturbation approach to steady-state Poisson-Nernst-Planck systems. SIAM J. Appl. Math. 2005, 65(3):754-766.
55. Liu X.D., Fedkiw R.P., Kang M. A boundary condition capturing method for Poisson's equation on irregular domains. J. Comput. Phys. 2000, 160:151-178.
56. MacKerell A.D., Bashford D., Bellot M., Dunbrack R.L., Evanseck J.D., Field M.J., Fischer S., Gao J., Guo H., Ha S., Joseph-McCarthy D., Kuchnir L., Kuczera K., Lau F.T.K., Mattos C., Michnick S., Ngo T., Nguyen D.T., Prodhom B., Reiher W.E., Roux B., Schlenkrich M., Smith J.C., Stote R., Straub J., Watanabe M., Wiórkiewicz-Kuczera J., Yin D., Karplus M. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 1998, 102(18):3586-3616.
57. Mamonov A.B., Kurnikova M.G., Coalson R.D. Diffusion constant of K+ inside Gramicidin A: a comparative study of four computational methods. Biophys. Chem. 2006, 124:268-278.
58. Marx D., Hutter J. Ab initio molecular dynamics: theory and implementation. NIC Series 2000, vol. 3:329-477. John von Neumann Institute for Computing, Jülich. J. Grotendorst (Ed.).
59. Mathur S.R., Murthy J.Y. A multigrid method for the Poisson-Nernst-Planck equations. SIAM J. Appl. Math. 2009, 52(17-18):4031-4039.
60. Mayo A. The fast solution of Poisson's and the biharmonic equations on irregular regions. SIAM J. Numer. Anal. 1984, 21:285-299.
61. Oevermann M., Klein R. A cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces. J. Comput. Phys. 2006, 219:749-769.
62. Peskin C.S. Numerical analysis of blood flow in the heart. J. Comput. Phys. 1977, 25(3):220-252.
63. Purves D., Augustine G.J., Fitzpatrick D., Katz L.C., LaMantia A.S., McNamara J.O., Williams S.M. Neuroscience 2001, Sinauer Associates, Sunderland, MA. second ed.
64. Richards F.M. Areas, volumes, packing, and protein structure. Annu. Rev. Biophys. Bioeng. 1977, 6(1):151-176.
65. Roux B., Allen T., Berneche S., Im W. Theoretical and computational models of biological ionchannels. Q. Rev. Biophys. 2004, 7(1):1-103.
66. Rubinstein I. Electro-diffusion of Ions 1990, SIAM, Philadelphia.
67. Sanner M.F., Olson A.J., Spehner J.C. Reduced surface: an efficient way to compute molecular surfaces. Biopolymers 1996, 38:305-320.
68. Schuss Z., Nadler B., Eisenberg R. Derivation of Poisson and Nernst-Planck equations in a bath and channel from a molecular model. Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 2001, 64:036116.
69. Swanson J.M.J., Mongan J., McCammon J.A. Limitations of atom-centered dielectric functions in implicit solvent models. J. Phys. Chem. B 2005, 109(31):14769-14772.
70. Tyagi S., Suzen M., Sega M., Barbosa M., Kantorovich S.S., Holm C. An iterative, fast, linear-scaling method for computing induced charges on arbitrary dielectric boundaries. J. Chem. Phys. 2010, 132(15):154112-154119.
71. Wei G.W. Differential geometry based multiscale models. Bull. Math. Biol. 2010, 72:1562-1622.
72. Xie D., Zhou S.Z. A new minimization protocol for solving nonlinear Poisson-Boltzmann mortar finite element equation. BIT Numer. Math. 2007, 47:853-871.
73. Yu S., Geng W., Wei G.W. Treatment of geometric singularities in implicit solvent models. J. Chem. Phys. 2007, 126:244108.
74. Yu S., Wei G.W. Three-dimensional matched interface and boundary (MIB) method for treating geometric singularities. J. Comput. Phys. 2007, 227:602-632.
75. Yu S., Zhou Y., Wei G.W. Matched interface and boundary (MIB) method for elliptic problems with sharp-edged interfaces. J. Comput. Phys. 2007, 224(2):729-756.
76. Zhao S., Wei G.W. High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces. J. Comput. Phys. 2004, 200(1):60-103.
77. Zhao S., Wei G.W. Matched interface and boundary (MIB) for the implementation of boundary conditions in high order central finite differences. Int. J. Numer. Methods Eng. 2009, 77:1690-1730.
78. Zhou Y.C., Feig M., Wei G.W. Highly accurate biomolecular electrostatics in continuum dielectric environments. J. Comput. Chem. 2008, 29:87-97.
79. Zhou Y.C., Wei G.W. On the fictitious-domain and interpolation formulations of the matched interface and boundary (MIB) method. J. Comput. Phys. 2006, 219(1):228-246.
80. Zhou Y.C., Zhao S., Feig M., Wei G.W. High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources. J. Comput. Phys. 2006, 213(1):1-30.
81. Zhou Z., Payne P., Vasquez M., Kuhn N., Levitt M. Finite-difference solution of the Poisson-Boltzmann equation: complete elimination of self-energy. J. Comput. Chem. 1996, 17:1344-1351.