The Gaussian curvature elastic energy of intermediates in membrane fusion

Comparing Multifunctional Viral and Eukaryotic Proteins for Generating Scission Necks in Membranes

Deterministic formation of membrane scission necks by protein machinery with multiplexed functions is critical in biology. A microbial example is M2 viroporin, a proton pump from the influenza A virus that is multiplexed with membrane remodeling activity to induce budding and scission in the host membrane during viral maturation. In comparison, the dynamin family constitutes a class of eukaryotic proteins implicated in mitochondrial fission, as well as various budding and endocytosis pathways. In the case of Dnm1, the mitochondrial fission protein in yeast, the membrane remodeling activity is multiplexed with mechanoenzyme activity to create fission necks. It is not clear why these functions are combined in these scission processes, which occur in drastically different compositions and solution conditions. In general, direct experimental access to changing neck sizes induced by individual proteins or peptide fragments is challenging due to the nanoscale dimensions and influence of thermal fluctuations. Here, we use a mechanical model to estimate the size of scission necks by leveraging small-angle X-ray scattering structural data of protein-lipid systems under different conditions. The influence of interfacial tension, lipid composition, and membrane budding morphology on the size of the induced scission necks is systematically investigated using our data and molecular dynamic simulations. We find that the M2 budding protein from the influenza A virus has robust pH-dependent membrane activity that induces nanoscopic necks within the range of spontaneous hemifission for a broad range of lipid compositions. In contrast, the sizes of scission necks generated by mitochondrial fission proteins strongly depend on lipid composition, which suggests a role for mechanical constriction.

A Coronin 1-Derived Peptide Inhibits Membrane Fusion by Modulating Membrane Organization and Dynamics

Membrane fusion is considered the first step in the entry of enveloped viruses into the host cell. Several targeted strategies have been implemented to block viral entry by limiting the fusion protein to form a six-helix bundle, which is a prerequisite for fusion. Nonetheless, the development of broad-spectrum fusion inhibitors is essential to combat emerging and re-emerging viral infections. TG-23, a coronin 1, a tryptophan-aspartate-rich phagosomal protein-derived peptide, demonstrated inhibition of fusion between small unilamellar vesicles (SUVs) by modulating the membrane's physical properties. However, its inhibitory efficacy reduces with an increasing concentration of membrane cholesterol. The present work aims to develop a fusion inhibitor whose efficacy would be unaltered in the presence of membrane cholesterol. A stretch of the tryptophan-aspartic acid-containing peptide with a similar secondary structure and hydrophobicity profile of TG-23 from coronin 1 was synthesized, and its ability to inhibit SUV-SUV fusion with varying concentrations of membrane cholesterol was evaluated. Our results demonstrate that the GG-21 peptide inhibits fusion irrespective of the cholesterol content of the membrane. We have further evaluated the peptide-induced change in the membrane organization and dynamics utilizing arrays of steady-state and time-resolved fluorescence measurements and correlated these results with their effect on fusion. Interestingly, GG-21 displays inhibitory efficacy in a wide variety of lipid compositions despite having a secondary structure and physical properties similar to those of TG-23. Overall, our results advocate that the secondary structure and physical properties of the peptide may not be sufficient to predict its inhibitory efficacy.

Synaptotagmin 7 C2 domains induce membrane curvature stress via electrostatic interactions and the wedge mechanism

Synaptotagmin 7 (Syt-7) is part of the synaptotagmin protein family that regulates exocytotic lipid membrane fusion. Among the family, Syt-7 stands out by its membrane binding strength and stabilization of long-lived membrane fusion pores. Given that Syt-7 vesicles form long-lived fusion pores, we hypothesize that its interactions with the membrane stabilize the specific curvatures, thicknesses, and lipid compositions that support a metastable fusion pore. Using all-atom molecular dynamics simulations and FRET-based assays of Syt-7's membrane-binding C2 domains (C2A and C2B), we found that Syt-7 C2 domains sequester anionic lipids, are sensitive to cholesterol, thin membranes, and generate lipid membrane curvature by two competing, but related mechanisms. First, Syt-7 forms strong electrostatic contacts with the membrane, generating negative curvature stress. Second, Syt-7's calcium binding loops embed in the membrane surface, acting as a wedge to thin the membrane and induce positive curvature stress. These curvature mechanisms are linked by the protein insertion depth as well as the resulting protein tilt. Simplified quantitative models of the curvature-generating mechanisms link simulation observables to their membrane-reshaping effectiveness.

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.

Lipid nanoparticle topology regulates endosomal escape and delivery of RNA to the cytoplasm

RNA therapeutics have the potential to resolve a myriad of genetic diseases. Lipid nanoparticles (LNPs) are among the most successful RNA delivery systems. Expanding their use for the treatment of more genetic diseases hinges on our ability to continuously evolve the design of LNPs with high potency, cellular-specific targeting, and low side effects. Overcoming the difficulty of releasing cargo from endocytosed LNPs remains a significant hurdle. Here, we investigate the fundamental properties of nonviral RNA nanoparticles pertaining to the activation of topological transformations of endosomal membranes and RNA translocation into the cytosol. We show that, beyond composition, LNP fusogenicity can be prescribed by designing LNP nanostructures that lower the energetic cost of fusion and fusion-pore formation with a target membrane. The inclusion of structurally active lipids leads to enhanced LNP endosomal fusion, fast evasion of endosomal entrapment, and efficacious RNA delivery. For example, conserving the lipid make-up, RNA-LNPs having cuboplex nanostructures are significantly more efficacious at endosomal escape than traditional lipoplex constructs.

Simulated dynamic cholesterol redistribution favors membrane fusion pore constriction

