Molecular dynamics simulations of rhodopsin in different one-component lipid bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.     

Effect of ions on a dipalmitoyl phosphatidylcholine bilayer. a molecular dynamics simulation study

The effect of physiological concentrations of different chlorides on the structure of a dipalmitoyl phosphatidylcholine (DPPC) bilayer has been investigated through atomistic molecular dynamics simulations. These calculations provide support to the concept that Li+, Na+, Ca2+, Mg2+, Sr2+, Ba2+, and Ac3+, but not K+, bind to the lipid-head oxygens. Ion binding exhibits an influence on lipid order, area per lipid, orientation of the lipid head dipole, the charge distribution in the system, and therefore the electrostatic potential across the head-group region of the bilayer. These structural effects are sensitive to the specific characteristics of each cation, i.e., radius, charge, and coordination properties. These results provide evidence aimed at shedding some light into the apparent contradictions among different studies reported recently regarding the ordering effect of ions on zwitterionic phosphatidylcholine lipid bilayers.     

Fusogenic tilted peptides induce nanoscale holes in supported phosphatidylcholine bilayers

Tilted peptides are known to insert in lipid bilayers with an oblique orientation, thereby destabilizing membranes and facilitating membrane fusion processes. Here, we report the first direct visualization of the interaction of tilted peptides with lipid membranes using in situ atomic force microscopy (AFM) imaging. Phase-separated supported dioleoylphosphatidylcholine/dipalmitoylphosphatidylcholine (DOPC/DPPC) bilayers were prepared by fusion of small unilamellar vesicles and imaged in buffer solution, in the absence and in the presence of the simian immunodeficiency virus (SIV) peptide. The SIV peptide was shown to induce the rapid appearance of nanometer scale bilayer holes within the DPPC gel domains, while keeping the domain shape unaltered. We attribute this behavior to a local weakening and destabilization of the DPPC domains due to the oblique insertion of the peptide molecules. These results were directly correlated with the fusogenic activity of the peptide as determined using fluorescently labeled DOPC/DPPC liposomes. By contrast, the nontilted ApoE peptide did not promote liposome fusion and did not induce bilayer holes but caused slight erosion of the DPPC domains. In conclusion, this work provides the first direct evidence for the production of stable, well-defined nanoholes in lipid bilayer domains by the SIV peptide, a behavior that we have shown to be specifically related to the tilted character of the peptide. A molecular mechanism underlying spontaneous insertion of the SIV peptide within lipid bilayers and the subsequent removal of bilayer patches is proposed, and its relevance to membrane fusion processes is discussed.     

Material properties of matrix lipids determine the conformation and intermolecular reactivity of diacetylenic phosphatidylcholine in the lipid bilayer

Photopolymerizable phospholipid DC(8,9)PC (1,2-bis-(tricosa-10,12-diynoyl)-sn-glycero-3-phosphocholine) exhibits unique assembly characteristics in the lipid bilayer. Because of the presence of the diacetylene groups, DC(8,9)PC undergoes polymerization upon UV (254 nm) exposure and assumes chromogenic properties. DC(8,9)PC photopolymerization in gel-phase matrix lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) monitored by UV-vis absorption spectroscopy occurred within 2 min after UV treatment, whereas no spectral shifts were observed when DC(8,9)PC was incorporated into liquid-phase matrix 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Liquid chromatography-tandem mass spectrometry analysis showed a decrease in the amount of DC(8,9)PC monomer in both DPPC and POPC environments without any change in the matrix lipids in UV-treated samples. Molecular dynamics (MD) simulations of DPPC/DC(8,9)PC and POPC/DC(8,9)PC bilayers indicate that the DC(8,9)PC molecules adjust to the thickness of the matrix lipid bilayer. Furthermore, the motions of DC(8,9)PC in the gel-phase bilayer are more restricted than in the fluid bilayer. The restricted motional flexibility of DC(8,9)PC (in the gel phase) enables the reactive diacetylenes in individual molecules to align and undergo polymerization, whereas the unrestricted motions in the fluid bilayer restrict polymerization because of the lack of appropriate alignment of the DC(8,9)PC fatty acyl chains. Fluorescence microscopy data indicates the homogeneous distribution of lipid probe 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-lissamine rhodamine B sulfonyl ammonium salt (N-Rh-PE) in POPC/DC(8,9)PC monolayers but domain formation in DPPC/DC(8,9)PC monolayers. These results show that the DC(8,9)PC molecules cluster and assume the preferred conformation in the gel-phase matrix for the UV-triggered polymerization reaction.     

