Sensitivity of hydrogen bond lifetime dynamics to the presence of ethanol at the interface of a phospholipid bilayer

Atomistic molecular dynamics simulations of a fully hydrated liquid crystalline lamellar phase of a dimyrystoylphosphatidylcholine lipid bilayer containing ethanol at 1:1 composition as well as of the pure lamellar phase of the bilayer have been performed. Detailed analyses have been carried out to investigate the effects of ethanol, if any, on the lifetime dynamics of lipid-water and water-water hydrogen bonds in the hydration layer of the lipid headgroups. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed the identification of the bound states at the bilayer interface and the quantification of the dynamic equilibrium between the bound and the free water molecules, in terms of time-dependent relaxation rates. The calculations show that the overall relaxation of phosphate-water hydrogen bonds is faster in the presence of ethanol. Studies of the residence time and the number fluctuation of the hydration layer water molecules reveal that the presence of ethanol molecules decreases the rigidity of the lipid hydration layer.     

Inverse-phosphocholine lipids: a remix of a common phospholipid

Zwitterionic inverse-phosphocholine (iPC) lipids contain headgroups with an inverted charge orientation relative to phosphocholine (PC) lipids. The iPC lipid headgroup has a quaternary amine adjacent to the bilayer interface and a phosphate that extends into the aqueous phase. Neutral iPC lipids with ethylated phosphate groups (CPe) and anionic iPC lipids nonethylated phosphate groups (CP) were synthesized. The surface potential of CPe liposomes remains negative across a broad pH range and in the presence of up to 10 mM Ca(2+). CP liposomes aggregate in the presence of Ca(2+), but at a slower rate than other anionic lipids. Hydrolysis of CP lipids by alkaline phosphatases generates a cationic lipid. CPe liposomes release encapsulated anionic carboxyfluorescein (CF) 20 times faster than PC liposomes and release uncharged glucose twice as fast as PC liposomes. As such, iPC lipids afford a unique opportunity to investigate the biophysical and bioactivity-related ramifications of a charge inversion at the bilayer surface.     

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.     

Hydration layer of a cationic micelle, C(10)TAB: structure, rigidity, slow reorientation, hydrogen bond lifetime, and solvation dynamics

We report a theoretical study of the structure and dynamics of the water layer (the hydration layer) present at the surface of the cationic micelle decyltrimethylammonium bromide (DeTAB) by using atomistic molecular dynamics simulations. The simulated micelle consisted of 47 surfactant molecules (and an equal number of bromide ions), in good agreement with the pioneering light scattering experiments by Debye which found an aggregation number of 50. In this micelle, three partially positively charged methyl groups of each surfactant headgroup face the surrounding water. The nature of the cationic micellar surface is found to play an important role in determining the arrangement of water which is quite different from that in the bulk or on the surface of an anionic micelle, like cesium perfluorooctanoate. Water molecules present in the hydration layer are found to be preferentially distributed in the region between the three partially charged methyl headgroups. It is found that both the translational and rotational motions of water exhibit appreciably slower dynamics in the layer than those in the bulk. The solvation time correlation function (TCF) of bromide ions exhibits a long time component which is found to originate primarily from the interaction of the probe with the micellar headgroups. Thus, the decay of the solvation TCF is controlled largely by the residence time of the probe in the surface. The residence time distribution of the water molecules also exhibits a slow time component. We also calculate the collective number density fluctuation in the layer and find a prominent slow component compared to the similar quantity in the bulk. This slow component demonstrates that water structure in the hydration layer is more rigid than that in the bulk. These results demonstrate that the slow dynamics of hydration layer water is generic to macromolecular surfaces of either polarity.     

Lipids and lipid mixtures in boundary layers: From hydration lubrication to osteoarthritis

The hydration layer surrounding the phosphocholine headgroups of single-component phosphatidylcholine lipids, or of lipid-mixtures, assembled at an interface greatly modifies the interfacial properties and interactions. As water molecules within the hydration layer are held tightly by the headgroup but are nonetheless very fluid on shear, the boundary lipid layers, exposing the highly hydrated headgroup arrays, can provide efficient boundary lubrication when sliding against an opposing surface, at physiologically high contact pressures. In addition, any free lipids in the surrounding liquid can heal defects which may form during sliding on the boundary phosphatidylcholine layer. Similar boundary lipid layers contribute to the lubricating, pressure-bearing, and wear-protection functions of healthy articular joints. This review presents a survey of the relationship between the molecular composition of the interfacial complex and the lubrication behavior of the lipid-based boundary layers, which could be beneficial for designing boundary lubricants for intra-articular injection for the treatment of early osteoarthritis.     

