Spin-labeled alkylphospholipids in a dipalmitoylphosphatidylcholine bilayer: molecular dynamics simulations

Molecular dynamics simulation has been performed to investigate the structural properties of perifosine and its synthetic spin-labeled alkylphospholipid analogues. The conformations adopted by these compounds in water and in a dipalmitoylphosphatidylcholine bilayer as a function of the presence and position of the N-oxyl-4',4'-dimethyloxazolidine ring (doxyl group) have been investigated by all-atom molecular dynamics. No predominant conformation was observed in water, but the molecules adopt specific orientations and conformations in the lipid bilayer. As is expected, alkyl chains tend to insert into the hydrophobic core, while charged groups stay at the lipid-water interface. A doxyl group in the middle of the alkyl chain moves up to the interface region, thus preventing adoption of the extended conformation. Compounds with a doxyl group close to the polar head group adopt conformations similar to that of unlabeled perifosine within the first nanoseconds of simulation. When the doxyl group is at the end of alkyl chain, the spin-labeled molecule needs more time to reach equilibrium. These results indicate a considerable effect of the doxyl position within the alkyl chain on its localization in the lipid bilayer and can be extended further to other similar spin probes used in the electron paramagnetic resonance spectroscopy of biological membranes.     

Structure of a gel phase lipid bilayer prepared by the Langmuir-Blodgett/Langmuir-Schaefer method characterized by sum-frequency vibrational spectroscopy

The structure of a planar supported lipid bilayer (PSLB) prepared by the Langmuir-Blodgett (LB)/Langmuir-Schaefer (LS) method was investigated by sum-frequency vibrational spectroscopy (SFVS). By using asymmetric lipid bilayers composed of selectively deuterated 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) lipids, the orientation of the fatty acid chains and phosphocholine headgroups has been determined independently for both leaflets of the bilayer. The alkyl chains of the lipids were found to be orientated approximately 13 degrees +/- 4 degrees from the surface normal for both leaflets. The lipid chains in both leaflets also contain some gauche content, which is consistent with previous NMR and FTIR studies of similar lipid systems. More importantly, the relative number of gauche defects does not seem to be influenced by the deposition method, LB versus LS. The headgroup orientation for the lipid film in contact with the silica support was determined to be 69 degrees +/- 3 degrees , whereas that in contact with the aqueous phase was 66 degrees +/- 4 degrees from the surface normal. The SFVS results indicate that the structure of the DSPC lipid film in contact with the solid support and the film adjacent to the aqueous phase are nearly identical in structure. These results suggesting the LB/LS deposition method do indeed produce symmetric lipid bilayers. These studies further add to the growing information on the efficacy of PSLBs as suitable models for biological membrane studies.     

3,4,5-Trimethoxybenzoate of Catechin,an Anticarcinogenic Semisynthetic Catechin,Modulates the Physical Properties of Anionic Phospholipid Membranes

3,4,5-Trimethoxybenzoate of catechin (TMBC) is a semisynthetic catechin which shows strong antiproliferative activity against malignant melanoma cells. The amphiphilic nature of the molecule suggests that the membrane could be a potential site of action, hence the study of its interaction with lipid bilayers is mandatory in order to gain information on the effect of the catechin on the membrane properties and dynamics. Anionic phospholipids, though being minor components of the membrane, possess singular physical and biochemical properties that make them physiologically essential. Utilizing phosphatidylserine biomimetic membranes, we study the interaction between the catechin and anionic bilayers, bringing together a variety of experimental techniques and molecular dynamics simulation. The experimental data suggest that the molecule is embedded into the phosphatidylserine bilayers, where it perturbs the thermotropic gel to liquid crystalline phase transition. In the gel phase, the catechin promotes the formation of interdigitation, and in the liquid crystalline phase, it decreases the bilayer thickness and increases the hydrogen bonding pattern of the interfacial region of the bilayer. The simulation data agree with the experimental ones and indicate that the molecule is located in the interior of the anionic bilayer as monomer and small clusters reaching the carbonyl region of the phospholipid, where it also disturbs the intermolecular hydrogen bonding between neighboring lipids. Our observations suggest that the catechin incorporates well into phosphatidylserine bilayers, where it produces structural changes that could affect the functioning of the membrane.     

