A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Ku?erka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (D_B), hydrophobic thickness (D_C), and lipid volume (V_L). The only minor discrepancy between simulation and experiment is a measure of (D_B − D_HH) / 2 where D_HH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.     

Structure and dynamics of the two amphipathic arginine-rich peptides RW9 and RL9 in a lipid environment investigated by solid-state NMR and MD simulations

Cell penetrating peptides (CPPs) are able to cross membranes without using receptors but only little information about the underlying mechanism is available. In this work, we investigate the interaction of the two arginine-rich CPPs RW9 and RL9 with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), and POPC/POPG membranes with varying POPG content using isothermal titration calorimetry (ITC), solid-state nuclear magnetic resonance (NMR) spectroscopy, and molecular dynamics (MD) simulations. Both peptides were derived from the known CPP penetratin and it was shown previously that RW9 is able to penetrate membranes better than RL9. Overall, the results show that both RW9 and RL9 have a relatively small influence on the membrane. They increase the order of the lipids in the headgroup region and reduce order in the acyl chains indicating that they are located in the lipid/water interface. In addition, the flexibility of the membrane is slightly increased by both peptides but RW9 has a larger influence than RL9. The differences observed in the influences on POPC and POPG as well as MD simulations on the mixed POPC/POPG bilayers of 850 ns length each show that both peptides preferentially associate with and enrich the charged PG lipids almost 2fold in an area of 12 Å around the peptides. As expected, we could not observe any membrane crossing on the simulation time scale of 850 ns but observed that some peptides flipped their orientation during binding to the membrane. Interestingly, all observed flips coincided with structural changes in the peptides indicating that structural changes or flexibility might play a role during the binding of arginine-rich CPPs to membranes.     

Perturbation of DPPC/POPG bilayers by the N-terminal helix of lung surfactant protein SP-B: a <Superscript>2</Superscript>H NMR study

SP-B_8–25 is a synthetic peptide comprising the N-terminal helix of the essential lung surfactant protein SP-B. Rat lung oxygenation studies have shown that SP-B_8–25 retains some of the function of full-length SP-B. We have used deuterium nuclear magnetic resonance (²H-NMR) to examine the influence of SP-B_8–25 on the mixing properties of saturated PC and unsaturated PG lipids in model mixed lipid bilayers containing dipalmitoylphosphatidylcholine (DPPC) and palmitoyl-oleoyl-phosphatidylglycerol (POPG), in a molar ratio of 7:3. In the absence of the peptide, ²H-NMR spectra of DPPC/POPG mixtures, with one or the other lipid component deuterated, indicate coexistence of large liquid crystal and gel domains over a range of about 10°C through the liquid crystal to gel transition of the bilayer. Addition of SP-B_8–25 has little effect on the width of the transition but the spectra through the transition range cannot be resolved into distinct liquid crystal and gel spectral components suggesting that the peptide interferes with the tendency of the DPPC and POPG lipid components in this mixture to phase separate near the bilayer transition temperature. Quadrupole echo decay observations suggest that the peptide may also reduce differences in the correlation times for local reorientation of the two lipids. These observations suggest that SP-B_8–25 promotes a more thorough mixing of saturated PC and unsaturated PG components and may be relevant to understanding the behaviour of lung surfactant material under conditions of lateral compression which might be expected to enhance the propensity for saturated and unsaturated surfactant lipid components to segregate.     

Exploring membrane selectivity of the antimicrobial peptide KIGAKI using solid-state NMR spectroscopy

