Depth of water penetration into lipid bilayers

The Chemical Reactivity of Membrane Lipids

It is well-known that aqueous dispersions of phospholipids spontaneously assemble into bilayer structures. These structures have numerous applications across chemistry and materials science and form the fundamental structural unit of the biological membrane. The particular environment of the lipid bilayer, with a water-poor low dielectric core surrounded by a more polar and better hydrated interfacial region, gives the membrane particular biophysical and physicochemical properties and presents a unique environment for chemical reactions to occur. Many different types of molecule spanning a range of sizes, from dissolved gases through small organics to proteins, are able to interact with membranes and promote chemical changes to lipids that subsequently affect the physicochemical properties of the bilayer. This Review describes the chemical reactivity exhibited by lipids in their membrane form, with an emphasis on conditions where the lipids are well hydrated in the form of bilayers. Key topics include the following: lytic reactions of glyceryl esters, including hydrolysis, aminolysis, and transesterification; oxidation reactions of alkenes in unsaturated fatty acids and sterols, including autoxidation and oxidation by singlet oxygen; reactivity of headgroups, particularly with reactive carbonyl species; and E/Z isomerization of alkenes. The consequences of reactivity for biological activity and biophysical properties are also discussed.

The Design and Optimization of Ceramide NP-Loaded Liposomes to Restore the Skin Barrier

The impairment of skin integrity derived from derangement of the orthorhombic lateral organization is mainly caused by dysregulation of ceramide amounts in the skin barrier. Ceramides, fatty acids, and cholesterol-containing nano-based formulations have been used to impair the skin barrier. However, there is still a challenge to formulate novel formulations consisting of ceramides due to their chemical structure, poor aqueous solubility, and high molecular weight. In this study, the design and optimization of Ceramide 3 (CER-NP)-loaded liposomes are implemented based on response surface methodology (RSM). The optimum CER-NP-loaded liposome was selected based on its particle size (PS) and polydispersity index (PDI). The optimum CER-NP-loaded liposome was imagined by observing the encapsulation by using a confocal laser scanning microscope (CLSM) within fluorescently labeled CER-NP. The characteristic liquid crystalline phase and lipid chain conformation of CER-NP-loaded liposomes were determined using attenuated total reflectance infrared spectroscopy (ATR-IR). The CER-NP-loaded liposomes were imagined using a field emission scanning electron microscope (FE-SEM). Finally, the in vitro release of CER-NP from liposomes was examined using modified Franz Cells. The experimental and predicted results were well correlated. The CLSM images of optimized liposomes were conformable with the other studies, and the encapsulation efficiency of CER-NP was 93.84 ± 0.87%. ATR-IR analysis supported the characteristics of the CER-NP-loaded liposome. In addition, the lipid chain conformation shows similarity with skin barrier lipid organization. The release pattern of CER-NP liposomes was fitted with the Korsmeyer-Peppas model. The cytotoxicity studies carried out on HaCaT keratinocytes supported the idea that the liposomes for topical administration of CER-NP could be considered relatively safe. In conclusion, the optimized CER-NP-loaded liposomes could have the potential to restore the skin barrier function.

Effects of sugar concentration on the electroporation, size distribution and average size of charged giant unilamellar vesicles

We investigated the effects of sugar concentration on the electroporation, size distribution and average size of giant unilamellar vesicles (GUVs). GUVs were prepared from 40 mol% of 1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DOPG) and 60 mol% of 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipids. Pulsed electric field was applied to the 40%DOPG/60%DOPC-GUVs and it induced lateral electric tension (σ_c) in the membranes of vesicles. The σ_c-induced probability of rupture (P_pore) and the rate constant of rupture (k_p) of GUVs under the sugar concentration, c = 40, 100 and 300 mM, were determined. Both the P_pore and k_p increased with the increase of σ_c, but higher tension was required to generate the same values of P_pore and k_p with increasing c. We also investigated average sizes of GUVs from the size distribution of vesicles under various sugar concentrations. With the increase of c, the peak of the size distribution histograms shifted to the region of smaller vesicles. The average size decreased 1.6-fold when c increased from 10 to 300 mM. These investigations help to understand various biomedical, biophysical, and biochemical processes in vesicles and cells. Electroporation, size distribution and average size of charged GUVs were investigated under various sugar concentrations. The sugar concentration influences the electroporation of vesicles and the average size of GUVs.

Surface Characterization of Lipid Biomimetic Systems

