Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration

We have applied a hybrid equilibration and sampling procedure for the atomic level simulation of a hydrated lipid bilayer to systems consisting of dipalmitoyl phosphatidylcholine (DPPC) and cholesterol, and palmitoyl-oleyl phosphatidylcholine (POPC) at low (approximately 6%) cholesterol concentration. The procedure is applied to bilayers of 94 molecules of DPPC, 6 molecules of cholesterol, and 3205 water molecules, and to bilayers of 120 molecules of POPC, 8 molecules of cholesterol, and 4268 water molecules, at a temperature of 325 K. After equilibration, three separate 400-ps continuous molecular dynamics runs, separated by 10,000 configurational bias Monte Carlo 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 DPPC and cholesterol, and of pure DPPC, and pure POPC. Certain properties of the bilayers are indistinguishable from cholesterol-free bilayers, including lateral diffusion and electron density. Other properties, most notably the order parameter profile, show the effect of cholesterol even at low concentrations.     

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.     

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.     

Plasmon-waveguide resonance studies of lateral segregation of lipids and proteins into microdomains (rafts) in solid-supported bilayers

Plasmon-waveguide resonance (PWR) spectroscopy has been used to examine solid-supported lipid bilayers consisting of dioleoylphosphatidylcholine (DOPC), palmitoyloleoylphosphatidylcholine (POPC), sphingomyelin (SM), and phosphatidylcholine/SM binary mixtures. Spectral simulation of the resonance curves demonstrated an increase in bilayer thickness, long-range order, and molecular packing density in going from DOPC to POPC to SM single component bilayers, as expected based on the decreasing level of unsaturation in the fatty acyl chains. DOPC/SM and POPC/SM binary mixtures yielded PWR spectra that can be ascribed to a superposition of two resonances corresponding to microdomains (rafts) consisting of phosphatidylcholine- and SM-rich phases coexisting within a single bilayer. These were formed spontaneously over time as a consequence of lateral phase separation. Each microdomain contained a small proportion (<20%) of the other lipid component, which increased their kinetic and thermodynamic stability. Incorporation of a glycosylphosphatidylinositol-linked protein (placental alkaline phosphatase) occurred within each of the single component bilayers, although the insertion was less efficient into the DOPC bilayer. Incorporation of placental alkaline phosphatase into a DOPC/SM binary bilayer occurred with preferential insertion into the SM-rich phase, although the protein incorporated into both phases at higher concentrations. These results demonstrate the utility of PWR spectroscopy to provide insights into raft formation and protein sorting in model lipid membranes.     

Influence of stigmastanol and stigmastanyl-phosphorylcholine, two plasma cholesterol lowering substances, on synthetic phospholipid membranes. A 2H- and 31P-NMR study.

Cholesterol, stigmastanol, and stigmastanyl-phosphorylcholine (ST-PC) were incorporated into model membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). POPC and ST-PC were deuterated at the lipid headgroup, DOPC at the cis-double bonds. The influence of the three sterols on the motion and conformation of the lipid headgroups and the hydrocarbon chains was monitored with 2H- and 31P-NMR. All three sterols were freely miscible with the lipid matrix in concentrations of up to 50 mol% without inducing phase separations or nonbilayer structures. However, the molecules exert quite different effects on the phospholipid bilayer. Cholesterol and stigmastanol are largely buried in the hydrocarbon part of the membrane, distinctly restricting the flexing motions of the fatty acyl chains whereas the conformation of the phospholipid headgroups is little affected. In contrast, ST-PC is anchored with its headgroup in the layer of phospholipid dipoles, preventing an extensive penetration of the sterol ring into the hydrocarbon layer. Hence ST-PC has almost no effect on the hydrocarbon chains but induces a characteristic conformational change of the phospholipid headgroups. The 2H- and 31P-NMR spectra of mixed phospholipid/ST-PC membranes further demonstrate that the PC headgroup of ST-PC has a similar orientation as the surrounding phosphatidylcholine headgroups. For both types of molecules the -P-N+ dipole is essentially parallel to the membrane surface. Addition of ST-PC induces a small rotation of the POPC headgroup towards the water phase.     

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 study of interaction of dimyristoyl phosphotidyl choline with water