Endo- and exocytosis proceed through a highly strained membrane fusion pore topology regardless of the aiding protein machinery. The membrane's lipid components bias fusion pores toward expansion or closure, modifying the necessary work done by proteins. Cholesterol, a key component of plasma membranes, promotes both inverted lipid phases with concave leaflets (i.e., negative total curvature, which thins the leaflet) and flat bilayer phases with thick, ordered hydrophobic interiors. We demonstrate by theory and simulation that both leaflets of nascent catenoidal fusion pores have negative total curvature. Furthermore, the hydrophobic core of bilayers with strong negative Gaussian curvature is thinned. Therefore, it is an open question whether cholesterol will be enriched in these regions because of the negative total curvature or depleted because of the membrane thinning. Here, we compare all-atom molecular dynamics simulations (built using a procedure to create specific fusion pore geometries) and theory to understand the underlying reasons for lipid redistribution on fusion pores. Our all-atom molecular dynamics simulations resolve this question by showing that cholesterol is strongly excluded from the thinned neck of fusion and fission pores, revealing that thickness (and/or lipid order) influences cholesterol distributions more than curvature. The results imply that cholesterol exclusion can drive fusion pore closure by creating a small, cholesterol-depleted zone in the neck. This model agrees with literature evidence that membrane reshaping is connected to cholesterol-dependent lateral phase separation.

Vesicle formation induced by thermal fluctuations

The process of fission and vesicle formation depends on the geometry of the membrane that will split. For instance, a flat surface finds it difficult to form vesicles because of the lack of curved regions where to start the process. Here we show that vesicle formation can be promoted by temperature, by using a membrane phase field model with Gaussian curvature. We find a phase transition between fluctuating and vesiculation phases that depends on temperature, spontaneous curvature, and the ratio between bending and Gaussian moduli. We analysed the energy dynamical behaviour of these processes and found that the main driving ingredient is the Gaussian energy term, although the curvature energy term usually helps with the process as well. We also found that the chemical potential can be used to investigate the temperature of the system. Finally we address how temperature changes the condition for spontaneous vesiculation for all geometries, making it happen in a wider range of values of the Gaussian modulus.

Serinc5 Restricts HIV Membrane Fusion by Altering Lipid Order and Heterogeneity in the Viral Membrane

The host restriction factor, Serinc5, incorporates into budding HIV particles and inhibits their infection by an incompletely understood mechanism. We have previously reported that Serinc5 but not its paralogue, Serinc2, blocks HIV cell entry by membrane fusion, specifically by inhibiting fusion pore formation and dilation. A body of work suggests that Serinc5 may alter the conformation and clustering of the HIV fusion protein, Env. To contribute an additional perspective to the developing model of Serinc5 restriction, we assessed Serinc2 and Serinc5's effects on HIV pseudoviral membranes. By measuring pseudoviral membrane thickness via cryo-electron microscopy and order via the fluorescent dye, FLIPPER-TR, Serinc5 was found to increase membrane heterogeneity, skewing the distribution toward a larger fraction of the viral membrane in an ordered phase. We also directly observed for the first time the coexistence of membrane domains within individual viral membrane envelopes. Using a total internal reflection fluorescence-based single particle fusion assay, we found that treatment of HIV pseudoviral particles with phosphatidylethanolamine (PE) rescued HIV pseudovirus fusion from restriction by Serinc5, which was accompanied by decreased membrane heterogeneity and order. This effect was specific for PE and did not depend on acyl chain length or saturation. Together, these data suggest that Serinc5 alters multiple interrelated properties of the viral membrane─lipid chain order, rigidity, line tension, and lateral pressure─which decrease the accessibility of fusion intermediates and disfavor completion of fusion. These biophysical insights into Serinc5 restriction of HIV infectivity could contribute to the development of novel antivirals that exploit the same weaknesses.

Light-Responsive Supramolecular Nanotubes-Based Chiral Plasmonic Assemblies

We describe the fabrication of dual-responsive (thermo/light) chiral plasmonic films. The idea is based on using photoswitchable achiral liquid crystal (LCs) forming chiral nanotubes for templating helical assemblies of Au NPs. Circular dichroism spectroscopy (CD) confirms chiroptical properties coming from the arrangement of organic and inorganic components, with up to 0.2 dissymmetry factor (g-factor). Upon exposure to UV light, organic molecules isomerize, resulting in controlled melting of organic nanotubes and/or inorganic nanohelices. The process can be reversed using visible light and further modified by varying the temperature, offering a control of chiroptical response of the composite material. These properties can play a key role in the future development of chiral plasmonics, metamaterials, and optoelectronic devices.

Effect of monolayer spontaneous curvature on constant tension-induced pore formation in lipid bilayers

The methodology of constant tension-induced rupture of giant unilamellar vesicles (GUVs) has provided information on tension-induced pore formation. This method was used to investigate the effect of spontaneous curvature (H₀) for a lipid monolayer on the rate constant (k_r) for constant tension (σ)-induced rupture, which originates from pore formation in lipid bilayers. Lipids were incorporated with different H₀ values into GUV membranes to change the overall H₀ value for the GUV monolayer. The dioleoylphosphatidylglycerol (DOPG)/dioleoylphosphatidylethanolamine (DOPE) (4/6, molar ratio, here and elsewhere) monolayer has a negative H₀, whereas the DOPG/dioleoylphosphatidylcholine (DOPC) (4/6) monolayer has an essentially zero H₀. A higher tension was required to induce the rupture of DOPG/DOPE (4/6)-GUVs compared with DOPG/DOPC (4/6)-GUVs. The line tension (Γ) for a pre-pore in DOPG/DOPE (4/6)-GUVs, determined by the analysis of the tension dependence of k_r, was 1.5 times larger than that in DOPG/DOPC (4/6)-GUVs. The k_r values for GUVs comprising DOPG/DOPC/18:1 lysophosphatidylcholine (LPC) (40/55/10), which has a positive H₀, were larger than those for DOPG/DOPC (4/6)-GUVs under the same tension. The Γ value for DOPG/DOPC/LPC (40/55/10)-GUVs was almost half that for DOPG/DOPC (4/6)-GUVs. These results indicate that Γ decreases with increasing H₀, which results in an increase in k_r. Based on these results, the effect of H₀ on k_r and Γ is discussed.

Determination of elastic parameters of lipid membranes from simulation under varied external pressure

