Originally published as Biophys J. BioFAST on August 31, 2004.
doi:10.1529/biophysj.104.046102
Biophysical Journal 87:3323-3335 (2004)
© 2004 The Biophysical Society
Dynamic Structure of Vesicle-Bound Melittin in a Variety of Lipid Chain Lengths by Solid-State NMR
Shuichi Toraya,
Katsuyuki Nishimura and
Akira Naito
Faculty of Engineering, Yokohama National University, Yokohama, Japan
Correspondence: Address reprint requests to Akira Naito, E-mail address: naito{at}ynu.ac.jp.
 |
ABSTRACT
|
|---|
Solid-state 31P- and 13C-NMR spectra were recorded in melittin-lecithin vesicles composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) or 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Highly ordered magnetic alignments were achieved with the membrane surface parallel to the magnetic field above the gel-to-liquid crystalline phase transition temperature (Tc). Using these magnetically oriented vesicle systems, dynamic structures of melittin bound to the vesicles were investigated by analyzing the 13C anisotropic and isotropic chemical shifts of selectively 13C-labeled carbonyl carbons of melittin under the static and magic-angle spinning conditions. These results indicate that melittin molecules adopt an
-helical structure and laterally diffuse to rotate rapidly around the membrane normal with tilt angles of the N-terminal helix being 33° and 36° and those of the C-terminal helix being 21° and 25° for DLPC and DPPC vesicles, respectively. The rotational-echo double-resonance method was used to measure the interatomic distance between [1-13C]Val8 and [15N]Leu13 to further identify the bending
-helical structure of melittin to possess the interhelical angles of 126° and 119° in DLPC and DPPC membranes, respectively. These analyses further lead to the conclusion that the
-helices of melittin molecules penetrate the hydrophobic cores of the bilayers incompletely as a pseudo-trans-membrane structure and induce fusion and disruption of vesicles.
 |
INTRODUCTION
|
|---|
Melittin is a major component of venom of a honeybee, Apis mellifera, consisting of 26 amino acid residues with the primary structure of GIGAVLKVLTTGLPALISWIKRKRQQ-NH2 (Habermann and Jentsch, 1967
). Melittin disrupts biomembranes and artificial lipid membranes (Habermann, 1972
; Sessa et al., 1969
) and promotes phospholipid hydrolysis catalyzed by phospholipase A2 (Mollay and Kreil, 1974
; Mollay et al., 1976
). Furthermore, melittin induces fusion of phospholipid vesicles (Morgan et al., 1983
; Eytan and Almary, 1983
). It was shown by electron microscopy with a freeze-fracture method that melittin induces fusion of membranes above the gel-to-liquid crystalline phase-transition temperature, Tc, and membrane disruption with forming discoidal membrane fragments whose edges are surrounded by melittin molecules below Tc in a melittin-lecithin membrane system (Dufourcq et al., 1986
). As acyl chain length is increased, the disk tends to be destabilized (Faucon et al., 1995
). Melittin has a voltage-dependent ion channel activity across a lipid bilayer (Tosteson and Tosteson, 1981
; Kempf et al., 1982
), so that melittin was occasionally regarded as a pore-forming peptide by oligomerization.
Three-dimensional structures of melittin have so far been studied under various conditions. In an aqueous medium, monomeric melittin takes random-coil conformation, whereas it mainly adopts an
-helical structure with increasing concentration of NaCl, thereby associating with the peptides to form a tetramer (Talbot et al., 1979
). The helicity and tetramerization of melittin in an aqueous solution depend on the concentration of melittin, the ionic strength, and pH (Bello et al., 1982
; Quay and Condie, 1983
). In a crystalline state, two crystal polymorphs with the space groups of P6122 and C2221 were reported for melittin (Anderson et al., 1980
). X-ray structural analyses showed that melittin adopts a similar
-helical structure and bends at the region of T11G12 with the kink angle of
120° in both of the crystals: the peptides form a tetramer, where the hydrophilic side chains extend mainly toward the outside of the bend and the hydrophobic sides face the center (Terwilliger and Eisenberg, 1982
; Terwilliger et al., 1982
). The kink angle between two helical rods, I2-T11 and L13-Q26, of melittin in methanol, is larger than the kink angles found in the crystals (Bazzo et al., 1988
). In membrane environments, melittin adopts an
-helical structure when it binds to a sodium dodecylsulfate micelle as shown by circular-dichroism spectra (Dawson et al., 1978
; Knöppel et al., 1979
), and binds to dodecylphosphocholine micelles with the
-helical axis parallel to the micelle/water interface (Inagaki et al., 1989
). Transferred nuclear Overhauser enhancement analysis indicated divergent conformations for the region R22-Q26 of melittin bound to lipid bilayers (Okada et al., 1994
).
To reveal the mechanisms of the actions of melittin upon lipid bilayers, it has been discussed that the
-helical axis of the peptide orients perpendicularly (Vogel et al., 1983
; Brauner et al., 1987
) or parallel (Altenbach et al., 1989
) to the bilayer plane depending on the types of experiments. The structure of melittin bound to mechanically oriented ditetradecylphosphatidylcholine (DTPC) membrane has been examined by observing the 13C-NMR spectra. The results indicated that melittin takes trans-membrane helix and the structure is found to be similar to that found in methanol with the interhelical angle of 160° (Smith et al., 1994
). On the other hand, a "toroidal" model was proposed for the pore-forming mechanism, based on information that both the orientations exist for membrane-bound melittin molecules (Yang et al., 2001
). Many attempts have been performed to elucidate the mechanism of the interaction between melittin and the membranes. However, information obtained from most of the spectroscopic studies in the structural biology has been limited to analyze structures of melittin bound to membranes under the conditions where the peptide does not exhibit fusion or disruption activities on membranes partly because the samples were less hydrated. It was shown by optical microscopic observations under the excess hydration condition that giant vesicles are formed in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer systems containing melittin or dynorphin A(117) at a temperature higher than Tc, where the peptides induces membrane fusion (Naito et al., 2000
, 2002
).
In the presence of the strong magnetic field, the membrane plane of the giant vesicle spontaneously orients parallel to the magnetic field above Tc because the spontaneous orientation is stabilized by the interaction energy between anisotropy of the magnetic susceptibility of the phospholipid molecule and the magnetic field. The vesicle thus undergoes elongation along the magnetic field to increase the orientation energy. Since the hydrated membrane is flexible, melittin molecules can have a dynamic structure. Using solid-state NMR with the oriented melittin-DMPC bilayer system, it was elucidated that melittin adopts a bent trans-membrane
-helical structure and laterally diffuses in the membrane with tilting the helical axes to the membrane normal by 30° and 10° for the N- and C-terminal helical axes, respectively. However, the interhelical angle was not determined to be either 160° or 140° above Tc because the signs of the tilt angles could not be determined (Naito et al., 2000
).
In this study, the influence of altering length of the lipid acyl chain on the membrane-bound structure of melittin was investigated under the fusion condition above Tc to understand the action of the melittin-inducing membrane fusion and disruption. We analyzed the dynamic structure of melittin bound to vesicles composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (12:0/12:0-phosphatidylcholine; DLPC) or 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (16:0/16:0-phosphatidylcholine; DPPC) using solid-state 31P- and 13C-NMR spectroscopies. We also report a novel approach to determine a high-resolution dynamic structure of melittin bound to the membranes under the fusion condition (above the Tc) in view of the 13C chemical shift tensors of the carbonyl carbons.
 |
