Membrane dipole potentials, hydration forces, and the ordering of water at membrane surfaces

Many-body polarization effects and the membrane dipole potential

Molecular dynamics simulations of a lipid monolayer at a water-air interface are used to investigate the dipole potential that arises at the water-lipid interface. One simulation explicitly accounts for many-body polarization effects by using a model based on classical Drude oscillators. The dipole potential of the Drude model monolayer is 0.35V in excellent agreement with experimental estimates that range between 0.3 and 0.4V, whereas, a simulation using a nonpolarizable model significantly overestimates the potential with a calculated value of 0.8V. Induced polarization effects in the nonpolar region of the monolayer are found to buffer the residual positive lipid potential that results from competing polarization effects at the polar water/monolayer interface. These results, indicate the utility of the inclusion of many-body polarization effects in empirical force field models of lipids.

Membrane transport in primitive cells

Although model protocellular membranes consisting of monoacyl lipids are similar to membranes composed of contemporary diacyl lipids, they differ in at least one important aspect. Model protocellular membranes allow for the passage of polar solutes and thus can potentially support cell-to functions without the aid of transport machinery. The ability to transport polar molecules likely stems from increased lipid dynamics. Selectively permeable vesicle membranes composed of monoacyl lipids allow for many lifelike processes to emerge from a remarkably small set of molecules.

Gadolinium ions block mechanosensitive channels by altering the packing and lateral pressure of anionic lipids

Effects of polyvalent ions on the lateral packing of phospholipids have been known for decades, but the physiological consequences have not been systematically studied. Gd(3+) is a relatively nonspecific agent that blocks mechano-gated channels with a variable affinity. In this study, we show that the large mechanosensitive channel MscL of Escherichia coli is effectively blocked by Gd(3+) only when reconstituted with negatively charged phospholipids (e.g., PS). Taking this lead, we studied effects of Gd(3+) on monolayers and unilamellar vesicles made of natural brain PS, DMPS, and its mixtures with DMPC. In monolayer experiments, we found that muM Gd(3+) present in the subphase leads to approximately 8% lateral compaction of brain PS (at 35 mN/m). Gd(3+) more strongly shrinks and rigidifies DMPS films causing a spontaneous liquid expanded-to-compact transition to the limiting 40 A(2)/mol. Pressure-area isotherms of uncharged DMPC were unaffected by Gd(3+), and neutralization of DMPS surface by low pH did not produce strong compaction. Upshifts of surface potential isotherms of DMPS monolayers reflected changes in the diffuse double layer due to neutralization of headgroup charges by Gd(3+), whereas the increased packing density produced up to a 200 mV change in the interfacial dipole potential. The slopes of surface potential versus reciprocal area predicted that Gd(3+) induced a modest ( approximately 18%) increase in the magnitude of the individual lipid dipoles in DMPS. Isothermal titration calorimetry indicated that binding of Gd(3+) to DMPS liposomes in the gel state is endothermic, whereas binding to liquid crystalline liposomes produces heat consistent with the isothermal liquid-to-gel phase transition induced by the ion. Both titration curves suggested a K(b) of approximately 10(6) M(-1). We conclude that anionic phospholipids serve as high-affinity receptors for Gd(3+) ions, and the ion-induced compaction generates a lateral pressure increase estimated as tens of mN/m. This pressure can "squeeze" the channel and shift the equilibrium toward the closed state.

Properties of hydrated excess protons near phospholipid bilayers

The behavior of the hydrated excess proton near different lipid membranes is studied with the third generation of the multistate empirical valence bond (MS-EVB3) model [Wu, Y. J.; Chen, H. N.; Wang, F.; Paesani, F.; Voth, G. A. J. Phys. Chem. B 2008, 112, 467]. Dioleoylphosphatidylcholine (DOPC), dioleoylphosphatidylethanolamine (DOPE), and dioleoylphosphatidylglycerol (DOPG) are selected as example lipids. In spite of the differences of the head groups, the molecular dynamics simulations show that all the lipid membranes have a proton-collecting antenna effect with no free energy barrier between the bulk water and interface regions. By comparison with classical hydronium model simulations, it is found that an appropriate description of proton Grotthuss shuttling and associated charge defect delocalization are necessary to obtain the correct free energy profile for the hydrated excess proton. In addition, nanosecond time scale sampling is essential to evaluate the free energy profiles, because certain slow motions are needed to stabilize the excess proton in the deep membrane interface region. It is also found that the lateral diffusion coefficients are 1 order of magnitude smaller in the interface region than in bulk water for all the lipids. These coefficients are almost the same as those of the lipid head groups. Finally, since the lipid phosphates may possibly be protonated due to the proton antenna effect of the membrane, phosphate group protonation is investigated and discussed within the MS-EVB framework.

The fusion peptide domain is the primary membrane-inserted region and enhances membrane interaction of the ectodomain of HIV-1 gp41

