Pore formation and rupture in fluid bilayers

Encapsulation and release of calcein from herceptin-conjugated eLiposomes

Achieving an optimal therapeutic level is crucial in effectively eradicating cancer cells during treatment. However, conventional chemotherapy-associated systemic administration of anticancer agents leads to many side effects. To achieve the desired control over the target site, active targeting of HER2-positive breast cancer cells can be achieved by conjugating liposomal vesicles with Human Epidermal growth factor Receptor 2 (HER2) and inducing release of the encapsulated drug using ultrasound. To further enhance the delivery efficiency, nanoemulsion droplets exhibiting responsiveness to low-frequency ultrasound are encapsulated within these lipid vesicles. In this study, we prepared four different liposomal formulations, namely pegylated liposomes, emulsion liposomes (eLiposomes), HER-conjugated liposomes, and HER-conjugated eLiposomes, each loaded with calcein and subjected to a thorough characterization process. Their sizes, phospholipid concentration, and amount of antibody conjugation were compared and analyzed. Cryogenic transmission electron microscopy was used to confirm the encapsulation of nanoemulsion droplets within the liposomes. The drug-releasing performance of Herceptin-conjugated eLiposomes was found to surpass that of other liposomal formulations with a notably higher calcein release and established it as a highly effective nanocarrier. The study showcases the efficacy of calcein-loaded and Herceptin-conjugated eLiposomes, which demonstrate rapid and efficient drug release among other liposomal formulations when subjected to ultrasound. This discovery paves the way for a more targeted, efficient, and humane approach to cancer therapy.

Supramolecular Engineering of Alkylated, Fluorinated, and Mixed Amphiphiles

The rational design of perfluorinated amphiphiles to control the supramolecular aggregation in aqueous medium is still a key challenge for the engineering of supramolecular architectures. Here we present the synthesis and physical properties of six novel non-ionic amphiphiles. We also studied the effect of mixed alkylated and perfluorinated segments in a single amphiphile and compared it with only alkylated and perfluorinated units. To explore their morphological behavior in aqueous medium, we used dynamic light scattering (DLS) and cryo-TEM/EM measurements. We further confirmed their assembly mechanisms with theoretical investigations, using the Martini model to perform large-scale coarse-grained molecular dynamics simulations. These novel synthesized amphiphiles offer a greater and more systematic understanding of how perfluorinated systems assemble in aqueous medium and suggest new directions for rational designing of new amphiphilic systems and interpreting their assembly process. This article is protected by copyright. All rights reserved.

Selecting a Cell Engineering Methodology During Cell Therapy Product Development

When considering the development pathway for a genetically modified cell therapy product, it is critically important that the product is engineered consistent with its intended human use. For scientists looking to develop and commercialize a new technology, the decision to select a genetic modification method depends on several practical considerations. Whichever path is chosen, the developer must understand the key risks and potential mitigations of the cell engineering approach. The developer should also understand the clinical implications: permanent/memory establishment versus transient expression, and clinical manufacturing considerations when dealing with transplantation of genetically engineered cells. This review covers important topics for mapping out a strategy for developers of new cell-based therapeutics. Biological, technological, manufacturing, and clinical considerations are all presented to map out development lanes for the initiation and risk management of new gene-based cell therapeutic products for human use.

Harnessed Dopant Block Copolymers Assist Decorating Membrane Pores: A Dissipative Particle Dynamics Study

Self-assembly of asymmetric block copolymers (BCPs) around active pore edges has emerged as an important strategy to produce smart membranes with tunable pathways for solute transport. However, thus far, it is still challenging to manipulate pore shape and functionality for directional transformation under external stimuli. Here, a versatile strategy by mesoscale simulations to design stimuli-responsive pores with various edge decorations in hybrid membranes is reported. Dopant BCPs are used as decorators to stabilize pore edges and extend their function in reconfiguring pores in response to repeated membrane stretching/shrinking caused by external stimuli. The decoration morphologies are predictable since the assemblies of dopant BCPs around pore edges are closely related to their self-assemblies in solution. The coassembly between different BCPs in the hybrid membrane for the control of pore morphology is featured, and the parameter settings, including block incompatibility and molecular architecture for the construction of a specific pore, are determined. Results show that harnessed dopant BCPs in the hybrid membrane can enhance pore formation and induce directional pore shape and functionality transformation. Diversified pore decorations exhibit potential that can be further explored in selective solute transport and the design of stimuli-responsive smart nanodevices.