Molecular characterization of gel and liquid-crystalline structures of fully hydrated POPC and POPE bilayers

Molecular dynamics simulations were used for a comprehensive study of the structural properties of monounsaturated POPC and POPE bilayers in the gel and liquid-crystalline state at a number of temperatures, ranging from 250 to 330 K. Though the chemical structures of POPC and POPE are largely similar (choline versus ethanolamine headgroup), their transformation processes from a gel to a liquid-crystalline state are contrasting. In the similarities, the lipid tails for both systems are tilted below the phase transition and become more random above the phase transition temperature. The average area per lipid and bilayer thickness were found less sensitive to phase transition changes as the unsaturated tails are able to buffer reordering of the bilayer structure, as observed from hysteresis loops in annealing simulations. For POPC, changes in the structural properties such as the lipid tail order parameter, hydrocarbon trans-gauche isomerization, lipid tail tilt-angle, and level of interdigitation identified a phase transition at about 270 K. For POPE, three temperature ranges were identified, in which the lower one (270-280 K) was associated with a pre-transition state and the higher (290-300 K) with the post-transition state. In the pre-transition state, there was a significant increase in the number of gauche arrangements formed along the lipid tails. Near the main transition (280-290 K), there was a lowering of the lipid order parameters and a disappearance of the tilted lipid arrangement. In the post-transition state, the carbon atoms along the lipid tails became less hindered as their density profiles showed uniform distributions. This study also demonstrates that atomistic simulations of current lipid force fields are capable of capturing the phase transition behavior of lipid bilayers, providing a rich set of molecular and structural information at and near the main transition state.     

Modification of Lipid Bilayer Structure by Diacylglycerol: A Comparative Study of Diacylglycerol and Cholesterol

Diacylglycerols (DAGs) are important second messengers in biomembranes, and they can activate protein kinase C and many other enzymes and receptors. However, their interactions with cholesterol and other lipids have not been previously studied using molecular dynamics (MD) simulation. In this study, nine independent atomistic MD simulations were performed to specifically investigate the interactions between di16:0DAG, 16:0,18:1-phosphatidylcholine (POPC), and cholesterol. Despite of their substantial differences in chemical structure, DAG and cholesterol produce some very similar effects in POPC bilayers: increasing acyl chain order and bilayer thickness, reducing volume-per-lipid, and decreasing lateral diffusion of molecules. More significantly, DAG also produces a strong "condensing effect" in PC bilayers. In comparison, cholesterol is more effective than DAG in producing the above effects. The driving force for the condensing effect is their molecular shape: DAG and cholesterol both have small polar headgroups and large hydrophobic bodies. In a lipid bilayer, in order to avoid the unfavorable exposure of their hydrophobic parts to water, neighboring phospholipid headgroups move toward cholesterol or DAG to provide cover. Thus, seemingly complex interactions between DAG, cholesterol and phospholipid can be clearly explained using the Umbrella Model. Our simulations confirmed the hypothesis that DAG increases the spacing between phospholipid headgroups, which is important for activating protein kinase C and other enzymes. Interestingly, our simulations also show that the conventional wisdom that the spacing created by a DAG is directly above the DAG molecule is incorrect; instead, the largest spacing usually occurs between the first and the second nearest-neighbor PC headgroups from a DAG, due to the umbrella effect.     

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.     

Freezing point depression of water in phospholipid membranes: a solid-state NMR study

