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

Ratiometric fluorescence nanoscopy and lifetime imaging of novel Nile Red analogs for analysis of membrane packing in living cells

Subcellular membranes have complex lipid and protein compositions, which give rise to organelle-specific membrane packing, fluidity, and permeability. Due to its exquisite solvent sensitivity, the lipophilic fluorescence dye Nile Red has been used extensively to study membrane packing and polarity. Further improvement of Nile Red can be achieved by introducing electron-donating or withdrawing functional groups. Here, we compare the potential of derivatives of Nile Red with such functional substitutions for super-resolution fluorescence microscopy of lipid packing in model membranes and living cells. All studied Nile Red derivatives exhibit cholesterol-dependent fluorescence changes in model membranes, as shown by spectrally resolved stimulated emission depletion (STED) microscopy. STED imaging of Nile Red probes in cells reveals lower membrane packing in fibroblasts from healthy subjects compared to those from patients suffering from Niemann Pick type C1 (NPC1) disease, a lysosomal storage disorder with accumulation of cholesterol and sphingolipids in late endosomes and lysosomes. We also find small but consistent changes in the fluorescence lifetime of the Nile Red derivatives in NPC1 cells, suggesting altered hydrogen-bonding capacity in their membranes. All Nile Red derivatives are essentially non-fluorescent in water but increase their brightness in membranes, allowing for their use in MINFLUX single molecule tracking experiments. Our study uncovers the potential of Nile Red probes with functional substitutions for nanoscopic membrane imaging.

The cellular zeta potential: cell electrophysiology beyond the membrane

The standard model of the cell membrane potential Vm describes it as arising from diffusion currents across a membrane with a constant electric field, with zero electric field outside the cell membrane. However, the influence of Vm has been shown to extend into the extracellular space where it alters the cell's ζ-potential, the electrical potential measured a few nm from the cell surface which defines how the cell interacts with charged entities in its environment, including ions, molecules, and other cells. The paradigm arising from surface science is that the ζ-potential arises only from fixed membrane surface charge, and has consequently received little interest. However, if the ζ-potential can mechanistically and dynamically change by alteration of Vm, it allows the cell to dynamically alter cell-cell and cell-molecule interactions and may explain previously unexplained electrophysiological behaviours. Whilst the two potentials Vm and ζ are rarely reported together, they are occasionally described in different studies for the same cell type. By considering published data on these parameters across multiple cell types, as well as incidences of unexplained but seemingly functional Vm changes correlating with changes in cell behaviour, evidence is presented that this may play a functional role in the physiology of red blood cells, macrophages, platelets, sperm, ova, bacteria and cancer. Understanding how these properties will improve understanding of the role of electrical potentials and charges in the regulation of cell function and in the way in which cells interact with their environment. Insight  The zeta (ζ) potential is the electrical potential a few nm beyond the surface of any suspensoid in water. Whilst typically assumed to arise only from fixed charges on the cell surface, recent and historical evidence shows a strong link to the cell's membrane potential Vm, which the cell can alter mechanistically through the use of ion channels. Whilst these two potentials have rarely been studied simultaneously, this review collates data across multiple studies reporting Vm, ζ-potential, electrical properties of changes in cell behaviour. Collectively, this points to Vm-mediated ζ-potential playing a significant role in the physiology and activity of blood cells, immune response, developmental biology and egg fertilization, and cancer among others.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

Chromone-Containing Allylmorpholines Influence Ion Channels in Lipid Membranes via Dipole Potential and Packing Stress