Zeta potential and dipole potential measures are direct operational methodologies to determine the adsorption, insertion and penetration of ions, amphipathic and neutral compounds into the membranes of cells and model systems. From these results, the contribution of charged and dipole groups can be deduced. However, although each method may give apparent affinity or binding constants, care should be taken to interpret them in terms of physical meaning because they are not independent properties. On the base of a recent model in which the lipid bilayer is considered as composed by two interphase regions at each side of the hydrocarbon core, this review describes how dipole potential and zeta potential are correlated due to water reorganization. From this analysis, considering that in a cell the interphase region the membrane extends to the cell interior or overlaps with the interphase region of another supramolecular structure, the correlation of dipole and electrostatic forces can be taken as responsible of the propagation of perturbations between membrane and cytoplasm and vice versa. Thus, this picture gives the membrane a responsive character in addition to that of a selective permeability barrier when integrated to a complex system.

Effects of osmotic pressure on the irreversible electroporation in giant lipid vesicles

Irreversible electroporation (IRE) is a nonthermal tumor/cell ablation technique in which a series of high-voltage short pulses are used. As a new approach, we aimed to investigate the rupture of giant unilamellar vesicles (GUVs) using the IRE technique under different osmotic pressures (Π), and estimated the membrane tension due to Π. Two categories of GUVs were used in this study. One was prepared with a mixture of dioleoylphosphatidylglycerol (DOPG), dioleoylphosphatidylcholine (DOPC) and cholesterol (chol) for obtaining more biological relevance while other with a mixture of DOPG and DOPC, with specific molar ratios. We determined the rate constant (kp) of rupture of DOPG/DOPC/chol (46/39/15)-GUVs and DOPG/DOPC (40/60)-GUVs induced by constant electric tension (σc) under different Π. The σc dependent kp values were fitted with a theoretical equation, and the corresponding membrane tension (σoseq) at swelling equilibrium under Π was estimated. The estimated membrane tension agreed well with the theoretical calculation within the experimental error. Interestingly, the values of σoseq were almost same for both types of synthesized GUVs under same osmotic pressure. We also examined the sucrose leakage, due to large osmotic pressure-induced pore formation, from the inside of DOPG/DOPC/chol(46/39/15)-GUVs. The estimated membrane tension due to large Π at which sucrose leaked out was very similar to the electric tension at which GUVs were ruptured without Π. We explained the σc and Π induced pore formation in the lipid membranes of GUVs.

An investigation into the critical tension of electroporation in anionic lipid vesicles

Irreversible electroporation (IRE) is a technique for the disruption of localized cells or vesicles by a series of short and high-frequency electric pulses which has been used for tissue ablation and treatment in certain diseases. It is well reported that IRE induces lateral tension in the membranes of giant unilamellar vesicles (GUVs). The GUVs are prepared by a mixture of anionic lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. Here the influence of DOPG mole fraction, XDOPG, on the critical tension of electroporation in GUVs has been investigated in sodium chloride-containing PIPES buffer. The critical tension decreases from 9.0 ± 0.3 to 6.0 ± 0.2 mN/m with the increase of XDOPG from 0.0 to 0.60 in the membranes of GUVs. Hence an increase in XDOPG greatly decreases the mechanical stability of membranes. We develop a theoretical equation that fits the XDOPG dependent normalized critical tension, and obtain a binding constant for the lipid-ion interaction of 0.75 M-1. The decrease in the energy barrier for formation of the nano-size nascent or prepore state, due to the increase in XDOPG, is the main factor explaining the decrease in critical tension of electroporation in vesicles.

Electrostatic effects on the electrical tension-induced irreversible pore formation in giant unilamellar vesicles

Irreversible electroporation (IRE) is a new technique in which a series of short pulses with high frequency electrical energy is applied on the targeted regions of cells or vesicles for their destruction or rupture formation. IRE induces lateral tension in the membranes of vesicles. We have investigated the electrostatic interaction effects on the constant electrical tension-induced rate constant of irreversible pore formation in the membranes of giant unilamellar vesicles (GUVs). The electrostatic interaction has been varied by changing the salt concentration in buffer and the surface charge density of membranes. The membranes of GUVs are synthesized by a mixture of negatively charged lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. The rate constant of pore formation increases with the decrease of salt concentration in buffer along with the increase of surface charge density of membranes. The tension dependent probability of pore formation and the rate constant of pore formation are fitted to the theoretical equation, and obtained the line tension of membranes. The decrease in energy barrier of a prepore due to electrostatic interaction is the key factor causing an increase of rate constant of pore formation.