Metastable Prepores in Tension-Free Lipid Bilayers

The formation and closure of aqueous pores in lipid bilayers is a key step in various biophysical processes. Large pores are well described by classical nucleation theory, but the free-energy landscape of small, biologically relevant pores has remained largely unexplored. The existence of small and metastable "prepores" was hypothesized decades ago from electroporation experiments, but resolving metastable prepores from theoretical models remained challenging. Using two complementary methods-atomistic simulations and self-consistent field theory of a minimal lipid model-we determine the parameters for which metastable prepores occur in lipid membranes. Both methods consistently suggest that pore metastability depends on the relative volume ratio between the lipid head group and lipid tails: lipids with a larger head-group volume fraction (or shorter saturated tails) form metastable prepores, whereas lipids with a smaller head-group volume fraction (or longer unsaturated tails) form unstable prepores.

Direct Measurement of the Critical Pore Size in a Model Membrane

We study pore nucleation in a model membrane system, a freestanding polymer film. Nucleated pores smaller than a critical size close, while pores larger than the critical size grow. Holes of varying size were purposefully prepared in liquid polymer films, and their evolution in time was monitored using optical and atomic force microscopy to extract a critical radius. The critical radius scales linearly with film thickness for a homopolymer film. The results agree with a simple model which takes into account the energy cost due to surface area at the edge of the pore. The energy cost at the edge of the pore is experimentally varied by using a lamellar-forming diblock copolymer membrane. The underlying molecular architecture causes increased frustration at the pore edge resulting in an enhanced cost of pore formation.

Line tension of multicomponent bilayer membranes

The line tension or edge energy of bilayer membranes self-assembled from binary amphiphilic molecules is studied using self-consistent-field theory (SCFT). Specifically, solutions of the SCFT equations corresponding to an infinite membrane with a circular pore, or an open membrane, are obtained for a coarse-grained model in which the amphiphilic species and hydrophilic solvents are represented by ABandED diblock copolymers and C homopolymers, respectively. The edge energy of the membrane is extracted from the free energy of the open membranes. Results for membranes composed of mixtures of symmetric and cone- or inverse cone-shaped amphiphilic molecules with neutral and/or repulsive interactions are obtained and analyzed. It is observed that an increase in the concentration of the cone-shaped species leads to a decrease of the line tension. In contrast, adding inverse cone-shaped copolymers results in an increase of the line tension. Furthermore, the density profile of the copolymers reveals that the line tension is regulated by the distribution of the amphiphiles at the bilayer edge.

Ultrasound-induced calcein release from eLiposomes

Ultrasound is explored as a method of inducing the release of encapsulated materials from eLiposomes, defined as liposomes containing emulsion droplets. Emulsions were formed using perfluorohexane and perfluoropentane. eLiposomes were formed by folding interdigitated lipid sheets into closed vesicles around the emulsion droplets. Cryogenic transmission electron microscopy was used to verify droplet encapsulation. Self-quenched calcein was also encapsulated inside the vesicles. A fluorometer was used to measure baseline fluorescence, calcein release after ultrasound exposure, and total release from the vesicles. eLiposome samples released 3 to 5 times more of the encapsulated calcein than did controls when exposed to 20-kHz ultrasound. Calcein release increased with exposure time and intensity of ultrasound. eLiposomes with large (400 nm) droplets produced more calcein release than small (100 nm) droplets. These observations suggest that the emulsions are vaporized by ultrasound and that the Laplace pressure in the emulsions has an effect on droplet vaporization.

Encapsulating nanoemulsions inside eLiposomes for ultrasonic drug delivery

An eLiposome is a liposome encapsulating an emulsion nanodroplet and can be used for drug delivery. For example, therapeutic agents are encapsulated inside the eLiposomes, and the application of ultrasound can cause the emulsion droplet to change from liquid to gas, thus increasing the volume inside the vesicle and causing rupture and the release of the drug. In this research, two different methods were used to prepare eLiposomes. In the first method, emulsion droplets were made of perfluorohexane or perfluoropentane and stabilized with 1,2-dipalmitoyl-sn-glycero-3-phosphate. A layer of 1,2-dimyristoyl-sn-glycero-3-phosphocholine was dried in a round-bottomed flask. Then the emulsion suspension was added to the flask. As the suspension hydrated the phospholipids, they formed liposomes around the emulsions. In the second method, emulsions and liposomes were made separately, and then they were mixed using ultrasound. The advantage of this second method compared to the previous one is that eLiposomes can be made with fewer restrictions because of incompatible combinations of surfactants. Dynamic light scattering and transmission electron microscopy were used to measure the size of the emulsions, liposomes, and eLiposomes. The size of eLiposomes is appropriate for extravasation into tumors with malformed capillary beds. We hypothesize that ultrasound breaks open these eLiposomes. Both types of eLiposomes were constructed with folate attached via a poly(ethylene glycol) tether to induce endocytosis of the eLiposome. The latter eLiposomes were successfully used to deliver calcein as a model drug to HeLa cells.