Herein, we report that chromone-containing allylmorpholines can affect ion channels formed by pore-forming antibiotics in model lipid membranes, which correlates with their ability to influence membrane boundary potential and lipid-packing stress. At 100 µg/mL, allylmorpholines 1, 6, 7, and 8 decrease the boundary potential of the bilayers composed of palmitoyloleoylphosphocholine (POPC) by about 100 mV. At the same time, the compounds do not affect the zeta-potential of POPC liposomes, but reduce the membrane dipole potential by 80-120 mV. The allylmorpholine-induced drop in the dipole potential produce 10-30% enhancement in the conductance of gramicidin A channels. Chromone-containing allylmorpholines also affect the thermotropic behavior of dipalmytoylphosphocholine (DPPC), abolishing the pretransition, lowering melting cooperativity, and turning the main phase transition peak into a multicomponent profile. Compounds 4, 6, 7, and 8 are able to decrease DPPC's melting temperature by about 0.5-1.9 °C. Moreover, derivative 7 is shown to increase the temperature of transition of palmitoyloleoylphosphoethanolamine from lamellar to inverted hexagonal phase. The effects on lipid-phase transitions are attributed to the changes in the spontaneous curvature stress. Alterations in lipid packing induced by allylmorpholines are believed to potentiate the pore-forming ability of amphotericin B and gramicidin A by several times.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Curvature increases permeability of the plasma membrane for ions, water and the anti-cancer drugs cisplatin and gemcitabine

In this work the permeability of a model asymmetric plasma membrane, for ions, water and the anti-cancer drugs cisplatin and gemcitabine is studied by means of all-atom molecular dynamics simulations. It is shown for the first time that permeability of the highly curved membrane increases from one to three orders of magnitude upon membrane bending depending on the compound and the sign of curvature. Our results suggest that the membrane curvature could be an important factor of drug translocation through the membrane.

Blocking ion channels induced by antifungal lipopeptide syringomycin E with amide-linked local anesthetics

The effects of the amide-linked (lidocaine (LDC), mepivacaine (MPV), prilocaine (PLC)) and ester-bound local anesthetics (benzocaine (BZC), procaine (PRC), and tetracaine (TTC)) on the pore-forming activity of the antifungal lipopeptide syringomycin E (SRE) in lipid bilayers were studied. Independently on electrolyte concentration in the membrane bathing solution the observed changes in conductance of SRE channels agreed with the altered membrane dipole potential under the action of ester-bound local anesthetics. Effects of aminoamides in diluted and concentrated solutions were completely different. At 0.1 M KCl (pH 7.4) the effects of amide-linked anesthetics were in accordance with changes in the membrane surface potential, while at 2 M KCl aminoamides blocked ion passage through the SRE channels, leading to sharp reductions in pore conductance at negative voltages and 100-fold decreases in the channel lifetimes. The effects were not practically influenced by the membrane lipid composition. The interaction cooperativity implied the existence of specific binding sites for amide-bound anesthetics in SRE channels.

The influence of curvature on the properties of the plasma membrane. Insights from atomistic molecular dynamics simulations

In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.

Changes in lipid membrane mechanics induced by di- and tri-phenyltins

Organotin compounds, being biologically active, affect a variety of cellular functions due to their ability to accumulate in and penetrate biological membranes. These compounds influence the distribution of electrostatic charges, alter organization, disrupt molecular dynamics and change mechanical properties of biological membranes. It was found that the membrane/water partition coefficient equals 4, a value significantly higher than octanol/water partition coefficient. In addition, the effect of di- and tri-phenyltin chlorides on the mechanics of model lipid membranes was measured for the first time. It has been determined that phenyltins affect the global model lipid bilayer properties by reducing the membrane expansion modulus, when measured using micromanipulation technique, and elevating the bending rigidity coefficient of the lipid bilayer, as determined with the flickering noise spectroscopy. In addition, the elevated water permeability shows that phenyltins also cause the local defects formation in the lipid bilayer, i.e. lipid pores. These data shows that phenyltins may interfere indirectly with variety cellular processes by altering non-specifically the entire cellular membrane system. Accordingly, when phenyltins are added to macrophages in culture, they inflict massive alterations of cell morphology and interfere with membrane-associated processes, as visualized using fluorescence labelling of selected subcellular compartments.

Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study

Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.

The plasma membrane as a capacitor for energy and metabolism

