Covalently linked gramicidin channels: effects of linker hydrophobicity and alkaline metals on different stereoisomers

The direct role of the dioxolane group on the gating and single-channel conductance of different stereoisomers of the dioxolane-linked gramicidin A (gA) channels reconstituted in planar lipid bilayers was investigated. Four different covalently linked gA dimers were synthesized. In two of them, the linker was the conventional dioxolane described previously (SS and RR channels). Two gAs were covalently linked with a novel modified dioxolane group containing a retinal attachment (ret-SS and ret-RR gA dimers). These proteins also formed ion channels in lipid bilayers and were selective for monovalent cations. The presence of the bulky and hydrophobic retinal group immobilizes the dioxolane linker in the bilayer core preventing its rotation into the hydrophilic lumen of the pore. In 1 M HCl the gating kinetics of the SS or RR dimers were indistinguishable from their retinal counterparts; the dwell-time distributions of the open and closed states in the SS and ret-SS were basically the same. In particular, the inactivation of the RR was not prevented by the presence of the retinal group. It is concluded that neither the fast closing events in the SS or RR dimers nor the inactivation of the RR are likely to be a functional consequence of the flipping of the dioxolane inside the pore of the channel. On the other hand, the inactivation of the RR dimer was entirely eliminated when alkaline metals (Cs(+) or K(+)) were the permeating cations in the channel. In fact, the open state of the RR channel became extremely stable, and the gating characteristics of both the SS and RR channels were different from what was seen before with permeating protons. As in HCl, the presence of a retinal in the dioxolane linker did not affect the gating behavior of the SS and RR in Cs(+)- or K(+)-containing solutions. Alternative hypotheses concerning the gating of linked gA dimers are discussed.     

Proton conduction in gramicidin A and in its dioxolane-linked dimer in different lipid bilayers.

Gramicidin A (gA) molecules were covalently linked with a dioxolane ring. Dioxolane-linked gA dimers formed ion channels, selective for monovalent cations, in planar lipid bilayers. The main goal of this study was to compare the functional single ion channel properties of natural gA and its covalently linked dimer in two different lipid bilayers and HCl concentrations (10-8000 mM). Two ion channels with different gating and conductance properties were identified in bilayers from the product of dimerization reaction. The most commonly observed and most stable gramicidin A dimer is the main object of this study. This gramicidin dimer remained in the open state most of the time, with brief closing flickers (tau(closed) approximately 30 micros). The frequency of closing flickers increased with transmembrane potential, making the mean open time moderately voltage dependent (tau(open) changed approximately 1.43-fold/100 mV). Such gating behavior is markedly different from what is seen in natural gA channels. In PEPC (phosphatidylethanolamine-phosphatidylcholine) bilayers, single-channel current-voltage relationships had an ohmic behavior at low voltages, and a marked sublinearity at relatively higher voltages. This behavior contrasts with what was previously described in GMO (glycerylmonooleate) bilayers. In PEPC bilayers, the linear conductance of single-channel proton currents at different proton concentrations was essentially the same for both natural and gA dimers. g(max) and K(D), obtained from fitting experimental points to a Langmuir adsorption isotherm, were approximately 1500 pS and 300 mM, respectively, for both the natural gA and its dimer. In GMO bilayers, however, proton affinities of gA and the dioxolane-dimer were significantly lower (K(D) of approximately 1 and 1.5 M, respectively), and the g(max) higher (approximately 1750 and 2150 pS, respectively) than in PEPC bilayers. Furthermore, the relationship between single-channel conductance and proton concentration was linear at low bulk concentrations of H+ (0.01-2 M) and saturated at concentrations of more than 3 M. It is concluded that 1) The mobility of protons in gramicidin A channels in different lipid bilayers is remarkably similar to proton mobilities in aqueous solutions. In particular, at high concentrations of HCl, proton mobilities in gramicidin A channel and in solution differ by only 25%. 2) Differences between proton conductances in gramicidin A channels in GMO and PEPC cannot be explained by surface charge effects on PEPC membranes. It is proposed that protonated phospholipids adjacent to the mouth of the pore act as an additional source of protons for conduction through gA channels in relation to GMO bilayers. 3) Some experimental results cannot be reconciled with simple alterations in access resistance to proton flow in gA channels. Said differences could be explained if the structure and/or dynamics of water molecules inside gramicidin A channels is modulated by the lipid environment and by modifications in the structure of gA channels. 4) The dioxolane ring is probably responsible for the closing flickers seen in the dimer channel. However, other factors can also influence closing flickers.     