Many cellular processes such as endocytosis, exocytosis, and vesicle trafficking involve membrane deformations, which can be analyzed in the framework of the elastic theories of lipid membranes. These models operate with phenomenological elastic parameters. A connection between these parameters and the internal structure of lipid membranes can be provided by three-dimensional (3D) elastic theories. Considering a membrane as a 3D layer, Campelo et al. [F. Campelo et al., Adv. Colloid Interface Sci. 208, 25 (2014)10.1016/j.cis.2014.01.018] developed a theoretical basis for the calculation of elastic parameters. In this work we generalize and improve this approach by considering a more general condition of global incompressibility instead of local incompressibility. Crucially, we find an important correction to the theory of Campelo et al., which if not taken into account leads to a significant miscalculation of elastic parameters. With the total volume conservation taken into account, we derive an expression for the local Poisson's ratio, which determines how the local volume changes upon stretching and permits a more precise determination of elastic parameters. Also, we substantially simplify the procedure by calculating the derivatives of the moments of the local tension with respect to stretching instead of calculating the local stretching modulus. We obtain a relation between the Gaussian curvature modulus as a function of stretching and the bending modulus, showing that these two elastic parameters are not independent, as was previously assumed. The proposed algorithm is applied to membranes composed of pure dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and their mixture. The following elastic parameters of these systems are obtained: the monolayer bending and stretching moduli, spontaneous curvature, neutral surface position, and local Poisson's ratio. It is shown that the bending modulus of the DPPC/DOPC mixture follows a more complex trend than predicted by the classical Reuss averaging, which is often employed in theoretical frameworks.

Interferon-induced transmembrane protein 3 (IFITM3) and its antiviral activity

Infections caused by enveloped viruses require fusion with cellular membranes for viral genome entry. Viral entry occurs following an interaction of viral and cellular membranes allowing the formation of fusion pores, by which the virus accesses the cytoplasm. Here, we focus on interferon-induced transmembrane protein 3 (IFITM3) and its antiviral activity. IFITM3 is predicted to block or stall viral fusion at an intermediate state, causing viral propagation to fail. After introducing IFITM3, we describe the generalized lipid membrane fusion pathway and how it can be stalled, particularly with respect to IFITM3, and current questions regarding IFITM3's topology, with specific emphasis on IFITM3's amphipathic α-helix (AAH) ₅₉V-₆₈M, which is necessary for the antiviral activity. We report new hydrophobicity and hydrophobic moment calculations for this peptide and a variety of active site peptides from known membrane-remodeling proteins. Finally, we discuss the effects of posttranslational modifications and localization, how IFITM3's AAH may block viral fusion, and possible ramifications of membrane composition.

Determination of Elastic Parameters of Lipid Membranes with Molecular Dynamics: A Review of Approaches and Theoretical Aspects

Lipid membranes are abundant in living organisms, where they constitute a surrounding shell for cells and their organelles. There are many circumstances in which the deformations of lipid membranes are involved in living cells: fusion and fission, membrane-mediated interaction between membrane inclusions, lipid-protein interaction, formation of pores, etc. In all of these cases, elastic parameters of lipid membranes are important for the description of membrane deformations, as these parameters determine energy barriers and characteristic times of membrane-involved phenomena. Since the development of molecular dynamics (MD), a variety of in silico methods have been proposed for the determination of elastic parameters of simulated lipid membranes. These MD methods allow for the consideration of details unattainable in experimental techniques and represent a distinct scientific field, which is rapidly developing. This work provides a review of these MD approaches with a focus on theoretical aspects. Two main challenges are identified: (i) the ambiguity in the transition from the continuum description of elastic theories to the discrete representation of MD simulations, and (ii) the determination of intrinsic elastic parameters of lipid mixtures, which is complicated due to the composition-curvature coupling effect.

The interaction of viral fusion peptides with lipid membranes

In this paper, we studied fusogenic peptides of class I-III fusion proteins, which are relevant to membrane fusion for certain enveloped viruses, in contact with model lipid membranes. We resolved the vertical structure and examined the adsorption or penetration behavior of the fusogenic peptides at phospholipid Langmuir monolayers with different initial surface pressures with x-ray reflectometry. We show that the fusion loops of tick-borne encephalitis virus (TBEV) glycoprotein E and vesicular stomatitis virus (VSV) G-protein are not able to insert deeply into model lipid membranes, as they adsorbed mainly underneath the headgroups with only limited penetration depths into the lipid films. In contrast, we observed that the hemagglutinin 2 fusion peptide (HA2-FP) and the VSV-transmembrane domain (VSV-TMD) can penetrate deeply into the membranes. However, in the case of VSV-TMD, the penetration was suppressed already at low surface pressures, whereas HA2-FP was able to insert even into highly compressed films. Membrane fusion is accompanied by drastic changes of the membrane curvature. To investigate how the peptides affect the curvature of model lipid membranes, we examined the effect of the fusogenic peptides on the equilibration of cubic monoolein structures after a phase transition from a lamellar state induced by an abrupt hydrostatic pressure reduction. We monitored this process in presence and absence of the peptides with small-angle x-ray scattering and found that HA2-FP and VSV-TMD drastically accelerate the equilibration, while the fusion loops of TBEV and VSV stabilize the swollen state of the lipid structures. In this work, we show that the class I fusion peptide of HA2 penetrates deeply into the hydrophobic region of membranes and is able to promote and accelerate the formation of negative curvature. In contrast, we found that the class II and III fusion loops of TBEV and VSV tend to counteract negative membrane curvature.

Aspects of Biological Replication and Evolution Independent of the Central Dogma: Insights from Protein-Free Vesicular Transformations and Protein-Mediated Membrane Remodeling

