Gramicidin A Channel as a Test Ground for Molecular Dynamics Force Fields
We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyrist...
Saved in:
Published in | Biophysical journal Vol. 84; no. 4; pp. 2159 - 2168 |
---|---|
Main Authors | , , , |
Format | Journal Article |
Language | English |
Published |
United States
Elsevier Inc
01.04.2003
Biophysical Society |
Subjects | |
Online Access | Get full text |
Cover
Loading…
Summary: | We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K
+ ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K
+ ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60
ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose. |
---|---|
Bibliography: | ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23 Toby Allen's present address is Dept. of Biochemistry, Weill Medical College of Cornell University, 1300 York Ave., New York, NY 10021. Address reprint requests to Dr. Serdar Kuyucak, Dept. of Theoretical Physics, Research School of Physical Sciences, Australian National University, Canberra, A.C.T. 0200, Australia. Tel.: 61-2-6125-2969; Fax: 61-2-6125-4676; E-mail: serdar.kuyucak@anu.edu.au. |
ISSN: | 0006-3495 1542-0086 |
DOI: | 10.1016/S0006-3495(03)75022-X |