耗散粒子动力学研究双层膜线张力与抗弯刚度的关系

利用耗散力粒子动力学模拟方法研究了拉伸状态下双层膜的抗弯刚度与膜孔线张力之间的关系.通过对双层膜在不同投影面积约束下的系统模拟,观察到3个区域:自由振动膜、伸展膜和穿孔膜.由前两个区域的应力张量计算得到的膜的面张力(σΣ)与拟合膜波动性质得到的面张力(σfluc)吻合的很好,除去在两区域的转变点附近σfluc略大于σΣ.当考虑在经典Helfrich弹性膜模型中被忽略的膜厚度时,线张力可以和抗弯刚度用一个简单的模型联系起来.通过对穿孔膜区域的数据进行分析,证明由本模型得到的抗弯刚度与拟合膜波动性质得到的抗弯刚度符合的很好.由此提出一种简便的测量方法,通过计算拉伸膜孔的面张力和统计膜厚度,拟合这个简单模型来测量膜的抗弯刚度.     

Size-dependent properties of small unilamellar vesicles formed by model lipids

The size-dependent behavior of small unilamellar vesicles is explored by dissipative particle dynamics, including the membrane characteristics and mechanical properties. The spontaneously formed vesicles are in the metastable state and the vesicle size is controlled by the concentration of model lipids. As the vesicle size decreases, the bilayer gets thinner and the area density of heads declines. Nonetheless, the area density in the inner leaflet is higher than that in the outer. The packing parameters are calculated for both leaflets. The result indicates that the shape of lipid in the outer leaflet is like a truncated cone but that in the inner leaflet resembles an inverted truncated cone. Based on a local order parameter, our simulations indication that the orientation order of lipid molecules decreases as the size of the vesicle reduces and this fact reveals that the bilayer becoming thinner for smaller vesicle is mainly attributed to the orientation disorder of the lipids. The membrane tension can be obtained through the Young-Laplace equation. The tension is found to grow with reducing vesicle size. Therefore, small vesicles are less stable against fusion. Using the inflation method, the area stretching and bending moduli can be determined and those moduli are found to grow with reducing size. Nonetheless, a general equation with a single numerical constant can relate bending modulus, area stretching modulus, and bilayer thickness irrespective of the vesicle size. Finally, a simple metastable model is proposed to explain the size-dependent behavior of bilayer thickness, orientation, and tension.     

Lipid bilayer elasticity measurements in giant liposomes in contact with a solubilizing surfactant

A new method to probe the modification of the elasticity of phospholipid bilayers is presented. The purpose here concerns the action of a solubilizing surfactant on a vesicle bilayer. This method is based on the measure of the under-field elongation of giant magnetic-fluid-loaded liposomes. The addition of the nonionic surfactant octyl-beta-d-glucopyranoside (OG) to vesicles at sublytic levels increases the elasticity of the membrane, as shown by the value of the bending modulus K(b), which decreases. K(b) measured around 20 kT for a pure 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayer indeed reaches a few kT in the case of the mixed OG-DOPC bilayer. The purpose and interest of this study are to allow the determination of the membrane bending modulus before and after the addition of OG on the same magnetic liposome. Moreover, the experimental conditions used in this work allow the control of lipid and surfactant molar fractions in the mixed aggregates. Then, optical microscopy observation can be performed on samples in well-defined regions of the OG-phospholipid state diagram.     

Dynamics and instabilities of lipid bilayer membrane shapes

Biological membranes undergo constant shape remodeling involving the formation of highly curved structures. The lipid bilayer represents the fundamental architecture of the cellular membrane with its shapes determined by the Helfrich curvature bending energy. However, the dynamics of bilayer shape transitions, especially their modulation by membrane proteins, and the resulting shape instabilities, are still not well understood. Here, we review in a unifying manner several theories that describe the fluctuations (i.e. undulations) of bilayer shapes as well as their local coupling with lipid or protein density variation. The coupling between local membrane curvature and lipid density gives rise to a ‘slipping mode’ in addition to the conventional ‘bending mode’ for damping the membrane fluctuation. This leads to a number of interesting experimental phenomena regarding bilayer shape dynamics. More importantly, curvature-inducing proteins can couple with membrane shape and eventually render the membrane unstable. A criterion for membrane shape instability is derived from a linear stability analysis. The instability criterion reemphasizes the importance of membrane tension in regulating the stability and dynamics of membrane geometry. Recent progresses in understanding the role of membrane tension in regulating dynamical cellular processes are also reviewed. Protein density is emphasized as a key factor in regulating membrane shape transitions: a threshold density of curvature coupling proteins is required for inducing membrane morphology transitions.     

