Water structure in the Gramicidin A transmembrane channel

Voltage-gated proton channels and other proton transfer pathways

Proton channels exist in a wide variety of membrane proteins where they transport protons rapidly and efficiently. Usually the proton pathway is formed mainly by water molecules present in the protein, but its function is regulated by titratable groups on critical amino acid residues in the pathway. All proton channels conduct protons by a hydrogen-bonded chain mechanism in which the proton hops from one water or titratable group to the next. Voltage-gated proton channels represent a specific subset of proton channels that have voltage- and time-dependent gating like other ion channels. However, they differ from most ion channels in their extraordinarily high selectivity, tiny conductance, strong temperature and deuterium isotope effects on conductance and gating kinetics, and insensitivity to block by steric occlusion. Gating of H(+) channels is regulated tightly by pH and voltage, ensuring that they open only when the electrochemical gradient is outward. Thus they function to extrude acid from cells. H(+) channels are expressed in many cells. During the respiratory burst in phagocytes, H(+) current compensates for electron extrusion by NADPH oxidase. Most evidence indicates that the H(+) channel is not part of the NADPH oxidase complex, but rather is a distinct and as yet unidentified molecule.

Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.

Hydrodynamic model of temperature change in open ionic channels

Most theories of open ionic channels ignore heat generated by current flow, but that heat is known to be significant when analogous currents flow in semiconductors, so a generalization of the Poisson-Nernst-Planck theory of channels, called the hydrodynamic model, is needed. The hydrodynamic theory is a combination of the Poisson and Euler field equations of electrostatics and fluid dynamics, conservation laws that describe diffusive and convective flow of mass, heat, and charge (i.e., current), and their coupling. That is to say, it is a kinetic theory of solute and solvent flow, allowing heat and current flow as well, taking into account density changes, temperature changes, and electrical potential gradients. We integrate the equations with an essentially nonoscillatory shock-capturing numerical scheme previously shown to be stable and accurate. Our calculations show that 1) a significant amount of electrical energy is exchanged with the permeating ions; 2) the local temperature of the ions rises some tens of degrees, and this temperature rise significantly alters for ionic flux in a channel 25 A long, such as gramicidin-A; and 3) a critical parameter, called the saturation velocity, determines whether ionic motion is overdamped (Poisson-Nernst-Planck theory), is an intermediate regime (called the adiabatic approximation in semiconductor theory), or is altogether unrestricted (requiring the full hydrodynamic model). It seems that significant temperature changes are likely to accompany current flow in the open ionic channel.

Sodium in gramicidin: an example of a permion

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.

The double pi pi 5.6 helix of gramicidin A predominates in unsaturated lipid membranes

The structure of the channel-forming polypeptide gramicidin A (GA) incorporated into phosphatidyl-choline (PC) liposomes has been studied as a function of the degree of unsaturation of the acyl chains of PC. The initial conformational state of GA in reconstituted bilayers is determined by the solvent in which the peptide and the lipid are initially co-dissolved, whereas the equilibrium conformational state (after heat incubation) is affected by the lipid structure rather than by the nature of the solvent. The conformational equilibrium of GA has been studied in liposomes prepared from PC having a variable number of double bonds in the fatty acid moiety, by circular dichroism and Fourier transform infrared. Liposomes were prepared from trifluoroethanol or ethanol solutions and incubated at 68 degrees C. GA was shown to retain the conformation of the right-handed pi-->6.3 pi<--6.3 helix in PC with saturated acyl chains and with one double bond, whereas in dilinoleoyl-PC, having two double bond in each chain, the thermodynamically preferred structures are left-handed antiparallel and parallel double pi pi 5.6 helices. Natural soybean PC also favours left-handed pi pi 5.6 helical structures of GA (approximately 75%). This finding is discussed in terms of the role of PC unsaturation in the dynamic properties of the lipid matrix. Differences between observed FTIR spectra of the increases decreases pi pi 5.6 helix in solution (and to a larger extent in the membrane) and the calculated IR spectra can be interpreted as resulting from deviation of the real structure from the theoretically derived ideal helix. The data obtained provide grounds for better understanding of a GA channel functioning in lipids of variable degrees of unsaturation.

Ion-water and ion-polypeptide correlations in a gramicidin-like channel. A molecular dynamics study

