Stability of Asymmetric Lipid Bilayers Assessed by Molecular Dynamics Simulations

Self-Assembly of Mycolic Acid in Water: Monolayer or Bilayer

The enduring pathogenicity of Mycobacterium tuberculosis can be attributed to its lipid-rich cell wall, with mycolic acids (MAs) being a significant constituent. Different MAs' fluidity and structural adaptability within the bacterial cell envelope significantly influence their physicochemical properties, operational capabilities, and pathogenic potential. Therefore, an accurate conformational representation of various MAs in aqueous media can provide insights into their potential role within the intricate structure of the bacterial cell wall. We have carried out MD simulations of MAs in an aqueous solution and shed light on various structural properties such as thickness, order parameters, area-per-MAs, conformational changes, and principle component (PC) in the single-component and mixture MAs monolayer. The different conformational populations in the monolayer were estimated using the distance-based analysis between the function groups represented as W, U, and Z conformations that lead to the fold of the MAs chain in the monolayer. Additionally, we have also simulated the mixture of alpha-MA (α-MA or AMA), methoxy-MA (MMA), and keto-MA (KMA) with 50.90% AMA, 36.36% MMA, and 12.72% KMA composition. The thickness of the MAs monolayer was observed to range from 5 to 7 nm with an average 820 kg/m3 density for α-MA, MMA, and KMA quantitative agreement with experimental results. The mero chain (long chain), consisting of a functional group at the proximal and distal positions, tends to fold and exhibit a more disordered phase than the short chain. The keto-MA showed the greatest WUZ total conformations (35.32%) with decreasing trend of eZ > eU > aU > aZ folds in both single component and mixture. Our results are in quantitative agreement with the experimental observations. The sZ folds show the lowest conformational probability in monolayer assembly (0.75% in a single component and 1.1% in a mixture). However, eU and aU folds are most probable for AMA and MMA. One striking observation is the abundance of MA conformers beyond the known WUZ convention because of the wide range distribution of intramolecular distances and change in dihedral angles. From a thermodynamic perspective, all mycolic acid monolayers in this study within the microsecond-long simulation, MA molecules self-assembled, and the self-assembled monolayer was found to be stable. The conformation of MAs corresponding to lower free energy minima in the monolayer gives rise to tighter packing and a highly dense self-assembly. Such a highly packed assembly shows higher resistance for drug permeability. Therefore, we concluded that the monolayer formed by AMA will be more densely packed and may cause more resistance for the drug molecules.

Molecular Dynamics Investigation into CerENP's Effect on the Lipid Matrix of Stratum Corneum

The extracellular lipid matrix in the stratum corneum (SC) plays a critical role in skin barrier functionality, comprising three primary components: ceramides, cholesterol, and free fatty acids. The diverse ceramides, differentiated by molecular structures such as hydroxylations and varying chain lengths, are essential for the lipid matrix's structural integrity. Recently, a new subclass of ceramide, 1-O-acylceramide NP (CerENP), has been identified; however, its precise role in the lipid matrix of the SC is still elusive. Herein, we investigate the role of CerENP on the structure and permeability of the SC using molecular dynamics simulations. Our findings indicate that CerENP contributes to a compact lipid matrix in the lateral dimension of our SC model with a repeat distance of about 13 nm. Additionally, ethanol permeability assessments show that CerENP effectively reduces molecular penetration through the lipid matrix. This study provides an insight into the role of a new subclass of ceramide in the SC, enhancing our understanding of skin structure and the mechanisms behind barrier dysfunction in skin diseases.

Methyl-β-cyclodextrin asymmetrically extracts phospholipid from bilayers, granting tunable control over differential stress in lipid vesicles

Lipid asymmetry - that is, a nonuniform lipid distribution between the leaflets of a bilayer - is a ubiquitous feature of biomembranes and is implicated in several cellular phenomena. Differential tension - that is, unequal lateral monolayer tensions comparing the leaflets of a bilayer- is closely associated with lipid asymmetry underlying these varied roles. Because differential tension is not directly measurable in combination with the fact that common methods to adjust this quantity grant only semi-quantitative control over it, a detailed understanding of lipid asymmetry and differential tension are impeded. To overcome these challenges, we leveraged reversible complexation of phospholipid by methyl-β-cyclodextrin (mbCD) to tune the direction and magnitude of lipid asymmetry in synthetic vesicles. Lipid asymmetry generated in our study induced (i) vesicle shape changes and (ii) gel-liquid phase coexistence in 1-component vesicles. By applying mass-action considerations to interpret our findings, we discuss how this approach provides access to phospholipid thermodynamic potentials in bilayers containing lipid asymmetry (which are coupled to the differential tension of a bilayer). Because lipid asymmetry yielded by our approach is (i) tunable and (ii) maintained over minute to hour timescales, we anticipate that this approach will be a valuable addition to the experimental toolbox for systematic investigation into the biophysical role(s) of lipid asymmetry (and differential tension).

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

The rise of FTIR spectroscopy in the characterization of asymmetric lipid membranes

In contrast to symmetric unilamellar liposomes (sLUVs) prepared from a mixture of different lipids, asymmetric ones (aLUVs) with different lipid composition in the inner and outer membrane leaflets are more suitable model systems of eukaryotic plasma membranes. However, apart from the challenging preparation of asymmetric liposomes and small amounts of obtained asymmetric unilamellar liposomes (aLUVs), a major drawback is the qualitative characterization of asymmetry, as each of the techniques used so far has certain limitations. In this regard, we prepared aLUVs composed dominantly of DPPC(out)/DPPS(in) lipids and, along with 1H NMR and DSC characterization, we showed for the first time how FTIR spectroscopy can be used in the presence of (a)symmetry between DPPC/DPPS lipid bilayers. Using second derivative FTIR spectra we demonstrated not only that the hydration of lipids glycerol backbone and choline moiety of DPPC differs in s/aLUVs, but in addition that the lateral interactions between hydrocarbon chains during the phase change display different trend in s/aLUVs. Molecular dynamics simulations confirmed different chain ordering and packing between s/a bilayers, with a significant influence of temperature, i.e. membrane phase.