Groove-spanning behavior of lipid membranes on microfabricated silicon substrates

We report on a spreading behavior of phospholipid membranes that arise from a lump of phospholipid (a lipid source) on topographically patterned substrates immersed in an aqueous solution. Microgrooves with well-defined shapes were prepared on Si111 surfaces by anisotropic etching in an alkaline solution. A spreading front that consists of membrane lobes and a single lipid bilayer was observed on the patterned silicon substrates by utilizing fluorescence interference contrast (FLIC) microscopy. FLIC images indicate that the membrane lobes span the microgrooves, while the underlying single lipid bilayer spread along the surface of the microgrooves. In fact, fluorescent polystyrene nanoparticles could be encapsulated in the microgrooves that were completely covered with the membrane lobes. The groove-spanning behavior of membrane lobes is discussed in terms of a balance between adhesion and bending energies of lipid bilayers.     

Mechanical properties of surfactant bilayer membranes from atomistic and coarse-grained molecular dynamics simulations

We use simulations to predict the stability and mechanical properties of two amphiphilic bilayer membranes. We carry out atomistic MD simulations and investigate whether it is possible to use an existing coarse-grained (CG) surfactant model to map the membrane properties. We find that certain membranes can be represented well by the CG model, whereas others cannot. Atomistic MD simulations of the erucate membrane yield a headgroup area per surfactant a(0) of 0.26 nm(2), an elastic modulus K(A) of 1.7 N/m, and a bending rigidity kappa of 5 k(B)T. We find that the CG model, with the right choice for the size and potential well depth of the head, correctly reproduces a(0), kappa, as well as the fluctuation spectrum over the whole range of q values. Atomistic MD simulations of EHAC, on the other hand, suggest that this membrane is unstable. This is indicated by the fact that kappa is of the order of k(B)T, which means that the interface is extremely flexible and diffuse, and K(A) is close to zero, which means that the surface tension is zero. We argue that the CG model can be used if the headgroups are uncharged, dipolar, or effectively dipolar due to headgroup charge screening induced by counterion condensation.     

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.     

A novel method for measuring the bending rigidity of model lipid membranes by simulating tethers

The tensile force along a cylindrical lipid bilayer tube is proportional to the membrane's bending modulus and inversely proportional to the tube radius. We show that this relation, which is experimentally exploited to measure bending rigidities, can be applied with even greater ease in computer simulations. Using a coarse-grained bilayer model we efficiently obtain bending rigidities that compare very well with complementary measurements based on an analysis of thermal undulation modes. We furthermore illustrate that no deviations from simple quadratic continuum theory occur up to a radius of curvature comparable to the bilayer thickness.     

The cooperative role of membrane skeleton and bilayer in the mechanical behaviour of red blood cells

Red blood cell (RBC) shape, behaviour and deformability can be consistently accounted for by a model for the elastic properties of the RBC membrane that includes the elasticity of the membrane skeleton in dilation and shear, and the local and nonlocal resistance of the bilayer to bending. The role of the corresponding energy terms in different RBC shape and deformation situations is analyzed. RBC shape transformations are compared to the shape transformations of phospholipid vesicles that are driven by the difference between the equilibrium areas of the bilayer leaflets (DeltaA0). It is deduced that the skeleton energy contributions play a crucial role in the formation of an echinocyte. The effect of a transformation of the natural biconcave RBC shape into an echinocyte on its resistance to entry into capillary-sized cylindrical tubes is analyzed. It is shown that, during the aspiration of an echinocyte into a pipette, there are two competing skeleton deformation effects, which arise due to skeleton density changes, one due to spicule formation and the other due to deformation induced by micropipette aspiration. Furthermore, the shift of the observed dependence of the projection length on the aspiration pressure of more crenated cells towards higher aspiration pressures can be accounted for by an increase of the equilibrium area difference DeltaA0 and consequent modification of the nonlocal contribution to the cell elastic energy.     

Composition dependence of bilayer elasticity

