BioMOCA—a Boltzmann transport Monte Carlo model for ion channel simulation

With the recent availability of high-resolution structural information for several key ion channel proteins and large-scale computational resources, Molecular Dynamics has become an increasingly popular tool for ion channel simulation. However, the CPU requirements for simulating ion transport on time scales relevant to conduction still exceed the resources presently available. To address this problem, we have developed Biology Monte Carlo (BioMOCA), a three-dimensional (3D) coarse-grained particle ion channel simulator based on the Boltzmann Transport Monte Carlo (BTMC) methodology. Although this approach is widely employed in the engineering community to study charge transport in electron devices, its application to molecular biology and electrolytes in general is new and hence must be validated. The pair correlation function, which is a measure of the microscopic structure of matter, provides a suitable benchmark to compare the BTMC method against the well-established Equilibrium Monte Carlo (EMC) approach. For validation purposes BioMOCA is used to simulate several simple homogeneous equilibrium electrolytes at concentrations of physiological interest. The ion–ion pair correlation functions computed from these simulations compare very well with those obtained from EMC simulations. We also demonstrate several performance-improving techniques that result in a several-fold speed-up without compromising the pair correlation function. BioMOCA is then used to perform full 3D simulations of ion transport in the gramicidin A channel in situ in a membrane environment, as well as to study the link between the electrostatic and dielectric properties of the protein and the channel's selectivity.

[1]  W. Im,et al.  Ions and counterions in a biological channel: a molecular dynamics simulation of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution. , 2002, Journal of molecular biology.

[2]  M. Kurnikova,et al.  Three-dimensional Poisson-Nernst-Planck theory studies: influence of membrane electrostatics on gramicidin A channel conductance. , 2000, Biophysical journal.

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

[4]  Umberto Ravaioli.  Vectorization of Monte Carlo Algorithms for Semiconductor Simulation , 1991 .

[5]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[6]  Randal R Ketchem,et al.  High-resolution conformation of gramicidin A in a lipid bilayer by solid-state NMR. , 1993, Science.

[7]  Michael J. Holst,et al.  The adaptive multilevel finite element solution of the Poisson-Boltzmann equation on massively parallel computers , 2001, IBM J. Res. Dev..

[8]  S. Hladky,et al.  Ion transfer across lipid membranes in the presence of gramicidin A. I. Studies of the unit conductance channel. , 1972, Biochimica et biophysica acta.

[9]  Robert S. Eisenberg,et al.  Two- and Three-Dimensional Poisson–Nernst–Planck Simulations of Current Flow Through Gramicidin A , 2002, J. Sci. Comput..

[10]  C. D. Cole,et al.  Noncontact dipole effects on channel permeation. I. Experiments with (5F-indole)Trp13 gramicidin A channels. , 1998, Biophysical journal.

[11]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[12]  R. Eisenberg,et al.  Three-Dimensional Continuum Simulations of Ion Transport Through Biological Ion Channels: Effect of Charge Distribution in the Constriction Region of Porin , 2002 .

[13]  P. Phale,et al.  Brownian dynamics simulation of ion flow through porin channels. , 1999, Journal of molecular biology.

[14]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

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

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

[17]  A. Nitzan,et al.  A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel. , 1999, Biophysical journal.

[18]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[19]  S. Selberherr.  Analysis and simulation of semiconductor devices , 1984 .

[20]  G. R. Smith,et al.  The dielectric properties of water within model transbilayer pores. , 1997, Biophysical journal.

[21]  W. Im,et al.  Theoretical and computational models of biological ion channels , 2004, Quarterly Reviews of Biophysics.

[22]  Benoît Roux,et al.  Molecular mechanism of H+ conduction in the single-file water chain of the gramicidin channel. , 2002, Biophysical journal.

[23]  Grace Jordison.  Molecular Biology of the Gene , 1965, The Yale Journal of Biology and Medicine.

[24]  E. Jakobsson,et al.  Brownian dynamics study of a multiply-occupied cation channel: application to understanding permeation in potassium channels. , 1994, Biophysical journal.

[25]  W. L. Jorgensen.  OPLS Force Fields , 2002 .

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

[27]  R. S. Eisenberg,et al.  Computing the Field in Proteins and Channels , 2010, 1009.2857.

[28]  T. Schlick.  Molecular modeling and simulation , 2002 .

[29]  Umberto Ravaioli,et al.  Integration of a particle-particle-particle-mesh algorithm with the ensemble Monte Carlo method for the simulation of ultra-small semiconductor devices , 2000 .

[30]  C. Brooks.  Computer simulation of liquids , 1989 .

[31]  Peter C. Jordan.  Microscopic approaches to ion transport through transmembrane channels: the model system gramicidin , 1987 .

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

[33]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[34]  M. Saraniti,et al.  Computational Issues in Modeling Ion Transport in Biological Channels: Self-Consistent Particle-Based Simulations , 2003 .

[35]  J. Simonin,et al.  Real Ionic Solutions in the Mean Spherical Approximation. 3. Osmotic and Activity Coefficients for Associating Electrolytes in the Primitive Model , 1998 .

[36]  W. Im,et al.  Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics, Brownian dynamics, and continuum electrodiffusion theory. , 2002, Journal of molecular biology.

