Interaction of fatty acids with lipid bilayers

We present a method by which it is possible to describe the binding of fatty acids to phospholipid bilayers. Binding constants for oleic acid and a number of fatty acids used as spectroscopic probes are deduced from electrophoresis measurements. There is a large shift in $\text{p}\text{K}$ value for the fatty acids on binding to the phospholipid bilayers, consistent with stronger binding of the uncharged form of the fatty acid. For dansylundecanoic acid, fluorescence titrations are consistent with the binding constants derived from the electrophoresis experiments. For 12-(9-anthroyloxy)stearic acid, fluorescence and electrophoresis data are inconsistent, and we attribute this to quenching of fluorescence at high molar ratios of 12-anthroylstearic acid to phospholipid in the bilayer.     

Interaction of small peptides with lipid bilayers.

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.     

Molecular dynamics simulations of indolicidin association with model lipid bilayers

Identifying the mechanisms responsible for the interaction of peptides with cell membranes is critical to the design of new antimicrobial peptides and membrane transporters. We report here the results of a computational simulation of the interaction of the 13-residue peptide indolicidin with single-phase lipid bilayers of dioleoylphosphatidylcholine, distearoylphosphatidylcholine, dioleoylphosphatidylglycerol, and distearoylphosphatidylglycerol. Ensemble analysis of the membrane-bound peptide revealed that, in contrast to the extended, linear backbone structure reported for indolicidin in sodium dodecyl sulphate detergent micelles, the peptide adopts a boat-shaped conformation in both phosphatidylglycerol and phosphatidylcholine lipid bilayers, similar to that reported for dodecylphosphocholine micelles. In agreement with fluorescence and NMR experiments, simulations confirmed that the peptide localizes in the membrane interface, with the distance between phosphate headgroups of each leaflet being reduced in the presence of indolicidin. These data, along with a concomitant decrease in lipid order parameters for the upper-tail region, suggest that indolicidin binding results in membrane thinning, consistent with recent in situ atomic force microscopy studies.     

Interaction of nonionic PEO-PPO diblock copolymers with lipid bilayers

The relationship between the molecular architecture of a series of poly(ethylene oxide)-b-poly(propylene oxide) (PEO-PPO) diblock copolymers and the nature of their interactions with lipid bilayers has been studied using small- and wide-angle X-ray scattering (SAXS and WAXS) and differential scanning calorimetry (DSC). The number of molecular repeat units in the hydrophobic PPO block has been found to be a critical determinant of the nature of diblock copolymer-lipid bilayer association. For dimyristoyl-sn-glycero-3-phosphocholine (DMPC)-based biomembrane structures, polymers whose PPO chain length approximates that of the acyl chains of the lipid bilayer yield highly ordered, expanded lamellar structures consistent with well-integrated (into the lipid bilayer) PPO blocks. Shorter diblock copolymers produce mixed lamellar and nonlamellar mesophases. The thermotropic phase behavior of the polymer-doped membrane systems is highly influenced by the presence and molecular architecture of the diblock copolymer, as evidenced by shifting of the main phase transition to higher temperatures, broadening of the main transition, and the appearance of other features. Studies of temperature-induced changes in the mesophase structure for compositions prepared with well-integrated PEO-PPO polymers indicate that they undergo reversible changes to a nonlamellar structure as the temperature is lowered. Increasing either the number of repeat units in the PEO block or the polymer concentration promotes a greater degree of structural ordering.     

Interaction of synaptotagmin with lipid bilayers, analyzed by single-molecule force spectroscopy

Synaptotagmin I is the major Ca²⁺ sensor for membrane fusion during neurotransmitter release. The cytoplasmic domain of synaptotagmin consists of two C2 domains, C2A and C2B. On binding Ca²⁺, the tips of the two C2 domains rapidly and synchronously penetrate lipid bilayers. We investigated the forces of interaction between synaptotagmin and lipid bilayers using single-molecule force spectroscopy. Glutathione-S-transferase-tagged proteins were attached to an atomic force microscope cantilever via a glutathione-derivatized polyethylene glycol linker. With wild-type C2AB, the force profile for a bilayer containing phosphatidylserine had both Ca²⁺-dependent and Ca²⁺-independent components. No force was detected when the bilayer lacked phosphatidylserine, even in the presence of Ca²⁺. The binding characteristics of C2A and C2B indicated that the two C2 domains cooperate in binding synaptotagmin to the bilayer, and that the relatively weak Ca²⁺-independent force depends only on C2A. When the lysine residues K189-192 and K326, 327 were mutated to alanine, the strong Ca²⁺-dependent binding interaction was either absent or greatly reduced. We conclude that synaptotagmin binds to the bilayer via C2A even in absence of Ca²⁺, and also that positively charged regions of both C2A and C2B are essential for the strong Ca²⁺-dependent binding of synaptotagmin to the bilayer.     