Estimation of negative membrane tension in lipid bilayers and its effect on antimicrobial peptide magainin 2-induced pore formation

Positive membrane tension in the stretched plasma membrane of cells and in the stretched lipid bilayer of vesicles has been well analyzed quantitatively, whereas there is limited quantitative information on negative membrane tension in compressed plasma membranes and lipid bilayers. Here, we examined negative membrane tension quantitatively. First, we developed a theory to describe negative membrane tension by analyzing the free energy of lipid bilayers to obtain a theoretical equation for negative membrane tension. This allowed us to obtain an equation describing the negative membrane tension (σosm) for giant unilamellar vesicles (GUVs) in hypertonic solutions due to negative osmotic pressure (Π). Then, we experimentally estimated the negative membrane tension for GUVs in hypertonic solutions by measuring the rate constant (kr) of rupture of the GUVs induced by the constant tension (σex) due to an external force as a function of σex. We found that larger σex values were required to induce the rupture of GUVs under negative Π compared with GUVs in isotonic solution and quantitatively determined the negative membrane tension induced by Π (σosm) by the difference between these σex values. At small negative Π, the experimental values of negative σosm agree with their theoretical values within experimental error, but as negative Π increases, the deviation increases. Negative tension increased the stability of GUVs because higher tensions were required for GUV rupture, and the rate constant of antimicrobial peptide magainin 2-induced pore formation decreased.

Seeing the Membrane from Both Sides Now: Lipid Asymmetry and Its Strange Consequences

Almost all biomembranes are constructed as lipid bilayers and, in almost all of these, the two opposing monolayers (leaflets) have distinct lipid compositions. This lipid asymmetry arises through the concerted action of a suite of energy-dependent enzymes that maintain living bilayers in a far-from-equilibrium steady-state. Recent discoveries reveal that lipid compositional asymmetry imparts biophysical asymmetries and that this dualistic organization may have major consequences for cellular physiology. Importantly, while transbilayer asymmetry appears to be an essential, near-ubiquitous characteristic of biological membranes, it has been challenging to reproduce in reconstituted or synthetic systems. Although recent methodological developments have overcome some critical challenges, it remains difficult to extrapolate results from available models to biological systems. Concurrently, there are few experimental approaches for targeted, controlled manipulation of lipid asymmetry in living cells. Thus, the biophysical and functional consequences of membrane asymmetry remain almost wholly unexplored. This perspective summarizes the current state of knowledge and highlights emerging themes that are beginning to make inroads into the fundamental question of why life tends toward asymmetry in its bilayers.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Complex electrostatic effects on the selectivity of membrane-permeabilizing cyclic lipopeptides

Cyclic lipopeptides (CLiPs) have many biological functions, including the selective permeabilization of target membranes, and technical and medical applications. We studied the anionic CLiP viscosin from Pseudomonas along with a neutral analog, pseudodesmin A, and the cationic viscosin-E2K to better understand electrostatic effects on target selectivity. Calcein leakage from liposomes of anionic phosphatidylglycerol (PG) and phosphatidylethanolamine (PE) is measured in comparison with net-neutral phosphatidylcholine by time-resolved fluorescence. By contrast to the typical selectivity of cationic peptides against anionic membranes, we find viscosin more active against PG/PE at 30 μM lipid than viscosin-E2K. At very low lipid concentration, the selectivity is reversed. An equi-activity analysis reveals the reciprocal partition coefficients, 1/K, and the CLiP-to-lipid mole ratio within the membrane as leakage after 1 h reaches 50%, Re50. As expected, 1/K to PG/PE is much lower (higher affinity) for viscosin-E2K (3 μM) than viscosin (15 μM). However, the local damage to the PG/PE membrane caused by a viscosin molecule is much stronger than that of viscosin-E2K. This can be explained by the strong membrane expansion due to PG/viscosin repulsion inducing asymmetry stress between the two leaflets and, ultimately, transient limited leakage at Re50 = 0.08. PG/viscosin-E2K attraction opposes expansion and leakage starts only as the PG charges in the outer leaflet are essentially compensated by the cationic peptide (Re50 = 0.32). In the high-lipid regime (at lipid concentrations cL ≫ 1/K), virtually all CLiP is membrane bound anyway and Re50 governs selectivity, favoring viscosin. In the low-lipid regime at cL ≪ 1/K, virtually all CLiP is in solution, 1/K becomes important and the "cation attacks anionic membrane" selectivity gets restored. Overall, activity and selectivity data can only properly be interpreted if the lipid regime is known and predictions for other lipid concentrations or cell counts require knowledge of 1/K and Re50.

Leaky membrane fusion: an ambivalent effect induced by antimicrobial polycations

Both antimicrobial peptides and their synthetic mimics are potential alternatives to classical antibiotics. They can induce several membrane perturbations including permeabilization. Especially in model studies, aggregation of vesicles by such polycations is often reported. Here, we show that unintended vesicle aggregation or indeed fusion can cause apparent leakage in model studies that is not possible in most microbes, thus potentially leading to misinterpretations. The interactions of a highly charged and highly selective membrane-active polycation with negatively charged phosphatidylethanolamine/phosphatidylglycerol (PE/PG) vesicles are studied by a combination of biophysical methods. At low polycation concentrations, apparent vesicle aggregation was found to involve exchange of lipids. Upon neutralization of the negatively charged vesicles by the polycation, full fusion and leakage occurred and leaky fusion is suspected. To elucidate the interplay of leakage and fusion, we prevented membrane contacts by decorating the vesicles with PEG-chains. This inhibited fusion and also leakage activity. Leaky fusion is further corroborated by increased leakage with increasing likeliness of vesicle-vesicle contacts. Because of its similar appearance to other leakage mechanisms, leaky fusion is difficult to identify and might be overlooked and more common amongst polycationic membrane-active compounds. Regarding biological activity, leaky fusion needs to be carefully distinguished from other membrane permeabilization mechanisms, as it may be less relevant to bacteria, but potentially relevant for fungi. Furthermore, leaky fusion is an interesting effect that could help in endosomal escape for drug delivery. A comprehensive step-by-step protocol for membrane permeabilization/vesicle leakage using calcein fluorescence lifetime is provided in the ESI.