Structure and thermotropic phase behavior of fluorinated phospholipid bilayers: a combined attenuated total reflection FTIR spectroscopy and imaging ellipsometry study

Lipid bilayers consisting of lipids with terminally perfluoroalkylated chains have remarkable properties. They exhibit increased stability and phase-separated nanoscale patterns in mixtures with nonfluorinated lipids. In order to understand the bilayer properties that are responsible for this behavior, we have analyzed the structure of solid-supported bilayers composed of 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC) and of a DPPC analogue with 6 terminal perfluorinated methylene units (F6-DPPC). Polarized attenuated total reflection Fourier-transform infrared spectroscopy indicates that for F6-DPPC, the tilt of the lipid acyl chains to the bilayer normal is increased to 39 degrees as compared to 21 degrees for native DPPC, for both lipids in the gel phase. This substantial increase of the tilt angle is responsible for a decrease of the bilayer thickness from 5.4 nm for DPPC to 4.5 nm for F6-DPPC, as revealed by temperature-controlled imaging ellipsometry on microstructured lipid bilayers and solution atomic force microscopy. During the main phase transition from the gel to the fluid phase, both the relative bilayer thickness change and the relative area change are substantially smaller for F6-DPPC than for DPPC. In light of these structural and thermotropic data, we propose a model in which the higher acyl-chain tilt angle in F6-DPPC is the result of a conformational rearrangement to minimize unfavorable fluorocarbon-hydrocarbon interactions in the center of the bilayer due to chain staggering.     

Engineered lipids that cross-link the inner and outer leaflets of lipid bilayers

The application of supported lipid bilayer systems as molecular sensors, diagnostic devices, and medical implants is limited by their lack of stability. In an effort to enhance the stability of supported lipid bilayers, three pairs of phosphatidylcholine lipids were designed to cross-link at the termini of their 2-position acyl chain upon the formation of lipid bilayers. The cross-linked lipids span the lipid bilayer, resembling naturally occurring bolaamphiphiles that stabilize archaebacterial membranes against high temperatures. The three reactions investigated here include the acyl chain cross-linking between thiol and bromine groups, thiol and acryloyl groups, and cyclopentadiene and acryloyl groups. All three reactive lipid pairs were found to cross-link in liposomal membranes, as determined by thin-layer chromatography, ion-spray mass spectrometry, and 1H NMR. The monolayer film properties of the reactive amphiphiles were characterized by surface pressure-area isotherms and showed that stable monolayers formed at the air-water interface with limiting molecular areas comparable to that of pure saturated phosphatidylcholine lipids. Langmuir-Blodgett bilayers of dimyristoylphosphatidylcholine incorporating 15 mol % of the reactive thiol and acryloyl lipids had diffusion coefficients comparable with pure dimyristoylphosphatidylcholine, while bilayers with more than 25 mol % of the reactive lipids were immobile, suggesting that interleaflet cross-linking of the lipids inhibited membrane diffusion. Our results show that the reactive lipids can cross-link within a lipid bilayer and are suitable for assembling supported lipid bilayers using Langmuir-Blodgett deposition. By using terminally reactive amphiphiles to build up supported lipid bilayers with cross-linked leaflets, bolaamphiphiles can be incorporated into asymmetric solid supported membranes to increase their stability in biosensor and medical implant applications.     

Molecular dynamics study on the stabilization of dehydrated lipid bilayers with glucose and trehalose