Bridging microscopic and mesoscopic simulations of lipid bilayers

A lipid bilayer is modeled using a mesoscopic model designed to bridge atomistic bilayer simulations with macro-scale continuum-level simulation. Key material properties obtained from detailed atomistic-level simulations are used to parameterize the meso-scale model. The fundamental length and time scale of the meso-scale simulation are at least an order of magnitude beyond that used at the atomistic level. Dissipative particle dynamics cast in a new membrane formulation provides the simulation methodology. A meso-scale representation of a dimyristoylphosphatidylcholine membrane is examined in the high and low surface tension regimes. At high surface tensions, the calculated modulus is found to be slightly less than the atomistically determined value. This result agrees with the theoretical prediction that high-strain thermal undulations still persist, which have the effect of reducing the value of the atomistically determined modulus. Zero surface tension simulations indicate the presence of strong thermal undulatory modes, whereas the undulation spectrum and the calculated bending modulus are in excellent agreement with theoretical predictions and experiment.     

Constant surface tension molecular dynamics simulations of lipid bilayers with trehalose

Surface areas and fluctuations evaluated from 50 ns molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers in a 1:2 trehalose:lipid ratio carried out at surface tensions 10, 17 and 25 dyn/cm/leaflet are compared with those of pure bilayers under the same conditions. Trehalose increases the surface area, as consistent with the surface tension lowering observed in simulations at constant area. The system bulk elastic modulus K _b  = 1.5 ± 0.3 × 10¹⁰ dyn/cm². It is independent of bilayer surface area and trehalose content within statistical error. In contrast, the area elastic modulus K _a shows a strong area dependence. At 64 Å²/lipid (the experimental surface area), K _a  = 138 ± 26 dyn/cm for a pure DPPC bilayer and 82 ± 10 dyn/cm for one with trehalose; i.e. trehalose increases fluidity of the bilayer surface at this area per lipid.     

Coupling field theory with mesoscopic dynamical simulations of multicomponent lipid bilayers

A method for simulating a two-component lipid bilayer membrane in the mesoscopic regime is presented. The membrane is modeled as an elastic network of bonded points; the spring constants of these bonds are parameterized by the microscopic bulk modulus estimated from earlier atomistic nonequilibrium molecular dynamics simulations for several bilayer mixtures of DMPC and cholesterol. The modulus depends on the composition of a point in the elastic membrane model. The dynamics of the composition field is governed by the Cahn-Hilliard equation where a free energy functional models the coupling between the composition and curvature fields. The strength of the bonds in the elastic network are then modulated noting local changes in the composition and using a fit to the nonequilibrium molecular dynamics simulation data. Estimates for the magnitude and sign of the coupling parameter in the free energy model are made treating the bending modulus as a function of composition. A procedure for assigning the remaining parameters in the free energy model is also outlined. It is found that the square of the mean curvature averaged over the entire simulation box is enhanced if the strength of the bonds in the elastic network are modulated in response to local changes in the composition field. We suggest that this simulation method could also be used to determine if phase coexistence affects the stress response of the membrane to uniform dilations in area. This response, measured in the mesoscopic regime, is already known to be conditioned or renormalized by thermal undulations.     

Interaction of tea catechins with lipid bilayers investigated by a quartz-crystal microbalance analysis

The quartz-crystal microbalance (QCM) technique was applied to investigate the interaction of tea catechins with lipid bilayers. The association constants obtained from the frequency changes of QCM revealed that (-)epicatechin gallate and (-)epigallocatechin gallate interacted with 1,2-dimyristoyl-sn-glycero-3-phosphocholine ca. 1000 times more strongly than (-)epicatechin and (-)epigallocatechin. The results exhibited good correlation with the strength of biological activity.     

Interaction of glycosylated human myelin basic protein with lipid bilayers

