Membrane-bound ARF1 peptide: interpretation of neutron diffraction data by molecular dynamics simulation methods.

Adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We have previously presented four possible conformations of the membrane bound, human ARF1 N-terminal peptide in planar lipid bilayers of DOPC and DOPG (7:3 molar ratio), determined from lamellar neutron diffraction and circular dichroism data. In this paper we analyse the four possible conformations by molecular dynamics simulations. The aim of these simulations was to use MD to distinguish which of the four possible membrane bound structures was the most likely. The most likely conformation was determined according to the following criteria: (a) location of label positions on the peptide in relation to the bilayer, (b) lowest mean square displacement from the initial structure, (c) lowest system energy, (d) most peptide-lipid headgroup hydrogen bonding, (e) analysis of phi/psi angles of the peptide. These findings demonstrate the application of molecular dynamics simulations to explore neutron diffraction data.

Finite-size effects do not reduce the repeat spacing of phospholipid multibilayer stacks on a rigid substrate.

Finite-size effects in stacks of phospholipid bilayers, in the fluid L alpha phase, are investigated using samples oriented on silicon substrates. Recently in this journal, such effects have been suggested as the probable cause of reduced lamellar repeat spacings in very thin samples made up of a few (<10) bilayers. Our systematic studies on samples of different thicknesses do not support this conclusion. At full hydration all samples are found to have the same repeat spacing, irrespective of their thickness. At lower hydrations, on the other hand, very thin samples, consisting of only a few bilayers, have a slightly larger spacing.

Relationship between the unbinding and main transition temperatures of phospholipid bilayers under pressure.

Using neutron diffraction and a specially constructed high pressure cell suitable for aligned multibilayer systems, we have studied, as a function of pressure, the much observed anomalous swelling regime in dimyristoyl- and dilauroyl-phosphatidylcholine bilayers, DMPC and DLPC, respectively. We have also reanalyzed data from a number of previously published experiments and have arrived at the following conclusions. (a). The power law behavior describing anomalous swelling is preserved in all PC bilayers up to a hydrostatic pressure of 240 MPa. (b). As a function of increasing pressure there is a concomitant decrease in the anomalous swelling of DMPC bilayers. (c). For PC lipids with hydrocarbon chains >or=13 carbons the theoretical unbinding transition temperature T small star, filled is coupled to the main gel-to-liquid crystalline transition temperature T(M). (d). DLPC is intrinsically different from the other lipids studied in that its T small star, filled is not coupled to T(M). (e). For DLPC bilayers we predict a hydrostatic pressure (>290 MPa) where unbinding may occur.

Spontaneously formed monodisperse biomimetic unilamellar vesicles: the effect of charge, dilution, and time.

Using small-angle neutron scattering and dynamic light scattering, we have constructed partial structural phase diagrams of lipid mixtures composed of the phosphatidylcholines dimyristoyl and dihexanoyl doped with calcium ions (Ca2+) and/or the negatively charged lipid, dimyristoyl phosphatidylglycerol (DMPG). For dilute solutions (lipid concentration < or =1 wt %), spontaneously forming unilamellar vesicles (ULVs) were found, and their polydispersity was determined to be approximately 20%. The stability of the Ca2+- or DMPG-doped ULVs was monitored over a period of 4 days and their structural parameters (e.g., average outer radius, ) were found to be insensitive to the lipid concentration (Clp). However, doping the dimyristoyl/dihexanoyl system with both Ca2+ and DMPG resulted in ULVs whose was found to be Clp dependent. The of DMPG-doped ULVs remained unchanged over an extended period of time (at least 4 days), a good indication of their stability.

Finite-size effects in biomimetic smectic films.

Thin stacks of lipid multibilayers supported on rigid silicon and mica substrates are found to exhibit finite-size effects. Using neutron diffraction we find that the repeat spacing (d) of stacks containing up to a few tens of bilayers depends on their thickness (D), with d increasing with decreasing D. Differences in d are larger in the low-temperature Lbeta' phase consisting of rigid bilayers than in the high-temperature Lalpha phase where the bilayers are more flexible. Various scenarios that may be responsible for this counterintuitive observation are discussed.

Concentration-independent spontaneously forming biomimetric vesicles.