To execute the membrane fusion function, it is necessary for the fusion protein of the virus to penetrate into the hydrophobic milieu of membrane bilayer. Hence identification of the region(s) of the ectodomain of viral fusion proteins involved in the membrane insertion and their interaction with the rest of the fusion protein in the membrane would be important for the mechanistic study of membrane fusion. To this end, we examined membrane activity of the fusion peptide, and the ectodomain protein with or without the fusion peptide domain of HIV-1 gp41 by several biophysical measurements. The results revealed that the ectodomain protein containing the fusion peptide domain had higher membrane-perturbing activity and deeper membrane insertion, while the construct lacking the fusion peptide domain had much lower membrane activity. Strikingly, the N-terminal heptad repeat region was found to be induced deeper into the membrane by the fusion peptide, consistent with the role of the latter in the membrane penetration. We concluded that the fusion peptide is the only stretch of gp41 ectodomain that embeds deeply in the membrane interior in the prefusion stage. The function of fusion peptide in terms of membrane interaction and the implications of its interplay with other domains of gp41 on the membrane fusion cascade were discussed.

Mitochondrial membranes with mono- and divalent salt: changes induced by salt ions on structure and dynamics

We employ atomistic simulations to consider how mono- (NaCl) and divalent (CaCl(2)) salt affects properties of inner and outer membranes of mitochondria. We find that the influence of salt on structural properties is rather minute, only weakly affecting lipid packing, conformational ordering, and membrane electrostatic potential. The changes induced by salt are more prominent in dynamical properties related to ion binding and formation of ion-lipid complexes and lipid aggregates, as rotational diffusion of lipids is slowed down by ions, especially in the case of CaCl(2). In the same spirit, lateral diffusion of lipids is slowed down rather considerably for increasing concentration of CaCl(2). Both findings for dynamic properties can be traced to the binding of ions with lipid head groups and the related changes in interaction patterns in the headgroup region, where the binding of Na(+) and Ca(2+) ions is clearly different. The role of cardiolipins in these phenomena turns out to be important.

On the Validation of Molecular Dynamics Simulations of Saturated and cis-Monounsaturated Phosphatidylcholine Lipid Bilayers: A Comparison with Experiment

Molecular dynamics simulations of fully hydrated pure bilayers of four widely studied phospholipids, 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) using a recent revision of the GROMOS96 force field are reported. It is shown that the force field reproduces the structure and the hydration of bilayers formed by each of the four lipids with high accuracy. Specifically, the solvation and the orientation of the dipole of the phosphocholine headgroup and of the ester carbonyls show that the structure of the primary hydration shell in the simulations closely matches experimental findings. This work highlights the need to reproduce a broad range of properties beyond the area per lipid, which is poorly defined experimentally, and to consider the effect of system size and sampling times well beyond those commonly used.

Screening for the drug-phospholipid interaction: correlation to phospholipidosis

Phospholipid bilayers represent a complex, anisotropic environment fundamentally different from bulk oil or octanol, for instance. Even "simple" drug association to phospholipid bilayers can only be fully understood if the slab-of-hydrocarbon approach is abandoned and the complex, anisotropic properties of lipid bilayers reflecting the chemical structures and organization of the constituent phospholipids are considered. The interactions of drugs with phospholipids are important in various processes, such as drug absorption, tissue distribution, and subcellular distribution. In addition, drug-lipid interactions may lead to changes in lipid-dependent protein activities, and further, to functional and morphological changes in cells, a prominent example being the phospholipidosis (PLD) induced by cationic amphiphilic drugs. Herein we briefly review drug-lipid interactions in general and the significance of these interactions in PLD in particular. We also focus on a potential causal connection between drug-induced PLD and steatohepatitis, which is induced by some cationic amphiphilic drugs.

Cholesterol, regulated exocytosis and the physiological fusion machine

Exocytosis is a highly conserved and essential process. Although numerous proteins are involved throughout the exocytotic process, the defining membrane fusion step appears to occur through a lipid-dominated mechanism. Here we review and integrate the current literature on protein and lipid roles in exocytosis, with emphasis on the multiple roles of cholesterol in exocytosis and membrane fusion, in an effort to promote a more molecular systems-level view of the as yet poorly understood process of Ca2+-triggered membrane mergers.

Fusion peptide from influenza hemagglutinin increases membrane surface order: an electron-spin resonance study

A spin-labeling study of interactions of a fusion peptide from the hemagglutinin of the influenza virus, wt20, and a fusion-inactive mutant DeltaG1 with dimyristoylphosphatidylcholine (DMPC) and 1-palmitoyl-2-oleoyl-phosphatdylcholine bilayers was performed. We found that upon binding of wt20, the ordering of headgroups and the ordering of acyl chains near the headgroup increased significantly, in a manner consistent with a cooperative phenomenon. However, changes in the order at the end of the acyl chains were negligible. The ordering effect of wt20 on the headgroup was much stronger at pH 5 than at pH 7. No effect of DeltaG1 binding on the order of bilayers was evident. We also found that 1-palmitoyl-2-hydroxyl phosphatidylcholine, a membrane-fusion inhibitor, decreased the ordering of DMPC headgroups, whereas arachidonic acid, a membrane-fusion promoter, increased the ordering of DMPC headgroups. These results suggest that increases in headgroup ordering may be important for membrane fusion. We propose that upon binding of wt20, which is known to affect only the outer leaflet of the bilayer, this outer leaflet becomes more ordered, and thus more solid-like. Then the coupling between the hardened outer leaflet and the softer inner leaflet generates bending stresses in the bilayer, which tend to increase the negative curvature of the bilayer. We suggest that the increased ordering in the headgroup region enhances dipolar interactions and lowers electrostatic energy, which may provide an energy source for membrane fusion. Possible roles of bending stresses in promoting membrane fusion are discussed.

Configuration of carbonyl groups at the lipid interphases of different topological arrangements of lipid dispersions

