Cascades of transient pores in giant vesicles: line tension and transport

Under ordinary circumstances, the membrane tension of a giant unilamellar vesicle is essentially nil. Using visible light, we stretch the vesicles, increasing the membrane tension until the membrane responds by the sudden opening of a large pore (several micrometers in size). Only a single pore is observed at a time in a given vesicle. However, a cascade of transient pores appear, up to 30-40 in succession, in the same vesicle. These pores are transient: they reseal within a few seconds as the inner liquid leaks out. The membrane tension, which is the driving force for pore opening, is relaxed with the opening of a pore and the leakage of the inner liquid; the line tension of the pore's edge is then able to drive the closure of a pore. We use fluorescent membrane probes and real-time videomicroscopy to study the dynamics of the pores. These can be visualized only if the vesicles are prepared in a viscous solution to slow down the leakout of the internal liquid. From measurements of the closure velocity of the pores, we are able to infer the line tension,. We have studied the effect of the shape of inclusion molecules on. Cholesterol, which can be modeled as an inverted cone-shaped molecule, increases the line tension when incorporated into the bilayers. Conversely, addition of cone-shaped detergents reduces. The effect of some detergents can be dramatic, reducing by two orders of magnitude, and increasing pore lifetimes up to several minutes. We give some examples of transport through transient pores and present a rough measurement of the leakout velocity of the inner liquid through a pore. We discuss how our results can be extended to less viscous aqueous solutions which are more relevant for biological systems and biotechnological applications.     

Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

Dynamics of fusion pores connecting membranes of different tensions

The energetics underlying the expansion of fusion pores connecting biological or lipid bilayer membranes is elucidated. The energetics necessary to deform membranes as the pore enlarges, in some combination with the action of the fusion proteins, must determine pore growth. The dynamics of pore growth is considered for the case of two homogeneous fusing membranes under different tensions. It is rigorously shown that pore growth can be quantitatively described by treating the pore as a quasiparticle that moves in a medium with a viscosity determined by that of the membranes. Motion is subject to tension, bending, and viscous forces. Pore dynamics and lipid flow through the pore were calculated using Lagrange's equations, with dissipation caused by intra- and intermonolayer friction. These calculations show that the energy barrier that restrains pore enlargement depends only on the sum of the tensions; a difference in tension between the fusing membranes is irrelevant. In contrast, lipid flux through the fusion pore depends on the tension difference but is independent of the sum. Thus pore growth is not affected by tension-driven lipid flux from one membrane to the other. The calculations of the present study explain how increases in tension through osmotic swelling of vesicles cause enlargement of pores between the vesicles and planar bilayer membranes. In a similar fashion, swelling of secretory granules after fusion in biological systems could promote pore enlargement during exocytosis. The calculations also show that pore expansion can be caused by pore lengthening; lengthening may be facilitated by fusion proteins.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Transbilayer Pores Induced by Thickness Fluctuations

Thermally-induced fluctuations of individual phospholipids in a bilayer lipid membrane (BLM) are converted into collective motions due to the intermolecular interactions. Here, we demonstrate that transbilayer stochastic pores can be generated via collective thermal movements (CTM). Using the elastic theory of continuous media applied to smectic-A liquid crystals, we estimate the pore radius and the energetic requirements for pore appearance. Three types of thermally-induced transbilayer pores could be formed through BLMs: open and stable, open and unstable, and closed. In most of the situations, two open and stable pores with different radii could be generated. Notably, the two pores have the same generation probability. Unstable pores are possible to appear across thin bilayers that contain phospholipids with a large polar headgroup. Closed pores are present throughout the cases that we have inspected. The effects of hydrophobic thickness, polar headgroup size of phospholipids, temperature, surface tension, and elastic compression on the pore formation and pore stability have been examined as well.     

The exocytotic fusion pore modeled as a lipidic pore.

