Effects of the Lipid Bilayer Phase State on the Water Membrane Interface

Atomic-scale structure of interfacial water on gel and liquid phase lipid membranes

Hydration of biological membranes is essential to a wide range of biological processes. In particular, it is intrinsically linked to lipid thermodynamic properties, which in turn influence key cell functions such as ion permeation and protein mobility. Experimental and theoretical studies of the surface of biomembranes have revealed the presence of an interfacial repulsive force, which has been linked to hydration or steric effects. Here, we directly characterise the atomic-scale structure of water near supported lipid membranes of 1,2-dimyristoyl-sn-glycero-3-phosphocholine in their gel and liquid phase through three-dimensional atomic force microscopy (3D AFM). First, we demonstrate the ability to probe the morphology of interfacial water of lipid bilayers in both phases with sub-molecular resolution by using ultrasharp tips. We then visualise the molecular arrangement of water at the lipid surface at different temperatures. Our experiments reveal that water is organised in multiple hydration layers on both the solid-ordered and liquid-disordered lipid phases. Furthermore, we observe a monotonic repulsive force, which becomes relevant only in the liquid phase. These results offer new insights into the water structuring near soft biological surfaces, and demonstrate the importance of investigating it with vertical and lateral sub-molecular resolution.

l-Phenylalanine Partitioning Mechanisms in Model Biological Membranes

Time-resolved fluorescence spectroscopy in combination with differential scanning calorimetry (DSC) was used to study the chemical interactions that occur when l-phenylalanine is introduced to solutions containing phosphatidylcholine vesicles. Studies reported in this work address open questions about l-Phe's affinity for lipid vesicle bilayers, the effects of l-Phe partitioning on bilayer properties, l-Phe's solvation within a lipid bilayer, and the amount of l-Phe within that local solvation environment. DSC data show that l-Phe reduces the amount of heat necessary to melt saturated phosphatidylcholine bilayers from their gel to liquid-crystalline state but does not change the transition temperature (T_gel-lc). Time-resolved emission shows only a single l-Phe lifetime at low temperatures corresponding to l-Phe remaining solvated in aqueous solution. At temperatures close to T_gel-lc, a second, shorter lifetime appears that is assigned to l-Phe already embedded within the membrane that becomes hydrated as water starts to permeate the lipid bilayer. This new lifetime is attributed to a conformationally restricted rotamer in the bilayer's polar headgroup region and accounts for up to 30% of the emission amplitude. Results reported for dipalmitoylphosphatidylcholine (DPPC, 16:0) lipid vesicles prove to be general, with similar effects observed for dimyristoylphosphatidylcholine (DMPC, 14:0) and distearoylphosphatidylcholine (DSPC, 18:0) vesicles. Taken together, these results create a complete and compelling picture of how l-Phe associates with model biological membranes. Furthermore, this approach to examining amino acid partitioning into membranes and the resulting solvation forces points to new strategies for studying the structure and chemistry of membrane-soluble peptides and selected membrane proteins.

Cooperativity between sodium ions and water molecules facilitates lipid mobility in model cell membranes

Cellular membranes are surrounded by an aqueous buffer solution containing various ions, which influence the hydration layer of the lipid head groups. At the same time, water molecules hydrating the lipids play a major role in facilitating the organisation and dynamics of membrane lipids. Employing fluorescence microscopy imaging and fluorescence recovery after photobleaching measurements, we demonstrate that the cooperativity between water and sodium (Na+) ions is crucial to maintain lipid mobility upon the removal of the outer hydration layer of the lipid membrane. Under similar hydration conditions, lipid diffusion ceases in the absence of Na+ ions. We find that Na+ ions (and similarly K+ ions) strengthen the water clathrate cage around the lipid phosphocholine headgroup and thus prevent its breaking upon removal of bulk water. Intriguingly, Ca2+ and Mg2+ do not show this effect. In this article, we provide a detailed molecular-level picture of ion specific dependence of lipid mobility and membrane hydration properties.

Effect of Cholesterol Versus Ergosterol on DPPC Bilayer Properties: Insights from Atomistic Simulations