In this Letter we present small-angle neutron scattering data from a biomimetic system composed of the phospholipids dimyristoyl and dihexanoyl phosphorylcholine (DMPC and DHPC, respectively). Doping DMPC-DHPC multilamellar vesicles with either the negatively charged lipid dimyristoyl phosphorylglycerol (DMPG, net charge -1) or the divalent cation, calcium (Ca2+), leads to the spontaneous formation of energetically stabilized monodisperse unilamellar vesicles whose radii are concentration independent and in contrast with previous experimental observations.

Molecular dynamics simulations of a mixed DOPC/DOPG bilayer.

We have constructed a mixed dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylglycerol bilayer (DOPG) bilayer utilizing MD simulations. The aim was to develop an explicit molecular model of biological membranes as a complementary technique to neutron diffraction studies that are well established within the group. A monolayer was constructed by taking a previously customized PDB file of each molecule and arranging them in a seven rows of ten molecules and duplicated and rotated to form a bilayer. The 140-molecule bilayer contained 98 DOPC molecules and 42 DOPG molecules, in a 7:3 ratio in favour of DOPC. Sodium counter ions were placed near the phosphate moiety of DOPG to counteract the negative charge of DOPG. This was representative of the lipid ratio in a sample used for neutron diffraction. The MD package GROMACS was used for confining the bilayer in a triclinic box, adding Simple Polar Charge water molecules, energy minimization (EM). The bilayer/solvent system was subjected to EM using the steepest descent method to nullify bad contacts and reduce the potential energy of the system. Subsequent MD simulation using an initial NVT (constant number of particles, volume and temperature) for a 20 ps MD run followed by a NPT (constant number of particles, pressure and temperature) was performed. Structural parameters including volume of lipid, area of lipid, order parameter of the fatty acyl carbons and electron density profiles generated by the MD simulation were verified with values obtained from experimental data of DOPC, as there are no comparable experimental data available for the mixed bilayer.

Inhibitors can arrest the membrane activity of human islet amyloid polypeptide independently of amyloid formation.

Human islet amyloid polypeptide (hIAPP), co-secreted with insulin from pancreatic beta cells, misfolds to form amyloid deposits in non-insulin-dependent diabetes mellitus (NIDDM). Like many amyloidogenic proteins, hIAPP is membrane-active: this may be significant in the pathogenesis of NIDDM. Non-fibrillar hIAPP induces electrical and physical breakdown in planar lipid bilayers, and IAPP inserts spontaneously into lipid monolayers, markedly increasing their surface area and producing Brewster angle microscopy reflectance changes. Congo red inhibits these activities, and they are completely arrested by rifampicin, despite continued amyloid formation. Our results support the idea that non-fibrillar IAPP is membrane-active, and may have implications for therapy and for structural studies of membrane-active amyloid.

Barrel-stave model or toroidal model? A case study on melittin pores.

Transmembrane pores induced by amphiphilic peptides, including melittin, are often modeled with the barrel-stave model after the alamethicin pore. We examine this assumption on melittin by using two methods, oriented circular dichroism (OCD) for detecting the orientation of melittin helix and neutron scattering for detecting transmembrane pores. OCD spectra of melittin were systematically measured. Melittin can orient either perpendicularly or parallel to a lipid bilayer, depending on the physical condition and the composition of the bilayer. Transmembrane pores were detected when the helices oriented perpendicularly to the plane of the bilayers, not when the helices oriented parallel to the bilayers. The evidence that led to the barrel-stave model for alamethicin and that to the toroidal model for magainin were reviewed. The properties of melittin pores are closely similar to that of magainin but unlike that of alamethicin. We conclude that, among naturally produced peptides that we have investigated, only alamethicin conforms to the barrel-stave model. Other peptides, including magainins, melittin and protegrins, all appear to induce transmembrane pores that conform to the toroidal model in which the lipid monolayer bends continuously through the pore so that the water core is lined by both the peptides and the lipid headgroups.

Oblique membrane insertion of viral fusion peptide probed by neutron diffraction.

