Computing induced charges in inhomogeneous dielectric media: application in a Monte Carlo simulation of complex ionic systems.

The efficient calculation of induced charges in an inhomogeneous dielectric is important in simulations and coarse-grained models in molecular biology, chemical physics, and electrochemistry. We present the induced charge computation (ICC) method for the calculation of the polarization charges based on the variational formulation of Allen et al. [Phys. Chem. Chem. Phys. 3, 4177 (2001)]. We give a different solution for their extremum condition that produces a matrix formulation. The induced charges are directly calculated by solving the linear matrix equation Ah=c, where h contains the discretized induced charge density, c depends only on the source charges-the ions moved in the simulation-and the matrix A depends on the geometry of dielectrics, which is assumed to be unchanged during the simulation. Thus, the matrix need be inverted only once at the beginning of the simulation. We verify the efficiency and accuracy of the method by means of Monte Carlo simulations for two special cases. In the simplest case, a single sharp planar dielectric boundary is present, which allows comparison with exact results calculated using the method of electrostatic images. The other special case is a particularly simple case where the matrix A is not diagonal: a slab with two parallel flat boundaries. Our results for electrolyte solutions in these special cases show that the ICC method is both accurate and efficient.

[1]  Shin-Ho Chung,et al.  Study of ionic currents across a model membrane channel using Brownian dynamics. , 1998, Biophysical journal.

[2]  D. Busath,et al.  Monte Carlo Study of the Effect of Ion and Channel Size on the Selectivity of a Model Calcium Channel , 2001 .

[3]  Shin-Ho Chung,et al.  Tests of continuum theories as models of ion channels. II. Poisson-Nernst-Planck theory versus brownian dynamics. , 2000, Biophysical journal.

[4]  Minoru Sakurai,et al.  Medium effects on the molecular electronic structure. I. The formulation of a theory for the estimation of a molecular electronic structure surrounded by an anisotropic medium , 1987 .

[5]  Rudolph A. Marcus,et al.  On the Theory of Oxidation‐Reduction Reactions Involving Electron Transfer. I , 1956 .

[6]  D. Levitt,et al.  Electrostatic calculations for an ion channel. I. Energy and potential profiles and interactions between ions. , 1978, Biophysical journal.

[7]  Dirk Gillespie,et al.  Ion Accumulation in a Biological Calcium Channel: Effects of Solvent and Confining Pressure , 2001 .

[8]  W. Im,et al.  Brownian dynamics simulations of ions channels: A general treatment of electrostatic reaction fields for molecular pores of arbitrary geometry , 2001 .

[9]  Uwe Hollerbach,et al.  Dielectric boundary force and its crucial role in gramicidin. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Peter C. Jordan Electrostatic modeling of ion pores. Energy barriers and electric field profiles. , 1982, Biophysical journal.

[11]  B. Nadler,et al.  Derivation of Poisson and Nernst-Planck equations in a bath and channel from a molecular model. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  D. Busath,et al.  Monte Carlo study of the selectivity of calcium channels: improved geometrical model , 2002 .

[13]  Daniel Borgis,et al.  Electrostatics on particles: Phenomenological and orientational density functional theory approach , 2002 .

[14]  J. Hansen,et al.  Nonlinear counterion screening in colloidal suspensions , 1993 .

[15]  S. Chung,et al.  Mechanisms of permeation and selectivity in calcium channels. , 2001, Biophysical journal.

[16]  G. Torrie,et al.  Electrical double layers. II. Monte Carlo and HNC studies of image effects , 1982 .

[17]  Shin-Ho Chung,et al.  Solutions of Poisson's equation in channel-like geometries , 1998 .

[18]  Shin-Ho Chung,et al.  COMPUTER SIMULATION OF ION CONDUCTANCE IN MEMBRANE CHANNELS , 1998 .

[19]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[20]  B. Eisenberg,et al.  Binding and selectivity in L-type calcium channels: a mean spherical approximation. , 2000, Biophysical journal.

[21]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[22]  D. Busath,et al.  Monte Carlo Simulations of the Mechanism for Channel Selectivity: The Competition between Volume Exclusion and Charge Neutrality , 2000 .

[23]  J. Tomasi,et al.  Description of the solvent effects for large molecules: a linear scaling procedure , 2001 .

[24]  P. Läuger,et al.  Nonlinear electrical effects in lipid bilayer membranes. II. Integration of the generalized Nernst-Planck equations. , 1969, Biophysical journal.

[25]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[26]  Peter C. Jordan Electrostatic modeling of ion pores. II. Effects attributable to the membrane dipole potential. , 1983, Biophysical journal.

[27]  B. U. Felderhof Fluctuations of polarization and magnetization in dielectric and magnetic media , 1977 .

[28]  Rudolph A. Marcus,et al.  Electrostatic Free Energy and Other Properties of States Having Nonequilibrium Polarization. I , 1956 .

[29]  S. Kuyucak,et al.  Role of the dielectric constants of membrane proteins and channel water in ion permeation. , 2003, Biophysical journal.

[30]  Intermittent permeation of cylindrical nanopores by water. , 2002, Physical review letters.

[31]  J. Hansen,et al.  Molecular dynamics investigation of water permeation through nanopores , 2003 .

[32]  Christian Silvio Pomelli,et al.  An improved iterative solution to solve the electrostatic problem in the polarizable continuum model , 2001 .

[33]  Shin-Ho Chung,et al.  Permeation of ions across the potassium channel: Brownian dynamics studies. , 1999, Biophysical journal.

[34]  C. Kelley Iterative Methods for Linear and Nonlinear Equations , 1987 .

[35]  R. Fox,et al.  Classical Electrodynamics, 3rd ed. , 1999 .

[36]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[37]  D. Levitt Electrostatic calculations for an ion channel. II. Kinetic behavior of the gramicidin A channel. , 1978, Biophysical journal.

[38]  Bob Eisenberg,et al.  Monte Carlo simulations of ion selectivity in a biological Na channel: Charge–space competition , 2002 .

[39]  Richard J. Sadus,et al.  Molecular Simulation of Fluids: Theory, Algorithms and Object-Orientation , 1999 .

[40]  G. R. Smith,et al.  Simulation approaches to ion channel structure–function relationships , 2001, Quarterly Reviews of Biophysics.

[41]  Jess P. Wilcoxon,et al.  Photooxidation of Organic Chemicals Catalyzed by Nanoscale MoS2 , 1999 .

[42]  P. Lugli,et al.  The Monte Carlo Method for Semiconductor Device Simulation , 1990 .

[43]  P. Attard Variational formulation for the electrostatic potential in dielectric continua , 2003 .

[44]  Dirk Gillespie,et al.  Monte Carlo Simulation Study of a System with a Dielectric Boundary: Application to Calcium Channel Selectivity , 2004 .

[45]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[46]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[47]  M. Marchi,et al.  A dielectric continuum molecular dynamics method , 2001 .

[48]  Peter C. Jordan Energy barriers for passage of ions through channels. Exact solution of two electrostatic problems. , 1981, Biophysical chemistry.

[49]  A. Parsegian,et al.  Energy of an Ion crossing a Low Dielectric Membrane: Solutions to Four Relevant Electrostatic Problems , 1969, Nature.

[50]  Simone Melchionna,et al.  Electrostatic potential inside ionic solutions confined by dielectrics: a variational approach , 2001 .

[51]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.