The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Solvent history dependence of gramicidin A conformations in hydrated lipid bilayers.

Reconstituted lipid bilayers of dimyristoylphosphatidylcholine (DMPC) and gramicidin A' have been prepared by cosolubilizing gramicidin and DMPC in one of three organic solvent systems followed by vacuum drying and hydration. The conformational state of gramicidin as characterized by 23Na NMR, circular dichroism, and solid state 15N NMR is dependent upon the cosolubilizing solvent system. In particular, two conformational states are described; a state in which Na+ has minimal interactions with the polypeptide, referred to as a nonchannel state, and a state in which Na+ interacts very strongly with the polypeptide, referred to as the channel state. Both of these conformations are intimately associated with the hydrophobic core of the lipid bilayer. Furthermore, both of these states are stable in the bilayer at neutral pH and at a temperature above the bilayer phase transition temperature. These results with gramicidin suggest that the conformation of membrane proteins may be dictated by the conformation before membrane insertion and may be dependent upon the mechanism by which the insertion is accomplished.     

A Molecular Dynamics Study of DMPC Lipid Bilayers Interacting with Dimethylsulfoxide–Water Mixtures

The diffusion process of dimethylsulfoxide (DMSO) through zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer was studied by means of molecular dynamics (MD) simulations. To account for the cryoprotectant concentration difference between the inside and the outside of the cell, dual DMPC lipid bilayers which separate two aqueous reservoirs with and without DMSO were modeled. The initial configuration of the simulation model had DMSO molecules present in one of the aqueous phases (outside the cell) at two different concentrations of ~3 and ~6?mol%. MD simulations were performed on the systems for 50?ns at 323?K and 1?bar. Although the simulation time considered in the study was insufficient for the DMSO molecules to reach the other aqueous phase and equilibrium, early stages of the diffusion process indicated that DMSO molecules had a tendency to diffuse towards the other aqueous phase. The effects of DMSO on bilayer structural characteristics during the diffusion process were investigated. Simulations were analyzed to correlate the following properties of lipid bilayers in the presence of two different aqueous phases: area per lipid, lipid thickness, mass density profiles, lipid tail order parameter and water dipole orientation. Area per lipid calculated for the leaflet facing the aqueous DMSO?Cwater mixture did not show any significant difference compared to area per lipid for the DMSO-free pure DMPC bilayer. Mass density profiles revealed that DMSO molecules had a strong tendency to diffuse toward the aqueous phase with pure water. The lipid tail order parameter calculated for the sn-1 tail of the leaflet facing the aqueous DMSO?Cwater mixture showed that the ordering of lipid tails decreased compared to the leaflet exposed to pure water. However, the ordering of lipid tails in a system where a single bilayer is hydrated by an aqueous DMSO?Cwater mixture is far lower.     

Area per molecule and distribution of water in fully hydrated dilauroylphosphatidylethanolamine bilayers

The area per lipid molecule for fully hydrated dilauroylphosphatidylethanolamine (DLPE) has been obtained in both the gel and liquid-crystalline states by combining wide-angle X-ray diffraction, electron density profiles, and previously published dilatometry results [Wilkinson, D. A., & Nagle, J. F. (1981) Biochemistry 20, 187-192]. The molecular area increases from 41.0 +/- 0.2 to 49.1 +/- 1.2 A2 upon melting from the gel to liquid-crystalline phase. The thickness of the bilayer, as measured from the electron density profiles, decreases about 4 A upon melting, from 45.2 +/- 0.3 to 41.0 +/- 0.6 A. A somewhat unexpected result is that the fluid layer between fully hydrated bilayers is the same in both gel and liquid-crystalline phases and is only about 5 A thick. From these data, plus the volume of the anhydrous DLPE molecule, it is possible to determine the number of water molecules per lipid and their approximate distribution relative to the lipid molecule. Our analysis shows that there are about 7 and 9 waters per DLPE molecule in the gel and liquid-crystalline phases, respectively. About half of the water is located in the fluid space between adjacent bilayers, and the remaining waters are intercalated into the bilayer, presumably in the head group region. There are significantly fewer water molecules in the fluid spaces between DLPE bilayers than in the fluid spaces in gel- or liquid-crystalline-phase phosphatidylcholine bilayers. This small fluid space in PE bilayers could arise from interbilayer hydrogen bond formation through the water molecules or electrostatic interactions between the amine and phosphate groups on apposing bilayers.     