Tracking the Endosomal Escape: A Closer Look at Calcein and Related Reporters

Crossing the cellular membrane and delivering active pharmaceuticals or biologicals into the cytosol of cells is an essential step in the development of nanomedicines. One of the most important intracellular processes regarding the cellular uptake of biologicals is the endolysosomal pathway. Sophisticated nanocarriers have been developed overcoming a major hurdle, the endosomal entrapment, and delivering their cargo to the required site of action. In parallel, in vitro assays have been established analyzing the performance of these nanocarriers. Among them, the release of the membrane-impermeable dye calcein has become a popular and straightforward method. It is accessible for most researchers worldwide, allows for rapid conclusions about the release potential, and enables the study of release mechanisms. This review is intended to provide an overview and guidance for scientists applying the calcein release assay. It comprises a survey of several applications in the study of endosomal escape, considerations of potential pitfalls, challenges and limitations of the assay, and a brief summary of complementary methods. Based on this review, we hope to encourage further research groups to take advantage of the calcein release assay for their own purposes and help to create a database for more efficient cross-correlations between nanocarriers. This article is protected by copyright. All rights reserved.

Interplay of Fusion, Leakage, and Electrostatic Lipid Clustering: Membrane Perturbations by a Hydrophobic Antimicrobial Polycation

Membrane active compounds are able to induce various types of membrane perturbations. Natural or biomimetic candidates for antimicrobial treatment or drug delivery scenarios are mostly designed and tested for their ability to induce membrane permeabilization, also termed leakage. Furthermore, the interaction of these usually cationic amphiphiles with negatively charged vesicles often causes colloidal instability leading to vesicle aggregation or/and vesicle fusion. We show the interplay of these modes of membrane perturbation in mixed phosphatidyl glycerol (PG)/phosphatidyl ethanolamine (PE) by the statistical copolymer MM:CO comprising, both, charged and hydrophobic subunits. MM:CO is a representative of partially hydrophobic, highly active, but less selective antimicrobial polycations. Cryo-electron microscopy indicates vesicle fusion rather than vesicle aggregation upon the addition of MM:CO to negatively charged PG/PE (1:1) vesicles. In a combination of fluorescence-based leakage and fusion assays, there is support for membrane permeabilization and pronounced vesicle fusion activity as distinct effects. To this end, membrane fusion and aggregation were prevented by including lipids with polyethylene glycol attached to their head groups (PEG-lipids). The leakage activity of MM:CO is very similar in the absence and presence of PEG-lipids. Vesicle aggregation and fusion however are largely suppressed. This strongly suggests that MM:CO induces leakage by asymmetric packing stress because of hydrophobically driven interactions which could lead to leakage. As a further membrane perturbation effect, MM:CO causes lipid clustering in model vesicles. We address potential artifacts and misinterpretations of experiments characterizing leakage and fusion. Additional to the leakage activity, the pronounced fusogenic activity of the polymer and potentially of many other similar compounds likely has implications for antimicrobial activity and beyond.

Leaflet Asymmetry Modeling in the Lipid Composition of Escherichia coli Cytoplasmic Membranes

Lipid composition asymmetry between leaflets is important to cell function and plays a key role in the "positive inside" rule in transmembrane proteins. In this work, Escherichia coli inner plasma membrane models reflecting this asymmetry have been investigated at the early-log and stationary stages during the bacterial lifecycle using all-atom molecular dynamics simulations. The CHARMM36 lipid force field is used, and selected membrane properties are tested for variations between two leaflets and whole membranes. Our models include bacterial lipids with a cyclopropane moiety on the sn-2 acyl chain in the stationary membrane model. The PE/PG ratio for two leaflets reflects the "positive inside" rule of membrane proteins, set to 6.8 and 2.8 for the inner and outer leaflets of the two models, respectively. We are the first to model leaflet asymmetry in the lipid composition of E. coli cytoplasmic membranes and observe the effect on membrane properties in leaflets and whole membranes. Specifically, our results show that for the stationary phase bilayer, the surface area per lipid (SA/lipid) is larger, the thickness (2D_C and D_B) is smaller, the tilt angle is larger, the tilt modulus is smaller, and the deuterium order parameters (S_CD) of sn-1 and sn-2 tails are lower, compared to the early-log stage. Moreover, the stationary stage bilayer has a positive spontaneous curvature, while the early-log stage has a near flat spontaneous curvature. For leaflet asymmetry, the inner leaflet has a larger SA/lipid, a smaller thickness, a smaller elastic tilt modulus (a larger tilt angle), and lower S_CD, compared to the outer leaflet in both stages. Moreover, an asymmetric membrane involves a lipid tilt and a lateral extension, varying from a reference state of a pre-equilibrium membrane. This work encourages a more profound exploration of leaflet asymmetry in various other membrane models and how this might affect the structure and function of membrane-associated peptides and proteins.

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

The Optimal Lipid Chain Length of a Membrane-Permeabilizing Lipopeptide Results From the Balance of Membrane Partitioning and Local Damage

