The influence of proteins and peptides on the phase properties of lipids

This paper reviews model membrane studies on the modulation of the macroscopic structure of lipids by lipid-protein interactions, with particular emphasis on the gramicidin molecule. This hydrophobic peptide has three main effects on lipid polymorphism: (1) in lysophosphatidylcholine it triggers a micellar to bilayer transition, (2) in phosphatidylethanolamine it lowers the bilayer to hexagonal HII phase transition temperature and (3) in phosphatidylcholine and other bilayer preferring lipids it is able to induce the formation of an HII phase. From experiments in which the gramicidin molecule was chemically modified it can be concluded that the tryptophan residues play a determining role in the peptide-induced changes in polymorphism. The experimental data lead to the proposal that gramicidin molecules have a tendency to self-associate, possibly mediated by tryptophan-tryptophan interactions and organize into tubular structures such as found in the HII phase.

Permeability of lipid bilayers to water and ionic solutes

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

Gramicidin induced aggregation and size increase of phosphatidylcholine vesicles

To investigate the role of membrane proteins in the fusion process, linear hydrophobic polypeptide gramicidin was used as fusogenic agent in small unilamellar vesicles (SUV) constituted of saturated lecithins. It was found that gramicidin, externally added to a suspension of vesicles, induces a reversible vesicles aggregation. When incorporated into the bilayer, gramicidin induces increase in vesicle size. The vesicle size increase was monitored by column chromatography and transmission electron microscopy. The process of vesicle size increase occurs only when the lipid membrane is in the gel state. A maximum is observed in the kinetics at a temperature of approx. 25 degrees C lower than the phase transition temperature of lipids. Higher rates of vesicle size increase are obtained as the lipid chain length increases. The process is accompanied by a release of internal vesicle content and by membrane lipid mixing.

The gramicidin A channel: energetics and structural characteristics of the progression of a sodium ion in the presence of water

The distribution of water molecules in the Gramicidin A (GA) channel is determined by theoretical computations, and the role of this water on the energetics of the system upon progression of a sodium cation through the channel is investigated. In the absence of the ion, water molecules form a chain along the channel, hydrogen bonded to one another and to the L carbonyl oxygens, while others stay at the entrances of the channel, hydrogen-bonded to the free carbonyl oxygens of the L-Tryptophan residues. According to the definition adopted for the "inside" and the "outside" of the channel, it is found to contain at most 7 or 9 water molecules. When a hydrated sodium cation approaches and enters the channel, the structural properties corresponding to the minimized total energy of the system GA-water-Na+ indicate a reorganization, but not a destruction, of the chain of water molecules. The "energy profile" for the system GA-Na+-(22 waters) is analyzed in terms of its components and in comparison to the corresponding intrinsic profile computed earlier in vacuo. It appears that the presence of water does not unduely modify the pathway or the qualitative features of the energetics of the cation passage, except at the entrance, where the partial and progressive dehydration of the cation plays an important role. The presence and characteristics of the minimum found earlier at 10.5 A from the center are conserved.

Comparative 2H- and 31P-NMR study on the properties of palmitoyllysophosphatidylcholine in bilayers with gramicidin, cholesterol and dipalmitoylphosphatidylcholine

The stoichiometric palmitoyllysophosphatidylcholine (lysoPC)/gramicidin (4:1, mol/mol) lamellar complex (Killian, J.A., De Kruijff, B., Van Echteld, C.J.A., Verkleij, A.J., Leunissen-Bijvelt, J. and De Gier, J. (1983) Biochim. Biophys. Acta 728, 141-144) is a useful model system to investigate the various aspects of lipid protein interactions. To study the effect of gramicidin on local order and motion of 1-palmitoyl-sn-glycero-3-phosphocholine (lysoPC) we employed 31P and 2H nuclear magnetic resonance (NMR) using selectively deuterated lysoPC's and we compared the results to those obtained for lysoPC in bilayers with cholesterol (1:1, mol/mol) and dipalmitoylphosphatidylcholine (DPPC) (1:4, mol/mol). 2H-NMR experiments on acyl chain deuterated lysoPC showed similar quadrupole splittings in the liquid crystalline state for the lysoPC/DPPC and the lysoPC/gramicidin samples. In the lysoPC/cholesterol sample an increase of the quadrupole splitting was found. T1 measurements showed that gramicidin decreases the lysoPC acyl chain motion, especially at the C12 position. In the lysoPC/cholesterol sample an increase of motion was observed as compared to lysoPC in fluid bilayers of DPPC. 31P-NMR and 2-H-NMR measurements of lysoPC, deuterated at the alpha- and beta-position of the choline moiety, indicated an increase in headgroup flexibility in all samples as compared to the parent compound DPPC. In addition, a change in headgroup conformation was observed. The alpha- and beta-segments in all samples exhibited concerted motion. It was found that also in the polar headgroup gramicidin induces a decrease of the rate of motion.

