Voltar
Carta Acesso aberto Revisado por pares
BIBTEX RIS

On the Potential Functions used in Molecular Dynamics Simulations of Ion Channels

2002; Elsevier BV; Volume: 82; Issue: 3 Linguagem: Inglês

10.1016/s0006-3495(02)75520-3

ISSN

1542-0086

Autores

Benoı̂t Roux, Simon Bernèche,

Tópico(s)

Nanopore and Nanochannel Transport Studies

Resumo

The determination of the structure of the KcsA K+ channel represents an extraordinary opportunity for understanding biological ion channels at the atomic level. In principle, molecular dynamics (MD) simulations based on detailed atomic models can complement the experimental data and help to characterize the microscopic factors that ultimately determine the permeation of ions through KcsA. A number of MD studies, broadly aimed at analyzing the dynamical motions of water molecules and ions in the KcsA channel, have now been reported (Guidoni et al., 1999Guidoni L. Torre V. Carloni P. Potassium and sodium binding to the outer mouth of the K+ channel.Biochemistry. 1999; 38: 8599-8604Crossref PubMed Scopus (139) Google Scholar; Allen et al., 1999Allen T.W. Kuyucak S. Chung S.H. Molecular dynamics study of the KcsA potassium channel.Biophys. J. 1999; 77: 2502-2516Abstract Full Text Full Text PDF PubMed Scopus (137) Google Scholar; Shrivastava and Sansom, 2000Shrivastava I.H. Sansom M.S. Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer.Biophys. J. 2000; 78: 557-570Abstract Full Text Full Text PDF PubMed Scopus (235) Google Scholar; Åqvist and Luzhkov, 2000Åqvist J. Luzhkov V. Ion permeation mechanism of the potassium channel.Nature. 2000; 404: 881-884Crossref PubMed Scopus (373) Google Scholar; Bernèche and Roux, 2000Bernèche S. Roux B. Molecular dynamics of the KcsA K(+) channel in a bilayer membrane.Biophys. J. 2000; 78: 2900-2917Abstract Full Text Full Text PDF PubMed Scopus (280) Google Scholar; Biggin et al., 2001Biggin P.C. Smith G.R. Shrivastava I. Choe S. Sansom M.S. Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations.Biochim. Biophys. Acta. 2001; 1510: 1-9Crossref PubMed Scopus (80) Google Scholar; Luzhkov and Åqvist, 2001Luzhkov V.B. Åqvist J. K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations.Biochim. Biophys. Acta. 2001; 1548: 194-202Crossref PubMed Scopus (101) Google Scholar; Crouzy et al., 2001Crouzy S. Bernèche S. Roux B. Extracellular blockade of K+ channels by TEA: results from molecular dynamics simulations of the KcsA channel.J. Gen. Phys. 2001; 118: 207-217Crossref PubMed Scopus (64) Google Scholar). The potential functions that were used to calculate the microscopic interatomic forces and generate the dynamical trajectory are listed in Table 1, where they can be seen to differ significantly. In particular, the atomic partial charges and the Lennard–Jones radii, which are at the heart of the potential function, varied widely. Furthermore, some include all atoms (AMBER and CHARMM PARAM22), whereas others are extended-atom models that treat only the polar hydrogens able to form hydrogen bonds explicitly (CHARMM PARAM19 and GROMOS). How these differences affect the results of MD calculations is an important concern of all scientists involved in investigations of ion channels, theoreticians and experimentalists alike. It is the goal of this short letter to discuss important aspects of potential functions related to MD studies of ion permeation.Table 1Potential energy function and MD simulations of KcsAMD SimulationsProtein and LipidType of Force FieldWaterIonsGuidoni et al., 1999Guidoni L. Torre V. Carloni P. Potassium and sodium binding to the outer mouth of the K+ channel.Biochemistry. 1999; 38: 8599-8604Crossref PubMed Scopus (139) Google ScholarAMBER*Cornell et al. (1995).All atomsTIP3†Jorgensen et al. (1983).Åqvist, 1990Åqvist J. Ion water interaction potential derived from free energy perturbation simulations.J. Phys. Chem. 1990; 94: 8021-8024Crossref Scopus (1350) Google ScholarBernèche and Roux, 2000Bernèche S. Roux B. Molecular dynamics of the KcsA K(+) channel in a bilayer membrane.Biophys. J. 2000; 78: 2900-2917Abstract Full Text Full Text PDF PubMed Scopus (280) Google ScholarCHARMM PARAM22‡Schlenkrich et al. (1996) for lipids and MacKerell et al. (1998) for proteins.All atomsTIP3†Jorgensen et al. (1983).Beglov and Roux, 1994Beglov D. Roux B. Finite representation of an infinite bulk system: solvent boundary potential for computer simulations.J. Chem. Phys. 1994; 100: 9050-9063Crossref Scopus (726) Google ScholarAllen et al., 1999Allen T.W. Kuyucak S. Chung S.H. Molecular dynamics study of the KcsA potassium channel.Biophys. J. 1999; 77: 2502-2516Abstract Full Text Full Text PDF PubMed Scopus (137) Google ScholarCHARMM PARAM19§Brooks et al. (1983).Extended atomsST2¶Stillinger and Rahman (1974).Heinzinger, 1985Heinzinger K. Computer simulations of aqueous electrolyte solutions.Physica. 1985; 131B: 196-216Google ScholarShrivastava and Sansom, 2000Shrivastava I.H. Sansom M.S. Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer.Biophys. J. 2000; 78: 557-570Abstract Full Text Full Text PDF PubMed Scopus (235) Google ScholarGROMOS∥Hermans et al. (1984).Extended atomsSPC**Berendsen et al. (1981).Straatsma et al., 1988Straatsma T.P. Berendsen H.J.C. Postma J.P.M. Free energy of ionic hydration: analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations.J. Chem. Phys. 1988; 89: 5876-5886Crossref Scopus (400) Google ScholarÅqvist and Luzhkov, 2000Åqvist J. Luzhkov V. Ion permeation mechanism of the potassium channel.Nature. 2000; 404: 881-884Crossref PubMed Scopus (373) Google ScholarGROMOS∥Hermans et al. (1984).Extended atomsSPC**Berendsen et al. (1981).Åqvist, 1990Åqvist J. Ion water interaction potential derived from free energy perturbation simulations.J. Phys. Chem. 1990; 94: 8021-8024Crossref Scopus (1350) Google Scholar* Cornell et al., 1995Cornell W.D. Cieplak P. Bayly C.I. Gould I.R. Merz Jr., K.M. Ferguson D.M. Spellmeyer D.C. Fox T. Caldwell J.W. Kollman P.A. A second generation force field for the simulation of proteins and nucleic acids.J. Am. Chem. Soc. 1995; 117: 5179-5197Crossref Scopus (10933) Google Scholar.† Jorgensen et al., 1983Jorgensen W.L. Chandrasekhar J. Madura J.D. Impey R.W. Klein M.L. Comparison of simple potential functions for simulating liquid water.J. Chem. Phys. 1983; 79: 926-935Crossref Scopus (26693) Google Scholar.‡ Schlenkrich et al., 1996Schlenkrich M.J. Brickmann J. MacKerell Jr., A.D. Karplus M. An empirical potential energy function for phospholipids: criteria for parameters optimization and applications.in: Merz K.M. Roux B. Biological Membranes. A Molecular Perspective from Computation and Experiment. Birkhauser, Boston, MA1996: 31-81Crossref Google Scholar for lipids and MacKerell et al., 1998MacKerell Jr., A.D. Bashford D. Bellot M. Dunbrack R.L. Evanseck J.D. Field M.J. Fischer S. Gao J. Guo H. Joseph-McCarthy D. Ha S. Kuchnir L. Kuczera K. Lau F.T.K. Mattos C. Michnick S. Ngo T. Nguyen D.T. Prodhom B. Reiher III, W.E. Roux B. Schlenkrich M. Smith J. Stote R. Straub J. Watanabe M. Wiorkiewicz-Kuczera J. Karplus M. All-atom empirical potential for molecular modeling and dynamics studies of proteins.J. Phys. Chems. B. 1998; 102: 3586-3616Crossref Scopus (10704) Google Scholar for proteins.§ Brooks et al., 1983Brooks B.R. Bruccoleri R.E. Olafson B.D. States D.J. Swaminathan S. Karplus M. CHARMM: a program for macromolecular energy minimization and dynamics calculations.J. Comput. Chem. 1983; 4: 187-217Crossref Scopus (13371) Google Scholar.¶ Stillinger and Rahman, 1974Stillinger F.H. Rahman A. Improved simulation of liquid water by molecular dynamics.J. Chem. Phys. 1974; 60: 1545-1557Crossref Scopus (950) Google Scholar.∥ Hermans et al., 1984Hermans J. Berendsen H.J.C. van Gunsteren W.F. Postma J.P.M. A consistent empirical potential for water–protein interactions.Biopol. 1984; 23: 1513-1518Crossref Scopus (679) Google Scholar.** Berendsen et al., 1981Berendsen H.J.C. Postma J.P.M. van Gunsteren W.F. Hermans J. Interaction models for water in relation to proteins hydration.in: Pullman B. Intermolecular Forces. Reidel, Dordrecht, The Netherlands1981: 331-342Crossref Google Scholar. Open table in a new tab For meaningful theoretical studies of permeation, it is necessary to have a potential energy function providing a realistic and accurate representation of the microscopic interactions. In practice, this presents a difficult challenge. The permeation process through KcsA involves the partial dehydration of a K+ ion, followed by the translocation through the interior of a narrow pore of 12-Å-long, lined by backbone carbonyl oxygens, which acts as a selectivity filter (Doyle et al., 1998Doyle D.A. Cabral J.M. Pfuetzner R.A. Kuo A. Gulbis J.M. Cohen S.L. Chait B.T. MacKinnon R. The structure of the potassium channel: molecular basis of K+ conduction and selectivity.Science. 1998; 280: 69-77Crossref PubMed Scopus (5468) Google Scholar). Thus, the conductance and selectivity of the KcsA channel results from a delicate balance of very strong microsopic interactions, the large energetic loss of dehydration being roughly compensated by coordination with main chain carbonyl oxygens. Gas phase experiments on model systems provide the most direct information concerning the individual microscopic interactions (Džidić and Kebarle, 1970Džidi I. Kebarle P. Hydration of the alkali ions in the gas phase. Enthalpies and entropies of reactions M+(H2O)n−1 + H2O = M+(H2O)n.J. Phys. Chem. 1970; 74: 1466-1474Crossref Scopus (668) Google Scholar; Klassen et al., 1996Klassen J.S. Anderson S.G. Blades A.T. Kebarle P. Reaction enthalpies for M(+)L = M(+) + L, where M(+) = Na+ and K+ and L equals acetamide, N-methylacetamide, N,N-dimethylacetamide, glycine, and glycylglycine, from determinations of the collision-induced dissociation thresholds.J. Phys. Chem. 1996; 100: 14218-14227Crossref Scopus (148) Google Scholar). High-level quantum-mechanical ab initio calculations can also be used to supplement the (often scarce) information available from experiments (Roux and Karplus, 1995Roux B. Karplus M. Potential energy function for cations–peptides interactions: an ab initio study.J. Comp. Chem. 1995; 16: 690-704Crossref Scopus (62) Google Scholar). The interaction of ions with a single water molecule, or with a single isolated N-methylacetamide (NMA) molecule, an excellent model of the backbone carbonyl of proteins, is of particular interest. The most important microscopic interactions energies for ion permeation through the K+ channel are given in Table 2. Despite the considerable uncertainty in the experimental data and the ab initio calculations, both clearly indicate that the interaction of cations with a single NMA is substantially larger than with a single water molecule. The binding enthalpy of K+ with a water molecule is 17.9 kcal/mole, whereas it is roughly 25–30 kcal/mole with NMA. The interactions are even larger in the case of Na+. This general trend is generally reproduced by all the potential functions, with the exception of GROMOS (Hermans et al., 1984Hermans J. Berendsen H.J.C. van Gunsteren W.F. Postma J.P.M. A consistent empirical potential for water–protein interactions.Biopol. 1984; 23: 1513-1518Crossref Scopus (679) Google Scholar). In this case, the interaction of K+ and Na+ with a single NMA is actually smaller than the interaction with a single water molecule. The difference in the interaction energy is directly related to the atomic charges assigned to the peptide backbone, i.e., the atomic charges from GROMOS (Hermans et al., 1984Hermans J. Berendsen H.J.C. van Gunsteren W.F. Postma J.P.M. A consistent empirical potential for water–protein interactions.Biopol. 1984; 23: 1513-1518Crossref Scopus (679) Google Scholar) are about 60% to 75% relative to those from AMBER (Cornell et al., 1995Cornell W.D. Cieplak P. Bayly C.I. Gould I.R. Merz Jr., K.M. Ferguson D.M. Spellmeyer D.C. Fox T. Caldwell J.W. Kollman P.A. A second generation force field for the simulation of proteins and nucleic acids.J. Am. Chem. Soc. 1995; 117: 5179-5197Crossref Scopus (10933) Google Scholar), CHARMM PARAM19 (Brooks et al., 1983Brooks B.R. Bruccoleri R.E. Olafson B.D. States D.J. Swaminathan S. Karplus M. CHARMM: a program for macromolecular energy minimization and dynamics calculations.J. Comput. Chem. 1983; 4: 187-217Crossref Scopus (13371) Google Scholar), or CHARMM PARAM22 (MacKerell et al., 1998MacKerell Jr., A.D. Bashford D. Bellot M. Dunbrack R.L. Evanseck J.D. Field M.J. Fischer S. Gao J. Guo H. Joseph-McCarthy D. Ha S. Kuchnir L. Kuczera K. Lau F.T.K. Mattos C. Michnick S. Ngo T. Nguyen D.T. Prodhom B. Reiher III, W.E. Roux B. Schlenkrich M. Smith J. Stote R. Straub J. Watanabe M. Wiorkiewicz-Kuczera J. Karplus M. All-atom empirical potential for molecular modeling and dynamics studies of proteins.J. Phys. Chems. B. 1998; 102: 3586-3616Crossref Scopus (10704) Google Scholar).Table 2Microscopic interactions (kcal/mol)*The interactions energies based on the different force fields were calculated by us assuming a rigid geometry of the water or NMA molecule. When unavailable for a given potential function, the atomic partial charges of NMA were deduced from those of a glycine dipeptide.K+Na+WaterNMAWaterNMAReference17.928.3–32.324.033.7–39.0Gas phase exp (Džidić and Kebarle, 1970Džidi I. Kebarle P. Hydration of the alkali ions in the gas phase. Enthalpies and entropies of reactions M+(H2O)n−1 + H2O = M+(H2O)n.J. Phys. Chem. 1970; 74: 1466-1474Crossref Scopus (668) Google Scholar; Klassen et al., 1996Klassen J.S. Anderson S.G. Blades A.T. Kebarle P. Reaction enthalpies for M(+)L = M(+) + L, where M(+) = Na+ and K+ and L equals acetamide, N-methylacetamide, N,N-dimethylacetamide, glycine, and glycylglycine, from determinations of the collision-induced dissociation thresholds.J. Phys. Chem. 1996; 100: 14218-14227Crossref Scopus (148) Google Scholar)15.9–17.624.8–31.724.0–25.838.4–40.4Ab initio (Roux and Karplus, 1995Roux B. Karplus M. Potential energy function for cations–peptides interactions: an ab initio study.J. Comp. Chem. 1995; 16: 690-704Crossref Scopus (62) Google Scholar)18.223.723.229.5Guidoni et al., 1999Guidoni L. Torre V. Carloni P. Potassium and sodium binding to the outer mouth of the K+ channel.Biochemistry. 1999; 38: 8599-8604Crossref PubMed Scopus (139) Google Scholar18.924.125.530.1Bernèche and Roux, 2000Bernèche S. Roux B. Molecular dynamics of the KcsA K(+) channel in a bilayer membrane.Biophys. J. 2000; 78: 2900-2917Abstract Full Text Full Text PDF PubMed Scopus (280) Google Scholar18.321.924.827.9Allen et al., 1999Allen T.W. Kuyucak S. Chung S.H. Molecular dynamics study of the KcsA potassium channel.Biophys. J. 1999; 77: 2502-2516Abstract Full Text Full Text PDF PubMed Scopus (137) Google Scholar17.816.622.820.6Åqvist and Luzhkov, 2000Åqvist J. Luzhkov V. Ion permeation mechanism of the potassium channel.Nature. 2000; 404: 881-884Crossref PubMed Scopus (373) Google Scholar17.616.826.323.6Shrivastava and Sansom, 2000Shrivastava I.H. Sansom M.S. Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer.Biophys. J. 2000; 78: 557-570Abstract Full Text Full Text PDF PubMed Scopus (235) Google Scholar* The interactions energies based on the different force fields were calculated by us assuming a rigid geometry of the water or NMA molecule. When unavailable for a given potential function, the atomic partial charges of NMA were deduced from those of a glycine dipeptide. Open table in a new tab Conductance and selectivity are primarily governed by relative free energies. For this reason, it is essential to consider also thermodynamics properties in the parametrization of the potential function in addition to the microscopic interactions. The solvation free energy of cations in liquid water and liquid NMA are particularly important for calibrating a potential function. In the case of water, it is possible to reproduce both the microsopic interactions and the solvation free energy of ions with the current potential functions (Straatsma et al., 1988Straatsma T.P. Berendsen H.J.C. Postma J.P.M. Free energy of ionic hydration: analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations.J. Chem. Phys. 1988; 89: 5876-5886Crossref Scopus (400) Google Scholar; Åqvist, 1990Åqvist J. Ion water interaction potential derived from free energy perturbation simulations.J. Phys. Chem. 1990; 94: 8021-8024Crossref Scopus (1350) Google Scholar; Beglov and Roux, 1994Beglov D. Roux B. Finite representation of an infinite bulk system: solvent boundary potential for computer simulations.J. Chem. Phys. 1994; 100: 9050-9063Crossref Scopus (726) Google Scholar). For example, the solvation free energy of K+ in liquid water is ∼80 kcal/mol (Dorman et al., 1996Dorman V. Partenskii M.B. Jordan P.C. A semi-microscopic Monte Carlo study of permeation energetics in a gramicidin-like channel: the origin of cation selectivity.Biophys. J. 1996; 70: 121-134Abstract Full Text PDF PubMed Scopus (56) Google Scholar) (though there is considerable uncertainty, see Pliego and Riveros, 2000Pliego J.R. Riveros J.M. New values for the absolute solvation free energy of univalent ions in aqueous solution.Chem. Phys. Lett. 2000; 332: 597-602Crossref Scopus (106) Google Scholar). Such a value can be reproduced quite well with a potential function yielding a microsopic interaction with a single water molecule on the order of 17–18 kcal/mol (Straatsma et al., 1988Straatsma T.P. Berendsen H.J.C. Postma J.P.M. Free energy of ionic hydration: analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations.J. Chem. Phys. 1988; 89: 5876-5886Crossref Scopus (400) Google Scholar; Åqvist, 1990Åqvist J. Ion water interaction potential derived from free energy perturbation simulations.J. Phys. Chem. 1990; 94: 8021-8024Crossref Scopus (1350) Google Scholar; Beglov and Roux, 1994Beglov D. Roux B. Finite representation of an infinite bulk system: solvent boundary potential for computer simulations.J. Chem. Phys. 1994; 100: 9050-9063Crossref Scopus (726) Google Scholar). In contrast, MD free energy calculations indicate that it is very difficult to reproduce both the cation–NMA microscopic energy and the solvation free energy in liquid NMA with current biomolecular potential functions. For example, the CHARMM PARAM22 potential function, which gives an interaction energy of 24.1 kcal/mol with a single NMA, yields a free energy of ∼88 kcal/mol in liquid NMA (S. Bernèche and B. Roux, unpublished results). Although the solvation free energy of K+ in liquid NMA is not known experimentally, data from other liquid amides suggests that such a large value is unrealistic and that a reasonable estimate should be ∼80–82 kcal/mol (Cox et al., 1974Cox B.G. Hedwig G.R. Parker A.J. Watts D.W. Solvation of ions XIX. Thermodynamic properties for transfer of single ions between protic and dipolar aprotic solvents.Aust. J. Chem. 1974; 27: 477-501Crossref Scopus (420) Google Scholar). For a given potential function, the calculated ion solvation free energy in liquid NMA is expected to be reflected directly upon the stability of K+ in the selectivity filter during MD simulations of the KcsA channel. Therefore, the present analysis suggests that the K+ ions bind too strongly to KcsA by ∼5–10 kcal/mol in MD simulations based on the all-atoms potential function AMBER and CHARMM PARAM22, such as used by Guidoni et al., 1999Guidoni L. Torre V. Carloni P. Potassium and sodium binding to the outer mouth of the K+ channel.Biochemistry. 1999; 38: 8599-8604Crossref PubMed Scopus (139) Google Scholar and Bernèche and Roux, 2000Bernèche S. Roux B. Molecular dynamics of the KcsA K(+) channel in a bilayer membrane.Biophys. J. 2000; 78: 2900-2917Abstract Full Text Full Text PDF PubMed Scopus (280) Google Scholar, respectively. In contrast, because the microscopic interaction energy of K+ with a single NMA is only on the order of 16–17 kcal/mol (see Table 2), the K+ ions bind probably too weakly to KcsA by as much as 20 kcal/mol in MD simulations based on the extended-atom GROMOS potential function such as used by Åqvist (Åqvist and Luzhkov, 2000Åqvist J. Luzhkov V. Ion permeation mechanism of the potassium channel.Nature. 2000; 404: 881-884Crossref PubMed Scopus (373) Google Scholar; Luzhkov and Åqvist, 2001Luzhkov V.B. Åqvist J. K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations.Biochim. Biophys. Acta. 2001; 1548: 194-202Crossref PubMed Scopus (101) Google Scholar) and Sansom (Shrivastava and Sansom, 2000Shrivastava I.H. Sansom M.S. Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer.Biophys. J. 2000; 78: 557-570Abstract Full Text Full Text PDF PubMed Scopus (235) Google Scholar; Biggin et al., 2001Biggin P.C. Smith G.R. Shrivastava I. Choe S. Sansom M.S. Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations.Biochim. Biophys. Acta. 2001; 1510: 1-9Crossref PubMed Scopus (80) Google Scholar). To obtain a free energy of ∼80 kcal/mol in liquid NMA, one can adjust the Lennard–Jones parameters of the cation–carbonyl oxygen pairs and reduce the microscopic cation–NMA interaction energy to ∼21.6 kcal/mol (S. Bernèche and B. Roux, unpublished results). This is one way to parametrize and calibrate the potential function for theoretical studies of ion permeation through KcsA. Clearly, if the potential function was an exact representation of the Born–Oppenheimer energy surface, success in reproducing the microsopic interactions would automatically lead to accurate thermodynamic properties. But current biomolecular potential functions try to account for many-body polarization effects in an average way using an effective parametrization of the atomic partial charges. Because of this approximation, the optimal parametrization is the result of a compromise between an accurate representation of the microscopic energies and bulk solvation properties. We believe that such potential functions can yield meaningful results of semi-quantitative accuracy. Recently, we have taken these factors into consideration in calibrating the potential function for a calculation of the free energy surface governing conduction of K+ ions through the selectivity filter of the KcsA K+ channel (Bernèche and Roux, 2001Bernèche S. Roux B. Energetics of ion conduction through the k+ channel.Nature. 2001; 414: 73-77Crossref PubMed Scopus (620) Google Scholar). In the particular case of this study, it should be stressed that meaningful results were not obtained until the potential function was adjusted to reproduce the correct free energies of K+ in liquid water and liquid NMA. In general, it ought to be possible to calibrate any potential function to reproduce solvation free energies using a similar approach (though the significantly underestimated ion–NMA interaction energy based on the GROMOS force field might require some modifications of the atomic charges). Further analysis suggest that the situation might be more difficult in the case of a small cation such as Na+ (Roux, 1993Roux B. Nonadditivity in cation–peptide interactions: a molecular dynamics and ab initio study of Na+ in the gramicidin channel.Chem. Phys. Lett. 1993; 212: 231-240Crossref Scopus (43) Google Scholar), suggesting that a quantitative simulation of the microscopic factors governing ion selectivity is probably beyond the ability of current biomolecular potential function. Ultimately, the influence of nonadditive many-body polarization should be viewed in a wider perspective. At the present time, computational chemists and theoreticians are actively pursuing the development of a new generation of force fields that will include induced polarization for computational studies of biological systems (Halgren and Damm, 2001Halgren T.A. Damm W. Polarizable force fields.Curr. Opin. Struct. Biol. 2001; 11: 236-242Crossref PubMed Scopus (409) Google Scholar). But much more work is needed before such potential functions are ready to be used in simulations of biological ion channels. Meanwhile, we believe that MD studies of ion channels can still yield meaningful results, as long as they are based on effective potential functions that have been calibrated to correctly reproduce solvation free energies.

Referência(s)