Molecular dynamics study of surfactant monolayers adsorbed at the oil/water and air/water interfaces

Atomistic molecular dynamics (MD) simulations have been carried out to investigate the physical properties of monolayers of monododecyl diethylene glycol (C(12)E(2)) surfactants adsorbed at the oil/water and air/water interfaces. The study shows that the surfactant molecules exhibit more extended conformations with a consequent increase of the thickness of the monolayer in the presence of the oil medium. It is noticed that the hydrocarbon tails of the surfactants are more vertically oriented at the oil/water interface. Interestingly, we notice that the presence of the oil medium has a strong influence in restricting both the translational and reorientational motions of the water molecules present in the hydration layer close to the surfactant headgroups.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Dynamics and state of lipid bilayer-internal water unraveled with solution state 1H dynamic nuclear polarization

The dynamics and state of lipid bilayer-internal hydration water of unilamellar lipid vesicles dispersed in solutions is characterized. This study was enabled by a recently developed technique based on Overhauser dynamic nuclear polarization (DNP)-driven amplification of (1)H nuclear magnetic resonance (NMR) signal of hydration water. This technique can, in the full presence of bulk water, selectively quantify the translational dynamics of hydration water within ～10 ? around spin labels that are specifically introduced to the local volume of interest within the lipid bilayer. With this approach, the local apparent diffusion coefficients of internal water at different depths of the lipid bilayer were determined. The modulation of these values as a response to external stimuli, such as the addition of sodium chloride or ethanol and the lipid phase transitions, that alter the fluctuations of bilayer interfaces together with the activation energy values of water diffusivity shows that water is not individually and homogeneously solvating lipid's hydrocarbon tails in the lipid bilayer. We provide experimental evidence that instead, water and the lipid membrane comprise a heterogeneous system whose constituents include transient hydrophobic water pores or water structures traversing the lipid bilayer. We show how these transient pore structures, as key vehicles for passive water transport can better reconcile our experimental data with existing literature data on lipid bilayer hydration and dynamics.     

Ultrafast hydration dynamics in the lipidic cubic phase: Discrete water structures in nanochannels

We report here our studies of hydration dynamics of confined water in aqueous nanochannels (approximately 50 A) of the lipidic cubic phase. By systematically anchoring the hydrocarbon tails of a series of tryptophan-alkyl ester probes into the lipid bilayer, we mapped out with femtosecond resolution the profile of water motions across the nanochannel. Three distinct time scales were observed, revealing discrete channel water structures. The interfacial water at the lipid surface is well-ordered, and the relaxation dynamics occurs in approximately 100-150 ps. These dynamically rigid water molecules are crucial for global structural stability of lipid bilayers and for stabilization of anchored biomolecules in membranes. The adjacent water layers near the lipid interface are hydrogen-bonded networks and the dynamical relaxation takes 10-15 ps. This quasi-bound water motion, similar to the typical protein surface hydration relaxation, facilitates conformation flexibility for biological recognition and function. The water near the channel center is bulklike, and the dynamics is ultrafast in less than 1 ps. These water molecules freely transport biomolecules near the channel center. The corresponding orientational relaxation at these three typical locations is well correlated with the hydration dynamics and local dynamic rigidity. These results reveal unique water structures and dynamical motions in nanoconfinements, which is critical to the understanding of nanoscopic biological activities and nanomaterial properties.     

Structure and dynamics of water at the interface with phospholipid bilayers

We have performed two molecular-dynamics simulations to study the structural and dynamical properties of water at the interface with phospholipid bilayers. In one of the simulations the bilayer contained neutral phospholipid molecules, dioleoylphosphatidylcholine (DOPC); in the second simulation the bilayer contained charged lipid molecules, dioleoylphosphatidylserine (DOPS). From the density profile of water we observe that water next to the DOPS bilayer is more perturbed as compared to water near the DOPC bilayer. Using an energetic criterion for the determination of hydrogen bonding we find that water molecules create strong hydrogen bonds with the headgroups of the phospholipid molecules. Due to the presence of these bonds and also due to the confinement of water, the translational and orientational dynamics of water at the interface are slowed down. The degree of slowing down of the dynamics depends upon the location of water molecules near a lipid headgroup.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Interaction of zervamicin IIB with lipid bilayers. Molecular dynamics study

In this work we have studied the interaction of zervamicin IIB (ZrvIIB) with the model membranes of eukaryotes and prokaryotes using all-atom molecular dynamics. In all our simulations zervamicin molecule interacted only with lipid headgroups but did not penetrate the hydrophobic core of the bilayers. During the interaction with the prokaryotic membrane zervamicin placed by its N-termini towards the lipids and rotated at an angle of 40° relatively to the bilayer surface. In the case of eukaryotic membrane zervamicin stayed in the water and located parallel to the membrane surface. We compared hydrogen bonds between peptide and lipids and concluded that interactions of ZrvIIB with prokaryotic membrane are stronger than those with eukaryotic one. Also it was shown that two zervamicin molecules formed dimer and penetrated deeper in the area of lipid headgroups.     

