Permeability of lipid bilayers to water and ionic solutes

The lipid bilayer moiety of biological membranes is considered to be the primary barrier to free diffusion of water and solutes. This conclusion arises from observations of lipid bilayer model membrane systems, which are generally less permeable than biological membranes. However, the nature of the permeability barrier remains unclear, particularly with respect to ionic solutes. For instance, anion permeability is significantly greater than cation permeability, and permeability to proton-hydroxide is orders of magnitude greater than other monovalent inorganic ions. In this review, we first consider bilayer permeability to water and discuss proposed permeation mechanisms which involve transient defects arising from thermal fluctuations. We next consider whether such defects can account for ion permeation, including proton-hydroxide flux. We conclude that at least two varieties of transient defects are required to explain permeation of water and ionic solutes.     

Electric Field-Driven Water Dipoles: Nanoscale Architecture of Electroporation

Electroporation is the formation of permeabilizing structures in the cell membrane under the influence of an externally imposed electric field. The resulting increased permeability of the membrane enables a wide range of biological applications, including the delivery of normally excluded substances into cells. While electroporation is used extensively in biology, biotechnology, and medicine, its molecular mechanism is not well understood. This lack of knowledge limits the ability to control and fine-tune the process. In this article we propose a novel molecular mechanism for the electroporation of a lipid bilayer based on energetics analysis. Using molecular dynamics simulations we demonstrate that pore formation is driven by the reorganization of the interfacial water molecules. Our energetics analysis and comparisons of simulations with and without the lipid bilayer show that the process of poration is driven by field-induced reorganization of water dipoles at the water-lipid or water-vacuum interfaces into more energetically favorable configurations, with their molecular dipoles oriented in the external field. Although the contributing role of water in electroporation has been noted previously, here we propose that interfacial water molecules are the main players in the process, its initiators and drivers. The role of the lipid layer, to a first-order approximation, is then reduced to a relatively passive barrier. This new view of electroporation simplifies the study of the problem, and opens up new opportunities in both theoretical modeling of the process and experimental research to better control or to use it in new, innovative ways.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions

Anionic palmitoyloleoylphosphatidylglycerol (POPG) is one of the most abundant lipids in nature, yet its atomic-scale properties have not received significant attention. Here we report extensive 150-ns molecular dynamics simulations of a pure POPG lipid membrane with sodium counterions. It turns out that the average area per lipid of the POPG bilayer under physiological conditions is approximately 19% smaller than that of a bilayer built from its zwitterionic phosphatidylcholine analog, palmitoyloleoylphosphatidylcholine. This suggests that there are strong attractive interactions between anionic POPG lipids, which overcome the electrostatic repulsion between negative charges of PG headgroups. We demonstrate that interlipid counterion bridges and strong intra- and intermolecular hydrogen bonding play a key role in this seemingly counterintuitive behavior. In particular, the substantial strength and stability of ion-mediated binding between anionic lipid headgroups leads to complexation of PG molecules and ions and formation of large PG-ion clusters that act in a concerted manner. The ion-mediated binding seems to provide a possible molecular-level explanation for the low permeability of PG-containing bacterial membranes to organic solvents: highly polar interactions at the water/membrane interface are able to create a high free energy barrier for hydrophobic molecules such as benzene.     

On the large error introduced in the estimate of the density of membrane pores from permeability measurements when diffusion in "unstirred layer" within the cells is disregarded