Lipid-water interaction plays an important role in the properties of lipid bilayers, cryoprotectants, and membrane-associated peptides and proteins. The temperature at which water bound to lipid bilayers freezes is lower than that of free water. Here, we report a solid-state NMR investigation on the freezing point depression of water in phospholipid bilayers in the presence and absence of cholesterol. Deuterium NMR spectra at different temperatures ranging from -75 to + 10 degrees C were obtained from fully (2)H2O-hydrated POPC (1-palmitoyl-2-oleoylphosphatidylcholine) multilamellar vesicles (MLVs), prepared with and without cholesterol, to determine the freezing temperature of water and the effect of cholesterol on the freezing temperature of water in POPC bilayers. Our 2H NMR experiments reveal the motional behavior of unfrozen water molecules in POPC bilayers even at temperatures significantly below 0 degrees C and show that the presence of cholesterol further lowered the freezing temperature of water in POPC bilayers. These results suggest that in the presence of cholesterol the fluidity and dynamics of lipid bilayers can be retained even at very low temperatures as exist in the liquid crystalline phase of the lipid. Therefore, bilayer samples prepared with a cryoprotectant like cholesterol should enable the performance of multidimensional solid-state NMR experiments to investigate the structure, dynamics, and topology of membrane proteins at a very low temperature with enhanced sample stability and possibly a better sensitivity. Phosphorus-31 NMR data suggest that lipid bilayers can be aligned at low temperatures, while 15N NMR experiments demonstrate that such aligned samples can be used to enhance the signal-to-noise ratio of is 15N chemical shift spectra of a 37-residue human antimicrobial peptide, LL-37.     

Monte Carlo study of lipid membranes: Simulation of diparmitoylphosphatidylcholine bilayers in gel and liquid-crystalline phases

The statistical properties of the bilayer membranes of diparmitoylphosphatidylcholine (DPPC) in the gel and liquid-crystal phases were studied by Monte Carlo (MC) simulation using potential functions of the Lennard-Jones, the simple Coulomb, and the bond torsion. The simulation was undertaken on a two-dimensional periodic condition imposed on the bilayer model consisting of faithfully described molecules. The structure and ordering of the model bilayers accorded well with experiments, and the segment order parameters were in agreement with those of the nuclear magnetic resonance (NMR) experiments. The two kinds of lipid chains of DPPC do not equivalently behave in the bilayers, and chain 2 has lower ordering than chain 1. The order parameters of the first eight segments of chain 2 in the liquid-crystal model are particularly small and are roughly constant. From electron density analysis, it has been observed that the liquid-crystal bilayer has about one excess water molecule per one lipid molecule in comparison with the gel bilayer. The energy difference between the two bilayer models, taking account of the water contribution, is consistent with the latent heat of the phase transition. © 1995 by John Wiley & Sons, Inc.     

Molecular dynamics simulations of ternary membrane mixture: phosphatidylcholine, phosphatidic acid, and cholesterol

A ternary mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and cholesterol (CHOL) works effectively for a functional conformation of nicotinic acetylcholine receptor (nAChR) that can undergo agonist-induced conformation changes, but POPC alone can stabilize only a desensitized state of nAChR. To gain insights into the lipid mixture that has strong impact to nAChR functions, we performed more than 50 ns all atom molecular dynamic (MD) simulations at 303 K on a fully hydrated bilayer consisting of 240 POPC, 80 POPA, and 80 CHOL (3:1:1). The MD simulation revealed various interactions between different types of molecular pairs that ultimately regulated lipid organization. The heterogeneous interactions among three different constituents resulted in a broad spectrum of lipid properties, including extensive distributions of average area per lipid and varied lipid ordering as a function of lipid closeness to CHOL. Higher percentage of POPA than POPC had close association with CHOL, which coincided with relatively higher ordering of POPA molecules in their acyl chains near lipid head groups. Lower fraction of gauche dihedrals was also found in the same region of POPA. Although the CHOL molecules had the effects on the enhancement of surrounding lipid order, relatively low lipid order parameters and high fraction of gauche bonds were observed in the ternary mixture. Collectively, these results suggest that the dynamical structure of the ternary system could be determinant for a functional nAChR.     