On the origin of closing flickers in gramicidin channels: a new hypothesis

The submillisecond closing events (flickers) and the single channel conductances to protons (g(H)) were studied in native gramicidin A (gA) and in the SS and RR diastereoisomers of dioxolane-linked gA channels in planar bilayers. Bilayers were formed from glycerylmonooleate (GMO) in various solvents. In GMO/decane (thick) bilayers, the largest flicker frequency occurred in the SS channel (39 s(-1)), followed by the RR (4 s(-1)) and native gA channels (3 s(-1)). These frequencies were attenuated in GMO/squalene (thin) bilayers by 100-, 30-, and 70-fold in the SS, RR, and native gA channels, respectively. In thin bilayers, the average burst duration of native gA channels was 30-fold longer than in thick bilayers. The RR dioxolane-linked gA dimer "inactivated" in GMO/decane but not in squalene-containing bilayers. The mean closed time of flickers (approximately 0.12 ms) was essentially the same in various gA channels. In thin bilayers, g(H) values were larger by approximately 10% (SS), 30% (RR), and 20% (native gA) in relation to thick bilayers. It is concluded that flickers are not related to pre-dissociation or dissociation states of gA monomers, and do not seem to be caused by intrinsic conformational changes of channel proteins. It is proposed that flickers are caused by undulations of the bilayer that obliterate the openings of gA channels. Differences between flicker frequencies in various gA channels are likely to result from variations in channel geometries at the bilayer/channel interface. The smaller g(H) in thick bilayers suggests that the deformation of these bilayers around the gA channel creates a diffusional pathway next to the mouths of the channel that is longer and more restrictive than in thin GMO bilayers. A possible molecular interpretation for these effects is attempted.     

Proton transfer in gramicidin channels is modulated by the thickness of monoglyceride bilayers

The thickness of monoglyceride planar bilayers has significant effects on the transfer of protons in both native gramicidin A (gA) and in covalently linked SS- and RR-dioxolane-linked gA proteins. Planar bilayers with various thicknesses were formed from an appropriate combination of monoglyceride with various fatty acid lengths and solvent. Bilayer thicknesses ranged from 25 A (monoolein in squalene) to 54 A (monoeicosenoin in decane). Single-channel conductances to protons (g(H)) were measured in the concentration range of 10-5000 mM HCl. In native gA as well as in RR channels, the shape of the log(g(H))-log([H(+)]) relationships was nonlinear and remained basically unaltered in monoglyceride bilayers with various thicknesses. For both native gA and RR channels, g(H) values were systematically and significantly larger in thin than in thick bilayers. By contrast, the shape of the log(g(H))-log([H(+)]) relationships in the SS channel was linear (with a slope considerably smaller than 1) in thick (>37 A) bilayers. However, in thin (<37 A) bilayers these plots became nonlinear and g(H) values approached those obtained in native gA channels. The linearization of the log-log plots in the SS channel in thick bilayers is a consequence of a dramatic increase (instead of a decrease as in native gA and RR channels) of g(H) in these bilayers in [H(+)] <1 M. The gating characteristics of the various gA channels as a function of bilayer thickness followed the same pattern as described previously. It was noticed, however, that in the thickest monoglyceride bilayer used in this study, both the SS- and RR-dioxolane-linked channels opened in a mode of bursting activity instead of remaining in the open state as in thin bilayers. It is proposed that the thickness of monoglyceride bilayers modulates proton transfer in native gA channels by a combination of factors including the access resistances of channels to H(+), and fluctuations in both the structure of the lipid bilayer and in the distance between gA monomers. The differential effects of relatively thick monoglyceride bilayers on proton transfer in both dioxolane-linked gA channels must relate to distinct interactions between the bilayers and the SS and RR dioxolanes.     

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.     

Gating and Permeation in Ion Channels Formed by Gramicidin A and Its Dioxolane-linked Dimer in Na+ and Cs+ Solutions