MATERIALS AND METHODS
|
|---|
Sample preparation
Syntheses of 9-fluorenylmethoxycarbonyl (Fmoc) [1-13C]-L-amino acids were carried out with n-(9-fluorenylmethoxycarbonyl)succinimide (Watanabe Chemical Industries, Hiroshima, Japan) and [1-13C]-labeled amino acids (Cambridge Isotope Laboratories, Andover, MA), following the method reported by Paquet (1982)
. Selectively 13C-labeled melittin molecules at the carbonyl carbon of Gly3, Ala4, Val5, Leu16, Ile17, or Ile20 were synthesized using an Applied Biosystems (Foster City, CA) 431A peptide synthesizer with Fmoc solid phase chemistry. The crude peptide yielded by a cleavage reaction of the peptide-resin with Reagent K (Kings et al., 1990
) was purified by a reverse-phase high-performance liquid chromatography.
DLPC and DPPC were purchased from Sigma (St. Louis, MO) and used without further purification. Powder of the total amount of 50 mg comprising melittin and DLPC or DPPC with the peptide/lipid molar ratio of 1:10 was dissolved in 9 ml of chloroform with a small amount of methanol. The solvent was subsequently evaporated in vacuo to prepare a homogeneous film, followed by removal of the residual solvent under high vacuum. A freeze-and-thaw cycle was repeated 30 times after the film was swelled with 300 µl of Tris buffer (20 mM Tris, 100 mM NaCl, pH 7.5). Five-millimeter O.D. glass and zirconia tubes were filled with the sample and sealed hermetically. These samples were used for NMR measurements of the melittin-lecithin bilayer systems in the highly hydrated states.
The hydrated membrane dispersion was allowed to stand in the liquid-crystalline state for 1 h. Subsequently, the sample was frozen rapidly so that membrane disruption was not induced by melittin. This frozen sample was lyophilized to retain the structure in the lipid bilayer. The lyophilized powder was placed into a 5-mm O.D. zirconia tube and subsequently used for measuring the principal values of the 13C chemical shift tensor of the 13C-labeled carbonyl carbon and the 13C-15N interatomic distance in the immobile state.
NMR measurements
The 31P- and 13C-NMR measurements were performed on a Chemagnetics CMX infinity-400 NMR spectrometer (Chemagnetics, Fort Collins, CO) operating at the 31P and 13C resonance frequencies of 161.15 and 100.11 MHz, respectively. In the 31P- and 13C-NMR experiments of the hydrated melittin-lecithin membrane dispersions, the free induction decay signals were obtained after 90°-excitation pulses of 5.0- and 5.5-µs widths under the presence of high-power proton decoupling pulses of 50- and 45-kHz amplitudes with the repetition times of 2 and 5 s, respectively. Lorentzian line broadenings of 60, 100, and 30 Hz were applied to the free-induction decay to obtain 13C-NMR spectra under the static, slow magic-angle spinning (slow MAS, spinning frequency of 100 Hz) and magic-angle spinning (MAS, spinning frequency of 2 kHz) conditions before Fourier transformations, respectively. To determine the 13C chemical shift tensors and the interatomic distances, 13C-NMR spectra of the lyophilized powder samples were measured at 0°C using cross-polarization (CP) and magic-angle spinning (MAS) with the contact time of 1 ms and the spinning frequencies of 2 kHz. The principal values of 13C chemical shift tensors of the carbonyl carbons were determined by comparing the sideband patterns obtained in the CP-MAS experiments with the simulated spectra. 31P and 13C chemical shift values were externally referred to 0 ppm for the phosphorus of 85% H3PO4 and 176.03 ppm for the carboxyl carbon of glycine from that for tetramethylsilane, respectively. The 15N rotational-echo double-resonance (REDOR) spectra were measured using a xy-4 irradiation pulse to compensate the errors of flip angle, off-resonance effect, and fluctuation of RF field for 13C nuclei. The lengths of
-pulses for 13C and 15N nuclei were 12.3 and 13.5 µs, respectively, and the proton decoupling amplitude was 65 kHz. The rotor frequency was controlled to 4000 ±2 Hz. REDOR and full-echo spectra were recorded at various NcTr values from 6 to 24 ms, where Nc and Tr are the rotor cycle and rotor period numbers, respectively. The REDOR differences were evaluated as
 | (1) |
where SREDOR and Sfull echo are the peak intensities of REDOR and full-echo spectra, respectively. These REDOR differences were plotted against the NcTr values to fit the theoretically obtained curves to determine the 13C15N interatomic distance.
 |
RESULTS
|
|---|
Magnetic orientation of melittin-DLPC and melittin-DPPC vesicles
Fig. 1 shows temperature variations of 31P-NMR spectra of the melittin-DLPC and melittin-DPPC bilayers in the hydrated dispersion systems. The 31P-NMR signals for both systems appeared at 0 ppm below Tc of pure lecithins (Tc = 0°C for DLPC, Tc = 41.5°C for DPPC; see Van Dijck et al., 1976
), indicating that the bilayer systems were disrupted to form small discoidal particles as reported by electron microscopy (Dufourcq et al., 1986
). When the temperature is increased to
Tc, powder patterns appeared in the 31P-NMR spectrum of the melittin-DLPC bilayers, indicating that the small particles fuse to form larger vesicles. At temperatures >Tc, signals became sharp and the chemical shifts were the same as those of the perpendicular components of the axially symmetric 31P chemical-shift powder patterns. It was therefore found that the membrane planes of the vesicles composed of respective DLPC and DPPC orient parallel to the static magnetic field above Tc, to form elongated vesicles in the same manner as melittin-DMPC bilayer systems (Naito et al., 2000
). We refer to these magnetically oriented vesicle systems as MOVS.