Infrared spectroscopic studies on gramicidin ion-channels: relation to the mechanisms of anesthesia

Fourier transform infrared spectroscopic studies are reported on gramicidin ion-channels in phospholipid bilayers and the effects on the spectra of the anesthetics and related compounds (methoxyflurane, halothane, chloroform, carbon tetrachloride, n-pentane and n-decane) have been determined. The addition of anesthetics containing the 'acidic hydrogen' caused unique changes particularly on the amide I bands at 1639 cm-1 and 1670 cm-1. The 1639 cm-1 band became more intense while the intensity near 1670 cm-1 decreased dramatically. These effects were not observed with carbon tetrachloride, n-pentane and n-decane. The 1670 cm-1 band is interpreted as arising from the carbonyls involved in the head-to-head hydrogen-bonded dimerization where the relationship between chains is analogous to that of the antiparallel beta-pleated sheet structure and the anesthetics with 'acidic hydrogens' are considered to disrupt the hydrogen-bonded dimerization by competitive hydrogen bonding to the carbonyls at the head-to-head junction. As the dimer-monomer equilibrium is the 'on-off' mechanism for gramicidin ion-channel conductance, the results are considered in terms of the mechanism of action of anesthetics and are taken to suggest, for certain anesthetics, a hydrogen-bonding role to protein ion-channel components.

Nuclear magnetic resonance of 23Na ions interacting with the gramicidin channel

Basic nuclear magnetic resonance (NMR) features of 23Na ions bound to the gramicidin channel (packaged into lecithin liposomes) were studied. The first binding constant K1 of Na+ was not significantly dependent on channel models employed. With the two-identical-site model (Model I), K1 was 13.7 (+/- 1.4) molal-1 (in the activity basis) at 25 degrees C; when the binding of a third ion was included (Model II), it was 13.0 (+/- 2.0) molal-1. The second binding constant K2 was model dependent; it was 1.6 (+/- 0.2) and 3-4 molal-1 for Models I and II, respectively. The rate constants, k-1 and k-2, of Na+ for exit from singly and doubly loaded channels, respectively, were 8 X 10(5) s-1 less than or equal to k-1 less than or equal to 3 X 10(6) s-1 and 8 X 10(5) s-1 less than or equal to k-2 less than or equal to 1.0 X 10(7) s-1 at 25 degrees C; the lower bound represents a rough approximation of k-1. The ratio k-2/k-1 was greater than one and did not greatly exceed 20. From the competition experiment, K1 of T1+ was 5.7 (+/- 0.6) X 10(2) molal-1. The longitudinal relaxation time T1 of bound 23Na in the state of single occupancy (T 1B sing) was virtually independent of models, 0.56 (+/- 0.03) and 0.55 (+/- 0.04) ms at 25 degrees C for Models I and II, respectively. For the state of double occupancy, T1 of bound 23Na (T 1B doub) was model dependent: 0.27 (+/- 0.01) and 0.4-0.6 ms for Models I and II. The correlation time tau c of bound 23Na was 2.2 (+/- 0.2) ns at 25 degrees C for single occupancy; tau c for double occupancy was not significantly different from this value. The estimated tau c was found to involve no appreciable contribution of the exchange of 23Na between the channel and the bulk solution. Thé quadrupole coupling constant chi was 1.0 (+/- 0.1) MHz for 23Na in single occupancy; chi for double occupancy was 0.9-1.4 MHz, depending on models. A lower bound of the average quadrupole coupling constant chi alpha was 0.13-0.26 MHz at 25 degrees C for 23Na in single occupancy; this value represents a rough approximation of chi alpha at this temperature. An argument based on the estimated chi alpha and the known conformation of the gramicidin channel suggests that the binding site is a small domain near the channel end. Within the framework of Model I, Tb was faster than Tljn; this inequality was attributed to an increased chi in the presence ofa second cation, which was not explained in terms of electrostatic interactions between bound cations, implying a conformation change upon binding of cations.