Myelin basic protein (MBP), isolated from normal human myelin, was glycosylated with UDP-N-acetyl-D-galactosamine and a glycosyltransferase isolated from porcine submaxillary glands. MBP containing 0.85 mol of N-acetyl-D-galactosamine per mole of protein was oxidized at carbon 6 by galactose oxidase and complexed with a spin-label, Tempoamine, in order to study its interactions with lipids. When the spin-labeled MBP was reacted with lipid vesicles consisting of DSPG, DPPG, and DMPG, most of the spin-label was motionally restricted in the gel phase, with a correlation time greater than 10(-8)s. The motion increased with increasing temperature and was sensitive to the lipid phase transition. Interaction with the gel phase of DPPA caused much less motional restriction of the probe. However, melting of the lipid allowed increased interaction and motional restriction of the probe, which was only partially reversed on cooling back to the gel phase. The motional restriction of the probe in these lipids is attributed to its penetration partway into the lipid bilayer in both the gel and liquid-crystalline phases. The fact that the probe bound to the protein can penetrate partway into the bilayer suggests that other hydrophobic side chains and residues of the protein can similarly penetrate into the bilayer. Additional evidence for penetration was provided by digestion of the lipid-bound protein with endoproteinase Lys-C. When nonglycosylated and glycosylated MBP in solution was treated with Lys-C, extensive digestion occurred. A single radioactive peptide which eluted at 25 min was identified as residues 92-105.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interaction of tea catechins with lipid bilayers investigated with liposome systems

Interaction of tea catechins with lipid bilayers was investigated with liposome systems, which enabled us to separate liposomes from the external medium by centrifugation. We found that epicatechin gallate had the highest affinity for lipid bilayers, followed by epigallocatechin gallate, epicatechin, and epigallocatechin. Epicatechin gallate and epigallocatechin gallate in the surface of lipid bilayer perturbed the membrane structure.     

Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

Determination of component volumes of lipid bilayers from simulations.

An efficient method for extracting volumetric data from simulations is developed. The method is illustrated using a recent atomic-level molecular dynamics simulation of L alpha phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer. Results from this simulation are obtained for the volumes of water (VW), lipid (V1), chain methylenes (V2), chain terminal methyls (V3), and lipid headgroups (VH), including separate volumes for carboxyl (Vcoo), glyceryl (Vgl), phosphoryl (VPO4), and choline (Vchol) groups. The method assumes only that each group has the same average volume regardless of its location in the bilayer, and this assumption is then tested with the current simulation. The volumes obtained agree well with the values VW and VL that have been obtained directly from experiment, as well as with the volumes VH, V2, and V3 that require certain assumptions in addition to the experimental data. This method should help to support and refine some assumptions that are necessary when interpreting experimental data.     

Interaction of protegrin-1 with lipid bilayers: membrane thinning effect

Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their beta-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The beta-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged beta-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.     

Interaction of beta-lactoglobulin with phospholipid bilayers: a molecular level elucidation as revealed by infrared spectroscopy

Fourier transform infrared (FTIR) spectroscopy has been used to study, at a molecular level, the interactions between beta-lactoglobulin (BLG), the most abundant globular protein in milk, and some lipids (sphingomyelin, SM; dimyristoylphosphatidylcholine, DMPC; dipalmytoylphosphatidylcholine, DPPC; dimyristoylphosphatidylserine-sodium salt, DMPS; dipalmitoylphosphatidylserine-sodium salt, DPPS) constituting the milk fat globule membrane (MFGM). The interactions were monitored with respect to alteration in the secondary structure of BLG, as registered by the amide I' band, and phospholipid conformation, as revealed by the acyl chain and carbonyl bands. The results show that neither the conformation nor the thermotropism of neutral bilayers containing DMPC or DPPC is affected by BLG. Reciprocally, the secondary structure and thermal behaviour of pure BLG remain the same in the presence of PC. These results suggest that no interaction occurs between PC and BLG, in agreement with previous studies. However, it is found that BLG interacts with neutral bilayers constituted by milk SM lipids, increasing gauche conformers and thus conformational disorder of the lipid acyl chains. This perturbing effect has been attributed to a partial penetration of BLG into the hydrophobic core of the bilayer, which allows hydrophobic interactions between BLG and SM. Moreover, the fact that SM possesses the same headgroup of PC implies that the head group does not prevent the occurrence of BLG-lipid interactions and other lipid regions can control the binding of BLG to lipids. Furthermore, BLG was found to interact electrostatically with charged bilayers containing PS, leading to a rigidification of the lipid hydrocarbon chains and a dehydration of the interfacial region. This last effect suggests that the protein limits the accessibility of water molecules to the interfacial region of the phospholipids by its presence at the membrane surface.     

