Thermodynamically reversible paths of the first fusion intermediate reveal an important role for membrane anchors of fusion proteins

Cell-inspired, massive electromodulation of friction via transmembrane fields across lipid bilayers

Transient electric fields across cell bilayer membranes can lead to electroporation and cell fusion, effects crucial to cell viability whose biological implications have been extensively studied. However, little is known about these behaviours in a materials context. Here we find that transmembrane electric fields can lead to a massive, reversible modulation of the sliding friction between surfaces coated with lipid-bilayer membranes-a 200-fold variation, up to two orders of magnitude greater than that achieved to date. Atomistic simulations reveal that the transverse fields, resembling those at cell membranes, lead to fully reversible electroporation of the confined bilayers and the formation of inter-bilayer bridges analogous to the stalks preceding intermembrane fusion. These increase the interfacial dissipation through reduced hydration at the slip plane, forcing it to revert in part from the low-dissipation, hydrated lipid-headgroup plane to the intra-bilayer, high-dissipation acyl tail interface. Our results demonstrate that lipid bilayers under transmembrane electric fields can have striking materials modification properties.

Novel non-helical antimicrobial peptides insert into and fuse lipid model membranes

This research addresses the growing menace of antibiotic resistance by exploring antimicrobial peptides (AMPs) as alternatives to conventional antibiotics. Specifically, we investigate two linear amphipathic AMPs, LE-53 (12-mer) and LE-55 (16-mer), finding that the shorter LE-53 exhibits greater bactericidal activity against both Gram-negative (G(-)) and Gram-positive (G(+)) bacteria. Remarkably, both AMPs are non-toxic to eukaryotic cells. The heightened effectiveness of LE-53 is attributed to its increased hydrophobicity (H) compared to LE-55. Circular dichroism (CD) reveals that LE-53 and LE-55 both adopt β-sheet and random coil structures in lipid model membranes (LMMs) mimicking G(-) and G(+) bacteria, so secondary structure is not the cause of the potency difference. X-ray diffuse scattering (XDS) reveals increased lipid chain order in LE-53, a potential key distinction. Additionally, XDS study uncovers a significant link between LE-53's upper hydrocarbon location in G(-) and G(+) LMMs and its efficacy. Neutron reflectometry (NR) confirms the AMP locations determined using XDS. Solution small angle X-ray scattering (SAXS) demonstrates LE-53's ability to induce vesicle fusion in bacterial LMMs without affecting eukaryotic LMMs, offering a promising strategy to combat antibiotic-resistant strains while preserving human cell integrity, whereas LE-55 has a smaller ability to induce fusion.

Myomerger Induces Membrane Hemifusion and Regulates Fusion Pore Expansion

Membrane fusion is a crucial mechanism in a wide variety of important events in cell biology from viral infection to exocytosis. However, despite many efforts and much progress, cell-cell fusion has remained elusive to our understanding. Along the life of the fusion pore, large conformational changes take place from the initial lipid bilayer bending, passing through the hemifusion intermediates, and ending with the formation of the first nascent fusion pore. In this sense, computer simulations are an ideal technique for describing such complex lipid remodeling at the molecular level. In this work, we studied the role played by the muscle-specific membrane protein Myomerger during the formation of the fusion pore. We have conducted μs length atomistic and coarse-grained molecular dynamics, together with free-energy calculations using ad hoc collective variables. Our results show that Myomerger favors the hemifusion diaphragm-stalk transition, reduces the nucleation-expansion energy difference, and promotes the formation of nonenlarging fusion pores.

Mechanisms of SNARE proteins in membrane fusion

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) are a family of small conserved eukaryotic proteins that mediate membrane fusion between organelles and with the plasma membrane. SNAREs are directly or indirectly anchored to membranes. Prior to fusion, complementary SNAREs assemble between membranes with the aid of accessory proteins that provide a scaffold to initiate SNARE zippering, pulling the membranes together and mediating fusion. Recent advances have enabled the construction of detailed models describing bilayer transitions and energy barriers along the fusion pathway and have elucidated the structures of SNAREs complexed in various states with regulatory proteins. In this Review, we discuss how these advances are yielding an increasingly detailed picture of the SNARE-mediated fusion pathway, leading from first contact between the membranes via metastable non-bilayer intermediates towards the opening and expansion of a fusion pore. We describe how SNARE proteins assemble into complexes, how this assembly is regulated by accessory proteins and how SNARE complexes overcome the free energy barriers that prevent spontaneous membrane fusion.

The local variation of the Gaussian modulus enables different pathways for fluid lipid vesicle fusion

Viral infections, fertilization, neurotransmission, and many other fundamental biological processes rely on membrane fusion. Straightforward calculations based on the celebrated Canham-Helfrich elastic model predict a large topological energy barrier that prevents the fusion process from being thermally activated. While such high energy is in accordance with the physical barrier function of lipid membranes, it is difficult to reconcile with the biological mechanisms involved in fusion processes. In this work, we use a Ginzburg-Landau type of free energy that recovers the Canham-Helfrich model in the limit of small width-to-vesicle-extension ratio, with the additional ability to handle topological transitions. We show that a local modification of the Gaussian modulus in the merging region both dramatically lowers the elastic energy barrier and substantially changes the minimal energy pathway for fusion, in accordance with experimental evidence. Therefore, we discuss biological examples in which such a modification might play a crucial role.

Free Energy Surface and Molecular Mechanism of Slow Structural Transitions in Lipid Bilayers