Pseudodesmin A (PSD) is a cyclic lipodepsipeptide produced by Pseudomonas that kills certain bacteria at MIC_1/2 in the single micromolar range, probably by permeabilizing their cellular membranes. Synthetic PSD variants, where the native decanoic (C10) acyl chain is varied in length from C4 to C8 and C12 to C14 carbons, were described to be not or less active against a panel of gram-positive strains, as compared to native PSD-C10. Here, we test the membrane-permeabilizing activity of PSD-C4 through PSD-C14 in terms of calcein release from liposomes, which is characterized in detail by the fluorescence-lifetime based leakage assay. Antagonistic concentrations and their chain length dependence agree well for liposome leakage and antimicrobial activity. The optimal chain length is governed by a balance between membrane partitioning (favoring longer chains) and the local perturbation or “damage” inflicted by a membrane-bound molecule (weakening for longer chains). Local perturbation, in turn, may involve at least two modes of action. Asymmetry stress between outer and inner leaflet builds up as the lipopeptides enter the outer leaflet and when it reaches a system-specific stability threshold, it causes a transient membrane failure that allows for the flip of some molecules from the outer to the inner leaflet. This cracking-in may be accompanied by transient, incomplete leakage from the aqueous cores of the liposomes observed, typically, for some seconds or less. The mismatch of the lipopeptide with the lipid leaflet geometry, expressed for example in terms of a spontaneous curvature, has two effects. First, it affects the threshold for transient leakage as described. Second, it controls the rate of equilibrium leakage proceeding as the lipopeptide has reached sufficient local concentrations in both leaflets to form quasi-toroidal defects or pores. Both modes of action, transient and equilibrium leakage, synergize for intermediate chain lengths such as the native, i.e., for PSD-C10. These mechanisms may also account for the reported chain-length dependent specificities of antibiotic action against the target bacteria.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

Activating mechanosensitive channels embedded in droplet interface bilayers using membrane asymmetry

Droplet microcompartments linked by lipid bilayers show great promise in the construction of synthetic minimal tissues. Central to controlling the flow of information in these systems are membrane proteins, which can gate in response to specific stimuli in order to control the molecular flux between membrane separated compartments. This has been demonstrated with droplet interface bilayers (DIBs) using several different membrane proteins combined with electrical, mechanical, and/or chemical activators. Here we report the activation of the bacterial mechanosensitive channel of large conductance (MscL) in a dioleoylphosphatidylcholine:dioleoylphosphatidylglycerol DIB by controlling membrane asymmetry. We show using electrical measurements that the incorporation of lysophosphatidylcholine (LPC) into one of the bilayer leaflets triggers MscL gating in a concentration-dependent manner, with partial and full activation observed at 10 and 15 mol% LPC respectively. Our findings could inspire the design of new minimal tissues where flux pathways are dynamically defined by lipid composition.

Hidden complexity in membrane permeabilization behavior of antimicrobial polycations

There are diverse membrane permeabilization behaviors of antimicrobial polycations in zwitterionic or charged vesicles; different mechanisms may occur over time.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Simulations of Asymmetric Membranes Illustrate Cooperative Leaflet Coupling and Lipid Adaptability

Biological membranes are composed of lipid bilayers that are often asymmetric with regards to the lipid composition and/or aqueous solvent they separate. Studying lipid asymmetry both experimentally and computationally is challenging. Molecular dynamics simulations of lipid bilayers with asymmetry are difficult due to finite system sizes and time scales accessible to simulations. Due to the very slow flip-flop rate for phospholipids, one must first choose how many lipids are on each side of the bilayer, but the resulting bilayer may be unstable (or metastable) due to differing tensile and compressive forces between leaflets. Here we use molecular dynamics simulations to investigate a number of different asymmetric membrane systems, both with atomistic and coarse-grained models. Asymmetries studied include differences in number of lipids, lipid composition (unsaturated and saturated tails and different headgroups), and chemical gradients between the aqueous phases. Extensive analysis of the bilayers' properties such as area per lipid, density, and lateral pressure profiles are used to characterize bilayer asymmetry. We also address how cholesterol (which flip-flops relatively quickly) influences membrane asymmetries. Our results show how each leaflet is influenced by the other and can mitigate the structural changes to the bilayer overall structure. Cholesterol can respond to changes in bilayer asymmetry to alleviate some of the effect on the bilayer structure, but that will alter its leaflet distribution, which in turn affects its chemical potential. Ionic imbalances are shown to have a modest change in bilayer structure, despite large changes in the electrostatic potential. Bilayer asymmetry can also induce a modest electrostatic potential across the membrane. Our results highlight the importance of membrane asymmetry on bilayer properties, the influence of lipid headgroups, tails and cholesterol on asymmetry, and the ability of lipids to adapt to different environments.

Spontaneous Curvature, Differential Stress, and Bending Modulus of Asymmetric Lipid Membranes

Lipid bilayers can exhibit asymmetric states, in which the physical characteristics of one leaflet differ from those of the other. This most visibly manifests in a different lipid composition, but it can also involve opposing lateral stresses in each leaflet that combine to an overall vanishing membrane tension. Here, we use theoretical modeling and coarse-grained simulation to explore the interplay between a compositional asymmetry and a nonvanishing differential stress. Minimizing the total elastic energy leads to a preferred spontaneous curvature that balances torques due to both bending moments and differential stress, with sometimes unexpected consequences. For instance, asymmetric flat bilayers, whose specific areas in each leaflet are matched to those of corresponding tensionless symmetric flat membranes, still exhibit a residual differential stress because the conditions of vanishing area strain and vanishing bending moment differ. We also measure the curvature rigidity of asymmetric bilayers and find that a sufficiently strong differential stress, but not compositional asymmetry alone, can increase the bending modulus. The likely cause is a stiffening of the compressed leaflet, which appears to be related to its gel transition but not identical with it. We finally show that the impact of cholesterol on differential stress depends on the relative strength of elastic and thermodynamic driving forces: if cholesterol solvates equally well in both leaflets, it will redistribute to cancel both leaflet tensions almost completely, but if its partitioning free energy prefers one leaflet over the other, the resulting distribution bias may even create differential stress. Because cells keep most of their lipid bilayers in an asymmetric nonequilibrium steady state, our findings suggest that biomembranes are elastically more complex than previously thought: besides a spontaneous curvature, they might also exhibit significant differential stress, which could strongly affect their curvature energetics.