After the development of the "black lipid membrane" techniques, studies of the permeability of labeled water and nonelectrolytes across these artificial membranes have yielded permeability constants comparable in magnitude to those obtained from tracer studies of living cell membranes. This general agreement has affirmed the belief that the living cell membranes are indeed closely similar to these bilayer phospholipid membranes. In this report, we draw attention to a hidden assumption behind such comparisons made: the assumption that labeled material passing through the cell membrane barriers instantly reaches diffusion equilibrium inside the cell. The permeability constants to labeled water (and nonelectrolytes) across lipid layers were obtained using setups in which the lipid membrane was sandwiched between aqueous compartments both of which were vigorously stirred. In studies of permeability of living cell membranes only the outside solution was stirred, the intracellular water remained stationary. Yet the calculations of permeability constants of the cell membrane were made with the tacit assumption, that once the labeled materials pass through the cell membrane, they were instantly mixed with the entire cell contents as if a stirrer operating at infinite speed had been present inside the cells. Ignoring this unstirred condition of the intracellular water, in fact, lumped all the real-life delay due to diffusion in the cytoplasm and added it to the resistance to diffusion of the membrane barrier. The result is an estimated membrane permeability to labeled water (and nonelectrolytes) many times slower than it actually is. The present report begins with a detailed analysis of a specific case: tritiated water diffusion from giant barnacle muscle fibers and two non-living models, one real, one imagined.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effect of membrane tension on the electric field and dipole potential of lipid bilayer membrane

The dipole potential of lipid bilayer membrane controls the difference in permeability of the membrane to oppositely charged ions. We have combined molecular dynamics (MD) simulations and experimental studies to determine changes in electric field and electrostatic potential of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer in response to applied membrane tension. MD simulations based on CHARMM36 force field showed that electrostatic potential of DOPC bilayer decreases by ~45mV in the physiologically relevant range of membrane tension values (0 to 15dyn/cm). The electrostatic field exhibits a peak (~0.8×10(9)V/m) near the water/lipid interface which shifts by 0.9? towards the bilayer center at 15dyn/cm. Maximum membrane tension of 15dyn/cm caused 6.4% increase in area per lipid, 4.7% decrease in bilayer thickness and 1.4% increase in the volume of the bilayer. Dipole-potential sensitive fluorescent probes were used to detect membrane tension induced changes in DOPC vesicles exposed to osmotic stress. Experiments confirmed that dipole potential of DOPC bilayer decreases at higher membrane tensions. These results are suggestive of a potentially new mechanosensing mechanism by which mechanically induced structural changes in the lipid bilayer membrane could modulate the function of membrane proteins by altering electrostatic interactions and energetics of protein conformational states.     

An NMR study of the temperature dependence of the coefficient of water self-diffusion through lipid bilayer membranes

The temperature dependence of the coefficient of water self-diffusion through plane-parallel lipid multilayers of the phospholipid dioleoylphosphatidylcholine oriented on a glass support has been studied in the temperature range of 20-60 degrees C by the method of NMR with magnetic field pulse gradient. The values of the coefficients of transbilayer water diffusion are by four orders of magnitude less than for bulky water and ten times less than the coefficients of lateral diffusion of the lipid under the same conditions. The temperature dependence of the coefficient of water diffusion is described by the Arrhenius law with an apparent activation energy of about 41 kJ/mol, which far exceeds the activation energy for the diffusion of bulky water (18 kJ/mol). The experimental data were analyzed using a "dissolving-diffusion" model, by simulating the passage of water through membrane channels, and by analyzing the exchange of water molecules in states with different modes of translation mobility, including pore channels and bilayer "defects". Each of the approaches used made it possible to take the significance of bilayer permeability for the apparent energy of activation of water diffusion into account and estimate the energies of activation of water diffusion in the hydrophobic moiety of the bilayer, which were found to be close to the values for bulky water. The coefficients of water diffusion in the system under examination and the coefficients of permeation of water through the bilayer were estimated, and the effect of bilayer "defects" on the coefficients of water diffusion along and across bilayers was studied.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

Contribution of aqueous interphases to the permeability barrier of lipid bilayers for non-electrolytes

Equilibrium dialysis experiments are used to measure excluded volumes for the non-electrolyte permeant [U-¹⁴C] erythritol in lipid bilayer systems. The data indicate amounts of water associated with the lipid membranes which correspond with amounts calculated from calorimetric measurements.The membrane systems can be described as composite elements consisting of the lipid bilayers and adjacent water layers on both sides. The finding that the permeant is excluded indicates that the water layers contribute to the permeability barrier.The mean thickness of the water layers is about 6 Å for planar bilayers in multilayered liposomes and 10 Å for curved bilayers in sonicated vesicles. Next to the difference in thickness of the water layers differences in interfacial adsorption between the two systems are apparent.     

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells’ membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane’s mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.     

