Free energy landscape of rim-pore expansion in membrane fusion

Transient pores in hemifusion diaphragms

Exchange of material across two membranes, as in the case of synaptic neurotransmitter release from a vesicle, involves the formation and poration of a hemifusion diaphragm (HD). The nontrivial geometry of the HD leads to environment-dependent control, regarding the stability and dynamics of the pores required for this kind of exocytosis. This work combines particle simulations, field-based calculations, and phenomenological modeling to explore the factors influencing the stability, dynamics, and possible control mechanisms of pores in HDs. We find that pores preferentially form at the HD rim, and that their stability is sensitive to a number of factors, including the three line tensions, membrane tension, HD size, and the ability of lipids to "flip-flop" across leaflets. Along with a detailed analysis of these factors, we discuss ways that vesicles or cells may use them to open and close pores and thereby quickly and efficiently transport material.

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.

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.

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.

Biomimetic Curvature and Tension-Driven Membrane Fusion Induced by Silica Nanoparticles

Fusion events in living cells are intricate phenomena that require the coordinate action of multicomponent protein complexes. However, simpler synthetic tools to control membrane fusion in artificial cells are highly desirable. Native membrane fusion machinery mediates fusion, driving a delicate balance of membrane curvature and tension between two closely apposed membranes. Here, we show that silica nanoparticles (SiO₂ NPs) at a size close to the cross-over between tension-driven and curvature-driven interaction regimes initiate efficient fusion of biomimetic model membranes. Fusion efficiency and mechanisms are studied by Förster resonance energy transfer and confocal fluorescence microscopy. SiO₂ NPs induce a slight increase in lipid packing likely to increase the lateral tension of the membrane. We observe a connection between membrane tension and fusion efficiency. Finally, real-time confocal fluorescence microscopy reveals three distinct mechanistic pathways for membrane fusion. SiO₂ NPs show significant potential for inclusion in the synthetic biology toolkit for membrane remodeling and fusion in artificial cells.

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.

Sequential Water and Headgroup Merger: Membrane Poration Paths and Energetics from MD Simulations

Membrane topology changes such as poration, stalk formation, and hemifusion rupture are essential to cellular function, but their molecular details, energetics, and kinetics are still not fully understood. Here, we present a unified energetic and mechanistic picture of metastable pore defects in tensionless lipid membranes. We used an exhaustive committor analysis to test and select optimal reaction coordinates and also to determine the nucleation mechanism. These reaction coordinates were used to calculate free-energy landscapes that capture the full process and end states. The identified barriers agree with the committor analysis. To enable sufficient sampling of the complete transition path for our molecular dynamics simulations, we developed a “gizmo” potential biasing scheme. The simulations suggest that the essential step in the nucleation is the initial merger of lipid headgroups at the nascent pore center. To facilitate this event, an indentation pathway is energetically preferred to a hydrophobic defect. Continuous water columns that span the indentation were determined to be on-path transients that precede the nucleation barrier. This study gives a quantitative description of the nucleation mechanism and energetics of small metastable pores and illustrates a systematic approach to uncover the mechanisms of diverse cellular membrane remodeling processes.

SNARE-Mediated Exocytosis in Neuronal Development

The formation of the nervous system involves establishing complex networks of synaptic connections between proper partners. This developmental undertaking requires the rapid expansion of the plasma membrane surface area as neurons grow and polarize, extending axons through the extracellular environment. Critical to the expansion of the plasma membrane and addition of plasma membrane material is exocytic vesicle fusion, a regulated mechanism driven by soluble N-ethylmaleimide-sensitive factor attachment proteins receptors (SNAREs). Since their discovery, SNAREs have been implicated in several critical neuronal functions involving exocytic fusion in addition to synaptic transmission, including neurite initiation and outgrowth, axon specification, axon extension, and synaptogenesis. Decades of research have uncovered a rich variety of SNARE expression and function. The basis of SNARE-mediated fusion, the opening of a fusion pore, remains an enigmatic event, despite an incredible amount of research, as fusion is not only heterogeneous but also spatially small and temporally fast. Multiple modes of exocytosis have been proposed, with full-vesicle fusion (FFV) and kiss-and-run (KNR) being the best described. Whereas most in vitro work has reconstituted fusion using VAMP-2, SNAP-25, and syntaxin-1; there is much to learn regarding the behaviors of distinct SNARE complexes. In the past few years, robust heterogeneity in the kinetics and fate of the fusion pore that varies by cell type have been uncovered, suggesting a paradigm shift in how the modes of exocytosis are viewed is warranted. Here, we explore both classic and recent work uncovering the variety of SNAREs and their importance in the development of neurons, as well as historical and newly proposed modes of exocytosis, their regulation, and proteins involved in the regulation of fusion kinetics.

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.

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.