Quantified Membrane Permeabilization Indicates the Lipid Selectivity of Membrane-Active Antimicrobials

Most antimicrobial peptides (AMPs) and their synthetic mimics (SMAMPs) are thought to act by permeabilizing cell membranes. For antimicrobial therapy, selectivity for pathogens over mammalian cells is a key requirement. Understanding membrane selectivity is thus essential for designing AMPs and SMAMPs to complement classical antibiotics in the future. This study focuses on membrane permeabilization induced by SMAMPs and their selectivity for membranes with different lipid compositions. We measure release and fluorescence lifetime of a self-quenching dye in lipid vesicles. Apart from the dose-response, we quantify the strength of individual leakage events, and, employing cumulative kinetics, categorize permeabilization behavior. We propose that differing selectivities in a series of SMAMPs arise from a combination of the effect of the antimicrobial agent and the susceptibility of the membrane (with a given lipid composition) for certain types of leakage behavior. The unselective and hemolytic SMAMP is found to act mainly by the asymmetry stress mechanism, mediated by hydrophobic insertion of SMAMPs into lipid layers. The more selective SMAMPs induced leakage events occurring stochastically over several hours. Lipid intrinsic properties might additionally amplify the efficiency of leakage events. Leakage behavior changes with both the design of the SMAMP and the lipid composition of the membrane. Understanding how leakage behavior contributes to the selectivity and activity of antimicrobial agents will aid the design and screening of antimicrobials. An understanding of the underlying processes facilitates the comparison of membrane permeabilization across in vitro and in vivo assays.

Bilayer Membranes with Frequent Flip-Flops Have Tensionless Leaflets

Biomembranes are built up from lipid bilayers with two leaflets that typically differ in their lipid composition. Each lipid molecule stays within one leaflet of the bilayer before it undergoes a transition, or flip-flop, to the other leaflet. The corresponding flip-flop times are very different for different lipid species and vary over several orders of magnitude. Here, we use molecular dynamics simulations to elucidate the consequences of this separation of time scales for compositionally asymmetric bilayers. We first study bilayers with two lipid components that do not undergo flip-flops on the accessible time scales. In such a situation, one must distinguish a bilayer state in which both leaflets have the same preferred area from another state in which each leaflet is tensionless. However, when we add a third lipid component that undergoes frequent flip-flops, the bilayer relaxes toward the state with tensionless leaflets, not to the state with equal preferred leaflet areas. Furthermore, we show that bilayers with compositional asymmetry acquire a significant spontaneous curvature even if both leaflets are tensionless. Our results can be extended to lipid bilayers with a large number of lipid components provided at least one of these components undergoes frequent flip-flops. For cellular membranes containing lipid pumps, the leaflet tensions also depend on the rates of protein-induced flip-flops.

Spontaneous and Stress-Induced Pore Formation in Membranes: Theory, Experiments and Simulations

The large plasticity, dynamics and adaptability of biological membranes allow different modes of intrinsic and inducible permeability. These phenomena are of physiological importance for a number of natural functions related to cell death and can also be manipulated artificially for practical purposes like gene transfer, drug delivery, prevention of infections or anticancer therapy. For these advances to develop in a controllable and specific way, we need a sufficient understanding of the membrane permeability phenomena. Since the formulation of early concepts of pore formation, there has been an enormous effort to describe membrane permeability by using theory, simulations and experiments. A major breakthrough has come recently through theoretical developments that allow building continuous trajectories of pore formation both in the absence and presence of stress conditions. The new model provides a coherent quantitative view of membrane permeabilization, useful to test the impact of known lipid properties, make predictions and postulate specific pore intermediates that can be studied by simulations. For example, this theory predicts unprecedented dependencies of the line tension on the pore radius and on applied lateral tension which explain previous puzzling results. In parallel, important concepts have also come from molecular dynamics simulations, of which the role of water for membrane permeabilization is of special interest. These advances open new challenges and perspectives for future progress in the study of membrane permeability, as experiments and simulations will need to test the theoretical predictions, while theory achieves new refinements that provide a physical ground for observations.

Improved Intracellular Delivery of Polyarginine Peptides with Cargoes

Complementary to endocytosis, cell-penetrating peptides (CPPs) at high concentrations can penetrate the cell membrane in a direct way, which further makes CPPs popular candidates for delivering therapeutic or diagnostic agents. Although featured as rapid uptake, the translocation efficiency and potential toxicity of the direct penetration are usually affected by cargoes, which is still unclear. Here, using coarse-grained molecular dynamics simulations, we show that the polyarginine (R₈) peptides penetrate the membrane through a water pore in the membrane, and the transmembrane efficiency is improved by conjugating to small nanoparticles (NPs) with proper linkers. It can be attributed to both the extension of the lifetime of the water pore by the NPs and outward diffusion of negative lipids in the asymmetry membrane, which induces the surrounding R₈-NP conjugates to the water pore before it is closed. The translocation efficiency is closely related to the length of the linkers, and it gets the maximum value when the length of the linkers is around half of the membrane thickness. Overlong linkers not only decrease the transmembrane efficiency because of the blockage of NPs in the water pore but may also cause cytotoxicity because of the unclosed water pore. The results provide insights into the internalization of CPPs and facilitate the design of CPP and drug conjugates with high efficiency and low toxicity.

Molecular Mechanism of Polycation-Induced Pore Formation in Biomembranes

