A stabilized finite element method for modified Poisson-Nernst-Planck equations to determine ion flow through a nanopore

Communications in Computational Physics

Volumn 15, Issue 1, 2014, Pages 93-125

  Chaudhry, Jehanzeb Hameed ^a   Comer, Jeffrey ^b   Aksimentiev, Aleksei ^b   Olson, Luke N ^b  

^a Department of Environmental and Radiological Health Sciences   (United States)

^b UNIVERSITY OF ILLINOIS AT URBANA CHAMPAIGN   (United States)

Author keywords

Continuum transport theory; Finite element method; Ion channels; Ionic current; Nucleic acids; PNP; Steric effects; SUPG

Indexed keywords

EID: 84887796020     PISSN: 18152406     EISSN: 19917120     Source Type: Journal    
DOI: 10.4208/cicp.101112.100413a     Document Type: Article

References (66)

1. E. Neher, B. Sakmann, J. Steinbach, The extracellular patch clamp: A method for resolving currents through individual open channels in biological membranes, Pflug. Arch. Eur. J. Physiol. 375 (1978) 219-228. (Pubitemid 8383955)
2. B. Sakmann, E.Neher, Patch clamp techniques for studying ionic channels in excitablemembranes, Ann. Rev. Physiol. 46 (1984) 455-472. (Pubitemid 14147809)
3. J. Mathé, H. Visram, V. Viasnoff, Y. Rabin, A. Meller, Nanopore unzipping of individual DNA hairpin molecules, Biophys. J. 87 (2004) 3205-3212. (Pubitemid 40468577)
4. Q. Zhao, G. Sigalov, V. Dimitrov, B. Dorvel, U.Mirsaidov, S. Sligar, A. Aksimentiev, G. Timp, Detecting SNPs using a synthetic nanopore, Nano Lett. 7 (2007) 1680-1685. (Pubitemid 47140444)
5. Q. Zhao, J. Comer, V. Dimitrov, A. Aksimentiev, G. Timp, Stretching and unzipping nucleic acid hairpins using a synthetic nanopore, Nucl. Acids Res. 36 (2008) 1532-1541. (Pubitemid 351426106)
6. X. F. Kang, S. Cheley, X. Y. Guan, H. Bayley, Stochastic detection of enantiomers, J. Am. Chem. Soc. 128 (2006) 10684-10685. (Pubitemid 44277749)
7. J. J. Kasianowicz, E. Brandin, D. Branton, D. W. Deamer, Characterization of individual polynucleotide molecules using a membrane channel, Proc. Natl. Acad. Sci. USA 93 (1996) 13770-13773. (Pubitemid 26424195)
8. M. Akeson, D. Branton, J. J. Kasianowicz, E. Brandin, D. W. Deamer, Microsecond timescale discrimination among polycytidylic acid, polyadenylic acid, and polyuridylic acid as homopolymers or as segments within singe RNAmolecules, Biophys. J. 77 (1999) 3227-3233.
9. D. Branton, D. Deamer, A. Marziali, H. Bayley, S. Benner, T. Butler, M. Di Ventra, S. Garaj, A. Hibbs, X. Huang, et al., The potential and challenges of nanopore sequencing, Nature Biotech. 26 (10) (2008) 1146-1153.
10. J. Clarke, H. Wu, L. Jayasinghe, A. Patel, S. Reid, H. Bayley, Continuous base identification for single-molecule nanopore DNA sequencing, Nature Nanotech. 4 (2009) 265-270.
11. I. Derrington, T. Butler, M. Collins, E. Manrao, M. Pavlenok, M. Niederweis, J. Gundlach, Nanopore DNA sequencing with MspA, Proc. Natl. Acad. Sci. USA 107 (2010) 16060.
12. W. Im, B. Roux, Ions and counterions in a biological channel: a molecular dynamics study of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution, J. Mol. Biol. 319 (2002) 1177-1197. (Pubitemid 34729427)
13. W. Im, B. Roux, Ion permeation and selectivity of OmpF porin: A theoretical study based on molecular dynamics, brownian dynamics, and continuum electrodiffusion theory, J. Mole. Bio. 322 (2002) 851-869.
14. S. Y. Noskov, W. Im, B. Roux, Ion permeation through the α-hemolysin channel: Theoretical studies based on Brownian Dynamics and Poisson-Nernst-Plank electrodiffusion theory, Biophys. J. 87 (2004) 2299-2309. (Pubitemid 41071331)
15. A. Aksimentiev, J. B. Heng, G. Timp, K. Schulten,Microscopic kinetics of DNA translocation through synthetic nanopores, Biophys. J. 87 (2004) 2086-2097. (Pubitemid 39172203)
16. A. Aksimentiev, K. Schulten, Imaging alpha-hemolysin withmolecular dynamics: Ionic conductance, osmotic permeability and the electrostatic potential map, Biophys. J. 88 (2005) 3745-3761. (Pubitemid 40991090)
17. J. Comer, V. Dimitrov, Q. Zhao, G. Timp, A. Aksimentiev, Microscopic mechanics of hairpin DNA translocation through synthetic nanopores, Biophys. J. 96 (2009) 593-608.
18. S. Pezeshki, C. Chimerel, A.N. Bessonov,M.Winterhalter, U. Kleinekathofer, Understanding ion conductance on a molecular level: An all-atom modeling of the bacterial porin OmpF, Biophys. J. 97 (2009) 1898-1906.
19. Y. Luo, B. Egwolf, D. Walters, B. Roux, Ion selectivity of α-hemolysin with a β-cyclodextrin adapter. I. Single ion potential ofmean force and diffusion coefficient, J. Phys. Chem. B (2009) 2035-2042.
20. A. Aksimentiev, Deciphering ionic current signatures of DNA transport through a nanopore, Nanoscale 2 (2010) 468-483.
21. C.Maffeo, S. Bhattacharya, J. Yoo, D.Wells, A. Aksimentiev,Modeling and simulation of ion channels, Chem. Rev. 112 (2012) 6250-6284.
22. M. P. Allen, D. J. Tildesley, Computer Simulation of Liquids, Oxford University Press, New York, 1987.
23. R. Carr, J. Comer, M. Ginsberg, A. Aksimentiev, Atoms-to-microns model for small solute transport through sticky nanochannels, Lab Chip 11 (2011) 3766-3773.
24. J. Comer, A. Aksimentiev, Predicting the DNA sequence dependence of nanopore ion current using atomic-resolution Brownian dynamics, J. Phys. Chem. C 116 (2012) 3376-3393.
25. M. E. Davis, J. A.McCammon, Electrostatics in biomolecular structure and dynamics, Chem. Rev. 90 (1990) 509-521.
26. V. Barcilon, D.-P. Chen and R. S. Eisenberg, Ion flow through narrow membrane channels. Part II, SIAMJ. Appl. Math. 52 (1992) 1405-1425. (Pubitemid 23589679)
27. M. G. Kurnikova, R. D. Coalson, P. Graf, A. Nitzan, A lattice relaxation algorithm for threedimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel, Biophys J. 76 (1999) 642-656. (Pubitemid 29264548)
28. B. Lu, Y. Zhou, G. A. Huber, S. D. Bond, M. J. Holst, J. A. McCammon, Electrodiffusion: A continuum modeling framework for biomolecular systems with realistic spatiotemporal resolution, J. Chem. Phys. 127 (2007) 135102 (17 pages).
29. M. S. Kilic, M. Z. Bazant, A. Ajdari, Steric effects in the dynamics of electrolytes at large applied voltages. II. Modified Poisson-Nernst-Planck equations, Phys. Rev. E 75 (2007) 021503.
30. D. S. Bolintineanu, A. Sayyed-Ahmad, H. T. Davis, Y. N. Kaznessis, Poisson-Nernst-Planck models of nonequilibrium ion electrodiffusion through a protegrin transmembrane pore, PLoS Comput. Biol. 5 (1) (2009) e1000277.
31. A. E. Cardenas, R. D. Coalson, M. G. Kurnikova, Three-dimensional Poisson-Nernst-Planck theory studies: Influence of membrane electrostatics on gramicidin A channel conductance, Biophys. J. 79 (1) (2000) 80-93. (Pubitemid 30436729)
32. H. Cohen, J. Cooley, The numerical solution of the time-dependentNernst- Planck equations, Biophys. J. 5 (1965) 145-162.
33. B. Lu, Y. C. Zhou, G. A. Huber, S. D. Bond, M. J. Holst, J. A. McCammon, Electrodiffusion: A continuum modeling framework for biomolecular systems with realistic spatiotemporal resolution, J. Chem. Phys. 127 (2007) 135102.
34. B. Lu, M. J. Holst, J. Andrew McCammon, Y. C. Zhou, Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes I: Finite element solutions, J. Comput. Phys. 229 (2010) 6979-6994.
35. B. Lu, Y. Zhou, Poisson-Nernst-Planck equations for simulating biomolecular diffusionreaction processes II: Size effects on ionic distributions and diffusion-reaction rates, Biophys. J. 100 (2011) 2475-2485.
36. Q. Zheng, D. Chen, G.-W.Wei, Second-order Poisson-Nernst-Planck solver for ion transport, J. Comput. Phys. 230 (2011) 5239-5262.
37. A. Logg, G. N. Wells, DOLFIN: Automated finite element computing, ACM Trans. Math. Softw. 37 (2010) 20:1-20:28.
38. D. Braess, Finite Elements: Theory, Fast Solvers, and Applications in Elasticity Theory, 3rd Edition, Cambridge University Press, 2007.
39. P. B. Bochev,M. D. Gunzburger, J.N. Shadid, Stability of the SUPG finite elementmethod for transient advection-diffusion problems, Comp. Meth. Appl. Mech. Engr. 193 (23-26) (2004) 2301-2323. (Pubitemid 38503358)
40. T. J. R. Hughes, Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations, Int. J. Numer. Meth. Fluids 7 (1987) 1261-1275.
41. T. J. R. Hughes, L. P. Franca, G. M. Hulbert, A new finite element formulation for computational fluid dynamics: VIII. The galerkin/least-squares method for advective-diffusive equations, Comp. Meth. Appl. Mech. Engr. 73 (1989) 173-189.
42. L. P. Franca, S. L. Frey, T. J. R. Hughes, Stabilized finite element methods: I. Application to the advective-diffusivemodel, Comp. Meth. Appl. Mech. Engr. 95 (1992) 253-276.
43. T. T. Eu, Stabilized finite element formulations for incompressible flow computations, Adv. Appl. Mech. 28 (1991) 1-44.
44. T. J. Hughes, M.Mallet,M. Akira, A new finite element formulation for computational fluid dynamics: II. Beyond SUPG, Comp. Meth. Appl. Mech. Engr. 54 (1986) 341-355.
45. M. K. Gilson, M. E. Davis, B. A. Luty, J. A. McCammon, Computation of electrostatic forces on solvatedmolecules using the poisson-boltzmann equation, J. Phys. Chem. 97 (1993) 3591-3600.
46. Z. Zhou, P. Payne,M. Vasquez, N. Kuhn,M. Levitt, Finite-difference solution of the Poisson- Boltzmann equation: Complete elimination of self-energy, J. Comput. Chem. 17 (1996) 1344-1351. (Pubitemid 126576751)
47. J. Chaudhry, S. Bond, L. Olson, Finite element approximation to a finite-size modified Poisson-Boltzmann equation, J. Sci. Comput. 47 (2011) 347-364, 10.1007/s10915-010-9441-7.
48. S. D. Bond, J. H. Chaudhry, E. C. Cyr, L. N. Olson, A first-order system least-squares finite element method for the Poisson-Boltzmann equation, J. Comput. Chem. 31 (8) (2010) 1625-1635.
49. I.-L. Chern, J.-G. Liu, W.-C. Wang, Accurate evaluation of electrostatics for macromolecules in solution, Meth. Appl. Anal. 10 (2003) 309-328. (Pubitemid 38825116)
50. B. Aksoylu, S. D. Bond, E. C. Cyr, M. J. Holst, Goal-oriented adaptivity and multilevel preconditioning for the Poisson-Boltzmann equation, J. Sci. Comput. 52 (2012) 202-225.
51. M. Z. Bazant, M. S. Kilic, B. D. Storey, A. Ajdari, Towards an understanding of inducedcharge electrokinetics at large applied voltages in concentrated solutions, Adv. Colloid Interface Sci. 152 (1-2) (2009) 48-88.
52. J. Heng, A. Aksimentiev, C. Ho, P. Marks, Y. Grinkova, S. Sligar, K. Schulten, G. Timp, The electromechanics of dna in a synthetic nanopore, Biophys. J. 90 (2006) 1098-1106.
53. H. R. Drew, R.M.Wing, T. Takano, C. Broka, S. Tanaka, K. Itakura, R. E. Dickerson, Structure of a b-dna dodecamer: Conformation and dynamics, Proc. Natl. Acad. Sci. USA 78 (1981) 2179-2183.
54. URL http://www.rcsb.org/pdb/
55. T. Dolinsky, P. Czodrowski, H. Li, J. Nielsen, J. Jensen, G. K. N. Baker, Pdb2pqr: Expanding and upgrading automated preparation of biomolecular structures formolecular simulations, Nucleic Acids Res. 35 (2007)W522-W525.
56. T. Dolinsky, J. Nielsen, J. McCammon, N. Baker, Pdb2pqr: An automated pipeline for the setup, execution, and analysis of Poisson-Boltzmann electrostatics calculations, Nucleic Acids Res. 35 (2004)W665-W667. (Pubitemid 38997419)
57. Z. Yu,M. J.Holst, Y. Cheng, J. A.McCammon, Feature-preserving adaptivemesh generation for molecular shape modeling and simulation, J. Mol. Graph. Model. 26 (2008) 1370-1380.
58. C. Geuzaine, J.-F. Remacle, Gmsh: A 3-d finite element mesh generator with built-in preand post-processing facilities, Int. J. Numer. Meth. Engr. 79 (2009) 1309-1331.
59. R.M.M. Smeets, U. F. Keyser, D. Krapf,M.-Y.Wu, N. H. Dekker, C. Dekker, Salt dependence of ion transport and dna translocation through solid-state nanopores, Nano Letters 6 (2006) 89-95. (Pubitemid 43166107)
60. M. Young, B. Jayaram, D. Beveridge, Local dielectric environment of B-DNA in solution: Results from a 14 ns molecular dynamics trajectory, J. Phys. Chem. B 102 (1998) 7666-7669. (Pubitemid 128578404)
61. B. Roux, The calculation of the potential of mean force using computer simulations, Comp. Phys. Commun. 91 (1995) 275-282.
62. J. C. Phillips, R. Braun,W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, K. Schulten, Scalable molecular dynamics with NAMD, J. Comp. Chem. 26 (2005) 1781-1802. (Pubitemid 43078511)
63. A. D. MacKerell, Jr., D. Bashford, M. Bellott, R. L. Dunbrack, Jr., J. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, I. W. E. Reiher, B. Roux, M. Schlenkrich, J. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, D. Yin, M. Karplus, Allatom empirical potential for molecular modeling and dynamics studies of proteins, J. Phys. Chem. B 102 (1998) 3586-3616. (Pubitemid 128576688)
64. D. Beglov, B. Roux, Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations, J. Chem. Phys. 100 (1994) 9050-9063.
65. W. Im, B. Roux, Brownian dynamics simulations of ions channels: A general treatment of electrostatic reaction fields for molecular pores of arbitrary geometry, J. Chem. Phys. 115 (2001) 4580.
66. D. B. Wells, V. Abramkina, A. Aksimentiev, Exploring transmembrane transport through α-hemolysin with grid-steered molecular dynamics, J. Chem. Phys. 127 (2007) 125101.