Fusion Pores Live on the Edge

Biological transmission of vesicular content occurs by opening of a fusion pore. Recent experimental observations have illustrated that fusion pores between vesicles that are docked by an extended flat contact zone are located at the edge (vertex) of this zone. We modeled this experimentally observed scenario by coarse-grained molecular simulations and elastic theory. This revealed that fusion pores experience a direct attraction toward the vertex. The size adopted by the resulting vertex pore strongly depends on the apparent contact angle between the adhered vesicles even in the absence of membrane surface tension. Larger contact angles substantially increase the equilibrium size of the vertex pore. Because the cellular membrane fusion machinery actively docks membranes, it facilitates a collective expansion of the contact zone and increases the contact angle. In this way, the fusion machinery can drive expansion of the fusion pore by free energy equivalents of multiple tens of k_BT from a distance and not only through the fusion proteins that reside within the fusion pore.

Phosphatidic acid in membrane rearrangements

Phosphatidic acid (PA) is the simplest cellular glycerophospholipid characterized by unique biophysical properties: a small headgroup; negative charge; and a phosphomonoester group. Upon interaction with lysine or arginine, PA charge increases from -1 to -2 and this change stabilizes protein-lipid interactions. The biochemical properties of PA also allow interactions with lipids in several subcellular compartments. Based on this feature, PA is involved in the regulation and amplification of many cellular signalling pathways and functions, as well as in membrane rearrangements. Thereby, PA can influence membrane fusion and fission through four main mechanisms: it is a substrate for enzymes producing lipids (lysophosphatidic acid and diacylglycerol) that are involved in fission or fusion; it contributes to membrane rearrangements by generating negative membrane curvature; it interacts with proteins required for membrane fusion and fission; and it activates enzymes whose products are involved in membrane rearrangements. Here, we discuss the biophysical properties of PA in the context of the above four roles of PA in membrane fusion and fission.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

SNARE-mediated membrane fusion arrests at pore expansion to regulate the volume of an organelle

Constitutive membrane fusion within eukaryotic cells is thought to be controlled at its initial steps, membrane tethering and SNARE complex assembly, and to rapidly proceed from there to full fusion. Although theory predicts that fusion pore expansion faces a major energy barrier and might hence be a rate-limiting and regulated step, corresponding states with non-expanding pores are difficult to assay and have remained elusive. Here, we show that vacuoles in living yeast are connected by a metastable, non-expanding, nanoscopic fusion pore. This is their default state, from which full fusion is regulated. Molecular dynamics simulations suggest that SNAREs and the SM protein-containing HOPS complex stabilize this pore against re-closure. Expansion of the nanoscopic pore to full fusion can thus be triggered by osmotic pressure gradients, providing a simple mechanism to rapidly adapt organelle volume to increases in its content. Metastable, nanoscopic fusion pores are then not only a transient intermediate but can be a long-lived, physiologically relevant and regulated state of SNARE-dependent membrane fusion.

Rate of hemifusion diaphragm dissipation and ability to form three-junction bound HD determined by lipid composition

Though the hemifusion diaphragm (HD) is widely accepted as an intermediate in bilayer membrane fusion, lipid contributions toward HD stability and dynamics are still not fully understood. In this paper, we study large, binary, protein-free HD systems at varying compositions of negative intrinsic curvature (NIC) lipids using molecular dynamics (MD) simulations of a solvent-free coarse-grained lipid model. Under MD, initially created HDs are found to relax to three major end states depending on the composition and lipid intrinsic curvature. Low compositions of NIC lipids or weak intrinsic curvature result in double-bilayer end states, and moderate compositions of moderate to strong NIC lipids result in metastable fusion pores. Importantly, high compositions of moderate NIC lipids result in a metastable HD that persists beyond μs time scales. NIC lipids stabilize the HD by filling the junction core around the HD. Sorting of NIC lipids toward the three-junction region occurs in fused-endpoint systems, but no significant sorting was seen in systems that end in a double bilayer indicating that high line tension at the triple junction drives HD dissipation faster than sorting can enrich that junction enough to lower that line tension. The appearance of three end states dependent on the NIC lipid composition highlights the necessity of NIC lipids for non-leaky fusion.

Hemagglutinin-Mediated Membrane Fusion: A Biophysical Perspective