A previously developed molecular level model for homogeneous lipid bilayers [Brannigan and Brown, J. Chem. Phys 120, 1059 (2004)] is extended to allow for multiple lipid species. Monte Carlo simulations (including species exchange moves for efficient sampling) reveal a variety of mixing behaviors in binary systems. Two species are identified that maintain stable, randomly mixed fluid membranes at vanishing tension over all possible binary compositions. The thermal and elastic properties of membranes formed by these lipids are characterized over the full composition range. Equilibrium area at constant tension is nonmonotonic with respect to composition, but consistent with that of a quadratic mixture. In the constant tension ensemble, the bending rigidity of the bilayer is minimized at an intermediate composition. The observed functional form of bending rigidity vs composition is fit to a simple expression motivated by linear elasticity theory; this expression accounts for membrane heterogeneity through a single parameter.     

Membrane orientation of MSI-78 measured by sum frequency generation vibrational spectroscopy

Antimicrobial peptides (AMPs) selectively disrupt bacterial cell membranes to kill bacteria whereas they either do not or weakly interact with mammalian cells. The orientations of AMPs in lipid bilayers mimicking bacterial and mammalian cell membranes are related to their antimicrobial activity and selectivity. To understand the role of AMP-lipid interactions in the functional properties of AMPs better, we determined the membrane orientation of an AMP (MSI-78 or pexiganan) in various model membranes using sum frequency generation (SFG) vibrational spectroscopy. A solid-supported single 1,2-dipalmitoyl-an-glycero-3-[phospho-rac-(1-glycerol)] (DPPG) bilayer or 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (POPG) bilayer was used as a model bacterial cell membrane. A supported 1,2-dipalmitoyl-an-glycero-3-phosphocholine (DPPC) bilayer or a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer was used as a model mammalian cell membrane. Our SFG results indicate that the helical MSI-78 molecules are associated with the bilayer surface with ～70° deviation from the bilayer normal in the negatively charged gel-phase DPPG bilayer at 400 nM peptide concentration. However, when the concentration was increased to 600 nM, MSI-78 molecules changed their orientation to make a 25° tilt from the lipid bilayer normal whereas multiple orientations were observed for an even higher peptide concentration in agreement with toroidal-type pore formation as reported in a previous solid-state NMR study. In contrary, no interaction between MSI-78 and a zwitterionic DPPC bilayer was observed even at a much higher peptide concentration (～12,000 nM). These results demonstrate that SFG can provide insights into the antibacterial activity and selectivity of MSI-78. Interestingly, the peptide exhibits a concentration-dependent membrane orientation in the lamellar-phase POPG bilayer and was also found to induce toroidal-type pore formation. The deduced lipid flip-flop from SFG signals observed from lipids also supports MSI-78-induced toroidal-type pore formation.     

Insertion and assembly of membrane proteins via simulation

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.     

Electrostatic properties of membrane lipids coupled to metarhodopsin II formation in visual transduction

Changes in lipid composition have recently been shown to exert appreciable influences on the activities of membrane-bound proteins and peptides. We tested the hypothesis that the conformational states of rhodopsin linked to visual signal transduction are related to biophysical properties of the membrane lipid bilayer. For bovine rhodopsin, the meta I-meta II conformational transition was studied in egg phosphatidylcholine (PC) recombinants versus the native rod outer segment (ROS) membranes by means of flash photolysis. Formation of metarhodopsin II was observed by the change in absorbance at 478 nm after a single actinic flash was delivered to the sample. The meta I/meta II ratio was investigated as a function of both temperature and pH. The data clearly demonstrated thermodynamic reversibility of the transition for both the egg PC recombinants and the native ROS membranes. A significant shift of the apparent pK(a) for the acid-base equilibrium to lower values was evident in the egg PC recombinant, with little meta II produced under physiological conditions. Calculations of the membrane surface pH using a Poisson-Boltzmann model suggested the free energies of the meta I and meta II states were significantly affected by electrostatic properties of the bilayer lipids. In the ROS membranes, phosphatidylserine (PS) is needed for full formation of meta II, in combination with phosphatidylethanolamine (PE) and polyunsaturated docosahexaenoic acid (DHA; 22:6omega3) chains. We propose that the PS surface potential leads to an accumulation of hydronium ions, H(3)O(+), in the electrical double layer, which drive the reaction together with the large negative spontaneous curvature (H(0)) conferred by PE plus DHA chains. The elastic stress/strain of the bilayer arises from an interplay of the approximately zero H(0) from PS and the negative H(0) due to the PE headgroups and polyunsaturated chains. The lipid influences are further explained in terms of matching of the bilayer spontaneous curvature to the curvature at the lipid/rhodopsin interface, as formulated by the Helfrich bending energy. These new findings guide current ideas as to how bilayer properties govern the conformational energetics of integral membrane proteins. Moreover, they yield knowledge of how membrane lipid-protein interactions involving acidic phospholipids such as PS and neutral polyunsaturated DHA chains are implicated in key biological functions such as vision.     