Lipid membrane remodeling, crucial for many cellular processes, is governed by the coupling of membrane structure and shape fluctuations. Given the importance of the ∼ nm length scale, details of the transition intermediates for conformational change are not fully captured by a continuum-mechanical description. Slow dynamics and the lack of knowledge of reaction coordinates (RCs) for biasing methods pose a challenge for all-atom (AA) simulations. Here, we map system dynamics on Langevin dynamics in a normal mode space determined from an elastic network model representation for the lipid-water Hamiltonian. AA molecular dynamics (MD) simulations are used to determine model parameters, and Langevin dynamics predictions for bilayer structural, mechanical, and dynamic properties are validated against MD simulations and experiments. Transferability to describe the dynamics of a larger lipid bilayer and a heterogeneous membrane-protein system is assessed. A set of generic RCs for pore formation in two tensionless bilayers is obtained by coupling Langevin dynamics to the underlying energy landscape for membrane deformations. Structure evolution is carried out by AA MD, wherein the generic RCs are used in a path metadynamics or an umbrella sampling simulation to determine the thermodynamics of pore formation and its molecular determinants, such as the role of distinct bilayer motions, lipid solvation, and lipid packing.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Membrane Tension Inhibits Lipid Mixing by Increasing the Hemifusion Stalk Energy

Fusion of biological membranes is fundamental in various physiological events. The fusion process involves several intermediate stages with energy barriers that are tightly dependent on the mechanical and physical properties of the system, one of which is membrane tension. As previously established, the late stages of fusion, including hemifusion diaphragm and pore expansions, are favored by membrane tension. However, a current understanding of how the energy barrier of earlier fusion stages is affected by membrane tension is lacking. Here, we apply a newly developed experimental approach combining micropipette-aspirated giant unilamellar vesicles and optically trapped membrane-coated beads, revealing that membrane tension inhibits lipid mixing. We show that lipid mixing is 6 times slower under a tension of 0.12 mN/m compared with tension-free membranes. Furthermore, using continuum elastic theory, we calculate the dependence of the hemifusion stalk formation energy on membrane tension and intermembrane distance and find the increase in the corresponding energy barrier to be 1.6 kBT in our setting, which can explain the increase in lipid mixing time delay. Finally, we show that tension can be a significant factor in the stalk energy if the pre-fusion intermembrane distance is on the order of several nanometers, while for membranes that are tightly docked, tension has a negligible effect.

Computational modeling of membrane trafficking processes: From large molecular assemblies to chemical specificity

In the last decade, molecular dynamics (MD) simulations have become an essential tool to investigate the molecular properties of membrane trafficking processes, often in conjunction with experimental approaches. The combination of MD simulations with recent developments in structural biology, such as cryo-electron microscopy and artificial intelligence-based structure determination, opens new, exciting possibilities for future investigations. However, the full potential of MD simulations to provide a molecular view of the complex and dynamic processes involving membrane trafficking can only be realized if certain limitations are addressed, and especially those concerning the quality of coarse-grain models, which, despite recent successes in describing large-scale systems, still suffer from far-from-ideal chemical accuracy. In this review, we will highlight recent success stories of MD simulations in the investigation of membrane trafficking processes, their implications for future research, and the challenges that lie ahead in this specific research domain.

Cholesterol-Containing Liposomes Decorated With Au Nanoparticles as Minimal Tunable Fusion Machinery

Membrane fusion is essential for the basal functionality of eukaryotic cells. In physiological conditions, fusion events are regulated by a wide range of specialized proteins, operating with finely tuned local lipid composition and ionic environment. Fusogenic proteins, assisted by membrane cholesterol and calcium ions, provide the mechanical energy necessary to achieve vesicle fusion in neuromediator release. Similar cooperative effects must be explored when considering synthetic approaches for controlled membrane fusion. We show that liposomes decorated with amphiphilic Au nanoparticles (AuLips) can act as minimal tunable fusion machinery. AuLips fusion is triggered by divalent ions, while the number of fusion events dramatically changes with, and can be finely tuned by, the liposome cholesterol content. We combine quartz-crystal-microbalance with dissipation monitoring (QCM-D), fluorescence assays, and small-angle X-ray scattering (SAXS) with molecular dynamics (MD) at coarse-grained (CG) resolution, revealing new mechanistic details on the fusogenic activity of amphiphilic Au nanoparticles (AuNPs) and demonstrating the ability of these synthetic nanomaterials to induce fusion regardless of the divalent ion used (Ca²⁺ or Mg²⁺ ). The results provide a novel contribution to developing new artificial fusogenic agents for next-generation biomedical applications that require tight control of the rate of fusion events (e.g., targeted drug delivery).

Two modes of fusogenic action for influenza virus fusion peptide

The entry of influenza virus into the host cell requires fusion of its lipid envelope with the host membrane. It is catalysed by viral hemagglutinin protein, whose fragments called fusion peptides become inserted into the target bilayer and initiate its merging with the viral membrane. Isolated fusion peptides are already capable of inducing lipid mixing between liposomes. Years of studies indicate that upon membrane binding they form bend helical structure whose degree of opening fluctuates between tightly closed hairpin and an extended boomerang. The actual way in which they initiate fusion remains elusive. In this work we employ atomistic simulations of wild type and fusion inactive W14A mutant of influenza fusion peptides confined between two closely apposed lipid bilayers. We characterise peptide induced membrane perturbation and determine the potential of mean force for the formation of the first fusion intermediate, an interbilayer lipid bridge called stalk. Our results demonstrate two routes through which the peptides can lower free energy barrier towards fusion. The first one assumes peptides capability to adopt transmembrane configuration which subsequently promotes the creation of a stalk-hole complex. The second involves surface bound peptide configuration and proceeds owing to its ability to stabilise stalk by fitting into the region of extreme negative membrane curvature resulting from its formation. In both cases, the active peptide conformation corresponds to tight helical hairpin, whereas extended boomerang geometry appears to be unable to provide favourable thermodynamic effect. The latter observation offers plausible explanation for long known inactivity of boomerang-stabilising W14A mutation.

Intrinsic Disorder in α-Synuclein Regulates the Exocytotic Fusion Pore Transition

Today, it is widely accepted that intrinsic disorder is strongly related to the cell cycle, during mitosis, differentiation, and apoptosis. Of particular interest are hybrid proteins possessing both structured and unstructured domains that are critical in human health and disease, such as α-synuclein. In this work, we describe how α-synuclein interacts with the nascent fusion pore as it evolves toward expansion. We unveil the key role played by its intrinsically disordered region as a thermodynamic regulator of the nucleation-expansion energy barrier. By analyzing a truncated variant of α-synuclein that lacks the disordered region, we find that the landscape of protein interactions with PIP₂ and POPS lipids is highly altered, ultimately increasing the energy cost for the fusion pore to transit from nucleation to expansion. We conclude that the intrinsically disordered region in full-length α-synuclein recognizes and allocates pivotal protein:lipid interactions during membrane remodeling in the first stages of the fusion pore.