Fusion peptides mimic the membrane fusion activities of the larger viral proteins from which they derive their sequences. A possible mode of activity involves their oblique insertion into lipid bilayers, causing membrane disruption by promoting highly curved hemifusion intermediates, leading to fusion. We have determined the location and orientation of the simian immunodeficiency virus (SIV) fusion peptide in planar lipid bilayers using neutron lamellar diffraction. The helical axis of the peptide adopts an angle of 55 degrees relative to the membrane normal, while it positions itself nearest the lipid bilayer surface. This is the first direct observation of the structural interaction between a fusion peptide and a phospholipid bilayer.

Membrane thinning effect of the beta-sheet antimicrobial protegrin.

Lipid bilayers containing the antimicrobial peptide protegrin-1 (PG-1) were studied by lamellar X-ray diffraction. Previously, we have shown that the peptide exists in two distinct states when associated with lipid bilayers depending on the peptide concentration [Heller, W. T., Waring, A. J., Lehrer, R. I., and Huang, H. W. (1998) Biochemistry 37, 17331-17338]. For concentrations below a lipid-dependent threshold, PG-1 exhibits a unique oriented circular dichroism spectrum called the S state. X-ray experiments show that in this state PG-1 decreases the thickness of the lipid bilayer in proportion to the peptide concentration, similar to alamethicin's membrane thinning effect. This indicates that the S state is adsorbed in the headgroup region of the lipid bilayer, where the peptide is in an inactive state. For PG-1 above the threshold concentration, X-ray diffraction shows that the interaction between the peptide and the bilayer changes significantly. These results suggest that PG-1 has the same concentration-gated mechanism of action as alamethicin.

Supramolecular structures of peptide assemblies in membranes by neutron off-plane scattering: method of analysis.

In a previous paper (Yang et al., Biophys. J. 75:641-645, 1998), we showed a simple, efficient method of recording the diffraction patterns of supramolecular peptide assemblies in membranes where the samples were prepared in the form of oriented multilayers. Here we develop a method of analysis based on the diffraction theory of two-dimensional liquids. Gramicidin was used as a prototype model because its pore structure in membrane in known. At full hydration, the diffraction patterns of alamethicin and magainin are similar to gramicidin except in the scale of q (the momentum transfer of scattering), clearly indicating that both alamethicin and magainin form pores in membranes but of different sizes. When the hydration of the multilayer samples was decreased while the bilayers were still fluid, the in-plane positions of the membrane pores became correlated from one bilayer to the next. We believe that this is a new manifestation of the hydration force. The effect is most prominent in magainin patterns, which are used to demonstrate the method of analysis. When magainin samples were further dehydrated or cooled, the liquid-like diffraction turned into crystal-like patterns. This discovery points to the possibility of investigating the supramolecular structures with high-order diffraction.

Theoretical analysis of hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

We present a quantitative analysis of the effects of hydrophobic matching and membrane-mediated protein-protein interactions exhibited by gramicidin embedded in dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) bilayers (Harroun et al., 1999. Biophys. J. 76:937-945). Incorporating gramicidin, at 1:10 peptide/lipid molar ratio, decreases the phosphate-to-phosphate (PtP) peak separation in the DMPC bilayer from 35.3 A without gramicidin to 32.7 A. In contrast, the same molar ratio of gramicidin in DLPC increases the PtP from 30.8 A to 32.1 A. Concurrently, x-ray in-plane scattering showed that the most probable nearest-neighbor separation between gramicidin channels was 26.8 A in DLPC, but reduced to 23.3 A in DMPC. In this paper we review the idea of hydrophobic matching in which the lipid bilayer deforms to match the hydrophobic surface of the embedded proteins. We use a simple elasticity theory, including thickness compression, tension, and splay terms to describe the membrane deformation. The energy of membrane deformation is compared with the energy cost of hydrophobic mismatch. We discuss the boundary conditions between a gramicidin channel and the lipid bilayer. We used a numerical method to solve the problem of membrane deformation profile in the presence of a high density of gramicidin channels and ran computer simulations of 81 gramicidin channels to find the equilibrium distributions of the channels in the plane of the bilayer. The simulations contain four parameters: bilayer thickness compressibility 1/B, bilayer bending rigidity Kc, the channel-bilayer mismatch Do, and the slope of the interface at the lipid-protein boundary s. B, Kc, and Do were experimentally measured; the only free parameter is s. The value of s is determined by the requirement that the theory produces the experimental values of bilayer thinning by gramicidin and the shift in the peak position of the in-plane scattering due to membrane-mediated channel-channel interactions. We show that both hydrophobic matching and membrane-mediated interactions can be understood by the simple elasticity theory.

