Ion permeation and selectivity of OmpF porin: A theoretical study based on molecular dynamics, brownian dynamics, and continuum electrodiffusion theory

被引:315
作者
Im, W [1 ]
Roux, B [1 ]
机构
[1] Cornell Univ, Dept Biochem, Weill Med Coll, New York, NY 10021 USA
基金
美国国家卫生研究院;
关键词
ion channel; membrane protein; diffusion; reversal potential; ion flux;
D O I
10.1016/S0022-2836(02)00778-7
中图分类号
Q5 [生物化学]; Q7 [分子生物学];
学科分类号
071010 ; 081704 ;
摘要
Three different theoretical approaches are used and compared to refine our understanding of ion permeation through the channel formed by OmpF porin from Escherichia coli. Those approaches are all-atom molecular dynamics (MD) in which ions, solvent, and lipids are represented explicitly, Brownian dynamics (BD) in which ions are represented explicitly, while solvent and lipids are represented as featureless dielectrics, and Poisson-Nernst-Planck (PNP) electrodiffusion theory in which both solvent and local ion concentrations are represented as a continuum. First, the ability of the different theoretical approaches in reproducing the equilibrium average ion density distribution in OmpF porin bathed by a 1 M KCl symmetric salt solution is examined. Under those conditions the PNP theory is equivalent to the non-linear Poisson-Boltzmann (PB) theory. Analysis shows that all the three approaches are able to capture the important electrostatic interactions between ions and the charge distribution of the channel that govern ion permeation and selectivity in OmpF. The K+ and Cl- density distributions obtained from the three approaches are very consistent with one another, which suggests that a treatment on the basis of a rigid protein and continuum dielectric solvent is valid in the case of OmpF. Interestingly, both BD and continuum electrostatics reproduce the distinct left-handed twisted ion pathways for K+ and Cl- extending over the length of the pore which were observed previously in MD. Equilibrium BD simulations in the grand canonical ensemble indicate that the channel is very attractive for cations, particularly at low salt concentration. On an average there is 1.55 K+ inside the pore in 10 mM KCl. Remarkably, there is still 0.17 K+ on average inside the pore even at a concentration as low as 1 muM KCl. Secondly, non-equilibrium ion flow through OmpF is calculated using BD and PNP and compared with experimental data. The channel conductance in 0.2 M and 1 M KCl calculated using BD is in excellent accord with the experimental data. The calculations reproduce the experimentally well-known conductance-concentration relation and also reveal an asymmetry in the channel conductance (a larger conductance is observed under a positive transmembrane potential). Calculations of the channel conductance for three mutants (R168A, R132A, and K16A) in 1 M KCl suggest that the asymmetry in the channel conductance arises mostly from the permanent charge distribution of the channel rather than the shape of the pore itself. Lastly, the calculated reversal potential in a tenfold salt gradient (0.1:1 M KCl) is 27.4(+/-1.3) mV (BD) and 22.1(+/-0.6) mV (PNP), in excellent accord with the experimental value of 24.3 mV. Although most of the results from PNP are qualitatively reasonable, the calculated channel conductance is about 50% higher than that calculated from BID probably because of a lack of some dynamical ion-ion correlations. (C) 2002 Elsevier Science Ltd. All rights reserved.
引用
收藏
页码:851 / 869
页数:19
相关论文
共 82 条
[1]  
Allen M. P., 1989, COMPUTER SIMULATION
[2]   Molecular dynamics estimates of ion diffusion in model hydrophobic and KcsA potassium channels [J].
Allen, TW ;
Kuyucak, S ;
Chung, SH .
BIOPHYSICAL CHEMISTRY, 2000, 86 (01) :1-14
[3]   FINITE REPRESENTATION OF AN INFINITE BULK SYSTEM - SOLVENT BOUNDARY POTENTIAL FOR COMPUTER-SIMULATIONS [J].
BEGLOV, D ;
ROUX, B .
JOURNAL OF CHEMICAL PHYSICS, 1994, 100 (12) :9050-9063
[4]   PERMEATION OF HYDROPHILIC MOLECULES THROUGH THE OUTER-MEMBRANE OF GRAM-NEGATIVE BACTERIA - REVIEW ON BACTERIAL PORINS [J].
BENZ, R ;
BAUER, K .
EUROPEAN JOURNAL OF BIOCHEMISTRY, 1988, 176 (01) :1-19
[5]   ION SELECTIVITY OF GRAM-NEGATIVE BACTERIAL PORINS [J].
BENZ, R ;
SCHMID, A ;
HANCOCK, REW .
JOURNAL OF BACTERIOLOGY, 1985, 162 (02) :722-727
[6]   Alteration of pore properties of Escherichia coli OmpF induced by mutation of key residues in anti-loop 3 region [J].
Bredin, J ;
Saint, N ;
Malléa, M ;
Dé, E ;
Molle, G ;
Pagès, JM ;
Simonet, V .
BIOCHEMICAL JOURNAL, 2002, 363 :521-528
[7]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[8]   Three-dimensional Poisson-Nernst-Planck theory studies:: Influence of membrane electrostatics on gramicidin A channel conductance [J].
Cárdenas, AE ;
Coalson, RD ;
Kurnikova, MG .
BIOPHYSICAL JOURNAL, 2000, 79 (01) :80-93
[9]   Stochastic problems in physics and astronomy [J].
Chandrasekhar, S .
REVIEWS OF MODERN PHYSICS, 1943, 15 (01) :0001-0089
[10]   THE NATURE OF ION AND WATER BARRIER CROSSINGS IN A SIMULATED ION CHANNEL [J].
CHIU, SW ;
NOVOTNY, JA ;
JAKOBSSON, E .
BIOPHYSICAL JOURNAL, 1993, 64 (01) :98-108