The association of two gramicidin A (gA) peptides via H-bonds in lipid bilayers causes the formation of an ion channel that is selective for monovalent cations only. In this study, two gAs were covalently linked with a dioxolane group (SS dimer). Some functional properties of natural gA channels were compared to that synthetic dimer in Na⁺- or Cs⁺-containing solutions. The SS dimer remained in the open configuration most of the time, while natural gA channels had a relatively brief mean open time. Single channel conductances to Na⁺ (g _Na ) or Cs⁺ (g _Cs ) in the SS dimer were smaller than in natural gA. However, g _Na was considerably more attenuated than g _Cs . This probably results from a tight solvation of Na⁺ by the dioxolane linker in the SS channel. In Cs⁺ solutions, the SS had frequent closures. By contrast, in Na⁺ solutions the synthetic dimer remained essentially in the open state. The mean open times of SS channels in different solutions (T _open,Na > T _open,Cs > T _open,H ) were inversely proportional to the single channel conductances (g _H > g _Cs > g _Na ). This suggests that ion occupancy inside the pore stabilizes the open configuration of the gA dimer. The mean closed time of the SS dimer was longer in Cs⁺ than in H⁺ solutions. Possible mechanisms for these effects are discussed. Received: 17 September 1999/Revised: 21 December 1999     

Gramicidin channels in phospholipid bilayers with unsaturated acyl chains.

In organic solvents gramicidin A (gA) occurs as a mixture of slowly interconverting double-stranded dimers. Membrane-spanning gA channels, in contrast, are almost exclusively single-stranded beta(6,3)-helical dimers. Based on spectroscopic evidence, it has previously been concluded that the conformational preference of gA in phospholipid bilayers varies as a function of the degree of unsaturation of the acyl chains. Double-stranded pi pi(5,6)-helical dimers predominate (over single-stranded beta(6,3)-helical dimers) in lipid bilayer membranes with polyunsaturated acyl chains. We therefore examined the characteristics of channels formed by gA in 1-palmitoyl-2-oleoylphosphatidylcholine/n-decane, 1,2-dioleoylphosphatidylcholine/n-decane, and 1,2-dilinoleoylphosphatidylcholine/n-decane bilayers. We did not observe long-lived channels that could be conducting double-stranded pi pi(5,6)-helical dimers in any of these different membrane environments. We conclude that the single-stranded beta(6,3)-helical dimer is the only conducting species in these bilayers. Somewhat surprisingly, the average channel duration and channel-forming potency of gA are increased in dilinoleoylphosphatidylcholine/n-decane bilayers compared to 1-palmitoyl-2-oleoylphosphatidylcholine/n-decane and dioleoylphosphatidylcholine/n-decane bilayers. To test for specific interactions between the aromatic side chains of gA and the acyl chains of the bilayer, we examined the properties of channels formed by gramicidin analogues in which the four tryptophan residues were replaced with naphthylalanine (gN), tyrosine (gT), and phenylalanine (gM). The results show that all of these analogue channels experience the same relative stabilization when going from dioleoylphosphatidylcholine to dilinoleoylphosphatidylcholine bilayers.     

Modulation of proton transfer in the water wire of dioxolane-linked gramicidin channels by lipid membranes