Influence of cholesterol on electroporation in lipid membranes of giant vesicles

Irreversible electroporation (IRE) is primarily a nonthermal ablative technology that uses a series of high-voltage and ultra-short pulses with high-frequency electrical energy to induce cell death. This paper presents the influence of cholesterol on the IRE-induced probability of pore formation and the rate constant of pore formation in giant unilamellar vesicles (GUVs). The GUVs are prepared by a mixture of dioleoylphosphatidylglycerol (DOPG), dioleoylphosphatidylcholine (DOPC) and cholesterol using the natural swelling method. An IRE signal of frequency 1.1 kHz is applied to the membranes of GUVs. The probability of pore formation and the rate constant of pore formation events are obtained using statistical analysis from several single GUVs. The time-dependent fraction of intact GUVs among all those examined is fitted to a single exponential decay function from where the rate constant of pore formation is calculated. The probability of pore formation and the rate constant of pore formation decreases with an increase in cholesterol content in the membranes of GUVs. Theoretical equations are fitted to the tension-dependent rate constant of pore formation and to the probability of pore formation, which allows us to obtain the line tension of membranes. The obtained line tension increases with an increase in cholesterol in the membranes. The increase in the energy barrier of the prepore state, due to the increase of cholesterol in membranes, is the main factor explaining the decrease in the rate constant of pore formation.

Effects of electrically-induced constant tension on giant unilamellar vesicles using irreversible electroporation

Stretching in membranes of cells and vesicles plays important roles in various physiological and physicochemical phenomena. Irreversible electroporation (IRE) is the irreversible permeabilization of the membrane through the application of a series of electrical field pulses of micro- to millisecond duration. IRE induces lateral tension due to stretching in the membranes of giant unilamellar vesicles (GUVs). However, the effects of electrically induced (i.e., IRE) constant tension in the membranes of GUVs have not been investigated yet in detail. To explore the effects of electrically induced tension on GUVs, firstly a microcontroller-based IRE technique is developed which produces electric field pulses (332 V/cm) with pulse width 200 µs. Then the electrodeformation, electrofusion and membrane rupture of GUVs are investigated at various constant tensions in which the membranes of GUVs are composed of dioleoylphosphatidylglycerol (DOPG) and dioleoylphosphatidylcholine (DOPC). Stochastic electropore formation is observed in the membranes at an electrically induced constant tension in which the probability of pore formation is increased with the increase of tension from 2.5 to 7.0 mN/m. The results of pore formation at different electrically-induced constant tensions are in agreement with those reported for mechanically-induced constant tension. The decrease in the energy barrier of the pre-pore state due to the increase of electrically-induced tension is the main factor increasing the probability of electropore formation. These investigations help to provide an understanding of the complex behavior of cells/vesicles in electric field pulses and can form the basis for practical applications in biomedical technology.

Topically applied lipid- and surfactant-based nanoparticles in the treatment of skin disorders

INTRODUCTION

In the treatment of dermatological disorders, topical drug administration is a mainstay. However, nanoparticle-based carrier systems could improve and expand the current therapeutic range via localized delivery of active ingredients. Areas covered: This review gives a detailed description of lipid- and surfactant-based drug delivery systems which have been explored for topical drug administration. To guide researchers in their choice of delivery system, an informative decision tree is included. Moreover, this review provides a complete overview of the topical or transdermal drug products, currently on the market or under clinical investigation, delivered via the discussed carriers, in the treatment of skin disorders. Expert opinion: Conventional liposomes are still popular in the domain of topical or transdermal drug delivery and dominate the market landscape. However, several other carriers, such as exosomes and niosomes, are being explored which offer distinct advantages over liposomes and should therefore not be disregarded when selecting a proper drug delivery system.

Computer modelling studies of the bilayer/water interface

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

An elastic liposomal formulation for RNAi-based topical treatment of skin disorders: Proof-of-concept in the treatment of psoriasis