Effect of phosphatidylethanolamine on the properties of phospholipid-apolipoprotein complexes

Plasma high density lipoproteins (HDL) are synthesized in intestinal mucosal cells and hepatocytes and are secreted into the blood. Factors influencing the structure and function of these HDL, such as lipid and protein composition, are poorly understood. It appears, however, that intracellular, discoidal HDL are enriched, relative to plasma HDL, in phosphatidylethanolamine (PE), a phospholipid known to generate unusual, nonbilayer structures of putative physiological significance. Although incubation of dimyristoylphosphatidylcholine (DMPC) with apolipoprotein A-I at the gel-liquid crystalline phase transition temperature results in the spontaneous formation of lipid-protein complexes, the presence of proportionately small amounts of PE prevents the formation of such complexes, suggesting that PE profoundly alters the phase properties of the phospholipid bilayers. However, by using a detergent-mediated method for the formation of PE-rich model nascent HDL from phospholipids and apolipoprotein A-I, lipid-protein complexes containing as much as 75% DLPE could be formed, thus demonstrating that the presence of PE causes a kinetic, rather than a thermodynamic, barrier to spontaneous complex formation. The products contained a DLPE:DMPC molar ratio similar to that of the initial incubation mixture; however, as the mole percentage of DLPE increased, the products became less heterogeneous, the buoyant density of the products increased, and the Stokes diameter of the products decreased. Similar results were obtained when dimyristoylphosphatidylethanolamine (DMPE) and dipalmitoylphosphatidylethanolamine (DPPE) were employed in lieu of DLPE. Electron microscopy of complexes containing DLPE and DMPC at a 1:1 molar ratio showed that these particles possessed a discoidal, bilayer morphology similar to that seen with complexes containing only phosphatidylcholine.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interpretation of small angle X-ray measurements guided by molecular dynamics simulations of lipid bilayers

Reconstruction and interpretation of lipid bilayer structure from X-ray scattering often rely on assumptions regarding the molecular distributions across the bilayer. It is usually assumed that changes in head-head spacings across the bilayer, as measured from electron density profiles, equal the variations in hydrocarbon thicknesses. One can then determine the structure of a bilayer by comparison to the known structure of a lipid with the same headgroup. Here we examine this procedure using simulated electron density profiles for the benchmark lipids DMPC and DPPC. We compare simulation and experiment in both real and Fourier space to address two main aspects: (i) the measurement of head-head spacings from relative electron density profiles, and (ii) the determination of the absolute scale for these profiles. We find supporting evidence for the experimental procedure, thus explaining the robustness and consistency of experimental structural results derived from electron density profiles. However, we also expose potential pitfalls in the Fourier reconstruction that are due to the limited number of scattering peaks. Volumetric analysis of simulated bilayers allows us to propose an improved, yet simple method for scale determination. In this way we are able to remove some of the restrictions imposed by limited scattering data in constructing reliable electron density profiles.     

