Meet me on the other side: trans-bilayer modulation of a model voltage-gated ion channel activity by membrane electrostatics asymmetry

While it is accepted that biomembrane asymmetry is generated by proteins and phospholipids distribution, little is known about how electric changes manifested in a monolayer influence functional properties of proteins localized on the opposite leaflet. Herein we used single-molecule electrophysiology and investigated how asymmetric changes in the electrostatics of an artificial lipid membrane monolayer, generated oppositely from where alamethicin--a model voltage-gated ion channel--was added, altered peptide activity. We found that phlorizin, a membrane dipole potential lowering amphiphile, augmented alamethicin activity and transport features, whereas the opposite occurred with RH-421, which enhances the monolayer dipole potential. Further, the monolayer surface potential was decreased via adsorption of sodium dodecyl sulfate, and demonstrated that vectorial modification of it also affected the alamethicin activity in a predictive manner. A new paradigm is suggested according to which asymmetric changes in the monolayer dipole and surface potential extend their effects spatially by altering the intramembrane potential, whose gradient is sensed by distantly located peptides.     

Voltage-dependent lipid flip-flop induced by alamethicin.

Alamethicin appears to allow voltage-dependent lipid exchange ("flip-flop") between leaflets of a planar bilayer. In membranes with one leaflet of phosphatidyl serine and one of phosphatidyl ethanolamine, the shape of the nonactin current-voltage curve accurately reports the difference in surface potential between the two sides of the membrane. The surface potential is itself a good measure of membrane asymmetry. Alamethicin added to the bathing solutions of an asymmetric membrane does not per se reduce the membrane asymmetry, but turning on the alamethicin conductance by application of a voltage pulse does. Immediately after application of a voltage pulse, large enough to turn on the alamethicin conductance, the asymmetry of the nonactin-K+ current voltage curve decreases, in some cases, nearly to zero. During the pulse, the alamethicin conductance activates if a decrease in surface potential favors turn-on of the alamethicin conductance or inactivates if a decrease in surface potential favors turn-off of the alamethicin conductance. After the pulse, the nonactin-K+ asymmetry returns to its original value if the alamethicin conductance is not turned on. The time-course of this return allows an estimate of the diffusion constant of lipid in the planar bilayer. The value obtained is 5.1 x 10(-8) cm2/s.     

A Coarse Grained Model for a Lipid Membrane with Physiological Composition and Leaflet Asymmetry

The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.     

The Influence of Natural Lipid Asymmetry upon the Conformation of a Membrane-inserted Protein (Perfringolysin O)

Eukaryotic membrane proteins generally reside in membrane bilayers that have lipid asymmetry. However, in vitro studies of the impact of lipids upon membrane proteins are generally carried out in model membrane vesicles that lack lipid asymmetry. Our recently developed method to prepare lipid vesicles with asymmetry similar to that in plasma membranes and with controlled amounts of cholesterol was used to investigate the influence of lipid composition and lipid asymmetry upon the conformational behavior of the pore-forming, cholesterol-dependent cytolysin perfringolysin O (PFO). PFO conformational behavior in asymmetric vesicles was found to be distinct both from that in symmetric vesicles with the same lipid composition as the asymmetric vesicles and from that in vesicles containing either only the inner leaflet lipids from the asymmetric vesicles or only the outer leaflet lipids from the asymmetric vesicles. The presence of phosphatidylcholine in the outer leaflet increased the cholesterol concentration required to induce PFO binding, whereas phosphatidylethanolamine and phosphatidylserine in the inner leaflet of asymmetric vesicles stabilized the formation of a novel deeply inserted conformation that does not form pores, even though it contains transmembrane segments. This conformation may represent an important intermediate stage in PFO pore formation. These studies show that lipid asymmetry can strongly influence the behavior of membrane-inserted proteins.     

Calcium-induced inactivation of alamethicin in asymmetric lipid bilayers

This paper discusses a calcium-dependent inactivation of alamethicin- induced conductance in asymmetric lipid bilayers. The bilayers used were formed with one leaflet of phosphatidyl ethanolamine (PE) and one of phosphatidyl serine (PS). Calcium, initially confined to the neutral lipid (PE) side, can pass through the open alamethicin channel to the negative lipid (PS) side, where it can bind to the negative lipid and reduce the surface potential. Under appropriate circumstances, the voltage-dependent alamethicin conductance is thereby inactivated. We have formulated a model for this process based on the diffusion of calcium in the aqueous phases and we show that the model describes the kinetic properties of the alamethicin conductance under various circumstances. EGTA on the PS side of the membrane reduces the effects of calcium dramatically as predicted by the model.     