The designed antimicrobial peptide KIGAKIKIGAKIKIGAKI possesses enhanced membrane selectivity for bacterial lipids, such as phosphatidylethanolamine and phosphatidylglycerol. The perturbation of the bilayer by the peptide was first monitored using oriented bilayer samples on glass plates. The alignment of POPE/POPG model membranes with respect to the bilayer normal was severely altered at 4 mol% KIGAKI while the alignment of POPC bilayers was retained. The interaction mechanism between the peptide and POPE/POPG bilayers was investigated by carefully comparing three bilayer MLV samples (POPE bilayers, POPG bilayers, and POPE/POPG 4/1 bilayers). KIGAKI induces the formation of an isotropic phase for POPE/POPG bilayers, but only a slight change in the ³¹P NMR CSA line shape for both POPE and POPG bilayers, indicating the synergistic roles of POPE and POPG lipids in the disruption of the membrane structure by KIGAKI. ²H NMR powder spectra show no reduction of the lipid chain order for both POPG and POPE/POPG bilayers upon peptide incorporation, supporting the evidence that the peptide acts as a surface peptide. ³¹P longitudinal relaxation studies confirmed that different dynamic changes occurred upon interaction of the peptide with the three different lipid bilayers, indicating that the strong electrostatic interaction between the cationic peptide KIGAKI and anionic POPG lipids is not the only factor in determining the antimicrobial activity. Furthermore, ³¹P and ²H NMR powder spectra demonstrated a change in membrane characteristics upon mixing of POPE and POPG lipids. The interaction between different lipids, such as POPE and POPG, in the mixed bilayers may provide the molecular basis for the KIGAKI carpet mechanism in the permeation of the membrane.     

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na⁺ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the ³¹P chemical shift anisotropy and ¹⁴N and CD₃ NMR quadrupole splittings.     

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.     

Molecular Dynamics Simulations of the Bacterial UraA H+-Uracil Symporter in Lipid Bilayers Reveal a Closed State and a Selective Interaction with Cardiolipin

The Escherichia coli UraA H⁺-uracil symporter is a member of the nucleobase/ascorbate transporter (NAT) family of proteins, and is responsible for the proton-driven uptake of uracil. Multiscale molecular dynamics simulations of the UraA symporter in phospholipid bilayers consisting of: 1) 1-palmitoyl 2-oleoyl-phosphatidylcholine (POPC); 2) 1-palmitoyl 2-oleoyl-phosphatidylethanolamine (POPE); and 3) a mixture of 75% POPE, 20% 1-palmitoyl 2-oleoyl-phosphatidylglycerol (POPG); and 5% 1-palmitoyl 2-oleoyl-diphosphatidylglycerol/cardiolipin (CL) to mimic the lipid composition of the bacterial inner membrane, were performed using the MARTINI coarse-grained force field to self-assemble lipids around the crystal structure of this membrane transport protein, followed by atomistic simulations. The overall fold of the protein in lipid bilayers remained similar to the crystal structure in detergent on the timescale of our simulations. Simulations were performed in the absence of uracil, and resulted in a closed state of the transporter, due to relative movement of the gate and core domains. Anionic lipids, including POPG and especially CL, were found to associate with UraA, involving interactions between specific basic residues in loop regions and phosphate oxygens of the CL head group. In particular, three CL binding sites were identified on UraA: two in the inner leaflet and a single site in the outer leaflet. Mutation of basic residues in the binding sites resulted in the loss of CL binding in the simulations. CL may play a role as a “proton trap” that channels protons to and from this transporter within CL-enriched areas of the inner bacterial membrane.     

Binding,folding and insertion of a β-hairpin peptide at a lipid bilayer surface: Influence of electrostatics and lipid tail packing

