Interface water dynamics and porating electric fields for phospholipid bilayers

Lipid bilayers, normally a barrier to charged species and large molecules, are permeabilized by electric fields, a phenomenon exploited by cell biologists and geneticists for porating and transfecting cells and tissues. Recent molecular simulation studies have advanced our understanding of electroporation, but the relative contributions of atomically local details (interface water and headgroup dipole and counterion configurations) and medium- and long-range electrostatic gradients and changes in membrane structural shifts to the initiating conditions and mechanisms of pore formation remain unclear. Molecular dynamics simulations of electroporation in several lipid systems presented here reveal the effects of lipid hydrocarbon tail length and composition on the magnitude of the field required for poration and on the location of the initial sites of field-driven water intrusion into the bilayer. Minimum porating external fields of 260 mV nm(-1), 280 mV nm(-1), 320 mV nm(-1), and 380 mV nm(-1) were found for 1,2-dilauroyl-sn-glycero-3-phosphatidylcholine (DLPC), 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), and 1,2-dioleoyl- sn-glycero-3-phosphatidylcholine (DOPC) bilayers, respectively, and correlated most strongly with the bilayer thickness. These phospholipid systems share several common features including a wide, dynamic distribution of the headgroup dipole angle with the bilayer normal ranging from 0 to 155 degrees that is only slightly shifted in a porating electric field, and similar electric field-induced shifts in water dipole orientation, although the mean water dipole moment profile at the aqueous-membrane interface is more sensitive to the electric field for DOPC than for the other phospholipids. The location of pore initiation, at the anode- or cathode-facing leaflet, varies with the composition of the bilayer and correlates with a change in the polarity of the localized electric field at the interface.     

Phase behavior of polarizable spherocylinders in external fields

Applied electric fields are known to induce significant changes in the properties of systems of polarizable molecules or particles. For rod-shaped molecules, the field-induced behavior can be rather surprising, as in the case of the negative electric birefringence of concentrated solutions of rodlike polyelectrolytes. We have investigated the interplay of shape anisotropy and field-induced anisotropy in molecular dynamics simulations of systems of polarizable soft spherocylinders in an electric field, in the limit of infinitely anisotropic polarizability, taking full account of mutual induction effects. We find a novel crystalline structure (K(2)) in the high-field limit, whose formation is driven by interactions between induced dipoles. For high pressures, the phase diagram exhibits a polar nematic phase between the hexagonal close-packed crystal phase and the K(2) phase. We also compare this system with an analogous system of spherocylinders with permanent electric dipoles and find that qualitatively similar behavior is obtained in the limit of strong coupling of the permanent dipoles to the external field.     

Electrophoretic transport and dynamic deformation of bio‐vesicles

Study of the deformation dynamics of cells and other sub‐micron vesicles, such as virus and neurotransmitter vesicles are necessary to understand their functional properties. This mechanical characterization can be done by submerging the vesicle in a fluid medium and deforming it with a controlled electric field, which is known as electrodeformation. Electrodeformation of biological and artificial lipid vesicles is directly influenced by the vesicle and surrounding media properties and geometric factors. The problem is compounded when the vesicle is naturally charged, which creates electrophoretic forcing on the vesicle membrane. We studied the electrodeformation and transport of charged vesicles immersed in a fluid media under the influence of a DC electric field. The electric field and fluid‐solid interactions are modeled using a hybrid immersed interface‐immersed boundary technique. Model results are verified with experimental observations for electric field driven translocation of a virus through a nanopore sensor. Our modeling results show interesting changes in deformation behavior with changing electrical properties of the vesicle and the surrounding media. Vesicle movement due to electrophoresis can also be characterized by the change in local conductivity, which can serve as a potential sensing mechanism for electrodeformation experiments in solid‐state nanopore setups.     

外电场作用下水化聚全氟磺酸钾膜中水分子的电渗迁移

用分子动力学模拟方法研究外电场(简称电场)作用下水化聚全氟磺酸钾膜中水分子的电渗迁移运动,并分析探讨膜的结构与水分子的电渗迁移特性的关系.结果表明,无外加电场时水分子和K+的速度都服从麦克斯韦分布;施加适当电场时,水分子和K+在垂直电场方向上的速度分量仍服从麦克斯韦分布,但平行电场方向上的速度分量则服从峰值漂移的麦克斯韦分布.并且,峰值漂移速度可作为水分子和K+的平均迁移速度的近似值,从而计算得到水分子的电渗系数.结果还显示,K+第一配位层内平均含有约4.04个水分子,它们的平均迁移速度只有K+的57%.这部分水分子贡献的电渗迁移系数为总电渗迁移系数(2.97)的77%.     