The purpose of this work is to analyze the conformation of the carbonyl groups of acyl phospholipids at the hydrocarbon-water interphase in different topological ensembles and phase states, such as micelles and bilayers. The separation of the band components in lipids dispersed in D(2)O is compared with that of PCs in a low hydrated state. When hydrated, the differences in the frequencies of the band components corresponding to the carbonyl groups identified as low hydrated and hydrated populations increase when dimyristoylphosphatidylcholine (DMPC) bilayers go from the lamellar gel to the ripple corrugated phase at the pretransition temperature. Below the pretransition, at which the membrane in the gel state is planar, the two components overlap making the deconvolution unreliable. A further analysis shows that the frequency of the highly hydrated population increases more noticeable than that corresponding to the low hydrated one following the sequence: micelles, fluid phase, ripple gel phase, and lamellar gel phase. This is confirmed by the increase in the separation of the band components when the liposomes are subjected to an osmotic dehydration suggesting that the hydrated population loses water and the dehydrated one partially hydrates. It is concluded that this behavior is a feature conferred by hydration of the different topological arrangements. The relevance of these results on the interphase properties of lipid membranes is discussed.

Monitoring biophysical properties of lipid membranes by environment-sensitive fluorescent probes

We review the main trends in the development of fluorescence probes to obtain information about the structure, dynamics, and interactions in biomembranes. These probes are efficient for studying the microscopic analogs of viscosity, polarity, and hydration, as well as the molecular order, environment relaxation, and electrostatic potentials at the sites of their location. Progress is being made in increasing the information content and spatial resolution of the probe responses. Multichannel environment-sensitive probes that can distinguish between different membrane physicochemical properties through multiple spectroscopic parameters show considerable promise.

Charge equilibration force fields for lipid environments: applications to fully hydrated DPPC bilayers and DMPC-embedded gramicidin A

Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

Nanoscopic description of biomembrane electrostatics: results of molecular dynamics simulations and fluorescence probing

Electrostatic fields generated on and inside biological membranes are recognized to play a fundamental role in key processes of cell functioning. Their understanding requires an adequate description on the level of elementary charges and the reconstruction of electrostatic potentials by integration over all elementary interactions. Out of all the available research tools, only molecular dynamics simulations are capable of this, extending from the atomic to the mesoscopic level of description on the required time and space scale. A complementary approach is that offered by molecular probe methods, with the application of electrochromic dyes. Highly sensitive to intermolecular interactions, they generate integrated signals arising from electric fields produced by elementary charges at the sites of their location. This review is an attempt to provide a critical analysis of these two approaches and their present and potential applications. The results obtained by both methods are consistent in that they both show an extremely complex profile of the electric field in the membrane. The nanoscopic view, with two-dimensional averaging over the bilayer plane and formal separation of the electrostatic potential into surface (Psi(s)), dipole (Psi(d)) and transmembrane (Psi(t)) potentials, is constructive in the analysis of different functional properties of membranes.

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Effects of ether vs. ester linkage on lipid bilayer structure and water permeability