Effects of a carane derivative local anesthetic on a phospholipid bilayer studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of two hydrated palmitoyloleoylphosphatidylcholine (POPC) bilayers each containing eight carane derivative (KP-23) local anesthetic (LA) molecules in neutral (POPC-LA) or protonated (POPC-LAH) forms were carried out to investigate the effect of KP-23 and its protonation on the bilayer. 3-ns trajectories were used for analyses. A pure POPC bilayer was employed as a reference system. In both POPC-LA and POPC-LAH systems a few KP-23 molecules intercalated into the bilayer and moved near the bilayer/water interface. They were located on the hydrophobic core side of the interface in the POPC-LA bilayer, but on the water phase side in the POPC-LAH bilayer. The order of the POPC chains was higher in the POPC-LA bilayer than in the pure POPC bilayer and was lower in the POPC-LAH bilayer. Interactions between polar groups of KP-23 and POPC or water were responsible for a lower hydration of POPC headgroups in POPC bilayers containing KP-23 than in the pure POPC bilayer. KP-23 molecules were found to form aggregates both in POPC-LA and POPC-LAH bilayers. Due to higher amphiphilicity of LAH, the LAH aggregate was more micelle-like and larger than the LA one. The results demonstrate the rapid timescales of the initial processes that take place at and near the bilayer interface as well as details of the atomic level interactions between local anesthetic and the lipid matrix of a cell membrane.     

Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation

Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a “ring” of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7–0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

The influence of 1-alkanols and external pressure on the lateral pressure profiles of lipid bilayers

The suggestion by Robert Cantor, that drug-induced pressure changes in lipid bilayers can change the conformational equilibrium between open and closed states of membrane proteins and thereby cause anesthesia, attracted much attention lately. Here, we studied the effect of both large external pressure and of 1-alkanols of different chain lengths—some of them anesthetics, others not—on the lateral pressure profiles across dimyristoylphosphatidylcholine (DMPC) bilayers by molecular dynamics simulations. For a pure DMPC bilayer, high pressure both reduced and broadened the tension at the interface hydrophobic/hydrophilic and diminished the repulsion between the phospholipid headgroups. Whereas the effect of ethanol on the lateral pressure profile was similar to the effect of a large external pressure on a DMPC bilayer, long-chain 1-alkanols significantly amplified local maxima and minima in the lateral pressure profile. For most 1-alkanols, external pressure had moderate effects and did not reverse the changes 1-alkanols exerted on the pressure profile. Nevertheless, assuming the bent helix model as a simple geometric model for the transmembrane region of a membrane protein, protein conformational equilibria were shifted in opposite directions by addition of 1-alkanols and additional application of external pressure.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

A phospholipid bilayer supported under a polymerized Langmuir film

Phospholipid single bilayers supported on a hydrophilic solid substrate are extensively used in the study of the interaction between model membranes and proteins or polypeptides. In this article, the formation of a single dimyristoylphosphatidylcholine (DMPC) bilayer under an octadecyltrimethoxysilane (OTMS) polymerized Langmuir monolayer at the air-water interface is followed by Brewster angle microscopy (BAM) and polarization-modulated infrared reflection absorption spectroscopy (PM-IRRAS). The formation of the bilayer is initiated by injection of dimyristoylphosphatidylcholine small unilamellar vesicles into the aqueous subphase. Brewster angle microscopy allows visualization of the kinetics of formation and the homogeneity of the bilayer. Spectral simulations of the polarization-modulated infrared reflection absorption spectroscopy spectra reveal that the bilayer thickness is 39 +/- 5 A. This system constitutes the first example of a phospholipid bilayer on a "nanoscopic" support and opens the way to studies involving supported bilayers using powerful experimental techniques such as x-ray reflectivity, vibrational spectroscopies, or Brewster angle microscopy.     

Study of curcumin behavior in two different lipid bilayer models of liposomal curcumin using molecular dynamics simulation