Effect of chain length and asymmetry on material properties of bilayer membranes

Dissipative particle dynamics is used to extract the material parameters (bending and area stretch moduli) of a bilayer membrane patch. Some experiments indicate that the area stretch modulus of lipid vesicles varies little as the chain length of the lipids composing the bilayer increases. Here we show that making the interactions between the hydrophilic head groups of the model amphiphiles proportional to the hydrophobic tail length reproduces the above result for the area stretch modulus. We also show that the area stretch modulus of bilayers composed of amphiphiles with the same number of tail beads but with asymmetric chains is less than that of bilayers with symmetric chains. The effects on the bilayer density and lateral stress profiles of changes to the amphiphile architecture are also presented.     

Simulations of edge behavior in a mixed-lipid bilayer: fluctuation analysis

Coarse-grained molecular dynamics simulations of lipid bilayer ribbons consisting of a mixture of lipids of different tail lengths have been performed to gain insight into bicelle mixtures. The line tension of the bilayer edge decreases as the mole fraction of short-chain lipids in the system is increased, dropping below zero between 30% and 35%. The mole fraction of short-chain lipids in the ribbon interior is lower than the total mole fraction, as the short-chain lipids segregate towards the edge, but continues to rise even after the line tension vanishes, in contrast to predictions of a two-component two-phase model. The fluctuations of the bilayer edge in both high and low line tension regimes have been analyzed to extract information about the factors that influence the length and shape of the edge. At high line tension the wavelength-dependent in-plane fluctuations of the edge are predicted quantitatively using a simple analytical model using only the line tension as input. Where line tension is vanishing, the fluctuations can be modeled as arising from a combination of harmonic fluctuations around a minimum energy contour length and an in-plane bending elasticity. The estimated value of the in-plane bending modulus is of order 10(-29) J m, placing the intrinsic persistence length for the edge near the bilayer thickness of 4 nm.     

单烷基乙二胺稀水溶液头基相互作用自组装形成 双分子膜

提出普通表面活性剂(单链两亲分子)亲水头基相互作用诱导疏水尾链平行聚集形成双分子膜的新机制。设计和合成了系列单烷基取代乙二胺C~nH~2~n~+~1NHC~2H~4NH~2(n=8,12,14,16,18)。通过电镜形态,分散液凝胶/液晶相变和对应铸膜的二维双层结构,表明单链两亲分子头基相互作用和脂链引入刚性片断一样,两者形成的双分子膜具有类似的结构和性能;展示了各体系取代乙二胺双层结构和性能的密切联系。指出了广泛认同的单链两亲分子形成双分子膜必须引入刚性片断的单一成膜机制的片面性,为组装新一类功能头基表面活性剂双分子膜独辟蹊径。     

Agar-supported lipid bilayers — basic structures for biosensor design. Electrical and mechanical properties

The lipid bilayer is widely accepted as the basic structure of all biological membranes. Known as BLM (bilayer lipid membrane), it can be prepared artificially. Suitably modified, the BLM serves as a very appropriate model for biological membranes. Recent investigations have verified the high analytical potential of artificial lipid membranes. With a structure and composition almost identical to the lipid moiety of biomembranes, the BLM may serve as an ideal host for receptor molecules of biological origin, thus becoming a transducer which could “see” the environment the way the living cell does. For the construction of lipid bilayer based biosensors; however, stable, easy to prepare and long-lasting lipid membranes are required. With this aim in mind, we have prepared lipid bilayer membranes which use an agar gel as support. This as-BLM (agar-supported BLM) has been shown to possess the same electrical, mechanical and dynamic properties the conventional BLM is famous for, along with the benefits of long-term stability and considerably elevated breakdown voltages. Its preparation on the tip of an agar-filled Teflon tube of 0.5 mm diameter is easy and can be performed even by less-skilled personnel.