Freeze-fracture electron micrographs from degranulating cells show that the lumen of the secretory granule is connected to the extracellular compartment via large (20 to 150 nm diameter) aqueous pores. These exocytotic fusion pores appear to be made up of a highly curved bilayer that spans the plasma and granule membranes. Conductance measurements, using the patch-clamp technique, have been used to study the fusion pore from the instant it conducts ions. These measurements reveal the presence of early fusion pores that are much smaller than those observed in electron micrographs. Early fusion pores open abruptly, fluctuate, and then either expand irreversibly or close. The molecular structure of these early fusion pores is unknown. In the simplest extremes, these early fusion pores could be either ion channel like protein pores or lipidic pores. Here, we explored the latter possibility, namely that of the early exocytotic fusion pore modeled as a lipid-lined pore whose free energy was composed of curvature elastic energy and work done by tension. Like early exocytotic fusion pores, we found that these lipidic pores could open abruptly, fluctuate, and expand irreversibly. Closure of these lipidic pores could be caused by slight changes in lipid composition. Conductance distributions for stable lipidic pores matched those of exocytotic fusion pores. These findings demonstrate that lipidic pores can exhibit the properties of exocytotic fusion pores, thus providing an alternate framework with which to understand and interpret exocytotic fusion pore data.     

The exocytotic fusion pore interface: a model of the site of neurotransmitter release

infrastructurel techniques have shown that an early event in the exocytotic fusion of a secretory vesicle is the formation of a narrow, water-filled pore spanning both the vesicle and plasma membranes and connecting the lumen of the secretory vesicle to the extracellular environment. Smaller precursors of the exocytotic fusion pore have been detected using electrophysio-logical techniques, which reveal a dynamic fusion pore that quickly expands to the size of the pores seen with electron microscopy. While it is clear that in the latter stages of expansion, when the size of the fusion pore is several orders of magnitude bigger than any known macromolecule, the fusion pore must be mainly made of lipids, the structure of the smaller precursors is unknown. Patch-clamp measurements of the activity of individual fusion pores in mast cells have shown that the fusion pore has some unusual and unexpected properties, namely that there is a large flux of lipid through the pore and the rate of pore closure has a discontinuous temperature dependency, suggesting a purely lipidic fusion pore. Moreover, comparisons of experimental data with theoretical fusion pores and with breakdown pores support the view that the fusion pore is initially a pore through a single bilayer, as would be expected for membrane fusion proceeding through a hemifusion mechanism. Based on these observations we present a model where the fusion pore is initially a pore through a single bilayer. Fusion pore formation is regulated by a macromolecular scaffold of proteins that is responsible for bringing the plasma membrane into a highly curved dimple very close to a tense secretory granule membrane, creating the architecture where the strongly attractive hydrophobic force causes the membranes to form a ‘hemifusion’ intermediate. Membrane fusion is completed by the formation of an aqueous pore after rupture of the shared bilayer. We also propose that the microenvironment of the interface when the pore first opens, dominated by the charged groups on the secretory vesicle matrix and phospholipids, will greatly influence the release of secretory products.     

PB1-F2, an influenza A virus-encoded proapoptotic mitochondrial protein, creates variably sized pores in planar lipid membranes

A frameshifted region of the influenza A virus PB1 gene encodes a novel protein, termed PB1-F2, a mitochondrial protein that can induce cell death. Many proapoptotic proteins are believed to act at the mitochondrial outer membrane to form an apoptotic pore with lipids. We studied the interaction of isolated, synthetic PB1-F2 (sPB1-F2) peptide with planar phospholipid bilayer membranes. The presence of nanomolar concentrations of peptide in the bathing solution induced a transmembrane conductance that increased in a potential-dependent manner. Positive potential on the side of protein addition resulted in a severalfold increase in the rate of change of membrane conductance. sPB1-F2-treated membranes became permeable to monovalent cations, chloride, and to a lesser extent, divalent ions. Despite various experimental conditions, we did not detect the distinctive conductance levels typical of large, stable pores, protein channels, or even pores that are partially proteinaceous. Rather, membrane conductance induced by sPB1-F2 fluctuated and visited almost all conductance values. sPB1-F2 also dramatically decreased bilayer stability in an electric field, consistent with a decrease in the line tension of a lipidic pore. Since similar membrane-destabilizing profiles are seen with proapoptotic proteins (e.g., Bax) and the cytoplasmic helix of human immunodeficiency virus gp41, we suggest that the basis for sPB1-F2-induced cell death may be the permeabilization and destabilization of mitochondrial membranes, leading to macromolecular leakage and apoptosis.     