Proton conductance (g(H)) in single SS stereoisomers of dioxolane-linked gramicidin A (gA) channels were measured in different phospholipid bilayers at different HCl concentrations. In particular, measurements were obtained in bilayers made of 1,2-diphytanoyl 3-phosphocholine (DiPhPC) or its ethylated derivative 1,2-diphytanoyl 3-ethyl-phosphocholine (et-DiPhPC,). The difference between these phospholipids is that in et-DiPhPC one of the phosphate oxygens is covalently linked to an ethyl group and cannot be protonated. In relatively dilute acid solutions, g(H) in DiPhPC is significantly higher than in et-DiPhPC. At high acid concentrations, g(H) is the same in both diphytanoyl bilayers. Such differences in g(H) can be accounted for by surface charge effects at the membrane/solution interfaces. In the linear portion of the log g(H)-log [H] relationship, g(H) values in diphytanoyl bilayers were significantly larger (approximately 10-fold) than in neutral glyceryl monooleate (GMO) membranes. The slopes of the linear log-log relationships between g(H) and [H] in diphytanoyl and GMO bilayers are essentially the same (approximately 0.76). This slope is significantly lower than the slope of the log-log plot of proton conductivity versus proton concentration in aqueous solutions (approximately 1.00). Because the chemical composition of the membrane-channel/solution interface is strikingly different in GMO and diphytanoyl bilayers, the reduced slope in g(H)-[HCl] relationships may be a characteristic of proton transfer in the water wire inside the SS channel. Values of g(H) in diphytanoyl bilayers were also significantly larger than in membranes made of the more common biological phospholipids 1-palmitoyl 2-oleoyl phosphocholine (POPC) or 1-palmitoyl 2-oleoyl phosphoethanolamine (POPE). These differences, however, cannot be accounted for by different surface charge effects or by different internal dipole potentials. On the other hand, maximum g(H) measured in the SS channel does not depend on the composition of the bilayer and is determined essentially by the reduced mobility of protons in concentrated acid solutions. Finally, no experimental evidence was found in support of a lateral proton movement at the phospholipid/solution interface contributing to g(H) in single SS channels. Protein-lipid interactions are likely to modulate g(H) in the SS channel.     

Kinetic isotope effects of proton transfer in aqueous and methanol containing solutions,and in gramicidin A channels

The electrochemical conductivities of HCL and DCI were measured in: H(2)O and D(2)O; in methanol and fully deuterated methanol; and in water-methanol solutions. The single channel conductances to H(+) (g(H)) and D(+) (g(D)) in various gramicidin A (gA) ion channels incorporated in glycerylmonooleate planar bilayers were also measured. Kinetic isotope effects (KIE) were estimated from the ratio of conductivity measurements. In 1 and 5 M HCl aqueous solutions and in 1 M HCl+3.7 M methanol, the KIE ( approximately 1.35) is not different from values previously determined in dilute acid solutions. This suggests that the mobility of protons in those solutions is largely determined by proton transfer. In 10 M HCl, however, where the mobility of protons is likely to be determined by hydrodynamic diffusion, the measured KIE is considerably larger (1.47). Possible causes for this effect are discussed. The KIE of proton conductivities in 5 and 50 mM HCl in methanol and d-methanol is approximately 1.15. This is considerably smaller than the ratio between conductivities of 5 mM KCl in methanol and d-methanol (1.24). The KIE values (1.22-1.37) for g(H) in gA channels in 1 M HCl are significantly larger than for other monovalent cations and consistent with H(+) transfer. Methanol reduces g(H) in gA channels. The KIE of this effect is not different from the one measured in the absence of methanol. Possible mechanisms for the methanol-induced block of H(+) conductivities in solution and gA channels are discussed.     

Proton transfer in gramicidin water wires in phospholipid bilayers: attenuation by phosphoethanolamine

The transfer of protons in water wires was studied in native gramicidin A (gA), and in the SS- and RR-diastereoisomers of dioxolane-linked gA channels (SS and RR channels). These peptides were incorporated into membranes comprised of distinct combinations of phospholipid headgroups and acyl chains. Quantitative relationships between single channel conductances to H+ (g(H)) and [H+] were determined in distinct phospholipid membranes, and are in remarkable contrast with results previously obtained in monoglyceride membranes. In particular: 1), g(H)-[H+] relationships for the various gA channels in distinct phospholipid membranes are well fitted by single adsorption isotherms. A simple kinetic model assuming mono-occupancy of channels by protons fits said relationships. This does not occur with monoglyceride membranes. 2), Under nonsaturating [H+], g(H) is approximately 1 order of magnitude larger in phospholipid than in monoglyceride membranes. 3), Differences between rates of H+ transfer in various gA channels are still present but considerably attenuated in phospholipid relative to monoglyceride membranes. 4), Charged phospholipid headgroups affect g(H) via changes in [H+] at the membrane/solution interfaces. 5), Phosphoethanolamine groups caused a marked attenuation of g(H) relative to membranes with other phospholipid headgroups. This attenuation is voltage-dependent and tends to saturate H+ currents at voltages larger than 250 mV. This effect is likely to occur by limiting the access and exit of H+ in and out of the channel due to relatively strong oriented H-bonds between waters and phosphoethanolamine groups at channel interfaces. The differential effects of phospholipids on proton transfer could be reasoned by considering solvation effects of side chain residues of gramicidin channels by double acyl chains and by the presence of polar headgroups facilitating the entrance/exit of protons through the channel mouths.     