Biological membrane remodeling is central to living systems. In spite of serving as “containers” of whole-living systems and functioning as dynamic compartments within living systems, biological membranes still find a “blue collar” treatment compared to the “white collar” nucleic acids and proteins in biology. This may be attributable to the fact that scientific literature on biological membrane remodeling is only 50 years old compared to ~ 150 years of literature on proteins and a little less than 100 years on nucleic acids. However, recently, evidence for symbiotic origins of eukaryotic cells from data only on biological membranes was reported. This, coupled with appreciation of reproducible amphiphilic self-assemblies in aqueous environments (mimicking replication), has already initiated discussions on origins of life beyond nucleic acids and proteins. This work presents a comprehensive compilation and meta-analyses of data on self-assembly and vesicular transformations in biological membranes—starting from model membranes to establishment of Influenza Hemagglutinin-mediated membrane fusion as a prototypical remodeling system to a thorough comparison between enveloped mammalian viruses and cellular vesicles. We show that viral membrane fusion proteins, in addition to obeying “stoichiometry-driven protein folding”, have tighter compositional constraints on their amino acid occurrences than general-structured proteins, regardless of type/class. From the perspective of vesicular assemblies and biological membrane remodeling (with and without proteins) we find that cellular vesicles are quite different from viruses. Finally, we propose that in addition to pre-existing thermodynamic frameworks, kinetic considerations in de novo formation of metastable membrane structures with available “third-party” constituents (including proteins) were not only crucial for origins of life but also continue to offer morphological replication and/or functional mechanisms in modern life forms, independent of the central dogma.

Bicontinuous inverted cubic phase stabilization as an index of antimicrobial and membrane fusion peptide activity

Some antimicrobial peptides (AMPs) and membrane fusion-catalyzing peptides (FPs) stabilize bicontinuous inverted cubic (QII) phases. Previous authors proposed a topological rationale: since AMP-induced pores, fusion intermediates, and QII phases all have negative Gaussian curvature (NGC), peptides which produce NGC in one structure also do it in another. This assumes that peptides change the curvature energy of the lipid membranes. Here I test this with a Helfrich curvature energy model. First, experimentally, I show that lipid systems often used to study peptide NGC have NGC without peptides at higher temperatures. To determine the net effect of an AMP on NGC, the equilibrium phase behavior of the host lipids must be determined. Second, the model shows that AMPs must make large changes in the curvature energy to stabilize AMP-induced pores. Peptide-induced changes in elastic constants affect pores and QII phase differently. Changes in spontaneous curvature affect them in opposite ways. The observed correlation between QII phase stabilization and AMP activity doesn't show that AMPs act by lowering pore curvature energy. A different rationale is proposed. In theory, AMPs could simultaneously stabilize QII phase and pores by drastically changing two particular elastic constants. This could be tested by measuring AMP effects on the individual constants. I propose experiments to do that. Unlike AMPs, FPs must make only small changes in the curvature energy to catalyze fusion. It they act in this way, their fusion activity should correlate with their ability to stabilize QII phases.

Self-consistent field modeling of mesomorphic phase changes of monoolein and phospholipids in response to additives

Mapping the topological phase behaviour of lipids in aqueous solution is time consuming and finding the ideal lipid system for a desired application is often a matter of trial and error. Modelling techniques that can accurately predict the mesomorphic phase behaviour of lipid systems are therefore of paramount importance. Here, the self-consistent field theory of Scheutjens and Fleer (SF-SCF) in which a lattice refinement has been implemented, is used to scrutinize how various additives modify the self-assembled phase behaviour of monoolein (MO) and 1,2-dioleoyl-phosphatidylcholine (DOPC) lipids in water. The mesomorphic behaviour is inferred from trends in the mechanical properties of equilibrium lipid bilayers with increasing additive content. More specifically, we focus on the Helfrich parameters, that is, the mean and Gaussian bending rigidities (κ and [small kappa, Greek, macron], respectively) supplemented with the spontaneous curvature of the monolayer (Jm0). We use previously established interaction parameters that position the unperturbed DOPC system in the lamellar Lα phase ([small kappa, Greek, macron] < 0, κ > 0 and Jm0 ≈ 0). Similar interaction parameters position the MO system firmly in a bicontinuous cubic phase ([small kappa, Greek, macron] > 0). In line with experimental data, a mixture of MO and DOPC tends to be in one of these two phases, depending on the mixing ratio. Moreover we find good correlations between predicted trends and experimental data concerning the phase changes of MO in response to a wide range of additives. These correlations give credibility to the use of SF-SCF modelling as a valuable tool to quickly explore the mesomorphic phase space of (phospho)lipid bilayer systems including additives.

Bridging the Antimicrobial Activity of Two Lactoferricin Derivatives in E. coli and Lipid-Only Membranes

We coupled the antimicrobial activity of two well-studied lactoferricin derivatives, LF11-215 and LF11-324, in Escherichia coli and different lipid-only mimics of its cytoplasmic membrane using a common thermodynamic framework for peptide partitioning. In particular, we combined an improved analysis of microdilution assays with ζ-potential measurements, which allowed us to discriminate between the maximum number of surface-adsorbed peptides and peptides fully partitioned into the bacteria. At the same time, we measured the partitioning of the peptides into vesicles composed of phosphatidylethanolamine (PE), phosphatidylgylcerol (PG), and cardiolipin (CL) mixtures using tryptophan fluorescence and determined their membrane activity using a dye leakage assay and small-angle X-ray scattering. We found that the vast majority of LF11-215 and LF11-324 readily enter inner bacterial compartments, whereas only 1-5% remain surface bound. We observed comparable membrane binding of both peptides in membrane mimics containing PE and different molar ratios of PG and CL. The peptides' activity caused a concentration-dependent dye leakage in all studied membrane mimics; however, it also led to the formation of large aggregates, part of which contained collapsed multibilayers with sandwiched peptides in the interstitial space between membranes. This effect was least pronounced in pure PG vesicles, requiring also the highest peptide concentration to induce membrane permeabilization. In PE-containing systems, we additionally observed an effective shielding of the fluorescent dyes from leakage even at highest peptide concentrations, suggesting a coupling of the peptide activity to vesicle fusion, being mediated by the intrinsic lipid curvatures of PE and CL. Our results thus show that LF11-215 and LF11-324 effectively target inner bacterial components, while the stored elastic stress makes membranes more vulnerable to peptide translocation.

Ebola virus glycoprotein interacts with cholesterol to enhance membrane fusion and cell entry