Membrane-forming properties of pseudoglyceryl backbone based gemini lipids possessing oxyethylene spacers

Five pseudoglyceryl backbone based gemini lipids possessing varying lengths of oxyethylene [(-CH2-O-CH2-)n] spacers between cationic ammonium head groups have been synthesized, where n varies from 1 to 5. The membrane-forming properties of these gemini cationic lipids have been investigated. All the gemini lipids formed stable suspensions in water. The presence of membranous aggregates in such lipid suspensions was evidenced by transmission electron microscopy. The membrane-forming characteristics of these gemini lipids were compared with those of the corresponding monomeric lipid with one head group to understand the effect of lipid dimerization. The lipid suspensions were further characterized by dynamic light scattering and zeta potential measurements. Except for the gemini lipid with -CH2-CH2-O-CH2-CH2- spacer (2a), zeta potential of aggregates of all other gemini lipids were significantly greater than that of monomeric lipid suspensions. X-ray diffraction studies with lipid cast films revealed the increase in membrane bilayer width with increase in the length of the spacer (-CH2-O-CH2-)n. Clear thermotropic phase transitions typical of membranous assemblies were observed for all the lipid suspensions by high sensitivity differential scanning calorimetry. Aggregates of gemini lipid 2a bearing one oxyethylene [-(CH2-CH2-O-CH2-CH2)-] unit between headgroups manifested the highest phase transition temperature as compared to other gemini analogues as well as that of monomeric lipid 1. The phase transitions were reversible and exhibited large hysteresis, indicating that the observed phase transitions were of first order. To probe the surface hydration of these membranous aggregates, Paldan fluorescence studies were performed. These studies indicated the high polarity of the vesicular surface of gemini lipid 2a both in the gel and fluid melted phase as compared to vesicles of other gemini lipids.     

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.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Interaction of the sugars trehalose, maltose and glucose with a phospholipid bilayer: a comparative molecular dynamics study

Molecular dynamics simulations are used to investigate the interaction of the sugars trehalose, maltose, and glucose with a phospholipid bilayer at atomic resolution. Simulations of the bilayer in the absence or in the presence of sugar (2 molal concentration for the disaccharides, 4 molal for the monosaccharide) are carried out at 325 and 475 K. At 325 K, the three sugars are found to interact directly with the lipid headgroups through hydrogen bonds, replacing water at about one-fifth to one-quarter of the hydrogen-bonding sites provided by the membrane. Because of its small size and of the reduced topological constraints imposed on the hydroxyl group locations and orientations, glucose interacts more tightly (at identical effective hydroxyl group concentration) with the lipid headgroups when compared to the disaccharides. At high temperature, the three sugars are able to prevent the thermal disruption of the bilayer. This protective effect is correlated with a significant increase in the number of sugar-headgroups hydrogen bonds. For the disaccharides, this change is predominantly due to an increase in the number of sugar molecules bridging three or more lipid molecules. For glucose, it is primarily due to an increase in the number of sugar molecules bound to one or bridging two lipid molecules.     

A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer

A series of molecular dynamics simulations of dimyristoylphosphatidylcholine bilayers, with different levels of hydration and temperature, were performed to examine the influence of hydration on properties of lipid membranes. Structural and dynamical properties such as area per lipid, electron densities, order parameters for all CH bonds and water, diffusion, and reorientation autocorrelation functions were determined and were all found to be affected by changes in the hydration level. The simulations give an overall picture of the bilayer going to a more ordered state when the hydration level is reduced. Lipid headgroups were found to adopt an orientation more parallel to the membrane plane when the water content was decreased. Dynamical properties such as lipid diffusion and relaxation of reorientation time correlation functions were found to become slower with the removal of water. Our simulation results generally agree with experimental data in cases where such data are available. One important conclusion drawn is that while structural properties are affected only moderately dynamical properties are affected very strongly by a decrease of water content.     

Hydration dependent studies of highly aligned multilayer lipid membranes by neutron scattering