Electrochemical and PM-IRRAS studies of the effect of cholesterol on the properties of the headgroup region of a DMPC bilayer supported at a Au(111) electrode

Polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS) was employed to investigate the interaction of cholesterol with the headgroups of dimyristoylphosphatidycholine (DMPC) molecules under a static electric field. DMPC/cholesterol (7:3 molar ratio) mixtures form a bilayer on a Au(111) electrode surface by fusion and spreading of small unilamellar vesicles. PM-IRRAS experiments provided detailed information concerning the conformation and hydration of headgroups of DMPC bilayers in the presence and absence of 30% cholesterol. The presence of 30% cholesterol increases the space between the headgroups of DMPC molecules and hence increases the hydration of the DMPC/cholesterol mixed bilayer. The conformational state of the headgroups of DMPC molecules in the mixed bilayer is also significantly changed. The phosphate group is closer to the surface compared with the pure DMPC bilayer. The conformation of the -O-C-C-N moiety changes from gauche to trans in the presence of cholesterol.     

Modulating membrane properties: the effect of trehalose and cholesterol on a phospholipid bilayer

The protective properties of trehalose on cholesterol-containing lipid dipalmitoylphosphatidylcholine (DPPC) bilayers are studied through molecular simulations. The ability of the disaccharide to interact with the phospholipid headgroups and stabilize the membrane persists even at high cholesterol concentrations and restricts some of the changes to the structure that would otherwise be imposed by cholesterol molecules. Predictions of bilayer properties such as area per lipid, tail ordering, and chain conformation support the notion that the disaccharide decreases the main melting transition in these multicomponent model membranes, which correspond more closely to common biological systems than pure bilayers. Molecular simulations indicate that the membrane dynamics are slowed considerably by the presence of trehalose, indicating that high sugar concentrations would serve to avert possible phase separations that could arise in mixed phospholipid systems. Various time correlation functions suggest that the character of the modifications in lipid dynamics induced by trehalose and cholesterol is different in the hydrophilic and hydrophobic regions of the membrane.     

Using spin polarised positive muons for studying guest molecule partitioning in soft matter structures

Fully polarised positive muons substituted for protons in organic free radicals can be used as spin labels which reveal information about the structure, dynamics and environment of these radicals. In applications via the technique of avoided-level-crossing muon spin resonance (ALC-microSR), the positive muon has been used to study the partitioning of phenyl alcohols in lamellar phase colloidal dispersions of a cationic dichain surfactant. Here we describe the experimental technique which permits highly sensitive spectroscopy as previously demonstrated for surfactant mixtures. We also demonstrate its capability in the study of partitioning of cosurfactant molecules in surfactant bilayers in order to elucidate the main factors which contribute to cosurfactant ordering at interfaces. The technique takes advantage of the positive muon combining with an electron to a hydrogen-like atom that is called muonium. This atom attaches to a phenyl group, forming a cyclohexadienyl-type radical that contains the muon as a polarised spin label, providing an excellent probe even for very low phenyl alcohol concentrations. The position of one type of resonance, which on the basis of spectroscopic selection rules is denoted as Delta(0), is related to the solvent polarity of the radicals' environment. The results derived from Delta(0) measurements reveal a systematic trend where the increasing chain length of the phenyl alcohol results in a deeper immersion of the phenyl ring of the alcohol into the surfactant bilayer with the OH group anchored at the interface. In addition, the data suggest partial penetration of water molecules into the bilayer. Furthermore, data ensuing from a second resonance (called Delta(1), which is dependent upon the degree of confinement of the radical within the surfactant aggregate structure) indicates not only that the phenyl alcohol resides in an anisotropic environment, (i.e. that the host molecule is unable to undergo full 3-D reorientation on a timescale of 50 ns), but the resonance line widths indicate that the radicals are undergoing fast rotation about a particular axis, in this instance about the first C-C substituent bond at the phenyl ring. Detailed analysis of these Delta(1) line shapes suggests that other types of motion involving reorientation of the above rotation axis are also present. At room temperature, the hydrocarbon chains of the double layers form an aggregate state commonly referred to as the L(beta) phase, where the motions of surfactant alkyl chains are effectively frozen out. These chains melt on heating over a temperature range which is solution composition dependent (ca. 51 to 67 degrees C), but in all cases leading to a liquid-like disordered hydrocarbon regime whilst retaining the overall lamellar structure (and in this state is termed L(alpha)). Above the L(alpha)/L(beta) chain ordering phase transition the tracer molecules reside within the bilayer, but below this transition (and depending on their water-oil solubility) they are completely or partly expelled. This interpretation is further supported by Heisenberg spin exchange experiments. The water-bilayer partitioning reflects both typical classical and nonclassical hydrophobic solvation depending on temperature and chain length of phenyl alcohols.     