Cholesterol serves critical roles in enveloped virus fusion by modulating membrane properties. The glycoprotein (GP) of Ebola virus (EBOV) promotes fusion in the endosome, a process that requires the endosomal cholesterol transporter NPC1. However, the role of cholesterol in EBOV fusion is unclear. Here we show that cholesterol in GP-containing membranes enhances fusion and the membrane-proximal external region and transmembrane (MPER/TM) domain of GP interacts with cholesterol via several glycine residues in the GP2 TM domain, notably G660. Compared to wild-type (WT) counterparts, a G660L mutation caused a more open angle between MPER and TM domains in an MPER/TM construct, higher probability of stalling at hemifusion for GP2 proteoliposomes and lower cell entry of virus-like particles (VLPs). VLPs with depleted cholesterol show reduced cell entry, and VLPs produced under cholesterol-lowering statin conditions show less frequent entry than respective controls. We propose that cholesterol-TM interactions affect structural features of GP2, thereby facilitating fusion and cell entry.

Mechanistic insights of host cell fusion of SARS-CoV-1 and SARS-CoV-2 from atomic resolution structure and membrane dynamics

The emerging and re-emerging viral diseases are continuous threats to the wellbeing of human life. Previous outbreaks of Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS had evidenced potential threats of coronaviruses in human health. The recent pandemic due to SARS-CoV-2 is overwhelming and has been going beyond control. Vaccines and antiviral drugs are ungently required to mitigate the pandemic. Therefore, it is important to comprehend the mechanistic details of viral infection process. The fusion between host cell and virus being the first step of infection, understanding the fusion mechanism could provide crucial information to intervene the infection process. Interestingly, all enveloped viruses contain fusion protein on their envelope that acts as fusion machine. For coronaviruses, the spike or S glycoprotein mediates successful infection through receptor binding and cell fusion. The cell fusion process requires merging of virus and host cell membranes, and that is essentially performed by the S2 domain of the S glycoprotein. In this review, we have discussed cell fusion mechanism of SARS-CoV-1 from available atomic resolution structures and membrane binding of fusion peptides. We have further discussed about the cell fusion of SARS-CoV-2 in the context of present pandemic situation.

Rate of Entropy Production in Evolving Interfaces and Membranes under Astigmatic Kinematics: Shape Evolution in Geometric-Dissipation Landscapes

This paper presents theory and simulation of viscous dissipation in evolving interfaces and membranes under kinematic conditions, known as astigmatic flow, ubiquitous during growth processes in nature. The essential aim is to characterize and explain the underlying connections between curvedness and shape evolution and the rate of entropy production due to viscous bending and torsion rates. The membrane dissipation model used here is known as the Boussinesq-Scriven fluid model. Since the standard approaches in morphological evolution are based on the average, Gaussian and deviatoric curvatures, which comingle shape with curvedness, this paper introduces a novel decoupled approach whereby shape is independent of curvedness. In this curvedness-shape landscape, the entropy production surface under constant homogeneous normal velocity decays with growth but oscillates with shape changes. Saddles and spheres are minima while cylindrical patches are maxima. The astigmatic flow trajectories on the entropy production surface, show that only cylinders and spheres grow under the constant shape. Small deviations from cylindrical shapes evolve towards spheres or saddles depending on the initial condition, where dissipation rates decrease. Taken together the results and analysis provide novel and significant relations between shape evolution and viscous dissipation in deforming viscous membrane and surfaces.

Vesicle formation mechanisms: an overview

Vesicle structures primarily embody spherical capsules composed of a single or multiple bilayers, entrapping a pool of aqueous solution in their interior. The bilayers can be synthesised by phospholipids or other amphiphiles (surfactants, block copolymers, etc.). Vesicles with broad-spectrum applications in numerous scientific disciplines, including biochemistry, biophysics, biology, and various pharmaceutical industries, have attracted widespread attention. Consequently, a multitude of protocols have been devised and proposed for their fabrication. In this review, with a motivation to derive the basic conditions for the formation of vesicles, the associated thermodynamic and kinetic aspects are comprehensively appraised. Contextually, an all-purpose overview of the underlying thermodynamics of bilayer/membrane generation and deformation, including the chemical potential of aggregates, geometric packing and the concept of elastic properties, is presented. Additionally, the current review highlights the probable, inherent mechanisms of vesicle formation under distinct modes of manufacturing. We lay focus on vesicle formation from pre-existing bilayers, as well as from bilayers, which form when lipids from an organic solvent are transferred into an aqueous medium. Furthermore, we outline the kinetic effects on vesicle formation from the lamellar phase, with and without the presence of shearing force. Wherever required, the experimental and/or theoretical outcomes, the driving forces for vesicle size selection, and various scaling laws are also reviewed, all of which facilitate an overall improved understanding of the vesicle formation mechanisms.

Mechanical properties of lipid bilayers: a note on the Poisson ratio

We investigate the Poisson ratio ν of fluid lipid bilayers, i.e., the question how area strains compare to the changes in membrane thickness (or, equivalently, volume) that accompany them. We first examine existing experimental results on the area- and volume compressibility of lipid membranes. Analyzing them within the framework of linear elasticity theory for homogeneous thin fluid sheets leads us to conclude that lipid membrane deformations are to a very good approximation volume-preserving, with a Poisson ratio that is likely about 3% smaller than the common soft matter limit . These results are fully consistent with atomistic simulations of a DOPC membrane at varying amount of applied lateral stress, for which we instead deduce ν by directly comparing area- and volume strains. To assess the problematic assumption of transverse homogeneity, we also define a depth-resolved Poisson ratio ν(z) and determine it through a refined analysis of the same set of simulations. We find that throughout the membrane's thickness, ν(z) is close to the value derived assuming homogeneity, with only minor variations of borderline statistical significance.

Boundary behaviour of open vesicles in axisymmetric case

