Atomistic Simulations of Bicelle Mixtures

Mixtures of long- and short-tail phosphatidylcholine lipids are known to self-assemble into a variety of aggregates combining flat bilayerlike and curved micellelike features, commonly called bicelles. Atomistic simulations of bilayer ribbons and perforated bilayers containing dimyristoylphosphatidylcholine (DMPC, di-C₁₄ tails) and dihexanoylphosphatidylcholine (DHPC, di-C₆ tails) have been carried out to investigate the partitioning of these components between flat and curved microenvironments and the stabilization of the bilayer edge by DHPC. To approach equilibrium partitioning of lipids on an achievable simulation timescale, configuration-bias Monte Carlo mutation moves were used to allow individual lipids to change tail length within a semigrand-canonical ensemble. Since acceptance probabilities for direct transitions between DMPC and DHPC were negligible, a third component with intermediate tail length (didecanoylphosphatidylcholine, di-C₁₀ tails) was included at a low concentration to serve as an intermediate for transitions between DMPC and DHPC. Strong enrichment of DHPC is seen at ribbon and pore edges, with an excess linear density of ∼3 nm⁻¹. The simulation model yields estimates for the onset of edge stability with increasing bilayer DHPC content between 5% and 15% DHPC at 300 K and between 7% and 17% DHPC at 323 K, higher than experimental estimates. Local structure and composition at points of close contact between pores suggest a possible mechanism for effective attractions between pores, providing a rationalization for the tendency of bicelle mixtures to aggregate into perforated vesicles and perforated sheets.     

Mechanism of structural transformations induced by antimicrobial peptides in lipid membranes

It has long been suggested that pore formation is responsible for the increase in membrane permeability by antimicrobial peptides (AMPs). To better understand the mechanism of AMP activity, the disruption of model membrane by protegrin-1 (PG-1), a cationic antimicrobial peptide, was studied using atomic force microscopy. We present here the direct visualization of the full range of structural transformations in supported lipid bilayer patches induced by PG-1 on zwitterionic 1,2-dimyristoyl-snglycero-phospho-choline (DMPC) membranes. When PG-1 is added to DMPC, the peptide first induces edge instability at low concentrations, then pore-like surface defects at intermediate concentrations, and finally wormlike structures with a specific length scale at high concentrations. The formation of these structures can be understood using a mesophase framework of a binary mixture of lipids and peptides, where PG-1 acts as a line-active agent. Atomistic molecular dynamics simulations on lipid bilayer ribbons with PG-1 molecules placed at the edge or interior positions are carried out to calculate the effect of PG-1 in reducing line tension. Further investigation of the placement of PG-1 and its association with defects in the bilayer is carried out using unbiased assembly of a PG-1 containing bilayer from a random mixture of PG-1, DMPC, and water. A generalized model of AMP induced structural transformations is also presented in this work. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Molecular dynamics simulation of lipid reorientation at bilayer edges

Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps—a common process in bilayer patterning—and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.     

A comparison of DMPC- and DLPE-based lipid bilayers.

A 250 ps molecular dynamics simulation of the dimyristoylphosphatidylcholine (DMPC)-based lipid bilayer, including explicit water molecules, is reported. The solvent environment of the head groups and other structural properties of the bilayer have been analyzed and compared with experimental results as well as our previous simulation of the dilauroylphosphatidylethanolamine (DLPE)-based bilayer. From this comparison we find that the solvent structure around the DMPC head group (clathrate shell) is significantly different than that around the DLPE head group (typical hydrogen bonding interactions). We have modeled the probable relationship between the different solvent environments around the R-N(CH3)3+ (DMPC) and R-NH3+ (DLPE) head groups and the different interlammelar distances in these systems by performing potential of mean force (PMF) simulations on two N(CH3)4+ and NH4+ ions in water. From the PMF simulations it appears that the differences in the hydration of the DMPC and DLPE head groups is not responsible for the differences in the hydration force observed for these systems. We also find that the orientational polarization of DLPE and DMPC is similar, which suggests that solvent polarization is not responsible for the differences in the hydration repulsion behavior observed in these systems. We also examined the order parameters for DMPC and found them to be in reasonable agreement with experiment. Given the different characteristics of the DLPE and DMPC head groups, we suggest an explanation of the differences in the interlammellar spacings of bilayers composed of these like-charged lipids. From our DLPE simulations we find that the R-NH3+ head groups can interact with the nonesterified oxygens of the phosphate group in an intraleaflet or an interleaflet manner. For the latter a "cross link" between two leaflets can be formed, which causes a stabilization of the interlamellar spacings at fairly short distances. Moreover, due to the strong intraleaflet interaction we find that the DLPE interface is relatively "flat" (as opposed to DMPC-based bilayers), which results in a surface that has regions of positive and negative charge that reside in the same plane along the bilayer normal. Based on this we propose that the DLPE bilayer interface can correlate itself with another DLPE interface by alignment of the regions of positive (or negative) charge on one leaflet with the opposite charges on the opposing leaflet.     