RNA interference (RNAi) is a rapidly emerging approach for targeted gene silencing to alleviate disease pathology. However, lack of efficient carriers for targeted delivery delays the clinical translation of RNAi. An interesting target for local RNAi therapeutics is the skin as it allows direct access to target cells. Still, applications are limited due to the effective skin barrier which hinders penetration. Herein, a description is given of a liposomal carrier, called 'DDC642', capable of delivering RNAi molecules to the epidermis of impaired and intact human skin, without targeting the dermis or circulatory system. In a psoriasis tissue model, down-regulation of the psoriasis marker human beta-defensin 2 by DDC642-delivered siRNA was confirmed, providing proof-of-concept. These liposomes thus hold great potential as topical delivery system for RNAi therapeutics in the treatment of numerous skin diseases.

Oral Delivery of Bovine Lactoferrin Using Pectin- and Chitosan-Modified Liposomes and Solid Lipid Particles: Improvement of Stability of Lactoferrin

A critical problem associated with delivery of bovine lactoferrin (bLf) by the oral route is low bioavailability, which is derived from the enzymatic degradation in the gastrointestinal tract and poor permeation across the intestinal epitheliums. Particulate carrier systems have been identified to protect bLf against proteolysis via encapsulation. This study aimed to evaluate the physico-chemical stability of bLf-loaded liposomes and solid lipid particles (SLPs) modified by pectin and chitosan when exposed to various stress conditions. Transmission electron microscopy results showed liposomes and SLPs had a classic shell-core structure with polymer layers surrounded on surface, but the structure appeared to be partially broken after digestion in simulated intestinal fluid (SIF). Although HPLC and sodium dodecyl sulphate-polyacrylamide gel electrophoresis methods qualitatively and quantitatively described either liposomes or SLPs could retain intact bLf against proteolysis in SIF to some extent, all liposome formulations showed rapid rate of lipolysis mediated by pancreatic enzymes. On the other hand, all SLP formulations showed higher heat resistance and greater electrolyte tolerance compared to liposome formulations. After 180 days storage time, liposome-loaded bLf was completely degraded, whereas almost 30% of intact bLf still remained in SLP formulations. Overall, SLPs are considered as primary choice for oral bLf delivery.

Hydration in Lipid Monolayers: Correlation of Water Activity and Surface Pressure

In order to give a physical meaning to each region of the membrane we define the interphase as the region in a lipid membrane corresponding to the polar head groups imbibed in water with net different properties than the hydrocarbon region and the water phase. The interphase region is analyzed under the scope of thermodynamics of surface and solutions based on the definition of Defay-Prigogine of an interphase and the derivation that it has in the understanding of membrane processeses in the context of biological response. In the view of this approach, the complete monolayer is considered as the lipid layer one molecule thick plus the bidimensional solution of the polar head groups inherent to it (the interphase region). Surface water activity appears as a common factor for the interaction of several aqueous soluble and surface active proteins with lipid membranes of different composition. Protein perturbation can be measured by changes in the surface pressure of lipid monolayers at different initial water surface activities. As predicted by solution chemistry, the increase of surface pressure is independent of the particle nature that dissolves. Therefore, membranes give a similar response in terms of the determined surface states given by water activity independent of the protein or peptide.

Membrane Hydration: A Hint to a New Model for Biomembranes

The classical view of a biological membrane is based on the Singer-Nicholson mosaic fluid model in which the lipid bilayer is the structural backbone. Under this paradigm, many studies of biological processes such as, permeability, active transport, enzyme activity and adhesion and fusion processes have been rationalized considering the lipid membrane as a low dielectric slab of hydrocarbon chains with polar head groups exposed to water at each side in which oil/water partition prevails. In spite of several analyses and evidence available in relation to membrane hydration, water is not taken into account as a functional component. For this purpose, new insights in the water organization in restricted environments and the thermodynamical and mechanical properties emerging from them are specifically analysed and correlated.This chapter summarizes the progress of the studies of water in membranes along the book in order to give a more realistic structural and dynamical picture accounting for the membrane functional properties.

Structural and thermodynamic properties of water-membrane interphases: significance for peptide/membrane interactions

Water appears as a common intermediary in the mechanisms of interaction of proteins and polypeptides with membranes of different lipid composition. In this review, how water modulates the interaction of peptides and proteins with lipid membranes is discussed by correlating the thermodynamic response and the structural changes of water at the membrane interphases. The thermodynamic properties of the lipid-protein interaction are governed by changes in the water activity of monolayers of different lipid composition according to the lateral surface pressure. In this context, different water populations can be characterized below and above the phase transition temperature in relation to the CH₂ conformers' states in the acyl chains. According to water species present at the interphase, lipid membrane acts as a water state regulator, which determines the interfacial water domains in the surface. It is proposed that those domains are formed by the contact between lipids themselves and between lipids and the water phase, which are needed to trigger adsorption-insertion processes. The water domains are essential to maintain functional dynamical properties and are formed by water beyond the hydration shell of the lipid head groups. These confined water domains probably carries information in local units in relation to the lipid composition thus accounting for the link between lipidomics and aquaomics. The analysis of these results contributes to a new insight of the lipid bilayer as a non-autonomous, responsive (reactive) structure that correlates with the dynamical properties of a living system.