Sterols have been ascribed a major role in the organization of biological membranes, in particular for the formation of liquid ordered domains in complex lipid mixtures. Here, we employed molecular dynamics simulations to compare the effects of cholesterol and ergosterol as the major sterol of mammalian and fungal cells, respectively, on binary mixtures with 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) as a proxy for saturated lipids. In agreement with previous work, we observe that the addition of sterol molecules modifies the order of DPPC both in the gel phase and in the liquid phase. When disentangling the overall tilt angle and the structure of the tail imposed by trans/gauche configurations of torsion angles in the tail, respectively, a more detailed picture of the impact of sterols can be formulated, revealing, for example, an approximate temperature-concentration superposition ranging from the liquid to the gel phase. Furthermore, a new quantitative measure to identify the presence of collective sterol effects is discussed. Moreover, when comparing both types of sterols, addition of cholesterol has a noticeably stronger impact on phospholipid properties than that of ergosterol. The observed differences can be attributed to higher planarity of the cholesterol ring system. This planarity combined with an inherent asymmetry in its molecular interactions leads to better alignment and hence stronger interaction with saturated acyl chains. Our results suggest that the high order demonstrated for ergosterol in fungal plasma membranes must therefore be generated via additional mechanisms.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE–CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Ion-induced modification of the sucrose network and its impact on melting of freeze-dried liposomes. DSC and molecular dynamics study

Disaccharides play an important role in survival of anhydrobiotic organisms during extreme environmental conditions. A key protection feature is their capability to form the hydrogen bond (HB) network in a similar fashion as the one made by water. Since various ions also affect the HB network in completely hydrated systems, it is of a great interest to understand how they impact preservation when incorporated in a disaccharide network. To address this, we employ a combination of experimental and modeling techniques to study behavior of multilamellar 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes freeze-dried with sucrose in presence of NaCl or NaH₂PO₄·H₂O at various concentrations (0.01-1M). Differential scanning calorimetry (DSC) was employed in order to determine the cooperative unit size (CUS), the number of lipid molecules that constitute a domain of cooperative motion in the liposome, and the melting temperature (T_m). In the absence of salt CUS was estimated to be 122±12, whereas in the presence of NaCl CUS increases more (347±34 for c=1M) than for NaH₂PO₄·H₂O (193±26 for 1M). When it comes to T_m, the situation is reversed; NaCl induces increase by about 1K, while NaH₂PO₄·H₂O by about 10K. These findings clearly demonstrate how different interaction forces-hydrogen bonding, charge pairing, and van der Waals interactions between acyl chains-affect CUS and T_m. Their interplay and contribution of particular interaction was further analyzed with molecular dynamics (MD) simulations. This analysis demonstrated that the HB network of DMPC and sucrose is partially disrupted in the presence of NaCl ions, and even to a greater extent in the case of NaH₂PO₄·H₂O ions. Notably, H₂PO₄^- ions outcompete and replace the sucrose molecules at the DMPC surface, which in turn alters the nature of the DMPC-surrounding interactions, from a weaker HB-dominated to a stronger CP-dominated interaction network.

Investigating the Structure of Multicomponent Gel-Phase Lipid Bilayers

Single- and multicomponent lipid bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphatidylcholine (DSPC), isostearyl isostearate, and heptadecanoyl heptadecanoate in the gel phase are studied via molecular dynamics simulations. It is shown that the structural properties of multicomponent bilayers can deviate strongly from the structures of their single-component counterparts. Specifically, the lipid mixtures are shown to adopt a compact packing by offsetting the positioning depths at which different lipid species are located in the bilayer. This packing mechanism affects the area per lipid, the bilayer height, and the chain tilt angles and has important consequences for other bilayer properties, such as interfacial hydrogen bonding and bilayer permeability. In particular, the simulations suggest that bilayers containing isostearyl isostearate or heptadecanoyl heptadecanoate are less permeable than pure 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine or DSPC bilayers. Furthermore, hydrogen-bond analysis shows that the residence times of lipid-water hydrogen bonds depend strongly on the bilayer composition, with longer residence times for bilayers that have a higher DSPC content. The findings illustrate and explain the fundamental differences between the properties of single- and multicomponent bilayers.

Calcium Assists Dopamine Release by Preventing Aggregation on the Inner Leaflet of Presynaptic Vesicles