Binding of Aβ peptide creates lipid density depression in DMPC bilayer

Using isobaric–isothermal replica exchange molecular dynamics and all-atom explicit water model we study the impact of Aβ monomer binding on the equilibrium properties of DMPC bilayer. We found that partial insertion of Aβ peptide into the bilayer reduces the density of lipids in the binding “footprint” and indents the bilayer thus creating a lipid density depression. Our simulations also reveal thinning of the bilayer and a decrease in the area per lipid in the proximity of Aβ. Although structural analysis of lipid hydrophobic core detects disordering in the orientations of lipid tails, it also shows surprisingly minor structural perturbations in the tail conformations. Finally, partial insertion of Aβ monomer does not enhance water permeation through the DMPC bilayer and even causes considerable dehydration of the lipid–water interface. Therefore, we conclude that Aβ monomer bound to the DMPC bilayer fails to perturb the bilayer structure in both leaflets. Limited scope of structural perturbations in the DMPC bilayer caused by Aβ monomer may constitute the molecular basis of its low cytotoxicity.     

Molecular dynamics generation of nonarbitrary membrane models reveals lipid orientational correlations

This report addresses the following problems associated with the generation of computer models of phospholipid bilayer membranes using molecular dynamics simulations: arbitrary initial structures and short equilibration periods, an Ewald-induced strong coupling of phospholipids, uncertainty regarding which value should be used for surface tension to alleviate the problem of the small size of the membrane, and simultaneous realization of both order parameters and the surface area. We generated a computer model of the liquid-crystalline L-alpha-dimyristoylphosphatidylcholine (DMPC) bilayer, starting from a configuration based on a crystal structure (rather than from an arbitrary structure). To break the crystalline structure, a 20-ps high-temperature pulse of 510 K (but not 450 or 480 K) was effective. The system finally obtained is an all-atom model, with Ewald summation to evaluate Coulombic interactions and a constant surface tension of 35 dynes/cm/water-membrane interface, equilibrated for 12 ns (over 50 ns total calculation time), which reproduces all of the experimentally observed parameters examined in this work. Furthermore, this model shows the presence of significant orientational correlations between neighboring alkyl chains and between shoulder vectors (which show the orientations of the lipids about their long axes) of neighboring DMPCs.     

Cholesterol effects on the phospholipid condensation and packing in the bilayer: a molecular simulation study

A 15-ns molecular dynamics simulation of the fully hydrated liquid-crystalline dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. The generated trajectory was analysed to investigate the mechanism of the Chol condensing effect on DMPC hydrocarbon chains and the influence of Chol on the chain packing in the membrane. Chol was found to induce stronger van der Waals interactions among the chains, whereas its interactions with the chains were weak. In the DMPC-Chol bilayer, as in the DMPC bilayer, DMPC chains were regularly packed around a chosen chain but around a Chol molecule they were not. DMPC gamma chains made closer contacts with Chol than the beta chains.     

Molecular Dynamics Simulation of Conformational Flexibility of Alamethicin Fragments in Aqueous and Membranous Environment

Abstract