In this study, we use molecular dynamics simulations to investigate and compare the interactions of DPPC bilayers with and without saccharides (glucose or trehalose) under dehydrated conditions. Results from the simulations indicate that unilamellar bilayers lose their structural integrity under dehydrated conditions in the absence of saccharides; however, in the presence of either glucose or trehalose, the bilayers maintain their stability. Hydrogen bond analysis shows that the saccharide molecules displace a significant amount of water surrounding the lipid headgroups. At the same time, the additional hydrogen bonds formed between water and saccharide molecules help to maintain a hydration layer on the lipid bilayer interface. On the basis of the hydrogen bond distributions, trehalose forms more hydrogen bonds with the lipids than glucose, and it is less likely to interact with neighboring saccharide molecules. These results suggest that the interaction between the saccharide and lipid molecules through hydrogen bonds is an essential component of the mechanism for the stabilization of lipid bilayers.     

Molecular dynamics studies of the molecular structure and interactions of cholesterol superlattices and random domains in an unsaturated phosphatidylcholine bilayer membrane

The effect of the molecular organization of lipid components on the properties of the bilayer membrane has been a topic of increasing interest. Several experimental and theoretical studies have suggested that cholesterol is not randomly distributed in the fluid-state lipid bilayer but forms nanoscale domains. Several cholesterol-enriched nanodomain structures have been proposed, including rafts, regular or maze arrays, complexes, and superlattices. At present, the molecular mechanisms by which lipid composition influences the formation and stability of lipid nanodomains remain unclear. In this study, we have used molecular dynamics (MD) simulations to investigate the effects of the molecular organization of cholesterol--superlattice versus random--on the structure of and interactions between lipids and water in lipid bilayers of cholesterol and 1-palmitoyl-2-oleoylphosphatidylcholine (cholesterol/POPC) at a fixed cholesterol mole fraction of 0.40. On the basis of four independent replicates of 200-ns MD simulations for a superlattice or random bilayer, statistically significant differences were observed in the lipid structural parameters, area per lipid, density profile, and acyl chain order profile, as well as the hydrogen bonding between various pairs (POPC and water, cholesterol and water, and POPC and cholesterol). The time evolution of the radial distribution of the cholesterol hydroxy oxygen suggests that the lateral distribution of cholesterol in the superlattice bilayer is more stable than that in the random bilayer. Furthermore, the results indicate that a relatively long simulation time, more than 100 ns, is required for these two-component bilayers to reach equilibrium and that this time is influenced by the initial lateral distribution of lipid components.     

Cholesterol packing around lipids with saturated and unsaturated chains: a simulation study

The fundamental role of cholesterol in the regulation of eukaryotic membrane structure is well-established. However the manner in which atomic level interactions between cholesterol and lipids, with varying degrees of chain unsaturation and polar groups, affect the overall structure and organization of the bilayer is only beginning to be understood. In this paper we describe a series of Molecular Dynamics simulations designed to provide new insights into lipid-cholesterol interactions as a function of chain unsaturation. We have run simulations of varying concentrations of cholesterol in dipalmitoyl phosphatidylcholine (DPPC), palmitoyl-oleyol phosphatidylcholine (POPC), and dioleyol phosphatidylcholine (DOPC) bilayers. Structural analysis of the simulations reveals both atomistic and systemic details of the interactions and are presented here. In particular, we find that the minimum partial molecular area of cholesterol occurs in POPC-Chol mixtures implying the most favorable packing. Physically, this appears to be related to the fact that the two faces of the cholesterol molecule are different from each other and that the steric cross section of cholesterol molecules drops sharply near the small chain tails.     

Structural studies of the monolayers and bilayers formed by a novel cholesterol-phospholipid chimera