Antimicrobial peptides (AMPs) act as host defenses against microbial pathogens. Here we investigate the interactions of SVS-1 (KVKVKVKV^dP^lPTKVKVKVK), an engineered AMP and anti-cancer β-hairpin peptide, with lipid bilayers using spectroscopic studies and atomistic molecular dynamics simulations. In agreement with literature reports, simulation and experiment show preferential binding of SVS-1 peptides to anionic over neutral bilayers. Fluorescence and circular dichroism studies of a Trp-substituted SVS-1 analogue indicate, however, that it will bind to a zwitterionic DPPC bilayer under high-curvature conditions and folds into a hairpin. In bilayers formed from a 1:1 mixture of DPPC and anionic DPPG lipids, curvature and lipid fluidity are also observed to promote deeper insertion of the fluorescent peptide. Simulations using the CHARMM C36m force field offer complementary insight into timescales and mechanisms of folding and insertion. SVS-1 simulated at an anionic mixed POPC/POPG bilayer folded into a hairpin over a microsecond, the final stage in folding coinciding with the establishment of contact between the peptide's valine sidechains and the lipid tails through a “flip and dip” mechanism. Partial, transient folding and superficial bilayer contact are seen in simulation of the peptide at a zwitterionic POPC bilayer. Only when external surface tension is applied does the peptide establish lasting contact with the POPC bilayer. Our findings reveal the influence of disruption to lipid headgroup packing (via curvature or surface tension) on the pathway of binding and insertion, highlighting the collaborative effort of electrostatic and hydrophobic interactions on interaction of SVS-1 with lipid bilayers.     

NMR investigations of interactions between anesthetics and lipid bilayers

Interactions between anesthetics (lidocaine and short chain alcohols) and lipid membranes formed by dimyristoylphosphatidylcholine (DMPC) were studied using NMR spectroscopy. The orientational order of lidocaine was investigated using deuterium NMR on a selectively labelled compound whereas segmental ordering in the lipids was probed by two-dimensional ¹H-¹³C separated local field experiments under magic-angle spinning conditions. In addition, trajectories generated in molecular dynamics (MD) computer simulations were used for interpretation of the experimental results. Separate simulations were carried out with charged and uncharged lidocaine molecules. Reasonable agreement between experimental dipolar interactions and the calculated counterparts was observed. Our results clearly show that charged lidocaine affects significantly the lipid headgroup. In particular the ordering of the lipids is increased accompanied by drastic changes in the orientation of the P-N vector in the choline group.     

NaCl interactions with phosphatidylcholine bilayers do not alter membrane structure but induce long-range ordering of ions and water

It is generally accepted that ions interact directly with lipids in biological membranes. Decades of biophysical studies on pure lipid bilayer systems have shown that only certain types of ions, most significantly large anions and multivalent cations, can fundamentally alter the structure and dynamics of lipid bilayers. It has long been accepted that at physiological concentrations NaCl ions do not alter the physical behavior or structure of bilayers composed solely of zwitterionic phosphatidylcholine (PC) lipids. Recent X-ray scattering experiments have reaffirmed this dogma, showing that below 1 M concentration, NaCl does not significantly alter bilayer structure. However, despite this history, there is an ongoing controversy within the molecular dynamics (MD) simulation community regarding NaCl/PC interactions. In particular, the CHARMM and GROMOS force fields show dramatically different behavior, including the effect on bilayer structure, surface potential, and the ability to form stable, coordinated ion–lipid complexes. Here, using long-timescale, constant-pressure simulations under the newest version of the CHARMM force field, we find that Na⁺ and Cl⁻ associate with PC head groups in a POPC bilayer with approximately equal, though weak, affinity, and that the salt has a negligible effect on bilayer structure, consistent with earlier CHARMM results and more recent X-ray data. The results suggest that interpretation of simulations where lipids interact with charged groups of any sort, including charged proteins, must be carefully scrutinized.     

Extracting Curvature Preferences of Lipids Assembled in Flat Bilayers Shows Possible Kinetic Windows for Genesis of Bilayer Asymmetry and Domain Formation in Biological Membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures—shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35–40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.     

Effect of Membrane Composition on Antimicrobial Peptides Aurein 2.2 and 2.3 From Australian Southern Bell Frogs