View larger version (17K):
[in this window]
[in a new window]
|
FIGURE 1 Temperature variations of 31P-NMR spectra of the melittin-DLPC and melittin-DPPC bilayers in the hydrated dispersion systems. 1000 transients were accumulated for each spectrum.
|
|
Temperature variation of 13C-NMR spectra of melittin-DLPC and melittin-DPPC bilayer systems
Fig. 2 shows temperature variations of 13C-NMR spectra of hydrated melittin-DPPC bilayers recorded under the static condition. Melittin molecules were singly 13C-labeled at carbonyl carbons of Ala4 and Leu16 in the N- and C-terminal helices, respectively. Spectra were recorded under the static condition except for those shown at the bottom of Fig. 2. Sanders et al., have assigned 13C-NMR signals of bicelles consisting of DMPC and DHPC in the magnetic field. They found that the signal of CO1 (carbonyl carbon in the 1-position of the glycerole of DMPC) shifts upfield from that of CO2 (carbonyl carbon in the 2-position of the glycerol of DMPC) as the acyl chains orients perpendicularly to the magnetic field (Sanders, 1993
). In the melittin-DPPC bilayer system at 50°C, the 13C-NMR signals observed at 174.0 and 168.0 ppm were therefore attributed to CO2 (denoted
) and CO1 (denoted
) of DPPC, respectively, the same as the case of magnetically oriented bicelle. This result as well as the 31P-NMR spectrum obtained at 50°C indicates that melittin-DPPC membrane planes orient parallel to the magnetic field. The signal intensity of CO1 (167.3 ppm) decreased at 40°C, and CO1 and CO2 signals merged into one peak at 173.8 ppm (denoted *) below 30°C, indicating that the lipid bilayers were disrupted.
At 50°C, relatively narrow signals of [1-13C]Ala4 and [1-13C]Leu16 appeared at 167.7 and 178.5 ppm, respectively. The signals of both [1-13C]Ala4 and [1-13C]Leu16 were broadened at 40°C, and at 35°C, co-existing narrow signals of [1-13C]Ala4 and [1-13C]Leu16 were observed at 177.9 and 175.9 ppm, respectively. At 30°C, the narrow signals of [1-13C]Ala4 and [1-13C]Leu16 appeared at 178.0 and 175.7 ppm, respectively, which were different from the chemical shift values of the melittin molecules bound to the oriented membranes at 50°C. These results indicate that melittin molecules also orient along the magnetic field by strongly binding to the lipid bilayers. In the 13C-NMR spectra of the melittin-DPPC bilayer systems measured at 20°C under the MAS condition (spinning frequency of 2 kHz), the isotropic signals of [1-13C]Ala4 and [1-13C]Leu16 appeared at 177.9 and 175.7 ppm, respectively, as shown in the bottom of Fig. 2. Since the isotropic values under the MAS condition were nearly identical to the chemical-shift values of the narrow peak obtained below 30°C under the static condition, it is apparent that melittin molecules undergo an isotropic motion as a result of formation of small membrane fragments below 30°C. 13C isotropic chemical shift values of [1-13C]Ala4 and [1-13C]Leu16 were 177.3 and 175.8 ppm at 50°C under the MAS condition, respectively (see Table 1). These values are the same as those in the disruption state (below Tc). Therefore, these results suggest that melittin is strongly bound to the membrane fragment in the disruption state with adopting structure similar to that in the vesicle.
Fig. 3 shows temperature variation of 13C-NMR spectra of hydrated melittin-DLPC bilayers recorded under the static condition. Melittin molecules were singly 13C-labeled at the carbonyl carbons of Ala4 and Leu16 in the N- and C-terminal helices, respectively. In the melittin-DLPC bilayer system at 30°C, the signals observed at 174.2 and 168.3 ppm were attributed to CO2 (denoted
) and CO1 (denoted
) of DLPC, respectively, following the results for bicelles (Sanders, 1993
). This result as well as the 31P-NMR spectrum obtained at 30°C indicates that melittin-DLPC membrane planes orient parallel to the magnetic field. For melittin, the signals of the carbonyl carbons of [1-13C]Ala4 and [1-13C]Leu16-melittin molecules appeared at 169.6 and 175.4 ppm at 30°C, respectively. The signals of both [1-13C]Ala4 and [1-13C]Leu16 were broadened as the temperature was decreased. At 0°C, the signals of the lipid appeared at 173.8 ppm, which is identical to the isotropic chemical shift value, although the signals of [1-13C]Ala4 and [1-13C]Leu16 were significantly broadened.
Since the spectra of [1-13C]Ala4 and [1-13C]Leu16-melittin molecules in both the membrane systems behave in a similar manner of motional average on the temperature variations across the Tc under the static condition, both the N- and C-terminal regions of melittin molecules strongly interact with the vesicles and the membrane fragments. It is also noticed that free melittin molecules do not exist in the lecithin-melittin dispersion systems because sharp signals of melittin do not appear in the spectra at temperatures higher than Tc.
13C-NMR spectra of melittin bound to magnetically oriented DLPC and DPPC vesicles in the magnetic field
Fig. 4 shows 13C-NMR spectra of a variety of singly [1-13C]-labeled melittin molecules bound to hydrated DLPC bilayers. All of the 13C-NMR measurements in the left column were performed at 30°C under the static condition where the bilayer planes oriented parallel to the magnetic field. 13C chemical shift values,
obs, of the signals from each carbonyl carbon of the oriented membrane-bound melittin were investigated under the static condition. Thus
obs at Gly3, Ala4, Val5, Leu16, Ile17, and Ile20 were determined as summarized in Table 1. It is noticed that large differences are seen among the values of each residues. The right column of Fig. 4 demonstrates the 13C-NMR spectra at 30°C under the MAS condition with the spinning frequency of 2 kHz. Conformation-dependent 13C isotropic chemical shift values have been well correlated to the secondary structures in model peptides (Saitô and Ando, 1989
). Using the relation, the isotropic chemical shifts,
iso, revealed that the vicinity of Gly3, Ala4, Val5, Leu16, Ile17, and Ile20 adopts
-helical conformations, and hence the overall structure of melittin bound to the DLPC membrane was found to be
-helical (see Table 1).

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 4 13C-NMR spectra of hydrated melittin-DLPC vesicles at 30°C under the static (left) and the MAS (right) conditions. A variety of carbonyl carbons were labeled with 13C nuclei. The symbols for the lipid signals are identical in meaning with those defined in Fig. 2. 790027,852 transients were accumulated for the spectra.
|
|
The left column of Fig. 5 shows 13C-NMR spectra of a variety of singly [1-13C]-labeled melittin molecules bound to the magnetically oriented DPPC bilayers at 50°C. Under the static condition,
obs values at Gly3, Ala4, Val5, Leu16, Ile17, and Ile20 were determined as summarized in Table 1. It is noticed that these values are significantly different from those at respective residues of the melittin in DLPC bilayers. This result indicates that the orientation of melittin helix in DPPC vesicles is slightly different from that in DLPC vesicles. Under the MAS condition,
iso at each residue was nearly equal to the isotropic chemical shift value obtained in the melittin-DLPC bilayer system. It was therefore found that melittin molecules bound to the DPPC bilayers also adopt an
-helical structure over the entire region.

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 5 13C-NMR spectra of hydrated melittin-DPPC vesicles at 50°C under the static (left) and the MAS (right) conditions. A variety of carbonyl carbons were labeled with 13C nuclei. The symbols for the lipid signals are identical in meaning with those defined in Fig. 2. 700030,000 transients were accumulated for the spectra.
|
|
As demonstrated in our previous study, the carbonyl peaks of melittin bound to magnetically oriented membranes coincide with the perpendicular components of axially symmetric powder patterns obtained under the slow MAS (spinning frequency of 100 Hz) condition. This experiment was performed to disturb the magnetic orientation to observe powder patterns in the NMR spectra. The results indicate that the chemical shifts under the static condition are the same as those of the perpendicular components of the axially symmetric powder pattern obtained from the slow MAS experiments. Therefore, it is classified that melittin laterally diffuses in the membrane with rotating about an axis that is perpendicular to the membrane plane, implying that the rotation axis is perpendicular to the magnetic field (Naito et al., 2000
).
In the present experiments, the carbonyl peaks under the oriented conditions in both the DLPC and DPPC membrane systems coincided with the perpendicular components, resonating at 
, of the axially symmetric powder patterns as shown in Fig. 6. Consequently, it was found that melittin molecules bound to DLPC and DPPC membranes rotate about the axis perpendicular to the membrane plane. It is stressed that the significant difference in
obs between the melittin-DLPC and melittin-DPPC bilayer systems at each residue were detected as summarized in Table 1. Since it turned out that
obs is equal to 
, the change of 13C chemical shift anisotropies, 
=
||
obs = 3(
iso
obs), of melittin bound to the differently oriented membranes, actually depends on the orientation of melittin molecules in membranes with the different acyl chain lengths.

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 6 13C-NMR spectra of hydrated [1-13C]Ile17-melittin-DLPC (left) and [1-13C]Val5-melittin DPPC (right) vesicle systems under the static (a and d; oriented), slow MAS (b and e; spinning frequency of 100 Hz), and MAS (c and f; spinning frequency of 2 kHz) conditions. 15,85631,744 transients were accumulated for the spectra.
|
|
13C-NMR spectra of melittin in the lyophilized lecithin bilayer system
To determine the principal values of the 13C chemical shift tensor without motional averaging for each carbonyl carbon of melittin bound to the lecithin bilayer, 13C-NMR spectra of the lyophilized powder samples of the melittin-lecithin bilayer systems were measured using cross-polarization and magic-angle spinning (CP-MAS) with the spinning frequency of 2 kHz (Fig. 7). Since membrane disruption can be induced in the processes of freezing and melting, lyophilization was carried out after rapid freezing by use of the liquid nitrogen. The isotropic chemical shift values,
obtained from the lyophilized powder samples were almost identical with
iso obtained from the hydrated bilayer dispersion systems at the respective carbonyl carbons, implying that the
-helical conformations are retained in the lyophilized powder samples. This suggests that this lyophilization process does not yield significant differences in the principal values of the 13C chemical shift tensors. Chemical-shift sideband patterns were analyzed by fitting with the simulated spectra to determine the 13C chemical shift tensors. Table 1 summarizes the 13C chemical shift values and the anisotropies obtained from this study. It is noticed that substantially different principal values were obtained depending on the amino acid residues. Although the orientations of the principal axes were not determined, it is probable that the orientation of
22 is close to the C=O bond direction as observed in simple peptides (Separovic et al., 1990
).

