Calculation of deformation energies and conformations in lipid membranes containing gramicidin channels

Attenuation of channel kinetics and conductance by cholesterol: an interpretation using structural stress as a unifying concept

The ubiquity of cholesterol in cell membranes and changes in its concentration during development, aging and in various diseases suggest that it plays an important role in modulating cell function. We examined this possibility by monitoring the effects of cholesterol on the activity of the calcium-activated potassium (BK) channel reconstituted into lipid bilayers from rat brain homogenates. Increasing the cholesterol concentration to 11% of total lipid weight resulted in a 70% reduction in channel mean open time and a reduction of the open probability of the channel by 80%. Channel conductance was reduced by 7%. Cholesterol is known to change the order state and the modulus of compressibility of bilayers. These physico-chemical changes may be translated into an overall increase in the structural stress in the bilayer, and this force may be transmitted to proteins residing therein. By examining the characteristics of the BK channel as a function of temperature, in the presence and absence of cholesterol, we were able to estimate the activation energy based on Arrhenius plots of channel kinetics. Cholesterol reduced the activation energy of the BK channel by 50% for the open to closed transition. This result is consistent with an increased stress energy in the bilayer and favors the channel moving into the closed state. Taken together, these data are consistent with a model in which cholesterol induces structural stress which enhances the transition from the open to the closed state of the channel. We suggest that this is an important mechanism for regulating the activity of membrane-integral proteins and therefore membrane function, and that the concept of structural stress may be relevant to understanding the modulation of ion channel activity in cell membranes.

Lysophospholipids modulate channel function by altering the mechanical properties of lipid bilayers