Temperature dependence of water self-diffusion through lipid bilayers assessed by NMR

The temperature dependence of the coefficient of water self-diffusion across plane-parallel multib-ilayers of dioleoylphosphatidylcholine oriented on a glass support was studied in the 20–60°C range by pulsed field gradient NMR. The coefficient for transbilayer diffusion of water proved almost four orders of magnitude smaller than for bulk water, and 10 times smaller than that for lateral diffusion of lipid under the same conditions. The temperature dependence obeyed the Arrhenius law with apparent activation energy of 41 kJ/mol, much higher than that for bulk water (18 kJ/mol). The experimental data were analyzed using the “dissolution-diffusion” model, by simulating water passage through membrane channels, and by examining water exchange in states with different modes of translational mobility, including pore channels and bilayer defects. Each approach could take into account the role of bilayer permeability and assess the apparent activation energy for water diffusion in the hydrophobic part of the bilayer, which proved close to the value for bulk water. Estimates were obtained for water diffusion coefficients in the system, coefficients of bilayer permeability for water, and the influence of bilayer defects on the lateral and transverse diffusion coefficients.     

Exploring gas permeability of cellular membranes and membrane channels with molecular dynamics

Aquaporins are a family of membrane proteins specialized in rapid water conduction across biological membranes. Whether these channels also conduct gas molecules and the physiological significance of this potential function have not been well understood. Here we report 140 ns of molecular dynamics simulations of membrane-embedded AQP1 and of a pure POPE bilayer addressing these questions. The permeability of AQP1 to two types of gas molecules, O2 and CO2, was investigated using two complementary methods, namely, explicit gas diffusion simulation and implicit ligand sampling. The simulations show that the central (tetrameric) pore of AQP1 can be readily used by either gas molecule to permeate the channel. The two approaches produced similar free energy profiles associated with gas permeation through the central pore: a -0.4 to -1.7 kcal/mol energy well in the middle, and a 3.6-4.6 kcal/mol energy barrier in the periplasmic vestibule. The barrier appears to be mainly due to a dense cluster of water molecules anchored in the periplasmic mouth of the central pore by four aspartate residues. Water pores show a very low permeability to O2, but may contribute to the overall permeation of CO2 due to its more hydrophilic nature. Although the central pore of AQP1 is found to be gas permeable, the pure POPE bilayer provides a much larger cross-sectional area, thus exhibiting a much lower free energy barrier for CO2 and O2 permeation. As such, gas conduction through AQP1 may only be physiologically relevant either in membranes of low gas permeability, or in cells where a major fraction of the cellular membrane is occupied by AQPs.     

Peptide Aggregation and Pore Formation in a Lipid Bilayer: A Combined Coarse-Grained and All Atom Molecular Dynamics Study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-μs CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 μs is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the α-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.     

Direct membrane effects of morphine and endorphins on Amoeba proteus.

Morphine, leu-enkephalinamide, met-enkephalin, alpha-neoendorphin and its Arg8 1-8 fragment increase contractile vacuole output in the freshwater Amoeba proteus at 18 microM. Significant effects of leu-enkephalin and naloxone are obtained at 180 microM. All compounds have reached their maximal activity at 720 microM. Alpha-neoendorphin and leu-enkephalin are inactive in the presence of isotonic, non-penetration sucrose, hence these compounds increase plasma membrane permeability to water. Results from molecular modeling show a clear correlation of activity with amphiphilicity, charge distribution and general flexibility of molecules. We conclude that, like previously-studied vasopressin analogues and non-hormonal amphiphilic peptides, active opioids embed themselves into the Amoeba plasma membrane, disrupting the lipid bilayer and increasing its permeability. In our Amoeba system, naloxone, a general morphine-like inhibitor, blocks active opioids as well as a vasopressin analogue. Naloxone, being less active than other tested amphiphiles, acts as a membrane stabilizer, protecting the lipid bilayer against the disruption action of more active compounds.     

Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map