View larger version (31K):
[in this window]
[in a new window]
|
FIGURE 7 13C-NMR spectra of lyophilized melittin-DLPC and melittin-DPPC bilayer systems using CP-MAS with the spinning frequency of 2 kHz. 478833,628 transients were accumulated for the spectra.
|
|
Measurement of interatomic distance between [1-13C]Val8 and [15N]Leu13
Interatomic distance of carbonyl carbon of Val8 and amide nitrogen of Leu13 was measured by means of the REDOR method (Gullion and Schaefer, 1989
) for the lyophilized powder sample as shown in Fig. 8. In this experiment, we have observed 15N-REDOR and full echo spectra, because the background signals due to natural abundance nuclei can be neglected (Naito et al., 1994
). After the experimentally obtained
S/S0 values were plotted as shown in Fig. 8 b, the best-fit curve was calculated by considering the effect of finite pulse length to obtain the interatomic distance accurately (Naito et al., 1996
). Finally the interatomic distance was determined to be 4.8 ± 0.2 Å. As is discussed later, this value provides information for a bending
-helical structure of melittin bound to lecithin bilayers.

View larger version (14K):
[in this window]
[in a new window]
|
FIGURE 8 (a) 15N-REDOR and full echo spectra of [1-13C]Val8, [15N]Leu-melittin in a lyophilized DPPC bilayer system at NcTr of 18 ms (9000 transients). (b) A plot of S/S0 against NcTr values and the best fit curve indicating the interatomic distance to be 4.8 Å.
|
|
 |