Exploration of the shapes of double-walled vesicles with a confined inner membrane

We investigate double-walled vesicles as a simple model system for multi-vesicular structures, where the inner membrane is confined within the outer membrane. Various shapes of double-walled vesicles in two dimensions are obtained by means of our recently-developed discrete space variation method, and the shapes of each layer are found to be interdependent. Confined within the outer membrane, an inner membrane with a larger surface area always shows a cristae shape. As previous simulations and theoretical analyses of a single-walled vesicle have been done before, the geometric properties of double-walled vesicles, including the mean square radius of gyration and volume within the vesicle membrane, are studied in detail as functions of the pressure and surface area. It is found that due to the inter-space restriction of each layer, double-walled vesicles exhibit different behaviors compared with the previously-observed scaling laws of single-walled vesicles. It is straightforward to extend this study to more complicated and realistic biological systems, such as those including electrostatic interactions between membranes and solvent, phase separation, and cooperative interactions between multicomponent membranes.

Minimum energy path to membrane pore formation and rupture

We combine dynamic self-consistent field theory with the string method to calculate the minimum energy path to membrane pore formation and rupture. In the regime where nucleation can occur on experimentally relevant time scales, the structure of the critical nucleus is between a solvophilic stalk and a locally thinned membrane. Classical nucleation theory fails to capture these molecular details and significantly overestimates the free energy barrier. Our results suggest that thermally nucleated rupture may be an important factor for the low rupture strains observed in lipid membranes.

Interaction of lipid membrane with nanostructured surfaces

Tiny details of the phospholipid (DMPC) membrane morphology in close vicinity to nanostructured silica surfaces have been discovered in the atomic force microscopy experiments. The structural features of the silica surface were varied in the experiments by the deposition of silica nanoparticles of different diameter on plane and smooth silica substrates. It was found that, due to the barrier function of the lipid membrane, only particles larger than 22 nm in diameter with a smooth surface were completely enveloped by the lipid membrane. However, nanoparticles with bumpy surfaces (curvature diameter of bumps as that of particles <22 nm) were only partially enveloped by the lipid bilayer. For the range of nanostructure dimensions between 1.2 and 22 nm, the lipid membrane underwent structural rearrangements by forming pores (holes). The nanoparticles were accommodated into the pores but not enveloped by the lipid bilayer. The study also found that the lipid membrane conformed to the substrate with surface structures of dimensions less than 1.2 nm without losing the membrane integrity. The experimental results are in accord with the analytical free energy model, which describes the membrane coverage, and numerical simulations which evaluate adhesion of the membrane and dynamics as a function of surface topology. The results obtained in this study are useful for the selection of dimensions and shapes for drug-delivery cargo and for the substrate for supported lipid bilayers. They also help in qualitative understanding the role of length scales involved in the mechanisms of endocytosis and cytotoxicity of nanoparticles. These findings provide a new approach for patterning supported lipid membranes with well-defined features in the 1.2-22 nm range.

Stability and rupture of archaebacterial cell membrane: a model study

It is known that the thermoacidophilic archaebacterium Sulfolobus acidocaldarius can grow in hot springs at 65-80 degrees C and live in acidic environments (pH 2-3); however, the origin of its unusual thermal stability remains unclear. In this work, using a vesicle as a model, we study the thermal stability and rupture of archaebacterial cell membrane. We perform a simulation investigation of the structure-property relationship of monolayer membrane formed by bolaform lipids and compare it with that of bilayer membrane formed by monopolar lipids. The origin of the unusually thermal stability of archaebacterial cell and the mechanism for its rupture are presented in molecular details.

Interaction of nanoparticles with lipid membrane