Mechanism of pore formation in stearoyl-oleoyl-phosphatidylcholine membranes subjected to lateral tension

Theoretical model of a through pore formation in lipid bilayer membrane under applied lateral tension was developed. In the framework of elastic theory of liquid crystals adapted to lipid membranes, we calculated a continuous trajectory from intact bilayer through a hydrophobic defect to a through pore. It was shown that the major energetic characteristic of membrane stability with respect to the pore formation, i. e., line tension, depends both on the pore radius and on the value of the applied lateral tension. This leads to a non-monotonous dependence of the average waiting time of the pore formation on the lateral tension: at low tensions the waiting time was large, then there was a local minimum, after which the average waiting time was increasing again. For membranes formed from stearoyl oleoyl phosphatidylcholine, the local minimum corresponded to the lateral tension of 7 mN/m; the calculated value of the edge line tension of a large pore was 16.5 pN. These results are consistent with available experimental data.     

Membrane mechanics can account for fusion pore dilation in stages.

Once formed, fusion pores rapidly enlarge to semi-stable conductance values. The membranes lining the fusion pore are continuous bilayer structures, so variations of conductance in time reflect bending and stretching of membranes. We therefore modeled the evolution of fusion pores using the theory of the mechanics of deforming homogeneous membranes. We calculated the changes in length and width of theoretical fusion pores according to standard dynamical equations of motion. Theoretical fusion pores quickly achieve semi-stable dimensions, which correspond to energy minima located in a canyon between energy barriers. The height of the barrier preventing pore expansion diminishes along the dimensions of length and width. The bottom of the canyon slopes gently downward along increasing length. As a consequence, theoretical fusion pores slowly lengthen and widen as the dimensions migrate along the bottom of the canyon, until the barrier vanishes and the pore rapidly enlarges. The dynamics of growth is sensitive to tension, spontaneous curvature, bending elasticity, and mobilities. This sensitivity can account for the quantitative differences in pore evolution observed in two experimental systems: HA-expressing cells fusing to planar bilayer membranes and beige mouse mast cell degranulation. We conclude that the mechanics of membranes could cause the phenomenon of stagewise growth of fusion pores.     

Interaction of hagfish cathelicidin antimicrobial peptides with model lipid membranes

Hagfish intestinal antimicrobial peptides (HFIAPs) are a family of polycationic peptides exhibiting potent, broad-spectrum bactericidal activity. In an attempt to unravel the mechanism of action of HFIAPs, we have studied their interaction with model membranes. Synthetic HFIAPs selectively bound to liposomes mimicking bacterial membranes, and caused the release of vesicle-encapsulated fluorescent markers in a size-dependent manner. In planar lipid bilayer membranes, HFIAPs induced erratic current fluctuations and reduced membrane line tension according to a general theory for lipidic pores, suggesting that HFIAP pores contain lipid molecules. Consistent with this notion, lipid transbilayer redistribution accompanied HFIAP pore formation, and membrane monolayer curvature regulated HFIAP pore formation. Based on these studies, we propose that HFIAPs kill target cells, at least in part, by interacting with their plasma membrane to induce formation of lipid-containing pores. Such a membrane-permeabilizing function appears to be an evolutionarily conserved host-defense mechanism of antimicrobial peptides.     

Protein Arcs May Form Stable Pores in Lipid Membranes