Attenuation of proton currents by methanol in a dioxolane-linked gramicidin A channel in different lipid bilayers.

The mobility of protons in a dioxolane-linked gramicidin A channel (D1) is comparable to the mobility of protons in aqueous solutions (Cukierman, S., E. P. Quigley, and D. S. Crumrine. 1997. Biophys. J. 73:2489-2502). Aliphatic alcohols decrease the mobility of H+ in aqueous solutions. In this study, the effects of methanol on proton conduction through D1 channels were investigated in different lipid bilayers and at different HCl concentrations. Methanol attenuated H+ currents in a voltage-independent manner. Attenuation of proton currents was also independent of H+ concentrations in solution. In phospholipid bilayers, methanol decreased the single channel conductance to protons without affecting the binding affinity of protons to bilayers. In glycerylmonooleate membranes, the attenuation of single channel proton conductances qualitatively resembled the decrease of conductivities of HCl solutions by methanol. However, in both types of lipid bilayers, single channel proton conductances through D1 channels were considerably more attenuated than the conductivities of different HCl solutions. This suggests that methanol modulates single proton currents through D1 channels. It is proposed that, on average, one methanol molecule binds to a D1 channel, and attenuates H+ conductance. The Gibbs free energy of this process (DeltaG0) is approximately 1.2 kcal/mol, which is comparable to the free energy of decrease of HCl conductivity in methanol solutions (1.6 kcal/mol). Apolar substances like urea and glucose that do not transport protons in HCl solutions and do not permeate D1 channels decreased solution conductivity and single channel conductance by a considerably larger proportion than methanol. Cs+ currents through D1 channels were considerably less (fivefold) attenuated by methanol than proton currents. It is proposed that methanol partitions inside the pore of gramicidin channels and delays the transfer of protons between water and methanol molecules, causing a significant attenuation of the single channel proton conductance. Gramicidin channels offer an interesting experimental model to study proton hopping along a single chain of water molecules interrupted by a single methanol molecule.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

Effect of Linker Length on Avidin Binding to Biotinylated Gramicidin A

Biotinylated gramicidins are an important component of the AMBRI^® “ion channel switch™” biosensor. These gramicidin A (gA) analogues have a biotin attached to the C-terminus of gA via a number of aminocaproyl linker groups (X). The structure of gA5XB has been determined in deuterated sodium dodecyl sulfate micelles and is similar to native gA and other modified gA analogues. The biotin and aminocaproyl groups were mobile and located in the aqueous phase and when avidin was added, NMR and MS studies showed that gA5XB bound more effectively to avidin than gA2XB. The length and flexibility of the linker appears to be important for biotin–avidin binding and, in the AMBRI^® biosensor, gA5XB is a more effective gated ion channel than gA2XB. The conformation and dynamics of the aminocaproyl linker groups were investigated using ²H solid-state NMR. Deuterated aminocaproyl linkers were coupled to gA and incorporated into oriented bilayers in order to analyse the order and dynamics of the aminocaproyl linker. The small ²H splittings and the T ₁ relaxation times indicated that the aminocaproyl linker is undergoing fast rotation in phospholipid bilayers. Native d ₄ -gA as well as d ₄ -gA2XB, where the ethanolamine has been deuterated, were also incorporated into oriented bilayers. Solid-state ²H NMR data showed that the addition of the linker group restricted the mobility of the ethanolamine. However, these modifications to the C-terminus of gA did not interfere with ion channel function and clarify how the biotinylated gA analogues perform in the lipid bilayer as part of the AMBRI^® biosensor.Australian Peptide Conference Issue.     

Genistein can modulate channel function by a phosphorylation-independent mechanism: importance of hydrophobic mismatch and bilayer mechanics