Langmuir isotherm, neutron reflectivity, and small angle neutron scattering studies have been conducted to characterize the monolayers and vesicular bilayers formed by a novel chimeric phospholipid, ChemPPC, that incorporates a cholesteryl moeity and a C-16 aliphatic chain, each covalently linked via a glycerol backbone to phosphatidylcholine. The structures of the ChemPPC monolayers and bilayers are compared against those formed from pure dipalmitoylphoshatidylcholine (DPPC) and those formed from a 60:40 mol % mixture of DPPC and cholesterol. In accord with previous findings showing that very similar macroscopic properties were exhibited by ChemPPC and 60:40 mol % DPPC/cholesterol vesicles, it is found here that the chimeric lipid and lipid/sterol mixture have very similar monolayer structures (each having a monolayer thickness of ～26 ?), and they also form vesicles with similar lamellar structure, each having a bilayer thickness of ～50 ? and exhibiting a repeat spacing of ～65 ?. The interfacial area of ChemPPC, however, is around 10 ?(2) greater than that of the combined DPPC/cholesterol unit in the mixed lipid monolayer (viz., 57 ± 1 vs 46 ± 1 ?(2), at 35 mN·m(-1)), and this difference in area is attributed to the succinyl linkage which joins the ChemPPC steroid and glyceryl moieties. The larger area of the ChemPPC is reflected in a slightly thicker monolayer solvent distribution width (9.5 vs 9 ? for the DPPC/cholesterol system) and by a marginal increase in the level of lipid headgroup hydration (16 vs 13 H(2)O per lipid, at 35 mN·m(-1)).     

The Role of Lateral Tension in Calcium-Induced DPPS Vesicle Rupture

We assess the role of lateral tension in rupturing anionic dipalmitoylphosphatidyserine (DPPS), neutral dipalmitoylphosphatidylcholine (DPPC), and mixed DPPS-DPPC vesicles. Binding of Ca(2+) is known to have a significant impact on the effective size of DPPS lipids and little effect on the size of DPPC lipids in bilayer structures. In the present work we utilized laser transmission spectroscopy (LTS) to assess the effect of Ca(2+)-induced stress on the stability of the DPPS and DPPC vesicles. The high sensitivity and resolution of LTS has permitted the determination of the size and shape of liposomes in solution. The results indicate a critical size after which DPPS single shell vesicles are no longer stable. Our measurements indicate Ca(2+) promotes bilayer fusion up to a maximum diameter of ca. 320 nm. These observations are consistent with a straightforward free-energy-based model of vesicle rupture involving lateral tension between lipids regulated by the binding of Ca(2+). Our results support a critical role of lateral interactions within lipid bilayers for controlling such processes as the formation of supported bilayer membranes and pore formation in vesicle fusion. Using this free energy model we are able to infer a lower bound for the area dilation modulus for DPPS (252 pN/nm) and demonstrate a substantial free energy increase associated with vesicle rupture.     

Using bicellar mixtures to form supported and suspended lipid bilayers on silicon chips

Bicellar mixtures, planar lipid bilayer assemblies comprising long- and short-chain phosphatidylcholine lipids in suspension, were used to form supported lipid bilayers on flat silicon substrate and on nanotextured silicon substrates containing arrays of parallel troughs (170 nm wide, 380 nm deep, and 300 nm apart). Confocal fluorescence and atomic force microscopies were used to characterize the resulting lipid bilayer. Formation of a continuous biphasic undulating lipid bilayer membrane, where the crests and troughs corresponded to supported and suspended lipid bilayer regions, is demonstrated. The use of interferometric lithography to fabricate nanotexured substrates provides an advantage over other nanotextured substrates such as nanoporous alumina by offering flexibility in designing different geometries for suspending lipid bilayers.     

Atomistic simulation of lipid and DiI dynamics in membrane bilayers under tension