In this study, the dopamine-lipid bilayer interactions were probed with three physiologically relevant ion compositions using atomistic molecular dynamics simulations and free energy calculations. The in silico results indicate that calcium is able to decrease significantly the binding of dopamine to a neutral (zwitterionic) phosphatidylcholine lipid bilayer model mimicking the inner leaflet of a presynaptic vesicle. We argue that the observed calcium-induced effect is likely in crucial role in the neurotransmitter release from the presynaptic vesicles docked in the active zone of nerve terminals. The inner leaflets of presynaptic vesicles, which are responsible for releasing neurotransmitters into the synaptic cleft, are mainly composed of neutral lipids such as phosphatidylcholine and phosphatidylethanolamine. The neutrality of the lipid head group region, enhanced by a low pH level, should limit membrane aggregation of transmitters. In addition, the simulations suggest that the high calcium levels inside presynaptic vesicles prevent even the most lipophilic transmitters such as dopamine from adhering to the inner leaflet surface, thus rendering unhindered neurotransmitter release feasible.

A computational study suggests that replacing PEG with PMOZ may increase exposure of hydrophobic targeting moiety

In a previous study we showed that the cause of failure of a new, proposed, targeting ligand, the AETP moiety, when attached to a PEGylated liposome, was occlusion by the poly(ethylene glycol) (PEG) layer due to its hydrophobic nature, given that PEG is not entirely hydrophilic. At the time we proposed that possible replacement with a more hydrophilic protective polymer could alleviate this problem. In this study we have used computational molecular dynamics modelling, using a model with all atom resolution, to suggest that a specific alternative protective polymer, poly(2-methyloxazoline) (PMOZ), would perform exactly this function. Our results show that when PEG is replaced by PMOZ the relative exposure to the solvent of AETP is increased to a level even greater than that we found in previous simulations for the RGD peptide, a targeting moiety that has previously been used successfully in PEGylated liposome based therapies. While the AETP moiety itself is no longer under consideration, the results of this computational study have broader significance: the use of PMOZ as an alternative polymer coating to PEG could be efficacious in the context of more hydrophobic targeting ligands. In addition to PMOZ we studied another polyoxazoline, poly(2-ethyloxazoline) (PEOZ), that has also been mooted as a possible alternate protective polymer. It was also found that the RDG peptide occlusion was significantly greater for the case of both oxazolines as opposed to PEG and that, unlike PEG, neither oxazoline entered the membrane. As far as we are aware this is the first time that polyoxazolines have been studied using molecular dynamics simulation with all atom resolution.

How does ytterbium chloride interact with DMPC bilayers? A computational and experimental study

Lanthanide salts have been studied for many years, primarily in Nuclear Magnetic Resonance (NMR) experiments of mixed lipid–protein systems and more recently to study lipid flip-flop in model membrane systems.

Water populations in restricted environments of lipid membrane interphases

We employ molecular dynamics simulations to study the hydration properties of Dipalmitoylphosphatidylcholine (DPPC) bilayers, both in the gel and the liquid crystalline states. We show that while the tight hydration centers (PO and CO moieties) are significantly hydrated in both phases, the gel-fluid transition involves significant changes at the second hydration shell, particularly at the buried region between the hydrocarbon tails. Thus, while almost no buried water population exists in the gel state below the carbonyls, this hydrophobic region becomes partially water accesible in the liquid crystalline state. We shall also show that such water molecules present a lower H-bond coordination as compared to the molecules at the primary hydration shell. This means that, while the latter are arranged in relatively compact nanoclusters (as already proposed), the buried water molecules tend to organize themselves in less compact structures, typically strings or branched strings, with a scarce population of isolated molecules. This behavior is similar to that observed in other hydration contexts, like water penetrating carbon nanotubes or model hydrophobic channels or pores, and reflects the reluctance of water to sacrifice HB coordination.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Mixed Mechanism of Lubrication by Lipid Bilayer Stacks

Although the key role of lipid bilayer stacks in biological lubrication is generally accepted, the mechanisms underlying their extreme efficiency remain elusive. In this article, we report molecular dynamics simulations of lipid bilayer stacks undergoing load and shear. When the hydration level is reduced, the velocity accommodation mechanism changes from viscous shear in hydration water to interlayer sliding in the bilayers. This enables stacks of hydrated lipid bilayers to act as efficient boundary lubricants for various hydration conditions, structures, and mechanical loads. We also propose an estimation for the friction coefficient; thanks to the strong hydration forces between lipid bilayers, the high local viscosity is not in contradiction with low friction coefficients.