When considering which components of the cell are the most critical to function and physiology, we naturally focus on the nucleus, the mitochondria that regulate energy and apoptotic signaling, or other organelles such as the endoplasmic reticulum, Golgi, ribosomes, etc. Few people will suggest that the membrane is the most critical element of a cell in terms of function and physiology. Those that consider the membrane critical will point to its obvious barrier function regulated by the lipid bilayer and numerous ion channels that regulate homeostatic gradients. What becomes evident upon closer inspection is that not all membranes are created equal and that there are lipid-rich microdomains that serve as platforms of signaling and a means of communication with the intracellular environment. In this review, we explore the evolution of membranes, focus on lipid-rich microdomains, and advance the novel concept that membranes serve as "capacitors for energy and metabolism." Within this framework, the membrane then is the primary and critical regulator of stress and disease adaptation of the cell.

Effects of Dipole Potential Modifiers on Heterogenic Lipid Bilayers

In this work, we examine the ability of dipole modifiers, flavonoids, and RH dyes to affect the dipole potential (φ d) and phase separation in membranes composed of ternary mixtures of POPC with different sphingolipids and sterols. Changes in the steady-state conductance induced by cation-ionophore complexes have been measured to evaluate the changes in dipole potential of planar lipid bilayers. Confocal fluorescence microscopy has been employed to investigate lipid segregation in giant unilamellar vesicles. The effects of flavonoids on φ d depend on lipid composition and dipole modifier type. The effectiveness of RH dyes to increase φ d depends on sphingolipid type but is not influenced by sterol content. Tested modifiers lead to partial or complete disruption of gel domains in bilayers composed of POPC, sphingomyelin, and cholesterol. Substitution of cholesterol to ergosterol or 7-dehydrocholesterol leads to a loss of fluidizing effects of modifiers except phloretin. This may be due to various compositions of gel domains. The lack of influence of modifiers on phase scenario in vesicles composed of ternary mixtures of POPC, cholesterol, and phytosphingosine or sphinganine is related to an absence of gel-like phase. It was concluded that the membrane lateral heterogeneity affects the dipole-modifying abilities of the agents that influence the magnitude of φ d by intercalation into the bilayer and orientation of its own large dipole moments (phloretin and RH dyes). The efficacy of modifiers that do not penetrate deeply and affect φ d through water adsorption (phlorizin, quercetin, and myricetin) is not influenced by lateral heterogeneity of membrane.

Membrane permeation of a peptide: it is better to be positive

A joint experimental and computational study investigates the translocation of a tryptophan molecule through a phospholipid membrane. Time dependent spectroscopy of the tryptophan side chain determines the rate of permeation into 150 nm phospholipid vesicles. Atomically detailed simulations are conducted to calculate the free energy profiles and the permeation coefficient. Different charging conditions of the peptide (positive, negative, or zwitterion) are considered. Both experiment and simulation reproduce the qualitative trend and suggest that the fastest permeation is when the tryptophan is positively charged. The permeation mechanism, which is revealed by molecular dynamics simulations, is of a translocation assisted by a local defect. The influence of long-range electrostatic interactions, such as the membrane dipole potential on the permeation process, is not significant.

Modifiers of membrane dipole potentials as tools for investigating ion channel formation and functioning

Electrostatic fields generated on and within biological membranes play a fundamental role in key processes in cell functions. The role of the membrane dipole potential is of particular interest because of its powerful impact on membrane permeability and lipid-protein interactions, including protein insertion, oligomerization, and function. The membrane dipole potential is defined by the orientation of electric dipoles of lipid headgroups, fatty acid carbonyl groups, and membrane-adsorbed water. As a result, the membrane interior is several hundred millivolts more positive than the external aqueous phase. This potential decrease depends on the lipid, and especially sterol, composition of the membrane. The adsorption of certain electroneutral molecules known as dipole modifiers may also lead to significant changes in the magnitude of the potential decrease. These agents are widely used to study the effects of the dipole potential on membrane transport. This review presents a critical analysis of a variety of data from studies dedicated to ion channel formation and functioning in membranes with different dipole potentials. The types of ion channels found in cellular membranes and pores formed by antimicrobial agents and toxins in artificial lipid membranes are summarized. The mechanisms underlying the influence of the membrane dipole potential on ion channel activity, including dipole-dipole and charge-dipole interactions in the pores and in membranes, are discussed. A hypothesis, in which lipid rafts in both model and cellular membranes also modulate ion channel activity by virtue of an increased or decreased dipole potential, is also considered.