A nanoscale range of surface feature curvatures where lipid membranes lose integrity and form pores has been found experimentally. The pores were experimentally observed in the l-alpha-dimyristoyl phosphatidylcholine membrane around 1.2-22 nm polar nanoparticles deposited on mica surface. Lipid bilayer envelops or closely follows surface features with the curvatures outside of that region. This finding provides essential information for the understanding of nanoparticle-lipid membrane interaction, cytotoxicity, preparation of biomolecular templates and supported lipid membranes on rough and patterned surfaces.

Shear-induced intracellular loading of cells with molecules by controlled microfluidics

This study tested the hypothesis that controlled flow through microchannels can cause shear-induced intracellular loading of cells with molecules. The overall goal was to design a simple device to expose cells to fluid shear stress and thereby increase plasma membrane permeability. DU145 prostate cancer cells were exposed to fluid shear stress in the presence of fluorescent cell-impermeant molecules by using a cone-and-plate shearing device or high-velocity flow through microchannels. Using a syringe pump, cell suspensions were flowed through microchannels of 50-300 microm diameter drilled through Mylar sheets using an excimer laser. As quantified by flow cytometry, intracellular uptake and loss of viability correlated with the average shear stress. Optimal results were observed when exposing the cells to high shear stress for short durations in conical channels, which yielded uptake to over one-third of cells while maintaining viability at approximately 80%. This method was capable of loading cells with molecules including calcein (0.62 kDa), large molecule weight dextrans (150-2,000 kDa), and bovine serum albumin (66 kDa). These results supported the hypothesis that shear-induced intracellular uptake could be generated by flow of cell suspensions through microchannels and further led to the design of a simple, inexpensive, and effective device to deliver molecules into cells. Such a device could benefit biological research and the biotechnology industry.

Pore formation in a lipid bilayer under a tension ramp: modeling the distribution of rupture tensions

The rupture of fluid membrane vesicles with a steady ramp of micropipette suction has been shown to produce a distribution of breakage tensions, with a mean that rises rapidly with tension rate. Starting from a lattice model that incorporates the essential features of the lipid bilayers held together with hydrophobic forces, and developing it to handle varying tension rates, we reproduce the main features of the experimental results. In essence, we show that the rupture kinetics are driven by the nucleation and growth of pores, with two limiting kinetics-growth-limited and nucleation-limited. The model has been extended to address the role of peptides in solution that can adsorb and insert themselves into the bilayer. At concentrations below those required to spontaneously rupture the membrane, the effect of the peptides is to lower the rupture tensions systematically for all tension rates.

Nuclear pores and membrane holes: generic models for confined chains and entropic barriers in pore stabilization

The lumen of the nuclear pore complex is increasingly understood to be lined by a polymer brush that entropically regulates transport in and out of the nucleus-and it seems likely that similar effects probably arise with glycocalyx-lined holes in cell membranes. Here we mimic such pore-confined brushes with self-assembled polymer membranes imbued with nano-holes. Experiment and theory help elucidate the entropic origin and stabilization of the pores, which appear to have a similar basis as steric stabilization of colloids bearing polymer brushes. Free energies of interacting brushes reveal stable minima at pore sizes smaller than the classical metastable point, with little effect of the particular pore geometry. Such entropic forces have potential implications for lock and key mechanisms of nuclear pore assembly as well as transient poration of cells and synthetic nano-pores with regulatory mechanisms for transport.

Pores in bilayer membranes of amphiphilic molecules: coarse-grained molecular dynamics simulations compared with simple mesoscopic models

We investigate pores in fluid membranes by molecular dynamics simulations of an amphiphile-solvent mixture, using a molecular coarse-grained model. The amphiphilic membranes self-assemble into a lamellar stack of amphiphilic bilayers separated by solvent layers. We focus on the particular case of tensionless membranes, in which pores spontaneously appear because of thermal fluctuations. Their spatial distribution is similar to that of a random set of repulsive hard disks. The size and shape distribution of individual pores can be described satisfactorily by a simple mesoscopic model, which accounts only for a pore independent core energy and a line tension penalty at the pore edges. In particular, the pores are not circular: their shapes are fractal and have the same characteristics as those of two-dimensional ring polymers. Finally, we study the size-fluctuation dynamics of the pores, and compare the time evolution of their contour length to a random walk in a linear potential.

Simulations of stable pores in membranes: system size dependence and line tension