Water state and carbonyl distribution populations in confined regions of lipid bilayers observed by FTIR spectroscopy

It has been suggested that water in confined regions presents different properties than bulk water, mainly because of the changes in water population species that may be induced by the adjacent walls of different polarities in terms of hydrogen bond formation. In this context, it would be expected that lipids in the gel and the fluid states should offer different templates for water organization. The presence of water pockets or defects in lipid bilayers has been proposed to explain the insertion of charged and polar peptides and amino acids in membranes. In this work, we provide direct evidence by means of FTIR spectroscopy that water band profiles are changed whether lipids are in the solid state, in the gel state after heating and cooling across the phase transition, or in the fluid state. The different bands found in each case were assigned to different H-bonded water populations in agreement with the exposure of carbonyl groups.

Connected and isolated CH2 populations in acyl chains and its relation to pockets of confined water in lipid membranes as observed by FTIR spectrometry

Analysis of the band corresponding to the frequency of vibrational symmetric stretching mode of methylene groups in the lipid acyl chains and the bands of water below and above the phase transition of different lipids by Fourier transform infrared spectroscopy gives strong support to the formation of confined water pockets in between the lipid acyl chains. Our measures and analysis consolidate the mechanism early proposed by Traüble, in the sense that water is present in kinks formed by trans-gauche isomers along the hydrocarbon tails. The formation of these regions depends on the acyl lipid composition, which determines the presence of different populations of water species, characterized by its degree of H bond coordination in fluid saturated or unsaturated lipids. The free energy excess due to the reinforcement of the water structure along few water molecules in the adjacencies of exposed membrane residues near the phase transition is a reasonable base to explain the insertion and translocation of polar peptides and amino acid residues through the biomembrane on thermodynamic and structural grounds.

Effect of membrane tension on the physical properties of DOPC lipid bilayer membrane

Molecular dynamics simulations of a dioleoylphosphocholine (DOPC) lipid bilayer were performed to explore its mechanosensitivity. Variations in the bilayer properties, such as area per lipid, volume, thickness, hydration depth (HD), hydration thickness (HT), lateral diffusion coefficient, and changes in lipid structural order were computed in the membrane tension range 0 to 15dyn/cm. We determined that an increase in membrane tension results in a decrease in the bilayer thickness and HD of ~5% and ~5.7% respectively, whereas area per lipid, volume, and HT/HD increased by 6.8%, 2.4%, and 5% respectively. The changes in lipid conformation and orientation were characterized using orientational (S(2)) and deuterium (S(CD)) order parameters. Upon increase of membrane tension both order parameters indicated an increase in lipid disorder by 10-20%, mostly in the tail end region of the hydrophobic chains. The effect of membrane tension on lipid lateral diffusion in the DOPC bilayer was analyzed on three different time scales corresponding to inertial motion, anomalous diffusion and normal diffusion. The results showed that lateral diffusion of lipid molecules is anomalous in nature due to the non-exponential distribution of waiting times. The anomalous and normal diffusion coefficients increased by 20% and 52% when the membrane tension changed from 0 to 15dyn/cm, respectively. In conclusion, our studies showed that membrane tension causes relatively significant changes in the area per lipid, volume, polarity, membrane thickness, and fluidity of the membrane suggesting multiple mechanisms by which mechanical perturbation of the membrane could trigger mechanosensitive response in cells.

Derivation and systematic validation of a refined all-atom force field for phosphatidylcholine lipids