Membrane surface dynamics of DNA-threaded nanopores revealed by simultaneous single-molecule optical and ensemble electrical recording

We describe a method for simultaneous single-molecule optical and electrical characterization of membrane-based sensors that contain ion-channel nanopores. The technique is used to study the specific and nonspecific interactions of streptavidin-capped DNA polymers with lipid bilayers composed of diphytanoyl phosphatidylcholine and diphytanoyl phosphatidylglycerol. Biotinylated DNA that is bound to fluorescently labeled streptavidin is electrophoretically driven into, or away from, the lumen of alpha hemolysin (alphaHL) ion channels by an external electric field. Confocal microscopy simultaneously captures single-molecule fluorescence dynamics from the membrane interface at different applied potentials. Fluorescence correlation analysis is used to determine the surface number density and diffusion constant of membrane-associated complexes. The dual optical and electrical approach can detect membrane-associated species at a surface coverage below 10(-5) monolayers of streptavidin, a sensitivity that surpasses most other in vitro surface analysis techniques. By comparing the change in transmembrane current to the number of fluorescent molecules leaving the bilayer when the electrical potential is reversed, we demonstrate the general utility of the approach within the context of nanopore-based sensing and discuss a mechanism by which DNA-streptavidin complexes can be nonspecifically retained at the membrane interface.     

Molecular-dynamics simulation of the effect of ions on a liquid-liquid interface for a partially miscible mixture

Molecular-dynamics simulations were performed to model the effect of added salt ions on the liquid-liquid interface in a partially miscible system. Simulations of the interface between saturated phases of a model 1-hexanol+water system show a bilayer structure of 1-hexanol molecules at the interface with -OH heads of the first layer directed into the water phase and the opposite orientation for the second layer. The alignment of the polar -OH groups at the interface stabilizes a charge separation of sodium and chloride ions when salt is introduced into the aqueous phase, producing an electrical double layer. Chloride ions aggregate nearer the interface and sodium ions move toward the bulk water phase, consistent with the explanation that the -OH alignment presents a region of partial positive charges to which the hydrated chloride atoms are attracted. Ions near the interface were found to be less solvated than those in the bulk phase. An electric field was also applied to drive ions through the interface. Ions crossing the interface tended to shed water molecules as they entered the hexanol bilayer, leaving a trail of water molecules. Stabilization and facilitated transport of the ion by interactions with the second layer of hexanol molecules appeared to be an important step in the mechanism of sodium ion transport.     

Implicit-solvent mesoscale model based on soft-core potentials for self-assembled lipid membranes

An efficient implicit-solvent model for self-assembled lipid bilayers is presented and analyzed using Langevin molecular dynamics simulations. The model is based on soft interactions between particles and short-range attractive interaction between lipid tails, leading for the self-assembly of a lipid bilayer without an explicit solvent. This allows for efficient simulations of large membranes over long times. The model exhibits a fluid phase at high temperatures and a gel phase at low temperatures, identified with the Lbeta-phase. The melting transition is investigated via analysis of the diffusivity of the lipid molecules, the chain-orientational order parameter, the sixfold bond-orientational order parameter, and the positional and bond-orientational correlation functions. The analysis suggests the existence of a hexatic phase over a narrow range of temperatures around the melting transition. The elastic properties of the membrane in the fluid phase are also investigated.     

Comparative simulation study of nitrogen and ammonia adsorption on graphitized and nongraphitized carbon blacks