Poly(DL-proline), a synthetic polypeptide behaving as an ion channel across bilayer membranes

The synthesis and characterization of poly(DL-proline) are reported in relation with its predicted property of forming ion channels across membranes. The analysis of the conductance induced in synthetic bilayer membranes doped with poly(DL-proline) shows ionic permeoselectivity and the characteristic time course of fluctuations of ion channels, according to the similarity with the active structure of gramicidin A in membranes during the ion passage. An alternative mechanism of ion transport across bilayer membranes is also advanced.

Interaction of K+ ion with the solvated gramicidin A transmembrane channel

Using Urry's gramicidin A (GA) atomic coordinates and ab into calculations, the interaction energies of a K+ ion with GA are examined. From these energies the values of the fitting parameters are obtained for 6-12-1 atom-atom pair potentials. The potential of the GA channel as experienced by the ion is analyzed in detail. An energy profile of the K+ ion in the GA channel is obtained by analyzing iso-energy maps. Using Monte Carlo simulations, the energy profiles of the K+ ion with the solvated GA channel are analyzed and the hydration structures in the presence of the K+ ion are studied.

The effect of the amino-acid side chains on the energy profiles for ion transport in the gramicidin A channel

Computations on the energy profiles for Na+ in the gramicidin A (GA) channel have been extended by introducing the effect, previously neglected, of the amino acid side chains of GA, fixed in their most stable conformations. The calculations have been performed in two approximations: 1) with the ethanolamine tail fixed in its most stable conformation, 2) with the tail allowed to optimize its conformation upon the progression of the ion. In both approximations the overall shape of the energy profile is very similar to that obtained in the absence of the side chains. One observes, however, a general lowering of the profile upon the adjunction of the side chains. The analysis of the factors responsible for this energy lowering indicates that it is due essentially to the electrostatic and polarisation components of the interaction which interplay differently, however, in the different parts of the channel. A particular role is attributed in this respect to the tryptophan residues of GA. The role of the 4 tryptophans present, Trp 15, 13, 11 and 9, is individualized by stripping of one of them at a time. The strongest effect on the energy deepening is due to Trp 13 and is particularly prominent in the entrance zone at 14.5A from the center of the channel. The result indicates the possibility of investigating theoretically the effect on the energy profiles of the substitution of the "natural" side chain by others.

Molecular dynamics simulation of cation motion in water-filled gramicidinlike pores

A model calculation is carried out to study the potential energy profile of a sodium ion with several water molecules inside a simplified model of the gramicidin ion channel. The sodium ion is treated as a Lennard-Jones sphere with a point charge at its center. The Barnes polarizable water model is used to mimic the water molecules. A polarizable and deformable gramicidinlike channel is constructed based on the model obtained by Koeppe and Kimura. Potential minima and saddle points are located and the static energy barriers are computed. The potential minima at the two mouths of the channel exhibit an aqueous solvation structure very different from that at any of the interior minima. These sites are approximately 23.6 and 24.4 A apart for binding of a sodium ion and a cesium ion, respectively. Ionic motion from these exterior sites to the first interior minimum requires substantial rearrangement of the waters of solvation; this rearrangement may be the hydration/dehydration step in ionic permeation through the channel. Based on these results, a mechanism by which the sodium ion moves from the exterior binding site to the interior of the channel is proposed. Our model channel accommodates about eight water molecules and the transport of the ion and water within the channel is found to be single file. Results of less extensive calculations for Cs+ and Li+ ions in a channel with or without water are also reported.

Vibrational analysis of peptides, polypeptides and proteins. XXVII. Structure of gramicidin S from normal mode analyses of low-energy conformations

Normal mode calculations have been carried out on three low-energy structures of gramicidin S obtained from conformational energy calculations. When the results on the amide modes are compared with observed bands in the infrared and Raman spectra of crystalline gramicidin S and its N-deuterated derivative, one of the structures is clearly disfavored. Of the other two, one is slightly favored, and it corresponds to the lowest-energy structure obtained from the energy calculations. Spectra from solutions in DMSO and CH3 OH suggest that the molecular conformation is essentially retained in these solvents.

