Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Interactions of the Auxilin-1 PTEN-like Domain with Model Membranes Result in Nanoclustering of Phosphatidyl Inositol Phosphates

Auxilin-1 is a neuron-specific membrane-binding protein involved in a late stage of clathrin-mediated endocytosis. It recruits Hsc70, thus initiating uncoating of the clathrin-coated vesicles. Interactions of auxilin-1 with the vesicle membrane are crucial for this function and are mediated via an N-terminal PTEN-like domain. We have used multiscale molecular dynamics simulations to probe the interactions of the auxilin-1 PTEN-like domain with lipid bilayers containing differing phospholipid composition, including bilayers containing phosphatidyl inositol phosphates. Our results suggest a novel, to our knowledge, model for the auxilin/membrane encounter and subsequent interactions. Negatively charged lipids (especially PIP₂) enhance binding of auxilin to lipid bilayers and facilitate its correct orientation relative to the membrane. Mutations in three basic residues (R301E/R307E/K311E) of the C2 subdomain of the PTEN-like domain perturbed its interaction with the bilayer, changing its orientation. The interaction of membrane-bound auxilin-1 PTEN-like domain with negatively charged lipid headgroups results in nanoclustering of PIP₂ molecules in the adjacent bilayer leaflet.     

Lateral organization in lipid-cholesterol mixed bilayers

Interactions between lipid and cholesterol molecules in membranes play an important role in the structural and functional properties of cell membranes. Although structural properties of lipid-cholesterol mixtures have been extensively studied, an understanding of the role of cholesterol in the lateral organization of bilayers has been elusive. In this article, we propose a simple yet powerful model, based on self-consistent mean-field theory and molecular dynamics simulations, for lipid bilayers containing cholesterol. Properties predicted by our model are shown to be in excellent agreement with experimental data. Our model predicts that cholesterol induces structural changes in the bilayer through the formation of regions of ordered lipids surrounding each cholesterol molecule. We find that the "smooth" and "rough" sides of cholesterol play crucial roles in formation and distribution of the ordered regions. Our model is predictive in that input parameters are obtained from independent atomistic molecular dynamics simulations. The model and method are general enough to describe other heterogeneous lipid bilayers, including lipid rafts.     

Cholesterol surrogates: a comparison of cholesterol and 16:0 ceramide in POPC bilayers

Experimental evidence indicates that, under some circumstances, "surrogate" molecules may play the same role as cholesterol in ordering membrane lipids. The simplest molecule in this class is Ceramide. In this article, we describe atomic-level molecular dynamics simulations designed to shed light on this phenomenon. We run simulations of hydrated phosphoryl-oleoyl phosphatidylcholine (POPC) bilayers containing cholesterol, and containing ceramide, in concentrations ranging from 5% to 33%. We also perform a simulation of a pure POPC bilayer to verify the simulation force fields against experimental structural data for POPC. Our simulation data are in good agreement with experimental data for the partial molecular volumes, areas, form factors, and order parameters. These simulations suggest that ceramide and cholesterol have a very similar effect on the POPC bilayer, although ceramide is less effective in inducing order in the bilayer compared with cholesterol at the same concentrations.     

Molecular dynamics simulation studies of lipid bilayer systems

The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study required the construction of a model bacterial membrane which consisted of two types of phospholipids and counter ions. Whenever published experimental data were available, the results of the simulations were compared with them.     

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.     

The effect of surface charge density on valinomycin-K+ complex formation in model membranes

The model membrane approach was used to investigate the surface charge effect on the ion-antibiotic complexation process. Mixed monolayers of valinomycin and lipids were spread on subphases containing K+ or Na+. The surface charge density was modified by spreading ionizable valinomycin analogs on aqueous subphases of different pH or by changing the nature of the lipid (neutral, negatively charged) in the mixed film. Surface pressure and surface potential measurements demonstrated that a neutral lipid (phosphatidylcholine) or positively charged valinomycin analogs didn't enhance the anti-biotic complexing capacity. However, a maximal complexation is reached for a critical lipid concentration in the valinomycin-phosphatidylserine mixed film. The role of the surface charge on the valinomycin complexing properties was examined in terms of the Gouy-Chapman theory. As a consequence of the negative charge of the lipid monolayer, the K+ concentration near the surface is larger than the bulk concentration, by a Boltzmann factor. A good agreement was observed between the experimental results and the theoretical predictions. Conductance measurements of asymmetric bilayers containing a neutral lipid (egg lecithin) on one side and a negatively charged lipid (phosphatidyl-serine) on the other, confirm the role of the surface charge. Indeed, addition of K+ to the neutral side of the bilayer containing valinomycin had no effect on the conductance whereas addition of K+ to the charged side of the bilayer caused a 80-fold conductance increase.     