Lipid metabolites, free fatty acids and lysophospholipids, modify the function of membrane proteins including ion channels. Such alterations can occur through signal transduction pathways, but may also result from "direct" effects of the metabolite on the protein. To investigate possible mechanisms for such direct effects, we examined the alterations of gramicidin channel function by lysophospholipids (LPLs): lysophosphatidylcholine (LPC), lysophosphatidylethanolamine (LPE), lysophosphatidylserine (LPS), and lysophosphatidylinositol (LPI). The experiments were done on planar bilayers formed by diphytanoylphosphatidylcholine in n-decane a system where receptor-mediated effects can be excluded. At aqueous concentrations below the critical micelle concentration (CMC), LPLs can increase the dimerization constant for membrane-bound gramicidin up to 500-fold (at 2 microM). The relative potency increases as a function of the size of the polar head group, but does not seem to vary as a function of head group charge. The increased dimerization constant results primarily from an increase in the rate constant for channel formation, which can increase more than 100-fold (in the presence of LPC and LPI), whereas the channel dissociation rate constant decreases only about fivefold. The LPL effect cannot be ascribed to an increased membrane fluidity, which would give rise to an increased channel dissociation rate constant. The ability of LPC to decrease the channel dissociation rate constant varies as a function of channel length (which is always less than the membrane's equilibrium thickness): as the channel length is decreased, the potency of LPC is increased. LPC has no effect on membrane thickness or the surface tension of monolayers at the air/electrolyte interface. The bilayer-forming glycerolmonooleate does not decrease the channel dissociation rate constant. These results show that LPLs alter gramicidin channel function by altering the membrane deformation energy, and that the changes in deformation energy can be related to the molecular "shape" of the membrane-modifying compounds. Similar alterations in the mechanical properties of biological membranes may form a general mechanism by which one can alter membrane protein function.

Gramicidin A/short-chain phospholipid dispersions: chain length dependence of gramicidin conformation and lipid organization

Gramicidin-lipid interactions were investigated using diacylphosphatidylcholines that contained two identical acyl chains of varying length, between 6 and 14 carbons. The gramicidin A (gA) conformation was monitored by circular dichroism (CD) spectroscopy and high-performance size-exclusion chromatography, and the lipid organization was investigated using 31P and 1H NMR spectroscopy and negative-stain electron microscopy. Diacylphosphatidylcholine (PC) lipids with chain lengths between 4 and 8 carbons have been previously shown to have a micellar organization in aqueous solution [Lin, T.-L., et al. (1986) J. Am. Chem. Soc. 108, 3499-3507]. CD spectra of aqueous gA/lipid dispersions, at a ratio of 1:28, demonstrated that the channel conformation of gA can be readily obtained when the acyl chain length is > or = 10, but not when the chain length is < or = 7. Size-exclusion chromatography revealed that the fraction of gA that could easily be dissociated into monomers in the dispersions increased with increasing acyl chain length, in agreement with the CD results. For a chain length of 8, the results were intermediate. The formation of the channel structure was found to depend on the "solvent-history", the temperature, the gA and lipid concentrations, the gA:lipid ratio, and consequently on the method of sample preparation. 1H and 31P NMR results suggest that codispersed gA increases the size of dioctanoyl-PC aggregates, but not of dihexanoyl-PC micelles. Negative-stain electron microscopy directly supports these findings. Dihexanoyl-PC (28 mM) was able to solubilize 1 mM gA in H2O, but the gA was not in the "channel" conformation.(ABSTRACT TRUNCATED AT 250 WORDS)

Some effects of short-chain phospholipids and n-alkanes on a transient potassium current (IA) in identified Helix neurons

Many effects of short-chain phospholipids and n-alkanes on the squid axon sodium current (INa) are consistent with mechanisms involving changes in membrane thickness. Here, we suggest that the actions of short-chain phospholipids on an A-type potassium current (IA) in two-microelectrode voltage clamped Helix D1 and F77 neurons are incompatible with such simple mechanisms. Diheptanoyl phosphatidylcholine (diC7PC, 0.2 and 0.3 mM) caused substantial (58 and 79%), and in some cases partially reversible, increases in IA amplitude. These were correlated with hyperpolarizing shifts of up to -7 mV in the voltage dependence of current activation. The voltage dependence of steady-state inactivation was also moved in the hyperpolarizing direction. These effects are the opposite of those described for squid INa. 0.5 Saturated n-pentane and saturated n-hexane caused significant (-3 and -6 mV) hyperpolarizing shifts in the voltage dependence of IA inactivation, qualitatively consistent with their effects on squid INa, while the voltage dependence of activation was moved slightly to the left or unchanged. Hydrocarbons had variable effects on peak current amplitude, although saturated n-pentane produced a clear suppression. DiC7PC caused a 25% increase in the time constant of macroscopic IA inactivation (tau b) but 0.5 saturated n-pentane and saturated n-hexane reduced tau b by 40%. The effects of these agents on current-clamped cells were broadly consistent with their opposing actions on tau b--phospholipids tended to reduce excitability and n-alkanes tended to increase it. Possible mechanisms of IA perturbation are discussed.

A molecular model for lipid-protein interaction in membranes: the role of hydrophobic mismatch

The interaction free energy between a hydrophobic, transmembrane, protein and the surrounding lipid environment is calculated based on a microscopic model for lipid organization. The protein is treated as a rigid hydrophobic solute of thickness dP, embedded in a lipid bilayer of unperturbed thickness doL. The lipid chains in the immediate vicinity of the protein are assumed to adjust their length to that of the protein (e.g., they are stretched when dP > doL) in order to bridge over the lipid-protein hydrophobic mismatch (dP-doL). The bilayer's hydrophobic thickness is assumed to decay exponentially to its asymptotic, unperturbed, value. The lipid deformation free energy is represented as a sum of chain (hydrophobic core) and interfacial (head-group region) contributions. The chain contribution is calculated using a detailed molecular theory of chain packing statistics, which allows the calculation of conformational properties and thermodynamic functions (in a mean-field approximation) of the lipid tails. The tails are treated as single chain amphiphiles, modeled using the rotational isometric state scheme. The interfacial free energy is represented by a phenomenological expression, accounting for the opposing effects of head-group repulsions and hydrocarbon-water surface tension. The lipid deformation free energy delta F is calculated as a function of dP-doL. Most calculations are for C14 amphiphiles which, in the absence of a protein, pack at an average area per head-group ao approximately equal to 32 A2 (doL approximately 24.5 A), corresponding to the fluid state of the membrane. When dP = doL, delta F > 0 and is due entirely to the loss of conformational entropy experienced by the chains around the protein. When dP > doL, the interaction free energy is further increased due to the enhanced stretching of the tails. When dP < doL, chain flexibility (entropy) increases, but this contribution to delta F is overcounted by the increase in the interfacial free energy. Thus, delta F obtains a minimum at dP-doL approximately 0. These qualitative interpretations are supported by detailed numerical calculations of the various contributions to the interaction free energy, and of chain conformational properties. The range of the perturbation of lipid order extends typically over few molecular diameters. A rather detailed comparison of our approach to other models is provided in the discussion.

The permeation properties of small organic cations in gramicidin A channels

The conductance properties of organic cations in single gramicidin A channels were studied using planar lipid bilayers. From measurements at 10 mM and at 27 mV the overall selectivity sequence was found to be NH4+ > K+ > hydrazinium > formamidinium > Na+ > methylammonium, which corresponds to Eisenman polyatomic cation sequence X'. Methylammonium and formamidinium exhibit self block, suggesting multiple occupancy and single filing. Formamidinium has an apparent dissociation constant (which is similar to those of alkali metal cations) for the first ion being 22 mM from the Eadie-Hofstee plot (G0 vs. G0/C), 12 mM from the rate constants of a three-step kinetic model. The rate-limiting step for formamidinium is translocation judging from supralinear I-V relations at low concentrations. 1 M formamidinium solutions yields exceptionally long single channel lifetimes, 20-fold longer than methylammonium, which yields lifetimes similar to those found with alkali metal cations. The average lifetime in formamidinium solution significantly decreases with increasing voltage up to 100 mV but is relatively voltage independent between 100 and 200 mV. At lower voltages (< or = 100 mV), the temperature and concentration dependences of the average lifetime of formamidinium were steep. At very low salt concentrations (0.01 M, 100 mV), there was no significant difference in average lifetime from that formed with 0.01 M methylammonium or hydrazinium. We conclude that formamidinium very effectively stabilizes the dimeric channel while inside the channel and speculate that it does so by affecting tryptophan-reorientation or tryptophan-lipid interactions at binding sites.

Model ion channels: gramicidin and alamethicin

We have discussed in some detail a variety of experimental studies which were designed to elucidate the conformational and dynamic properties of gramicidin and alamethicin. Although the behavior of these peptides is by no means fully characterized, these studies have already permitted aspects of ion channel activity to be understood in molecular terms. Studies with gramicidin in a variety of organic solutions have revealed conformational heterogeneity of this peptide; at least five major isomers exist, several of which have been characterized in detail using NMR spectroscopy and X-ray crystallography. When added to lipid membranes gramicidin undergoes a further conformational conversion. Although the conformation of gramicidin in membranes is not as well characterized as the solution conformation(s) and an X-ray structure is not yet available, detailed data, particularly from solid-state NMR studies, continue to become available and a right-handed beta 6.3 helical conformation of the peptide backbone is now generally accepted. Two of these beta 6.3 helices joined at their N-termini are believed to form the conducting channel. The conformational behavior of the side-chains of gramicidin in the membrane-bound form is not well established and several NMR, CD, fluorescence and theoretical studies are now focussed on this. Although the side-chains do not directly contact the permeating ions, they can have distinct effects on conductance and selectivity by altering the electrostatic environment sensed by the ion. The dynamics of both side-chain and backbone conformations of gramicidin appear critical to a detailed understanding of the ion transport process in this channel. As the description of the membrane-bound conformation of gramicidin becomes more detailed, simulations of ion transport using computational methods are likely to improve and will further our understanding of the processes of ion transport. As well as internal motion of the backbone and side-chains, gramicidin undergoes rotational and translational motion in the plane of the membrane. These motions do not appear to be essential for the process of ion transport but can affect channel lifetime since lifetime is determined by the rate of association and dissociation of gramicidin monomers. Gramicidin-membrane interactions are also likely to be involved in the frequency of occurrence of channel subconductance states, the frequency of channel flickering and fundamentally in the stability of the membrane-bound gramicidin conformation. Alamethicin forms channels in membranes which are strongly voltage-dependent. The molecular origin of voltage-dependent conductances has been a fundamental problem in biophysics for many years.(ABSTRACT TRUNCATED AT 400 WORDS)

Influence of ion occupancy and membrane deformation on gramicidin A channel stability in lipid membranes

The average lifetime of gramicidin A channels in monoolein/decane bilayer membranes was measured. The results support the hypothesis of channel stabilization by ion occupancy. The effects of electric field and salt concentration are consistent with the expected effects on both occupancy and membrane compression. The lifetime in asymmetric solutions with divalent cation blockers on one side of the membrane shows a voltage dependence such that the lifetime decreases for positive voltages applied from the blocking side and increases for negative voltages. This result strongly supports the occupancy hypothesis. The lifetime increases with permeant ion concentration, and at the one molar level it also increases with voltage. The voltage dependence of lifetime for a low concentration of permeant ion depends on the total salt level. The results for these conditions are consistent with the assumption that membrane compression also influences the lifetime, even for the "soft" solvent-containing membrane considered here. It is proposed that the channel nearest neighbor lipids need not be fixed in a plane at the channel end. Using a liquid crystal model it may then be shown that surface tension is the major component of the membrane deformation free energy, which may explain the significant effects of the membrane compression on the lifetime.

On the supramolecular organization of gramicidin channels. The elementary conducting unit is a dimer

The question, whether the conducting channels formed by the linear gramicidins are dimers (as is generally believed) or tetramers (as has been recently proposed [Stark G., M. Strässle, and Z. Takacz. 1986. J. Membr. Biol. 89:23-37; Strässle, M., G. Stark, M. Wilhelm, P. Daumas, F. Heitz, and R. Lazaro. 1989. Biochim. Biophys. Acta. 980:305-314]) has been addressed in single-channel experiments. The experimental approach was based on the ability of electrophysiological (single-channel) experiments to resolve the number of hybrid channel types that could form between gramicidin A or C and O-pyromellityl-gramicidin A or C (in which a pyromellitic acid residue has been esterified to the ethanolamine-OH group [Apell, H.-J., E. Bamberg, H. Alpes, and P. Läuger. 1977. J. Membr. Biol. 31:171-188]). The presence of the bulky, negatively charged pyromellityl group at the channel entrances endows the hybrid channels with characteristically different features and thus facilitates the resolution of the different hybrid channel types. Only two hybrid channel types were detected, indicating that the conducting channels are membrane-spanning dimers. There was likewise no evidence for lateral association between conducting channels and nonconducting monomers. These results can be reconciled with those of Stark et al. (op. cit.) if gramicidin channel formation involves a (slow) folding into beta 6.3-helical monomers followed by the dimerization step.

Monitoring the surface tension of lipid membranes by a bubble method

A new method for the determination of lipid membrane surface tension was developed. Advantages of this method are that it allows for multiple measurements on a single membrane, is fast and direct not requiring empirical corrections, may be applied for dynamical surface-tension measurements and may be used with thin films and asymmetrical electrolytes. The pressure and radius of a bubble are measured. A piezoresistive sensor is used to minimize the transducer compliance. By moulding the sensor to a brass plate a resolution of 0.025 mm H2O (0.25 Pa) is obtained. The bubble is filmed using a videocamera and the radius of the bubble determined with the aid of a microcomputer. Data for monoolein/hexadecane in potassium chloride solutions and a cooling curve are presented and compared with previous results.

Effective pore radius of the gramicidin channel. Electrostatic energies of ions calculated by a three-dielectric model

Electrostatic calculation of the gramicidin channel is performed on the basis of a three-dielectric model in which the peptide backbone of the channel is added as a third dielectric region to the conventional two-dielectric channel model (whose pore radius is often referred to as the effective pore radius reff). A basic principle for calculating electrostatic fields in three-dielectric models is introduced. It is shown that the gramicidin channel has no unique value of reff. The reff with respect to the "self-image energy" (i.e., the image energy in the presence of a single ion) is 2.6-2.7 A, slightly depending upon the position of the ion (the least-square value over the whole length of the pore is 2.6 A). In contrast, the reff with respect to the electric potential due to an ion (and hence the reff with respect to the interaction energy between two ions) is dependent upon the distance s of separation; it ranges from 2.6 to greater than 5 A, increasing with an increase in s. However, for the purpose of rough estimation, the reff with respect to the self-image energy can also be used in calculating the electric potential and the interaction energy, because the error introduced by this approximation is an overestimation of the order of 30% at most. It is also shown that the apparent dielectric constant for the interaction between two charges depends markedly upon the positions of the charges. In the course of this study, the dielectric constant and polarizability of the peptide backbone in the beta-sheet structure is estimated to be 10 and 8.22 A3.