Amphiphilic bilayers with a pore were simulated using a coarse grained model. By stretching the bilayer to 70% beyond its equilibrium surface area, we established the phase diagram of pores, identifying regions where pores are stable, metastable, or unstable. A simple theoretical model is proposed to explain the phase diagram, and to calculate the critical and equilibrium relative stretches. Interestingly, these are found to scale with the inverse cubic root of the number of amphiphiles in the bilayer, thus explaining the order of magnitude difference between the simulated and the measured values. Three different methods are used to calculate a line tension coefficient of (3.5-4.0) x 10(-11) J/m, in good agreement with experimental data.

Pore nucleation in mechanically stretched bilayer membranes

We report a computer-simulation study of the free-energy barrier for the nucleation of pores in the bilayer membrane under constant stretching lateral pressure. We find that incipient pores are hydrophobic but as the lateral size of the pore nucleus becomes comparable with the molecular length, the pore becomes hydrophilic. In agreement with previous investigations, we find that the dynamical process of growth and closure of hydrophilic pores is controlled by the competition between the surface tension of the membrane and the line tension associated with the rim of the pore. We estimate the line tension of a hydrophilic pore from the shape of the computed free-energy barriers. The line tension thus computed is in a good agreement with available experimental data. We also estimate the line tension of hydrophobic pores at both macroscopic and microscopic levels. The comparison of line tensions at these two different levels indicates that the "microscopic" line tension should be carefully distinguished from the "macroscopic" effective line tension used in the theoretical analysis of pore nucleation. The overall shape of the free-energy barrier for pore nucleation shows no indication for the existence of a metastable intermediate during pore nucleation.

Pore formation in fluctuating membranes

We study the nucleation of a single pore in a fluctuating lipid membrane, specifically taking into account the membrane fluctuations, as well as the shape fluctuations of the pore. For large enough pores, the nucleation free energy is well-described by shifts in the effective membrane surface tension and the pore line tension. Using our framework, we derive the stability criteria for the various pore formation regimes. In addition to the well-known large-tension regime from the classical nucleation theory of pores, we also find a low-tension regime in which the effective line and surface tensions can change sign from their bare values. The latter scenario takes place at sufficiently high temperatures, where the opening of a stable pore of finite size is entropically favorable.

Direct imaging and mesoscale modelling of phase transitions in a nanostructured fluid

The kinetics of phase transitions is essential for understanding pattern formation in structured fluids. These fluids play a key role in the morphogenesis of biological cells, and they are very common in pharmaceutical products and plastic materials. Until now, it has not been possible to follow phase transitions in structured fluids experimentally in real time and with high spatial resolution. Previous work has relied on static images and indirect experimental evidence from spatially averaging scattering experiments. Simulating the processes with computer models is a further challenge because of the multiple time and length scales involved. Our movies based on in situ scanning force microscopy show the time sequence of the elementary steps of a phase transition in a fluid film of block copolymer from the cylinder to the perforated lamella phase. The movies validate a versatile simulation model that gives physical insight into the nature of the process. Our approach provides a means of improving the study and understanding of pattern formation processes in nanostructured fluids. We expect a significant impact on nanotechnology where block copolymers serve as self-organized templates for the synthesis of inorganic nanostructured materials.

The effect of thermal fluctuations on Schulman area elasticity

We study the elastic properties of a two-dimensional fluctuating surface whose area density is allowed to deviate from its optimal (Schulman) value. The behavior of such a surface is determined by an interplay between the area-dependent elastic energy, the curvature elasticity, and the entropy. We identify three different elastic regimes depending on the ratio Ap/As between the projected (frame) and the saturated areas. We show that thermal fluctuations modify the elastic energy of stretched surfaces (Ap/As>1), and dominate the elastic energy of compressed surfaces (Ap/As<1). When Ap approximately As the elastic energy is not much affected by the fluctuations; the frame area at which the surface tension vanishes becomes smaller than As and the area elasticity modulus increases.

An experimental and theoretical analysis of ultrasound-induced permeabilization of cell membranes