Experimental evidence for hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

Hydrophobic matching, in which transmembrane proteins cause the surrounding lipid bilayer to adjust its hydrocarbon thickness to match the length of the hydrophobic surface of the protein, is a commonly accepted idea in membrane biophysics. To test this idea, gramicidin (gD) was embedded in 1, 2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1, 2-myristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers at the peptide/lipid molar ratio of 1:10. Circular dichroism (CD) was measured to ensure that the gramicidin was in the beta6.3 helix form. The bilayer thickness (the phosphate-to-phosphate distance, or PtP) was measured by x-ray lamellar diffraction. In the Lalpha phase near full hydration, PtP is 30.8 A for pure DLPC, 32.1 A for the DLPC/gD mixture, 35.3 A for pure DMPC, and 32.7 A for the DMPC/gD mixture. Gramicidin apparently stretches DLPC and thins DMPC toward a common thickness as expected by hydrophobic matching. Concurrently, gramicidin-gramicidin correlations were measured by x-ray in-plane scattering. In the fluid phase, the gramicidin-gramicidin nearest-neighbor separation is 26.8 A in DLPC, but shortens to 23.3 A in DMPC. These experiments confirm the conjecture that when proteins are embedded in a membrane, hydrophobic matching creates a strain field in the lipid bilayer that in turn gives rise to a membrane-mediated attractive potential between proteins.

Neutron off-plane scattering of aligned membranes. I. Method Of measurement.

We describe a method of measuring neutron scattering of aligned membranes with the momentum transfer oriented parallel or partly perpendicular to the plane of the membranes. The method obtains the complete information for the structures within fluid membranes obtainable by scattering. Data from alamethicin- and magainin-induced pores are presented. Although the in-plane scattering curves of these two peptides are similar to each other, their off-plane scattering patterns are strikingly distinct. Magainin pores exhibit intermembrane correlations.

Effect of changing the size of lipid headgroup on peptide insertion into membranes.

Adsorption of amphiphilic peptides to the headgroup region of a lipid bilayer is a common mode of protein-membrane interactions. Previous studies have shown that adsorption causes membrane thinning. The degree of the thinning depends on the degree of the lateral expansion caused by the peptide adsorption. If this simple molecular mechanism is correct, the degree of lateral expansion and consequently the membrane thinning should depend on the size of the headgroup relative to the cross section of the hydrocarbon chains. Previously we have established the connection between the alamethicin insertion transition and the membrane thinning effect. In this paper we use oriented circular dichroism to study the effect of varying the size of the headgroup, while maintaining a constant cross section of the lipid chains, on the insertion transition. A simple quantitative prediction agrees very well with the experiment.

Membrane pores induced by magainin.

Magainin, found in the skin of Xenopus laevis, belongs to a broad class of antimicrobial peptides which kill bacteria by permeabilizing the cytoplasmic membrane but do not lyse eukaryotic cells. The 23-residue peptide has been shown to form an amphiphilic helix when associated with membranes. However, its molecular mechanism of action has been controversial. Oriented circular dichroism has detected helical magainin oriented perpendicular to the plane of the membrane at high peptide concentrations, but Raman, fluorescence, differential scanning calorimetry, and NMR all indicate that the peptide is associated with the head groups of the lipid bilayer. Here we show that neutron in-plane scattering detects pores formed by magainin 2 in membranes only when a substantial fraction of the peptide is oriented perpendicular to the membrane. The pores are almost twice as large as the alamethicin pores. On the basis of the in-plane scattering data, we propose a toroidal (or wormhole) model, which differs from the barrel-stave model of alamethicin in that the lipid bends back on itself like the inside of a torus. The bending requires a lateral expansion in the head group region of the bilayer. Magainin monomers play the role of fillers in the expansion region thereby stabilizing the pore. This molecular configuration is consistent with all published magainin data.