Interpretation of small angle X-ray measurements guided by molecular dynamics simulations of lipid bilayers

Reconstruction and interpretation of lipid bilayer structure from X-ray scattering often rely on assumptions regarding the molecular distributions across the bilayer. It is usually assumed that changes in head-head spacings across the bilayer, as measured from electron density profiles, equal the variations in hydrocarbon thicknesses. One can then determine the structure of a bilayer by comparison to the known structure of a lipid with the same headgroup. Here we examine this procedure using simulated electron density profiles for the benchmark lipids DMPC and DPPC. We compare simulation and experiment in both real and Fourier space to address two main aspects: (i) the measurement of head-head spacings from relative electron density profiles, and (ii) the determination of the absolute scale for these profiles. We find supporting evidence for the experimental procedure, thus explaining the robustness and consistency of experimental structural results derived from electron density profiles. However, we also expose potential pitfalls in the Fourier reconstruction that are due to the limited number of scattering peaks. Volumetric analysis of simulated bilayers allows us to propose an improved, yet simple method for scale determination. In this way we are able to remove some of the restrictions imposed by limited scattering data in constructing reliable electron density profiles.     

Interaction of fluorescently labeled pardaxin and its analogues with lipid bilayers.

Fluorescence measurements were used to monitor the interaction of the neurotoxin pardaxin and its analogues with membranes. Eight peptides were selectively labeled with the fluorophore 7-nitrobenz-2-oxa-1,3-diazole-4-yl, either at their N-terminal or at their C-terminal. No detectable changes in membrane permeability or hemolytic activity were observed upon modification. Upon the titration of solutions containing the different peptides with small unilamellar vesicles, the fluorescent emission spectra of 7-nitrobenz-2-oxa-1,3-diazole-4-yl-labeled pardaxin and its analogues, but not those of control peptides, displayed blue shifts in addition to enhanced intensities upon relocation of the probe to a more apolar environment. The results revealed that the N terminus of pardaxin is buried within the lipid bilayer while the C terminus is located at the bilayer's surface. Binding isotherms were obtained from the observed increases in the fluorescence emission yields, from which surface partition constants, in the range of 10(4) M-1, were in turn derived. The existence of an aggregation process was suggested by the shape of the binding isotherms. Furthermore, the results show good correlation between the incidence of aggregation and the ability of the different analogues to induce the release of relatively large molecules from vesicles. As such, our results suggest that the mechanism of pore formation employed by pardaxin and its analogues could be described by the "barrel stave" model.     

Interaction of sporidesmin,a mycotoxin fromPithomyces chartarum,with lipid bilayers

Sporidesmin, a mycotoxin fromPithomyces chartarum is a hydrophobic molecule. It can therefore be easily incorporated in the cell membrane, where it is likely to cause changes in the bilayer organization and the properties of membrane proteins. In order to understand the redox behaviour of sporidesmin in a hydrophobic environment, we have investigated the effects of oxidized and reduced sporidesmin on the phase transition properties of bilayers and on the susceptibility of bilayers to pancreatic phospholipase A₂ (PLA₂). The changes induced by sporidesmin in the thermotropic phase transition profiles of dimyristoyl-sn-3-phosphatidyl choline (DMPC) bilayers were similar to those caused by solutes known to localize in the glycerol-backbone region of the lipid bilayer, suggesting a similar localization for oxidized and reduced sporidesmin. Neither form of toxin disrupt the bilayer or membrane organization even at relatively high mole fractions. At concentrations <10 mole% both forms partitioned equally well in the gel and liquid-crystalline phases, whereas at higher concentrations (30 mole%) reduced sporidesmin is preferentially localized in the liquid-crystalline phase. These effects of sporidesmin on the phase properties of DMPC vesicles were also reported by the fluorescence behavior of 10-pyrenedecanoic acid (PDA). The effects of oxidized and reduced sporidesmins on PLA₂ kinetics are consistent with their ability to perturb bilayer organisation.     

Interaction of a two-transmembrane-helix peptide with lipid bilayers and dodecyl sulfate micelles

To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, F?rster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix-helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods.     

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.