Grand canonical Monte Carlo simulation is used to study the adsorption of nitrogen at 77 K and ammonia at 240 K to represent weakly polar and polar molecules, respectively, on infinite and finite graphite surfaces. These graphite surfaces were modeled with different percentages of carbons removed (defects) from the top graphite layer. Increasing the number of defects increases the adsorption and the isosteric heat of nitrogen at low pressure. At moderate pressures the amount adsorbed is less due to the disruption in the packing of the nitrogen in the first layer. In contrast, the adsorption of ammonia at all pressures is reduced as the percentage of defects is increased. This is due to the disruption in ammonia bonding caused by the defects. The condensation-like step change in the ammonia isotherm on the perfect graphite surface is not observed for any of these surfaces with defects even for the case of only 10% defects. At high percentage of defects the adsorption isotherm is close to Henry law behavior for much of the pressure range. The adsorption on finite surfaces shows that the amount adsorbed for both molecules decreases compared with that of the infinite surfaces, resulting from interaction potentials with the surface and other fluid molecules at the edge. The decrease is much greater for the ammonia adsorption because the bonding between ammonia molecules is disrupted, meaning that the adsorption cannot follow the mechanism of condensation seen for the infinite surface.     

Lipids and lipid mixtures in boundary layers: From hydration lubrication to osteoarthritis

The hydration layer surrounding the phosphocholine headgroups of single-component phosphatidylcholine lipids, or of lipid-mixtures, assembled at an interface greatly modifies the interfacial properties and interactions. As water molecules within the hydration layer are held tightly by the headgroup but are nonetheless very fluid on shear, the boundary lipid layers, exposing the highly hydrated headgroup arrays, can provide efficient boundary lubrication when sliding against an opposing surface, at physiologically high contact pressures. In addition, any free lipids in the surrounding liquid can heal defects which may form during sliding on the boundary phosphatidylcholine layer. Similar boundary lipid layers contribute to the lubricating, pressure-bearing, and wear-protection functions of healthy articular joints. This review presents a survey of the relationship between the molecular composition of the interfacial complex and the lubrication behavior of the lipid-based boundary layers, which could be beneficial for designing boundary lubricants for intra-articular injection for the treatment of early osteoarthritis.     

Effect of electric field during gel-layer controlled ultrafiltration of synthetic and fruit juice

Electric field enhanced ultrafiltration of pectin–sucrose mixture (synthetic juice) and mosambi (Citrus sinensis (L.) Osbeck) fruit juice using 50,000 (MWCO) polyerthersulfon membrane is studied in a cross-flow cell. Pectin, completely rejected by the membrane, forms a gel type layer over the membrane surface. Under the application of an external dc electric field across the membrane, gel-layer formation is restricted leading to an enhancement of permeate flux. During ultrafiltration of synthetic juice, application of dc electric field (800 V/m) increases the permeate flux to almost threefold compared to that with zero electric field. A theoretical model based on integral method assuming suitable concentration profile in the boundary layer is developed. The proposed model is used to predict the permeate flux in gel-layer governed electric field enhanced ultrafiltration. Predictions of the model are successfully compared with the experimental results under a wide range of operating conditions. Experiments with fruit juice also demonstrated significant increase in flux with the application of a suitable electric field.     

Poly(amidoamine) dendrimers on lipid bilayers II: Effects of bilayer phase and dendrimer termination

The molecular structures and enthalpy release of poly(amidoamine) (PAMAM) dendrimers binding to 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) bilayers were explored through atomistic molecular dynamics. Three PAMAM dendrimer terminations were examined: protonated primary amine, neutral acetamide, and deprotonated carboxylic acid. Fluid and gel lipid phases were examined to extract the effects of lipid tail mobility on the binding of generation-3 dendrimers, which are directly relevant to the nanoparticle interactions involving lipid rafts, endocytosis, lipid removal, and/or membrane pores. Upon binding to gel phase lipids, dendrimers remained spherical, had a constant radius of gyration, and approximately one-quarter of the terminal groups were in close proximity to the lipids. In contrast, upon binding to fluid phase bilayers, dendrimers flattened out with a large increase in their asphericity and radii of gyration. Although over twice as many dendrimer-lipid contacts were formed on fluid versus gel phase lipids, the dendrimer-lipid interaction energy was only 20% stronger. The greatest enthalpy release upon binding was between the charged dendrimers and the lipid bilayer. However, the stronger binding to fluid versus gel phase lipids was driven by the hydrophobic interactions between the inner dendrimer and lipid tails.     

Streaming potential in cylindrical pores of poly(ethylene terephthalate) track-etched membranes: variation of apparent zeta potential with pore radius