Interfacial water screens the protein-induced transmembrane voltage

Transmembrane proteins are crucial in cellular traffic, signal transduction, and energy storage in a form of transmembrane voltage. These proteins are stabilized by hydrophobic and hydrophilic interactions, in which cytoplasmic and exoplasmic water plays a special role. Water structural ordering generates the dipole potential that typically overcompensates for an intrinsic membrane-protein potential gradient, and thus it modifies and sustains an overall cellular electrostatics. Although the transmembrane voltage has been extensively studied, the dipole potential has attracted very little attention. Here, by using molecular dynamics, we examined water electrostatic response to the transmembrane charge, field, and potential asymmetry introduced by the presence of four integral membrane proteins: typical of inner (α-helix) and outer membrane (β-barrel). In all cases, the protein presence introduces electrostatic directionality in the transmembrane dipole field and voltage. In particular, water generates a deep potential sink if strongly polar residues are densely packed on one side of bilayer, as frequently occurs in a selectivity filter of the K(+) channel. We also found that protein secondary structure is less important than the polar residue distribution along the protein channel. Our findings are relevant for understanding the driving force behind biomembrane conductivity: the ability of biological water to electrostatically screen the transmembrane voltage.

The effects of solvent composition on the affinity of a peptide towards hair keratin: experimental and molecular dynamics data

Molecular dynamics simulations with a developed hair protofibril model demonstrated the ability to improve peptide uptake by hair shafts.

Cholesterol behavior in asymmetric lipid bilayers: insights from molecular dynamics simulations

Asymmetric lipid composition of the cell membranes plays an important role in the multitude of important biological functions. Much less is known, however, about the distribution and dynamics of cholesterol in asymmetric biological membranes. In this work we show how this issue could be addressed computationally by molecular dynamics simulations. The influence of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and the surface tension of the membrane on cholesterol distribution in asymmetric lipid bilayers is investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC, or DSPS lipids, were simulated on the time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with distribution ratio ~0.4-0.6. Multiple flip-flop transitions of cholesterol were observed directly and their mean times range from 350 to 2,000 ns. It was shown that the distribution of cholesterol in the asymmetric bilayer depends not only on the type of lipids but also on the local membrane curvature and the surface tension. The geometric shape of spontaneously curved asymmetric bilayer changes dramatically in the presence of cholesterol. The membrane curvature becomes less homogeneous with large patches of flattened regions interleaved by rather sharp bends.

Bringing bioelectricity to light

Any bilayer lipid membrane can support a membrane voltage. The combination of optical perturbation and optical readout of membrane voltage opens the door to studies of electrophysiology in a huge variety of systems previously inaccessible to electrode-based measurements. Yet, the application of optogenetic electrophysiology requires careful reconsideration of the fundamentals of bioelectricity. Rules of thumb appropriate for neuroscience and cardiology may not apply in systems with dramatically different sizes, lipid compositions, charge carriers, or protein machinery. Optogenetic tools are not electrodes; thus, optical and electrode-based measurements have different quirks. Here we review the fundamental aspects of bioelectricity with the aim of laying a conceptual framework for all-optical electrophysiology.

Lattice simulations of phase morphology on lipid bilayers: renormalization, membrane shape, and electrostatic dipole interactions