alpha-Hemolysin of Staphylococcus aureus is a self-assembling toxin that forms a water-filled transmembrane channel upon oligomerization in a lipid membrane. Apart from being one of the best-studied toxins of bacterial origin, alpha-hemolysin is the principal component in several biotechnological applications, including systems for controlled delivery of small solutes across lipid membranes, stochastic sensors for small solutes, and an alternative to conventional technology for DNA sequencing. Through large-scale molecular dynamics simulations, we studied the permeability of the alpha-hemolysin/lipid bilayer complex for water and ions. The studied system, composed of approximately 300,000 atoms, included one copy of the protein, a patch of a DPPC lipid bilayer, and a 1 M water solution of KCl. Monitoring the fluctuations of the pore structure revealed an asymmetric, on average, cross section of the alpha-hemolysin stem. Applying external electrostatic fields produced a transmembrane ionic current; repeating simulations at several voltage biases yielded a current/voltage curve of alpha-hemolysin and a set of electrostatic potential maps. The selectivity of alpha-hemolysin to Cl(-) was found to depend on the direction and the magnitude of the applied voltage bias. The results of our simulations are in excellent quantitative agreement with available experimental data. Analyzing trajectories of all water molecule, we computed the alpha-hemolysin's osmotic permeability for water as well as its electroosmotic effect, and characterized the permeability of its seven side channels. The side channels were found to connect seven His-144 residues surrounding the stem of the protein to the bulk solution; the protonation of these residues was observed to affect the ion conductance, suggesting the seven His-144 to comprise the pH sensor that gates conductance of the alpha-hemolysin channel.     

Effect of ionic polarizability on electrodiffusion in lipid bilayer membranes.

Ion-carrier complexes and organic ions of similar size and shape have mobilities in lipid bilayer membranes which span several orders of magnitude. In this communication, an examination is made of the hypothesis that the basis for this unusually wide range of ionic mobilities is the potential energy barrier arising from image forces which selectively act on ions according to their polarizability. Using Poisson's equation to evaluate the electrostatic interaction between an ion and its surroundings, the potential energy barrier to ion transport due to image effects is computed, with the result that the potential energy barrier height depends strongly on ionic polarizability. Theoretical membrane potential energy profile calculations are used in conjunction with Nernst-Planck electrodiffusion equation to analyze the available mobility data for several ion-carrier complexes and lipid-soluble ions in lipid bilayer membranes. The variation among the mobilities of different ions is shown to be in agreement with theoretical predictions based on ionic polarizability and size. Furthermore, the important influence exerted by image forces on ion transport in lipid bilayer membranes compared to the frictional effect of membrane viscosity is established by contrasting available data on the activation energy of ionic conductivity with that for membrane fluidity.     

Molecular dynamics simulation studies of lipid bilayer systems

The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study required the construction of a model bacterial membrane which consisted of two types of phospholipids and counter ions. Whenever published experimental data were available, the results of the simulations were compared with them.     

Does CO2 permeate through aquaporin-1?

Aquaporins facilitate water permeation across biological membranes. Additionally, glycerol and other small neutral solutes are permeated by related aquaglyceroporins. The role of aquaporins in gas permeation has been a long-standing and controversially discussed issue. We present an extensive set of atomistic molecular dynamics simulations that address the question of CO(2) permeation through human aquaporin-1. Free energy profiles derived from the simulations display a barrier of approximately 23 kJ/mol in the aromatic/arginine constriction region of the water pore, whereas a barrier of approximately 4 kJ/mol was observed for a palmitoyloleoylphosphatidylethanolamine lipid bilayer membrane. The results indicate that significant aquaporin-1-mediated CO(2) permeation is to be expected only in membranes with a low intrinsic CO(2) permeability.     

Synergistic effect of electric field and lipid oxidation on the permeability of cell membranes

BackgroundStrong electric fields are known to affect cell membrane permeability, which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on how this affects the cell membrane permeability.MethodWe here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.Results and conclusionsWe show how oxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.General significanceThis study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.