Polycations are an attractive class of macromolecules with promising applications as drug/gene carriers and biocides. The chemical structure and concentration of a polycation determine its interaction with cellular membranes and, hence, are crucial parameters for designing efficient nontoxic polycations. However, the interaction of polycations with biomembranes at the molecular level and the corresponding free-energy landscape is not well understood. In this work, we investigate the molecular mechanism of interaction between a strong polycation substituted with alkyl moieties and zwitterionic membranes via long-time-scale all-atom molecular dynamics simulations and free-energy calculations combined with Langmuir monolayer, atomic force microscopy, and calcein-release experimental measurements. We found that the membrane activity of the polycation and its ability to induce pores in the membranes can be attributed to the polycation-induced changes in the bilayer organization, such as reduced membrane thickness, increased disorder of the acyl chains, reduced packing, and electrostatic field gradients between membrane leaflets. These changes facilitate the penetration of water into the membrane and the formation of aqueous defects/pores. The calculated free-energy profiles indicate that the polycation lowers the nucleation barrier for pore opening and the free energy for pore formation in a concentration-dependent manner. Above the critical coverage of the membrane, the polycation nucleates spontaneous pores in zwitterionic membranes. Our work demonstrates the potential of combining enhanced sampling methods in MD simulations with experiments for a quantitative description of various events in the polycation-membrane interaction cycle, such as strong adsorption on the membrane due to hydrophobic and electrostatic interactions, and pore formation.

Accurate In Silico Modeling of Asymmetric Bilayers Based on Biophysical Principles

Technological advances in the last decade have enabled the study of ever more complex and physiologically relevant model membranes to help dispel the mystery surrounding the role of plasma membrane asymmetry in various cellular processes. The slowly accumulating body of experimental data is fueling renewed interest in and the need for computational methods to support interpretations and address a wide range of problems that are still not amenable to direct experimental study. The specific appeal of molecular dynamics simulations for this purpose lies in their ability to access information at atomic resolution, which is useful for the formulation of testable mechanistic hypotheses. But, the range of questions that can be addressed reliably with such simulations is determined by the appropriate construction and simulation of asymmetric bilayer models. One essential way to achieve this goal is to follow rigorous biophysical criteria and principles. In this context, we show that the requirement for a robust comparison between the properties of simulated asymmetric and symmetric model membranes is for the tension in each bilayer leaflet to be zero. Commonly used methods for constructing asymmetric bilayers, including matching the average areas of the leaflets from the corresponding symmetric systems, do not ensure zero leaflet tension, thus precluding physically realistic changes in the areas of the two leaflets. We present, to our knowledge, a new method for identifying the ideal lipid packing in bilayers with different leaflet compositions that achieves the zero-tension goal, and discuss the basic principles underlying the biophysically correct computational study of asymmetric membranes.

Intracellular Delivery by Membrane Disruption: Mechanisms, Strategies, and Concepts

Intracellular delivery is a key step in biological research and has enabled decades of biomedical discoveries. It is also becoming increasingly important in industrial and medical applications ranging from biomanufacture to cell-based therapies. Here, we review techniques for membrane disruption-based intracellular delivery from 1911 until the present. These methods achieve rapid, direct, and universal delivery of almost any cargo molecule or material that can be dispersed in solution. We start by covering the motivations for intracellular delivery and the challenges associated with the different cargo types-small molecules, proteins/peptides, nucleic acids, synthetic nanomaterials, and large cargo. The review then presents a broad comparison of delivery strategies followed by an analysis of membrane disruption mechanisms and the biology of the cell response. We cover mechanical, electrical, thermal, optical, and chemical strategies of membrane disruption with a particular emphasis on their applications and challenges to implementation. Throughout, we highlight specific mechanisms of membrane disruption and suggest areas in need of further experimentation. We hope the concepts discussed in our review inspire scientists and engineers with further ideas to improve intracellular delivery.

Mechanism of Initial Stage of Pore Formation Induced by Antimicrobial Peptide Magainin 2

Antimicrobial peptide magainin 2 forms pores in lipid bilayers, a property that is considered the main cause of its bactericidal activity. Recent data suggest that tension or stretching of the inner monolayer plays an important role in magainin 2-induced pore formation in lipid bilayers. Here, to elucidate the mechanism of magainin 2-induced pore formation, we investigated the effect on pore formation of asymmetric lipid distribution in two monolayers. First, we developed a method to prepare giant unilamellar vesicles (GUVs) composed of dioleoylphosphatidylglycerol (DOPG), dioleoylphosphatidylcholine (DOPC), and lyso-PC (LPC) in the inner monolayer and of DOPG/DOPC in the outer monolayer. We consider that in these GUVs, the lipid packing in the inner monolayer was larger than that in the outer monolayer. Next, we investigated the interaction of magainin 2 with these GUVs with an asymmetric distribution of LPC using the single GUV method, and found that the rate constant of magainin 2-induced pore formation, k_p, decreased with increasing LPC concentration in the inner monolayer. We constructed a quantitative model of magainin 2-induced pore formation, whereby the binding of magainin 2 to the outer monolayer of a GUV induces stretching of the inner monolayer, causing pore formation. A theoretical equation defining k_p as a function of magainin 2 surface concentration, X, reasonably explains the experimental relationship between k_p and X. This model quantitatively explains the effect on k_p of the LPC concentration in the inner monolayer. On the basis of these results, we discuss the mechanism of the initial stage of magainin 2-induced pore formation.

Systematic coarse-grained modeling of complexation between small interfering RNA and polycations

All-atom molecular dynamics simulations can provide insight into the properties of polymeric gene-delivery carriers by elucidating their interactions and detailed binding patterns with nucleic acids. However, to explore nanoparticle formation through complexation of these polymers and nucleic acids and study their behavior at experimentally relevant time and length scales, a reliable coarse-grained model is needed. Here, we systematically develop such a model for the complexation of small interfering RNA (siRNA) and grafted polyethyleneimine copolymers, a promising candidate for siRNA delivery. We compare the predictions of this model with all-atom simulations and demonstrate that it is capable of reproducing detailed binding patterns, charge characteristics, and water release kinetics. Since the coarse-grained model accelerates the simulations by one to two orders of magnitude, it will make it possible to quantitatively investigate nanoparticle formation involving multiple siRNA molecules and cationic copolymers.

In situ, Reversible Gating of a Mechanosensitive Ion Channel through Protein-Lipid Interactions