The structure and water permeability of bilayers composed of the ether-linked lipid, dihexadecylphosphatidylcholine (DHPC), were studied and compared with the ester-linked lipid, dipalmitoylphosphaditdylcholine (DPPC). Wide angle X-ray scattering on oriented bilayers in the fluid phase indicate that the area per lipid A is slightly larger for DHPC than for DPPC. Low angle X-ray scattering yields A=65.1A(2) for DHPC at 48 degrees C. LAXS data provide the bending modulus, K(C)=4.2x10(-13)erg, and the Hamaker parameter H=7.2x10(-14)erg for the van der Waals attractive interaction between neighboring bilayers. For the low temperature phases with ordered hydrocarbon chains, we confirm the transition from a tilted L(beta') gel phase to an untilted, interdigitated L(beta)I phase as the sample hydrates at 20 degrees C. Our measurement of water permeability, P(f)=0.022cm/s at 48 degrees C for fluid phase DHPC is slightly smaller than that of DPPC (P(f)=0.027cm/s) at 50 degrees C, consistent with our triple slab theory of permeability.

Ordered water structure at hydrophobic graphite interfaces observed by 4D, ultrafast electron crystallography

Interfacial water has unique properties in various functions. Here, using 4-dimensional (4D), ultrafast electron crystallography with atomic-scale spatial and temporal resolution, we report study of structure and dynamics of interfacial water assembly on a hydrophobic surface. Structurally, vertically stacked bilayers on highly oriented pyrolytic graphite surface were determined to be ordered, contrary to the expectation that the strong hydrogen bonding of water on hydrophobic surfaces would dominate with suppressed interfacial order. Because of its terrace morphology, graphite plays the role of a template. The dynamics is also surprising. After the excitation of graphite by an ultrafast infrared pulse, the interfacial ice structure undergoes nonequilibrium "phase transformation" identified in the hydrogen-bond network through the observation of structural isosbestic point. We provide the time scales involved, the nature of ice-graphite structural dynamics, and relevance to properties related to confined water.

Molecular dynamics simulations of a DMPC bilayer using nonadditive interaction models

We present a polarizable force field based on the charge-equilibration formalism for molecular dynamics simulations of phospholipid bilayers. We discuss refinement of headgroup dihedral potential parameters to reproduce ab initio conformational energies of dimethylphosphate calculated at the MP2/cc-pVTZ level of theory. We also address the refinement of electrostatic and Lennard-Jones (van der Waals) parameters to reproduce ab initio polarizabilities and water interaction energies of dimethylphosphate and tetramethylammonium. We present results of molecular dynamics simulations of a solvated dimyristoylphosphatidylcholine bilayer using this polarizable force field as well as the nonpolarizable, fixed-charge CHARMM27 and CHARMM27r force fields for comparison. Calculated atomic and electron-density profiles, deuterium order parameters, and headgroup orientations are found to be consistent with previous simulations and with experiment. Polarizable interaction models for solvent and lipid exhibit greater water penetration into the lipid interior; this is due to the variation of water molecular dipole moment from a bulk value of 2.6 Debye to a value of 1.9 Debye in the membrane interior. The reduction in the electrostatic component of the desolvation free-energy penalty allows for greater water density. The surface dipole potential predicted by the polarizable model is 0.95 V compared to the value of 0.8 V based on nonpolarizable force-field calculations. Effects of inclusion of explicit polarization are discussed in relation to water dipole moment and varying charge distributions. Dielectric permittivity profiles for polarizable and nonpolarizable interactions exhibit subtle differences arising from the nature of the individual component parameterizations; for the polarizable force field, we obtain a bulk dielectric permittivity of 79, whereas the nonpolarizable force field plateaus at 97 (the value for pure TIP3P water). In the membrane interior, both models predict unit permittivities, with the polarizable models contributing from one to two more units due to the optical dielectric (high-frequency dipole fluctuations). This contribution is a step toward the continuing development of a CHARMM (Chemistry at Harvard Molecular Mechanics) polarizable force field for simulations of biomacromolecular systems.

Effect of headgroup on the dipole potential of phospholipid vesicles

The dipole potentials, psi (d), of phospholipid vesicles composed of pure dimyristoylphosphatidylcholine (DMPC) or vesicles in which 50 mol% of the DMPC was substituted by dimyristoylphosphatidylserine (DMPS), dimyristoylphosphatidylglycerol (DMPG), dimyristoylethanolamine (DMPE), dimyristoylphosphatidic acid (DMPA) or monomyristoylphosphatidylcholine (MMPC) were measured via a fluorescent ratiometric method utilizing the probe di-8-ANEPPS. The PS and PG headgroups were found to cause only minor changes in psi (d). PE caused an increase in psi (d) of 51 mV. This could be explained by a decrease in the dielectric constant of the glycerol backbone region as well as a movement of the P(-)-N(+) dipole of the less bulky PE headgroup to a position more parallel to the membrane surface than in PC. The negatively charged PA headgroup increases psi (d) by 215 mV relative to PC alone. This indicates that the positive pole of the dipole predominantly responsible for the dipole potential is located at a position closer to the interior of the membrane than the phosphate group. The increase in the charge of the negative pole of the dipole by the phosphate group of PA increases the electrical potential drop across the lipid headgroup region. The incorporation of the single chain lipid MMPC into the membrane causes a decrease in psi (d) of 142 mV. This can be explained by a decrease in packing density within the membrane of carbonyl dipoles from the sn-2 chain of DMPC. The results presented should contribute to a better understanding of the electrical effect of lipid headgroups on the functioning of membrane proteins.

Structural and functional properties of hydration and confined water in membrane interfaces

The scope of the present review focuses on the interfacial properties of cell membranes that may establish a link between the membrane and the cytosolic components. We present evidences that the current view of the membrane as a barrier of permeability that contains an aqueous solution of macromolecules may be replaced by one in which the membrane plays a structural and functional role. Although this idea has been previously suggested, the present is the first systematic work that puts into relevance the relation water-membrane in terms of thermodynamic and structural properties of the interphases that cannot be ignored in the understanding of cell function. To pursue this aim, we introduce a new definition of interphase, in which the water is organized in different levels on the surface with different binding energies. Altogether determines the surface free energy necessary for the structural response to changes in the surrounding media. The physical chemical properties of this region are interpreted in terms of hydration water and confined water, which explain the interaction with proteins and could affect the modulation of enzyme activity. Information provided by several methodologies indicates that the organization of the hydration states is not restricted to the membrane plane albeit to a region extending into the cytoplasm, in which polar head groups play a relevant role. In addition, dynamic properties studied by cyclic voltammetry allow one to deduce the energetics of the conformational changes of the lipid head group in relation to the head-head interactions due to the presence of carbonyls and phosphates at the interphase. These groups are, apparently, surrounded by more than one layer of water molecules: a tightly bound shell, that mostly contributes to the dipole potential, and a second one that may be displaced by proteins and osmotic stress. Hydration water around carbonyl and phosphate groups may change by the presence of polyhydroxylated compounds or by changing the chemical groups esterified to the phosphates, mainly choline, ethanolamine or glycerol. Thus, surface membrane properties, such as the dipole potential and the surface pressure, are modulated by the water at the interphase region by changing the structure of the membrane components. An understanding of the properties of the structural water located at the hydration sites and the functional water confined around the polar head groups modulated by the hydrocarbon chains is helpful to interpret and analyze the consequences of water loss at the membranes of dehydrated cells. In this regard, a correlation between the effects of water activity on cell growth and the lipid composition is discussed in terms of the recovery of the cell volume and their viability. Critical analyses of the properties of water at the interface of lipid membranes merging from these results and others from the literature suggest that the interface links the membrane with the aqueous soluble proteins in a functional unit in which the cell may be considered as a complex structure stabilized by water rather than a water solution of macromolecules surrounded by a semi permeable barrier.

Interaction of peptides with biomembranes assessed by potential-sensitive fluorescent probes

Peptide-membrane interaction is an important step to be evaluated in a study of the activity and mode of action of several bioactive peptides. A variety of methods are available; however, few of them satisfy the criteria of being sensitive, biocompatible, versatile, easy to perform, and allowing real-time monitoring as the use of potential-sensitive fluorescent probes. Here we review methods for detecting the effects of membrane-active peptides, even those that are not intrinsically fluorescent, on the different types of membrane potentials, with a special emphasis on studies conducted with living cells. FPE is a probe sensitive to surface potential and detects electrostatic interactions at the water-lipid interface. Di-8-ANEPPS is sensitive to dipole potential and detects membrane incorporations. Transmembrane potential changes reveal major membrane destabilizations, such as in pore formation. The combination of the information obtained from the three potential variations can lead to a more elucidative picture of the mechanisms of the interaction of relevant peptides with biomembranes.

Effect of NaCl and KCl on phosphatidylcholine and phosphatidylethanolamine lipid membranes: insight from atomic-scale simulations for understanding salt-induced effects in the plasma membrane

To gain a better understanding of how monovalent salt under physiological conditions affects plasma membranes, we have performed 200 ns atomic-scale molecular dynamics simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers. These two systems provide representative models for the outer and inner leaflets of the plasma membrane, respectively. The implications of cation-lipid interactions in these lipid systems have been considered in two different aqueous salt solutions, namely NaCl and KCl, and the sensitivity of the results on the details of interactions used for ions is determined by repeating the simulations with two distinctly different force fields. We demonstrate that the main effect of monovalent salt on a phospholipid membrane is determined by cations binding to the carbonyl region of a membrane, while chloride anions mostly stay in the water phase. It turns out that the strength and character of the cation-lipid interactions are quite different for different types of lipids and cations. PC membranes and Na+ ions demonstrate strongest interactions, leading to notable membrane compression. This finding was confirmed by both force fields (Gromacs and Charmm) employed for the ions. The binding of potassium ions to PC membranes (and the overall effect of KCl), in turn, was found to be much weaker mainly due to the larger size of a K+ ion compared to Na+. Furthermore, the effect of KCl on PC membranes was found to be force-field sensitive: The binding of a potassium ion was not observed at all in simulations performed with the Gromacs force-field, which seems to exaggerate the size of a K+ ion. As far as PE lipid bilayers are concerned, they are found to be influenced by monovalent salt to a significantly lesser extent compared to PC bilayers, which is a direct consequence of the ability of PE lipids to form both intra- and intermolecular hydrogen bonds and hence to adopt a more densely packed bilayer structure. Whereas for NaCl we observed weak binding of Na+ cations to the PE lipid-water interface, in the case of KCl we witnessed almost complete lack of cation binding. Overall, our findings indicate that monovalent salt ions affect lipids in the inner and outer leaflets of plasma cell membranes in substantially different ways.

Statistical thermodynamics through computer simulation to characterize phospholipid interactions in membranes

This chapter describes the major issues that are involved in the statistical thermodynamics of phospholipid membranes at the atomic level. The ingredients going into models of lipid bilayers are summarized: force fields, representation of long-range interactions, and boundary conditions. Next, the choice of thermodynamic ensembles, and the two main options for the generation of a representative sample of configurations: molecular dynamics and Monte Carlo are discussed. The final issue that is dealt with describes the various ways the generated ensembles can be analyzed.

A quantitative coarse-grain model for lipid bilayers

A simplified particle-based computer model for hydrated phospholipid bilayers has been developed and applied to quantitatively predict the major physical features of fluid-phase biomembranes. Compared with available coarse-grain methods, three novel aspects are introduced. First, the main electrostatic features of the system are incorporated explicitly via charges and dipoles. Second, water is accurately (yet efficiently) described, on an individual level, by the soft sticky dipole model. Third, hydrocarbon tails are modeled using the anisotropic Gay-Berne potential. Simulations are conducted by rigid-body molecular dynamics. Our technique proves 2 orders of magnitude less demanding of computational resources than traditional atomic-level methodology. Self-assembled bilayers quantitatively reproduce experimental observables such as electron density, compressibility moduli, dipole potential, lipid diffusion, and water permeability. The lateral pressure profile has been calculated, along with the elastic curvature constants of the Helfrich expression for the membrane bending energy; results are consistent with experimental estimates and atomic-level simulation data. Several of the results presented have been obtained for the first time using a coarse-grain method. Our model is also directly compatible with atomic-level force fields, allowing mixed systems to be simulated in a multiscale fashion.

Effects of fluorescent probe NBD-PC on the structure, dynamics and phase transition of DPPC. A molecular dynamics and differential scanning calorimetry study

We present a combined theoretical (molecular dynamics, MD) and experimental (differential scanning calorimetry, DSC) study of the effect of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD) acyl chain-labeled fluorescent phospholipid analogs (C6-NBD-PC and C12-NBD-PC) on 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers. DSC measurements reveal that <1 mol% of NBD-PC causes elimination of the pre-transition and a large loss of cooperativity of the main transition of DPPC. Labeling with C6-NBD-PC or C12-NBD-PC shifts the main transition temperature to lower or higher values, respectively. Following our recent report on the location and dynamics of these probes (BBA 1768 (2007) 467-478) in fluid phase DPPC, we present a detailed analysis of 100-ns MD simulations of systems containing either C6-NBD-PC or C12-NBD-PC, focused on their influence on several properties of the host bilayer. Whereas most monitored parameters are not severely affected for 1.6 mol% of probe, for the higher concentration studied (6.2 mol%) important differences are evident. In agreement with published reports, we observed that the average area per phospholipid molecule increases, whereas DPPC acyl chain order parameters decrease. Moreover, we predict that incorporation of NBD-PC should increase the electrostatic potential across the bilayer and, especially for C12-NBD-PC, slow lateral diffusion of DPPC molecules and rotational mobility of DPPC acyl chains.