We present here results of molecular dynamics (MD) simulations on hydrated bilayers of 40 molecules of 1-2-dimyristoyl-sn-glycero-3-phosphatidyl choline (DMPC) in liquid crystalline (Lα) phase using two different models (i) with same (A) conformation for all DMPC molecules, (ii) with alternate rows having different (A and B reported in crystallographic studies on DMPC) conformations. The bilayers were hydrated using 776 and 1064 water molecules. Simulations have been carried out at 310K with AMBER 4.0 program, using united atom force field for 200 pico seconds (ps) after equilibration. During heating and equilibration constant pressure temperature (PT) conditions were maintained while in simulation of equillibrated bilayers constant volume temperature (VT) conditions were used. Subaveraged atomic coordinates were used to calculate geometric parameters of lipid molecules and lipid water interaction. Our results show larger flexibility of polar head group and glycerol region in Lα phase compared to gel or non-hydrated bilayers. Chain disorder was more towards end. Sn-2 chains were more disordered. Use of two types of starting conformations increased disorder. Trans fraction of chain torsional angle was higher in non-hydrated bilayer. However it was more disordered due to ‘swing’ movement of chains because of distortion in torsional angles α2 and 03 due to absence of water molecules. Trans fraction of the chains, order parameter and water penetration showed general agreement with the available experimental results. On the whole MD technique was found to be quite useful for depicting microscopic behaviour of liquid crystalline system and correlating the same with macroscopic changes observed experimentally.     

Diffraction-based density restraints for membrane and membrane-peptide molecular dynamics simulations

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.     

Interfacial tryptophan residues: a role for the cation-pi effect?

Integral membrane proteins are characterized by having a preference for aromatic residues, e.g., tryptophan (W), at the interface between the lipid bilayer core and the aqueous phase. The reason for this is not clear, but it seems that the preference is related to a complex interplay between steric and electrostatic forces. The flat rigid paddle-like structure of tryptophan, associated with a quadrupolar moment (aromaticity) arising from the pi-electron cloud of the indole, interacts primarily with moieties in the lipid headgroup region hardly penetrating into the bilayer core. We have studied the interaction between the nitrogen moiety of lipid molecule headgroups and the pi-electron distribution of gramicidin (gA) tryptophan residues (W9, W11, W13, and W15) using molecular dynamics (MD) simulations of gA embedded in two hydrated lipid bilayers composed of 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) and 1-palmitoyl-2-oleoylphosphatidyl-choline (POPC), respectively. We use a force field model for tryptophan in which polarizability is only implicit, but we believe that classical molecular dynamics force fields are sufficient to capture the most prominent features of the cation-pi interaction. Our criteria for cation-pi interactions are based on distance and angular requirements, and the results from our model suggest that cation-pi interactions are relevant for W(PE)1), W(PE)13, W(PE)15, and, to some extent, W(PC)11 and W(PC)13. In our model, W9 does not seem to engage in cation-pi interactions with lipids, neither in POPE nor POPC. The criteria for the cation-pi effect are satisfied more often in POPE than in POPC, whereas the H-bonding ability between the indole donor and the carbonyl acceptor is similar in POPE and POPC. This suggests an increased affinity for lipids with ethanolamine headgroups to transmembrane proteins enriched in interfacial tryptophans.     

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).     

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.     

Location of the TEMPO moiety of TEMPO-PC in phosphatidylcholine bilayers is membrane phase dependent

The (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) moiety tethered to the headgroup of phosphatidylcholine (PC) lipid is employed in spin labeling electron paramagnetic resonance spectroscopy to probe the water dynamics near lipid bilayer interfaces. Due to its amphiphilic character, however, TEMPO spin label could partition between aqueous and lipid phases, and may even be stabilized in the lipid phase. Accurate assessment of the TEMPO-PC configuration in bilayer membranes is essential for correctly interpreting the data from measurements. Here, we carry out all-atom molecular dynamics (MD) simulations of TEMPO-PC probe in single-component lipid bilayers at varying temperatures, using two standard MD force fields. We find that, for a dipalmitoylphosphatidylcholine (DPPC) membrane whose gel-to-fluid lipid phase transition occurs at 314 K, while the TEMPO spin label is stabilized above the bilayer interface in the gel phase, there is a preferential location of TEMPO below the membrane interface in the fluid phase. For bilayers made of unsaturated lipids, 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), which adopt the fluid phase at ambient temperature, TEMPO is unequivocally stabilized inside the bilayers. Our finding of membrane phase-dependent positioning of the TEMPO moiety highlights the importance of assessing the packing order and fluidity of lipids under a given measurement condition.     