In an attempt of further miniaturization the concept of the as-BLM was applied to thin-film micro-systems manufactured by standard micro-electronic techniques. The result is a lipid bilayer system, which, while preserving all the essential properties of the bilayer lipid membrane, can serve as a basic building block for cheap, disposable biosensoric systems.     

Voltammetric study of gold-supported lipid membranes in the presence of perfluorooctanesulphonic acid

The effect of perfluorooctanesulphonic acid (PFOS) on lipid membranes was studied using supported 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer as the model membrane. Phospholipid bilayer was deposited on gold electrode using a combination of the Langmuir–Blodgett and Langmuir–Schaefer (LB/LS) techniques. Electrodes were modified with two different types of membranes: DMPC bilayers initially containing PFOS and pure DMPC bilayers later exposed to the PFOS solutions. Such approach allowed studying both the changes in membrane characteristic imposed by the perfluorinated compound present in the model membrane and the process of its incorporation into the membrane. Studies with anticancer drug doxorubicin revealed that PFOS inhibits drug transport through the phospholipid bilayer and its effect can be compared to that of cholesterol. Moreover, the different trends observed in the changes in electron transfer rate constant (k_s) calculated for ferricyanides and in peak current of hexaamineruthenium chloride showed that electrostatic interactions between electroactive probes and PFOS molecules incorporating into phospholipid bilayers play an important role and should be taken into account while explaining the interactions of perfluorooctanesulphonic acid with model biological membranes.     

Effect of bilayer thickness on membrane bending rigidity

The bending rigidity k(c) of bilayer vesicles self-assembled from amphiphilic diblock copolymers has been measured using single- and dual-micropipet techniques. These copolymers are nearly a factor of 5 greater in hydrophobic membrane thickness d than their lipid counterparts and an order of magnitude larger in molecular weight M(n). The macromolecular structure of these amphiphiles lends insight into and extends relationships for traditional surfactant behavior. We find the scaling of k(c) with thickness to be nearly quadratic, in good agreement with existing theories for bilayer membranes. The results here are key to understanding and designing soft interfaces such as biomembrane mimetics.     

Surface response methodology for the study of supported membrane formation

We report on the investigations of the formation of the tethered lipid bilayer by vesicle deposition on amine-functionalized surfaces. The tethered bilayer was created by the deposition of egg-PC vesicles containing 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-poly-(ethyleneglycol)-N-hydroxysuccinimide as anchoring molecules on an amine-coated surface. This approach is an easy route for the formation of a biomimetic-supported membrane. A Doelhert experimental design was applied to determine the conditions leading to the formation of a continuous and defect-free tethered bilayer on different surfaces (gold and glass). Doehlert designs allow modeling of the experimental responses by second-order polynomial equations as a function of experimental factors. Four factors expected to influence bilayer formation were studied: the lipid concentration in the vesicle suspension, the mass percentage of anchoring molecules in the vesicles, the contact time between the vesicles and the surface, and the resting time of the membrane after buffer rinse. The optimization of the membrane preparation parameters was achieved by monitoring lipid assembly formation using surface plasmon resonance spectroscopy on gold and by fluorescence recovery after photobleaching on glass. Three characteristic responses were systematically measured: the bilayer thickness, the lipid diffusion coefficient, and the lipid mobile fraction. The simultaneous inspection of the three characteristics revealed that a restricted experimental domain leads to properties that are in accordance with a bilayer presence. The factors of this domain are a lipid concentration from 0.1 to 1 mg/mL, 4-8% of anchoring molecules in the vesicles, 1-4 h of contact time between vesicles and surface, and 21-24 h of resting time after buffer rinse. Under these conditions, a membrane having a lipid mass per surface between 545 +/- 5 and 590 +/- 10 ng/cm2, a diffusion coefficient of between 2.5 +/- 0.3 x 10(-8) and 3.60 +/- 0.5 x 10(-8) cm2/s, and a mobile fraction between 94 +/- 2 and 99 +/- 1% was formed. These findings were confirmed by atomic force microscopy observations, which showed the presence of a continuous and homogeneous bilayer in the determined experimental domain. This formation procedure presents many advantages; it provides an easily obtainable biomimetic membrane model for proteins studies and offers a versatile tethered bilayer because it can be adapted easily to various types of supports.