DISCUSSION
|
|---|
Analysis of dynamic structures of melittin molecules bound to magnetically oriented vesicles
It is important to reveal the membrane-bound structures of peptides that exhibit their activities by binding to lipid bilayers, if one is to understand the interactions between peptides and lipids. In solid-state NMR studies, membrane-bound structures of peptides have been mostly investigated using mechanically oriented bilayer systems where the bilayers are aligned between glass plates (Nicholson et al., 1987
). The structure of melittin was analyzed to be a trans-membrane
-helix using mechanically oriented bilayer systems (Smith et al., 1994
). However, the melittin molecule did not exhibit its activity as fusion, disruption, or channel formation on membrane in the systems. We believe that insufficient hydration restricts the mobility of the lipid on the glass plate. Experiments with magnetically oriented vesicles dispersed with excess hydration are thus reasonable to determine the active structure and to understand the mechanism of membrane fusion and disruption by melittin.
13C-NMR signals of the carbonyl carbons of the membrane-bound melittin molecules exhibit the chemical shift anisotropies of
150 ppm when the peptide is immobilized, as shown in Fig. 7 and Table 1. We have proposed that dynamic structures of melittin molecules bound to membranes under the fusion condition can be revealed by analyzing how these large anisotropies are averaged by the motion (Naito et al., 2000
). Since anisotropy of magnetic susceptibility of a lecithin is negative, the acyl chains orient perpendicularly to the magnetic field (Boroske and Helfrich, 1978
). Thus, when giant vesicles are formed as a result of fusion induced by melittin in the magnetic field, they spontaneously orient the membrane planes parallel to the magnetic field by forming elongated vesicles. Therefore, determining the direction of the rotation axis of melittin with respect to the oriented vesicle makes it possible to analyze a high-resolution dynamic structure of the peptide bound to the membrane in the mobile system as described below. We therefore developed an approach of analyses of the anisotropic 13C chemical shift interactions for backbone carbonyl carbons of peptides bound to membranes. As shown in Fig. 9, Euler rotations following the convention of Rose (1957)
were utilized to transform the 13C chemical shift tensors from the principal-axis frame (PAF) to the helical molecular frame (HMF), and consecutively to the molecular rotating frame (MRF) and finally to the laboratory frame (LF) using the three-step transformation for analyses of the dynamic structures of melittin bound to membranes. Assuming that the
22 axis of the carbonyl carbon is parallel to the helical axis, the Z axis, of melittin (Fig. 9 a) for a transformation from the PAF (
11,
22,
33) to the HMF (X, Y, Z), the matrix R1 of the Euler rotation is given by
 | (2) |