Application of ultrasound transiently permeabilizes cell membranes and offers a nonchemical, nonviral, and noninvasive method for cellular drug delivery. Although the ability of ultrasound to increase transmembrane transport has been well demonstrated, a systematic dependence of transport on ultrasound parameters is not known. This study examined cell viability and cellular uptake of calcein using 3T3 mouse cell suspension as a model system. Cells were exposed to varying acoustic energy doses at four different frequencies in the low frequency regime (20-100 kHz). At all frequencies, cell viability decreased with increasing acoustic energy dose, while the fraction of cells exhibiting uptake of calcein showed a maximum at an intermediate energy dose. Acoustic spectra under various ultrasound conditions were also collected and assessed for the magnitude of broadband noise and subharmonic peaks. While the cell viability and transport data did not show any correlation with subharmonic (f/2) emission, they correlated with the broadband noise, suggesting a dominant contribution of transient cavitation. A theoretical model was developed to relate reversible and irreversible membrane permeabilization to the number of transient cavitation events. The model showed that nearly every stage of transient cavitation, including bubble expansion, collapse, and subsequent shock waves may contribute to membrane permeabilization. For each mechanism, the volume around the bubble within which bubbles induce reversible and irreversible membrane permeabilization was determined. Predictions of the model are consistent with experimental data.

Lattice model for the kinetics of rupture of fluid bilayer membranes

We have constructed a model for the kinetics of rupture of membranes under tension, applying physical principles relevant to lipid bilayers held together by hydrophobic interactions. The membrane is characterized by the bulk compressibility (for expansion) K, the thickness 2h(t) of the hydrophobic part of the bilayer, the hydrophobicity sigma, and a parameter gamma characterizing the tail rigidity of the lipids. The model is a lattice model which incorporates strain relaxation, and considers the nucleation of pores at constant area, constant temperature, and constant particle number. The particle number is conserved by allowing multiple occupancy of the sites. An equilibrium "phase diagram" is constructed as a function of temperature and strain with the total pore surface and distribution as the order parameters. A first-order rupture line is found with increasing tension, and a continuous increase in protopore concentration with rising temperature till instability. The model explains current results on saturated and unsaturated phosphatidylcholine lipid bilayers and thicker artificial bilayers made of diblock copolymers. Pore size distributions are presented for various values of area expansion and temperature, and the fractal dimension of the pore edge is evaluated.

Adhesion of nanoparticles to vesicles: a Brownian dynamics simulation

We studied the interaction of bilayer vesicles and adhesive nanoparticles using a Brownian dynamics simulation. The nanoparticles are simple models of proteins or colloids. The adhering nanoparticle induces the morphological change of the vesicle: budding, formation of two vesicles in which only outer monolayers are connected, and fission. We also show that the nanoparticle promotes the fusion process: fusion-pore opening from a stalk intermediate, a neck-like structure that only connects outer monolayers of two vesicles. The nanoparticle bends the stalk, and induces the pore opening.

Polymersomes: tough vesicles made from diblock copolymers

Vesicles were made from amphiphilic diblock copolymers and characterized by micromanipulation. The average molecular weight of the specific polymer studied, polyethyleneoxide-polyethylethylene (EO40-EE37), is several times greater than that of typical phospholipids in natural membranes. Both the membrane bending and area expansion moduli of electroformed polymersomes (polymer-based liposomes) fell within the range of lipid membrane measurements, but the giant polymersomes proved to be almost an order of magnitude tougher and sustained far greater areal strain before rupture. The polymersome membrane was also at least 10 times less permeable to water than common phospholipid bilayers. The results suggest a new class of synthetic thin-shelled capsules based on block copolymer chemistry.

Thermally induced proliferation of pores in a model fluid membrane

The growth of thermally induced pores in a two-dimensional model fluid membrane is investigated by Monte Carlo simulation. Holes appear in the membrane via an activated process, and their subsequent growth is controlled by an edge energy per unit length or line tension. The barrier height and line tension, together with a lateral tension, are the independent parameters of the model. In the resulting phase diagram, a rupture transition separates an intact membrane from a disintegrated state. The approach to the ruptured state shows distinct regimes. Reducing the barrier height at large line tension produces multiple, quasi-independent, small holes whose behavior is dominated by their edge energy, whereas at lower line tensions shape fluctuations of the holes facilitate their coalescence into a single large hole. At a small value of line tension and large barrier height, a single hole spontaneously permeabilizes the membrane in an entropically driven phase transition. Entropy dominates pore growth for line tensions not far below those measured for artificial vesicles. Permeabilization of lipid bilayers by certain peptides involves perturbing lipid-lipid cohesive energies, and our simulations show that at small line tensions the entropy of hole shape fluctuations destroys the model membrane's stability.