Simulation studies of the interaction of antimicrobial peptides and lipid bilayers

Experimental studies of a number of antimicrobial peptides are sufficiently detailed to allow computer simulations to make a significant contribution to understanding their mechanisms of action at an atomic level. In this review we focus on simulation studies of alamethicin, melittin, dermaseptin and related antimicrobial, membrane-active peptides. All of these peptides form amphipathic alpha-helices. Simulations allow us to explore the interactions of such peptides with lipid bilayers, and to understand the effects of such interactions on the conformational dynamics of the peptides. Mean field methods employ an empirical energy function, such as a simple hydrophobicity potential, to provide an approximation to the membrane. Mean field approaches allow us to predict the optimal orientation of a peptide helix relative to a bilayer. Molecular dynamics simulations that include an atomistic model of the bilayer and surrounding solvent provide a more detailed insight into peptide-bilayer interactions. In the case of alamethicin, all-atom simulations have allowed us to explore several steps along the route from binding to the membrane surface to formation of transbilayer ion channels. For those antimicrobial peptides such as dermaseptin which prefer to remain at the surface of a bilayer, molecular dynamics simulations allow us to explore the favourable interactions between the peptide helix sidechains and the phospholipid headgroups.     

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

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

Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models

Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts toward evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that 1) leaflet compositional asymmetry can induce curvature asymmetry, 2) lipids are sorted by curvature to different extents, and 3) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations and potentially interpreting the functional roles of specific lipids.     

Transbilayer effects of raft-like lipid domains in asymmetric planar bilayers measured by single molecule tracking

Cell membranes have complex lipid compositions, including an asymmetric distribution of phospholipids between the opposing leaflets of the bilayer. Although it has been demonstrated that the lipid composition of the outer leaflet of the plasma membrane is sufficient for the formation of raft-like liquid-ordered (l(o)) phase domains, the influence that such domains may have on the lipids and proteins of the inner leaflet remains unknown. We used tethered polymer supports and a combined Langmuir-Blodgett/vesicle fusion (LB/VF) technique to build asymmetric planar bilayers that mimic plasma membrane asymmetry in many ways. We show that directly supported LB monolayers containing cholesterol-rich l(o) phases are inherently unstable when exposed to water or vesicle suspensions. However, tethering the LB monolayer to the solid support with the lipid-anchored polymer 1,2-dimyristoyl phophatidylethanolamine-N-[poly(ethylene glycol)-triethoxysilane] significantly improves stability and allows for the formation of complex planar-supported bilayers that retain >90% asymmetry for 1-2 h. We developed a single molecule tracking (SPT) system for the study of lipid diffusion in asymmetric bilayers with coexisting liquid phases. SPT allowed us to study in detail the diffusion of individual lipids inside, outside, or directly opposed to l(o) phase domains. We show here that l(o) phase domains in one monolayer of an asymmetric bilayer do not induce the formation of domains in the opposite leaflet when this leaflet is composed of palmitoyl-oleoyl phosphatidylcholine and cholesterol but do induce domains when this leaflet is composed of porcine brain phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and cholesterol. The diffusion of lipids is similar in l(o) and liquid-disordered phase domains and is not affected by transbilayer coupling, indicating that lateral and transverse lipid interactions that give rise to the domain structure are weak in the biological lipid mixtures that were employed in this work.     

Alamethicin-induced current-voltage curve asymmetry in lipid bilayers.

We have examined the causes of the asymmetry of the current-voltage curve induced by addition of alamethicin to one side of a black lipid membrane. We find that the alamethicin-induced current-voltage (I-V) curve has an inherent asymmetry. If it were possible to confine all alamethicin molecules to one side of the membrane, the I-V curve would exhibit a positive branch (voltage being measured with respect to the side of the membrane trans to the alamethicin addition) of steeper logarithmic slope than the negative branch and at a lower absolute value of potential. This condition is not usually realized, however, because alamethicin can leak through the membrane, so that, except at very high alamethicin concentrations and in certain kinds of membranes, the positive branch of the current-voltage curve has the same logarithmic slope as the negative branch and appears to arise from alamethicin which diffuses from the cis to the trans side of the membrane. We develop simple quantitative models for these two cases.     

Lipid asymmetry in DLPC/DSPC-supported lipid bilayers: a combined AFM and fluorescence microscopy study