High curvature promotes fusion of lipid membranes: Predictions from continuum elastic theory

The fusion of lipid membranes progresses through a series of hemifusion intermediates with two significant energy barriers related to the formation of stalk and fusion pore, respectively. These energy barriers determine the speed and success rate of many critical biological processes, including the fusion of highly curved membranes, for example synaptic vesicles and enveloped viruses. Here we use continuum elastic theory of lipid monolayers to determine the relationship between membrane shape and energy barriers to fusion. We find that the stalk formation energy decreases with curvature by up to 31 kBT in a 20-nm-radius vesicle compared with planar membranes and by up to 8 kBT in the fusion of highly curved, long, tubular membranes. In contrast, the fusion pore formation energy barrier shows a more complicated behavior. Immediately after stalk expansion to the hemifusion diaphragm, the fusion pore formation energy barrier is low (15-25 kBT) due to lipid stretching in the distal monolayers and increased tension in highly curved vesicles. Therefore, the opening of the fusion pore is faster. However, these stresses relax over time due to lipid flip-flop from the proximal monolayer, resulting in a larger hemifusion diaphragm and a higher fusion pore formation energy barrier, up to 35 kBT. Therefore, if the fusion pore fails to open before significant lipid flip-flop takes place, the reaction proceeds to an extended hemifusion diaphragm state, which is a dead-end configuration in the fusion process and can be used to prevent viral infections. In contrast, in the fusion of long tubular compartments, the surface tension does not accumulate due to the formation of the diaphragm, and the energy barrier for pore expansion increases with curvature by up to 11 kBT. This suggests that inhibition of polymorphic virus infection could particularly target this feature of the second barrier.

Membrane permeabilization can be crucially biased by a fusogenic lipid composition - leaky fusion caused by antimicrobial peptides in model membranes

Induced membrane permeabilization or leakage is often taken as an indication for activity of membrane-active molecules, such as antimicrobial peptides (AMPs). The exact leakage mechanism is often unknown, but important, because certain mechanisms might actually contribute to microbial killing, while others are unselective, or potentially irrelevant in an in vivo situation. Using an antimicrobial example peptide (cR₃W₃), we illustrate one of the potentially misleading leakage mechanisms: leaky fusion, where leakage is coupled to membrane fusion. Like many others, we examine peptide-induced leakage in model vesicles consisting of binary mixtures of anionic and zwitterionic phospholipids. In fact, phosphatidylglycerol and phosphatidylethanolamine (PG/PE) are supposed to reflect bacterial membranes, but exhibit a high propensity for vesicle aggregation and fusion. We describe the implications of this vesicle fusion and aggregation for the reliability of model studies. The ambiguous role of the relatively fusogenic PE-lipids becomes clear as leakage decreases significantly when aggregation and fusion are prevented by sterical shielding. Furthermore, the mechanism of leakage changes if PE is exchanged for phosphatidylcholine (PC). We thus point out that the lipid composition of model membranes can be biased towards leaky fusion. This can lead to discrepancies between model studies and activity in true microbes, because leaky fusion is likely prevented by bacterial peptidoglycan layers. In conclusion, choosing the model membrane might implicate the type of effect (here leakage mechanism) that is observed. In the worst case, as with leaky fusion of PG/PE vesicles, this is not directly relevant for the intended antimicrobial application.

On the role of membrane embedding, protein rigidity and transmembrane length in lipid membrane fusion

The fusion of biological membranes is ubiquitous in natural processes like exo- and endocytosis, intracellular trafficking and viral entry. Membrane fusion is also utilized in artificial biomimetic fusion systems, e.g. for drug delivery. Both the natural and the biomimetic fusion systems rely on a wide range of (artificial) proteins mediating the fusion process. Although the exact mechanisms of these proteins differ, clear analogies in their general behavior can be observed in bringing the membranes in close proximity and mediating the fusion reaction. In our study, we use molecular dynamics simulations with coarse grained models, mimicking the general behavior of fusion proteins (spikes), to systematically examine the effects of specific characteristics of these proteins on the fusion process. The protein characteristics considered are (i) the type of membrane embedding, i.e., either transmembrane or not, (ii) the rigidity, and (iii) the transmembrane domain (TMD) length. The results show essential differences in fusion pathway between monotopic and transmembrane spikes, in which transmembrane spikes seem to inhibit the formation of hemifusion diaphragms, leading to a faster fusion development. Furthermore, we observed that an increased rigidity and a decreased TMD length both proved to contribute to a faster fusion development. Finally, we show that a single spike may suffice to successfully induce a fusion reaction, provided that the spike is sufficiently rigid and attractive. Not only does this shed light on biological fusion of membranes, it also provides clear design rules for artificial membrane fusion systems.

Tuning the Extracellular Vesicles Membrane through Fusion for Biomedical Applications

Membrane fusion is one of the key phenomena in the living cell for maintaining the basic function of life. Extracellular vesicles (EVs) have the ability to transfer information between cells through plasma membrane fusion, making them a promising tool in diagnostics and therapeutics. This study explores the potential applications of natural membrane vesicles, EVs, and their fusion with liposomes, EVs, and cells and introduces methodologies for enhancing the fusion process. EVs have a high loading capacity, bio-compatibility, and stability, making them ideal for producing effective drugs and diagnostics. The unique properties of fused EVs and the crucial design and development procedures that are necessary to realize their potential as drug carriers and diagnostic tools are also examined. The promise of EVs in various stages of disease management highlights their potential role in future healthcare.

Computational Advances of Protein/Neurotransmitter-membrane Interactions Involved in Vesicle Fusion and Neurotransmitter Release