When liquid phases coexist at equilibrium but are not driven to minimize domain interfacial contact energy, the resulting patterns of phase domains can have important implications for living cells. In this study we explore some of the interactions and conditions that produce the stable patterned phases that are observed in model lipid mixtures. By use of Monte Carlo simulations we find that background curvature is important for the formation of patterned (modulated) phases. The interactions that stabilize nanoscopic phase separation are still not well understood. We show that inclusion of an electrostatic dipole repulsion with decay lengths as short as two to four lipid diameters can break up domains at the nanometer scale and that the location of the miscibility critical point is sensitive to this interaction. The use of a coarse-grained simulation raises questions about comparing parameters in simulations performed at different length scales. Using renormalization group techniques we show how to reconcile this problem, treating line tension as a running coupling constant.

The alteration of lipid bilayer dynamics by phloretin and 6-ketocholestanol

Lipid bilayer properties are quantified with a variety of arbitrary selected parameters such as molecular packing and dynamics, electrostatic potentials or permeability. In the paper we determined the effect of phloretin and 6-ketocholestanol (dipole potential modifying agents) on the membrane hydration and efficiency of the trans-membrane water flow. The dynamics of water molecules within the lipid bilayer interface was evaluated using solvent relaxation method, whereas the osmotically induced trans-membrane water flux was estimated with the stopped-flow method using the liposome shrinkage kinetics. The presence of phloretin or 6-ketocholestanol resulted in a change of both, the interfacial hydration level and osmotically driven water fluxes. Specifically, the presence of 6-ketocholestanol reduced the amount and mobility of water in the membrane interface. It also slows the osmotically induced water flow. The interfacial hydration change caused by phloretin was much smaller and the effect on osmotically induced water flow was opposite to that of 6-ketocholestanol.

A comparison of the electromechanical properties of structurally diverse proteins by molecular dynamics simulation

Proteins are subjected to electric fields both within the cell and during routine biochemical analysis. We have used atomistic molecular dynamics simulations to study conformational changes within three structurally diverse proteins subjected to high electric fields. At electric fields in excess of .5 V/nm, major structural changes were observed in all three proteins due to charge redistribution within the biomolecule. However, the electromechanical resilience was found to be highly dependent on the protein secondary structure, with α-helices showing a particularly high susceptibility to deformation by the applied electric field.

Effects of oriented surface dipole on photoconversion efficiency in an alkane/lipid-hybrid-bilayer-based photovoltaic model system

When a phospholipid monolayer containing a zinc-coordinated porphyrin species formed atop a self-assembled monolayer of heptadecafluoro-1-decanethiol (CF3(CF2)7(CH2)2SH) is subjected to photoelectrochemical current generation, a significant modulation effect is observed. Compared with devices that contain similar photoactive lipid monolayers but formed on 1-dodecanethiol SAMs, these fluorinated hybrid bilayers produce a >60% increase in cathodic currents and a similar decrease in anodic currents. Photovoltages recorded from these hybrid bilayers are found to vary in the same fashion. The modulation of photovoltaic responses in these hybrid-bilayer-based devices is explained by the opposite surface dipoles associated with the thiols employed in this study, which in one case (fluorothiol) increase and in another (alkanethiol) decrease the work function of the underlying gold substrates. A similar trend of photovoltage/photocurrent modulation is also observed if fullerene is used as the photoagent in these devices. Our results reveal the intricacy of orientated surface dipole in influencing the photovoltaic processes, and its subtle interplay with other factors related to the photoagents, such as their location and orientation within the organic matrix.

Excited-state proton coupled charge transfer modulated by molecular structure and media polarization

Charge and proton transfer reactions in the excited states of organic dyes can be coupled in many different ways. Despite the complementarity of charges, they can occur on different time scales and in different directions of the molecular framework. In certain cases, excited-state equilibrium can be established between the charge-transfer and proton-transfer species. The interplay of these reactions can be modulated and even reversed by variations in dye molecular structures and changes of the surrounding media. With knowledge of the mechanisms of these processes, desired rates and directions can be achieved, and thus the multiple emission spectral features can be harnessed. These features have found versatile applications in a number of cutting-edge technological areas, particularly in fluorescence sensing and imaging.

Cholesterol induces uneven curvature of asymmetric lipid bilayers

