Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. II. Rates and mechanisms of water transport

Proton mobilities in water and in different stereoisomers of covalently linked gramicidin A channels

Proton conductivities in bulk solution (lambda(H)) and single-channel proton conductances (g(H)) in two different stereoisomers of the dioxolane-linked gramicidin A channel (the SS and RR dimers) were measured in a wide range of bulk proton concentrations ([H], 0.1-8000 mM). Proton mobilities (micro(H)) in water as well as in the SS and RR dimers were calculated from the conductivity data. In the concentration range of 0.1-2000 mM, a straight line with a slope of 0.75 describes the log (g(H))-log ([H]) relationship in the SS dimer. At [H] > 2000 mM, saturation is followed by a decline in g(H). The g(H)-[H] relationship in the SS dimer is qualitatively similar to the [H] dependence of lambda(H). However, the slope of the straight line in the log(lambda(H))-log([H]) plot is 0.96, indicating that the rate-limiting step for proton conduction through the SS dimer is not the diffusion of protons in bulk solution. The significant difference between the slopes of those linear relationships accounts for the faster decline of micro(H) as a function of [H] in the SS dimer in relation to bulk solution. In the high range of [H], saturation and decline of g(H) in the SS dimer can be accounted for by the significant decrease of micro(H) in bulk solution. At any given [H], g(H) in the RR dimer is significantly smaller than in the SS. Moreover, the g(H)-[H] relationship in the RR stereoisomer is qualitatively different from that in the SS. Between 1 and 50 mM [H], g(H) can be fitted with an adsorption isotherm, suggesting the presence of a proton-binding site inside the pore (pK(a) approximately 2), which limits proton exit from the channel. At 100 mM < [H] < 3000 mM, g(H) increases linearly with [H]. The distinctive shape of the g(H)-[H] relationship in the RR dimer suggests that the channel can be occupied simultaneously by more than one proton. At higher [H], the saturation and decline of g(H) in the RR dimer reflect the properties of micro(H) in bulk solution. In the entire range of [H], protons seem to cross the SS and RR channels via a Grotthuss-like mechanism. The rate-limiting step for proton transfer in the SS dimer is probably the membrane-channel/bulk solution interface. It is also proposed that the smaller g(H) in the RR dimer is the consequence of a different organization and dynamics of the H-bonded network of water molecules inside the pore of the channel, resulting in a slower proton transfer and multiple pore occupancy by protons.

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.

Can we use rate constants and state models to describe ion transport through gramicidin channels?

Can we use rate constants and state models to describe ion transport through gramicidin channels? Maybe, but only if rate constants are just proportionality constants between rates and probabilities of observing states of the channel. This approach is natural if the system of channel plus ions (plus water) is almost always in one or another of a small number of identifiable states. Many features of ion transport through gramicidin, including the conductance-concentration relationship, concentration-dependent permeability ratios, anomalous mole fraction effect and to some extent flux ratio exponents, are consistent with a description in which there are four occupation 'states' of the pore: only water; an ion at one end; an ion at the other; and ions at both ends. Current-voltage relationships can (and must) also be fitted, but until there is a theory to predict the potential dependence of the rate constants this success will remain hollow. Other features have resisted interpretation. These include the failures to determine 'binding constants' consistent with all the data; the variation of flux ratio exponents with ion type; and, probably, the variation of the currents with asymmetrical ion concentrations. Nevertheless, state models still have one attractive feature, they allow consideration of the effects that one ion within the pore has on the movements of another.