Superficially active water in lipid membranes and its influence on the interaction of an aqueous soluble protease

The purpose of this paper is to demonstrate that the interaction of an aqueous soluble enzyme with lipid membranes is influenced by the lipid composition of the interphase. The results show that the interaction of an aqueous soluble protease, Rennet from Mucor miehei, depends on the exposure of the carbonyl and phosphate groups at the membrane interphase. The changes produced by the protease on the surface pressure of monolayers of dimyristoylphosphatidylcholine (DMPC); dioleoylphosphatidylcholine (DOPC); diphytanoylphosphatidylcholine (DPhPC); dipalmitoylphosphatidylcholine (DPPC); di-O-tetradecylphosphatidyl-choline [D(ether)PC]; dimyristoylphosphatidylethanolamine (DMPE); di-O-tetradecyl-phosphatidylethanolamine [D(ether)PE] were measured at different initial surface pressures. The meaning of the DeltaPi vs. Pi curves was interpreted in the light of the concept of interphase given by Defay and Prigogine [R. Defay, I. Prigogine, Surface Tension and Adsorption, John Wiley & Sons, New York, 1966, pp. 273-277] considering the interphase as a bidimensional solution of polar head groups. With this approach, and based on reported evidences that carbonyls and phosphates are the main hydration sites of the lipid membranes, it is suggested that the mechanism of interaction of aqueous soluble protein involves water beyond the hydration shell. At high surface pressure, only water strongly bound to carbonyl and phosphate groups is present and the interaction is not occurring. In contrast, at low surface pressures, the protease-membrane interaction is a function of acyl chain for different polar groups. This is interpreted, as a consequence of the changes in the interfacial tension produced by the displacement of water confined between the hydrated head groups.