We present here results on molecular dynamics (MD) simulation on two fragments of channel forming antibiotic peptide Alamethicin, containing isoamino butyric acid (Aib). Simulations are carried out in aqueous and membranous environment in a bilayer of 39 molecules of Dimyristoyl phosphatidyl choline (DMPC). The peptides Boc—;Pro-Aib-Ala-Aib- OBzl (Alam 1) and Boc-Leu-Aib-Pro-OBzl (Alam 2) were simulated from their crystallography coordinates. The bilayers were built from two different conformations (A and B) of DMPC reported in crystal data. The P-N dipoles were arranged hexagonally with surface area per lipid molecule 66.5 A°² and P-P separation across the bilayer 34 A°. They were hydrated by 28.6 and 25.5 water molecules per DMPC molecule. Simulations are done using AMBER 4.0 package in constant number volume temperature (NVT) condition for 100 pico seconds (ps) in aqueous environment and 250 ps of equilibrated bilayer. Geometric parameters of lipids as: bilayer thickness, order parameter of the chains, transfraction of chain torsional angles were monitored. We also monitored geometric parameters of the peptides as backbone torsional angles, distances amongst Ca atoms, angles between Cα atoms, movement of center of gravity (CG) along and perpendicular to bilayer normal. We find that membrane bilayer is slightly disturbed due to the presence of peptides. In case of alam 2 in water angles ψ1 and ψ3 showed larger variation in water compared to same in the bilayer. The peptide conformation is more stable in DMPC bilayer. However the peptides showed movement along and perpendicular to bilayer normal. This we believe is due to hydrophobic nature of these peptides.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Opsin stability and folding: modulation by phospholipid bicelles

Integral membrane proteins do not fare well when extracted from biological membranes and are unstable or lose activity in detergents commonly used for structure and function investigations. We show that phospholipid bicelles provide a valuable means of preserving alpha-helical membrane proteins in vitro by supplying a soluble lipid bilayer fragment. Both 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)/3-[(cholamidopropyl)dimethyl-ammonio]-1-propane sulfonate (Chaps) and DMPC/l-α-1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) bicelles dramatically increase the stability of the mammalian vision receptor rhodopsin as well as its apoprotein, opsin. Opsin is particularly unstable in detergent solution but can be directly purified into DMPC/Chaps. We show that opsin can also be directly purified in DMPC/DHPC bicelles to give correctly folded functional opsin, as shown by the ability to regenerate rhodopsin to  70% yield. These well-characterised DMPC/DHPC bicelles enable us to probe the influence of bicelle properties on opsin stability. These bicelles are thought to provide DMPC bilayer fragments with most DHPC capping the bilayer edge, giving a soluble bilayer disc. Opsin stability is shown to be modulated by the q value, the ratio of DMPC to DHPC, which reflects changes in the bicelle size and, thus, proportion of DMPC bilayer present. The observed changes in stability also correlate with loss of opsin secondary structure as determined by synchrotron far-UV circular dichroism spectroscopy; the most stable bicelle results in the least helix loss. The inclusion of Chaps rather than DHPC in the DMPC/Chaps bicelles, however, imparts the greatest stability. This suggests that it is not just the DMPC bilayer fragment in the bicelles that stabilises the protein, but that Chaps provides additional stability either through direct interaction with the protein or by altering the DMPC/Chaps bilayer properties within the bicelle. The significant stability enhancements and preservation of secondary structure reported here in bicelles are pertinent to other membrane proteins, notably G-protein-coupled receptors, which are unstable in detergent solution.     

Conformation of nifedipine in hydrated 1,2-di-myristoyl-sn-glycero-3-phosphorylcholine bilayer molecular dynamics simulation

The conformation of nifedipine, a cardiac and smooth muscle calcium ion channel antagonist is studied in a hydrated bilayer of forty nine 1,2-di-myristoyl-sn-glycero-3-phosphorylcho-line (DMPC) molecules using molecular dynamics (MD) simulation technique. The simulation was carried out in conditions of constant number, volume and temperature (NVT) at 310 K, which is above the liquid crystalline (Lα) transition temperature of DMPC. The periodic boundary conditions were applied in three-dimensions. Thus the model represented an infinite bilayer. The important geometric parameters characteristic to DMPC and nifedipine molecules were calculated and compared with other theoretical and experimental results pertaining to nifedipine and other related dihydrophyridine (DHP) analogues. Our results suggest that conformational parameters required for antagonist activity are fairly conserved during the interaction of nifedipine with DMPC bilary and bilayer stabilizes the drug conformation in the bioactive form.     

Non-polar interactions between cholesterol and phospholipids: a molecular dynamics simulation study

A 15-ns molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. An 8-ns trajectory was analysed to investigate the effect of Chol on the chain packing in the bilayer core. While the packing of DMPC chains on the smooth alpha-face side of the Chol ring is similar to that in the pure DMPC bilayer, the packing on the rough beta-face side is less regular and less tight. Two methyl groups located on the Chol beta-face disturb the packing; in effect, van der Waals (vdW) interactions between Chol rings and DMPC chains are weaker than the ones between sole DMPC chains. VdW interactions between an alkyl chain of DMPC and an isooctyl tail of Chol are similarly strong as those between two DMPC chains.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

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.     

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.     

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.