Liposomal formulation of curcumin is an important therapeutic agent for the treatment of various cancers. Despite extensive studies on the biological effects of this formulation in cancer treatment, much remains unknown about curcumin–liposome interactions. Understanding how different lipid bilayers respond to curcumin molecule may help us to design more effective liposomal curcumin. Here, we used molecular dynamics simulation method to investigate the behavior of curcumin in two lipid bilayers commonly used in preparation of liposomal curcumin, namely dipalmitoylphosphatidylcholine (DPPC) and dimyristoylphosphatidylglycerol (DMPG). First, the free energy barriers for translocation of one curcumin molecule from water to the lipid bilayer were determined by using the potential of mean force (PMF). The computed free energy profile exhibits a global minimum at the solvent–headgroup interface (LH region) for both lipid membranes. We also evaluated the free energy difference between the equilibrium position of curcumin in the lipid bilayer and bulk water as the excess chemical potential. Our results show that curcumin has the higher affinity in DMPG compared to DPPC lipid bilayer (?8.39 vs. ?1.69 k_BT) and this is related to more hydrogen bond possibility for curcumin in DMPG lipid membrane. Next, using an unconstrained molecular dynamic simulation with curcumin initially positioned at the center of lipid bilayer, we studied various properties of each lipid bilayer system in the presence of curcumin molecule that was in full agreement with PMF and experimental data. The results of these simulation studies suggest that membrane composition could have a large effect on interaction of curcumin–lipid bilayer.     

Variation in hydration forces between neutral phospholipid bilayers: evidence for hydration attraction

It is now generally recognized that hydration forces dominate close interactions of lipid hydrophilic surfaces. The commonality of their characteristics has been reasonably established. However, differences in measured net repulsion, particularly evident when phosphatidylethanolamine (PE) and phosphatidylcholine (PC) bilayers are compared, suggest there exists a variety of behavior wider than expected from earlier models of hydration and fluctuation repulsion balanced by van der Waals attraction. To find a basis for this diverse behavior, we have looked more closely at measured structural parameters, degrees of hydration, and interbilayer repulsive forces for the lamellar phases of the following lipids: 1-palmitoyl-2-oleoyl-PE (POPE), egg PE, transphosphatidylated egg PE (egg PE-T), mono- and dimethylated egg PE-T (MMPE and DMPE), 1-stearoyl-2-oleoyl-PC (SOPC), and mixtures of POPE and SOPC. POPE and SOPC bilayers differ not only in their maximum degrees of hydration but also in the empirical hydration force coefficients and decay lengths that characterize their interaction. When mixed with POPE, SOPC effects sudden and disproportionate increases in hydration. POPE, egg PE, and egg PE-T differ in their degree of hydration, molecular area, and hydration repulsion. A single methylation of egg PE-T almost completely converts its hydration and bilayer repulsive properties to those of egg PC; little progression of hydration is seen with successive methylations. In order to reconcile these observations with the conventional scheme of balancing interbilayer hydration and fluctuation-enhanced repulsion with van der Waals attraction, it is necessary to relinquish the fundamental idea that the decay of hydration forces is a constant determined by the properties of the aqueous medium. Alternatively, one can retain that fundamental idea if one recognizes the possibility that polar group hydration has an attractive component to it. In the latter view, that attractive component originates from interbilayer hydrogen-bonded water bridges between apposing bilayer surfaces, arising from correlation of zwitterionic or other complementary polar groups or from factors that affect polar group solubility. The same Marcelja and Radic formalism that accounts so well for the repulsive component also leads to an estimate of the attractive one. We suggest that the full range of degrees of hydration and of interbilayer spacings observed for different neutral bilayers results in part from variable contributions of the attractive and repulsive hydration components.(ABSTRACT TRUNCATED AT 400 WORDS)     

Examining the Origins of the Hydration Force Between Lipid Bilayers Using All-Atom Simulations

Using 237 all-atom double bilayer simulations, we examined the thermodynamic and structural changes that occur as a phosphatidylcholine lipid bilayer stack is dehydrated. The simulated system represents a micropatch of lipid multilayer systems that are studied experimentally using surface force apparatus, atomic force microscopy and osmotic pressure studies. In these experiments, the hydration level of the system is varied, changing the separation between the bilayers, in order to understand the forces that the bilayers feel as they are brought together. These studies have found a curious, strongly repulsive force when the bilayers are very close to each other, which has been termed the “hydration force,” though the origins of this force are not clearly understood. We computationally reproduce this repulsive, relatively free energy change as bilayers come together and make qualitative conclusions as to the enthalpic and entropic origins of the free energy change. This analysis is supported by data showing structural changes in the waters, lipids and salts that have also been seen in experimental work. Increases in solvent ordering as the bilayers are dehydrated are found to be essential in causing the repulsion as the bilayers come together.     