The effects of hydrophobic thickness and the molar phosphatidylglycerol (PG) content of lipid bilayers on the structure and membrane interaction of three cationic antimicrobial peptides were examined: aurein 2.2, aurein 2.3 (almost identical to aurein 2.2, except for a point mutation at residue 13), and a carboxy C-terminal analog of aurein 2.3. Circular dichroism results indicated that all three peptides adopt an α-helical structure in the presence of a 3:1 molar mixture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (DMPC/DMPG), and 1:1 and 3:1 molar mixtures of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (POPC/POPG). Oriented circular dichroism data for three different lipid compositions showed that all three peptides were surface-adsorbed at low peptide concentrations, but were inserted into the membrane at higher peptide concentrations. The ³¹P solid-state NMR data of the three peptides in the DMPC/DMPG and POPC/POPG bilayers showed that all three peptides significantly perturbed lipid headgroups, in a peptide or lipid composition-dependent manner. Differential scanning calorimetry results demonstrated that both amidated aurein peptides perturbed the overall phase structure of DMPC/DMPG bilayers, but perturbed the POPC/POPG chains less. The nature of the perturbation of DMPC/DMPG bilayers was most likely micellization, and for the POPC/POPG bilayers, distorted toroidal pores or localized membrane aggregate formation. Calcein release assay results showed that aurein peptide-induced membrane leakage was more severe in DMPC/DMPG liposomes than in POPC/POPG liposomes, and that aurein 2.2 induced higher calcein release than aurein 2.3 and aurein 2.3-COOH from 1:1 and 3:1 POPC/POPG liposomes. Finally, DiSC₃5 assay data further delineated aurein 2.2 from the others by showing that it perturbed the lipid membranes of intact S. aureus C622 most efficiently, whereas aurein 2.3 had the same efficiency as gramicidin S, and aurein 2.3-COOH was the least efficient. Taken together, these data show that the membrane interactions of aurein peptides are affected by the hydrophobic thickness of the lipid bilayers and the PG content.     

Phase-transition properties of glycerol-monopalmitate lipid bilayers investigated by molecular dynamics simulation: influence of the system size and force-field parameters

Phase-transition properties of glycerol-1-monopalmitate (GMP) bilayers are investigated using explicit-solvent molecular dynamics (MD) simulations, initiated from structures appropriate for the gel (GL) or liquid crystal (LC) phases, and carried out at different hydration levels and temperatures. Building up on a previous study and based on 600 ns simulations, the influence of the system size and of the force field on the equilibrium thermodynamic and dynamic parameters of the bilayers in the GL and LC phases, as well as on the temperature T_m and properties of the GL ? LC phase transition, are analysed. Qualitatively speaking, the results agree with the available experimental data for the area per lipid in the two phases and for the phase-transition temperatures at the three hydration levels irrespective of the selected model parameters. They also suggest that the total number of hydrogen bonds formed between a lipid headgroup and its environment is essentially constant, amounting to about four in both the LC and the GL phases. Quantitatively speaking, the dependence of T_m on the hydration level is found to be non-systematic across the different combinations of model parameters. This results in part from a sensitivity of the results on the system size and force-field parameters but also from the limited accuracy of the bracketing approach employed here to estimate T_m. Finally, a simple kinetic model is proposed to account for the timescales of the transitions. This model involves enthalpy and entropy increases of about 26 kJ mol^? 1 and 83 J mol^? 1 K^? 1 per lipid, upon going from the GL to the LC phase. The transition state is associated with activation parameters corresponding to 13% and 11%, respectively, of these values along the GL → LC transition, resulting in an activation free energy of about 0.3 kJ mol^? 1 per lipid at T_m.     

Comparative molecular dynamics simulations of the antimicrobial peptide CM15 in model lipid bilayers