Vesicle fusion is of crucial importance to neuronal communication at neuron terminals. The exquisite but complex fusion machinery for neurotransmitter release is tightly controlled and regulated by protein/neurotransmitter-membrane interactions. Computational 'microscopies', in particular molecular dynamics simulations and related techniques, have provided notable insight into the physiological process over the past decades, and have made enormous contributions to fields such as neurology, pharmacology and pathophysiology. Here we review the computational advances of protein/neurotransmitter-membrane interactions related to presynaptic vesicle-membrane fusion and neurotransmitter release, and outline the in silico challenges ahead for understanding this important physiological process.

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.

Activation energy and force fields during topological transitions of fluid lipid vesicles

Topological transitions of fluid lipid membranes are fundamental processes for cell life. For example, they are required for endo- and exocytosis or to enable neurotransmitters to cross the neural synapses. Here, inspired by the idea that fusion and fission proteins could have evolved in Nature in order to carry out a minimal work expenditure, we evaluate the minimal free energy pathway for the transition between two spherical large unilamellar vesicles and a dumbbell-shaped one. To address the problem, we propose and successfully use a Ginzburg-Landau type of free energy, which allows us to uniquely describe without interruption the whole, full-scale topological change. We also compute the force fields needed to overcome the involved energy barriers. The obtained forces are in excellent agreement, in terms of intensity, scale, and spatial localization with experimental data on typical fission protein systems, whereas they suggest the presence of additional features in fusion proteins.

A Core-Shell Approach for Systematically Coarsening Nanoparticle-Membrane Interactions: Application to Silver Nanoparticles

The continuous release of engineered nanomaterial (ENM) into the environment may bring about health concerns following human exposure. One important source of ENMs are silver nanoparticles (NPs) that are extensively used as anti-bacterial additives. The introduction of ENMs into the human body can occur via ingestion, skin uptake or the respiratory system. Therefore, evaluating how NPs translocate over bio-membranes is essential in assessing their primary toxicity. Unfortunately, data regarding membrane-NP interaction is still scarce, as is theoretical and in silico insight into what governs adhesion and translocation for the most relevant NPs and membranes. Coarse-grained (CG) molecular descriptions have the potential to alleviate this situation, but are hampered by the absence of a direct link to NP materials and membrane adhesion mechanisms. Here, we interrogate the relationship between the most common NP representation at the CG level and the adhesion characteristics of a model lung membrane. We find that this representation for silver NPs is non-transferable, meaning that a proper CG representation for one size is not suited for other sizes. We also identify two basic types of primary adhesion-(partial) NPs wrapping by the membrane and NP insertion into the membrane-that closely relate to the overall NP hydrophobicity and significantly differ in terms of lipid coatings. The proven non-transferability of the standard CG representation with size forms an inspiration for introducing a core-shell model even for bare NPs that are uniform in composition. Using existing all-atom molecular dynamics (MD) data as a reference, we show that this extension does allow us to reproduce size-dependent NP adhesion properties and lipid responses to NP binding at the CG level. The subsequent CGMD evaluation for 10 nm Ag NPs provides new insight into membrane binding for relevant NP sizes and into the role of water in trapping NPs into defected mixed monolayer-bilayer states. This development will be instrumental for simulating NP-membrane adhesion towards more experimentally relevant length and time scales for particular NP materials.

Enhanced Expansion and Reduced Kiss-and-Run Events in Fusion Pores Steered by Synaptotagmin-1 C2B Domains

The fusion pore controls the release of exocytotic vesicle contents through a precise orchestration of lipids from the fusing membranes and proteins. There is a major lipid reorganization during the different stages in life of the fusion pore (membrane fusion, nucleation, and expansion) that can be scrutinized thermodynamically. In this work, using umbrella sampling simulations we describe the expansion of the fusion pore. We have calculated free energy profiles to drive a nascent, just nucleated, fusion pore to its expanded configuration. We have quantified the effects on the free energy of one and two Synaptotagmin-1 C2B domains in the cytosolic space. We show that C2B domains cumulatively reduce the cost for expansion, favoring the system to evolve toward full fusion. Finally, by conducting thousands of unbiased molecular dynamics simulations, we show that C2B domains significantly decrease the probability of kiss-and-run events.

Modeling biomolecular kinetics with large-scale simulation

The molecular details of biomolecular kinetics present a challenging estimation problem because the identities of relevant intermediates and the rates of exchange between them must be determined. These can be derived from prior knowledge, but in recent years, great advances have been made in the development and application of methods to systematically determine states and rates using biomolecular simulation. Doing this for biological systems of reasonable complexity requires substantial computational power, and contemporary methods leverage distributed computing or leadership-class computing resources to accomplish this. The result has been substantial insight into pressing contemporary problems, including structural activation of pandemic viruses. Here, we highlight recent developments in both methodology and exciting applications.

Cell-Penetrating Peptides

In this introductory chapter, we first define cell-penetrating peptides (CPPs), give short overview of CPP history and discuss several aspects of CPP classification. Next section is devoted to the mechanism of CPP penetration into the cells, where direct and endocytic internalization of CPP is explained. Kinetics of internalization is discussed more extensively, since this topic is not discussed in other chapters of this book. At the end of this section some features of the thermodynamics of CPP interaction with the membrane is also presented. Finally, we present different cargoes that can be transferred into the cells by CPPs and briefly discuss the effect of cargo on the rate and efficiency of penetration into the cells.

Synaptotagmin-1 C2B domains cooperatively stabilize the fusion stalk via a master-servant mechanism

Synaptotagmin-1 is a low-affinity Ca2+ sensor that triggers synchronous vesicle fusion. It contains two similar C2 domains (C2A and C2B) that cooperate in membrane binding, being the C2B domain the...

Free energies of membrane stalk formation from a lipidomics perspective