Electron microscopy and atomic force microscopy images of cholesterol-dependent cytolysins and related proteins that form large pores in lipid membranes have revealed the presence of incomplete rings, or arcs. Some evidence indicates that these arcs are inserted into the membrane and induce membrane leakage, but other experiments seem to refute that. Could such pores, only partially lined by protein, be kinetically and thermodynamically stable? How would the lipids be structured in such a pore? Using the antimicrobial peptide protegrin-1 as a model, we test the stability of pores only partially lined by peptide using all-atom molecular dynamics simulations in POPC and POPE/POPG membranes. The data show that, whereas pure lipid pores close rapidly, pores partially lined by protegrin arcs are stable for at least 300 ns. Estimates of the thermodynamic stability of these arcs using line tension data and implicit solvent calculations show that these arcs can be marginally stable in both zwitterionic and anionic membranes. Arcs provide an explanation for the observed ion selectivity in protegrin electrophysiology experiments and could possibly be involved in other membrane permeabilization processes where lipids are thought to participate, such as those induced by antimicrobial peptides and colicins, as well as the Bax apoptotic pore.     

Protein Arcs May Form Stable Pores in Lipid Membranes

Electron microscopy and atomic force microscopy images of cholesterol-dependent cytolysins and related proteins that form large pores in lipid membranes have revealed the presence of incomplete rings, or arcs. Some evidence indicates that these arcs are inserted into the membrane and induce membrane leakage, but other experiments seem to refute that. Could such pores, only partially lined by protein, be kinetically and thermodynamically stable? How would the lipids be structured in such a pore? Using the antimicrobial peptide protegrin-1 as a model, we test the stability of pores only partially lined by peptide using all-atom molecular dynamics simulations in POPC and POPE/POPG membranes. The data show that, whereas pure lipid pores close rapidly, pores partially lined by protegrin arcs are stable for at least 300 ns. Estimates of the thermodynamic stability of these arcs using line tension data and implicit solvent calculations show that these arcs can be marginally stable in both zwitterionic and anionic membranes. Arcs provide an explanation for the observed ion selectivity in protegrin electrophysiology experiments and could possibly be involved in other membrane permeabilization processes where lipids are thought to participate, such as those induced by antimicrobial peptides and colicins, as well as the Bax apoptotic pore.     

Pore Formation in a Binary Giant Vesicle Induced by Cone-Shaped Lipids

We have investigated shape deformations of binary giant unilamellar vesicles (GUVs) composed of cone- and cylinder-shaped lipids. By coupling the spontaneous curvature of lipids with the phase separation, we demonstrated pore opening and closing in GUVs. When the temperature was set below the chain melting transition temperature of the cylinder-shaped lipid, the GUVs burst and then formed a single large pore, where the cone shape lipids form a cap at the edge of the bilayer to stabilize the pore. The pore closed when we increased the temperature above the transition temperature. The pore showed three types of shapes depending on the cone-shaped lipid concentration: simple circular, rolled-rim, and wrinkled-rim pores. These pore shape changes indicate that the distribution of the cone- and cylinder-shaped lipids is asymmetric between the inner and outer leaflets in the bilayer. We have proposed a theoretical model for a two-component membrane with an edge of bilayer where lipids can transfer between two leaflets. Using this model, we have reproduced numerically the observed shape deformations at the rim of pore.     

Solubilization kinetics determines the pulsatory dynamics of lipid vesicles exposed to surfactant

We establish a biophysical model for the dynamics of lipid vesicles exposed to surfactants. The solubilization of the lipid membrane due to the insertion of surfactant molecules induces a reduction of membrane surface area at almost constant vesicle volume. This results in a rate-dependent increase of membrane tension and leads to the opening of a micron-sized pore. We show that solubilization kinetics due to surfactants can determine the regime of pore dynamics: either the pores open and reseal within a second (short-lived pore), or the pore stays open up to a few minutes (long-lived pore). First, we validate our model with previously published experimental measurements of pore dynamics. Then, we investigate how the solubilization kinetics and membrane properties affect the dynamics of the pore and construct a phase diagram for short and long-lived pores. Finally, we examine the dynamics of sequential pore openings and show that cyclic short-lived pores occur with a period inversely proportional to the solubilization rate. By deriving a theoretical expression for the cycle period, we provide an analytical tool to estimate the solubilization rate of lipid vesicles by surfactants. Our findings shed light on some fundamental biophysical mechanisms that allow simple cell-like structures to sustain their integrity against environmental stresses, and have the potential to aid the design of vesicle-based drug delivery systems. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.     

