Correspondence to: Benoît Roux, GRTM, Department of Chemistry, University of Montreal, CP 6128, Succ. Centre-Ville, Montreal, Quebec, Canada H3C 3J7. Fax: 514-343-7586; E-mail:rouxb{at}plgcn.umontreal.ca.
The Perspectives on Ion Permeation in the June issue of The Journal of General Physiology provided a unique exposition of contrasting views regarding kinetic rate models, electrodiffusions, Brownian dynamics, and Poisson-Nernst-Planck (PNP). The debate is timely, healthy, and desirable, because recent progress in the determination of the three-dimensional structure of biological ion channels at atomic resolution gives a fresh impetus to efforts directed at understanding the fundamental principles governing ion permeation (
Misconception About Molecular Dynamics Simulations
Several statements about molecular dynamics (MD) indicate that the approach is misunderstood. The approach consists of constructing detailed atomic models of the macromolecular system of interest and, having described the microscopic forces with a potential function, use Newton's classical equation, F = MA, to literally "simulate" the dynamical motions of all the atoms as a function of time (
It is wrong to state that numerical integration of Newton's equation of motions is not reliable for times longer than several picoseconds because the calculated trajectories are very sensitive to initial conditions and round-off errors, and therefore diverge exponentially. This question was debated among theoretical chemists and physicists in the late 1960's or so, and it was shown that these concerns are perfectly solvable (
Misconception About Kinetic Rate Models
Kinetic models are constructed on the basis of two assumptions: first, it is assumed that the total configurational space of the whole system comprises a complete collection of distinct subspaces (the states), and, second, it is assumed that the system possesses no dynamical memory when it leaves one state to enter another (the Markov assumption). While it is always possible to define a complete collection of discrete states for any system, the Markov assumption is not necessarily valid (e.g., if there are no free energy barriers between different regions and the movement is purely diffusive). A kinetic rate model is valid as long as the long-time behavior of a molecular system can be described by a finite number of states whose lifetimes are exponentially distributed. Equilibrium properties do not rely on the Markov assumption, as it can be shown rigorously that the equilibrium probabilities for occupancy of multiply occupied channels have the familiar algebraic form for saturation behavior that is obtained from kinetic models with discrete states corresponding to a specific number of ions inside the pore (
The concept of a rate constant has been used and discussed for a long time (
The best aspect of kinetic rate models is their flexibility; one can adjust the rate constants constituting the model to fit most complex observed behaviors. The worst aspect of kinetic rate models is also their flexibility. The lack of internal constraints between the various rate constants makes it virtually impossible to guarantee a unique interpretation of what is going on at the microscopic level. The most sensible choice is thus to keep the models as simple as possible, even though a complex reality is not necessarily best described by a simple model.
Misconception About the Transmembrane Potential
There are many misconceptions about the transmembrane potential and its role in ion permeation. For example, in rate models, the transmembrane potential is usually treated as an external (constant?) field independent of the free energy profile (
Misconception About Poisson-Nernst-Planck Electrodiffusion
PNP electrodiffusion is an approximate theory combining the diffusion equation under the influence of an electric field, which itself is evaluated based on a continuum electrostatic approximation using the average (mean) charge density of the diffusing ions. Hence its "mean-field" characteristic. Although the earlier form probably dates back to Planck, a more complete description was developed by L. Onsager in the 1940's (
The best aspect of PNP is that it aims at doing everything at once: ionion, ionchannel, ionwater interactions, and the transmembrane potential are all treated in a consistent way. The worst aspect of PNP is that, while it aims at doing everything at once, it leaves out much of the atomic reality that we know is important at the microscopic level (e.g., van der Waals interactions, core repulsion, induction, hydrogen bonding, solvation structure, and protein flexibility). In practice, PNP is based on several simplifications: rigid channel structure, structureless dielectric solvent, and mean-field ionion interactions. If one is to adopt a continuum electrodiffusion approach, such simplifications are necessary to have partial differential equations that can be solved numerically. Interestingly, the debate about PNP often hinges on the use of the mean-field approximation to represent ionion interactions. While this is a nontrivial approximation (see below), the most fundamental problems with PNP are related to the approximations about channel rigidity and the representation of the solvent in terms of continuum electrostatics. Let us examine the physical significance of those approximations.
The flexibility of ion channels, as any proteins, plays an important role in its function (12; the interaction energy between two anions or two cations in bulk water are different, while continuum electrostatics is unable to make that distinction. In single file channels, the deviations from the continuum behavior are expected to be even more significant; e.g., ionion interactions in the gramicidin channel are species dependent even at a distance of 20 Å (
PNP is a consistent but approximate theory. It may, or may not, provide a useful picture of ion permeation because it relies on several physical approximations (rigid channel, continuum electrostatics, and mean-field ionion interactions) that are of unknown validity in the context in which they are used. Ultimately, the significance of that picture should not be expected to exceed that of the physical approximations upon which it is built, as is the case of the Born model of solvation, Debye-Hückel, or Poisson-Boltzmann theories. It would be useful to determine the validity of those physical approximations, but that would require more than algebraic mathematical considerations. It requires a comparison between the results from atomic models with explicit molecules and experimental data. Prediction of experimental results alone cannot reveal the limitations of PNP at the microscopic level.
What Is the Role of Theory in Biology?
Rather than focussing on the narrow question, "Which is the best: kinetic rate models or electrodiffusion?," one should ask the deeper question, "Where do these theoretical models stand within modern biology and ion channel science?" It is necessary to step back and try to address more fundamental questions about the role of theory in biology.
Understanding the function of biological systems is one the greatest scientific challenges of our times. To this day, biology remains primarily an experimental science, as it should. Complex biological systems are broken down in isolated elementary constituents, which are then analyzed for their structural and functional properties. These efforts represent more or less what one might (rather pompously) call the enterprise of modern scientific reductionism (
When there is no three-dimensional atomic structure available, the goal of theory in biology is to help formulate plausible and reasonable models to help organize the information from experimental data (e.g., currentvoltage-concentration relation). Simple models with a limited number of adjustable parameters are most desirable. There is no reason to complicate theoretical models when the atomic structure is not known.
When a three-dimensional atomic-resolution structure is available, the goal of theory is to analyze all the details that play an important role (van der Waals, electrostatic and hydrogen bonding interactions, protein flexibility, and hydration/dehydration processes) on the observed properties (permeation, selectivity, and gating). Can one predict observed properties such as a channel conductance from the atomic structure using only the fundamental laws of physics? At the present time, this is not really possible. The complexity of biological systems requires a hierarchy of inter-related levels of descriptions, the laws governing each level emerging from the fundamental behavior of the lower level: electrons and atoms obey, more or less, the laws of quantum mechanics as described by the Schrödinger equation, the forces acting on atoms and molecules are, more or less, described by the Born-Oppenheimer approximation, the dynamical trajectory of molecules follows, more or less, the laws of classical mechanics, the complexity of dense molecular systems leads, more or less, to chaotic diffusive motions guided by some sort of free energy potential surface, variations in the local composition may (or may not) follow a Markovian kinetic rate process, and so on. One should be able to continue like this all the way to the macroscopic physiological level (
Useful contributions in theoretical biology should clearly acknowledge their relationship to the other levels of description (above and below). One level below connects with a more fundamental basis, one level above connects to the macroscopic laws. The need to consider and incorporate several levels of description is an important methodological aspect of theoretical biology. Good theory teaches us all (experimentalists and theoreticians) how to better understand the behavior of biological systems by providing an integration of the information from different levels of description. For example: "How does an atomic model of an ion channel relate to a continuous electrodiffusion or to the kinetic rate model?" or "Can one make such a relation?" Detailed computations based on atomic models can contribute by helping in assessing quantitatively the relative importance of microscopic factors.
Theoretical models, at any level, are approximations. Excessive criticism of a theoretical model therefore becomes irrelevant when it is taken out of context. For example, kinetic rate models have been criticized extensively by the proponents of continuum electrodiffusion. But many of the criticisms of kinetic rate theory are largely unjustified and reflect misunderstandings and misconceptions about fundamental molecular statistical physics. The danger becomes that such criticism deters experimentalists from doing quantitative analysis of experimental data because they fear that they cannot attain the required level of theoretical rigor. Such an outcome would be very unfortunate because a quantitative characterization of biological systems (a painstaking task that is not always glamorous) is very important.
The ultimate goal should be to understand mechanisms better, not to develop a black-box that spits out numbers. A useful calculation can reveal important aspects of the function of an ion channel while it fails to reproduce exactly the conductance of the channel. In that sense, accurate reproduction of experimental data is desirable, but not necessarily the absolute criteria because no one at this point can reproduce the observed macroscopic behavior of a biological system starting from the most fundamental level. One must be patient and not ask too much too soon.
References
Becker, M.D., Koeppe, R.E., II, Andersen, O.S. 1992. Amino acid substitutions and ion channel function. Model-dependent conclusions. Biophys. J. 62:25-27[Medline].
Bock, G.R., Goode, J.A. 1998. The limits of reductionism in biology. Chichester, England, John Wiley & Sons, Ltd.
Born, M. 1920. Volumen und hydratationswarme der ionen. Z. Phys. 1:45.
Brooks, C.L., III, Karplus, M., Pettitt, B.M. 1988. Proteins. A theoretical perspective of dynamics, structure and thermodynamics. In Prigogine I., Rice S.A., eds. Advances in Chemical Physics. Vol. LXXI. New York, NY, John Wiley & Sons.
Chandler, D. 1978. Statistical mechanics of isomerization dynamics in liquids and the transition state approximation. J. Chem. Phys. 68:2959-2970.
Chang, G., Spencer, R.H., Lee, A.T., Barclay, M.T., Rees, D.C. 1998. Structure of the mscl homolog from mycobacterium tuberculosis: a gated mechanosensitive ion channel. Science. 282:2220-2226
Chen, D., Lear, J., Eisenberg, B. 1997. Permeation through an open channel: Poisson-Nernst-Planck theory of a synthetic ionic channel. Biophys. J. 72:97-116[Abstract].
Cowan, S.W., Schirmer, T., Rummel, G., Steiert, M., Gosh, R., Pauptit, R.A., Jansonius, J.N. 1992. Crystal structures explain functional properties of two E. coli porins. Nature. 358:727-733[Medline].
Crouzy, S., Woolf, T.B., Roux, B. 1994. A molecular dynamics study of gating in dioxolane-linked gramicidin A channels. Biophys. J. 67:1370-1386[Abstract].
Debye, P., Hückel, E. 1923. Zur theorie der elektrolyte II. Das Grenzgesetz für die elektrische leifähigkeit. Phys. Z. 24:305-325.
Doyle, D.A., Cabral, J.M., Pfuetzner, R.A., Kuo, A., Gulbis, J.M., Cohen, S.L., Chait, B.T., MacKinnon, R. 1998. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science 280:69-77
Hand, A.R., Steinmetz, P.R., Kohn, O.F. 1996. Scales of urine acidification: apical membrane-associated particles in turtle bladder. Kidney Int. 49:1655-1659[Medline].
Honig, B., Nicholls, A. 1995. Classical electrostatics in biology and chemistry. Science. 268:1144[Medline].
Kramers, H.A. 1940. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica 7:284-304.
Kurnikova, M.G., Coalson, R.D., Graf, P., Nitzan, A. 1998. A lattice relaxation algorithm for 3D Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel. Biophys. J 76:642-656
Läuger, P. 1972. Ion transport through pores: a rate theory analysis. Biochim. Biophys. Acta. 311:423-441.
Nonner, W., Chen, D., Eisenberg, B. 1999. Progress and prospects in permeation. J. Gen. Physiol 113:773-782
Berne, B.J., and R. Pecora. 1976. Dynamic light scattering. In Electrolyte solutions. John Wiley & Sons, New York, NY. Chapter 9.
Pettitt, B.M., Rossky, P.J. 1986. Alkali halides in water: ionsolvent and ionion potential of mean force at infinite dilution. J. Chem. Phys. 84:5836-5844.
Roux, B. 1999. Statistical mechanical equilibrium theory of selective ion channels. Biophys. J. 77:139-153
Roux, B., Karplus, M. 1995. Potential energy function for cationspeptides interactions: an Ab Initio study. J. Comp. Chem 16:690-704.
Roux, B., Prod'hom, B., Karplus, M. 1995. Ion transport in the gramicidin channel: molecular dynamics study of single and double occupancy. Biophys. J. 68:876-892[Abstract].
Roux, B., Yu, H.A., Karplus, M. 1990. Molecular basis for the born model of ion solvation. J. Phys. Chem 94:4683-4688.
Zwolinski, B.I., Eyring, H., Resese, C.E. 1949. Diffusion and membrane permeability. J. Physiol. Colloid. Chem 53:1426-1453.
|
|