Combined Monte Carlo and molecular dynamics simulation of hydrated 18:0 sphingomyelin-cholesterol lipid bilayers

We have carried out atomic level molecular dynamics and Monte Carlo simulations of hydrated 18:0 sphingomyelin (SM)-cholesterol (CHOL) bilayers at temperatures of 20 and 50 degrees C. The simulated systems each contained 266 SM, 122 CHOL, and 11861 water molecules. Each simulation was run for 10 ns under semi-isotropic pressure boundary conditions. The particle-mesh Ewald method was used for long-range electrostatic interactions. Properties of the systems were calculated over the final 3 ns. We compare the properties of 20 and 50 degrees C bilayer systems with each other, with experimental data, and with experimental and simulated properties of pure SM bilayers and dipalmitoyl phospatidyl choline (DPPC)-CHOL bilayers. The simulations reveal an overall similarity of both systems, despite the 30 degrees C temperature difference which brackets the pure SM main phase transition. The area per molecule, lipid chain order parameter profiles, atom distributions, and electron density profiles are all very similar for the two simulated systems. Consistent with simulations from our lab and others, we find strong intramolecular hydrogen bonding in SM molecules between the phosphate ester oxygen and the hydroxyl hydrogen atoms. We also find that cholesterol hydroxyl groups tend to form hydrogen bonds primarily with SM carbonyl, methyl, and amide moieties and to a lesser extent methyl and hydroxyl oxygens.     

Computational modeling of poly(alkylthiophene) conductive polymer insertion into phospholipid bilayers

We have previously demonstrated that some poly(alkylthiophenes) (PATs) are able to increase the electrical conductance of unsupported phospholipid bilayers and have hypothesized that this effect is due to the ability of some PAT side chains to permit stable insertion into the bilayer. We have further proposed the development of long-term intracellular electrodes based on that phenomenon. In this article, we apply molecular dynamics techniques to study the insertion of two model PATs into a patch of a lipid bilayer. Steered molecular dynamics is used to obtain potential trajectories of insertion, followed by umbrella sampling to determine the free-energy change upon insertion. Our results indicate that both branched-side-chain poly(3-(2-ethylhexyl)thiophene) (EHPT) and straight-side-chain poly(3-hexylthiophene) (HPT) are able to enter the bilayer but only EHPT can cross the center of the membrane and establish an electrical bridge. HPT penetrates the head groups but is not able to enter the alkyl tail phase. These findings support the feasibility of our electrode concept and raise questions regarding the mechanisms by which branched side chains grant PATs greater solubility in a lipid bilayer environment. The parameters and methods used in this study establish a novel framework for studying these and similar systems, and the results hold promise for the use of EHPT in biosensing and neural interfacing.     

Construction of higher-ordered monolayer membranes derived from archaeal membrane lipid-inspired cyclic lipids with longer alkyl chains

