Molecular dynamics study of ion transport in transmembrane protein channels

Abstract

Ion transport through biological membranes often takes place via pore-like protein channels. The elementary process of this transport can be described as a motion of the ion in a quasi-periodic multi-well potential. In this study molecular dynamics simulations of ion transport in a model channel were performed in order to test the validity of reaction-rate theory for this process. The channel is modelled as a hexagonal helix of infinite length, and the ligand groups interacting with the ion are represented by dipoles lining the central hole of the channel. The dipoles interact electrostatically with each other and are allowed to oscillate around an equilibrium orientation. The coupled equations of motion for the ion and the dipoles were solved simultaneously with the aid of a numerical integration procedure. From the calculated ion trajectories it is seen that, particularly at low temperatures, the ion oscillates back and forth in the trapping site many times before it leaves the site and jumps over the barrier. The observed oscillation frequency was found to be virtually temperature-independent (nu 0 approximately equal to 2 X 10(12) s-1) so that the strong increase of transport rate with temperature results almost exclusively from the Arrhenius-type exponential dependence of jump probability w on 1/T. At higher temperatures simultaneous jumps over several barriers occasionally occur. Although the exponential form of w(T) was in agreement with the predictions of rate theory, the activation energy Ea as determined from w(T) was different from the barrier height which was calculated from the static potential of the ion in the channel; the actual transport rate was 1 X 10(3) times higher than the rate predicted from the calculated barrier height. This observation was interpreted by the notion that ion transport in the channel is strongly influenced by thermal fluctuations in the conformation of the ligand system which in turn give rise to fluctuations of barrier height.