We report altogether 3-μs molecular dynamics (MD) simulations of the antimicrobial peptide CM15 to systematically investigate its interaction with two model lipid bilayers, pure POPC and mixed POPG:POPC (1:2). Starting with either an α-helical or a random-coil conformation, CM15 is found to insert into both bilayers. Peptide-lipid interaction is stronger with the anionic POPG:POPC than the zwitterionic POPC, which is largely attributed to the electrostatic attraction between CM15 and the negatively charged POPG. Simulations initiated with CM15 as a random coil allowed us to study peptide folding at the lipid-water interface. Interestingly, CM15 folding appears to be faster in POPC than POPG:POPC, which may be explained by a lower activation energy barrier of structural rearrangement in the former system. Our data also suggest that compared with the random-coil conformation, CM15 in a pre-folded α-helix has significantly reduced interactions with the lipids, indicating that peptide initial structures may bias the simulation results considerably on the 100-ns timescale. The implications of this result should be considered when preparing and interpreting future AMP simulations.     

Comparative study of stability and transport of molecules through cyclic peptide nanotube and aquaporin: a molecular dynamics simulation approach

Abstract

The structural stability and transport properties of the cyclic peptide nanotube (CPN) 8?×?[Cys–Gly–Met–Gly]₂ in different phospholipid bilayers such as POPA (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) and POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine) with water have been investigated using molecular dynamics (MD) simulation. The hydrogen bonds and non-bonded interaction energies were calculated to study the stability in different bilayers. One µs MD simulation in POPA lipid membrane reveals the stability of the cyclic peptide nanotube, and the simulations at various temperatures manifest the higher stability of 8?×?[Cys–Gly–Met–Gly]₂. We demonstrated that the presence of sulphur-containing amino acids in CPN enhances the stability through disulphide bonds between the adjacent rings. Further, the water permeation coefficient of the CPN is calculated and compared with human aquaporin-2 (AQP2) channel protein. It is found that the coefficients are highly comparable to the AQP2 channel though the mechanism of water transport is not similar to AQP 2; the flow of water in the CPN is taking place as a two-line 1–2–1–2 file fashion. In addition to that, the transport behavior of Na⁺ and K⁺ ions, single water molecule, urea and anti-cancer drug fluorouracil were investigated using pulling simulation and potential of mean force calculation. The above transport behavior shows that Na⁺ is trapped in CPN for a longer time than other molecules. Also, the interactions of the ions and molecules in Cα and mid-Cα plane were studied to understand the transport behavior of the CPN. Abbreviations AQP2 Aquaporin-2

CPN Cyclic peptide nanotube

MD Molecular dynamics

POPA 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid

POPE 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine

POPG 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol

POPS 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine

Communicated by Ramaswamy H. Sarma     

Exploring membrane selectivity of the antimicrobial peptide KIGAKI using solid-state NMR spectroscopy

The designed antimicrobial peptide KIGAKIKIGAKIKIGAKI possesses enhanced membrane selectivity for bacterial lipids, such as phosphatidylethanolamine and phosphatidylglycerol. The perturbation of the bilayer by the peptide was first monitored using oriented bilayer samples on glass plates. The alignment of POPE/POPG model membranes with respect to the bilayer normal was severely altered at 4 mol% KIGAKI while the alignment of POPC bilayers was retained. The interaction mechanism between the peptide and POPE/POPG bilayers was investigated by carefully comparing three bilayer MLV samples (POPE bilayers, POPG bilayers, and POPE/POPG 4/1 bilayers). KIGAKI induces the formation of an isotropic phase for POPE/POPG bilayers, but only a slight change in the (31)P NMR CSA line shape for both POPE and POPG bilayers, indicating the synergistic roles of POPE and POPG lipids in the disruption of the membrane structure by KIGAKI. (2)H NMR powder spectra show no reduction of the lipid chain order for both POPG and POPE/POPG bilayers upon peptide incorporation, supporting the evidence that the peptide acts as a surface peptide. (31)P longitudinal relaxation studies confirmed that different dynamic changes occurred upon interaction of the peptide with the three different lipid bilayers, indicating that the strong electrostatic interaction between the cationic peptide KIGAKI and anionic POPG lipids is not the only factor in determining the antimicrobial activity. Furthermore, (31)P and (2)H NMR powder spectra demonstrated a change in membrane characteristics upon mixing of POPE and POPG lipids. The interaction between different lipids, such as POPE and POPG, in the mixed bilayers may provide the molecular basis for the KIGAKI carpet mechanism in the permeation of the membrane.     