An all-atomistic force field (FF) has been developed for fully saturated phospholipids. The parametrization has been largely based on high-level ab initio calculations in order to keep the empirical input to a minimum. Parameters for the lipid chains have been developed based on knowledge about bulk alkane liquids, for which thermodynamic and dynamic data are excellently reproduced. The FFs ability to simulate lipid bilayers in the liquid crystalline phase in a tensionless ensemble was tested in simulations of three lipids: 1,2-diauroyl-sn-glycero-3-phospocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1,2-dipalmitoyl-sn-glycero-3-phospcholine (DPPC). Computed areas and volumes per lipid, and three different kinds of bilayer thicknesses, have been investigated. Most importantly NMR order parameters and scattering form factors agree in an excellent manner with experimental data under a range of temperatures. Further, the compatibility with the AMBER FF for biomolecules as well as the ability to simulate bilayers in gel phase was demonstrated. Overall, the FF presented here provides the important balance between the hydrophilic and hydrophobic forces present in lipid bilayers and therefore can be used for more complicated studies of realistic biological membranes with protein insertions.

ELECTROFUSION OF MODEL LIPID MEMBRANES VIEWED WITH HIGH TEMPORAL RESOLUTION

The interaction of electric fields with lipid membranes and cells has been extensively studied in the last decades. The phenomena of electroporation and electrofusion are of particular interest because of their widespread use in cell biology and biotechnology. Giant vesicles, being of cell size and convenient for microscopy observations, are the simplest model of the cell membrane. However, optical microscopy observation of effects caused by electric DC pulses on giant vesicles is difficult because of the short duration of the pulse. Recently this difficulty has been overcome in our lab. Using a digital camera with high temporal resolution, we were able to access vesicle fusion dynamics on a sub-millisecond time scale. In this report, we present some observations on electrodeformation and –poration of single vesicles followed by an extensive study on the electrofusion of vesicle couples. Finally, we suggest an attractive approach for creating multidomain vesicles using electrofusion and present some preliminary results on the effect of membrane stiffness on the fusion dynamics.

Effect of temperature on the structure of charged membranes

Interactions between charged and neutral self-assembled phospholipid membranes are well understood and take into account temperature dependence. Yet, the manner in which the structure of the membrane is affected by temperature was hardly studied. Here we study the effect of temperature on the thickness, area per lipid, and volume per lipid of charged membranes. Two types of membranes were studied: membranes composed of charged lipids and dipolar (neutral) membranes that adsorbed divalent cations and became charged. Small-angle X-ray scattering data demonstrate that the thickness of charged membranes decreases with temperature. Wide-angle X-ray scattering data show that the area per headgroup increases with temperature. Intrinsically charged membranes linearly thin with temperature, whereas neutral membranes that adsorb divalent ions and become charged show an exponential decrease of their thickness. The data indicate that, on average, the tails shorten as the temperature rises. We attribute this behavior to higher lipid tail entropy and to the weaker electrostatic screening of the charged headgroups, by their counterions, at elevated temperatures. The latter effect leads to stronger electrostatic repulsion between the charged headgroups that increases the area per headgroup and decreases the bilayer thickness.

Lipid bilayer molecular dynamics study of lipid-derived agonists of the putative cannabinoid receptor, GPR55

Both L-α-lysophosphatidylinositol (LPI) and 2-arachidonoyl-sn-glycero-3-phosphoinositol (2-AGPI) have been reported to activate the putative cannabinoid receptor, GPR55. Recent microsecond time-scale molecular dynamics (MD) simulations and isothiocyanate covalent labeling studies have suggested that a transmembrane helix 6/7 (TMH6/7) lipid pathway for ligand entry may be necessary for interaction with cannabinoid receptors. Because LPI and 2-AGPI are lipid-derived ligands, conformations that each assumes in the lipid bilayer are therefore likely important for their interaction with GPR55. We report here the results of 70 ns NAMD molecular dynamics (MD) simulations of LPI and of 2-AGPI in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). These simulations are compared with a 70 ns simulation of the cannabinoid CB1 receptor endogenous ligand, N-arachidonoylethanolamine (anandamide, AEA) in a POPC bilayer. These simulations revealed that (1) LPI and 2-AGPI sit much higher in the bilayer than AEA, with inositol headgroups that can at times be solvated completely by water; (2) the behavior of the acyl chains of AEA and 2-AGPI are similar in their flexibilities in the bilayer, while the acyl chain of LPI has reduced flexibility; and (3) both 2-AGPI and LPI can adopt a tilted headgroup orientation by hydrogen bonding to the phospholipid phosphate/glycerol groups or via intramolecular hydrogen bonding. This tilted head group conformation (which represents over 40% of the conformer population of LPI (42.2 ± 3.3%) and 2-AGPI (43.7 ± 1.4%)) may provide a low enough profile in the lipid bilayer for LPI and 2-AGPI to enter GPR55 via the putative TMH6/7 entry port.