A fundamental attribute of cell membranes is transmembrane asymmetry, specifically the formation of ordered phase domains in one leaflet that are compositionally different from the opposing leaflet of the bilayer. Using model membrane systems, many previous studies have demonstrated the formation of ordered phase domains that display complete transmembrane symmetry; but there have been few reports on the more biologically relevant asymmetric membrane structures. Here we report on a combined atomic force microscopy and fluorescence microscopy study whereby we observe three different states of transmembrane symmetry in phase-separated supported lipid bilayers formed by vesicle fusion. We find that if the leaflets differ in gel-phase area fraction, then the smaller domains in one leaflet are in registry with the larger domains in the other leaflet and the system is dynamic. In a presumed lipid flip-flop process similar to Ostwald ripening, the smaller domains in one leaflet erode away whereas the large domains in the other leaflet grow until complete compositional asymmetry is reached and remains stable. We have quantified this evolution and determined that the lipid flip-flop event happens most frequently at the interface between symmetric and asymmetric DSPC domains. If both leaflets have identical area fraction of gel-phase, gel-phase domains are in registry and are static in comparison to the first state. The stability of these three DSPC domain distributions, the degree of registry observed, and the domain immobility have biological significance with regards to maintenance of lipid asymmetry in living cell membranes, communication between inner leaflet and outer leaflet, membrane adhesion, and raft mobility.     

Osmoporin OmpC forms a complex with MlaA to maintain outer membrane lipid asymmetry in Escherichia coli

Gram‐negative bacteria can survive in harsh environments in part because the asymmetric outer membrane (OM) hinders the entry of toxic compounds. Lipid asymmetry is established by having phospholipids (PLs) confined to the inner leaflet of the membrane and lipopolysaccharides (LPS) to the outer leaflet. Perturbation of OM lipid asymmetry, characterized by PL accumulation in the outer leaflet, disrupts proper LPS packing and increases membrane permeability. The multi‐component Mla system prevents PL accumulation in the outer leaflet of the OM via an unknown mechanism. Here, we demonstrate that in Escherichia coli, the Mla system maintains OM lipid asymmetry with the help of osmoporin OmpC. We show that the OM lipoprotein MlaA interacts specifically with OmpC and OmpF. This interaction is sufficient to localize MlaA lacking its lipid anchor to the OM. Removing OmpC, but not OmpF, causes accumulation of PLs in the outer leaflet of the OM in stationary phase, as was previously observed for MlaA. We establish that OmpC is an additional component of the Mla system; the OmpC‐MlaA complex may function to remove PLs directly from the outer leaflet to maintain OM lipid asymmetry. Our work reveals a novel function for the general diffusion channel OmpC in lipid transport.     

Investigations of antimicrobial peptides in planar film systems

Planar systems - monolayers and films - constitute a useful platform for studying membrane-active peptides. Here, we summarize varied approaches for studying peptide organization and peptide-lipid interactions at the air/water interface, and focus on three representative antimicrobial membrane-associated peptides—alamethicin, gramicidin, and valinomycin. Experimental data, specifically surface pressure/area isotherms and Brewster angle microscopy images, provided information on peptide association and the effects of the lipid monolayers on peptide surface organization. In general, film analysis emphasized the effects of lipid layers in promoting peptide association and aggregation at the air/water interface. Importantly, the data demonstrated that in many cases peptide domains are phase-separated within the phospholipid monolayers, suggesting that this behavior contributes to the biological actions of membrane-active antimicrobial peptides.     

Lipid chain-length dependence for incorporation of alamethicin in membranes: electron paramagnetic resonance studies on TOAC-spin labeled analogs

Alamethicin is a 19-residue hydrophobic peptide, which is extended by a C-terminal phenylalaninol but lacks residues that might anchor the ends of the peptide at the lipid-water interface. Voltage-dependent ion channels formed by alamethicin depend strongly in their characteristics on chain length of the host lipid membranes. EPR spectroscopy is used to investigate the dependence on lipid chain length of the incorporation of spin-labeled alamethicin in phosphatidylcholine bilayer membranes. The spin-label amino acid TOAC is substituted at residue positions n = 1, 8, or 16 in the sequence of alamethicin F50/5 [TOAC(n), Glu(OMe)(7,18,19)]. Polarity-dependent isotropic hyperfine couplings of the three TOAC derivatives indicate that alamethicin assumes approximately the same location, relative to the membrane midplane, in fluid diC(N)PtdCho bilayers with chain lengths ranging from N = 10-18. Residue TOAC(8) is situated closest to the bilayer midplane, whereas TOAC(16) is located farther from the midplane in the hydrophobic core of the opposing lipid leaflet, and TOAC(1) remains in the lipid polar headgroup region. Orientational order parameters indicate that the tilt of alamethicin relative to the membrane normal is relatively small, even at high temperatures in the fluid phase, and increases rather slowly with decreasing chain length (from 13 degrees to 23 degrees for N = 18 and 10, respectively, at 75 degrees C). This is insufficient for alamethicin to achieve hydrophobic matching. Alamethicin differs in its mode of incorporation from other helical peptides for which transmembrane orientation has been determined as a function of lipid chain length.     