Many biological membranes are asymmetric and exhibit complex lipid composition, comprising hundreds of distinct chemical species. Identifying the biological function and advantage of this complexity is a central goal of membrane biology. Here, we study how membrane complexity controls the energetics of the first steps of membrane fusions, that is, the formation of a stalk. We first present a computationally efficient method for simulating thermodynamically reversible pathways of stalk formation at coarse-grained resolution. The method reveals that the inner leaflet of a typical plasma membrane is far more fusogenic than the outer leaflet, which is likely an adaptation to evolutionary pressure. To rationalize these findings by the distinct lipid compositions, we computed ~200 free energies of stalk formation in membranes with different lipid head groups, tail lengths, tail unsaturations, and sterol content. In summary, the simulations reveal a drastic influence of the lipid composition on stalk formation and a comprehensive fusogenicity map of many biologically relevant lipid classes.

Fusion of cellular membranes begins with the formation of a stalk. Here, the authors develop a computationally efficient method for coarse-grained simulations of stalk formation and apply this approach to comprehensively analyse how stalk formation is influenced by the membrane lipid composition.

Cooperation of Conical and Polyunsaturated Lipids to Regulate Initiation and Processing of Membrane Fusion

The shape of lipids has long been suspected to be a critical determinant for the control of membrane fusion. To experimentally test this assertion, we used conical and malleable lipids and measured their influence on the fusion kinetics. We found that, as previously suspected, both types of lipids accelerate fusion. However, the implicated molecular mechanisms are strikingly different. Malleable lipids, with their ability to change shape with low energy cost, favor fusion by decreasing the overall activation energy. On the other hand, conical lipids, with their small polar head relative to the area occupied by the hydrophobic chains, tend to make fusion less energetically advantageous because they tend to migrate towards the most favorable lipid leaflet, hindering fusion pore opening. They could however facilitate fusion by generating hydrophobic defects on the membranes; this is suggested by the similar trend observed between the experimental rate of fusion nucleation and the surface occupied by hydrophobic defects obtained by molecular simulations. The synergy of dual-process, activation energy and nucleation kinetics, could facilitate membrane fusion regulation in vivo.

Artificial stabilization of the fusion pore by intra-organelle styrene-maleic acid copolymers

Styrene-maleic acid copolymers have become an advantageous detergent-free alternative for membrane protein isolation. Since their discovery, experimental membrane protein extraction and purification by keeping intact their lipid environment has become significantly easier. With the aim of identifying new applications of these interesting copolymers, their molecular binding and functioning mechanisms have recently become intense objects of study. In this work, we describe the use of styrene-maleic acid copolymers as an artificial tool to stabilize the fusion pore. We show that when these copolymers circumscribe the water channel that defines the fusion pore, they keep it from shrinking and closing. We describe how only intra-organelle copolymers have stabilizing capabilities while extra-organelle ones have negligible or even contrary effects on the fusion pore life-time.

Fast bilayer-micelle fusion mediated by hydrophobic dipeptides

To understand the transition from inanimate matter to life, we studied a process that directly couples simple metabolism to evolution via natural selection, demonstrated experimentally by Adamala and Szostak. In this process, dipeptides synthesized inside precursors of cells promote absorption of fatty acid micelles to vesicles, inducing their preferential growth and division at the expense of other vesicles. The process is explained on the basis of coarse-grained molecular dynamics simulations, each extending for tens of microseconds, carried out to model fusion between a micelle and a membrane, both made of fatty acids in the absence and presence of hydrophobic dipeptides. In all systems with dipeptides, but not in their absence, fusion events were observed. They involve the formation of a stalk made by hydrophobic chains from the micelle and the membrane, similar to that postulated for vesicle-vesicle fusion. The emergence of a stalk is facilitated by transient clusters of dipeptides, side chains of which form hydrophobic patches at the membrane surface. Committor probability calculations indicate that the size of a patch is a suitable reaction coordinate and allows for identifying the transition state for fusion. Free-energy barrier to fusion is greatly reduced in the presence of dipeptides to only 4-5 kcal/mol, depending on the hydrophobicity of side chains. The mechanism of mediated fusion, which is expected to apply to other small peptides and hydrophobic molecules, provides a robust means by which a nascent metabolism can confer evolutionary advantage to precursors of cells.

Joint Reaction Coordinate for Computing the Free-Energy Landscape of Pore Nucleation and Pore Expansion in Lipid Membranes

Topological transitions of membranes, such as pore formation or membrane fusion, play key roles in biology, biotechnology, and in medical applications. Calculating the related free-energy landscapes has been complicated by the fact that such processes involve a sequence of transitions along highly distinct directions in conformational space, making it difficult to define good reaction coordinates (RCs) for the overall process. In this study, a new RC capable of driving both pore nucleation and pore expansion in lipid membranes is presented. The potential of mean force (PMF) along the RC computed with molecular dynamics simulations provides a comprehensive view on the free-energy landscape of pore formation, including a barrier for pore nucleation; the size, free energy, and metastability of the open pore; and the energetic cost for further pore expansion against the line tension of the pore rim. The RC is illustrated by quantifying the effects of (i) simulation system size and (ii) the addition of dimethyl sulfoxide on the free-energy landscape of pore formation. PMF calculations along the RC provide mechanistic and energetic understanding of pore formation, hence they will be useful to rationalize the effects of membrane-active peptides, electric fields, and membrane composition on transmembrane pores.

How does curvature affect the free-energy barrier of stalk formation? Small vesicles vs apposing, planar membranes

Using molecular simulations of POPC lipids in conjunction with the calculation of the Minimum Free-Energy Path (MFEP), we study the effect of strong membrane curvature on the formation of the first fusion intermediate—the stalk between a vesicle and its periodic image. We find that the thermodynamic stability of this hourglass-shaped, hydrophobic connection between two vesicles is largely increased by the strong curvature of small vesicles, whereas the intrinsic barrier to form a stalk, i.e., associated with dimple formation and lipid tails protrusions, is similar to the case of two, apposing, planar membranes. A significant reduction of the barrier of stalk formation, however, stems from the lower dehydration free energy that is required to bring highly curved vesicle into a distance, at which stalk formation may occur, compared to the case of apposing, planar membranes.

Computational Insight Into the Mechanism of SARS-CoV-2 Membrane Fusion