Relaxation processes in the adsorption of surface layer proteins to lipid membranes

The present work evaluates the kinetics of the interaction of S-layer protein from Lactobacillus brevis with lipid monolayers by measuring the changes in the surface pressure as a function of time for different lipid compositions and at different lateral compressions. At high surface pressures, or at high cholesterol ratios, in which membrane rigidity and surface polarity are increased, the kinetics can be described by a pure diffusional process. At low pressures or in the absence of cholesterol, the kinetics of protein interaction can be interpreted as a consequence of a relaxation process of the membrane structure coupled to diffusion. As the less packed monolayers are more hydrated, the relaxation processes at low initial surface pressures could be ascribed to changes in water organization in the membrane. These observations denote that kinetic insertion of proteins can be modulated by components that modify the hydration state of the interface.

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.

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.

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.

Finite element analysis of microelectrotension of cell membranes

Electric fields can be focused by micropipette-based electrodes to induce stresses on cell membranes leading to tension and poration. To date, however, these membrane stress distributions have not been quantified. In this study, we determine membrane tension, stress, and strain distributions in the vicinity of a microelectrode using finite element analysis of a multiscale electro-mechanical model of pipette, media, membrane, actin cortex, and cytoplasm. Electric field forces are coupled to membranes using the Maxwell stress tensor and membrane electrocompression theory. Results suggest that micropipette electrodes provide a new non-contact method to deliver physiological stresses directly to membranes in a focused and controlled manner, thus providing the quantitative foundation for micreoelectrotension, a new technique for membrane mechanobiology.

Giant vesicles in electric fields

This review is dedicated to electric field effects on giant unilamellar vesicles, a cell-size membrane system. We summarize various types of behavior observed when vesicles are subjected either to weak AC fields at various frequency, or to strong DC pulses. Different processes such as electro-deformation, -poration and -fusion of giant vesicles are considered. We describe some recent developments, which allowed us to detect the dynamics of the vesicle response with a resolution below milliseconds for all of these processes. Novel aspects on electric field effects on vesicles in the gel phase are introduced.

Effects of atorvastatin and simvastatin on atrial plateau currents

Recent evidence has shown that the inhibitors of the 3-hydroxy-3-methylglutaryl coenzyme A reductase (statins) might exert antiarrhythmic effects both in experimental models and in humans. In this study we analyzed the effects of atorvastatin and simvastatin acid (SVA) on the currents responsible for the duration of the plateau of human atrial action potentials: hKv1.5, Kv4.3, and L-type Ca(2+) (I(Ca,L)). hKv1.5 and Kv4.3 currents were recorded in transfected Ltk(-) and Chinese hamster ovary cells, respectively, and I(Ca,L) in mouse ventricular myocytes, using whole-cell patch-clamp. Atorvastatin and SVA produced a concentration-dependent block of hKv1.5 channels (IC(50)=4.5+/-1.7 microM and 5.7+/-0.03 microM, respectively) and shifted the midpoint of the activation and inactivation curves to more negative potentials. Importantly, atorvastatin- and SVA-induced block was added to that produced by quinidine, a drug that blocks hKv1.5 channels by binding to their pore cavity. Atorvastatin and SVA blocked Kv4.3 channels in a concentration-dependent manner (IC(50)=13.9+/-3.6 nM and 7.0+/-0.8 microM, respectively). Both drugs accelerated the inactivation kinetics and shifted the inactivation curve to more negative potentials. SVA (10 nM), but not atorvastatin, also blocked I(Ca,L) producing a frequency-dependent block that, at 2 Hz, reached a 50.2+/-1.5%. As a consequence of these effects, at nanomolar concentrations, atorvastatin lengthened, whereas SVA shortened, the duration of mouse atrial action potentials. The results suggest that atorvastatin and SVA alter Kv1.5 and Kv4.3 channel activity following a complex mechanism that does not imply the binding of the drug to the channel pore.

Structure and dynamics of phospholipid bilayers using recently developed general all-atom force fields

Two fully hydrated pure-species phospholipids bilayers, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dioleoyl-sn-glycero-3-phosphorylcholine (DOPC), in the fluid phase and explicit solvent have been studied using molecular dynamics simulation. Atom interactions were modeled using recently developed force fields based on AMBER with full atomistic details. Several representative liquid phase properties for the structure and dynamics of lipids with different length of hydrocarbon chains and different level of saturation have been reproduced without artificially biasing the system in order to match experimental data. In particular, as the new GAFF (General Amber Force Field) has not been explicitly developed to reproduce lipid characteristics and is naturally compatible with standard AMBER nucleic acids and proteins parameters, it is here proven a promising tool to study mixed lipid-protein processes as protein activity dependence on membrane composition, permeation of solute across membranes, and other cellular processes.