A continuous transformation from a closed vesicle to an open vesicle requires that the area of open hole enlarges from zero. Since the shape equation and boundary conditions of lipid open vesicles with free edges have been obtained, we want to know whether this process can be achieved with valid parameters. By studying the boundary conditions in the axisymmetric case, the analytic expression of the boundary edges is obtained generally. It reveals that the radius and line tension of boundary edges are confined strongly by bending moduli. In some cases, there is the minimal nonzero boundary radius and the line tension needs to surmount the maxim following the increase of boundary radius. Without the spontaneous curvature, the line tension will trend to infinite when the boundary radius shrinks to zero. The continuous opening up process requires that the spontaneous curvature is nonzero and the ratio between the bending moduli of Gauss curvature and mean curvature satisfies [Formula: see text], which is smaller than the value [Formula: see text] from experiments and simulations. This result indicates that the opening up process is discontinuous.

All-Factor Analysis and Correlations on the Transmembrane Process for Arginine-Rich Cell-Penetrating Peptides

Currently, arginine-rich cell-penetrating peptides (CPPs), due to their little cytotoxicity and high transmembrane efficiency, are considered as one of the important intracellular carriers. Although the mechanism of the transmembrane process for arginine-rich CPPs was proposed, the quantitative correlations and the key factors involved in this process still deserve further investigation. In this study, all-atom molecular dynamics and the umbrella sampling technique were employed to study the arginine-rich CPPs transmembrane process. In the adsorption process of CPPs from solution to the surface of the lipid bilayer, the adsorption free energy (ΔG_A) is found to be linearly related to the interaction energy change (ΔE_A): ΔG_A = 0.0426ΔE_A + 36.7, R² = 0.92. In the CPPs transmembrane process, the transmembrane free energy barrier (ΔG_B) is roughly correlated with the corresponding interaction energy change (ΔE_B): ΔG_B = 0.108ΔE_B +135, R² = 0.73. The multiple salt bridges of guanidinium-PO₄ account for 65% of the overall interaction energy, so the increased negative charges of the lipid bilayer or more salt bridges would facilitate CPPs adsorption and transmembrane processes. Also, the increased negative charges of the lipid bilayer would reduce the amount of water to be carried into the pore and further reduce the ΔG_B. The peptide backbone would not have a direct impact on transmembrane efficiency. The ΔG_B is also found to be related to the length of the pore (L): ΔG_B = 46.2L - 31.3, R² = 0.92, which makes the transmembrane efficiency estimable. This work is expected to deliver an in-depth understanding and help the optimization of CPPs.

Simulation Best Practices for Lipid Membranes [Article v1.0]

We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.

Molecular mechanism of vesicle division induced by coupling between lipid geometry and membrane curvatures

We investigated the effects of lipid geometry on vesicle division using coarse grained molecular dynamics simulations. When the vesicle is composed of zero and negative spontaneous curvature lipids (ZSLs and NSLs), the difference in their molecular spontaneous curvatures destabilizes the neck of the limiting shape vesicle. In the vesicle division pathway, the neck developed into the stalk intermediates. The stalk was broken when the NSLs were expelled from the stalk. Free energy analysis shows that the coupling between the lipid geometry and the Gaussian rigidity is responsible for the observed vesicle division.

Computing Curvature Sensitivity of Biomolecules in Membranes by Simulated Buckling

Membrane curvature sensing, where the binding free energies of membrane-associated molecules depend on the local membrane curvature, is a key factor to modulate and maintain the shape and organization of cell membranes. However, the microscopic mechanisms are not well understood, partly due to absence of efficient simulation methods. Here, we describe a method to compute the curvature dependence of the binding free energy of a membrane-associated probe molecule that interacts with a buckled membrane, which has been created by lateral compression of a flat bilayer patch. This buckling approach samples a wide range of curvatures in a single simulation, and anisotropic effects can be extracted from the orientation statistics. We develop an efficient and robust algorithm to extract the motion of the probe along the buckled membrane surface, and evaluate its numerical properties by extensive sampling of three coarse-grained model systems: local lipid density in a curved environment for single-component bilayers, curvature preferences of individual lipids in two-component membranes, and curvature sensing by a homotrimeric transmembrane protein. The method can be used to complement experimental data from curvature partition assays and provides additional insight into mesoscopic theories and molecular mechanisms for curvature sensing.

Pulling force and surface tension drive membrane fusion

Despite catalyzed by fusion proteins of quite different molecular architectures, intracellular, viral, and cell-to-cell fusions are found to have the essential common features and the nearly same nature of transition states. The similarity inspires us to find a more general catalysis mechanism for membrane fusion that minimally depends on the specific structures of fusion proteins. In this work, we built a minimal model for membrane fusion, and by using dissipative particle dynamics simulations, we propose a mechanism that the pulling force generated by fusion proteins initiates the fusion process and the membrane tension regulates the subsequent fusion stages. The model shows different features compared to previous computer simulation studies: the pulling force catalyzes membrane fusion through lipid head overcrowding in the contacting region, leading to an increase in the head-head repulsion and/or the unfavorable head-tail contacts from opposing membranes, both of which destabilize the contacting leaflets and thus promote membrane fusion or vesicle rupture. Our simulations produce a variety of shapes and intermediates, closely resembling cases seen experimentally. Our work strongly supports the view that the tight pulling mechanism is a conserved feature of fusion protein-mediated fusion and that the membrane tension plays an essential role in fusion.

Lipid Phase Control and Secondary Structure of Viral Fusion Peptides Anchored in Monoolein Membranes

The fusion of lipid membranes involves major changes of the membrane curvatures and is mediated by fusion proteins that bind to the lipid membranes. For a better understanding of the way fusion proteins steer this process, we have studied the interaction of two different viral fusion peptides, HA2-FP and TBEV-FP, with monoolein mesophases as a function of temperature and pressure at limited hydration. The fusion peptides are derived from the influenza virus hemagglutinin fusion protein (HA2-FP) and from the tick-borne encephalitis virus envelope glycoprotein E (TBEV-FP). By using synchrotron X-ray diffraction, the changes of the monoolein phase behavior upon binding the peptides have been determined and the concomitant secondary structures of the peptides have been analyzed by FTIR spectroscopy. As main results we have found that the fusion peptides interact differently with monoolein and change the pressure and temperature dependent lipid phase behavior to different extents. However, they both destabilize the fluid lamellar phase and favor phases with negative curvature, i.e. inverse bicontinuous cubic and inverse hexagonal phases. These peptide-induced phase changes can partially be reversed by the application of high pressure, demonstrating that the promotion of negative curvature is achieved by a less dense packing of the monoolein membranes by the fusion peptides. Pressure jumps across the cubic-lamellar phase transition reveal that HA2-FP has a negligible effect on the rates of the cubic and the lamellar phase formation. Interestingly, the secondary structures of the fusion peptides appear unaffected by monoolein fluid-fluid phase transitions, suggesting that the fusion peptides are the structure dominant species in the fusion process of lipid membranes.