We investigated molecular motions on a picosecond timescale of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) model membranes as a function of hydration by using elastic and quasielastic neutron scattering. Two different hydrations corresponding to approximately nine and twelve water molecules per lipid were studied, the latter being the fully hydrated state. In our study, we focused on head group motions by using chain deuterated lipids. Information on in-plane and out-of-plane motions could be extracted by using solid supported DMPC multilayers. Our studies confirm and complete former investigations by Ko?nig et al. [J. Phys. II (France) 2, 1589 (1992)] and Rheinsta?dter et al. [Phys. Rev. Lett. 101, 248106 (2008)] who described the dynamics of lipid membranes, but did not explore the influence of hydration on the head group dynamics as presented here. From the elastic data, a clear shift of the main phase transition from the P(β) ripple phase to the L(α) liquid phase was observed. Decreasing water content moves the transition temperature to higher temperatures. The quasielastic data permit a closer investigation of the different types of head group motion of the two samples. Two different models are needed to fit the elastic incoherent structure factor and corresponding radii were calculated. The presented data show the strong influence hydration has on the head group mobility of DMPC.     

Interaction of cholesterol-conjugated ionizable amino lipids with biomembranes: lipid polymorphism, structure-activity relationship, and implications for siRNA delivery

Delivery of siRNA is a major obstacle to the advancement of RNAi as a novel therapeutic modality. Lipid nanoparticles (LNP) consisting of ionizable amino lipids are being developed as an important delivery platform for siRNAs, and significant efforts are being made to understand the structure-activity relationship (SAR) of the lipids. This article uses a combination of small-angle X-ray scattering (SAXS) and differential scanning calorimetry (DSC) to evaluate the interaction between cholesterol-conjugated ionizable amino lipids and biomembranes, focusing on an important area of lipid SAR--the ability of lipids to destabilize membrane bilayer structures and facilitate endosomal escape. In this study, cholesterol-conjugated amino lipids were found to be effective in increasing the order of biomembranes and also highly effective in inducing phase changes in biological membranes in vitro (i.e., the lamellar to inverted hexagonal phase transition). The phase transition temperatures, determined using SAXS and DSC, serve as an indicator for ranking the potency of lipids to destabilize endosomal membranes. It was found that the bilayer disruption ability of amino lipids depends strongly on the amino lipid concentration in membranes. Amino lipids with systematic variations in headgroups, the extent of ionization, tail length, the degree of unsaturation, and tail asymmetry were evaluated for their bilayer disruption ability to establish SAR. Overall, it was found that the impact of these lipid structure changes on their bilayer disruption ability agrees well with the results from a conceptual molecular "shape" analysis. Implications of the findings from this study for siRNA delivery are discussed. The methods reported here can be used to support the SAR screening of cationic lipids for siRNA delivery, and the information revealed through the study of the interaction between cationic lipids and biomembranes will contribute significantly to the design of more efficient siRNA delivery vehicles.     

Elastic deformation of membrane bilayers probed by deuterium NMR relaxation

In deuterium ((2)H) NMR spectroscopy of fluid lipid bilayers, the average structure is manifested in the segmental order parameters (S(CD)) of the flexible molecules. The corresponding spin-lattice relaxation rates (R(1Z) depend on both the amplitudes and the rates of the segmental fluctuations, and indicate the types of lipid motions. By combining (2)H NMR order parameter measurements with relaxation studies, we have obtained a more comprehensive picture of lipids in the liquid-crystalline (L(alpha)) state than formerly possible. Our data suggest that a lipid bilayer constitutes an ordered fluid, in which the phospholipids are grafted to the aqueous interface via their polar headgroups, whereas the fatty acyl chains are in effect liquid hydrocarbon. Studies of (2)H-labeled saturated lipids indicate their R(1Z) rates and S(CD) order parameters are correlated by a model-free, square-law functional dependence, signifying the presence of relatively slow bilayer fluctuations. A new composite membrane deformation model explains simultaneously the frequency (magnetic field) dependence and the angular anisotropy of the relaxation. The results imply the R(1Z) rates are due to a broad spectrum of 3-D collective bilayer excitations, together with effective axial rotations of the lipids. For the first time, NMR relaxation studies show that the viscoelastic properties of membrane lipids at megahertz frequencies are modulated by the lipid acyl length (bilayer thickness), polar headgroups (bilayer interfacial area), inclusion of a nonionic detergent (C(12)E(8)), and the presence of cholesterol, leading to a range of bilayer softness. Our findings imply the concept of elastic deformation is relevant on lengths approaching the bilayer thickness and less (the mesoscopic scale), and suggest that application of combined R(1Z) and S(CD) studies of phospholipids can be used as a simple membrane elastometer. Heuristic estimates of the bilayer bending rigidity kappa and the area elastic modulus K(a) enable comparison to other biophysical studies, involving macroscopic deformation of thin membrane lipid films. Finally, the bilayer softness may be correlated with the lipid diversity of biomembranes, for example, with regard to membrane curvature, repulsive interactions between bilayers, and lipid-protein interactions.