Computational electrophysiology: the molecular dynamics of ion channel permeation and selectivity in atomistic detail.

Presently, most simulations of ion channel function rely upon nonatomistic Brownian dynamics calculations, indirect interpretation of energy maps, or application of external electric fields. We present a computational method to directly simulate ion flux through membrane channels based on biologically realistic electrochemical gradients. In close analogy to single-channel electrophysiology, physiologically and experimentally relevant timescales are achieved. We apply our method to the bacterial channel PorB from pathogenic Neisseria meningitidis, which, during Neisserial infection, inserts into the mitochondrial membrane of target cells and elicits apoptosis by dissipating the membrane potential. We show that our method accurately predicts ion conductance and selectivity and elucidates ion conduction mechanisms in great detail. Handles for over- coming channel-related antibiotic resistance are identified.

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

[2]  Alan E Mark,et al.  Simulation of pore formation in lipid bilayers by mechanical stress and electric fields. , 2003, Journal of the American Chemical Society.

[3]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[4]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[5]  H. Berendsen,et al.  Molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine bilayer with different macroscopic boundary conditions and parameters , 1996 .

[6]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[7]  B. Roux.  The membrane potential and its representation by a constant electric field in computer simulations. , 2008, Biophysical journal.

[8]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[9]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[10]  Martin Neumann,et al.  Dipole moment fluctuation formulas in computer simulations of polar systems , 1983 .

[11]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[12]  Gerhard Hummer,et al.  Osmotic water transport through carbon nanotube membranes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Sanjay Ram,et al.  The role of porins in neisserial pathogenesis and immunity. , 2003, Trends in microbiology.

[14]  S. Hallermann,et al.  Short openings in high resolution single channel recordings of mouse nicotinic receptors , 2005, The Journal of physiology.

[15]  O. Berger,et al.  Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. , 1997, Biophysical journal.

[16]  M. Winterhalter,et al.  The porin and the permeating antibiotic: a selective diffusion barrier in Gram-negative bacteria , 2008, Nature Reviews Microbiology.

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

[18]  Camilo Aponte-Santamaría,et al.  Crystal Structure of a Yeast Aquaporin at 1.15 Å Reveals a Novel Gating Mechanism , 2009, PLoS biology.

[19]  I. Vattulainen,et al.  Calculation of the electrostatic potential of lipid bilayers from molecular dynamics simulations: methodological issues. , 2009, The Journal of chemical physics.

[20]  R. Nicholas,et al.  Porin-Mediated Antibiotic Resistance in Neisseria gonorrhoeae: Ion, Solute, and Antibiotic Permeation through PIB Proteins with penB Mutations , 2006, Journal of bacteriology.

[21]  O. McManus,et al.  Ion Channels as Drug Targets: The Next GPCRs , 2008, The Journal of general physiology.

[22]  Kresten Lindorff-Larsen,et al.  Principles of conduction and hydrophobic gating in K+ channels , 2010, Proceedings of the National Academy of Sciences.

[23]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[24]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[25]  M. Colombini,et al.  Successful recovery of the normal electrophysiological properties of PorB (class 3) porin from Neisseria meningitidis after expression in Escherichia coli and renaturation. , 1998, Biochimica et biophysica acta.

[26]  T. Iverson,et al.  Structural basis for solute transport, nucleotide regulation, and immunological recognition of Neisseria meningitidis PorB , 2010, Proceedings of the National Academy of Sciences.

[27]  Andrew J Ballard,et al.  Replica exchange with nonequilibrium switches , 2009, Proceedings of the National Academy of Sciences.

[28]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[29]  Thomas B Woolf,et al.  Atomistic simulations of biologically realistic transmembrane potential gradients. , 2004, The Journal of chemical physics.

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

[31]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[32]  Shin-Ho Chung,et al.  Reservoir boundaries in Brownian dynamics simulations of ion channels. , 2002, Biophysical Journal.

[33]  Helmut Grubmüller,et al.  Kinetics, statistics, and energetics of lipid membrane electroporation studied by molecular dynamics simulations. , 2008, Biophysical journal.

[34]  Heike Wulff,et al.  Voltage-gated potassium channels as therapeutic targets , 2009, Nature Reviews Drug Discovery.

[35]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[36]  Maarten G. Wolf,et al.  g_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbation , 2010, J. Comput. Chem..

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

[38]  Electric field-controlled water permeation coupled to ion transport through a nanopore. , 2003, The Journal of chemical physics.

[39]  B. L. de Groot,et al.  The mechanism of proton exclusion in the aquaporin-1 water channel. , 2003, Journal of molecular biology.

[40]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[41]  O. Kepp,et al.  Bacterial Porin Disrupts Mitochondrial Membrane Potential and Sensitizes Host Cells to Apoptosis , 2009, PLoS pathogens.

[42]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[43]  Sohail Murad,et al.  A computer simulation of the classic experiment on osmosis and osmotic pressure , 1993 .