where
is a constant and the helical axis is directed from the N- to C-terminus. When the helical axis rotates about the Z' axis at a constant tilt angle of
to the rotation axis (Fig. 9 b), for a transformation from the HMF to the MRF (X', Y', Z'), the matrix R2 of the Euler rotation is given by
 | (3) |
where
varies continuously with a rotational motion at the Z' axis. It is emphasized that the direction of the tilt angle is important to determine the interhelical angle in the case where more than two helices exist in the same molecule. If the rotation axis inclines
to the magnetic field, the z axis (Fig. 9 c), for a transformation from the MRF to the LF (x, y, z) the matrix R3 of the Euler rotation is given by
 | (4) |
Here,
is an arbitrary constant under the static condition, and
=
/2 is chosen accordingly
 | (4') |
Consequently, a rotation matrix of R = R3R2R1 transforms the chemical shift tensor from the PAF to the LF, as
 | (5) |
When the rotation axis inclines 90° to the magnetic field, the observed chemical shift value (
zz)
=
/2 = 
is expressed as
 | (6) |
On the other hand, when the rotation axis inclines 0° to the magnetic field, the observed chemical shift value (
zz)
=0 =
|| is expressed as
 | (7) |
The chemical shift anisotropy, 
=
||
, is obtained by combining Eqs. 6 and 7 as
 | (8) |
Here, functions of
can be averaged over a cycle when the helical axis rotates rapidly about the axis parallel to the membrane normal that corresponds to the Z' axis. In this case, we obtain an expression where the averaged anisotropy depends on the phase angle,
, of the peptide plane about the helical axis, and the tilt angle,
, of the helical axis from the rotation axis as
 | (9) |
where
 | (10) |
Since the C- and N-terminal regions of melittin adopt
-helical structures, the phase angle,
, varies by 100° per consecutive residue in the direction toward the C-terminus. Eq. 9 indicates that
of the carbonyl carbon can be oscillatorily changed as a function of
along the consecutive amino acid sequence in the
-helical region. The amplitude of the oscillation depends on the tilt angle
. We call this behavior a chemical-shift oscillation. This behavior can be compared with the PISA wheel experiments, where the 15N-1H dipolar and 15N chemical shift interactions are the functions of
and
(Marassi and Opella, 2000
; Wang et al., 2000
). To determine the tilt angle,
, and the phase angle,
, for
-helical rods of the N- and C-terminals individually, we take root mean-square deviations (RMSDs) of the 13C chemical shift anisotropies obtained from the three amino acid residues for the N-terminal helical rod (Gly3, Ala4, and Val5) and the C-terminal helical rod (Leu16, Ile17, and Ile20) as
 | (11) |
The least RMSD provides the phase angles,
G3 and
L16, of the peptide plane at Gly3 and Leu16 and the tilt angles,
N and
C, of the helical axes of the N- and C-termini from the membrane-normal, respectively. Once
G3 and
L16 are determined, the other phase angle in each rod is simultaneously determined by considering the structure of
-helix.

View larger version (21K):
[in this window]
[in a new window]
|
FIGURE 9 Euler rotations for description of the rotational motion of melittin bound to the membrane. (a) Rotation from the principal-axis frame (PAF) to the helical molecular frame (HMF). (b) Rotation from the HMF to the molecular rotation frame (MRF). (c) Rotation from the MRF to the laboratory frame (LF).
|
|
We also consider a case where the membrane-bound melittin also rotates about the helical axis. In this case, the Euler angle
'' is set to a time-dependent variable,
, instead of 0 in the rotation matrix transforming the HMF into the MRF, and the X,Y plane rotates rapidly around the helical axis. Because time-averaging of
-dependent terms simultaneously causes functions of
to be averaged,
becomes independent of
. In the 13C-NMR experiments,
exhibited such large changes that the sign was altered between the carbonyl carbons observed at neighboring residues as shown in Table 1. These results imply that the helices do not rotate about the helical axes but rotate about the membrane normal, and
'' must thus be a constant.
Fig. 10 shows RMSD contour plots of the 13C chemical shift anisotropies, where RMSD is plotted against
as abscissa and
as ordinate. Owing to the symmetry property of the Eq. 9, the relations
 | (12) |
are fulfilled. For the N- and C-terminal helices in the melittin-DLPC bilayer systems, the least RMSD indicates
 | (13) |
and
 | (14) |
respectively. For the N- and C-terminal helices in the melittin-DPPC bilayer systems, the least RMSD indicates
 | (15) |
and
 | (16) |
respectively. Since Thr10 cannot form a hydrogen bond to Pro14, and the MAS experiments (Figs. 4 and 5) and our previous study (Naito et al., 2000
) showed that Gly3, Ala4, Val5, Gly12, Ala15, Leu16, Ile17, and Ile20 adopt
-helical conformations, we assume that Thr11 locates as a boundary between Gly1-Thr10 and Gly12-Gln26, and that the two regions adopt typical
-helical structures. The two possible structures of melittin bound to the DPPC vesicles were shown in Fig. 11 out of 16 possible structures. In the case of extended structure shown in Fig. 11 a, the 13C15N interatomic distance between [1-13C]Val8 and [15N]Leu13 is expected to be longer than 6.9 Å. On the other hand, this distance is expected to be 4.8 Å in the case of bend structure shown in Fig. 11 b. In our previous report, NMR experiments could not distinguish such two structures as shown in Fig. 10, a and b (Naito et al., 2000
). We have determined the interatomic distance between [1-13C]Val8 and [15N]Leu13 for melittin bound to lyophilized DPPC bilayers by means of a REDOR method. The distance was determined to be 4.8 ± 0.2 Å, clearly indicating that the structure with
 | (17) |
as illustrated in Fig. 11 b, is the actual structure and the other possible structures were ruled out. When melittin adopts the structure shown in Fig. 11 b, it agrees well with the structure found by the x-ray diffraction studies (Terwilliger and Eisenberg, 1982
; Terwilliger et al., 1982
). Namely the carbonyl group at Val8 can form hydrogen bond to the amide group at Leu13, and this allows the molecule to behave as a rigid body in membrane. Also for the case of dynamic structure of melittin bound to the DLPC membrane, we finally obtain
 | (18) |
for the actual structure.

View larger version (41K):
[in this window]
[in a new window]
|
FIGURE 10 Contour plots of RMSD of experimental and theoretical 13C chemical shift anisotropies. (a, c) Contour plots for the N-terminal helix of melittin bound to DLPC and DPPC vesicles, respectively. (b, d) Contour plots for the C-terminal helix of melittin bound to DLPC and DPPC vesicles, respectively. RMSD values increase from blue to red.
|
|

View larger version (36K):
[in this window]
[in a new window]
|
FIGURE 11 Schematic representations of the possible structures of melittin bound to the DPPC vesicle from the analyses of the 13C chemical shift anisotropies. (a) ( G3, N) = (+76°, 36°) and ( L16, C) = (82°, 25°). (b) ( G3, N) = (+76°, 36°) and ( L16, C) = (82°, +25°). Structure b is proved to be the actual structure based on the interatomic distance of 4.8 ± 0.2 Å, between [1-13C]Val8 and [15N]Leu13 of melittin bound to the DPPC bilayers, determined by the REDOR measurements. The helical wheel representations were illustrated in the HMF, which is identical with that defined in Fig. 9 b. The Z' axis is the rotation axis of a melittin molecule, which is parallel to the membrane normal.
|
|
Interestingly, alteration of the lipid acyl-chain length brought about significant differences in the tilt angles,
N and
C, rather than the phase angles,
G3 and
L16. Since Gly12 of melittin bound to the DMPC membrane adopts a typical
-helical conformation under a fusion condition (Naito et al., 2000
), the kink angle of the helical molecule is attributed to the torsion angles of predominantly the vicinity of Thr11 that is the linking position of the two helical rods. The interhelical angle of melittin changed from 126° ± 8° in DLPC to 119° ± 6° in DPPC through the torsion angles of Thr11 as a result of altering the acyl groups of the lipid forming the bilayer from lauroyl to palmitoyl. We could detect this small change as significant differences between the chemical shift values,
obs, observed under the oriented conditions as summarized in Table 1. These angles are slightly smaller than those of melittin in the lyophilized powder and the hydrated gel phase of DTPC (Lam et al., 2001
). Conversely, the experiments indicate that our approach, as proposed in this study, is useful for analyses of structures and dynamics of peptides bound to dynamic membranes. It is also emphasized that this approach can determine not only the structure of membrane-bound molecules but also the orientation with respect to the membrane.
Mechanism of membrane fusion induced by melittin
It is demonstrated that alteration of lipid-acyl chain length causes a slight change in the interhelical angle of melittin (126° ± 8° and 119° ± 6° for DLPC and DPPC, respectively). Lis et al. (1982)
reported that an increase of acyl chain length of lipid molecule increases the bilayer thickness in the liquid crystalline states (Lis et al., 1982
). It is of interest to note that distances between the termini of melittin are estimated to be
34 and
33 Å in the DLPC and DPPC bilayer systems, respectively, where melittin adopts the helical structure with 5.4 Å pitch as shown in Fig. 11 b. On the other hand, bilayer thicknesses of DLPC and DPPC bilayers are 30.0 and 34.2 Å, respectively (Lis et al., 1982
). It is therefore possible for melittin molecules to take trans-membrane structures in DLPC and DPPC bilayers. However, the hydrophilic residues of K21-R24 may make one cycle of
-helix and locate in the hydrophilic region of lecithin bilayers by orienting parallel to the membrane plane as shown in Fig. 11 b. Tryptophan residue has also been reported to show an interfacial anchoring property (de Planque et al., 2003
). Thus, W19 may locate itself at the interfacial regions of lecithin bilayers. It is therefore reasonable to predict that the consecutive basic residues are located in the polar region. Consequently, the N-terminus of melittin cannot reach to the water-lipid interface, penetrating the hydrophobic core of the bilayers as a pseudo-trans-membrane structure. As length of the acyl chain increases, interactions between the hydrophobic bulky side chains of such residues as Ala, Ile, Leu, and Val inside of the bend and the hydrophobic core of lipid bilayer increases to cause the slight change in the interhelical angle. It is discussed that the incomplete insertion of the N-terminal
-helix of melittin may cause a great disorder in the lipid bilayer surface on the N-terminal side. The surface disorder induces lipid mixing between adjacent vesicles, thereby membrane fusion occurs.
 |
CONCLUSIONS
|
|---|
It was demonstrated that the magnetically oriented vesicle systems (MOVS) were formed in the melittin-lecithin bilayer systems. Using MOVS, dynamic structure of melittin bound to the lipid bilayer was found as a pseudo-trans-membrane bend structure, and melittin molecules rotate rapidly about the axis parallel to the bilayer normal. Furthermore, we could experimentally detect the small change in the tilt angles of the helical axes of the N- and C-terminal helices from the membrane normal as significant differences between the chemical shift values observed under the conditions where the membrane plane spontaneously orient parallel to the magnetic field. Consequently, detailed analyses of 13C chemical shift interactions of carbonyl carbons revealed the interhelical angles of melittin molecules bound to the DLPC and DPPC vesicles to be 126° and 119°, respectively. Although there are certain differences in the dynamic structures of melittin molecules bound to membranes depending on the membrane systems, melittin exhibits common activities on DLPC, DMPC, and DPPC bilayers such as membrane fusion (above Tc) and disruption (below Tc). In this study, the membrane fusion mechanism can be attributed to the fact that the melittin strongly binds to the vesicle with the pseudo-trans-membrane structure and orientation, which causes disorder of the lipid bilayer surface.
 |
ACKNOWLEDGEMENTS
|
|---|
This work was supported, in part, by a Grant-in-Aid for Scientific Research on Priority Areas (13024263) from the Ministry of Culture, Sports, Science and Technology of Japan.
Submitted on May 21, 2004;
accepted for publication August 20, 2004.
 |
REFERENCES
|
|---|
Altenbach, C., W. Froncisz, J. S. Hyde, and W. L. Hubbell. 1989. Conformation of spin-labeled melittin at membrane surfaces investigated by pulse saturation recovery and continuous wave power saturation electron paramagnetic resonance. Biophys. J. 56:11831191.[Abstract/Free Full Text]
Anderson, D., T. C. Terwilliger, W. Wickner, and D. Eisenberg. 1980. Melittin forms crystals which are suitable for high resolution x-ray structural analysis and which reveal a molecular twofold axis of symmetry. J. Biol. Chem. 255:25782582.[Abstract/Free Full Text]
Bazzo, R., M. J. Tappin, A. Pastore, T. S. Harvey, J. A. Carver, and I. D. Campbell. 1988. The structure of melittin. A 1H-NMR study in methanol. Eur. J. Biochem. 173:139146.[Medline]
Bello, J., H. R. Bello, and E. Granados. 1982. Conformation and aggregation of melittin: dependence on pH and concentration. Biochemistry. 21:461465.[CrossRef][Medline]
Boroske, E., and W. Helfrich. 1978. Magnetic anisotropy of egg lecithin membranes. Biophys. J. 24:863868.[Abstract/Free Full Text]
Brauner, J. W., R. Mendelsohn, and F. G. Prendergast. 1987. Attenuated total reflectance Fourier transform infrared studies of the interaction of melittin, two fragments of melittin, and
-hemolysin with phosphatidylcholines. Biochemistry. 26:81518158.[CrossRef][Medline]
Dawson, C. R., A. F. Drake, J. Helliwell, and R. C. Hider. 1978. The interaction of bee melittin with lipid bilayer membranes. Biochim. Biophys. Acta. 510:7586.[Medline]
de Planque, M. R. R., B. B. Bonev, J. A. A. Demmers, D. V. Greathouse, G. E. Koeppe, II, F. Separovic, A. Watts, and J. A. Killian. 2003. Interfacial anchor properties of tryptophan residues in trans-membrane peptides can dominate over hydrophobic matching effects in peptide-lipid interactions. Biochemistry. 47:53415348.
Dufourcq, J., J.-F. Faucon, G. Fourche, J.-L. Dasseux, M. Le Maire, and T. Gulik-Krzywicki. 1986. Morphological changes of phosphatidylcholine bilayers induced by melittin: vesicularization, fusion, discoidal particles. Biochim. Biophys. Acta. 859:3348.[Medline]
Eytan, G. D., and T. Almary. 1983. Melittin-induced fusion of acidic liposomes. FEBS Lett. 156:2932.[CrossRef][Medline]
Faucon, J.-F., J.-M. Bonmatin, J. Dufourcq, and E. J. Dufourc. 1995. Acyl-chain-length dependence in the stability of melittin-phosphatidylcholine complexes. A light scattering and 31P-NMR study. Biochim. Biophys. Acta. 1234:235243.[Medline]
Gullion, T., and J. Schaefer. 1989. Rotational echo double-resonance NMR. J. Magn. Reson. 81:196200.
Habermann, E. 1972. Bee and wasp venoms. Science. 177:314322.[Free Full Text]
Habermann, E., and J. Jentsch. 1967. Sequenzanalyse des melittins aus den tryptischen und peptischen spaltstücken. Hoppe-Seyler's Z. Physiol. Chem. 348:3750.
Inagaki, F., I. Shimada, K. Kawaguchi, M. Hirano, I. Terasawa, T. Ikura, and N. G
. 1989. Structure of melittin bound to perdeuterated dodecylphosphocholine micelles as studied by two-dimensional NMR and distance geometry calculations. Biochemistry. 28:59855991.[CrossRef]
Kempf, C., R. D. Klausner, J. N. Weinstein, J. Van Renswoude, M. Pincus, and R. Blumenthal. 1982. Voltage-dependent trans-bilayer orientation of melittin. J. Biol. Chem. 257:24692476.[Free Full Text]
Kings, D. S., C. G. Fields, and G. B. Fields. 1990. A cleavage method which minimizes side reactions following Fmoc solid phase peptide synthesis. Int. J. Pept. Protein Res. 36:255266.[Medline]
Knöppel, E., D. Eisenberg, and W. Wickner. 1979. Interactions of melittin, a protein model, with detergents. Biochemistry. 18:41774181.[CrossRef][Medline]
Lam, Y.-H., S. R. Wassall, C. I. Morton, R. Smith, and F. Separovic. 2001. Solid-state NMR structure determination of melittin in a lipid environment. Biophys. J. 81:27522764.[Abstract/Free Full Text]
Lis, L. J., M. McAlister, N. Fuller, R. P. Rand, and V. A. Parsegian. 1982. Interactions between neutral phospholipid bilayer membranes. Biophys. J. 37:657666.[Abstract/Free Full Text]
Marassi, F. M., and S. J. Opella. 2000. A solid-state NMR index of helical membrane protein structure and topology. J. Magn. Reson. 144:150155.[CrossRef][Medline]
Mollay, C., and G. Kreil. 1974. Enhancement of bee venom phospholipase A2 activity by melittin, direct lytic factor from cobra venom and polymyxin B. FEBS Lett. 46:141144.[CrossRef][Medline]
Mollay, C., G. Kreil, and H. Berger. 1976. Action of phospholipases on the cytoplasmic membrane of Escherichia coli. Biochim. Biophys. Acta. 426:317324.[Medline]
Morgan, C. G., H. Williamson, S. Fuller, and B. Hudson. 1983. Melittin induces fusion of unilamellar phospholipid vesicles. Biochim. Biophys. Acta. 732:668674.[Medline]
Naito, A., T. Nagao, K. Norisada, T. Mizuno, S. Tuzi, and H. Saitô. 2000. Conformation and dynamics of melittin bound to magnetically oriented lipid bilayers by solid-state 31P- and 13C-NMR spectroscopy. Biophys. J. 78:24052417.[Abstract/Free Full Text]
Naito, A., T. Nagao, M. Obata, Y. Shindo, M. Okamoto, S. Yokoyama, S. Tuzi, and H. Saitô. 2002. Dynorphin-induced magnetic ordering in lipid bilayers as studied by 31P-NMR spectroscopy. Biochim. Biophys. Acta. 1558:3444.[Medline]
Naito, A., K. Nishimura, S. Kimura, M. Aida, N. Yasuoka, S. Tuzi, and H. Saitô. 1996. Determination of the three-dimensional structure of a new crystalline form of n-acetyl-Pro-Gly-Phe as revealed by 13C REDOR, x-ray diffraction, and molecular dynamics calculation. J. Phys. Chem. B. 100:1499515004.
Naito, A., K. Nishimura, S. Tuzi, and H. Saitô. 1994. Inter- and intra-molecular contributions of neighboring dipolar pairs to the precise determination of interatomic distances in a simple [13C, 15N]-peptide by 13C, 15N-REDOR NMR spectroscopy. Chem. Phys. Lett. 229:506511.[CrossRef]
Nicholson, L. K., F. Moll, T. E. Mixon, P. V. LoGrasso, J. C. Lay, and T. A. Cross. 1987. Solid-state 15N-NMR of oriented lipid bilayer bound gramicidin A'. Biochemistry. 26:66216626.[CrossRef][Medline]
Okada, A., K. Wakamatsu, T. Miyazawa, and T. Higashijima. 1994. Vesicle bound conformation of melittin: transferred nuclear Overhauser enhancement analysis in the presence of perdeuterated phosphatidylcholine vesicles. Biochemistry. 33:94389446.[CrossRef][Medline]
Paquet, A. 1982. Introduction of 9-fluorenylmethoxycarbonyl, trichloroethoxycarbonyl, and benzyloxycarbonyl amine protecting groups into O-unprotected hydroxyamino acids using succinimidyl carbonates. Can. J. Chem. 60:976980.
Quay, S. C., and C. C. Condie. 1983. Conformational studies of aqueous melittin: thermodynamics parameters of the monomer-tetramer self-association reaction. Biochemistry. 22:695700.[CrossRef][Medline]
Rose, M. E. 1957. Elementary Theory of Angular Momentum. John Wiley and Sons, New York.
Saitô, H., and I. Ando. 1989. High-resolution solid-state NMR studies of synthetic and biological macromolecules. Annu. Rep. NMR Spectrosc. 21:209290.
Sanders, C. R. 1993. Solid state 13C-NMR of unlabeled phosphatidylcholine bilayers: spectral assignments and measurement of carbon-phosphorus dipolar couplings and 13C chemical shift anisotropies. Biophys. J. 64:171181.[Abstract/Free Full Text]
Separovic, F., R. Smith, C. S. Yannoni, and B. A. Cornell. 1990. Molecular sequence effect on the 13C carbonyl chemical shift shielding tensor. J. Am. Chem. Soc. 112:83248328.[CrossRef]
Sessa, G., J. H. Freer, G. Colaccico, and G. Weissmann. 1969. Interaction of a lytic polypeptide, melittin, with lipid membrane systems. J. Biol. Chem. 244:35753582.[Abstract/Free Full Text]
Smith, R., F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis. 1994. Structure and orientation of the pore-forming peptide, melittin, in lipid bilayers. J. Mol. Biol. 241:456466.[CrossRef][Medline]
Talbot, J. C., J. Dufourcq, J. de Bony, J. F. Faucon, and C. Lussan. 1979. Conformational change and self-association of monomeric melittin. FEBS Lett. 102:191193.[CrossRef][Medline]
Terwilliger, T. C., and D. Eisenberg. 1982. The structure of melittin. J. Biol. Chem. 257:60106015.[Abstract/Free Full Text]
Terwilliger, T. C., L. Weissman, and D. Eisenberg. 1982. The structure of melittin in the form I crystals and its implication for melittin's lytic and surface activities. Biophys. J. 37:353361.[Abstract/Free Full Text]
Tosteson, M. T., and D. C. Tosteson. 1981. Melittin forms channels in lipid bilayers. Biophys. J. 36:109116.[Abstract/Free Full Text]
Van Dijck, P. W. M., B. De Kruijff, L. L. M. Van Deenen, J. De Gier, and R. A. Demel. 1976. The preference of cholesterol for phosphatidylcholine in mixed phosphatidylcholine-phosphatidylethanolamine bilayers. Biochim. Biophys. Acta. 455:576587.[Medline]