Membrane fusion, a key step in the early stages of virus propagation, allows the release of the viral genome in the host cell cytoplasm. The process is initiated by fusion peptides that are small, hydrophobic components of viral membrane-embedded glycoproteins and are typically conserved within virus families. Here, we attempted to identify the correct fusion peptide region in the Spike protein of SARS-CoV-2 by all-atom molecular dynamics simulations of dual membrane systems with varied oligomeric units of putative candidate peptides. Of all of the systems tested, only a trimeric unit of a 40-amino-acid region (residues 816-855 of SARS-CoV-2 Spike) was effective in triggering the initial stages of membrane fusion, within 200 ns of simulation time. Association of this trimeric unit with dual membranes resulted in the migration of lipids from the upper leaflet of the lower bilayer toward the lower leaflet of the upper bilayer to create a structural unit reminiscent of a fusion bridge. We submit that residues 816-855 of Spike represent the bona fide fusion peptide of SARS-CoV-2 and that computational methods represent an effective way to identify fusion peptides in viral glycoproteins.

How proteins open fusion pores: insights from molecular simulations

Fusion proteins can play a versatile and involved role during all stages of the fusion reaction. Their roles go far beyond forcing the opposing membranes into close proximity to drive stalk formation and fusion. Molecular simulations have played a central role in providing a molecular understanding of how fusion proteins actively overcome the free energy barriers of the fusion reaction up to the expansion of the fusion pore. Unexpectedly, molecular simulations have revealed a preference of the biological fusion reaction to proceed through asymmetric pathways resulting in the formation of, e.g., a stalk-hole complex, rim-pore, or vertex pore. Force-field based molecular simulations are now able to directly resolve the minimum free-energy path in protein-mediated fusion as well as quantifying the free energies of formed reaction intermediates. Ongoing developments in Graphics Processing Units (GPUs), free energy calculations, and coarse-grained force-fields will soon gain additional insights into the diverse roles of fusion proteins.

Where are those lipid nano rings?

Highly curved toroidal micelles with diameters as small as 100 nm have been successfully constructed by self-assembly of amphiphilic block copolymers. These structures may have potential applications in gene or drug delivery. Experimental observations suggest that toroidal micelles likely originate from spherical or disc-like micelles which are tricked into forming toroidal micelles upon external stimuli ('smart' materials). Since self-assembly of polymeric and lipid surfactants is guided by the same physical principles, we hypothesize that 'smart' lipid surfactants can be equivalently tricked into forming highly curved toroidal micelles that are tenfold smaller (≃10 nm diameter). Paradoxically, these 'nano rings' have never been observed. Using coarse-grained molecular dynamics (MD) simulations in conjunction with a state-of-the-art free energy calculation method (a string method), we illustrate how a thermo-responsive lipid surfactant is able to form toroidal micelles. These micelles originate from disc-like micelles that are spontaneously perforated upon heat shocking, thereby supporting a longstanding hypothesis on the possible origin of polymeric toroidal micelle phases observed in experiments. We illustrate that kinetically stable 'nano rings' are substantially shorter lived than their tenfold larger polymeric analogs. The estimated life-time (milliseconds) is in fact similar to the characteristic breaking time of the corresponding worm-like micelle. Finally, we resolve the characteristic finger print which 'nano rings' leave in time-resolved X-ray spectra and illustrate how the uptake of small DNA fragments may enhance their stability. Despite a shared kinetics of self-assembly, length scale dependent differences in the life-time of surfactant phases can occur when phases are kinetically rather than thermodynamically stable. This results in the apparent absence or presence of toroidal micelle phases on different length scales. Our theoretical work precisely illustrates that the universality of surfactants nevertheless remains conserved even at different length scales.

Effect of pH on the influenza fusion peptide properties unveiled by constant-pH molecular dynamics simulations combined with experiment

The influenza virus fusion process, whereby the virus fuses its envelope with the host endosome membrane to release the genetic material, takes place in the acidic late endosome environment. Acidification triggers a large conformational change in the fusion protein, hemagglutinin (HA), which enables the insertion of the N-terminal region of the HA2 subunit, known as the fusion peptide, into the membrane of the host endosome. However, the mechanism by which pH modulates the molecular properties of the fusion peptide remains unclear. To answer this question, we performed the first constant-pH molecular dynamics simulations of the influenza fusion peptide in a membrane, extending for 40 µs of aggregated time. The simulations were combined with spectroscopic data, which showed that the peptide is twofold more active in promoting lipid mixing of model membranes at pH 5 than at pH 7.4. The realistic treatment of protonation introduced by the constant-pH molecular dynamics simulations revealed that low pH stabilizes a vertical membrane-spanning conformation and leads to more frequent contacts between the fusion peptide and the lipid headgroups, which may explain the increase in activity. The study also revealed that the N-terminal region is determinant for the peptide’s effect on the membrane.

Cell-cell interfaces as specialized compartments directing cell function

Cell-cell interfaces are found throughout multicellular organisms, from transient interactions between motile immune cells to long-lived cell-cell contacts in epithelia. Studies of immune cell interactions, epithelial cell barriers, neuronal contacts and sites of cell-cell fusion have identified a core set of features shared by cell-cell interfaces that critically control their function. Data from diverse cell types also show that cells actively and passively regulate the localization, strength, duration and cytoskeletal coupling of receptor interactions governing cell-cell signalling and physical connections between cells, indicating that cell-cell interfaces have a unique membrane organization that emerges from local molecular and cellular mechanics. In this Review, we discuss recent findings that support the emerging view of cell-cell interfaces as specialized compartments that biophysically constrain the arrangement and activity of their protein, lipid and glycan components. We also review how these biophysical features of cell-cell interfaces allow cells to respond with high selectivity and sensitivity to multiple inputs, serving as the basis for wide-ranging cellular functions. Finally, we consider how the unique properties of cell-cell interfaces present opportunities for therapeutic intervention.

Quantitative Synaptic Biology: A Perspective on Techniques, Numbers and Expectations