Loss of plasma membrane lipid asymmetry can induce ordered domain (raft) formation

In some cases, lipids in one leaflet of an asymmetric artificial lipid vesicle suppress the formation of ordered lipid domains (rafts) in the opposing leaflet. Whether this occurs in natural membranes is unknown. Here, we investigated this issue using plasma membrane vesicles (PMVs) from rat leukemia RBL-2H3 cells. Membrane domain formation and order was assessed by fluorescence resonance energy transfer and fluorescence anisotropy. We found that ordered domains in PMVs prepared from cells by N-ethyl maleimide (NEM) treatment formed up to ～37°C, whereas ordered domains in symmetric vesicles formed from the extracted PMV lipids were stable up to 55°C, indicating the stability of ordered domains was substantially decreased in intact PMVs. This behavior paralleled lesser ordered domain stability in artificial asymmetric lipid vesicles relative to the corresponding symmetric vesicles, suggesting intact PMVs exhibit some degree of lipid asymmetry. This was supported by phosphatidylserine mislocalization on PMV outer leaflets as judged by annexin binding, which indicated NEM-induced PMVs are much more asymmetric than PMVs formed by dithiothreitol/paraformaldehyde treatment. Destroying asymmetry by reconstitution of PMVs using detergent dilution also showed stabilization of domain formation, even though membrane proteins remained associated with reconstituted vesicles. Similar domain stabilization was observed in artificial asymmetric lipid vesicles after destroying asymmetry via detergent reconstitution. Proteinase K digestion of proteins had little effect on domain stability in NEM PMVs. We conclude that loss of PMV lipid asymmetry can induce ordered domain formation. The dynamic control of lipid asymmetry in cells may regulate domain formation in plasma membranes.     

Making a membrane on the other side of the wall

The outer membrane (OM) of Gram-negative bacteria is positioned at the frontline of the cell's interaction with its environment and provides a barrier against influx of external toxins while still allowing import of nutrients and excretion of wastes. It is a remarkable asymmetric bilayer with a glycolipid surface-exposed leaflet and a glycerophospholipid inner leaflet. Lipid asymmetry is key to OM barrier function and several different systems actively maintain this lipid asymmetry. All OM components are synthesized in the cytosol before being secreted and assembled into a contiguous membrane on the other side of the cell wall. Work in recent years has uncovered the pathways that transport and assemble most of the OM components. However, our understanding of how phospholipids are delivered to the OM remains notably limited. Here we will review seminal works in phospholipid transfer performed some 40 years ago and place more recent insights in their context. This article is part of a Special Issue entitled: Bacterial Lipids edited by Russell E. Bishop.     

Alamethicin interaction with lipid membranes: a spectroscopic study on synthetic analogues

Alamethicin (Alm) is one of the most extensively studied membrane-active antibiotic peptides, but several aspects of its mechanism of action are still under debate. In this study, synthetic analogues of natural Alm F50/5 (Alm-N), labeled with a 9H-fluoren-9-yl group at the N- (F-Alm) or C-terminus (Alm-F), were employed to investigate the position and orientation of this peptide in the membrane environment. Depth-dependent fluorescence quenching and polarized ATR-FT-IR experiments demonstrated that, in the absence of a transmembrane potential, Alm inserts its N-terminus into the membrane, while the C-terminus is exposed to the outer aqueous phase. We also found that the peptaibol populates different orientations with respect to the membrane normal. Furthermore, fluorescence resonance-energy transfer (FRET) indicated that no peptide translocation to the inner leaflet of lipid bilayers occurs. The mechanism of action of Alm is discussed on the basis of these findings. Two other Alm analogues, Alm-P and Alm-S, were exploited to investigate the role of specific Alm residues in terms of membrane-perturbing activity. Substitution of two or three Gln (E) residues (the only polar amino acids in the alamethicin sequence) by gamma-methyl glutamate (Glu(OMe)) residues induced marked variations in the aggregation and partition behaviors of the peptaibols, which, in turn, modulate their membrane activity. In particular, substitution of Gln(18) and Gln(19) caused a six-fold increase in membrane-perturbing activity, thus demonstrating that these residues are not essential for the stabilization of Alm pores.     