Streaming potential variation with pressure measured through poly(ethylene terephthalate) track-etched membranes of different pore sizes led to the determination of an apparent interfacial potential zetaa in the presence of 10-2 M KCl. The variation of zetaa with the pore radius r0 is interpreted by (i) the electric double layer overlap effect and (ii) the presence of a conductive gel layer. We propose a method which integrates both effects by assuming a simple model for the conductive gel at the pore wall. We observed the following three domains of pore size: (i) r0 > 70 nm, where surface effects are negligible; (ii) approximately 17 nm < r0 < 70 nm, where the pore/solution interface could be described as a conductive gel of thickness around 1 nm; (iii) r0 < approximately 17 nm, which corresponds to the region strongly damaged by the ion beam and is not analyzed here. The first one (zeta = -36.2 mV) corresponds to the raw material when etching has completely removed the ion beam predamaged region, which corresponds to the second intermediate domain (zeta = -47.3 mV). There the conductance of the gel layer deduced from the treatment of streaming potential data was found to be compatible with the number of ionic sites independently determined by the electron spin resonance technique.     

Nanopore formation and phosphatidylserine externalization in a phospholipid bilayer at high transmembrane potential

Atomic-resolution molecular dynamics simulations of lipid bilayers containing 7% phosphatidylserine (PS) on one leaflet are consistent with experimental observations of membrane poration and PS externalization in living cells exposed to nanosecond, megavolt-per-meter electric pulses. Nanometer-diameter aqueous pores develop within nanoseconds after application of an electric field of 450 mV/nm, and electrophoretic transport of the anionic PS headgroup along the newly constructed hydrophilic pore surface commences even while pore formation is still in progress.     

Forming Lipid Bilayer Membrane Arrays on Micropatterned Polyelectrolyte Film Surfaces

A novel method of forming lipid bilayer membrane arrays on micropatterned polyelectrolyte film surfaces is introduced. Polyelectrolyte films were fabricated by the layer‐by‐layer technique on a silicon oxide surface modified with a 3‐aminopropyltriethoxysilane (APTES) monolayer. The surface pK_a value of the APTES monolayer was determined by cyclic voltammetry to be approximately 5.61, on the basis of which a pH value of 2.0 was chosen for layer‐by‐layer assembly. Micropatterned polyelectrolyte films were obtained by deep‐UV (254 nm) photolysis though a mask. Absorbed fluorescent latex beads were used to visualize the patterned surfaces. Lipid bilayer arrays were fabricated on the micropatterned surfaces by immersing the patterned substrates into a solution containing egg phosphatidylcholine vesicles. Fluorescence recovery after photobleaching studies yielded a lateral diffusion coefficient for probe molecules of 1.31±0.17 μm² s^?1 in the bilayer region, and migration of the lipid NBD PE in bilayer lipid membrane arrays was observed in an electric field.     

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.     

Interactions of poly(amidoamine) dendrimers with Survanta lung surfactant: the importance of lipid domains

The interaction of generation 5 (G5) and 7 (G7) poly(amidoamine) (PAMAM) dendrimers with mica-supported Survanta bilayers is studied with atomic force microscopy (AFM). In these experiments, Survanta forms distinct gel and fluid domains with differing lipid composition. Nanoscale defects are induced by the PAMAM dendrimers. The positively charged dendrimers remove lipid from the fluid domains at a significantly greater rate than for the gel domains. Dendrimer accumulation on lipid edges and terraces preceding lipid removal has been directly imaged. Immediately following lipid removal, the mica surface is clean, indicating that lipid defects are not induced by dendrimers binding to the mica substrate and displacing the lipid.     

Elastic deformation of membrane bilayers probed by deuterium NMR relaxation