Combined monte carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.     

Molecular dynamics simulations of SOPS and sphingomyelin bilayers containing cholesterol

We performed a molecular dynamics simulation of an asymmetric bilayer that contained different lipid mixtures in its outer and inner leaflets. The outer leaflet contained a mixture of sphingomyelin (SM) with cholesterol and the inner leaflet a mixture of stearoyl-oleoyl-phosphatidylserine (SOPS) with cholesterol. For comparison purposes, we also performed two simulations on symmetric bilayers: the first simulation was performed on a bilayer containing a binary mixture of SOPS with cholesterol; the second contained a mixture of SM with cholesterol. We studied the hydrogen-bonding network of the bilayers in our simulations and the difference in the network properties in the monolayers either with SM or SOPS. We observed that in the asymmetric bilayer the properties of monolayers were the same as in the corresponding monolayers in the symmetric bilayers.     

Structure, dynamics, and hydration of POPC/POPS bilayers suspended in NaCl, KCl, and CsCl solutions

Effects of alkali metal chlorides on the properties of mixed negatively charged lipid bilayers are experimentally measured and numerically simulated. Addition of 20mol% of negatively charged phosphatidylserine to zwitterionic phosphatidylcholine strengthens adsorption of monovalent cations revealing their specificity, in the following order: Cs(+)相似文献   

Effects of a carane derivative local anesthetic on a phospholipid bilayer studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of two hydrated palmitoyloleoylphosphatidylcholine (POPC) bilayers each containing eight carane derivative (KP-23) local anesthetic (LA) molecules in neutral (POPC-LA) or protonated (POPC-LAH) forms were carried out to investigate the effect of KP-23 and its protonation on the bilayer. 3-ns trajectories were used for analyses. A pure POPC bilayer was employed as a reference system. In both POPC-LA and POPC-LAH systems a few KP-23 molecules intercalated into the bilayer and moved near the bilayer/water interface. They were located on the hydrophobic core side of the interface in the POPC-LA bilayer, but on the water phase side in the POPC-LAH bilayer. The order of the POPC chains was higher in the POPC-LA bilayer than in the pure POPC bilayer and was lower in the POPC-LAH bilayer. Interactions between polar groups of KP-23 and POPC or water were responsible for a lower hydration of POPC headgroups in POPC bilayers containing KP-23 than in the pure POPC bilayer. KP-23 molecules were found to form aggregates both in POPC-LA and POPC-LAH bilayers. Due to higher amphiphilicity of LAH, the LAH aggregate was more micelle-like and larger than the LA one. The results demonstrate the rapid timescales of the initial processes that take place at and near the bilayer interface as well as details of the atomic level interactions between local anesthetic and the lipid matrix of a cell 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.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Structure of Sphingomyelin Bilayers: A Simulation Study

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50°C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.     

Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.     

The effect of surface charge density on valinomycin-K+ complex formation in model membranes

The model membrane approach was used to investigate the surface charge effect on the ion-antibiotic complexation process. Mixed monolayers of valinomycin and lipids were spread on subphases containing K⁺ or Na⁺. The surface charge density was modified by spreading ionizable valinomycin analogs on aqueous subphases of different pH or by changing the nature of the lipid (neutral, negatively charged) in the mixed film. Surface pressure and surface potential measurements demonstrated that a neutral lipid (phosphatidylcholine) or positively charged valinomycin analogs didn't enhance the antibiotic complexing capacity. However, a maximal complexation is reached for a critical lipid concentration in the valinomycin-phosphatidylserine mixed film. The role of the surface charge on the valinomycin complexing properties was examined in terms of the Gouy-Chapman theory. As a consequence of the negative charge of the lipid monolayer, the K⁺ concentration near the surface is larger than the bulk concentration, by a Boltzmann factor. A good agreement was observed between the experimental results and the theoretical predictions. Conductance measurements of asymmetric bilayers containing a neutral lipid (egg lecithin) on one side and a negatively charged lipid (phosphatidylserine) on the other, confirm the role of the surface charge. Indeed, addition of K⁺ to the neutral side of the bilayer containing valinomycin had no effect on the conductance whereas addition of K⁺ to the charged side of the bilayer caused a 80-fold conductance increase.     

Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.     

Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.     

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.