Understanding the functioning of ion channels, as well as utilizing their properties for biochemical applications requires control over channel activity. Herein we report a reversible control over the functioning of a mechanosensitive ion channel by interfering with its interaction with the lipid bilayer. The mechanosensitive channel of large conductance from Escherichia coli is reconstituted into liposomes and activated to its different sub-open states by titrating lysophosphatidylcholine (LPC) into the lipid bilayer. Activated channels are closed back by the removal of LPC out of the membrane by bovine serum albumin (BSA). Electron paramagnetic resonance spectra showed the LPC-dose-dependent gradual opening of the channel pore in the form of incrementally increasing spin label mobility and decreasing spin-spin interaction. A method to reversibly open and close mechanosensitive channels to distinct sub-open conformations during their journey from the closed to the fully open state enables detailed structural studies to follow the conformational changes during channel functioning. The ability of BSA to revert the action of LPC opens new perspectives for the functional studies of other membrane proteins that are known to be activated by LPC.

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets?

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5–10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40–160 lipids per leaflet).

Data supporting beta-amyloid dimer structural transitions and protein-lipid interactions on asymmetric lipid bilayer surfaces using MD simulations on experimentally derived NMR protein structures

This data article supports the research article entitled “Maximally Asymmetric Transbilayer Distribution of Anionic Lipids Alters the Structure and interaction with Lipids of an Amyloidogenic Protein Dimer Bound to the Membrane Surface” [1]. We describe supporting data on the binding kinetics, time evolution of secondary structure, and residue-contact maps of a surface-absorbed beta-amyloid dimer protein on different membrane surfaces. We further demonstrate the sorting of annular and non-annular regions of the protein/lipid bilayer simulation systems, and the correlation of lipid-number mismatch and surface area per lipid mismatch of asymmetric lipid membranes.

Maximally asymmetric transbilayer distribution of anionic lipids alters the structure and interaction with lipids of an amyloidogenic protein dimer bound to the membrane surface

We used molecular dynamics simulations to explore the effects of asymmetric transbilayer distribution of anionic phosphatidylserine (PS) lipids on the structure of a protein on the membrane surface and subsequent protein-lipid interactions. Our simulation systems consisted of an amyloidogenic, beta-sheet rich dimeric protein (D42) absorbed to the phosphatidylcholine (PC) leaflet, or protein-contact PC leaflet, of two membrane systems: a single-component PC bilayer and double PC/PS bilayers. The latter comprised of a stable but asymmetric transbilayer distribution of PS in the presence of counterions, with a 1-component PC leaflet coupled to a 1-component PS leaflet in each bilayer. The maximally asymmetric PC/PS bilayer had a non-zero transmembrane potential (TMP) difference and higher lipid order packing, whereas the symmetric PC bilayer had a zero TMP difference and lower lipid order packing under physiologically relevant conditions. Analysis of the adsorbed protein structures revealed weaker protein binding, more folding in the N-terminal domain, more aggregation of the N- and C-terminal domains and larger tilt angle of D42 on the PC leaflet surface of the PC/PS bilayer versus the PC bilayer. Also, analysis of protein-induced membrane structural disruption revealed more localized bilayer thinning in the PC/PS versus PC bilayer. Although the electric field profile in the non-protein-contact PS leaflet of the PC/PS bilayer differed significantly from that in the non-protein-contact PC leaflet of the PC bilayer, no significant difference in the electric field profile in the protein-contact PC leaflet of either bilayer was evident. We speculate that lipid packing has a larger effect on the surface adsorbed protein structure than the electric field for a maximally asymmetric PC/PS bilayer. Our results support the mechanism that the higher lipid packing in a lipid leaflet promotes stronger protein-protein but weaker protein-lipid interactions for a dimeric protein on membrane surfaces.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Correlating antimicrobial activity and model membrane leakage induced by nylon-3 polymers and detergents

Most antimicrobial peptides act upon target microorganisms by permeabilizing their membranes. The mode of action is often assessed by vesicle leakage experiments that use model membranes, with the assumption that biological activity correlates with the permeabilization of the lipid bilayer. The current work aims to extend the interpretation of vesicle leakage results and examine the correlation between vesicle leakage and antimicrobial activity. To this end, we used a lifetime-based leakage assay with calcein-loaded vesicles to study the membrane permeabilizing properties of a novel antifungal polymer poly-NM, two of its analogs, and a series of detergents. In conjunction, the biological activities of these compounds against Candida albicans were assessed and correlated with data from vesicle leakage. Poly-NM induces all-or-none leakage in polar yeast lipid vesicles at the polymer's MIC, 3 μg mL(-1). At this and higher concentrations, complete leakage after an initial lag time was observed. Concerted activity tests imply that this polymer acts independently of the detergent octyl glucoside (OG) for both vesicle leakage and activity against C. albicans spheroplasts. In addition, poly-NM was found to have negligible activity against zwitterionic vesicles and red blood cells. Our results provide a consistent, detailed picture of the mode of action of poly-NM: this polymer induces membrane leakage by electrostatic lipid clustering. In contrast, poly-MM:CO, a nylon-3 polymer comprised of both cationic and hydrophobic segments, seems to act by a different mechanism that involves membrane asymmetry stress. Vesicle leakage for this polymer is transient (limited to <100%) and graded, non-specific among zwitterionic and polar yeast lipid vesicles, additive with detergent action, and correlates poorly with biological activity. Based on these results, we conclude that comprehensive leakage experiments can provide a detailed description of the mode of action of membrane permeabilizing compounds. Without this thorough approach, it would have been logical to assume that the two nylon-3 polymers we examined act via similar mechanisms; it is surprising that their mechanisms are so distinct. Some, but not all mechanisms of vesicle permeabilization allow for antimicrobial activity.

The activation mode of the mechanosensitive ion channel, MscL, by lysophosphatidylcholine differs from tension-induced gating