Effect of cholesterol on the structure and dynamic properties of unsaturated phospholipid bilayers

Molecular dynamics (MD) computer simulations of five different hydrated unsaturated phosphatidylcholine lipid bilayers built up by 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules with 40 mol% cholesterol, and the same five pure phosphatidylcholine bilayers have been performed at 303 K. The simulation box of a lipid bilayer contained 96 phosphatidylcholines, 64 cholesterols, and 3840 water molecules (48 phosphatidylcholine molecules and 32 cholesterols per layer and 24 water molecules per phospholipid or cholesterol in each case). The lateral self-diffusion coefficients of the lipids in these systems and mass density profiles with respect to the bilayer normal have been analyzed. It has been found that the lateral diffusion coefficients of phosphatidylcholine molecules increase with increasing number of double bonds in one of the lipid chains, both in pure bilayers and in bilayers with cholesterol. It has been found as well that the lateral diffusion coefficient of phosphatidylcholine molecules of a lipid bilayer with 40 mol% cholesterol is smaller than that for the corresponding pure phosphatidylcholine bilayer.     

Molecular interpretation of fluorescence solvent relaxation of Patman and 2H NMR experiments in phosphatidylcholine bilayers

The analysis of time-dependent fluorescence shifts of the bilayer probe 6-hexadecanoyl-2-(((2-(trimethylammonium)ethyl)methyl)amino)naphthalene chloride (Patman) offers valuable information on the hydration and dynamics of phospholipid headgroups. Quenching studies on vesicles composed of four phosphatidylcholines with different hydrocarbon chains (18:1c9/18:1c9, DOPC; 16:0/18:1c9, POPC; 18:1c9/16:0, OPPC; 18:1c6/18:1c6, PCDelta6) show that the chromophore of Patman is defined located at the level of the sn-1 ester-group in the phospholipid, which is invariant to the hydrocarbon chain. The so-called solvent relaxation (SR) approach as well as solid-state 2H NMR reveals that DOPC and PCDelta6 are more hydrated than POPC and OPPC. A strong dependence of SR kinetics on the position of double bond in the investigated fatty acid chains was observed. Apparently, the closer the double bond is located to the hydrated sn-1 ester-group, the more mobile this group becomes. This work demonstrates that the SR approach can report mobility changes within phospholipid bilayers with a remarkable molecular resolution.     

Phospholipase A2 hydrolysis of supported phospholipid bilayers: a neutron reflectivity and ellipsometry study

We have investigated the phospholipase A(2) catalyzed hydrolysis of supported phospholipid bilayers using neutron reflection and ellipsometry. At the hydrophilic silica-water interface, hydrolysis of phosphatidylcholine bilayers by phospholipase A(2) from Naja mossambica mossambica venom is accompanied by destruction of the bilayer at an initial rate, which is comparable for DOPC and DPPC but is doubled for POPC. The extent of bilayer destruction at 25 degrees C decreases from DOPC to POPC and is dramatically reduced for DPPC. Neutron reflectivity measurements indicate that the enzyme penetrates into the bilayers in increasing order for DOPC, POPC, and DPPC, while the amount of enzyme adsorbed at the interface is smallest for DPPC and exhibits a maximum for POPC. Penetration into the hydrophobic chain region in the bilayer is further supported by the fact that the enzyme adsorbs strongly and irreversibly to a hydrophobic monolayer of octadecyltrichlorosilane. These results are rationalized in terms of the properties of the reaction products and the effect of their accumulation in the membrane on the kinetics of enzyme catalysis.     

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.     

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25?ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.     

Membrane interface composition drives the structure and the tilt of the single transmembrane helix protein PMP1: MD studies

