A dynamic lattice Monte Carlo model of ion transport in inhomogeneous dielectric environments: Method and implementation

Journal of Physical Chemistry B

Volumn 104, Issue 51, 2000, Pages 12324-12338

  Graf, Peter ^a   Nitzan, Abraham ^a   Kurnikova, Maria G ^b   Coalson, Rob D ^b  

^a TEL AVIV UNIVERSITY   (Israel)

^b UNIVERSITY OF PITTSBURGH   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

BOUNDARY CONDITIONS; COMPUTER SIMULATION; CRYSTAL LATTICES; ELECTROSTATICS; IONS; MATHEMATICAL MODELS; MONTE CARLO METHODS; PROTEINS; WATER;

DYNAMIC LATTICE MONTE CARLO (DLMC) SIMULATION;

DIELECTRIC MATERIALS;

EID: 0034516442     PISSN: 15206106     EISSN: None     Source Type: Journal    
DOI: 10.1021/jp001282s     Document Type: Article

References (67)

1. Schmickler, W. Interfacial Electrochemistry; Oxford University Press: New York, 1996.
2. Bockris, J. O'M.; Khan, S. U. M. Surface Electrochemistry. A molecular level approach; Plenum Press: New York, 1993.
3. Gray, F. M. Solid Polymer Electrolytes. Fundamentals and technological applications; VCH: Weinheim, Germany, 1991.
4. Hille, B. Ionic Channels of Excitable Membranes, 2nd ed.; Sinauer Associates, Inc.: Sunderland, MA, 1992.
5. Levitt, D. G. J. Gen. Physiol. 1999, 113, 789.
6. Nicholls, A.; Sharp, K. A.; Honig, B. DelPhi V3.0; Columbia University, New York, 1990.
7. Honig, B.; Sharp, K.; Yana, A. S. J. Pliys. Client. 1993, 97, 1101.
8. Honig, B.; Nicholls, A. Science 1995, 268, 1144.
9. Eisenberg, R. S. J. Membr. Biol. 1996, 750, 1.
10. Eisenberg, R. S. Ace. Client. Res. 1998, 31, 117.
11. Eisenberg, R. S. J. Membr. Biol. 1999, 171, 1.
12. Nonner, W.; Chen, D. P.; Eisenberg, R. S. J. Gen. Physiol. 1999, 113, 773.
13. Syganow, A.; von Kitzing, E. J. Phys. Chem. 1995, 99, 12030.
14. Syganow, A.; von Kitzing, E. Biophys. J. 1999, 76, 768.
15. Syganow, A.; von Kitzing, E. Eur. Biophys. J. 1999, 28, 393.
16. Kurnikova, M. G.; Coalson, R. D.; Graf, P.; Nitzan, A. Biophys. J. 1999, 76, 642.
17. Cdrdenas, A. E.; Coalson R. D.; Kurnikova, M. G. Biophys. J. 2000, 79, 80.
18. Hollerbach, U.; Chen, D. P.; Busath, D. D.; Eisenberg, B. Langmuir 2000, 16, 5509.
19. See, e.g., Caillol, J. M. J. Chem. Phys. 1999, 111, 6528, and references therein.
20. Chung, S.-Ho; Hoyles, M.; Allen, T.; Kuyucak, S. Biophys. J. 1998, 75, 793.
21. Hoyles, M.; Kuyucak, S.; Chung, S.-Ho Phys. Rev. E 1998, 55, 3654.
22. Corry, B.; Kuyucak, S.; Chung, S.-Ho J. Gen. Physiol. 1999,114, 597.
23. Chung, S. H.; Allen, T. W.; Hoyles, M.; Kuyucak, S. Biophys. J. 1999, 77, 2517.
24. Allen, T. W.; Hoyles, M.; Kuyucak, S.; Chung, S. H. Chem. Phys. Lett. 1999, 313, 358.
25. Three relevant papers have appeared after our paper was submitted, (a) Moy, G.; Corry, B.; Kuyucak, S.; Chung, S.-Ho Biophys. J. 2000, 78, 2349-2363,2364-2381, describes Brownian dynamics simulations of ion transport through a pore with comparison to mean field (PB and PNP) calculations. Our methodology is different from the approach taken in these papers. Nevertheless our results are consistent with the observations made in these works.
26. (b) Im, W.; Seefeld, S.; Roux, B. Biophys. J. 2000, 79, 788, describes a Brownian dynamics simulation of ion motion in a Porin channel. In this paper dielectric effects on self-energy and interparticle interactions of mobile ions are disregarded.
27. Elber, R.; Chen, D. P.; Rojewska, D.; Eisenberg, R. Biophys. J. 1995, 68, 906.
28. Woolf, T. B. and Roux, B. Biophys. J. 1997, 72, 1930.
29. Tieleman, D. P.; Berendsen, H. J. C. Biophys. J. 1998, 74, 2786.
30. Knödler, D.; Pendzig, P.; Dieterich, W. 1994, 70, 356.
31. Knödler, D.; Dieterich, W.; Lonsky, C.; Nitzan, A. J. Client. Phys. 1995, 702, 465.
32. Pendzig, P.; Dieterich, W.; Nitzan, A. J. Noncryst. Solids. 1998, 235, 748.
33. For an impressive demonstration of such a program in the study of fracture dynamics in silicon see: Abraham, F. F.; Broughton, J. Q.; Bernstein, N.; Kaxiras, E. Comput. Phys. 1998, 12, 538.
34. A description of an equilibrium ionic channel in which the system is separated into an inner and an outer subsystems was recently presented by: Roux, B. Biophys. J. 1999, 77, 139.
35. See, e.g.: Barcilon, V.; Chen, D.; Eisenberg, R. S.; Ratner, M. A. J. Chem. Phys. 1993, 98, 1193.
36. Böttcher, C. J. F. Theory of Electric Polarization. Dielectrics in Static Fields, 2nd ed.; Eisevier: New York, 1973.
37. Shaw, P. B. Phys. Rev. A 1985, 32, 2476.
38. Newton, M. D.; Friedman, H. L. J. Chem. Phys. 1988, 88, 4460.
39. In the application of these periodic boundary conditions in the present work the corresponding interion interactions were computed with the minimum image approximation.
40. For example, composite grids could be used: one grid for the calculation of the electrostatic potential and a different, finer grid to represent ion configuration and ionic motion where the ion occupies more than a single grid cell. For the use of composite grids in electrostatic calculations see: Beck, T. L. J. Comput. Chem. 1999, 20, 1731-1739.
41. Binder, K.; Heermann, D. W. Monte Carlo simulation in statistical physics; Springer-Verlag: New York, 1988.
42. See, e.g.: Binder, K. Monte Carlo Investigations of Phase Transitions and Critical Phenomena. Phase Transitions and Critical Phenomena; Academic Press: New York, 1976; Vol. 5b, Chapter 1, pp 1-105.
43. 2/(2a))(1/ε - 1).
44. Using Coulomb potentials with a cutoff is a potential source of error. Alternatively, using Ewald summation in inhomogeneous restricted geometry systems (in our case a 2-dimensional Ewald summation is required) is computationally very demanding. Moreover, the inhomogeneous net charge distribution of mobile ions in our system would lead to unphysical interactions between periodic extensions of the system in the Ewald method. There are indications that the errors introduced by using Coulomb interactions with cutoff are not very large for many applications, see, e.g.: Schweighofer, K.; Benjamin, I. J. Chem. Phys. 2000,112, 1474-1482 and references therein. Controlling these errors is important in situations where system properties are dominated by the balance between opposing interactions, e.g., Coulombic repulsion vs dielectric attraction in solvated ions of the same charge sign
45. (Bader, J. S.; Chandler, D. J. Phys. Chem. 1992, 96, 6423).
46. c.
47. (0) reflects in the usual way the periodic boundary condition in the x- and y-direction and the finite system size in the z-direction.
48. 6
49. The unphysical image of polarization charges, induced by individual ions on the internal dielectric interfaces, on the imposed boundaries is not eliminated by this procedure. If, however, the imposed boundaries are far from the dielectric interfaces the resulting unphysical contribution to the total energy is small.
50. 2 = |r - (0, d)| are the distances between point r and the charges q and q′, respectively. (Note that we compute here only the polarization contribution to the potential Φ. These expressions, which diverge on the interface, become valid far enough (from distances of order a, or our grid spacing) from the interface.
51. See: Khurkats, Y. I.; Ullstrup, J. J. Electroanal. Chem. 1991, 308, 17-26 for a discussion of this point.
52. As in section 6b, a part that is not eliminated comes from the image of the polarization charges induced at the internal dielectric boundaries; however, this contribution can be made small (see footnote 47).
53. B = 80 is used in all calculations reported here, because in anticipated applications ions reside in a water environment.
54. An alternative procedure would be to calculate the contribution of the outer medium to the solvation energy of an ion near the boundary using, e.g., a Poisson-Boltzmann description of the outer medium (see, e.g.: Sadus, R. J. Molecular Simulations of Fluids; Elsevier: Amsterdam, 1999). This would, however, raise the issue of validity of the PB approximation for computing solvation energies near a system's boundary.
55. See for example: McQuarry, D. A. Statistical Mechanics: Harper & Row: New York, 1976.
56. Mafé, S.; Pellicer, J.; Aguilella, V. M. J. Phys. Chem. 1986, 90, 6045.
57. Aguilella, V. M.; Mafé, S.; Pellicer, J. Electrochim. Ada 1987, 32, 483.
58. Riveros, O. J.; Croxton, T. L.; Armstrong, W. McD. J. Theor. Biol. 1989, 140, 221.
59. In fact, the correspondence between the two approaches is not complete because the PNP theory disregards the retarding effect of the ionic solvation atmosphere on an ion's motion. This retardation affects the conductivity of a dilute homogeneous electrolyte solution by a few percent.
60. R = 0.096 M that correspond to the bulk concentration obtained in the equilibrium DLMC simulation. See Figure 5 and ensuing discussion.
61. Dieckmann, G. R.; Lear, J. D.; Zhong, Q.; Klein, M. L.; DeGrado, W. F.; Sharp, K. A. Biophys. J. 1999, 76, 618.
62. In realistic biological channels, which may contain internal structural barriers, blocking of ionic motion as a result of ions being trapped for a relatively long time inside the channel may actually lead to a sublinear I/ΔΦ characteristic. See, e.g.: Andersen, O. S. Biophys. J. 1983, 41, 135.
63. Sham, Y. Y.; Muegge, I.; Warshel, A. Biophys. J. 1998, 74, 1744-1753.
64. (b) Muegge, I.; Schweins, T.; Warshel, A. Proteins: Struct. Fund. Genet. 1998, 30, 407-423.
65. Simonson, T.; Brooks, C. L. J. Am. Chem. Soc. 1996, 118, 8452-8458.
66. See, e.g.: Press: W. H.; Teukolsky, S. A.; Vetterling, W. T.; Flannery, B. P. P. Numerical Recipes in FORTRAN: The An of Scientific Computing, 2nd ed.; Cambridge University Press: London, 1992.
67. Wesseling, P. An Introduction to Multigrid Methods; John Wiley & Sons Ltd.: Chichester, UK, 1992.