One of the best-studied mechanosensitive channels is the mechanosensitive channel of large conductance (MscL). MscL senses tension in the membrane evoked by an osmotic down shock and directly couples it to large conformational changes leading to the opening of the channel. Spectroscopic techniques offer unique possibilities to monitor these conformational changes if it were possible to generate tension in the lipid bilayer, the native environment of MscL, during the measurements. To this end, asymmetric insertion of l-α-lysophosphatidylcholine (LPC) into the lipid bilayer has been effective; however, how LPC activates MscL is not fully understood. Here, the effects of LPC on tension-sensitive mutants of a bacterial MscL and on MscL homologs with different tension sensitivities are reported, leading to the conclusion that the mode of action of LPC is different from that of applied tension. Our results imply that LPC shifts the free energy of gating by interfering with MscL-membrane coupling. Furthermore, we demonstrate that the fine-tuned addition of LPC can be used for controlled activation of MscL in spectroscopic studies.

Utilizing zeta potential measurements to study the effective charge, membrane partitioning, and membrane permeation of the lipopeptide surfactin

The effective charge of membrane-active molecules such as the fungicidal lipopeptide surfactin (SF) is a crucial property governing solubility, membrane partitioning, and membrane permeability. We present zeta potential measurements of liposomes to measure the effective charge as well as membrane partitioning of SF by utilizing what we call an equi-activity analysis of several series of samples with different lipid concentrations. We observe an effective charge of -1.0 for SF at pH8.5 and insignificantly lower at pH7.4, illustrating that the effective charge may deviate strongly from the nominal value (-2 for 1 Asp, 1 Glu). The apparent partition coefficient decreases from roughly 100 to 20/mM with increasing membrane content of SF in agreement with the literature. Finally, by comparing zeta potentials measured soon after the addition of peptide to liposomes with those measured after a heat treatment to induce transmembrane equilibration of SF, we quantified the asymmetry of partitioning between the outer and inner leaflets. At very low concentration, SF binds exclusively to the outer leaflet. The onset of partial translocation to the inner leaflet occurs at about 5mol-% SF in the membrane. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Freezing of stressed bilayers and vesicles

We show using a minimalistic theoretical framework that phase transition decoupling in lipid bilayers is caused by a nonuniform stress profile due to an asymmetric distribution of lipids between the two leaflets. Applying this framework to vesicles, we demonstrate that their anomalous freezing is also caused by a stress asymmetry, but that this is due to lipid tail extension on freezing. Finally, we predict a previously unknown dependence of surface tension on temperature, find the phenomenon also in coarse grained molecular dynamics simulations, and suggest that it might have relevance in thermosensitive protein gating.

Deformation of a Two-domain Lipid Bilayer due to Asymmetric Insertion of Lipid-modified Ras Peptides

Ras proteins are attached to the inner leaflet of the plasma membrane via a lipid-modified anchor. Membrane-bound Ras proteins laterally segregate into nanoscale signaling platforms called nanoclusters. It has been shown that the membrane domain preference of Ras nanoclusters varies with the nature of lipidation but their effect on the membrane has not been well understood. To investigate the effect of Ras insertion on membrane structure, we carried out numerous coarse-grained molecular dynamics (CGMD) simulations on a two-domain DPPC/DLiPC/cholesterol lipid bilayer in which different number and type of H-ras peptides were attached on one side. We have shown previously that this lipid mixture forms co-existing liquid-ordered/liquid-disordered (L o /L d ) domains and that different H-ras peptides form clusters that variously accumulate at the L o /L d regions or the boundary between them. Here we show that asymmetric insertion of each of these peptides induces a vertical relative displacement of the domains and deforms the bilayer, with the domain boundary serving as the center of deformation. The extent of the deformation, however, varies with the type and number of lipid modification. This is because the number and type of the Ras lipid tails determines the degree to which the stress caused by asymmetric peptide insertion is relieved by inter-leaflet cholesterol transfer and lipid tilt. In addition, we have characterized the mechanism of bilayer deformation based on the collective effect of the Ras peptides on inter-leaflet surface area, pressure profile and line tension differences. This allowed us to elucidate how Ras lipid modification affects membrane geometry and how a two-domain bilayer adjusts its shape through boundary deformation. The result contributes to a better understanding of Ras signaling platforms and highlights some of the mechanisms by which a multi-domain membrane responds to external perturbation.

Embedding Aβ42 in heterogeneous membranes depends on cholesterol asymmetries

Using a coarse-grained lipid and peptide model, we show that the free energy stabilization of amyloid-β in heterogeneous lipid membranes is predicted to have a dependence on asymmetric distributions of cholesterol compositions across the membrane leaflets. We find that a highly asymmetric cholesterol distribution that is depleted on the exofacial leaflet but enhanced on the cytofacial leaflet of the model lipid membrane thermodynamically favors membrane retention of a fully embedded Aβ peptide. However, in the case of cholesterol redistribution that increases concentration of cholesterol on the exofacial layer, typical of aging or Alzheimer's disease, the free energy favors peptide extrusion of the highly reactive N-terminus into the extracellular space that may be vulnerable to aggregation, oligomerization, or deleterious oxidative reactivity.

Lipid nanotechnology

Nanotechnology is a multidisciplinary field that covers a vast and diverse array of devices and machines derived from engineering, physics, materials science, chemistry and biology. These devices have found applications in biomedical sciences, such as targeted drug delivery, bio-imaging, sensing and diagnosis of pathologies at early stages. In these applications, nano-devices typically interface with the plasma membrane of cells. On the other hand, naturally occurring nanostructures in biology have been a source of inspiration for new nanotechnological designs and hybrid nanostructures made of biological and non-biological, organic and inorganic building blocks. Lipids, with their amphiphilicity, diversity of head and tail chemistry, and antifouling properties that block nonspecific binding to lipid-coated surfaces, provide a powerful toolbox for nanotechnology. This review discusses the progress in the emerging field of lipid nanotechnology.