Synapses play a central role for the processing of information in the brain and have been analyzed in countless biochemical, electrophysiological, imaging, and computational studies. The functionality and plasticity of synapses are nevertheless still difficult to predict, and conflicting hypotheses have been proposed for many synaptic processes. In this review, we argue that the cause of these problems is a lack of understanding of the spatiotemporal dynamics of key synaptic components. Fortunately, a number of emerging imaging approaches, going beyond super-resolution, should be able to provide required protein positions in space at different points in time. Mathematical models can then integrate the resulting information to allow the prediction of the spatiotemporal dynamics. We argue that these models, to deal with the complexity of synaptic processes, need to be designed in a sufficiently abstract way. Taken together, we suggest that a well-designed combination of imaging and modelling approaches will result in a far more complete understanding of synaptic function than currently possible.

Density Field Thermodynamic Integration (DFTI): A "Soft" Approach to Calculate the Free Energy of Surfactant Self-Assemblies

Thermodynamic integration is one of the most established methods to quantify excess free energies between different metastable states. Excess intermolecular interactions in surfactant assemblies are on the scale of the energy of thermal fluctuations. Therefore, these materials can be deformed and topologically altered via relatively small mechanical stresses. It is thus intuitive to design reaction paths and associated order parameters that exploit the "soft" nature of these materials to mechanically rather than alchemically morph surfactant assemblies from state to state. Here, we propose a novel method coined "density field thermodynamic integration" (DFTI) that adopts the universality and transferability of alchemical methods while simultaneously exploiting the soft excess interactions between surfactant molecules. DFTI was designed for a rapid quantification of the free energy differences between different metastable structures in soft fluid materials. The DFTI method uses an external field coupled to the local density to mechanically morph the system between metastable states of interest. Here, we explored the capability of the DFTI method to swiftly and accurately calculate free energy differences between states. To this aim, we studied two different coarse-grained lipidic surfactant systems: (i) a fusion stalk and (ii) a worm-like micelle. Our results illustrate that DFTI can provide an efficient, versatile, and rather reliable method to calculate the free energy differences between surfactant assemblies.

SNAREs, tethers and SM proteins: how to overcome the final barriers to membrane fusion?

Physiological membrane vesicles are built to separate reaction spaces in a stable manner, even when they accidentally collide or are kept in apposition by spatial constraints in the cell. This requires a natural resistance to fusion and mixing of their content, which originates from substantial energetic barriers to membrane fusion [1]. To facilitate intracellular membrane fusion reactions in a controlled manner, proteinaceous fusion machineries have evolved. An important open question is whether protein fusion machineries actively pull the fusion reaction over the present free energy barriers, or whether they rather catalyze fusion by lowering those barriers. At first sight, fusion proteins such as SNARE complexes and viral fusion proteins appear to act as nano-machines, which mechanically transduce force to the membranes and thereby overcome the free energy barriers [2,3]. Whether fusion proteins additionally alter the free energy landscape of the fusion reaction via catalytic roles is less obvious. This is a question that we shall discuss in this review, with particular focus on the influence of the eukaryotic SNARE-dependent fusion machinery on the final step of the reaction, the formation and expansion of the fusion pore.

Synergistic actions of v-SNARE transmembrane domains and membrane-curvature modifying lipids in neurotransmitter release

Vesicle fusion is mediated by assembly of SNARE proteins between opposing membranes. While previous work suggested an active role of SNARE transmembrane domains (TMDs) in promoting membrane merger (Dhara et al., 2016), the underlying mechanism remained elusive. Here, we show that naturally-occurring v-SNARE TMD variants differentially regulate fusion pore dynamics in mouse chromaffin cells, indicating TMD flexibility as a mechanistic determinant that facilitates transmitter release from differentially-sized vesicles. Membrane curvature-promoting phospholipids like lysophosphatidylcholine or oleic acid profoundly alter pore expansion and fully rescue the decelerated fusion kinetics of TMD-rigidifying VAMP2 mutants. Thus, v-SNARE TMDs and phospholipids cooperate in supporting membrane curvature at the fusion pore neck. Oppositely, slowing of pore kinetics by the SNARE-regulator complexin-2 withstands the curvature-driven speeding of fusion, indicating that pore evolution is tightly coupled to progressive SNARE complex formation. Collectively, TMD-mediated support of membrane curvature and SNARE force-generated membrane bending promote fusion pore formation and expansion.

A molecular view on the escape of lipoplexed DNA from the endosome

The use of non-viral vectors for in vivo gene therapy could drastically increase safety, whilst reducing the cost of preparing the vectors. A promising approach to non-viral vectors makes use of DNA/cationic liposome complexes (lipoplexes) to deliver the genetic material. Here we use coarse-grained molecular dynamics simulations to investigate the molecular mechanism underlying efficient DNA transfer from lipoplexes. Our computational fusion experiments of lipoplexes with endosomal membrane models show two distinct modes of transfection: parallel and perpendicular. In the parallel fusion pathway, DNA aligns with the membrane surface, showing very quick release of genetic material shortly after the initial fusion pore is formed. The perpendicular pathway also leads to transfection, but release is slower. We further show that the composition and size of the lipoplex, as well as the lipid composition of the endosomal membrane, have a significant impact on fusion efficiency in our models.

Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif

Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.

Influenza hemagglutinin drives viral entry via two sequential intramembrane mechanisms

Enveloped viruses enter cells via a process of membrane fusion between the viral envelope and a cellular membrane. For influenza virus, mutational data have shown that the membrane-inserted portions of the hemagglutinin protein play a critical role in achieving fusion. In contrast to the relatively well-understood ectodomain, a predictive mechanistic understanding of the intramembrane mechanisms by which influenza hemagglutinin drives fusion has been elusive. We used molecular dynamics simulations of fusion between a full-length hemagglutinin proteoliposome and a lipid bilayer to analyze these mechanisms. In our simulations, hemagglutinin first acts within the membrane to increase lipid tail protrusion and promote stalk formation and then acts to engage the distal leaflets of each membrane and promote stalk widening, curvature, and eventual fusion. These two sequential mechanisms, one occurring before stalk formation and one after, are consistent with our experimental measurements of single-virus fusion kinetics to liposomes of different sizes. The resulting model also helps explain and integrate previous mutational and biophysical data, particularly the mutational sensitivity of the fusion peptide N terminus and the length sensitivity of the transmembrane domain. We hypothesize that entry by other enveloped viruses may also use sequential processes of acyl tail exposure, followed by membrane curvature and distal leaflet engagement.