1H-13C selective NOE studies of the decapeptide gramicidin S

The cyclic decapeptide gramicidin S has been used as a model biopolymer to test the reliability of a structural method which is based on a relaxation analysis of heteronuclear selective NOEs. The observation of through-the-space dipolar couplings between intra- and inter residue amide protons and carbonyl carbons, perfectly consistent with the well established peptide solution conformation, confirms the effectiveness of this structural approach. As a corollary of the latter, carbonyl carbon resonances are unequivocally assigned. Moreover, a direct experimental proof of a Orn-NH2----Phe C = O hydrogen bonding is here given.

Temperature dependence of single channel currents and the peptide libration mechanism for ion transport through the gramicidin A transmembrane channel

A study of the temperature dependence of gramicidin A conductance of K+ in diphytanoyllecithin/n-decane membranes shows the plot of In (single channel conductance) as a function of reciprocal temperature to be nonlinear for the most probable set of conductance states. These results are considered in terms of a series of barriers, of the dynamics of channel conformation, vis-a-vis the peptide libration mechanism, and of the effect of lipid viscosity on side chain motions again as affecting the energetics of peptide libration.

Proton relaxation mechanisms and the measurements of r phi, r psi and transannular interproton distances in gramicidin S

Monoselective, Rio(SE), biselective, Rio(i,j), and nonselective proton spin-lattice relaxation rates have been measured for dilute solutions of gramicidin S in dimethyl sulfoxide and used to evaluate cross-relaxation rates (sigma ij = Rio(i,j)-Rio(SE)) and Fi ratios (Fi = Ri(NS)/Rio(SE)). The cross-relaxation parameters, sigma, and Fi ratios measured for backbone gramicidin S protons predict that the same correlation time, tau c = 1.2 X 10(-9)s, modulates all the dipolar proton-proton interactions and that these interactions represent the main source for the proton spin-lattice relaxation process. The larger relaxation rates for amide versus alpha-protons of the backbone are attributed to dipolar relaxation between 14N and its directly bonded protons and is an approximate measure of the extent of this. The intrabackbone proton-proton distances, evaluated from sigma values, were consistent with the antiparallel beta-plated sheet/beta II'-turn conformation previously proposed for gramicidin S in solution.

Ion currents through pores. The roles of diffusion and external access steps in determining the currents through narrow pores

External access steps, which may include restricted aqueous diffusion, are introduced into a kinetic model for ion transport through narrow pores. The conductance-concentration relation and the concentration dependence of the biionic permeability are calculated using two alternative assumptions: (a) access to the mouth of the pore is allowed only when no ion is within the lumen or at either mouth; (b) ions remain at the mouth only very transiently. With either assumption the concentration dependence of the fluxes is the same as in previous treatments in which all steps in access were lumped into a single process. Also as before, the biionic permeability ratio is independent of concentration so long as the lumen is never doubly occupied. For narrow pores, such as those formed by gramicidin A, the slowest external portion of the access process must occur close to the pore's mouth, and thus the region an ion must occupy to gain access is small. As a consequence, the probability of finding an ion within this region is also small. On this basis, it is argued that the second assumption is appropriate for these pores. The kinetic equations that result are identical to those used by Urban, B., S.B. Hladky, and D.A. Haydon (1980, Biochim. Biophys. Acta. 602:331-354).

Shortened analog of the gramicidin A channel argues for the doubly occupied channel as the dominant conducting state

A shortened analog of the gramicidin A transmembrane channel has been synthesized and its transport characterized in planar lipid bilayer membranes. General considerations of a shorter diffusional length and a shorter distance over which the voltage drop occurs (i.e., an increased electric field) would contribute to an increase in single-channel conductance. The finding of a decreased single-channel conductance supports the perspective that the dominant conducting state is the doubly occupied channel wherein distance-dependent repulsion due to the first ion in the channel impedes entry of the second ion in the shorter channel.

The gramicidin A channel: comparison of the energy profiles of Na+, K+ and Cs+. Influence of the flexibility of the ethanolamine end chain on the profiles