In deuterium ((2)H) NMR spectroscopy of fluid lipid bilayers, the average structure is manifested in the segmental order parameters (S(CD)) of the flexible molecules. The corresponding spin-lattice relaxation rates (R(1Z) depend on both the amplitudes and the rates of the segmental fluctuations, and indicate the types of lipid motions. By combining (2)H NMR order parameter measurements with relaxation studies, we have obtained a more comprehensive picture of lipids in the liquid-crystalline (L(alpha)) state than formerly possible. Our data suggest that a lipid bilayer constitutes an ordered fluid, in which the phospholipids are grafted to the aqueous interface via their polar headgroups, whereas the fatty acyl chains are in effect liquid hydrocarbon. Studies of (2)H-labeled saturated lipids indicate their R(1Z) rates and S(CD) order parameters are correlated by a model-free, square-law functional dependence, signifying the presence of relatively slow bilayer fluctuations. A new composite membrane deformation model explains simultaneously the frequency (magnetic field) dependence and the angular anisotropy of the relaxation. The results imply the R(1Z) rates are due to a broad spectrum of 3-D collective bilayer excitations, together with effective axial rotations of the lipids. For the first time, NMR relaxation studies show that the viscoelastic properties of membrane lipids at megahertz frequencies are modulated by the lipid acyl length (bilayer thickness), polar headgroups (bilayer interfacial area), inclusion of a nonionic detergent (C(12)E(8)), and the presence of cholesterol, leading to a range of bilayer softness. Our findings imply the concept of elastic deformation is relevant on lengths approaching the bilayer thickness and less (the mesoscopic scale), and suggest that application of combined R(1Z) and S(CD) studies of phospholipids can be used as a simple membrane elastometer. Heuristic estimates of the bilayer bending rigidity kappa and the area elastic modulus K(a) enable comparison to other biophysical studies, involving macroscopic deformation of thin membrane lipid films. Finally, the bilayer softness may be correlated with the lipid diversity of biomembranes, for example, with regard to membrane curvature, repulsive interactions between bilayers, and lipid-protein interactions.     

Lateral diffusion of molecules in two-component lipid bilayer: a Monte Carlo simulation study

Lateral diffusion of membrane components makes possible any in-plane membrane reaction and has a key role in signaling in cell membranes. In this report the equilibrium lateral diffusion of intrinsic molecules in an equimolar DMPC/DSPC mixture is simulated using a thoroughly tested two-state model of two-component phospholipid bilayers. The model has been successful in calculating the excess heat capacity function, the most frequent center-to-center distances between DSPC clusters, and the fractal dimensions of gel clusters (Sugar, I. P., Thompson, T. E., Biltonen, R. L. Biophys. J. 1999, 76, 2099-2110). In the gel/fluid mixed phase region, a diffusing intrinsic molecule may change its state from fluid to gel (or from gel to fluid) at any time. A common characterization of the diffusion of intrinsic molecules is given by the simulated average first-passage time curves. We find that these curves can be described as power functions containing two parameters, alpha and beta, except near the percolation threshold of gel/fluid or compositional clusters. We find also that the intrinsic molecules are involved in approximately normal diffusion, i.e., beta approximately 2 in the extreme gel and fluid phase regions, while in the gel/fluid and gel/gel mixed phase regions the diffusion is anomalous, i.e., beta not equal 2. In the mixed phase regions, when the initial local state of the diffusing molecule is not specified, each component is involved in sub-diffusion (beta > 2). In the gel/fluid mixed phase region molecules situated initially inside a fluid cluster are involved in sub-diffusion, but DMPC molecules situated initially inside a gel cluster are involved in super-diffusion (beta < 2). The possibility of anomalous diffusion in membranes apparently arises because the diffusing molecule visits a variety of different environments characterized by its relative proximity to various membrane components. The diffusion is actually anomalous when the components of the bilayer are nonrandomly distributed. The deviation from random distribution is strongly correlated with beta. Similar to the results of the NMR experiments, the calculated relative diffusion coefficient continuously decreases in the gel/fluid mixed phase region with decreasing temperature. In apparent contradiction, diffusion measured by fluorescence recovery after photobleaching (FRAP) demonstrates the existence of a threshold temperature, below which long-range diffusion of FRAP probe molecules is essentially blocked. This threshold temperature is highly correlated with the percolation temperature of gel clusters.     

Molecular dynamics simulation study of the interaction of trehalose with lipid membranes

The interactions of the cryoprotective agent trehalose with a lipid membrane made of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine at 323 K were studied by means of molecular dynamics simulations. It was observed that trehalose binds to the phospholipid headgroups with its main axis parallel to the membrane normal. Trehalose establishes hydrogen bonds with the carbonyl and phosphate groups and replaces water molecules from the lipid headgroup. Notably, the number of hydrogen bonds (HBs) that the membrane made with its environment was conserved after trehalose binding. The HBs between lipid and trehalose have a longer lifetime than those established between lipid and water. The binding of the sugar does not produce changes either in the lipid area or in the lipid order parameter. The effect of trehalose on the dipole potential is in agreement with experimental results. The contribution of the different components to the membrane dipole potential was analyzed. It was observed that the binding of trehalose produces changes in the different components and the sugar itself contributes to the surface potential due to the polarization of its hydroxyl in the interface.