Determining the Gaussian Modulus and Edge Properties of 2D Materials: From Graphene to Lipid Bilayers

The dominant deformation behavior of two-dimensional materials (bending) is primarily governed by just two parameters: bending rigidity and the Gaussian modulus. These properties also set the energy scale for various important physical and biological processes such as pore formation, cell fission and generally, any event accompanied by a topological change. Unlike the bending rigidity, the Gaussian modulus is, however, notoriously difficult to evaluate via either experiments or atomistic simulations. In this Letter, recognizing that the Gaussian modulus and edge tension play a nontrivial role in the fluctuations of a 2D material edge, we derive closed-form expressions for edge fluctuations. Combined with atomistic simulations, we use the developed approach to extract the Gaussian modulus and edge tension at finite temperatures for both graphene and various types of lipid bilayers. Our results possibly provide the first reliable estimate of this elusive property at finite temperatures and appear to suggest that earlier estimates must be revised. In particular, we show that, if previously estimated properties are employed, the graphene-free edge will exhibit unstable behavior at room temperature. Remarkably, in the case of graphene, we show that the Gaussian modulus and edge tension even change sign at finite temperatures.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

Achiral symmetry breaking and positive Gaussian modulus lead to scalloped colloidal membranes

In the presence of a nonadsorbing polymer, monodisperse rod-like particles assemble into colloidal membranes, which are one-rod-length-thick liquid-like monolayers of aligned rods. Unlike 3D edgeless bilayer vesicles, colloidal monolayer membranes form open structures with an exposed edge, thus presenting an opportunity to study elasticity of fluid sheets. Membranes assembled from single-component chiral rods form flat disks with uniform edge twist. In comparison, membranes composed of a mixture of rods with opposite chiralities can have the edge twist of either handedness. In this limit, disk-shaped membranes become unstable, instead forming structures with scalloped edges, where two adjacent lobes with opposite handedness are separated by a cusp-shaped point defect. Such membranes adopt a 3D configuration, with cusp defects alternatively located above and below the membrane plane. In the achiral regime, the cusp defects have repulsive interactions, but away from this limit we measure effective long-ranged attractive binding. A phenomenological model shows that the increase in the edge energy of scalloped membranes is compensated by concomitant decrease in the deformation energy due to Gaussian curvature associated with scalloped edges, demonstrating that colloidal membranes have positive Gaussian modulus. A simple excluded volume argument predicts the sign and magnitude of the Gaussian curvature modulus that is in agreement with experimental measurements. Our results provide insight into how the interplay between membrane elasticity, geometrical frustration, and achiral symmetry breaking can be used to fold colloidal membranes into 3D shapes.

Fusion pores and their control of neurotransmitter and hormone release

Chang et al. review fusion pore structure and dynamics and discuss the implications for hormone and neurotransmitter release

Ca²⁺-triggered exocytosis functions broadly in the secretion of chemical signals, enabling neurons to release neurotransmitters and endocrine cells to release hormones. The biological demands on this process can vary enormously. Although synapses often release neurotransmitter in a small fraction of a millisecond, hormone release can be orders of magnitude slower. Vesicles usually contain multiple signaling molecules that can be released selectively and conditionally. Cells are able to control the speed, concentration profile, and content selectivity of release by tuning and tailoring exocytosis to meet different biological demands. Much of this regulation depends on the fusion pore—the aqueous pathway by which molecules leave a vesicle and move out into the surrounding extracellular space. Studies of fusion pores have illuminated how cells regulate secretion. Furthermore, the formation and growth of fusion pores serve as a readout for the progress of exocytosis, thus revealing key kinetic stages that provide clues about the underlying mechanisms. Herein, we review the structure, composition, and dynamics of fusion pores and discuss the implications for molecular mechanisms as well as for the cellular regulation of neurotransmitter and hormone release.

Coarse-grained molecular dynamics study of membrane fusion: Curvature effects on free energy barriers along the stalk mechanism

The effects of membrane curvature on the free energy barrier for membrane fusion have been investigated using coarse-grained molecular dynamics (CG-MD) simulations, assuming that fusion takes place through a stalk intermediate. Free energy barriers were estimated for stalk formation as well as for fusion pore formation using the guiding potential method. Specifically, the three different geometries of two apposed membranes were considered: vesicle-vesicle, vesicle-planar, and planar-planar membranes. The free energy barriers for the resulting fusion were found to depend importantly on the fusing membrane geometries; the lowest barrier was obtained for vesicular membranes. Further, lipid sorting was observed in fusion of the mixed membranes of dimyristoyl phosphatidylcholine and dioleoyl phosphatidylethanolamine (DOPE). Specifically, DOPE molecules were found to assemble around the stalk to support the highly negative curved membrane surface. A consistent result for lipid sorting was observed when a simple continuum model (CM) was used, where the Helfrich energy and mixing entropy of the lipids were taken into account. However, the CM predicts a much higher free energy barrier than found using CG-MD. This discrepancy originates from the conformational changes of lipids, which were not considered in the CM. The results of the CG-MD simulations reveal that a large conformational change in the lipid takes place around the stalk region, which results in a reduction of free energy barriers along the stalk mechanism of membrane fusion.

Revisit the Correlation between the Elastic Mechanics and Fusion of Lipid Membranes