The energy profiles for single occupancy by Cs+, K+ and Na+ in the gramicidin A channel assumed to be in a head-to-head beta 6.3 3.3 helical dimeric structure, were computed: (A) allowing complete conformational freedom to the ethanolamine end, (B) constraining it to stay in its intrinsically preferred conformation. Whatever the constraint, both the entrance barrier and the central barrier appear in the order Cs+ less than K+ less than Na+. Introducing the flexibility of the tail modifies appreciably the profiles and the location of the extrema along it.

[Contribution of NMR spectroscopy to the study of structure-function relations of proteins and peptides]

NMR spectroscopy provides a unique means to study molecular conformation, mechanisms of action and structure-function relationships for peptides and proteins in solution under conditions approaching those of their physiological environment. Development of NMR techniques, especially directed to the peptide and protein conformational analysis, is considered under the topics of two-level signal assignment and structural significance of homo- and heteronuclear spin-spin couplings. The results of NMR conformational analysis are presented for solution spatial structure of valinomy cin and gramicidin A antibiotics, honey-bee neurotoxin apamin, scorpion insectotoxin I5A and snake venom neurotoxins of "short" and "long" types. The structure-function relationships are discussed for these biologically active molecules.

Structure and dynamics of ion transport through gramicidin A

Molecular dynamics calculations in which all atoms were allowed to move were performed on a water-filled ion channel of the polypeptide dimer gramicidin A (approximately 600 atoms total) in the head-to-head Urry model conformation. Comparisons were made among nine simulations in which four different ions (lithium, sodium, potassium, and cesium) were each placed at two different locations in the channel as well as a reference simulation with only water present. Each simulation lasted for 5 ps and was carried out at approximately 300 K. The structure and dynamics of the peptide and interior waters were found to depend strongly on the ion tested and upon its location along the pore. Speculations on the solution and diffusion of ions in gramicidin are offered based on the observations in our model that smaller ions tended to lie off axis and to distort the positions of the carbonyl oxygens more to achieve proper solvation and that the monomer-monomer junction was more distortable than the center of the monomer. With the potential energy surface used, the unique properties of the linear chain of interior water molecules were found to be important for optimal solvation of the various ions. Strongly correlated motions persisting over 25 A among the waters in the interior single-file column were observed.

[Acidity and solubility of gramicidin S in water]

The acid-base transformations of the gramicidin S molecule in water were studied. The protonization constants of the antibiotic amino group were calculated by the data of the potentiometric titration and the antibiotic distribution in the system of chloroform-water: K1 1.55 X 10(10), K2 1.38 X 10(6), the logarithm of the distribution coefficient of gramicidin S in the system of chloroform-water (1:1) lg alpha G 4.10. By the same data the constants of water solubility of gramicidin S base (1.02 X 10(7) mol/l), gramicidin S monohydrate (1.06 X 10(-4) mol/l) and gramicidin S dihydrochloride (2.08 X 10(-4) mol/l) were calculated.

Water structure in the Gramicidin A transmembrane channel

The interaction energy and the structure of water molecules either inside the Gramicidin A transmembrane channel or at its two extremities is examined with the use of iso-energy maps and Monte Carlo simulations. The shape of the channel as experienced by water is analyzed in detail. Variations in the hydration structure due to the presence of a sodium ion placed at several positions along the channel are simulated, analyzed and discussed. Preliminary data on Li+ and K+ interacting with Gramicidin A and the system of water molecules are reported. The Gramicidin A atomic coordinates have been taken from Urry's recent papers (Urry, D.W. (1971) Proc. Natl. Acad. Sci. U.S.A. 68, 672-676 and Urry, D.W., Trapane, T.L. and Prasad, K.U. (1982) Int. J. Quant. Chem. Quant. Biol. Symp. 9, 31-40).

Water structure in the Gramicidin A transmembrane channel

The interaction energy and the structure of water molecules either inside the Gramicidin A transmembrane channel or at its two extremities is examined with the use of iso-energy maps and Monte Carlo simulations. The shape of the channel as experienced by water is analyzed in detail. Variations in the hydration structure due to the presence of a Na+ ion placed at several positions along the channel are simulated, analyzed and discussed. Preliminary data on Li+ and K+ interacting with Gramicidin A and the system of water molecules are reported. The Gramicidin A atomic coordinates have been taken from Urry's recent papers.