Structural determinants of the packing and electrostatic behavior of unsaturated phosphoglycerides

Docosahexaenoic acid-containing phosphoglycerides accumulate preferentially in membranes of the retina, brain, and spermatozoa, but the functional significance of this largely remains to be determined. Previously we compared the physical properties of homogeneous monolayers of these and other phosphoglyceride species to obtain insights into their physiological roles. Particularly noteworthy were the unusually low dipole moments of species having sn-2-docosahexaenoyl chains. In this study, we have investigated the electrostatic and lateral packing properties of related phosphoglycerides and found that: 1), The dipole moment-lowering effect of the docosahexaenoyl group arises from its having a Z double bond at chain position n-3. 2), The large dipole moment-lowering effects at sn-1 of an ether bond to an alkyl or a 1Z alkenyl chain and that of a sn-2-esterified n-3 fatty acid are additive. 3), The 1Z double bond in an alkenyl chain lowers the molecular area of a phosphoglyceride and, concomitantly, makes it less compressible. 4), Ethanolamine-containing phosphoglycerides are generally less compressible than their corresponding choline analogs. Our data showing that relatively small lipid structural changes markedly alter lipid physical properties in fluid phases underscores the need to study the function of peripheral and integral membrane proteins in the presence of appropriate lipid species.

Effect of dipolar moments in domain sizes of lipid bilayers and monolayers

Lipid domains are found in systems such as multicomponent bilayer membranes and single component monolayers at the air-water interface. It was shown by Keller et al. [J. Phys. Chem. 91, 6417 (1987)] that in monolayers, the size of the domains results from balancing the line tension, which favors the formation of a large single circular domain, against the electrostatic cost of assembling the dipolar moments of the lipids. In this paper, we present an exact analytical expression for the electric potential, ion distribution, and electrostatic free energy for different problems consisting of three different slabs with different dielectric constants and Debye lengths, with a circular homogeneous dipolar density in the middle slab. From these solutions, we extend the calculation of domain sizes for monolayers to include the effects of finite ionic strength, dielectric discontinuities (or image charges), and the polarizability of the dipoles and further generalize the calculations to account for domains in lipid bilayers. In monolayers, the size of the domains is dependent on the different dielectric constants but independent of ionic strength. In asymmetric bilayers, where the inner and outer leaflets have different dipolar densities, domains show a strong size dependence with ionic strength, with molecular-sized domains that grow to macroscopic phase separation with increasing ionic strength. We discuss the implications of the results for experiments and briefly consider their relation to other two dimensional systems such as Wigner crystals or heteroepitaxial growth.

Orientation and motion of tryptophan interfacial anchors in membrane-spanning peptides