A series of artificial cyclic lipids that mimic archaeal membrane ones has been synthesized. The structural features of these molecules include a longer cyclic framework, in which the alkyl chain length ranges from 24 to 32 in carbon number, which is longer than our first analogous molecule with 20-carbon long alkyl chains [K. Miyawaki, T. Takagi, M. Shibakami, Synlett 8 (2002) 1326]. Microscopic observation reveals that these molecules have a self-assembling ability: hydration of the lipids yields multilamellar vesicles in aqueous solution and monolayer sheets on solid supports. High-sensitivity differential scanning calorimetry (24- and 28-carbon alkyl chain lipids) indicates that (i) the alkyl chain length affects their phase behavior and (ii) the enthalpies of endothermic peaks accompanied by phase transition were considerably lower than those of their monomeric phospholipid analogs. Fluorescence polarization measurements suggest that the membranes made from the 24-carbon alkyl chain lipid have a higher polarization factor than membranes composed of DMPC and DMPC plus cholesterol. These findings imply that the cyclic lipids containing 24- and 28-carbon alkyl chain construct well-organized monolayer membranes and, in particular, that the molecular order of the 24-carbon alkyl chain lipid is higher than that of bilayer membranes in the liquid-ordered phase.     

Heterogeneous molecular distribution in supported multicomponent lipid bilayers

Membrane domains contribute important structural and functional attributes to biological membranes. We describe the heterogeneous nanoscale distribution of lipid molecules within microscale membrane domains in multicomponent lipid bilayers composed of dipalmitoylphosphatidylcholine (DPPC), dilauroylphosphatidylcholine (DLPC), and cholesterol (chol). The lipids were labeled with the fluorescent lipid analogues Bodipy-PC and DiI-C20:0 to identify the distribution of individual membrane components. We used a near-field scanning optical microscope (NSOM) at room temperature to identify the nanoscale structures in the membrane. Simultaneous multicolor NSOM imaging at the emission maxima of the fluorescent analogues revealed a patchy distribution of Bodipy-PC and DiI-C20:0 indicative of phase separations in the bilayer. In a cholesterol-free system (DPPC/DLPC = 1:1), NSOM images proved that the two phosphatidylcholine molecules can coexist in domains at the micrometer level but form nanoscopic patches within the domains; DPPC occurs at the edge of the domains, whereas DLPC is present throughout the domains. In the presence of cholesterol (DPPC/DLPC = 7:3, chol = 18.9%), the two lipid molecules were more miscible but incomplete phase separations also occurred. The average domain sizes were 140-200 nm, well below the resolution capabilities of diffraction-limited light microscopy techniques; the domains were unresolvable by confocal microscopy. Our high-resolution NSOM studies of membrane domain behavior provide a better understanding of complex membrane phase phenomena in multicomponent biological membranes.     

Visualizing association of N-ras in lipid microdomains: influence of domain structure and interfacial adsorption

In this study, two-photon fluorescence microscopy on giant unilamellar vesicles and tapping-mode atomic force microscopy (AFM) are applied to follow the insertion of a fluorescently (4,4-difluoro-4-bora-3a,4a-diaza-s-indacene, BODIPY) labeled and completely lipidated (hexadecylated and farnesylated) N-Ras protein into heterogeneous lipid bilayer systems. The bilayers consist of the canonical raft mixture 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), sphingomyelin, and cholesterol, which-depending on the concentration of the constituents-separates into liquid-disordered (l(d)), liquid-ordered (l(o)), and solid-ordered (s(o)) phases. The results provide direct evidence that partitioning of N-Ras occurs preferentially into liquid-disordered lipid domains, which is also reflected in a faster kinetics of incorporation into the fluid lipid bilayers. The phase sequence of preferential binding of N-Ras to mixed-domain lipid vesicles is l(d) > l(o) > s(o). Intriguingly, we detect, using the better spatial resolution of AFM, also a large proportion of the lipidated protein located at the l(d)/l(o) phase boundary, thus leading to a favorable decrease in line tension that is associated with the rim of the demixed phases. Such an interfacial adsorption effect may serve as an alternative vehicle for association processes of signaling proteins in membranes.     