Genistein, a generic tyrosine kinase inhibitor, has been used extensively as a tool to investigate the possible regulation of membrane function by tyrosine phosphorylation. Genistein, in micromolar concentrations, alters the function of numerous ion channels and other membrane proteins, but only in few cases has it been demonstrated that the changes in membrane protein (ion channel) function are due to changes in a protein's phosphorylation status. The major common denominator characterizing proteins that are modulated by genistein seems to be that they are imbedded into, and span, the bilayer component of the plasma membrane. We therefore explored whether genistein could alter ion channel function by a bilayer-mediated mechanism and examined genistein's effect on gramicidin A (gA) channels in planar phospholipid bilayers. gA channels form by transmembrane dimerization of two nonconducting subunits, and genistein potentiates gA channel activity by increasing the appearance rate and prolonging the lifetime of bilayer-spanning gA dimers. That is, genistein shifts the equilibrium between nonconducting monomers and conducting dimers in favor of the bilayer-spanning dimers; the changes in channel activity therefore cannot be due to changes in bilayer fluidity. To obtain further insights into the mechanism underlying this modulation of gA channel function, we examined the effects of genistein on channels formed by gA analogues that differ in amino acid sequence. For a given channel length, the effects of genistein on gA dimerization do not depend on the specific sequence, or the chirality, of the channel-forming gA analogues. In contrast, when we change the channel length (by decreasing or increasing the number of amino acid residues in the sequence), or the bilayer thickness (by changing methylene groups in the acyl chains), the magnitude of genistein's effect increases with increasing hydrophobic mismatch between the channel length and the bilayer thickness. These results strongly suggest that genistein alters bilayer mechanical properties, which in turn modulates channel function. This bilayer-mediated mechanism is likely to apply to other pharmacological reagents and membrane proteins.     

Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A.

Proton transfer in biological systems is thought to often proceed through hydrogen-bonded chains of water molecules. The ion channel, gramicidin A (gA), houses within its helical structure just such a chain. Using the density functional theory based ab initio molecular dynamics Car-Parrinello method, the structure and dynamics of proton diffusion through a polyglycine analog of the gA ion channel has been investigated. In the channel, a proton, which is initially present as hydronium (H3O+), rapidly forms a strong hydrogen bond with a nearest neighbor water, yielding a transient H5O2+ complex. As in bulk water, strong hydrogen bonding of this complex to a second neighbor solvation shell is required for proton transfer to occur. Within gA, this second neighbor shell included not only a channel water molecule but also a carbonyl of the channel backbone. The present calculations suggest a transport mechanism in which a priori carbonyl solvation is a requirement for proton transfer.     

The membrane interface dictates different anchor roles for "inner pair" and "outer pair" tryptophan indole rings in gramicidin A channels

We investigated the effects of substituting two of the four tryptophans (the "inner pair" Trp(9) and Trp(11) or the "outer pair" Trp(13) and Trp(15)) in gramicidin A (gA) channels. The conformational preferences of the doubly substituted gA analogues were assessed using circular dichroism spectroscopy and size-exclusion chromatography, which show that the inner tryptophans 9 and 11 are critical for the gA's conformational preference in lipid bilayer membranes. [Phe(13,15)]gA largely retains the single-stranded helical channel structure, whereas [Phe(9,11)]gA exists primarily as double-stranded conformers. Within this context, the (2)H NMR spectra from labeled tryptophans were used to examine the changes in average indole ring orientations, induced by the Phe substitutions and by the shift in conformational preference. Using a method for deuterium labeling of already synthesized gAs, we introduced deuterium selectively onto positions C2 and C5 of the remaining tryptophan indole rings in the substituted gA analogues for solid-state (2)H NMR spectroscopy. The (least possible) changes in orientation and overall motion of each indole ring were estimated from the experimental spectra. Regardless of the mixture of backbone folds, the indole ring orientations observed in the analogues are similar to those found previously for gA channels. Both Phe-substituted analogues form single-stranded channels, as judged from the formation of heterodimeric channels with the native gA. [Phe(13,15)]gA channels have Na(+) currents that are ~50% and lifetimes that are ~80% of those of native gA channels. The double-stranded conformer(s) of [Phe(9,11)]gA do not form detectable channels. The minor single-stranded population of [Phe(9,11)]gA forms channels with Na(+) currents that are ~25% and single-channel lifetimes that are ~300% of those of native gA channels. Our results suggest that Trp(9) and Trp(11), when "reaching" for the interface, tend to drive both monomer folding (to "open" a channel) and dimer dissociation (to "close" a channel). Furthermore, the dipoles of Trp(9) and Trp(11) are relatively more important for the single-channel conductance than are the dipoles of Trp(13) and Trp(15).     