Long-Range Lipid-Water Interaction as Observed by ATR-FTIR Spectroscopy

It is commonly assumed that the structure of water at a lipid-water interface is influenced mostly in the first hydration layer. However, recent results from different experimental methods show that perturbation extends through several hydration layers. Due to its low light penetration depth, attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy is specifically suited to study interlamellar water structure in multibilayers. Results obtained by this technique confirm the long-range water structure disturbance. Consequently, in confined membrane environments nearly all water molecules can be perturbed. It is important to note that the behavior of confined water molecules differs significantly in samples prepared in excess water and in partially hydrated samples. We show in what manner the interlamellar water perturbation is influenced by the hydration level and how it is sequentially modified with a step-by-step dehydration of samples either by water evaporation or by osmotic pressure. Our results also indicate that besides different levels of hydration the lipid-water interaction is modulated by different lipid headgroups and different lipid phases as well. Therefore, modification of interlamellar water properties may clarify the role of water-mediated effects in biological processes.

Role of hydration in phosphatidylcholine reverse micelle structure and gelation in cyclohexane: a molecular dynamics study

Molecular simulations reveal lipid headgroup hydration is at key role in organogel transition in lecithin–water–cyclohexane systems.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. III. Water clustering and vibrational energy transfer

Water clustering and connectivity around lipid bilayers strongly influences the properties of membranes and is important for functions such as proton and ion transport. Vibrational anisotropic pump-probe spectroscopy is a powerful tool for understanding such clustering, as the measured anisotropy depends upon the time-scale and degree of intra- and intermolecular vibrational energy transfer. In this article, we use molecular dynamics simulations and theoretical vibrational spectroscopy to help interpret recent experimental measurements of the anisotropy of water in lipid multi-bilayers as a function of both lipid hydration level and isotopic substitution. Our calculations are in satisfactory agreement with the experiments of Piatkowski, Heij, and Bakker, and from our simulations we can directly probe water clustering and connectivity. We find that at low hydration levels, many water molecules are in fact isolated, although up to 70% of hydration water forms small water clusters or chains. At intermediate hydration levels, water forms a wide range of cluster sizes, while at higher hydration levels, the majority of water molecules are part of a large, percolating water cluster. Therefore, the size, number, and nature of water clusters are strongly dependent on lipid hydration level, and the measured anisotropy reflects this through its dependence on intermolecular energy transfer.

Cholesterol level affects surface charge of lipid membranes in saline solution

Cholesterol is an important component of all biological membranes as well as drug delivery liposomes. We show here that increasing the level of cholesterol in a phospholipid membrane decreases surface charge in the physiological environment. Through molecular dynamics simulation we have shown that increasing the level of cholesterol decreases Na+ ion binding. Complementary experimental ζ--potential measurements have shown a decreased ζ--potential with increasing cholesterol content, indicative of reduced surface charge. Both experiments and simulations have been carried out on both saturated 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and monounsaturated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes. This result is particularly important because membrane surface charge plays an important role in the interactions of biomembranes with peripheral membrane proteins and drug delivery liposomes with the immune system.

Pairing of cholesterol with oxidized phospholipid species in lipid bilayers

We claim that (1) cholesterol protects bilayers from disruption caused by lipid oxidation by sequestering conical shaped oxidized lipid species such as 1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine (PZPC) away from phospholipid, because cholesterol and the oxidized lipid have complementary shapes and (2) mixtures of cholesterol and oxidized lipids can self-assemble into bilayers much like lysolipid–cholesterol mixtures. The evidence for bilayer protection comes from molecular dynamics (MD) simulations and dynamic light scattering (DLS) measurements. Unimodal size distributions of extruded vesicles (LUVETs) made up of a mixture of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and PZPC containing high amounts of PZPC are only obtained when cholesterol is present in high concentrations. In simulations, bilayers containing high amounts of PZPC become porous, unless cholesterol is also present. The protective effect of cholesterol on oxidized lipids has been observed previously using electron paramagnetic resonance (EPR) and electron microscopy imaging of vesicles. The evidence for the pairing of cholesterol and PZPC comes mainly from correlated 2-D density and thickness plots from simulations, which show that these two molecules co-localize in bilayers. Further evidence that the two molecules can cohabitate comes from self-assembly simulations, where we show that cholesterol-oxidized lipid mixtures can form lamellar phases at specific concentrations, reminiscent of lysolipid–cholesterol mixtures. The additivity of the packing parameters of cholesterol and PZPC explains their cohabitation in a planar bilayer. Oxidized lipids are ubiquitously present in significant amounts in high- and low-density lipoprotein (HDL and LDL) particles, diseased tissues, and in model phospholipid mixtures containing polyunsaturated lipids. Therefore, our hypothesis has important consequences for cellular cholesterol trafficking; diseases related to oxidized lipids, and to biophysical studies of phase behaviour of cholesterol-containing phospholipid mixtures.