Comparative molecular dynamics study of ether- and ester-linked phospholipid bilayers

The lipid membranes found in archaea have high bilayer stability and low permeability. The molecular structure of their constituent lipids is characterized by ether-linked, branched hydrophobic chains, whereas the conventional lipids obtained from eukaryotic or eubacterial sources have ester linked straight chains. In order to elucidate the influence of the ether linkage, instead of an ester one, on the physical properties of the lipid bilayers, we have carried out comparative 10 ns molecular dynamics simulations of diphytanyl phosphatidylcholine (ether-DPhPC) and diphytanoyl phosphatidylcholine (ester-DPhPC) bilayers in water, respectively. We analyze bilayer structures, hydration of the lipids, membrane dipole potentials, and free energy profiles of water and oxygen across the bilayers. We observe that the membrane dipole potential for the ether-DPhPC bilayer, which arises mainly from the ether linkage, is about half of that of the ester-DPhPC. The calculated free energy barrier for a water molecule in the ether-DPhPC bilayer system is slightly higher than that in the ester-DPhPC counterpart, which is in accord with experimental data.     

Interleaflet interaction and asymmetry in phase separated lipid bilayers: molecular dynamics simulations

In order to investigate experimentally inaccessible, molecular-level detail regarding interleaflet interaction in membranes, we have run an extensive series of coarse-grained molecular dynamics simulations of phase separated lipid bilayers. The simulations are motivated by differences in lipid and cholesterol composition in the inner and outer leaflets of biological membranes. Over the past several years, this phenomenon has inspired a series of experiments in model membrane systems which have explored the effects of lipid compositional asymmetry in the two leaflets. The simulations are directed at understanding one potential consequence of compositional asymmetry, that being regions of bilayers where liquid-ordered (L(o)) domains in one leaflet are opposite liquid-disordered (L(d)) domains in the other leaflet (phase asymmetry). The simulated bilayers are of two sorts: 1) Compositionally symmetric leaflets where each of the two leaflets contains an identical, phase separated (L(o)/L(d)) mixture of cholesterol, saturated and unsaturated phospholipid; and 2) Compositionally asymmetric leaflets, where one leaflet contains a phase separated (L(o)/L(d)) mixture while the other contains only unsaturated lipid, which on its own would be in the L(d) phase. In addition, we have run simulations where the lengths of the saturated lipid chains as well as the mole ratios of the three lipid components are varied. Collectively, we report on three types of interleaflet coupling within a bilayer. First, we show the effects of compositional asymmetry on acyl chain tilt and order, lipid rotational dynamics, and lateral diffusion in regions of leaflets that are opposite L(o) domains. Second, we show substantial effects of compositional asymmetry on local bilayer curvature, with the conclusion that phase separated leaflets resist curvature, while inducing large degrees of curvature in an opposing L(d) leaflet. Finally, in compositionally symmetric, phase separated bilayers, we find phase asymmetry (domain antiregistration) between the two leaflets occurs as a consequence of mismatched acyl chain-lengths in the saturated and unsaturated lipids.     

Molecular dynamics simulations of lipid bilayers

Computer simulation methods are becoming increasingly widespread as tools for studying the structure and dynamics of lipid bilayer membranes. The length scale and time scale accessible to atomic-level molecular dynamics simulations are rapidly increasing, providing insight into the relatively slow motions of molecular reorientation and translation and demonstrating that effects due to the finite size of the simulation cell can influence simulation results. Additionally, significant advances have been made in the complexity of membrane systems studied, including bilayers with cholesterol, small solute molecules, and lipid-protein and lipid-DNA complexes. Especially promising is the progress that continues to be made in the comparison of simulation results with experiment, both to validate the simulation algorithms and to aid in the interpretation of existing experimental data.