Electric pulses induce cylindrical deformations on giant vesicles in salt solutions

In this article, we report for the first time unusual shape changes of vesicles subjected to strong electric pulses in salt solutions of low concentration. The electric field is created by two parallel electrodes between which the vesicle solution is located. Surprisingly, the vesicles assume cylindrical shapes during the pulse. These deformations are short-lived (their lifetime is approximately 1 ms) and occur only in the presence of salt outside the vesicles, irrespective of their content. When the solution conductivities inside and outside are the same, vesicles with square cross section are observed. Using a fast digital camera, we were able to record these deformations and study the vesicle shape dynamics. The aims of this article are to report the new vesicle morphologies and their dynamics and to provoke theoretical work in this direction.

A practical guide to giant vesicles. Probing the membrane nanoregime via optical microscopy

Research on giant vesicles is becoming increasingly popular. Giant vesicles provide model biomembrane systems for systematic measurements of mechanical and rheological properties of bilayers as a function of membrane composition and temperature, as well as hydrodynamic interactions. Membrane response to external factors (for example electric fields, ions and amphiphilic molecules) can be directly visualized under the microscope. In this paper we review our current understanding of lipid bilayers as obtained from studies on giant unilamellar vesicles. Because research on giant vesicles increasingly attracts the interest of scientists from various backgrounds, we also try to provide a concise introduction for newcomers in the field. Finally, we summarize some recent developments on curvature effects induced by polymers, domain formation in membranes and shape transitions induced by electric fields.

Cholesterol effect on the dipole potential of lipid membranes

The effect of cholesterol removal by methyl-beta-cyclodextrin on the dipole potential, psi(d), of membrane vesicles composed of natural membrane lipids extracted from the kidney and brain of eight vertebrate species was investigated using the voltage-sensitive fluorescent probe di-8-ANEPPS. Cyclodextrin treatment reduced cholesterol levels by on average 80% and this was associated with an average reduction in psi(d) of 50 mV. Measurements of the effect of a range of cholesterol derivatives on the psi(d) of DMPC lipid vesicles showed that the magnitude of the effect correlated with the component of the sterol's dipole moment perpendicular to the membrane surface. The changes in psi(d) observed could not be accounted for solely by the electric field originating from the sterols' dipole moments. Additional factors must arise from sterol-induced changes in lipid packing, which changes the density of dipoles in the membrane, and changes in water penetration into the membrane, which changes the effective dielectric constant of the interfacial region. In DMPC membranes, the cholesterol-induced change in psi(d) was biphasic, i.e., a maximum in psi(d) was observed at approximately 35-45 mol %, after which psi(d) started to decrease. We suggest that this could be associated with a maximum in the strength of DMPC-cholesterol intermolecular forces at this composition.

Ion permeation through a narrow channel: using gramicidin to ascertain all-atom molecular dynamics potential of mean force methodology and biomolecular force fields

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.

Electro-deformation and poration of giant vesicles viewed with high temporal resolution

Fast digital imaging was used to study the deformation and poration of giant unilamellar vesicles subjected to electric pulses. For the first time the dynamics of response and relaxation of the membrane at micron-scale level is revealed at a time resolution of 30 micros. Above a critical transmembrane potential the lipid bilayer ruptures. Formation of macropores (diameter approximately 2 microm) with pore lifetime of approximately 10 ms has been detected. The pore lifetime has been interpreted as interplay between the pore edge tension and the membrane viscosity. The reported data, covering six decades of time, show the following regimes in the relaxation dynamics of the membrane. Tensed vesicles first relax to release the acquired stress due to stretching, approximately 100 micros. In the case of poration, membrane resealing occurs with a characteristic time of approximately 10 ms. Finally, for vesicles with excess area an additional slow regime was observed, approximately 1 s, which we associate with relaxation of membrane curvature. Dimensional analysis can reasonably well explain the corresponding characteristic timescales. Being performed on cell-sized giant unilamellar vesicles, this study brings insight to cell electroporation. The latter is widely used for gene transfection and drug transport across the membrane where processes occurring at different timescales may influence the efficiency.