A remarkable flexibility is observed in biological membranes, which allows them to form the structures of different curvatures. We addressed the question of intrinsic ability of phospholipid membranes to form highly curved structures and the role of cholesterol in this process. The distribution of cholesterol in the highly curved asymmetric DOPC/DOPS lipid bilayer was investigated by the coarse-grained molecular dynamics simulations in the membrane patches with large aspect ratio. It is shown that cholesterol induces uneven membrane curvature promoting the formation of extended flattened regions of the membrane interleaved by sharp bends. It is shown that the affinity of cholesterol to anionic DOPS or neutral DOPC lipids is curvature dependent. The cholesterol prefers DOPS to DOPC in either planar or highly curved parts of the membrane. In contrast, in the narrow interval of moderate membrane curvatures this preference is inverted. Our data suggest that there is a complex self-consistent interplay between the membrane curvature and cholesterol distribution in the asymmetric lipid bilayers. The suggested new function of cholesterol may have a biological relevance.

Will C-Laurdan dethrone Laurdan in fluorescent solvent relaxation techniques for lipid membrane studies?

Development of fluorescence methods involves the necessity of understanding the fluorescent probes behavior in their ground and excited states. In the case of biological membranes, the position of the dyes in the lipid bilayer and their response after excitation are necessary parameters to correctly analyze the experiments. In the present work, we focus on two fluorescent markers, Laurdan (6-lauroyl-2-(N,N-dimethylamino)naphthalene) and its derivative C-Laurdan (6-dodecanoyl-2-[N-methyl-N-(carboxymethyl)amino]naphthalene), recently proposed for lipid raft visualization [Kim, H. M.; et al. ChemBioChem 2007, 8, 553]. C-Laurdan, by the presence of an additional carboxyl group, has an advantage over Laurdan since it has a higher sensitivity to the membrane polarity at the lipid headgroup region and a higher water solubility. This theoretical study, based on quantum mechanical (QM) and molecular dynamics (MD) simulations in a fully hydrated lipid membrane model, compare the equilibrium and dynamic properties of both probes taking into account their fluorescence lifetimes. C-Laurdan is found to be more stable than Laurdan in the headgroup region of the membrane and also much more aligned with the lipids. This study suggests that, besides the lipid raft imaging, the C-Laurdan marker can considerably extend the capabilities of fluorescent solvent relaxation method.

Beyond annexin V: fluorescence response of cellular membranes to apoptosis

Dramatic changes in the structure of cell membranes on apoptosis allow easy, sensitive and non-destructive analysis of this process with the application of fluorescence methods. The strong plasma membrane asymmetry is present in living cells, and its loss on apoptosis is commonly detected with the probes interacting strongly and specifically with phosphatidylserine (PS). This phospholipid becomes exposed to the cell surface, and the application of annexin V labeled with fluorescent dye is presently the most popular tool for its detection. Several methods have been suggested recently that offer important advantages over annexin V assay with the ability to study apoptosis by spectroscopy of cell suspensions, flow cytometry and confocal or two-photon microscopy. The PS exposure marks the integrated changes in the outer leaflet of cell membrane that involve electrostatic potential and hydration, and the attempts are being made to provide direct probing of these changes. This review describes the basic mechanisms underlying the loss of membrane asymmetry during apoptosis and discusses, in comparison with the annexin V-binding assay, the novel fluorescence techniques of detecting apoptosis on cellular membrane level. In more detail we describe the detection method based on smart fluorescent dye F2N12S incorporated into outer leaflet of cell membrane and reporting on apoptotic cell transformation by easily detectable change of the spectral distribution of fluorescent emission. It can be adapted to any assay format.

Differential effect of cholesterol and its biosynthetic precursors on membrane dipole potential