PMP1, a regulatory subunit of the yeast plasma membrane H⁺-ATPase, is a single transmembrane helix protein. Its cytoplasmic C-terminus possesses several positively charged residues and interacts with phosphatidylserine lipids as shown through both ¹H- and ²H-NMR experiments. We used all-atom molecular dynamics simulations to obtain atomic-scale data on the effects of membrane interface lipid composition on PMP1 structure and tilt. PMP1 was embedded in two hydrated bilayers, differing in the composition of the interfacial region. The neutral bilayer is composed of POPC (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylcholine) lipids and the negatively charged bilayer is composed of POPC and anionic POPS (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylserine) lipids. Our results were consistent with NMR data obtained previously, such as a lipid sn-2 chain lying on the W28 aromatic ring and in the groove formed on one side of the PMP1 helix. In pure POPC, the transmembrane helix is two residues longer than the initial structure and the helix tilt remains constant at 6 ± 3°. By contrast, in mixed POPC-POPS, the initial helical structure of PMP1 is stable throughout the simulation time even though the C-terminal residues interact strongly with POPS headgroups, leading to a significant increase of the helix tilt within the membrane to 20 ± 5°.     

Lipid bilayer molecular dynamics study of lipid-derived agonists of the putative cannabinoid receptor, GPR55

Both L-α-lysophosphatidylinositol (LPI) and 2-arachidonoyl-sn-glycero-3-phosphoinositol (2-AGPI) have been reported to activate the putative cannabinoid receptor, GPR55. Recent microsecond time-scale molecular dynamics (MD) simulations and isothiocyanate covalent labeling studies have suggested that a transmembrane helix 6/7 (TMH6/7) lipid pathway for ligand entry may be necessary for interaction with cannabinoid receptors. Because LPI and 2-AGPI are lipid-derived ligands, conformations that each assumes in the lipid bilayer are therefore likely important for their interaction with GPR55. We report here the results of 70 ns NAMD molecular dynamics (MD) simulations of LPI and of 2-AGPI in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). These simulations are compared with a 70 ns simulation of the cannabinoid CB1 receptor endogenous ligand, N-arachidonoylethanolamine (anandamide, AEA) in a POPC bilayer. These simulations revealed that (1) LPI and 2-AGPI sit much higher in the bilayer than AEA, with inositol headgroups that can at times be solvated completely by water; (2) the behavior of the acyl chains of AEA and 2-AGPI are similar in their flexibilities in the bilayer, while the acyl chain of LPI has reduced flexibility; and (3) both 2-AGPI and LPI can adopt a tilted headgroup orientation by hydrogen bonding to the phospholipid phosphate/glycerol groups or via intramolecular hydrogen bonding. This tilted head group conformation (which represents over 40% of the conformer population of LPI (42.2 ± 3.3%) and 2-AGPI (43.7 ± 1.4%)) may provide a low enough profile in the lipid bilayer for LPI and 2-AGPI to enter GPR55 via the putative TMH6/7 entry port.     

Partitioning of anesthetics into a lipid bilayer and their interaction with membrane-bound peptide bundles

Molecular dynamics simulations have been performed to investigate the partitioning of the volatile anesthetic halothane from an aqueous phase into a coexisting hydrated bilayer, composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipids, with embedded alpha-helical peptide bundles based on the membrane-bound portions of the alpha- and delta-subunits, respectively, of nicotinic acetylcholine receptor. In the molecular dynamics simulations halothane molecules spontaneously partitioned into the DOPC bilayer and then preferentially occupied regions close to lipid headgroups. A single halothane molecule was observed to bind to tyrosine (Tyr-277) residue in the alpha-subunit, an experimentally identified specific binding site. The binding of halothane attenuated the local loop dynamics of alpha-subunit and significantly influenced global concerted motions suggesting anesthetic action in modulating protein function. Steered molecular dynamics calculations on a single halothane molecule partitioned into a DOPC lipid bilayer were performed to probe the free energy profile of halothane across the lipid-water interface and rationalize the observed spontaneous partitioning. Partitioned halothane molecules affect the hydrocarbon chains of the DOPC lipid, by lowering of the hydrocarbon tilt angles. The anesthetic molecules also caused a decrease in the number of peptide-lipid contacts. The observed local and global effects of anesthetic binding on protein motions demonstrated in this study may underlie the mechanism of action of anesthetics at a molecular level.