[37]  W. M. Haynes.  CRC Handbook of Chemistry and Physics , 1990 .

[38]  O. Andersen,et al.  Importance of Tryptophan Dipoles for Protein Function: 5-Fluorination of Tryptophans in Gramicidin A Channels , 1998 .

[39]  W. Im,et al.  A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. , 2000, Biophysical journal.

[40]  Gail Cardew,et al.  Gramicidin and related ion channel-forming peptides , 1999 .

[41]  S. Chung,et al.  Brownian dynamics study of ion transport in the vestibule of membrane channels. , 1998, Biophysical journal.

[42]  BioMOCA: A Transport Monte Carlo Model for Ion Channels , 2003 .

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

[44]  Benoît Roux,et al.  Ion transport in a gramicidin-like channel: dynamics and mobility , 1991 .

[45]  O. Andersen,et al.  Gramicidin channels that have no tryptophan residues. , 1992, Biochemistry.

[46]  A. Warshel,et al.  Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor. , 1985, Journal of molecular biology.

[47]  Serge Durand-Vidal,et al.  Electrolytes at interfaces , 2000 .

[48]  E. Jakobsson,et al.  Stochastic theory of ion movement in channels with single-ion occupancy. Application to sodium permeation of gramicidin channels. , 1987, Biophysical journal.

[49]  Uwe Hollerbach,et al.  Predicting Function from Structure Using the Poisson−Nernst−Planck Equations: Sodium Current in the Gramicidin A Channel , 2000 .

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

[51]  B. Pettitt,et al.  Statistical Thermodynamic Approach to the Chemical Activities in Two-Component Solutions , 2004 .

[52]  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.

[53]  B. Eisenberg,et al.  Ion Channels as Devices , 2003, Bio-, Micro-, and Nanosystems (IEEE Cat. No.03EX733).

[54]  R. Eisenberg,et al.  Simulating Ion Permeation Through the ompF Porin Ion Channel Using Three-Dimensional Drift-Diffusion Theory , 2003 .

[55]  D. Haydon,et al.  Ion transfer across lipid membranes in the presence of gramicidin A. II. The ion selectivity. , 1972, Biochimica et biophysica acta.

[56]  R. Eisenberg,et al.  The effects of the antibiotics gramicidina, amphotericin B, and nystatin on the electrical properties of frog skeletal muscle , 1973 .

[57]  S. Chung,et al.  Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics. , 2000, Biophysical journal.

[58]  F. Reif,et al.  Fundamentals of Statistical and Thermal Physics , 1965 .

[59]  M B Jackson,et al.  Single‐Channel Recording , 1998, Current protocols in neuroscience.

[60]  R. W. Kelsall,et al.  The Monte Carlo method for semiconductor device simulation , 1995 .

[61]  J. Davies,et al.  Molecular Biology of the Cell , 1983, Bristol Medico-Chirurgical Journal.

[62]  Robert S. Eisenberg,et al.  Ion flow through narrow membrane channels: part II , 1992 .

[63]  P. Wolynes,et al.  The theory of ion transport through membrane channels. , 1985, Progress in biophysics and molecular biology.

[64]  O. Andersen,et al.  Molecular determinants of channel function. , 1992, Physiological reviews.

[65]  Shin-Ho Chung,et al.  Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations. , 2002, Biophysical journal.

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

[67]  G. Rummel,et al.  Crystal structures explain functional properties of two E. coli porins , 1992, Nature.

[68]  M. Muir.  Physical Chemistry , 1888, Nature.

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

[70]  Shin-Ho Chung,et al.  Continuum electrostatics fails to describe ion permeation in the gramicidin channel. , 2002, Biophysical journal.

[71]  J. Lundbaek,et al.  Gramicidin and Related Ion channel-Forming Peptides , 2000 .

[72]  R. Eisenberg,et al.  From Structure to Function in Open Ionic Channels , 1999, The Journal of Membrane Biology.

[73]  J. Lear,et al.  Permeation through an open channel: Poisson-Nernst-Planck theory of a synthetic ionic channel. , 1997, Biophysical journal.

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

[75]  A. Warshel,et al.  Electrostatic effects in macromolecules: fundamental concepts and practical modeling. , 1998, Current opinion in structural biology.

[76]  J. Barthel,et al.  Physical Chemistry of Electrolyte Solutions: Modern Aspects , 1998 .

[77]  Ernst Bamberg,et al.  ION TRANSPORT THROUGH THE GRAMICIDIN A CHANNEL , 1976 .

[78]  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.

[79]  C. D. Cole,et al.  Noncontact dipole effects on channel permeation. III. Anomalous proton conductance effects in gramicidin. , 1999, Biophysical journal.

[80]  R. Keynes.  The ionic channels in excitable membranes. , 1975, Ciba Foundation symposium.

[81]  B. Wallace,et al.  Gramicidin channels and pores. , 1990, Annual review of biophysics and biophysical chemistry.

[82]  S. Hladky,et al.  Discreteness of Conductance Change in Bimolecular Lipid Membranes in the Presence of Certain Antibiotics , 1970, Nature.

[83]  B. Roux,et al.  Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data. , 2003, Journal of the American Chemical Society.

[84]  Steven E. Laux.  On particle-mesh coupling in Monte Carlo semiconductor device simulation , 1996, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..