Pathways connecting two opposed bilayers with a fusion pore: a molecularly-informed phase field approach

A phase field model with two phase fields, representing the concentration and the head-tail separation of amphiphilic molecules, respectively, has been constructed using an extension of the Ohta-Kawasaki model (Macromolecules, 1986, 19, 2621-2632). It is shown that this molecularly-informed phase field model is capable of producing various self-assembled amphiphilic aggregates, such as bilayers, vesicles and micelles. Furthermore, pathways connecting two opposed bilayers with a fusion pore are obtained by using a combination of the phase field model and the string method. Multiple fusion pathways, including a classical pathway and a leaky pathway, have been obtained depending on the initial separation of the two bilayers. The study shed light on the understanding of the membrane fusion pathways and, more importantly, laid a foundation for further investigation of more complex membrane morphologies and transitions.

Structural Basis for Broad HIV-1 Neutralization by the MPER-Specific Human Broadly Neutralizing Antibody LN01

Potent and broadly neutralizing antibodies (bnAbs) are the hallmark of HIV-1 protection by vaccination. The membrane-proximal external region (MPER) of the HIV-1 gp41 fusion protein is targeted by the most broadly reactive HIV-1 neutralizing antibodies. Here, we examine the structural and molecular mechansims of neutralization by anti-MPER bnAb, LN01, which was isolated from lymph-node-derived germinal center B cells of an elite controller and exhibits broad neutralization breadth. LN01 engages both MPER and the transmembrane (TM) region, which together form a continuous helix in complex with LN01. The tilted TM orientation allows LN01 to interact simultaneously with the peptidic component of the MPER epitope and membrane via two specific lipid binding sites of the antibody paratope. Although LN01 carries a high load of somatic mutations, most key residues interacting with the MPER epitope and lipids are germline encoded, lending support for the LN01 epitope as a candidate for lineage-based vaccine development.

•

bNAb LN01 neutralizes 92% of a 118-strain virus panel

•

LN01 targets the HIV-1 gp41 MPER, the TM region, and lipids

•

LN01-complexed MPER forms a continuous helix with TM

•

Most LN01 paratope residues interacting with MPER-TM and lipids are germline encoded

Direct Derivation of Free Energies of Membrane Deformation and Other Solvent Density Variations From Enhanced Sampling Molecular Dynamics

We report a methodology to calculate the free energy of a shape transformation in a lipid membrane directly from a molecular dynamics simulation. The bilayer need not be homogeneous or symmetric and can be atomically detailed or coarse grained. The method is based on a collective variable that quantifies the similarity between the membrane and a set of predefined density distributions. Enhanced sampling of this "Multi-Map" variable re-shapes the bilayer and permits the derivation of the corresponding potential of mean force. Calculated energies thus reflect the dynamic interplay of atoms and molecules, rather than postulated effects. Evaluation of deformations of different shape, amplitude, and range demonstrates that the macroscopic bending modulus assumed by the Helfrich-Canham model is increasingly unsuitable below the 100-Å scale. In this range of major biological significance, direct free-energy calculations reveal a much greater plasticity. We also quantify the stiffening effect of cholesterol on bilayers of different composition and compare with experiments. Lastly, we illustrate how this approach facilitates analysis of other solvent reorganization processes, such as hydrophobic hydration. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Curvature induction and sensing of the F-BAR protein Pacsin1 on lipid membranes via molecular dynamics simulations

F-Bin/Amphiphysin/Rvs (F-BAR) domain proteins play essential roles in biological processes that involve membrane remodelling, such as endocytosis and exocytosis. It has been shown that such proteins transform the lipid membrane into tubes. Notably, Pacsin1 from the Pacsin/Syndapin subfamily has the ability to transform the membrane into various morphologies: striated tubes, featureless wide and thin tubes, and pearling vesicles. The molecular mechanism of this interesting ability remains elusive. In this study, we performed all-atom (AA) and coarse-grained (CG) molecular dynamics simulations to investigate the curvature induction and sensing mechanisms of Pacsin1 on a membrane. From AA simulations, we show that Pacsin1 has internal structural flexibility. In CG simulations with parameters tuned from the AA simulations, spontaneous assembly of two Pacsin1 dimers through lateral interaction is observed. Based on the complex structure, we show that the regularly assembled Pacsin1 dimers bend a tensionless membrane. We also show that a single Pacsin1 dimer senses the membrane curvature, binding to a buckled membrane with a preferred curvature. These results provide molecular insights into polymorphic membrane remodelling.

Peripheral Antimicrobial Peptide Gomesin Induces Membrane Protrusion, Folding, and Laceration

Optical microscopy shows that the peripheral antimicrobial peptide (AMP) gomesin does not disrupt the bacterial membrane by forming stable transmembrane pores but induces lipid accumulation domains, which is followed by a sudden burst near the domains. The molecular action mechanisms of gomesin on vesicle and planar bilayer membranes are investigated in this work using coarse-grained molecular dynamics simulations. By comparing the membrane morphology and property changes induced by gomesin and the pore-forming AMP melittin, we determined that the amphiphilic shape of the AMPs is a key factor affecting the mechanism of cell death. The binding of wedge-shaped gomesin, with a small hydrophobic surface, onto the membrane induces protrusion and folding of the outer monolayer followed by sudden membrane lacerations at the axillae of the protuberances. Alternatively, cylinder-shaped melittins with comparable hydrophilic and hydrophobic surfaces destroy membranes by forming stable pores coexisting with exocytosis-like buddings and endocytosis-like invaginations. The multiple actions of AMPs on the bacterial membrane suggest diverse paradigms for designing molecular carriers for delivering drugs to the cell.