Evidence of pores and thinned lipid bilayers induced in oriented lipid membranes interacting with the antimicrobial peptides, magainin-2 and aurein-3.3

Dynamic structures of supramolecular lipid assemblies, such as toroidal pores and thinned bilayers induced in oriented lipid membranes, which are interacting with membrane-acting antimicrobial peptides (AMPs), magainin-2 and aurein-3.3, were explored by ³¹P and ²H solid-state NMR (ssNMR) spectroscopy. Various types of phospholipid systems, such as POPC-d₃₁, POPC-d₃₁/POPG, and POPC-d₃₁/cholesterol, were investigated to understand the membrane disruption mechanisms of magainin-2 and aurein-3.3 peptides at various peptide-to-lipid (P:L) ratios. The experimental lineshapes of anisotropic ³¹P and ²H ssNMR spectra measured on these peptide-lipid systems were simulated reasonably well by assuming the presence of supramolecular lipid assemblies, such as toroidal pores and thinned bilayers, in membranes. Furthermore, the observed decrease in the anisotropic frequency span of either ³¹P or ²H ssNMR spectra of oriented lipid bilayers, particularly when anionic POPG lipids are interacting with AMPs at high P:L ratios, can directly be explained by a thinned membrane surface model with fast lateral diffusive motions of lipids. The spectral analysis protocol we developed enables extraction of the lateral diffusion coefficients of lipids distributed on the curved surfaces of pores and thinned bilayers on a few nanometers scale.     

Membrane models of E. coli containing cyclic moieties in the aliphatic lipid chain

Nearly all molecular dynamics simulations of bacterial membranes simplify the lipid bilayer by composing it of only one or two lipids. Previous attempts of developing a model E. coli membrane have used only 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) POPG lipids. However, an important constituent of bacterial membranes are lipids containing a cyclopropane ring within the acyl chain. We have developed a complex membrane that more accurately reflects the diverse population of lipids within E. coli cytoplasmic membranes, including lipids with a cyclic moiety. Differences between the deuterium order profile of cyclic lipids and monounsaturated lipids are observed. Furthermore, the inclusion of the cyclopropane ring decreases the surface density of the bilayer and produces a more rigid membrane as compared to POPE/POPG membranes. Additionally, the diverse acyl chain length creates a thinner bilayer which matches the hydrophobic thickness of E. coli transmembrane proteins better than the POPE/POPG bilayer. We believe that the complex lipid bilayer more accurately describes a bacterial membrane and suggest the use of it in molecular dynamic simulations rather than simple POPE/POPG membranes.     

Structure of molecular aggregates of 1-(3-sn-phosphatidyl)-l-myo-inositol 3,4-bis(phosphate) in water

The molecular organization of 1-(3-sn-phosphatidyl)-l-myo-inositol 3,4-bis-(phosphate)/water systems is investigated over a wide range of lipid concentrations using X-ray diffraction, calorimetry, analytical ultracentrifugation, densitometry and viscometry.At high lipid concentrations, the lipid molecules are found to form a lamellar phase. The repeat distance increases from 60 to 120 Å with increasing water content to 70 wt% and the surface area per lipid molecule increases from 41.7 Ų to a limiting value of 100 Ų.On the other hand, at very low lipid concentrations the molecules are found to form not vesicles but micelles, the total molecular weight of which takes a value of 93 000.This finding revises the prevalent view that lipids containing two (or more) hydrocarbon chains form extended bilayers or vesicles, whereas single chained lipids form micelles. (Tanford, C.(1972) J. Phys. Chem. 76, 3020–3024).