This work describes a molecular dynamics study of ion-water and ion-polypeptide correlation in a model gramicidin-like channel (the polyglycine analogue) based upon interaction between polarizable, multipolar groups. The model suggests that the vicinity of the dimer junction and of the ethanolamine tail are regions of unusual flexibility. Cs+ binds weakly in the mouth of the channel: there it coordinates five water molecules and the #11CO group with which it interacts strongly and is ideally aligned. In the channel interior it is generally pentacoordinate; at the dimer junction, because of increased channel flexibility, it again becomes essentially hexacoordinate. The ion is also strongly coupled to the #13 CO but not to either #9 or #15, consistent with 13C NMR data. Water in the channel interior is strikingly different from bulk water; it has a much lower mean dipole moment. This correlates with our observation (which differs from that of previous studies) that water-water angular correlations do not persist within the channel, a result independent of ion occupancy or ionic polarity. In agreement with streaming potential measurements, there are seven single file water molecules associated with Cs+ permeation; one of these is always in direct contact with bulk water. At the mouth of an ion-free channel, there is a pattern of dipole moment alteration among the polar groups. Due to differential interaction with water, exo-carbonyls have unusually large dipole moments whereas those of the endo-carbonyls are low. The computed potential of mean force for CS+ translocation is qualitatively reasonable. However, it only exhibits a weakly articulated binding site and it does not quantitatively account for channel energetics. Correction for membrane polarization reduces, but does not eliminate, these problems.

Gramicidin A--phospholipid model systems

Gramicidin A forms ion-conducting channels which can traverse the hydrocarbon core of lipid bilayer membranes. The structures formed by gramicidin A are among the best characterized of all membrane-bound polypeptides or proteins. In this review a brief summary is given of the occurrence, conformation, and synthesis of gramicidin A, and of its use as a model for ion transport and the interaction of proteins and lipids in biological membranes.

Proton permeation of lipid bilayers

Proton permeation of the lipid bilayer barrier has two unique features. First, permeability coefficients measured at neutral pH ranges are six to seven orders of magnitude greater than expected from knowledge of other monovalent cations. Second, proton conductance across planar lipid bilayers varies at most by a factor of 10 when pH is varied from near 1 to near 11. Two mechanisms have been proposed to account for this anomalous behavior: proton conductance related to contaminants of lipid bilayers, and proton translocation along transient hydrogen-bonded chains (tHBC) of associated water molecules in the membrane. The weight of evidence suggests that trace contaminants may contribute to proton conductance across planar lipid membranes at certain pH ranges, but cannot account for the anomalous proton flux in liposome systems. Two new results will be reported here which were designed to test the tHBC model. These include measurements of relative proton/potassium permeability in the gramicidin channel, and plots of proton flux against the magnitude of pH gradients. (1) The relative permeabilities of protons and potassium through the gramicidin channel, which contains a single strand of hydrogen-bonded water molecules, were found to differ by at least four orders of magnitude when measured at neutral pH ranges. This result demonstrates that a hydrogen-bonded chain of water molecules can provide substantial discrimination between protons and other cations. It was also possible to calculate that if approximately 7% of bilayer water was present in a transient configuration similar to that of the gramicidin channel, it could account for the measured proton flux. (2) The plot of proton conductance against pH gradient across liposome membranes was superlinear, a result that is consistent with one of three alternative tHBC models for proton conductance described by Nagle elsewhere in this volume.

Simulation of voltage-driven hydrated cation transport through narrow transmembrane channels

Molecular dynamics studies for the voltage-driven transport of the alkali metal ions lithium, sodium, and potassium through gramicidin A-type channels filled with water molecules are presented. The number of water molecules in the channel is obtained from a previous study (Skerra, A., and J. Brickmann, 1987, Biophys. J., 51:969-976). It is shown that the selectivity of the intrachannel ion diffusion through our model pore conforms to the experimentally observed selectivity of the gramicidin A channel. It is demonstrated that the number of water molecules in the channel plays a key role for the selectivity.

Why is gramicidin valence selective? A theoretical study

Calculations contrasting the channel solvation energy for cesium ions and chloride ions associated with water in gramicidin-like channels are presented. The energy profile for the cation exhibits a deep well at the channel entrance; within the single file region the solvation energy is roughly constant. The anion exhibits a totally different energy profile. There is an energy barrier at the channel entrance; if the ion could surmount this barrier, it would be quite stable within the channel. At the channel entrance, the calculated solvation energy difference between anion and cation is approximately 15 kcal mol-1. This is completely consistent with the observation that chloride neither permeates nor blocks the channel since the estimated rate of ion entry would be approximately 0.01-10(-5) s-1, far slower than the rate at which the channel dimer dissociates into monomers.

Permeability of lipid bilayers to water and ionic solutes

The lipid bilayer moiety of biological membranes is considered to be the primary barrier to free diffusion of water and solutes. This conclusion arises from observations of lipid bilayer model membrane systems, which are generally less permeable than biological membranes. However, the nature of the permeability barrier remains unclear, particularly with respect to ionic solutes. For instance, anion permeability is significantly greater than cation permeability, and permeability to proton-hydroxide is orders of magnitude greater than to other monovalent inorganic ions. In this review, we first consider bilayer permeability to water and discuss proposed permeation mechanisms which involve transient defects arising from thermal fluctuations. We next consider whether such defects can account for ion permeation, including proton-hydroxide flux. We conclude that at least two varieties of transient defects are required to explain permeation of water and ionic solutes.