Proton transfer in water wires in proteins: modulation by local constraint and polarity in gramicidin a channels

The transfer of protons in membrane proteins is an essential phenomenon in biology. However, the basic rules by which H(+) transfer occurs in water wires inside proteins are not well characterized. In particular, the effects of specific atoms and small groups of atoms on the rate of H(+) transfer in water wires are not known. In this study, new covalently linked gramicidin-A (gA) peptides were synthesized, and the effects of specific atoms and peptide constraints on the rate of H(+) transfer were measured in single molecules. The N-termini of two gA peptides were linked to various molecules: S,S-cyclopentane diacid, R,R-cyclopentane diacid, and succinic acid. Single-channel proton conductances (g(H)) were measured at various proton concentrations ([H(+)]) and compared to previous measurements obtained in the S,S- and R,R-dioxolane-linked as well as in native gA channels. Replacing the S,S-dioxolane by an S,S-cyclopentane had no effects on the g(H)-[H(+)] relationships, suggesting that the constrained and continuous transition between the two gA peptides via these S,S linkers is ultimately responsible for the two- to fourfold increase in g(H) relative to native gA channels. It is likely that constraining a continuous transition between the two gA peptides enhances the rate of H(+) transfer in water wires by decreasing the number of water wire configurations that do not transfer H(+) at higher rates as in native gA channels (a decrease in the activation entropy of the system). On the other hand, g(H) values in the R,R-cyclopentane are considerably larger than those in R,R-dioxolane-linked gA channels. One explanation would be that the electrostatic interactions between the oxygens in the dioxolane and adjacent carbonyls in the R,R-dioxolane-linked gA channel attenuate the rate of H(+) transfer in the middle of the pore. Interestingly, g(H)-[H(+)] relationships in the R,R-cyclopentane-linked gA channel are quite similar to those in native gA channels. g(H) values in succinyl-linked gA channels display a wide distribution of values that is well represented by a bigaussian. The larger peaks of these distributions are similar to g(H) values measured in native gA channel. This observation is also consistent with the notion that constraining the transition between the two beta-helical gA peptides enhances the rate of H(+) transfer in water wires by decreasing the activation entropy of the system.     

Membrane-Mediated Protein-Protein Interactions and Connection to Elastic Models: A Coarse-Grained Simulation Analysis of Gramicidin A Association

To further foster the connection between particle based and continuum mechanics models for membrane mediated biological processes, we carried out coarse-grained (CG) simulations of gramicidin A (gA) dimer association and analyzed the results based on the combination of potential of mean force (PMF) and stress field calculations. Similar to previous studies, we observe that the association of gA dimers depends critically on the degree of hydrophobic mismatch, with the estimated binding free energy of >10 kcal/mol in a distearoylphosphatidylcholine bilayer. Qualitative trends in the computed PMF can be understood based on the stress field distributions near a single gA dimer and between a pair of gA dimers. For example, the small PMF barrier, which is ∼1 kcal/mol independent of lipid type, can be captured nearly quantitatively by considering membrane deformation energy associated with the region confined by two gA dimers. However, the PMF well depth is reproduced poorly by a simple continuum model that only considers membrane deformation energy beyond the annular lipids. Analysis of lipid orientation, configuration entropy, and stress distribution suggests that the annular lipids make a significant contribution to the association of two gA dimers. These results highlight the importance of explicitly considering contributions from annular lipids when constructing approximate models to study processes that involve a significant reorganization of lipids near proteins, such as protein-protein association and protein insertion into biomembranes. Finally, large-scale CG simulations indicate that multiple gA dimers also form clusters, although the preferred topology depends on the protein concentration. Even at high protein concentrations, every gA dimer requires contact to lipid hydrocarbons to some degree, and at most three to four proteins are in contact with each gA dimer; this observation highlights another aspect of the importance of interactions between proteins and annular lipids.