Membrane tension modulates cellular processes by initiating changes in the dynamics of its molecular constituents. To quantify the precise relationship between tension, structural properties of the membrane, and the dynamics of lipids and a lipophilic reporter dye, we performed atomistic molecular dynamics (MD) simulations of DiI-labeled dipalmitoylphosphatidylcholine (DPPC) lipid bilayers under physiological lateral tensions ranging from -2.6 mN m(-1) to 15.9 mN m(-1). Simulations showed that the bilayer thickness decreased linearly with tension consistent with volume-incompressibility, and this thinning was facilitated by a significant increase in acyl chain interdigitation at the bilayer midplane and spreading of the acyl chains. Tension caused a significant drop in the bilayer's peak electrostatic potential, which correlated with the strong reordering of water and lipid dipoles. For the low tension regime, the DPPC lateral diffusion coefficient increased with increasing tension in accordance with free-area theory. For larger tensions, free area theory broke down due to tension-induced changes in molecular shape and friction. Simulated DiI rotational and lateral diffusion coefficients were lower than those of DPPC but increased with tension in a manner similar to DPPC. Direct correlation of membrane order and viscosity near the DiI chromophore, which was just under the DPPC headgroup, indicated that measured DiI fluorescence lifetime, which is reported to decrease with decreasing lipid order, is likely to be a good reporter of tension-induced decreases in lipid headgroup viscosity. Together, these results offer new molecular-level insights into membrane tension-related mechanotransduction and into the utility of DiI in characterizing tension-induced changes in lipid packing.     

Composition of supported model membranes determined by neutron reflection

We have investigated the formation of supported bilayers by coadsorption of dipalmitoyl phosphatidylcholine (DPPC) with the nonionic surfactant beta-D-dodecyl maltoside. The adsorption of mixed phospholipid-surfactant micelles on hydrophilic silica surfaces at 25 degrees C was followed as a function of bulk concentration by neutron reflection. Using chain-deuterated d(25)-beta-D-dodecyl maltoside and d(62)-DPPC, we demonstrate that it is possible to determine the composition of the bilayers at each stage of a sequential dilution process, which enriches the adsorbed layer in phospholipid and leads to complete elimination of the surfactant. The final supported bilayers have thicknesses of 51 +/- 3 A and are stable to heating to 37 degrees C once all surfactant has been removed, and the structures agree well with other published data on DPPC supported bilayers. The coadsorption of cholesterol in a DPPC-surfactant mixture was also achieved, and the location and volume fraction of cholesterol in the DPPC bilayer was determined. Cholesterol is located in a 18 +/- 1 A thick layer below the lipid headgroup region and leads to an increased bilayer thickness of 58 +/- 2 A at 26 mol % of cholesterol.     

Cholesterol induced suppression of large swelling of water layer in phosphocholine floating bilayers

The effect of the addition of 10 mol % cholesterol to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) floating bilayers has been investigated by means of neutron reflectivity measurements. The large swelling of the water layer between the two bilayers found in pure phosphocholine systems around the lipid main phase transition is greatly reduced in the 9:1 DPPC:cholesterol mixture. The analysis of the structure of the bilayer reveals that in the gel phase cholesterol induces the presence of a high rms roughness that disappears in the fluid phase.     

Atomistic simulation of cholesterol effects on miscibility of saturated and unsaturated phospholipids: implications for liquid-ordered/liquid-disordered phase coexistence

Mixed MD/MC simulation at fixed difference in chemical potential (Δμ) between two lipid types provides a computational indicator of the relative affinities of the two lipids for different environments. Applying this technique to ternary DPPC/DOPC/cholesterol bilayers yields a DPPC/DOPC ratio that increases with increasing cholesterol content at fixed Δμ, consistent with the known enrichment of DPPC and cholesterol-rich in liquid-ordered phase domains in the fluid-fluid coexistence region of the ternary phase diagram. Comparison of the cholesterol-dependence of PC compositions at constant Δμ with experimentally measured coexistence tie line end point compositions affords a direct test of the faithfulness of the atomistic model to experimental phase behavior. DPPC/DOPC ratios show little or no dependence on cholesterol content at or below 16% cholesterol in the DOPC-rich region of the composition diagram, indicating cooperativity in the favorable interaction between DPPC and cholesterol. The relative affinity of DPPC and DOPC for high cholesterol bilayer environments in simulations is explicitly shown to depend on the degree of cholesterol alignment with the bilayer normal, suggesting that a source of the cooperativity is the composition dependence of cholesterol tilt angle distributions.     

Micropatterned composite membranes of polymerized and fluid lipid bilayers