Dipole potential is the potential difference within the membrane bilayer, which originates due to the nonrandom arrangement of lipid dipoles and water molecules at the membrane interface. Cholesterol, a representative sterol in higher eukaryotic membranes, is known to increase membrane dipole potential. In this work, we explored the effects of immediate (7-DHC and desmosterol) and evolutionary (ergosterol) precursors of cholesterol on membrane dipole potential, monitored by the dual wavelength ratiometric approach utilizing the probe di-8-ANEPPS. Our results show that the effect of these precursors on membrane dipole potential is very different from that observed with cholesterol, although the structural differences among them are subtle. These results assume relevance, since accumulation of cholesterol precursors due to defective cholesterol biosynthesis has been reported to result in several inherited metabolic disorders such as the Smith-Lemli-Opitz syndrome. Interestingly, cholesterol (and its precursors) has a negligible effect on dipole potential in polyunsaturated membranes. We interpret these results in terms of noncanonical orientation of cholesterol in these membranes. Our results constitute the first report on the effect of biosynthetic and evolutionary precursors of cholesterol on dipole potential, and imply that a subtle change in sterol structure can significantly alter the dipolar field at the membrane interface.

Measurements and implications of the membrane dipole potential

There are three kinds of membrane potentials: the surface potentials, resulting from the accumulation of charges at the membrane surfaces; the transmembrane potential, determined by imbalance of charge in the aqueous solutions; and the dipole potential, a membrane-internal potential from the dipolar components of the phospholipids and interface water. The absolute value of the dipole potential has been very difficult to measure, although its value has been estimated to be in the range of 200-1,000 mV from ion translocation rates (determined by the planar lipid bilayer method), the surface potential of lipid monolayers (determined by the lipid monolayer method), molecular-dynamics calculations, and electron scattering using cryoelectron microscopy (cryo-EM). Spectroscopy methods have also been used to monitor the dipole potential changes on the basis of the observed fluorescence changes of voltage-sensitive probes. The dipole potential accounts for the much larger permeability of a bare phospholipid membrane to anions than cations and affects the conformation and function of membrane proteins.

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

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

Absorption and fluorescence of PRODAN in phospholipid bilayers: a combined quantum mechanics and classical molecular dynamics study

Absorption and fluorescence spectra of PRODAN (6-propionyl-2-dimethylaminonaphthalene) were studied by means of the time-dependent density functional theory and the algebraic diagrammatic construction method. The influence of environment, a phosphatidylcholine lipid bilayer and water, was taken into account employing a combination of quantum chemical calculations with empirical force-field molecular dynamics simulations. Additionally, experimental absorption and emission spectra of PRODAN were measured in cyclohexane, water, and lipid vesicles. Both planar and twisted configurations of the first excited state of PRODAN were taken into account. The twisted structure is stabilized in both water and a lipid bilayer, and should be considered as an emitting state in polar environments. Orientation of the excited dye in the lipid bilayer significantly depends on configuration. In the bilayer, the fluorescence spectrum can be regarded as a combination of emission from both planar and twisted structures.

Electrostatics of deformable lipid membranes

It was recently demonstrated that significant local deformations of biological membranes take place due to the fields of charged peptides and ions, challenging the standard model of membrane electrostatics. The ability of ions to retain their immediate hydration environment, combined with the lack of sensitivity of permeability to ion type or even ion pairs, led us to question the extent to which hydration energetics and electrostatics control membrane ion permeation. Using the arginine analog methyl-guanidinium as a test case, we find that although hydrocarbon electronic polarizability causes dramatic changes in ion solvation free energy, as well as a significant change (approximately 0.4 V) in the membrane dipole potential, little change in membrane permeation energetics occurs. We attribute this to compensation of solvation terms from polar and polarizable nonpolar components within the membrane, and explain why the dipole potential is not fully sensed in terms of the locally deformed bilayer interface. Our descriptions provide a deeper understanding of the translocation process and allow predictions for poly-ions, ion pairs, charged lipids, and lipid flip-flop. We also report simulations of large hydrophobic-ion-like membrane defects and the ionophore valinomycin, which exhibit little membrane deformation, as well as hydrophilic defects and the ion channel gramicidin A, to provide parallels to membranes deformed by unassisted ion permeation.