Effects of cyclosporin A on model lipid membranes

Cyclosporin A (CSA) is a widely used immunosuppressant drug for transplant therapy, however its limitation is its toxicity. The effect of CSA on model membranes such as dimyristoyl phosphatidylcholine (DMPC) bilayers was studied using small-angle X-ray diffraction and differential scanning calorimetry (DSC). CSA abolishes the pretransition and affects the transition of DMPC model membranes in a concentration-related manner as is shown by DSC. CSA induces a second peak at the high temperature side of the main transition, which is interpreted as a phase separation between areas rich and poor in CSA concentration. Small angle X-ray diffraction shows that the repeat distance of the DMPC bilayers in the lamellar Lalpha state increases as a function of concentration up to 10 mol% and remains constant thereafter. Furthermore, CSA affects the fatty acyl chains of the bilayer, especially the part of the chain proximal to the head group. In conclusion, CSA, as both small-angle X-ray diffraction and DSC show, affects in a concentration-wise manner the DMPC model membranes and perturbs the bilayer, in particular the acyl chain region.     

Human cells and cell membrane molecular models are affected in vitro by chlorpromazine

This study presents evidence that chlorpromazine (CPZ) affects human cells and cell membrane molecular models. Human SH-SY5Y neuroblastoma cells incubated with 0.1 mM CPZ suffered a decrease of cell viability. On the other hand, phase contrast microscopy observations of human erythrocytes indicated that they underwent a morphological alteration as 1 μM CPZ changed their discoid normal shape to stomatocytes, and to hemolysis with 1 mM CPZ. X-ray diffraction experiments performed on dimyristoylphosphatidylcholine (DMPC) and dimyristoylphosphatidylethanolamine (DMPE) bilayers, classes of the major phospholipids present in the outer and inner sides of the erythrocyte membrane, respectively showed that CPZ disordered the polar head and acyl chain regions of both DMPC and DMPE, where these interactions were stronger with DMPC bilayers. Fluorescence spectroscopy on DMPC LUV at 18 °C confirmed these results. In fact, the assays showed that CPZ induced a significant reduction of their generalized polarization (GP) and anisotropy (r) values, indicative of enhanced disorder at the polar head and acyl chain regions of the DMPC lipid bilayer.     

Lateral interactions among phosphatidylcholine and phosphatidylethanolamine head groups in phospholipid monolayers and bilayers

We develop theory for the lateral interactions among the zwitterionic head groups of phospholipids in monolayers and bilayers, particularly phosphatidylcholine (PC) and phosphatidylethanolamine (PE). With the P- end of the head group anchored at the water/hydrocarbon interface, a balance of two effects dictates the angle that the P--N+ dipole makes with respect to the plane of the bilayer: N+ is driven toward water due to the (Born) electrostatic free energy, but the hydrophobic effect drives the methyl and methylene groups around the N+ charge toward the hydrocarbon. The only adjustable parameter of the model is the average fluctuation of the oil/water interface or, alternatively, the dielectric constant of the hydrocarbon phase. The model predicts that at 5 degrees C the head group dipole should lie largely in the bilayer plane, in accord with X-ray, neutron diffraction, and NMR studies. The theory makes the novel prediction that the N+ end of the dipole becomes increasingly submerged in hydrocarbon with increasing temperature, leading to strongly enhanced lateral repulsion between PC head groups. This prediction is in good agreement with second and third viral coefficients of monolayer lateral pressures, and with the temperature dependence of the former. The theoretical model is consistent with head group fluctuations measured by neutron diffraction of PC and PE bilayers. Because PE has a smaller hydrophobic cluster near N+, its lateral repulsion should be much smaller and less temperature dependent than for PC, also in agreement with equation-of-state measurements. This suggests why at high density PE monolayers have higher melting temperatures than PC monolayers and more propensity for reversed curvature.     