Cholesterol asymmetry in synaptic plasma membranes

Lipids are essential for the structural and functional integrity of membranes. Membrane lipids are not randomly distributed but are localized in different domains. A common characteristic of these membrane domains is their association with cholesterol. Lipid rafts and caveolae are examples of cholesterol enriched domains, which have attracted keen interest. However, two other important cholesterol domains are the exofacial and cytofacial leaflets of the plasma membrane. The two leaflets that make up the bilayer differ in their fluidity, electrical charge, lipid distribution, and active sites of certain proteins. The synaptic plasma membrane (SPM) cytofacial leaflet contains over 85% of the total SPM cholesterol as compared with the exofacial leaflet. This asymmetric distribution of cholesterol is not fixed or immobile but can be modified by different conditions in vivo: (i) chronic ethanol consumption; (ii) statins; (iii) aging; and (iv) apoE isoform. Several potential candidates have been proposed as mechanisms involved in regulation of SPM cholesterol asymmetry: apoE, low-density lipoprotein receptor, sterol carrier protein-2, fatty acid binding proteins, polyunsaturated fatty acids, P-glycoprotein and caveolin-1. This review examines cholesterol asymmetry in SPM, potential mechanisms of regulation and impact on membrane structure and function.     

Antimicrobial and cell-penetrating peptides induce lipid vesicle fusion by folding and aggregation

According to their distinct biological functions, membrane-active peptides are generally classified as antimicrobial (AMP), cell-penetrating (CPP), or fusion peptides (FP). The former two classes are known to have some structural and physicochemical similarities, but fusogenic peptides tend to have rather different features and sequences. Nevertheless, we found that many CPPs and some AMPs exhibit a pronounced fusogenic activity, as measured by a lipid mixing assay with vesicles composed of typical eukaryotic lipids. Compared to the HIV fusion peptide (FP23) as a representative standard, all designer-made peptides showed much higher lipid-mixing activities (MSI-103, MAP, transportan, penetratin, Pep1). Native sequences, on the other hand, were less fusogenic (magainin 2, PGLa, gramicidin S), and pre-aggregated ones were inactive (alamethicin, SAP). The peptide structures were characterized by circular dichroism before and after interacting with the lipid vesicles. A striking correlation between the extent of conformational change and the respective fusion activities was found for the series of peptides investigated here. At the same time, the CD data show that lipid mixing can be triggered by any type of conformation acquired upon binding, whether α-helical, β-stranded, or other. These observations suggest that lipid vesicle fusion can simply be driven by the energy released upon membrane binding, peptide folding, and possibly further aggregation. This comparative study of AMPs, CPPs, and FPs emphasizes the multifunctional aspects of membrane-active peptides, and it suggests that the origin of a peptide (native sequence or designer-made) may be more relevant to define its functional range than any given name.     

Formation of asymmetric vesicles via phospholipase D-mediated transphosphatidylation

Most biomembranes have an asymmetric structure with regard to phospholipid distribution between the inner and outer leaflets of the lipid bilayers. Control of the asymmetric distribution plays a pivotal role in several cellular functions such as intracellular membrane fusion and cell division. The mechanism by which membrane asymmetry and its alteration function in these transformation processes is not yet clear. To understand the significance of membrane asymmetry on trafficking and metabolism of intracellular vesicular components, a system that experimentally reproduces the asymmetric nature of biomembranes is essential. Here, we succeeded in obtaining asymmetric vesicles by means of transphosphatidylation reactions with phospholipase D (PLD), which acts exclusively on phosphatidylcholine (PC) present in the outer leaflet of vesicles. By treating PC vesicles with PLD in the presence of 1.7 M serine and 0.3 M ethanolamine, we obtained asymmetric vesicles that are topologically similar to intracellular vesicles containing phosphatidylserine and phosphatidylethanolamine in the cytosolic leaflet. PLD and other unwanted compounds could be removed by trypsin digestion followed by dialysis. Our established technique has a great advantage over conventional methods in that asymmetric vesicles can be provided at high yield and high efficiency, which is requisite for most physicochemical assays.