| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Dept. of Physiology, Loyola University Medical Center, Maywood, Illinois 60153
Correspondence: Address reprint requests to Samuel Cukierman, Department of Physiology, Loyola University Medical Center, 2160 South First Avenue, Maywood, IL 60153. Tel.: 708-216-9471; Fax: 708-216-6308; E-mail: scukier{at}lumc.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Unidimensional chains of water molecules have been found in restricted spaces in various proteins that participate in bioenergetic processes (Baciou and Michel, 1995
; Branden et al., 2001
; Luecke et al., 1999
; Sass et al., 2000
; Zaslavsky and Gennis, 1998
). In virtually all cells, the production of ATP is ultimately triggered by the movement of protons across a membrane protein. The complexity of those proteins associated with a tight functional coupling between proton transfer and redox potentials (Trumpower and Gennis, 1994
), and the fact that proton transfer cannot be directly measured in a single molecule, makes it extremely difficult to analyze the fine features of proton transfer in bioenergetic proteins. On the other hand, proton transfer in a relatively simple structure such as gramicidin A (gA) channels can be directly measured at the single molecular level, and insights can be gained into the basic rules that govern proton transfer in proteins.
gA is a highly hydrophobic pentadecapeptide secreted by Bacillus brevis. In lipid bilayers, its primary structure (HCO-L-Val-Gly-L-Ala-D-Leu-L-Ala-D-Val-L-Val-D-Val-(L-Trp-D-Leu)3-L-Trp-NH-(CH2)2-OH) defines a right-hand ß6.3 helix in which the side chain residues are in contact with the lipid environment, and the carbonyl and amide groups line the pore of the protein (Arseniev et al., 1985
; Ketchem et al., 1993
, 1997
; Sarges and Witkopf, 1965
; Urry, 1971
). The association via six intermolecular H-bonds between the amino termini of two gA peptides, each located in a distinct monolayer, results in the formation of a water-filled ion channel that is selective for monovalent cations (Andersen, 1984
; Koeppe and Andersen, 1996
; Hladky and Haydon, 1972
). Disruption of intermolecular H-bonds results in the dissociation of gA monomers with the consequent loss of ion channel function.
Four distinct and direct experimental evidences support the notion that the single-channel conductance to protons (gH) in gA channels is determined by proton transfer. 1) Proton permeation in gA channels is 12 orders of magnitude larger than the second most permeable ionic species (Cs+; Armstrong et al., 2001
; Busath and Szabo, 1988
; Eisenman et al., 1980
; Hladky and Haydon, 1972
; Myers and Haydon, 1972
; Quigley et al., 2000
). 2) Levitt et al. (1978)
and Finkelstein (1987)
demonstrated that the permeation of protons through native gA channels is not accompanied by water movement as with other monovalent cations. 3) The kinetic isotope effect for H+ transfer in native gA (Akeson and Deamer, 1991
), and in both the SS- and RR-dioxolane-linked gA dimers (Chernyshev et al., 2003
) is larger than the kinetic isotope effect for other monovalent cations (Tredgold and Jones, 1979
), and consistent with a proton transfer mechanism (Chernyshev et al., 2003
). 4) Activation energies for proton transfer in native gA and in the SS- and RR-dioxolane-linked gA dimers (Akeson and Deamer, 1991
; Chernyshev and Cukierman, 2002
) are in general less than with single-channel conductances to alkalines (Chernyshev and Cukierman, 2002
; DeCoursey and Cherny, 1998
) and consistent with those of proton transfer in bulk solution (Chernyshev and Cukierman, 2002
).
Two desformylated gA peptides have been covalently linked with various chemical groups: malonyl (Bamberg and Janko, 1977
; Urry et al., 1971
), glutaryl (Rudnev et al., 1981
), and the SS and RR diacid dioxolane (Cukierman et al., 1997
; Quigley et al., 1999
; Stankovic et al., 1989
). Those covalently linked gA peptides form ion channels in lipid bilayers with single-channel conductance properties similar to native gA channels. As expected, the average lifetime of covalently linked gA channels is considerably longer than in native gA. In our laboratory, the SS- and RR-dioxolane-linked dimers have been used as experimental models to probe structurefunction relationships of proton transfer in proteins (Armstrong et al., 2001
; Armstrong and Cukierman, 2002
; Chernyshev and Cukierman, 2002
; Chernyshev et al., 2003
; Cukierman et al., 1997
; Cukierman, 1999
; 2000
; Godoy and Cukierman, 2001
; Quigley et al., 1999
). One of the most challenging questions that we pose is illustrated in Fig. 1
.
|
|
2 M. The bottom panel shows that the relationship between proton conductivity in water and [H+] is a straight line (up to a [H+] of
1 M) with a slope of 0.98 if proton concentrations are used or 1.00 if those concentrations are corrected for their thermodynamic activities (Cukierman, 2000
In addressing the molecular origin of the brief closures in various gA channels, Armstrong and Cukierman (2002)
demonstrated that in 1 M HCl, gH in GMO/squalene bilayers was considerably larger than in GMO/decane bilayers for native gA, SS, and RR channels. The significant difference between these bilayers is that a GMO/squalene bilayer is
50% thinner than a GMO/decane bilayer (see below). This result was intriguing and insightfulintriguing, because previous measurements by various investigators have demonstrated that the single-channel conductances to alkalines were not affected by manipulations of membrane thickness via solvent effects (see below, and Hladky and Haydon, 1972
; Kolb and Bamberg, 1977
; Rudnev et al., 1981
); and insightful, because it offered an opportunity to reinvestigate gH-[H+] relationships in bilayers with various thicknesses with the aim of identifying the molecular factors that modulate said relationships.
Our objective in this study was to perform a set of systematic measurements of gH-[H+] relationships in native gA, SS, and RR channels in monoglyceride bilayers with various thicknesses. It is now demonstrated that gH-[H+] relationships in these various gA channels channel are distinctly modulated by the thickness of lipid bilayers. This modulation is considerably more pronounced and distinct in the SS channel in the [H+] range of 0.011 M, which corresponds to the shoulder region in log-log plots of gH vs. [H+] (Eisenman et al., 1980
; Gowen et al., 2002
). We demonstrate that gH values in the SS channel in that [H+] range converge to those of native gA channels in thin bilayers. The shape of gH-[H+] relationships in distinct gA channels is determined by lipid-protein interactions.
| MATERIALS AND METHODS |
|---|
|
|
|---|
. It has been shown (see below) that the thickness of monoglyceride bilayers depends on the solvent used to form the bilayer. In this study, the monoglycerides contained cis-mono-unsaturated fatty acid chains. The composition of the bilayer will be referred to as GM-Cx/solvent, where x is the number of carbons in the fatty acid chain, and the solvent is decane, hexadecane, or squalene. Monoglycerides were purchased from Nu-Check Prep (Elysian, MN): Monoerucin-GM-C22 (
13 cis-monodocosenoin), Monoeicosenoin-GM-C20 (
11 cis-monoeicosenoin), Monoolein-GM-C18 (
9 cis-monoolein), and Monopalmitolein-GM-C16 (
9 cis-monopalmitolein). Decane, hexadecane, and squalene were obtained from Sigma (St. Louis, MO), and twice purified in a column containing (from bottom to top) neutral, basic, and acid chromatographic grade alumina (Sigma). This purification step is particularly important for squalene (White, 1978
60 mg of monoglyceride in 1 mL of solvent.
The thicknesses of planar bilayers formed by a combination of monoglycerides and solvents were measured by various investigators using several techniques (Benz et al., 1975
; Dilger, 1981
; Dilger and Benz, 1985
; Elliot et al., 1983
; Lewis and Engelman, 1983
; Requena et al., 1975
; Rudnev et al., 1981
; White, 1978
). Table 1 shows the thicknesses of the bilayers used in this study. The thickness of these planar bilayers is essentially determined by the distribution of solvent in the bilayer core. Relatively long chain alkanes (hexadecane) are mostly distributed between fatty acid chains of the monoglycerides whereas alkanes with shorter chains (decane, for example) partition between the monolayers of the bilayer (McIntosh et al., 1980
; Simon et al., 1977
). On the other hand, squalene seems to be virtually excluded from the lipid bilayer (White, 1978
).
|
) of various glycerides at the air/water interface was measured with the Wilhelmy plate method using a tensiometer (KSV Instruments, Helsinki, Finland). A droplet of lipid solution in chloroform was deposited on the water surface, and 1520 min were allowed for the chloroform to evaporate completely. These measurements were obtained at room temperature (24°C), and are also reported in Table 1.
Channels
The gA peptides were obtained from Fluka (Milwaukee, WI). The SS and RR stereoisomers of dioxolane-linked gA channels were synthesized, purified, and characterized as previously described (Cukierman et al., 1997
; Quigley et al., 1999
; Stankovic et al., 1989
).
Solutions
Most experiments were performed with symmetrical solutions of HCl of various concentrations across the bilayer. Some experiments were also performed in 1 or 0.25 M CsCl. All experiments were performed at room temperature (2224°C).
Electrical measurements
The planar bilayer was voltage-clamped at voltages where gH has an ohmic behavior, and single-channel currents were measured using a pair of 0.5 mm Ag/AgCl silver wires immersed in distinct solutions across the lipid bilayer. For a given experimental condition (one gA channel in one type of bilayer in a given solution), at least seven measurements of single channels were obtained from at least two (usually more than four) distinct planar bilayers. Single-channel conductances were measured in sectors of the recording where flickers were absent or occurred at a very low frequency. Single-channel conductances are expressed as mean ± SE.
Analysis
Single-channel recordings were digitized and analyzed using pClamp software (Axon Instruments, Union City, CA). Statistical analysis and graphs were done with Sigmaplot 2000 (SPSS Science, Chicago, IL).
| RESULTS |
|---|
|
|
|---|
|
|
|
Let us first concentrate on the open symbols of Fig. 5 (circles, native gA channels; squares, SS channels) and at [H+] < 2 M. In bilayers whose thicknesses were between 54 and 42 Å, the log[gH]-log[H+] relationships were linear for the SS channel in [H+]
2 M. Notice that the slopes of these linear relationships are slightly different (0.660.85; see ) in various bilayers. Perhaps the most interesting experimental finding in this article concerns the fact that in monoglyceride bilayers whose thicknesses are smaller than 37 Å, the relationships between gH and [H+] in the SS channel became essentially indistinguishable from those in native gA channels in loglog plots. In the thicker bilayers, loglog relationships became linear in the SS channel.
By contrast, the log(gH)-log([H+]) relationships for native gA channels display the usual shoulder-like appearance (Eisenman et al., 1980
; Gowen et al., 2002
) as previously mentioned in the Introduction. This shape was independent of the thickness of monoglyceride bilayers (see also Fig. 7, top panel). In thick (>37 Å) monoglyceride bilayers, gH values in native gA channels at [H+] < 1 M were significantly smaller than in SS channels. At [H+] > 1 M, however, gH values were larger in native gA channels.
|
1 M (see Fig. 3). This point is further documented in Fig. 6. In this figure, the ratios between gH values in GM-C18/decane (48 Å) and GM-C18/hexadecane (32 Å) for the SS (squares), and native gA (circles) and RR channels (triangles) were plotted as a function of [H+]. In native gA and RR channels, the ratios between gH values remained relatively constant or slightly less than 1 at various [H+] (see Figs. 2, 6, and 8). By contrast, these ratios could be twice as large in the SS channel (see also Fig. 3). Consequently, the convergence between the gH values of the SS and native gA channels in [gH]-[H+] plots in thin GM-C18/hexadecane bilayers (middle panel, right column of Fig. 4) is a consequence of a major attenuation of gH in the SS channel in thin bilayers at [H+] < 1 M.
|
|
; GM-C16/hexadecane,
). By contrast, gHmax values for the SS channel were measured in either GM-C18/decane () or GM-C16/decane (
) bilayers. Notice that for the SS channel, the solid symbols (thick bilayers) in Fig. 7 are consistently larger than the open (thin bilayers) symbols in various [H+].
, GM-C22/hexadecane,
). For the SS channel, however, gHmin values were systematically obtained in GM-C20/hexadecane (
, see Fig. 9 below and Discussion), followed by measurements in the thinnest bilayers (GM-C16/hexadecane,
; GM-C18/hexadecane,
; or GM-C18/squalene,
). Proton transfer in the SS channel in relation to native gA and RR channels is differentially modulated by bilayer thickness.
|
37 Å. In gA and RR channels, changes in gH were smaller, but definitely in the opposite direction of the SS. Significant attenuations of gH in GM-C20 (with either decane or hexadecane) bilayers for the SS channels, and in GM-C20/hexadecane for gA channels were measured (Fig. 8). It is likely that GM-C20 has a thickness-independent effect on the conformation of native gA and SS channels that disfavors proton transfer. At this point not much can be added to this observation.
It has been shown that within the range of 1 M < [H+] < 5 M, gH saturates in the SS channel (Cukierman et al., 1997
; Cukierman, 2000
; Godoy and Cukierman, 2001
). At these high [H+] levels, these observations also apply to bilayers with various thicknesses (Figs. 5 and 7). For the SS channel in various bilayers, gH saturated at values less than 2000 pS. For native gA channels, however, maximum gH values were close to 3000 pS, and saturated at larger [H+] than for the SS channel (Figs. 5 and 7).
So far, the thickness of the bilayer has been considered as the independent variable, and one of the differences between thick and thin bilayers concerned the solvent used (decane for thick and hexadecane for thin bilayers). The molecular environment between the monolayers in a bilayer changes with the solvent (McIntosh et al., 1980
; Simon et al., 1977
; White, 1978
). Thus, an alternative interpretation could be that the solvent itself was modulating proton transfer in SS channels. This possibility was addressed by measuring gH values in GM-C22/hexadecane bilayers (49 Å), and comparing them to those obtained in GM-C18/decane (48 Å). Another set of experiments was performed with GM-C18/squalene (25 Å), in which case the experimental measurements of gH were compared to those obtained with GM-C16/hexadecane (28 Å). These four sets of measurements were plotted as solid circles (native gA channels) and squares (SS channels) in the middle graph on the left column, and bottom graph on the right column of Fig. 5, respectively.
In the middle panel of the left column in Fig. 5, the relationship between log(gH) and log([H+]) for [H+]
2 M could still be represented by a straight line for the SS channel in GM-C22/hexadecane bilayers (solid squares). However, that straight line had a slope slightly larger (0.85) than in GM-C18/decane (0.75) due to a decreased gH at [H+]
0.1 M. Interestingly, this observation also applied to native gA channels. For native gA channels, there was a good agreement between gH values in GM-C18/decane and GM-C22/hexadecane bilayers. Most importantly. however, is the fact that gH values were still larger in the SS than in native gA channels. As for the comparison between gH-[H+] relationships in GM-C18/squalene and GM-C16/hexadecane bilayers (bottom panel on right column in Fig. 5), there was excellent agreement between the open and solid symbols (for both channels) in the log-log plots. It seems that the modulation of proton transfer in gA channels is a consequence of bilayer thickness.
Surface tension measurements at air/water interface were performed in the various lipids used in this study. No one-one relationships were found between gH (for the various gA) and surface tension measurements (see Table 1).
The single-channel conductances to several alkalines in native gA channels seem to be independent of the thickness of the bilayer as determined by the use of various solvents (Hladky and Haydon, 1972
; Kolb and Bamberg, 1977
; Rudnev et al., 1981
). These results contrast with gH measurements shown above for the SS and to a more limited extent in native gA and RR channels (Figs. 57). In view of the conceptual importance of this issue, gCs values were measured on both the SS and native gA channels in either GM-C18/decane or GM-C18/hexadecane bilayers. The data in Table 2 shows that bilayer thickness (as determined by the solvent) does not have a pronounced effect on gCs on both types of channels. In 0.25 M CsCl, the ratio between gCs values in the SS channel in thick and thin bilayers is
0.9. By contrast, the same ratio in 0.25 M HCl solutions is
2.0 (Figs. 5 and 6).
|
A unique observation was documented in the thickest bilayers used in this study (GM-C20/decane). This is shown in Fig. 9. In GM-C20/decane bilayers, it was common for both the SS and RR channels to display bursts of openings instead of very long openings (compare recordings in Figs. 3 and 9). In contrast to Fig. 3 in which the SS channel was open (with flickers) for more than 9 min, Fig. 9 shows that in GM-C20/decane bilayers, the channel opens in a bursting behavior (similar to native gA channels).
| DISCUSSION |
|---|
|
|
|---|
The hydrophobic length of gA channels is
22 Å (Elliot et al., 1983
). The function of gA channels (and of any membrane protein) depends on the appropriate interactions between the hydrophobic portions of gA and the membrane. In order for a gA channel to be functional in a bilayer with a hydrophobic thickness larger than 22 Å, the bilayer around the openings of a gA channel must adopt a conformation in which the openings of the channels are exposed to the outside solutions at the same time that the hydrophobic side chain residues of gA are shielded by the core of the bilayer (de Planque et al., 1998
; Hendry et al., 1978
; Hladky and Haydon, 1972
; Helfrich and Jakobsson, 1990
; Huang, 1986
; Elliot et al., 1983
; Killian et al., 1998
; Kolb and Bamberg, 1977
; Lundbæk and Andersen, 1994
; Lundbæk et al., 1996
; Ring, 1986
; van der Wel et al., 2000
). Lipid membranes are dynamical structures, and collective motions or undulations of planar lipid bilayers have been recently measured (Bayerl, 2000
; Hirn et al., 1998
). It is conceivable that the functionality of a native gA channel is a consequence of bilayer undulations. Two gA monomers located in opposite monolayers of a membrane will form an ion channel when thermal fluctuations in the two monolayers of the membrane bring their amino termini close enough to establish intermolecular H-bonds (see Helfrich and Jakobsson, 1990
, and their Fig. 6 in particular). Consequently, "deformations" of the bilayer around gA channels are not static structures but the result of intrinsic dynamic fluctuations in the lipid bilayer thickness (Armstrong and Cukierman, 2002
). In particular, it was proposed that undulations of the lipid bilayer underlie the fast closing flickers of various gA channels (Armstrong and Cukierman, 2002
).
In a previous study, Armstrong and Cukierman (2002)
have shown that proton transfer in native gA, SS, and RR channels was significantly larger in 1 M HCl in GM-C18/squalene bilayers than in GM-C18/decane bilayers. Considering that the resistance to proton transfer in the channel is comparable to that measured in bulk solution (Cukierman, 1999
, 2000
), it was reasoned that the lower gH in thick bilayers is a consequence of the presence of significant access resistances in series with the channel. These resistances would derive from long (proportional to the difference between the thickness of bilayer and channel's length) and diffusion-restricted spaces adjacent to the mouths of the pore in thick bilayers (Armstrong and Cukierman, 2002
). Although this hypothesis must be taken into consideration in attempting to explain the present experimental results, it is clear that the effects of bilayer thickness (or lipid/protein interactions) on the modulation of proton transfer in various gA channels are more complex than previously anticipated.
Proton transfer in monoglyceride bilayers
Experiments in various gA channels in [H+] > 1 M
Previous results (Cukierman, 2000
; Godoy and Cukierman, 2001
) and Figs. 5 and 7 in this article demonstrated that in [H+] > 1 M, gH saturates. It is likely that this saturation is a consequence of diffusion limitation in bulk solution (see bottom graph in Fig. 1; also see Cukierman, 2000
) at high [HCl]. As [HCl] increases, the mobility of protons in solution becomes increasingly determined by the hydrodynamic mobility of solvated H+ instead of H+ transfer (Agmon, 1998
; Cukierman, 2000
; Lengyel et al., 1962
; Owen and Sweeton, 1941
). This would attenuate the mobility of protons in solutions with the effect of limiting the access (exit) of protons to (from) the channel from (to) solution in concentrated acid solutions. An additional observation consistent with the mechanism described above concerns the behavior of gH in positively charged phospholipid bilayers (Godoy and Cukierman, 2001
). In those bilayers, saturation of gH values corrected for proton concentration at the membrane/solution interface occurred at [H+]
0.2 M which corresponds to a [H+]bulk of
2 M. In neutral monoglyceride bilayers in 0.2 M HCl solutions (Figs. 5 and 7), gH is still increasing with [H+]bulk. In relatively high [HCl], the effects of monoglyceride bilayer thickness on gH for the various gA channels were the same as in low [HCl] (see below).
Native gA channels
Small or no changes in the single-channel conductances to alkalines in native gA channels were measured in bilayers whose thicknesses were modified by the use of distinct solvents with the same monoglyceride (Hladky and Haydon, 1972
; Kolb and Bamberg, 1977
). This result has been confirmed (Table 2) for gCs in GM-C18 with either decane or hexadecane as the solvent. However, when the bilayer thickness was increased by using monoglycerides with longer fatty acid chains, a clear reduction in gCs, gNa, and gK (Hladky and Haydon, 1972
; Kolb and Bamberg, 1977
; Rudnev et al., 1981
) occurred in the respective 1 M chloride solutions. These experimental results are in a sense contradictory (see below).
Whereas the shape of the log(gH)-log([H+]) plots in native gA channels was maintained in the various monoglyceride bilayers (Fig. 7, top panel), at a given
and
had a relatively modest (compared to SS channels), but systematic, dependence on bilayer thickness. These extreme gH values were always measured in thin and thick monoglyceride bilayers, respectively. Contrary to previous results with alkalines (Hladky and Haydon, 1972
; Kolb and Bamberg, 1977
; Rudnev et al., 1981
; also see Table 2) there was a clear increase in proton transfer in bilayers made thinner by using hexadecane or squalene as substitutes for decane with the same monoglyceride (see Figs. 2, 6, and 8, for specific examples). On the other hand, in GM-C20 bilayers, significant attenuations of gH were measured. This is in agreement with results with alkalines (Kolb and Bamberg, 1977
; Rudnev et al., 1981
). In summary, in relatively thick bilayers (GM-C20/decane and GM-C22/hexadecane, our results; GM-C20/hexadecane, Kolb and Bamberg, 1977
; and GM-C22/heptane, Rudnev et al., 1981
) there is a general decrease in the ion-transport properties in gA channels, whereas in thin bilayers (with shorter fatty acid chains and/or by using hexadecane or squalene as solvents) there is a small but consistent increase of H+ transfer only (Figs. 2, 57, and Table 2).
We now examine a few distinct possibilities (still at the qualitative level) that could account for the apparently contradictory effects of monoglyceride bilayer thickness on the single-channel conductance to protons or alkalines in native gA channels.
One possibility concerns the development of an access resistance adjacent to the mouths of the channel in thick bilayers (Armstrong and Cukierman, 2002
). As discussed above, this would have the effect of attenuating gH by increasing the access resistance to protons in thick bilayers. At relatively high alkaline concentration, the rate-limiting step for alkaline permeation resides in the channel itself (Andersen, 1983
). Thus, and contrary to H+ transfer in the channel, this access resistance is not likely to contribute significantly to the total single-channel conductance in thick or thin bilayers. Consequently, this would explain the lack of effect of bilayer thickness on the single-channel conductance to alkalines but not to gH in native gA channels. It may well be that at low concentration of alkalines, in which diffusion limitation or access resistance plays a significant role and limits single-channel conductance (Andersen, 1983
), the single-channel conductance is larger in thin than in thick bilayers. Unfortunately, the small signal-to-noise ratio for the single-channel conductances in these experimental conditions makes the investigation of this phenomenon extremely difficult.
The general decrease in the transport properties of gA channels in the thickest bilayers may result from significant alterations of the channel's structure. Therefore, a second consideration is that the overall decline of single-channel conductances in thick bilayers could reflect an increased probability of mismatch between the hydrophobic portions of channel and the thick bilayer. In thick bilayers, one or both openings of gA channels could be obliterated by the bilayer. Undulations of thick lipid bilayers could be wide enough to obliterate the openings of the channel for longer times than in thin bilayers (Armstrong and Cukierman, 2002
; Bayerl, 2000
; Hirn et al., 1998
). Thus, the single-channel conductance to any cation in principle would represent an average value that reflects the relative times the channel spends in those conductive and nonconductive states. The fact that bursts of activity are present in dioxolane-linked gA channels (Fig. 9) provides some support to such an hypothesis. In bilayers whose thickness is
49 Å, the SS channel remains in the open state (with flickers) for tens of minutes or even hours (Cukierman et al., 1997
; Quigley et al., 1999
). However, in 54 Å bilayers, the SS (or RR) channels display the bursting pattern exemplified in Fig. 9. It seems reasonable to hypothesize that the intervals between bursts are a consequence of the electrical insulation of the mouths of covalently linked gA dimers from one or both solutions facing the mouths of the pore.
A third consideration is that significant fluctuations in the distance between the two gA monomers could occur in thick bilayers. These fluctuations would be a consequence of mobilities of gA monomers in monolayers in a direction perpendicular to the plane of the bilayer, and also to undulations of the membrane. Some of these gA conformations (including monomer dissociation) would disfavor ion permeation or proton transfer in the water wire in the channel, and as a consequence, single-channel conductances would be attenuated in general. In thin bilayers, fluctuations in gA intermonomeric distance would be more constrained favoring interactions that facilitate ion permeation or proton transfer. Because the transfer of protons in gA channels occurs with a time scale of 12 orders of magnitude faster than the permeation of alkalines (Grotthuss mechanism vs. hydrodynamic diffusion of water and alkalines), it is possible that such fluctuations would be more influential on proton transfer than on alkaline permeation. This could also explain why bilayers thinner than GM-C18/decane bilayers have a small or no effect on alkaline permeation but a significant effect on proton transfer.
RR- and SS-dioxolane-linked gA channels
The influence of monoglyceride bilayer thickness on proton transfer in RR channels is qualitatively the same as in native gA channels: the shoulder shape of the log-log plots is maintained in monoglyceride bilayers with various thicknesses, and gH varies inversely with bilayer thickness. Thus, most of the discussion above in the context of native gA channels could also apply to RR channels. It should be remarked, though, that as in our previous studies, gH values in the RR channel are systematically and significantly lower than in gA (or SS) channels.
Perhaps the most interesting and completely unexpected observation in this study was that gH in SS channels is considerably larger in thick (>37 Å) than in thin monoglyceride bilayers. This occurs more dramatically in [H+] < 1 M, and causes the linearization (slope <1) in the log(gH)-log[H+] plots. This effect was not measured with gCs (Table 2), and provides an interesting contrast to the experimental observations with the RR channel. Consequently, the enhancement of gH in thick bilayers cannot be related to the presence of the dioxolane linker per se.
In view of some recent findings in other laboratories and our experimental results in this study, it is of interest to address some experimental factors that regulate the linear vs. nonlinear behavior in log-log plots of gH vs. [H+] in gA channels.
Cukierman (2000)
demonstrated that not only gH values in the SS and RR channels have different magnitudes, but the shapes of the relationships between gH and [H+] were also distinct (Fig. 1). The shape of the log-log relationship between gH and [H+] in RR channels was similar to native gA channels (Eisenman et al., 1980
; Gowen et al., 2002
). This can be interpreted as if two (or more) H+ can simultaneously occupy the pore of RR channels (Hille and Schwarz, 1978
). However, this idea is not immediately evident or applicable for the SS channel, which has a linear (with slopes considerably less than 1.00) log-log relationship in thick monoglyceride bilayers over a wide range of proton concentrations. We have performed calculations on the dependence of gH on [H+] (results not shown) using Eyring rate theory (Eyring et al., 1949
; Hille and Schwarz, 1978
). The model for proton permeation inside the channel consisted of a sequence of three energy peaks separated by two wells. Allowance was given for double occupancy of the pore by protons. Calculations were performed with the energy peaks and valleys located at various electrical distances inside the channel. Using this model, it was not difficult to obtain log(gH)-log([H+]) linear relationships with a slope of 1.00 (see, for example, Fig. 3 in Hille and Schwarz, 1978
). However, linear relationships with slopes less than 1.00 in the log-log plots were not obtained. We have previously argued that the linear log-log plot in the SS channel could be a consequence of a rate-limiting step for H+ transfer in the solution/channel interfaces (Cukierman, 2000
) or inside the channel itself (Godoy and Cukierman, 2001
) but not to diffusion in bulk solution. We suspect that these distinct phenomena have different relative contributions to gH at various [H+]. Although the cause for the linearity (slope < 1) vs. nonlinearity of H+ transfer in gA channels in log-log plots remains unknown, in this study we demonstrated that the thickness of monoglyceride bilayers is definitely one parameter involved in this phenomenon in the SS channel.
Since our original observations (Fig. 1, Cukierman, 2000
), two recent publications have demonstrated that native gA channels can also display linear loglog relationships in gH-[H+] plots. In one of them, Rotikskaya et al. (2002)
showed that this happens in diphytanoylphosphatidylcholine/decane bilayers (slope = 0.51). We have preliminarily confirmed their results, albeit with a slope of
0.60 (Chernyshev and Cukierman, work in progress).
Gowen et al. (2002)
demonstrated that the replacement of the four Trp in gA channels by Phe (gM channels) caused linearization (slope = 1) of the log(gH)-log([H+]) plot. This linearization in gM channels is a consequence of a decrease and increase in gH within the [H+] range of 250 mM and 1001000 mM, respectively. In this study, the analysis of single-channel H+ currents was restricted to the low domain of [H+], and the marked changes in the gH of gM channels were attributed to a combination of three factors (gM relative to gA channels): 1) an increase in the energy of H+ inside the channel by
3 kcal/mol, which has the effect of causing a substantial decrease in the probability of finding the channel occupied simultaneously by two protons over a wide range of [H+]; 2) the exit rate of H+ is sped up 65-fold; and 3) the entrance rate of H+ in the channel is voltage-dependent in the gM channel (but not in gA), and threefold slower in gM channels. These alterations were discussed in terms of electrostatic interactions between protons and the indol groups of the Trp residues (Gowen et al., 2002
).
The thermodynamic and kinetic factors that may account for the linearization of log(gH)-log([H+]) plots in gM channels under monooccupancy conditions (Gowen et al., 2002
) do not seem to apply to a similar phenomenon caused by thick monoglyceride membranes in the SS channel: 1) the slopes of log-log plots for the SS channel are significantly smaller than 1 in thick bilayers (see also Rotikskaya et al., 2002
); 2) thick monoglyceride bilayers cause an overall increase of gH at [H+] < 1 M, an observation which contrasts with those in gM and gA channels (see above); and 3) it does not seem likely that different electrostatic interactions occur between protons and the Trp residues in native gA, SS, and RR channels. Nevertheless, a clear differential effect of the bilayer thickness on proton transfer in [H+] < 1 M was demonstrated for the SS channel.
The most prominent structural difference between the SS and RR channels is a distortion localized in the middle of the RR channel caused by an approximate 90° tilt of the RR in relation to the SS dioxolane (Quigley et al., 1999
; Stankovic et al., 1989
; Yu et al., 2003
). Considering this localized structural difference, distinct interactions in the middle of the SS and RR channels with thick monoglyceride bilayers could in principle account for differences in H+ transfer between these channels. Evidently, the open question concerns the atomic details of this phenomenon.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on July 1, 2002; accepted for publication August 29, 2002.
| REFERENCES |
|---|
|
|
|---|
Agmon, N. 1998. Structure of concentrated HCl solutions. J. Phys. Chem. 102:192199.
Akeson, M., and D. W. Deamer. 1991. Proton conductance by the gramicidin water wire. Model for proton conductance in the F0F1ATPases? Biophys. J. 60:101109.
Andersen, O. S. 1983. Ion movement through gramicidin A channels. Single channel measurements at very high potentials. Biophys. J. 41:119133.
Andersen, O. S. 1984. Gramicidin channels. Annu. Rev. Physiol. 46:531548.[Medline]
Armstrong, K. M., and S. Cukierman. 2002. On the origin of closing flickers in gramicidin channels: a new hypothesis. Biophys. J. 82:13291337.
Armstrong, K. M., E. P. Quigley, P. Quigley, D. S. Crumrine, and S. Cukierman. 2001. Covalently linked gramicidin channels: effects of linker hydrophobicity and alkaline metals on different stereoisomers. Biophys. J. 80:18101818.
Arseniev, A. S., I. L. Barsukov, V. F. Bystrov, A. L. Lonize, and Y. A. Ovchinikov. 1985. Proton NMR study of gramicidin A transmembrane ion channel: head-to-head right-handed, single-stranded helices. FEBS Lett. 186:168174.[Medline]
Baciou, L., and H. Michel. 1995. Interruption of the water chain in the reaction center from r. sphaeroides reduces the rate of the proton uptake and of the second electron transfer to Q_B. Biochemistry. 34:79677972.[Medline]
Bamberg, E., and K. Janko. 1977. The action of a carbonsuboxide dimerized gramicidin A on lipid bilayer membranes. Biochem. Biophys. Acta. 465:486499.[Medline]
Bayerl, T. M. 2000. Collective membrane motions. Curr. Opinion Coll. Interf. Sci. 5:232236.
Benz, R., O. Frohlich, P. Lauger, and M. Montal. 1975. Electrical capacity of black lipid films and of lipid bilayers made from monolayers. Biochem. Biophys. Acta. 394:323334.[Medline]
Bernal, J. D., and R. H. Fowler. 1933. A theory of water and ionic solution, with particular reference to hydrogen and hydroxyl ions. J. Chem. Phys. 1:515548.
Branden, M., H. Sigurdson, A. Namslauer, R. B. Gennis, P. Adelroth, and P. Brzezinski. 2001. On the role of the K-proton transfer pathway in cytochrome c oxidase. Proc. Natl. Acad. Sci. USA. 98:50135018.
Busath, D. D., and G. Szabo. 1988. Permeation and characteristics of gramicidin conformers. Biophys. J. 53:697707.
Chernyshev, A., and S. Cukierman. 2002. Thermodynamic view of activation energies of proton transfer in various gramicidin A channels. Biophys. J. 82:182192.
Chernyshev, A., R. Pomès, and S. Cukierman. 2003. Kinetic isotope effects of proton transfer in aqueous and methanol containing solutions, and in gramicidin channels. Biophys. Chem. In press.
Conway, B. E., J. O. M. Bockris, and H. Linton. 1956. Proton conductance and the existence of the H3O+ ion. J. Chem. Phys. 24:834852.
Cukierman, S. 1999. Flying protons in linked gramicidin A channels. Isr. J. Chem. 39:419426.
Cukierman, S. 2000. Proton mobilities in water and in different stereoisomers of covalently linked gramicidin A channels. Biophys. J. 78:18251834.
Cukierman, S., E. P. Quigley, and D. S. Crumrine. 1997. Proton Conduction in gramicidin A and in its dioxolane-linked dimer in different lipid bilayers. Biophys. J. 73:24892502.
Danneel, V. H. 1905. Notiz über ionengeschwindigkeiten. Z. Elektrochem. 11:249252.
Day, T. J. F., U. W. Schmitt, and G. A. Voth. 2000. The mechanism of hydrated proton transfer in water. J. Am. Chem. Soc. 122:1202712028.
De Planque, M. R. R., D. V. Greathouse, R. E. Koeppe II, H. Schäfer, D. Marsh, and J. A. Killian. 1998. Influence of lipid/peptide hydrophobic mismatch on the thickness of diacylphosphatidylcholine bilayers. A 2H NMR and ESR using designed transmembrane
-helical peptides and gramicidin A. Biochemistry. 37:93339345.[Medline]
DeCoursey, T. E., and V. V. Cherny. 1994. Voltage-activated hydrogen ion currents. J. Membr. Biol. 141:203223.[Medline]
DeCoursey, T. E., and V. V. Cherny. 1998. Temperature dependence of voltage-gated H+ currents in human neutrophils, rat alveolar epithelial cells, and mammalian phagocytes. J. Gen. Physiol. 112:503522.
Dilger, J. P. 1981. The thickness of monoolein lipid bilayers as determined from reflectance measurements. Biochim. Biophys. Acta. 645:357363.[Medline]
Dilger, J. P., and R. Benz. 1985. Optical and electrical properties of thin monoolein bilayers. 1985. J. Membr. Biol. 85:181189.[Medline]
Eisenman, G., B. Enos, J. Hägglund, and J. Sandblom. 1980. Gramicidin A as an example of a single filing ionic channel. Ann. N. Y. Acad. Sci. 329:820.
Elliot, J. R., D. Neddham, J. P. Dilger, and D. A. Haydon. 1983. The effects of bilayer thickness and tension on gramicidin single-channel lifetime. Biochim. Biophys. Acta. 735:95103.[Medline]
Eyring, H., R. Lumry, and J. W. Woodbury. 1949. Some applications of modern rate theory to physiological systems. Rec. Chem. Prog. 10:100114.
Finkelstein, A. 1987. Water Movement Through Lipid Bilayers, Pores, and Plasma Membrane. Theory and Reality. John Wiley & Sons, New York.
Godoy, C. M. G., and S. Cukierman. 2001. Modulation of proton transfer in the water wire of dioxolane-linked gramicidin channels by lipid membranes. Biophys. J. 81:14301438.
Gowen, J. A., J. C. Markham, S. E. Morrison, T. A. Cross, D. D. Busath, E. J. Mapes, and M. F. Schumaker. 2002. The role of Trp side chains in tuning single proton conduction through gramicidin channels. Biophys. J. 83:880898.
Helfrich, P., and E. Jakobsson. 1990. Calculation of deformation energies and conformations in lipid membranes con