The reaction path and free energy profile of Na+ were computed in the interior of the channel protein gramicidin, with the program MOIL. Gramicidin was represented in atomic detail, but surrounding water and lipid molecules were not included. Thus, only short range interactions were investigated. The permeation path of the ion was an irregular spiral, far from a straight line. Permeation cannot be described by motions of a single Na+ ion. The minimal energy path includes significant motion of water and channel atoms as well as motion of the permeating ion. We think of permeation as motion of a permion, a quasi-particle that includes the many body character of the permeation process, comparable with quasi-particles like holes, phonons, and electrons of solid-state physics. Na+ is accompanied by a plug of water molecules, and motions of water, Na+, and the atoms of gramicidin are highly correlated. The permion moves like a linear polymer made of waters and ion linked and moving coherently along a zigzag line, following the reptation mechanism of polymer transport. The effective mass, free energy, and memory kernel (of the integral describing time-dependent friction) of short range interactions were calculated. The effective mass of the permion (properly normalized) is much less than Na+. Friction varies substantially along the path. The free energy profile has two deep minima and several maxima. In certain regions, the dominant motions along the reaction path are those of the channel protein, not the permeating ion: there, ion waits while the other atoms move. At these waiting sites, the permion's motion along the reaction path is a displacement of the atoms of gramicidin that prepare the way for the Na+ ion.
[1]
R. Eisenberg,et al.
Diffusion as a chemical reaction: Stochastic trajectories between fixed concentrations
,
1995
.
[2]
E Jakobsson,et al.
Water and polypeptide conformations in the gramicidin channel. A molecular dynamics study.
,
1989,
Biophysical journal.
[3]
U. Singh,et al.
A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS
,
1984
.
[4]
A. Pullman.
Energy profiles in the gramicidin A channel
,
1987,
Quarterly Reviews of Biophysics.
[5]
M Karplus,et al.
Molecular dynamics simulations of the gramicidin channel.
,
1994,
Annual review of biophysics and biomolecular structure.
[6]
B. Berne,et al.
Spatial dependence of time‐dependent friction for pair diffusion in a simple fluid
,
1990
.
[7]
W. L. Jorgensen,et al.
Comparison of simple potential functions for simulating liquid water
,
1983
.
[8]
P. H. Berens,et al.
Structure and dynamics of ion transport through gramicidin A.
,
1984,
Biophysical journal.
[9]
R. Eisenberg,et al.
Charges, currents, and potentials in ionic channels of one conformation.
,
1993,
Biophysical journal.
[10]
G. N. Patey,et al.
The free energy of spheres with dipoles: Monte Carlo with multistage sampling
,
1973
.