Micropatterned composite membranes of polymerized and fluid lipid bilayers were constructed on solid substrates. Lithographic photopolymerization of a diacetylene-containing phospholipid, 1,2-bis(10,12-tricosadiynoyl)-sn-glycero-3-phosphocholine (DiynePC), and subsequent removal of nonreacted monomers by a detergent solution (0.1 M sodium dodecyl sulfate (SDS)) yielded a patterned polymeric bilayer matrix on the substrate. Fluid lipid bilayers of phosphatidylcholine from egg yolk (egg-PC) were incorporated into the lipid-free wells surrounded by the polymeric bilayers through the process of fusion and reorganization of suspended small unilamellar vesicles. Spatial distribution of the fluid bilayers in the patterned bilayer depended on the degree of photopolymerization that in turn could be modulated by varying the applied UV irradiation dose. The polymeric bilayer domains blocked lateral diffusion of the fluid lipid bilayers and confined them in the defined areas (corrals), if the polymerization was conducted with a sufficiently large UV dose. On the other hand, lipid molecules of the fluid bilayers penetrated into the polymeric bilayer domains, if the UV dose was relatively small. A direct correlation was observed between the applied UV dose and the lateral diffusion coefficient of fluorescent marker molecules in the fluid bilayers embedded within the polymeric bilayer domains. Artificial control of lateral diffusion by polymeric bilayers may lead to the creation of complex and versatile biomimetic model membrane arrays.     

Interaction of propyl paraben with dipalmitoyl phosphatidylcholine bilayer: a differential scanning calorimetry and nuclear magnetic resonance study

The influence of the preservative, propyl paraben (PPB) on the biophysical properties of dipalmitoyl phosphatidyl choline (DPPC) vesicles, both in multilamellar vesicle (MLV) and unilamellar vesicle (ULV) forms, has been studied using DSC and (¹H and ³¹P) NMR. The mechanism by which PPB interacts with DPPC bilayers was found to be independent of the morphological organization of the lipid bilayer. Incorporation of PPB in DPPC vesicles causes a significant depression in the transition temperature and enthalpy of both the pre-transition (PT) and the gel to liquid crystalline transition. The presence of the PPB also reduces the co-operativity of these transitions. However, at high PPB concentration the PT disappears. DSC and NMR findings indicate that: (i) PPB is bound strongly to the lipid bilayer leading to increased headgroup fluidity due to reduced headgroup–headgroup interaction and (ii) the PPB molecules are intercalated between the DPPC polar headgroups with its alkyl chain penetrate into the co-operative region. MLV incorporated with high PPB concentration shows additional transitions whose intensity increases with increasing PPB concentration. This phase segregation observed could probably be due to co-existence of PPB-rich and PPB-poor phospholipid domains within the bilayers. The effect of inclusion of cholesterol in the PPB-free and PPB-doped DPPC dispersion was also studied. Equilibration studies suggest that PPB molecules are very strongly bound and remain intercalated between the polar headgroup for prolonged time.     

Impact of cholesterol on voids in phospholipid membranes

Free volume pockets or voids are important to many biological processes in cell membranes. Free volume fluctuations are a prerequisite for diffusion of lipids and other macromolecules in lipid bilayers. Permeation of small solutes across a membrane, as well as diffusion of solutes in the membrane interior are further examples of phenomena where voids and their properties play a central role. Cholesterol has been suggested to change the structure and function of membranes by altering their free volume properties. We study the effect of cholesterol on the properties of voids in dipalmitoylphosphatidylcholine (DPPC) bilayers by means of atomistic molecular dynamics simulations. We find that an increasing cholesterol concentration reduces the total amount of free volume in a bilayer. The effect of cholesterol on individual voids is most prominent in the region where the steroid ring structures of cholesterol molecules are located. Here a growing cholesterol content reduces the number of voids, completely removing voids of the size of a cholesterol molecule. The voids also become more elongated. The broad orientational distribution of voids observed in pure DPPC is, with a 30% molar concentration of cholesterol, replaced by a distribution where orientation along the bilayer normal is favored. Our results suggest that instead of being uniformly distributed to the whole bilayer, these effects are localized to the close vicinity of cholesterol molecules.