Interaction of electric dipoles with phospholipid head groups. A 2H and 31P NMR study of phloretin and phloretin analogues in phosphatidylcholine membranes.

Phloretin, 4-hydroxyvalerophenone, and 2-hydroxy-omega-phenylpropiophenone are lipophilic dipolar substances that modify ionic conductances of bilayer membranes. The structural changes at the level of the head groups and the hydrocarbon chains as induced by the incorporation of phloretin and its analogues were investigated with deuterium and phosphorus nuclear magnetic resonance. Membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) were selectively deuterated at the choline head group and at the hydrocarbon chains, and 2H and 31P NMR spectra were recorded with varying concentrations of dipolar agents. Incorporation of phloretin leaves the bilayer structure intact, induces only a small disordering of the hydrocarbon chains, and has no significant effect on the head-group dynamics. On the other hand, quite distinct structural changes are observed for the phosphocholine head group. While the -P-N+ dipole is oriented approximately parallel to the membrane surface for pure POPC bilayers, addition of phloretin, and to a lesser extent 4-hydroxyvalerophenone and 2-hydroxy-omega-phenylpropiophenone, rotates the N+ end of the -P-N+ dipole closer to the hydrocarbon layer. The resulting normal component of the -P-N+ dipole partly compensates the electric field of the dipolar agents. In addition to this structural change, phloretin also modifies the hydration layer at the lipid-water interface. Much less 2H2O is adsorbed to the membrane surface when the bilayer contains phloretin, 4-hydroxyvalerophenone, or 2-hydroxy-omega-phenylpropiophenone. Moreover, a rather large change in the residual phosphorus chemical shielding anisotropy argues in favor of hydrogen-bond formation between the phosphate segment and the phloretin hydroxyl groups.     

Lipid hydration and mobility: an interplay between fluorescence solvent relaxation experiments and molecular dynamics simulations

Fluorescence solvent relaxation experiments are based on the characterization of time-dependent shifts in the fluorescence emission of a chromophore, yielding polarity and viscosity information about the chromophore’s immediate environment. A chromophore applied to a phospholipid bilayer at a well-defined location (with respect to the z-axis of the bilayer) allows monitoring of the hydration and mobility of the probed segment of the lipid molecules. Specifically, time-resolved fluorescence experiments, fluorescence quenching data and molecular dynamic (MD) simulations show that 6-lauroyl-2-dimethylaminonaphthalene (Laurdan) probes the hydration and mobility of the sn-1 acyl groups in a phosphatidylcholine bilayer. The time-dependent fluorescence shift (TDFS) of Laurdan provides information on headgroup compression and expansion induced by the addition of different amounts of cationic lipids to phosphatidylcholine bilayers. Those changes were predicted by previous MD simulations. Addition of truncated oxidized phospholipids leads to increased mobility and hydration at the sn-1 acyl level. This experimental finding can be explained by MD simulations, which indicate that the truncated chains of the oxidized lipid molecules are looping back into aqueous phase, hence creating voids below the glycerol level. Fluorescence solvent relaxation experiments are also useful in understanding salt effects on the structure and dynamics of lipid bilayers. For example, such experiments demonstrate that large anions increase hydration and mobility at the sn-1 acyl level of phosphatidylcholine bilayers, an observation which could not be explained by standard MD simulations. If polarizability is introduced into the applied force field, however, MD simulations show that big soft polarizable anions are able to interact with the hydrophilic/hydrophobic interface of the lipid bilayer, penetrating to the level probed by Laurdan, and that they expand and destabilize the bilayer making it more hydrated and mobile.