The tryptophans of integral membrane proteins have been suggested to play specific roles as "interfacial anchors", based on their preference for a location near the lipid head groups. Still, the underlying mechanism behind this behavior remains unclear. NMR experiments can provide an important tool to study this interaction in an actual bilayer environment. Here solid-state deuterium nuclear magnetic resonance was used to study the tryptophans in membrane-spanning model peptides from the WALP family (acetyl-GWW(LA)nWWA-ethanolamide with n = 5 and 6.5) in samples of mechanically aligned dimyristoylphosphatidylcholine (DMPC) bilayers. The data indicate that the tryptophans near the C-terminal end of the peptide display a significantly different behavior from those near the N-terminus. This is reflected prominently in a large difference in the motion experienced by the indoles at either end of the peptide, highlighting the directionality of the helix. Nevertheless, our observations indicate high levels of motional freedom for all tryptophans in these membrane spanning domains that exceed the dynamics for the helix itself. These observations signify that steric and dynamic features of the polypeptide context modulate the tryptophan anchoring in the membrane interface. Measurements of WALP19 in the ether-linked DMPC analogue ditetradecylphosphatidylcholine (missing the lipid carbonyls) show very similar Trp dynamics and suggest similar orientations for some or all of the tryptophans. This suggests that the lipid acyl chain carbonyls play at most a minor role in the anchoring interaction between these Trp residues and the DMPC interfacial region.

Hydration of POPC bilayers studied by 1H-PFG-MAS-NOESY and neutron diffraction

The stability of lipid bilayers is ultimately linked to the hydrophobic effect and the properties of water of hydration. Magic angle spinning (MAS) nuclear Overhauser enhancement spectroscopy (NOESY) with application of pulsed magnetic field gradients (PFG) was used to study the interaction of water with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers in the fluid phase. NOESY cross-relaxation between water and polar groups of lipids, but also with methylene resonances of hydrophobic hydrocarbon chains, has been observed previously. This observation led to speculations that substantial amounts of water may reside in the hydrophobic core of bilayers. Here, the results of a quantitative analysis of cross-relaxation in a lipid 1-palmitoyl-2-oleoyl-sn-glycero-3 phosphocholine (POPC)/water mixture are reported. Coherences were selected via application of pulsed magnetic field gradients. This technique shortens acquisition times of NOESY spectra to 20 min and reduces t (1)-spectral noise, enabling detection of weak crosspeaks, like those between water and lipids, with higher precision than with non-gradient NOESY methods. The analysis showed that water molecules interact almost exclusively with sites of the lipid-water interface, including choline, phosphate, glycerol, and carbonyl groups. The lifetime of lipid-water associations is rather short, on the order of 100 ps, at least one order of magnitude shorter than the lifetime of lipid-lipid associations. The distribution of water molecules over the lipid bilayer was measured at identical water content by neutron diffraction. Water molecules penetrate deep into the interfacial region of bilayers but water concentration in the hydrophobic core is below the detection limit of one water molecule per lipid, in excellent agreement with the cross-relaxation data.

Effect of arbutin on the dipole potential and area per lipid of ester and ether phosphatidylcholine and phosphatidyl ethanolamine monolayers

The present results report for the first time a systematic study of the effect of arbutin on the dipole potential of lipid membranes. The dipole potential and the area per lipid were measured in monolayers of dimyristoylphosphatidylcholine (DMPC), 1,2-di-O-tetradecyl-sn-glycero-3-phosphocholine (dietherPC), dimyristoylphosphatidylethanolamine (DMPE) and 1,2-di-O-tetradecyl-sn-glycero-3-phosphoethanolamine (dietherPE), spread on aqueous solutions of different concentrations of arbutin. The decrease of the dipole potential of DMPC, both in condensed and expanded monolayers, is parallel to an increase in the area per lipid. In contrast, for dietherPC, the area per lipid is not affected, in spite of the fact that arbutin is also able to decrease the dipole potential in a less drastic extent. In the case of DMPE, the response is similar to that observed with dietherPC: the dipole potential decreases, while the area per lipid remains unchanged. However, when the carbonyl groups are absent in phosphatidylethanolamine derivatives such as the dietherPE, the dipole potential is not affected by arbutin, with a small decrease in the area. The effect of arbutin on the dipole potential differs from that of sucrose, trehalose and phloretin and is congruent with previous results obtained by FTIR on its interaction with the CO groups. Arbutin binding is interpreted in terms of the exposure to water of the phosphate and carbonyl groups at the membrane interface of the different monolayers.

Fluorinated polar heads can strikingly increase or invert the dipole moments at the Langmuir monolayer-water boundary: possible effects from headgroup conformations