Influenza hemagglutinin (HA) is a viral membrane protein responsible for the initial steps of the entry of influenza virus into the host cell. It mediates binding of the virus particle to the host-cell membrane and catalyzes fusion of the viral membrane with that of the host. HA is therefore a major target in the development of antiviral strategies. The fusion of two membranes involves high activation barriers and proceeds through several intermediate states. Here, we provide a biophysical description of the membrane fusion process, relating its kinetic and thermodynamic properties to the large conformational changes taking place in HA and placing these in the context of multiple HA proteins working together to mediate fusion. Furthermore, we highlight the role of novel single-particle experiments and computational approaches in understanding the fusion process and their complementarity with other biophysical approaches.

Molecular dynamics simulations of membrane proteins and their interactions: from nanoscale to mesoscale

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics.

Protein-lipid interactions critical to replication of the influenza A virus

Influenza A virus (IAV) assembles on the plasma membrane where viral proteins localize to form a bud encompassing the viral genome, which ultimately pinches off to give rise to newly formed infectious virions. Upon entry, the virus faces the opposite task-fusion with the endosomal membrane and disassembly to deliver the viral genome to the cytoplasm. There are at least four influenza proteins-hemagglutinin (HA), neuraminidase (NA), matrix 1 protein (M1), and the M2 ion channel-that are known to directly interact with the cellular membrane and modify membrane curvature in order to both assemble and disassemble membrane-enveloped virions. Here, we summarize and discuss current knowledge of the interactions of lipids and membrane proteins involved in the IAV replication cycle.

Calculating Transition Energy Barriers and Characterizing Activation States for Steps of Fusion

We use continuum mechanics to calculate an entire least energy pathway of membrane fusion, from stalk formation, to pore creation, and through fusion pore enlargement. The model assumes that each structure in the pathway is axially symmetric. The static continuum stalk structure agrees quantitatively with experimental stalk architecture. Calculations show that in a stalk, the distal monolayer is stretched and the stored stretching energy is significantly less than the tilt energy of an unstretched distal monolayer. The string method is used to determine the energy of the transition barriers that separate intermediate states and the dynamics of two bilayers as they pass through them. Hemifusion requires a small amount of energy independently of lipid composition, while direct transition from a stalk to a fusion pore without a hemifusion intermediate is highly improbable. Hemifusion diaphragm expansion is spontaneous for distal monolayers containing at least two lipid components, given sufficiently negative diaphragm spontaneous curvature. Conversely, diaphragms formed from single-component distal monolayers do not expand without the continual injection of energy. We identify a diaphragm radius, below which central pore expansion is spontaneous. For larger diaphragms, prior studies have shown that pore expansion is not axisymmetric, and here our calculations supply an upper bound for the energy of the barrier against pore formation. The major energy-requiring deformations in the steps of fusion are: widening of a hydrophobic fissure in bilayers for stalk formation, splay within the expanding hemifusion diaphragm, and fissure widening initiating pore formation in a hemifusion diaphragm.

The Interaction between Influenza HA Fusion Peptide and Transmembrane Domain Affects Membrane Structure

Viral glycoproteins, such as influenza hemagglutinin (HA) and human immunodeficiency virus gp41, are anchored by a single helical segment transmembrane domain (TMD) on the viral envelope membrane. The fusion peptides (FP) of the glycoproteins insert into the host membrane and initiate membrane fusion. Our previous study showed that the FP or TMD alone perturbs membrane structure. Interaction between the influenza HA FP and TMD has previously been shown, but its role is unclear. We used PC spin labels dipalmitoylphospatidyl-tempo-choline (on the headgroup), 5PC and 14PC (5-C and 14-C positions on the acyl chain) to detect the combined effect of FP-TMD interaction by titrating HA FP to TMD-reconstituted 1,2-dimyristoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol)/cholesterol lipid bilayers using electron spin resonance. We found that the FP-TMD increases the lipid order at all positions, which has a greater lipid ordering effect than the sum of the FP or TMD alone, and this effect reaches deeper into the membranes. Although HA-mediated membrane fusion is pH dependent, this combined effect is observed at both pH 5 and pH 7. In addition to increasing lipid order, multiple components are found for 5PC at increased concentration of FP-TMD, indicating that distinct domains are induced. However, the mutation of Gly1 in the FP and L187 in the TMD eliminates the perturbations, consistent with their fusogenic phenotypes. Electron spin resonance on spin-labeled peptides confirms these observations. We suggest that this interaction may provide a driving force in different stages of membrane fusion: initialization, transition from hemifusion stalk to transmembrane contact, and fusion pore formation.