Molecular dynamics simulations of the bacterial ABC transporter SAV1866 in the closed form

The ATP binding cassette (ABC) transporter family of proteins contains members involved in ATP-mediated import or export of ligands at the cell membrane. For the case of exporters, the translocation mechanism involves a large-scale conformational change that involves a clothespin-like motion from an inward-facing open state, able to bind ligands and adenosine triphosphate (ATP), to an outward-facing closed state. Our work focuses on SAV1866, a bacterial member of the ABC transporter family for which the structure is known for the closed state. To evaluate the ability of this protein to undergo conformational changes at physiological temperature, we first performed conventional molecular dynamics (MD) on the cocrystallized adenosine diphosphate (ADP)-bound structure and on a nucleotide-free structure. With this assessment of SAV1866's stability, conformational changes were induced by steered molecular dynamics (SMD), in which the nucleotide binding domains (NBD) were pushed apart, simulating the ATP hydrolysis energy expenditure. We found that the transmembrane domain is not easily perturbed by large-scale motions of the NBDs.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. II. Two-dimensional infrared and peak shift observables within different theoretical approximations

In a previous report, we calculated the infrared absorption spectrum and both the isotropic and anisotropic pump-probe signals for the OD stretch of isotopically dilute water in dilauroylphosphatidylcholine (DLPC) multi-bilayers as a function of the lipid hydration level. These results were then compared to recent experimental measurements and are in generally good agreement. In this paper, we will further investigate the structure and dynamics of hydration water using molecular dynamics simulations and calculations of the two-dimensional infrared and vibrational echo peak shift observables for hydration water in DLPC membranes. These observables have not yet been measured experimentally, but future comparisons may provide insight into spectral diffusion processes and hydration water heterogeneity. We find that at low hydration levels the motion of water molecules inside the lipid membrane is significantly arrested, resulting in very slow spectral diffusion. At higher hydration levels, spectral diffusion is more rapid, but still slower than in bulk water. We also investigate the effects of several common approximations on the calculation of spectroscopic observables by computing these observables within multiple levels of theory. The impact of these approximations on the resulting spectra affects our interpretation of these measurements and reveals that, for example, the cumulant approximation, which may be valid for certain systems, is not a good approximation for a highly heterogeneous environment such as hydration water in lipid multi-bilayers.

The ELBA force field for coarse-grain modeling of lipid membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (ε(r) = 1). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities.

Ultrafast energy migration pathways in self-assembled phospholipids interacting with confined water

Phospholipids self-assembled into reverse micelles in benzene are introduced as a new model system to study elementary processes relevant for energy transport in hydrated biological membranes. Femtosecond vibrational spectroscopy gives insight into the dynamics of the antisymmetric phosphate stretching vibration ν(AS)(PO(2))(-), a sensitive probe of local phosphate-water interactions and energy transport. The decay of the ν(AS)(PO(2))(-) mode with a 300-fs lifetime transfers excess energy to a subgroup of phospholipid low-frequency modes, followed by redistribution among phospholipid vibrations within a few picoseconds. The latter relaxation is accelerated by adding a confined water pool, an efficient heat sink in which the excess energy induces weakening or breaking of water-water and water-phospholipid hydrogen bonds. In parallel to vibrational relaxation, resonant energy transfer between ν(AS)(PO(2))(-) oscillators delocalizes the initial excitation.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. I. Infrared spectra and ultrafast pump-probe observables

The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.