A lipocentric view of peptide-induced pores

Although lipid membranes serve as effective sealing barriers for the passage of most polar solutes, nonmediated leakage is not completely improbable. A high activation energy normally keeps unassisted bilayer permeation at a very low frequency, but lipids are able to self-organize as pores even in peptide-free and protein-free membranes. The probability of leakage phenomena increases under conditions such as phase coexistence, external stress or perturbation associated to binding of nonlipidic molecules. Here, we argue that pore formation can be viewed as an intrinsic property of lipid bilayers, with strong similarities in the structure and mechanism between pores formed with participation of peptides, lipidic pores induced by different types of stress, and spontaneous transient bilayer defects driven by thermal fluctuations. Within such a lipocentric framework, amphipathic peptides are best described as pore-inducing rather than pore-forming elements. Active peptides bound to membranes can be understood as a source of internal surface tension which facilitates pore formation by diminishing the high activation energy barrier. This first or immediate action of the peptide has some resemblance to catalysis. However, the presence of membrane-active peptides has the additional effect of displacing the equilibrium towards the pore-open state, which is then maintained over long times, and reducing the size of initial individual pores. Thus, pore-inducing peptides, regardless of their sequence and oligomeric organization, can be assigned a double role of increasing the probability of pore formation in membranes to high levels as well as stabilizing these pores after they appear.     

Pore formation by a Bax-derived peptide: effect on the line tension of the membrane probed by AFM

Bax is a critical regulator of physiological cell death that increases the permeability of the outer mitochondrial membrane and facilitates the release of the so-called apoptotic factors during apoptosis. The molecular mechanism of action is unknown, but it probably involves the formation of partially lipidic pores induced by Bax. To investigate the interaction of Bax with lipid membranes and the physical changes underlying the formation of Bax pores, we used an active peptide derived from helix 5 of this protein (Bax-alpha5) that is able to induce Bax-like pores in lipid bilayers. We report the decrease of line tension due to peptide binding both at the domain interface in phase-separated lipid bilayers and at the pore edge in atomic force microscopy film-rupture experiments. Such a decrease in line tension may be a general strategy of pore-forming peptides and proteins, as it affects the energetics of the pore and stabilizes the open state.     

Paradoxical lipid dependence of pores formed by the Escherichia coli alpha-hemolysin in planar phospholipid bilayer membranes

alpha-Hemolysin (HlyA) is an extracellular protein toxin (117 kDa) secreted by Escherichia coli that targets the plasma membranes of eukaryotic cells. We studied the interaction of this toxin with membranes using planar phospholipid bilayers. For all lipid mixtures tested, addition of nanomolar concentrations of toxin resulted in an increase of membrane conductance and a decrease in membrane stability. HlyA decreased membrane lifetime up to three orders of magnitude in a voltage-dependent manner. Using a theory for lipidic pore formation, we analyzed these data to quantify how HlyA diminished the line tension of the membrane (i.e., the energy required to form the edge of a new pore). However, in contrast to the expectation that adding the positive curvature agent lysophosphatidylcholine would synergistically lower line tension, its addition significantly stabilized HlyA-treated membranes. HlyA also appeared to thicken bilayers to which it was added. We discuss these results in terms of models for proteolipidic pores.     

Pores formed in lipid bilayer membranes by nystatin, Differences in its one-sided and two-sided action

Nystatin and amphotericin B induce a cation-selective conductance when added to one side of a lipid bilayer membrane and an anion-selective conductance when added to both sides. The concentrations of antibiotic required for the one-sided action are comparable to those employed on plasma membranes and are considerably larger than those required for the two-sided action. We propose that the two-sided effect results from the formation of aqueous pores formed by the hydrogen bonding in the middle of the bilayer of two "half pores," whereas the one-sided effect results from the half pores alone. We discuss, in terms of the flexibility of bilayer structure and its thickness, how it is possible to have conducting half pores and "complete pores" in the same membrane. The role of sterol (cholesterol and ergosterol) in pore formation is also examined.