The dipole potential of lipid monolayers and bilayers is positive toward their nonpolar moiety. In previous papers, we have shown that designed molecules with fluorinated polar heads can invert the polarity of un-ionized Langmuir films. Monolayers of long-chain trifluoroethyl ester RCOOCH2CF3 and trifluoroethyl ether ROCH2CF3 exhibit large negative DeltaV values, shifted by 150-200% from the positive dipole potentials of their non-fluorinated analogs (Petrov and Möhwald J. Phys. Chem. 1996, 100, 18458; Petrov et al. J. Phys. Chem. B 2005, 109, 14102). Here we report large positive surface (dipole) potentials of monolayers of N-trifluoroethyl docosanamide RCONHCH2CF3 and a 300% DeltaV shift with respect to the non-fluorinated N-ethyl docosanamide films. Comparing the dipole potentials and normal dipole moments of the RCONHCH2CF3 and RCOOCH2CF3 monolayers and the maps of the local electrostatic potential (MEP) and lipophilicity (MLP) of their molecules in vacuum, we conclude that the opposite DeltaV shifts and the difference of 1480 mV between the films of these structurally similar amphiphiles seem to be due to strongly different conformations of their heads. The large positive DeltaV values of the N-trifluoroethyl amide monolayer was related to the network of -NH...O=C- bonds fixing the orientation of the hydrophobic delta+C-F3delta- dipoles toward water. The trifluoroethyl ester heads do not form H-bonds and can adjust their energetically optimal conformation orienting the hydrophobic delta+C-F3delta- dipoles toward air. The opposite signs of the dipole potential and the apparent normal dipole moments of the trifluoroethyl ester and ethyl ester monolayers were explained via energy minimization of 36 upright closely packed molecules with "hook-like" heads. The equilibrium architecture of this ensemble shows statistical distribution of the headgroup conformations and a nano-rough monolayer-water boundary as known from X-ray reflectivity experiments and molecular dynamic simulations of phospholipid monolayers and bilayers. The average of the vertical molecular dipole moments at equilibrium agree fairly well with the measured values of mu perpendicular, and the mean molecular area in the ensemble 19.3 A2 matches the value of 18.9 +/- 0.2 A2 determined via X-ray diffraction at gracing incidence surprisingly well. These results reflect the balance of the attractive and repulsive forces between the closely packed "dry" amphiphilic molecules, but a more sophisticated molecular modeling explicitly including water would better serve to reveal the mechanism of the observed effects.

Measurement of lipid forces by X-ray diffraction and osmotic stress

Lipid suspensions in aqueous solutions most often form multilamellar vesicles of uniformly spaced bilayers. Interlamellar spacing is determined by the balance of attractive van der Waals (charge fluctuation) and repulsive forces. This balance of forces, as well as membrane elasticity, can be probed by applied osmotic stress. We describe how osmotic stress can be imposed on multilamellar lipid samples to study lipid interactions.

Comparison of excitation and emission ratiometric fluorescence methods for quantifying the membrane dipole potential

We are interested in developing fluorescence methods for quantifying lateral variations in the dipole potential across cell surfaces. Previous work in this laboratory showed that the ratio of fluorescence intensities of the voltage-sensitive dye di-8-ANEPPS using excitation wavelengths at 420 and 520 nm correlates well with measurements of the dipole potential. In the present work we evaluate the use of di-8-ANEPPS and an emission ratiometric method for measuring dipole potentials, as Bullen and Saggau (Biophys. J. 65 (1999) 2272-2287) have done to follow changes in the membrane potential in the presence of an externally applied field. Emission ratiometric methods have distinct advantages over excitation methods when applied to fluorescence microscopy because only a single wavelength is needed for excitation. We found that unlike the excitation ratio, the emission ratio does not correlate with the dipole potential of vesicles made from different lipids. A difference in the behaviour of the emission ratio in saturated compared to unsaturated lipid vesicles was noted. Furthermore, the emission ratio did not respond in the same way as the excitation ratio when cholesterol, 6-ketocholestanol, 7-ketocholesterol, and phloretin were added to dimyristoylphosphatidylcholine (DMPC) vesicles. We attribute the lack of correlation between the emission ratio and the dipole potential to simultaneous changes in membrane fluidity caused by changes in membrane composition, which do not occur when the electric field is externally applied as in the work of Bullen and Saggau. Di-8-ANEPPS can, thus, only be used via an excitation ratiometric method to quantify the dipole potential.

Human Opioid Peptide Met-Enkephalin Binds to Anionic Phosphatidylserine in High Preference to Zwitterionic Phosphatidylcholine: Natural-Abundance 13C NMR Study on the Binding State in Large Unilamellar Vesicles

A human opioid neuropeptide, Met-enkephalin (M-Enk: Tyr1-Gly2-Gly3-Phe4-Met5), having no net charge binds to anionic phosphatidylserine (PS) in high preference to zwitterionic phosphatidylcholine (PC). The binding mechanism in the PS and PC bilayers was studied on the basis of the inter- and intramolecular interaction data obtained by natural-abundance 13C nuclear magnetic resonance (NMR) of the peptide. Prominent upfield changes of the 13C resonance were observed in the C-terminal residue upon binding to PS, whereas no such marked change was observed upon binding to PC. The upfield chemical shift changes with their characteristic carbon site dependence are ascribed to the electrostatic binding between the peptide C-terminal CO2- and the PS headgroup NH3+. Despite the net negative charge of the PS bilayer surface, M-Enk thus anchors the negatively charged C-terminus. In the N-terminal residue, on the other hand, marked downfield chemical shift changes are observed upon binding to both the PS and PC bilayers, the magnitude of the changes being much larger in the PS system. The downfield changes with their characteristic carbon site dependence are ascribed to the electrostatic binding between the peptide N-terminal NH3+ and the lipid headgroup negative charge(s) (CO2- or PO4- in PS, PO4- in PC). Perturbation on the signal half-widths due to membrane binding also indicates the preferential and deeper binding of M-Enk on the PS membrane surface than on the PC membrane surface. Local charge cancellation takes place efficiently between M-Enk termini and the PS headgroups and compensates for the strong electrostatic hydration of the ionic groups. Distribution of the charged (positive and negative) and uncharged sites in the headgroups along the bilayer normal is responsible for the marked difference between PS and PC headgroups in controlling the binding state of the zwitterionic M-Enk.