Membrane fusion is a vital process in key cellular events. The fusion capability of a membrane depends on its elastic properties and varies with its lipid composition. It is believed that as the composition varies, the consequent change in C₀ (monolayer spontaneous curvature) is the major factor dictating fusion, owing to the associated variation in G_Es (elastic energies) of the fusion intermediates (e.g. stalk). By exploring the correlations among fusion, C₀ and K_cp (monolayer bending modulus), we revisit this long-held belief and re-examine the fusogenic contributions of some relevant factors. We observe that not only C₀ but also K_cp variations affect fusion, with depression in K_cp leading to suppression in fusion. Variations in G_Es and inter-membrane interactions cannot account for the K_cp-fusion correlation; fusion is suppressed even as the G_Es decrease with K_cp, indicating the presence of factor(s) with fusogenic importance overtaking that of G_E. Furthermore, analyses find that the C₀ influence on fusion is effected via modulating G_E of the pre-fusion planar membrane, rather than stalk. The results support a recent proposition calling for a paradigm shift from the conventional view of fusion and may reshape our understanding to the roles of fusogenic proteins in regulating cellular fusion machineries.

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.

Function Suggests Nano-Structure: Quantitative Structural Support for SNARE-Mediated Pore Formation

Granule secretory content is released in either basal or calcium-activated complete exocytosis mode. A vital element in these processes is the establishment of a fusion pore between the granule membrane and the plasma membrane, initiated by the formation of a circular rosette docking arrangement of SNARE protein complexes. The controversially disputed number of SNARE complexes needed for granule priming leading to the formation of the fusion pore, is granule-size dependent and varies between secretion modes. Resorting to a statistical mechanics approach that views SNARE complexes and Ca(2+) ions as interacting particles, we have developed a relationship that links secretion rate to SNARE rosette size, Ca(2+) concentration and Ca(2+) ion cooperativity. Data are presented and discussed which suggest this SNARE-dependent generalization of existing narrow-range biophysical models that correlate secretion rate with Ca(2+) concentration and maximal Ca(2+) ion cooperativity. Evidence from dozens of examples in the literature advocate for this relation, which holds through the entire biological range. The coalescence of so many areas of diverse research methodologies has greatly augmented our understanding of so many different sequences of granule life cycle. Accordingly, these new tools may become valuable in a variety of electrophysiological experiments.

How could SNARE proteins open a fusion pore?

The SNARE (Soluble NSF Attachment protein REceptor) complex, which in mammalian neurosecretory cells is composed of the proteins synaptobrevin 2 (also called VAMP2), syntaxin, and SNAP-25, plays a key role in vesicle fusion. In this review, we discuss the hypothesis that, in neurosecretory cells, fusion pore formation is directly accomplished by a conformational change in the SNARE complex via movement of the transmembrane domains.

Introduction to membrane lipids

Biological membranes are composed largely of lipids and proteins. The most common arrangement of lipids in biological membranes is as a bilayer. This arrangement spontaneously forms a barrier for the passage of polar materials. The bilayer is thin but can have a large area in the dimension perpendicular to its thickness. The physical nature of the bilayer membrane will vary according to the conditions of the environment as well as the chemical structure of the lipid constituents of the bilayer. These physical properties determine the function of the membrane together with specific structural features of the lipids that allow them to have signaling properties. The lipids of the membrane are not uniformly distributed. There is an intrinsic asymmetry between the two monolayers that constitute the bilayer. In addition, some lipids tend to be enriched in particular regions of the membrane, termed domains. There is evidence that certain domains recruit specific proteins into that domain. This has been suggested to be important for allowing interaction among different proteins involved in certain signal transduction pathways. Membrane lipids have important roles in determining the physical properties of the membrane, in modulating the activity of membrane-bound proteins and in certain cases being specific secondary messengers that can interact with specific proteins. A large variety of lipids present in biological membranes result in them possessing many functions.

Fast, Ca2+-dependent exocytosis at nerve terminals: shortcomings of SNARE-based models

Investigations over the last two decades have made major inroads in clarifying the cellular and molecular events that underlie the fast, synchronous release of neurotransmitter at nerve endings. Thus, appreciable progress has been made in establishing the structural features and biophysical properties of the calcium (Ca2+) channels that mediate the entry into nerve endings of the Ca2+ ions that trigger neurotransmitter release. It is now clear that presynaptic Ca2+ channels are regulated at many levels and the interplay of these regulatory mechanisms is just beginning to be understood. At the same time, many lines of research have converged on the conclusion that members of the synaptotagmin family serve as the primary Ca2+ sensors for the action potential-dependent release of neurotransmitter. This identification of synaptotagmins as the proteins which bind Ca2+ and initiate the exocytotic fusion of synaptic vesicles with the plasma membrane has spurred widespread efforts to reveal molecular details of synaptotagmin's action. Currently, most models propose that synaptotagmin interfaces directly or indirectly with SNARE (soluble, N-ethylmaleimide sensitive factor attachment receptors) proteins to trigger membrane fusion. However, in spite of intensive efforts, the field has not achieved consensus on the mechanism by which synaptotagmins act. Concurrently, the precise sequence of steps underlying SNARE-dependent membrane fusion remains controversial. This review considers the pros and cons of the different models of SNARE-mediated membrane fusion and concludes by discussing a novel proposal in which synaptotagmins might directly elicit membrane fusion without the intervention of SNARE proteins in this final fusion step.

Free energy analysis along the stalk mechanism of membrane fusion

The free energy profile of the stalk model of membrane fusion has been calculated using coarse-grained molecular dynamics simulations. The proposed method guides the lipid configuration using a guiding wall potential to make the transition from two apposed membranes to a stalk and a fusion pore. The free energy profile is obtained with a thermodynamic integration scheme using the mean force working on the guiding wall as a response of the system. We applied the method to two apposed flat bilayers composed of dioleoyl phosphatidylethanolamine/dioleoyl phosphatidylcholine expanding over the simulation box under the periodic boundary conditions. The two transition states are identified as pre-stalk and pre-pore states. The free energy barrier for the latter is confirmed to be in good agreement with that estimated by the pulling method. The present method provides a practical way to calculate the free energy profile along the stalk mechanism.