A new mechanism of model membrane fusion determined from Monte Carlo simulation

Membrane fission via transmembrane contact

Division of intracellular organelles often correlates with additional membrane wrapping, e.g., by the endoplasmic reticulum or the outer mitochondrial membrane. Such wrapping plays a vital role in proteome and lipidome organization. However, how an extra membrane impacts the mechanics of the division has not been investigated. Here we combine fluorescence and cryo-electron microscopy experiments with self-consistent field theory to explore the stress-induced instabilities imposed by membrane wrapping in a simple double-membrane tubular system. We find that, at physiologically relevant conditions, the outer membrane facilitates an alternative pathway for the inner-tube fission through the formation of a transient contact (hemi-fusion) between both membranes. A detailed molecular theory of the fission pathways in the double membrane system reveals the topological complexity of the process, resulting both in leaky and leakless intermediates, with energies and topologies predicting physiological events.

Two modes of fusogenic action for influenza virus fusion peptide

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

Effect of polyunsaturated free fatty acids on the membrane fusion mechanism

Membrane fusion is one of the important processes for the survival of eukaryotic cells and the entry of enveloped viruses into the host cells. Lipid composition plays a crucial role by modulating the organization and dynamics of the membrane, as well as the structure and conformation of membrane proteins. The diversity of the lipid acyl chain in its length and degree of unsaturation originates from the variation in free fatty acids (FFAs). We have studied the effect of linoleic (LA) and alpha-linolenic (ALA) acids on the depth-dependent organization, dynamics, and fusion of DOPC/DOPE (70/30 mol%) membranes utilizing steady-state and time-resolved fluorescence spectroscopic methods. Our results suggest that membranes with 5 mol% LA stabilize the stalk-intermediate and promote lipid mixing at the early stage of the process, i.e., the fusion follows the classical stalk model. Conversely, the extents of lipid and content mixing at the stalk intermediate are similar in the presence of 5 mol% of ALA, indicating the fusion mechanism as a nonclassical one like in the DOPC/DOPE (70/30 mol%) membranes. Our results provide an in-depth insight into the effect of the increasing degree of fatty acid tail unsaturation on membrane organization and dynamics and their impact on the membrane fusion mechanism.

Computational methods to study enveloped viral entry

The protein-membrane interactions that mediate viral infection occur via loosely ordered, transient assemblies, creating challenges for high-resolution structure determination. Computational methods and in particular molecular dynamics simulation have thus become important adjuncts for integrating experimental data, developing mechanistic models, and suggesting testable hypotheses regarding viral function. However, the large molecular scales of virus-host interaction also create challenges for detailed molecular simulation. For this reason, continuum membrane models have played a large historical role, although they have become less favored for high-resolution models of protein assemblies and lipid organization. Here, we review recent progress in the field, with an emphasis on the insight that has been gained using a mixture of coarse-grained and atomic-resolution molecular dynamics simulations. Based on successes and challenges to date, we suggest a multiresolution strategy that should yield the best mixture of computational efficiency and physical fidelity. This strategy may facilitate further simulations of viral entry by a broader range of viruses, helping illuminate the diversity of viral entry strategies and the essential common elements that can be targeted for antiviral therapies.

Spontaneous Formation and Fusion of Raspberry Vesicle Self-Assembled from Star Block Terpolymers in Aqueous Solution

The spontaneous formation and fusion of raspberry vesicles was studied using the dissipative particle dynamics (DPD) method. The vesicles were formed through the self-assembly of amphiphilic E₁₂O₆F₂ star terpolymers in selective solvent. E and F blocks are solvophobic and the O block is solvophilic. The shortest F block plays a major role in the formation of raspberry vesicles. Distinct vesicle formation mechanisms were observed at different polymer concentrations. At higher concentrations, vesicles form via the bending and closure of an oblate F-bump-E bilayer. At lower concentrations, the formation pathway contains: the initial formation of a vesicle with a core, the combination of such vesicles into cylindrical micelles, and the bending of the cylindrical micelles to form a hollow vesicle. In addition, raspberry vesicle fusion is regulated by F bumps through the continuous coalescence of them from apposed vesicle membranes. The contact area bends, followed by the formation of a fusion pore and a tilted inner layer. As the pore sealed, the hemifusion structure appears, which further restructures to form a vesicle. Our results provide guidance on understanding the dynamic processes of complex vesicles and biological membrane fusion.

Effect of Phosphatidylethanolamine and Oleic Acid on Membrane Fusion: Phosphatidylethanolamine Circumvents the Classical Stalk Model

Membrane fusion is one of the most important processes for the survival of eukaryotic cells and entry of enveloped viruses to the host cells. Lipid composition plays a crucial role in the process by modulating the organization and dynamics of the membrane, as well as the structure and conformation of membrane proteins. Phosphatidylethanolamine (PE), a lipid molecule with intrinsic negative curvature, promotes membrane fusion by stabilizing the non-lamellar intermediate structures in the fusion process. Conversely, oleic acid (OA), with intrinsic positive curvature, inhibits membrane fusion. The current study aimed to investigate polyethylene glycol-mediated lipid mixing, content mixing, content leakage, and depth-dependent membrane organization and dynamics, using arrays of steady-state and time-resolved fluorescence techniques, to determine the causative role of PE and OA in membrane fusion. The results demonstrated that the presence of 30 mol % PE in the membrane promotes membrane fusion through a mechanism that circumvents the classical stalk model. On the contrary, membranes containing OA showed reduced rate and extent of fusion, despite following the same mechanism. Collectively, our findings in terms of membrane organization and dynamics indicated a plausible role of PE and OA in membrane fusion.

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.

Continuum Models of Membrane Fusion: Evolution of the Theory

Starting from fertilization, through tissue growth, hormone secretion, synaptic transmission, and sometimes morbid events of carcinogenesis and viral infections, membrane fusion regulates the whole life of high organisms. Despite that, a lot of fusion processes still lack well-established models and even a list of main actors. A merger of membranes requires their topological rearrangements controlled by elastic properties of a lipid bilayer. That is why continuum models based on theories of membrane elasticity are actively applied for the construction of physical models of membrane fusion. Started from the view on the membrane as a structureless film with postulated geometry of fusion intermediates, they developed along with experimental and computational techniques to a powerful tool for prediction of the whole process with molecular accuracy. In the present review, focusing on fusion processes occurring in eukaryotic cells, we scrutinize the history of these models, their evolution and complication, as well as open questions and remaining theoretical problems. We show that modern approaches in this field allow continuum models of membrane fusion to stand shoulder to shoulder with molecular dynamics simulations, and provide the deepest understanding of this process in multiple biological systems.

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

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

Nanotube-Mediated Path to Protocell Formation

Cellular compartments are membrane-enclosed, spatially distinct microenvironments that confine and protect biochemical reactions in the biological cell. On the early Earth, the autonomous formation of compartments is thought to have led to the encapsulation of nucleotides, thereby satisfying a starting condition for the emergence of life. Recently, surfaces have come into focus as potential platforms for the self-assembly of prebiotic compartments, as significantly enhanced vesicle formation was reported in the presence of solid interfaces. The detailed mechanism of such formation at the mesoscale is still under discussion. We report here on the spontaneous transformation of solid-surface-adhered lipid deposits to unilamellar membrane compartments through a straightforward sequence of topological changes, proceeding via a network of interconnected lipid nanotubes. We show that this transformation is entirely driven by surface-free energy minimization and does not require hydrolysis of organic molecules or external stimuli such as electrical currents or mechanical agitation. The vesicular structures take up and encapsulate their external environment during formation and can subsequently separate and migrate upon exposure to hydrodynamic flow. This may link the self-directed transition from weakly organized bioamphiphile assemblies on solid surfaces to protocells with secluded internal contents.

Membrane Composition Modulates Fusion by Altering Membrane Properties and Fusion Peptide Structure

Membrane fusion, one of the most essential processes in the life of eukaryotes, occurs when two separate lipid bilayers merge into a continuous bilayer and internal contents of two separated membranes mingle. There is a certain class of proteins that assist the binding of the viral envelope to the target host cell and catalyzing fusion. All class I viral fusion proteins contain a highly conserved 20–25 amino-acid amphipathic peptide at the N-terminus, which is essential for fusion activity and is termed as the ‘fusion peptide’. It has been shown that insertion of fusion peptides into the host membrane and the perturbation in the membrane generated thereby is crucial for membrane fusion. Significant efforts have been given in the last couple of decades to understand the lipid-dependence of structure and function of the fusion peptide in membranes to understand the role of lipid compositions in membrane fusion. In addition, the lipid compositions further change the membrane physical properties and alter the mechanism and extent of membrane fusion. Therefore, lipid compositions modulate membrane fusion by changing membrane physical properties and altering structure of the fusion peptide.

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

Biological membrane fusion proceeds via an essential topological transition of the two membranes involved. Known players such as certain lipid species and fusion proteins are generally believed to alter the free energy and thus the rate of the fusion reaction. Quantifying these effects by theory poses a major challenge since the essential reaction intermediates are collective, diffusive and of a molecular length scale. We conducted molecular dynamics simulations in conjunction with a state-of-the-art string method to resolve the minimum free-energy path of the first fusion intermediate state, the so-called stalk. We demonstrate that the isolated transmembrane domains (TMDs) of fusion proteins such as SNARE molecules drastically lower the free energy of both the stalk barrier and metastable stalk, which is not trivially explained by molecular shape arguments. We relate this effect to the local thinning of the membrane (negative hydrophobic mismatch) imposed by the TMDs which favors the nearby presence of the highly bent stalk structure or prestalk dimple. The distance between the membranes is the most crucial determinant of the free energy of the stalk, whereas the free-energy barrier changes only slightly. Surprisingly, fusion enhancing lipids, i.e., lipids with a negative spontaneous curvature, such as PE lipids have little effect on the free energy of the stalk barrier, likely because of its single molecular nature. In contrast, the lipid shape plays a crucial role in overcoming the hydration repulsion between two membranes and thus rather lowers the total work required to form a stalk.

The role of scaffold reshaping and disassembly in dynamin driven membrane fission

The large GTPase dynamin catalyzes membrane fission in eukaryotic cells, but despite three decades of experimental work, competing and partially conflicting models persist regarding some of its most basic actions. Here we investigate the mechanical and functional consequences of dynamin scaffold shape changes and disassembly with the help of a geometrically and elastically realistic simulation model of helical dynamin-membrane complexes. Beyond changes of radius and pitch, we emphasize the crucial role of a third functional motion: an effective rotation of the filament around its longitudinal axis, which reflects alternate tilting of dynamin’s PH binding domains and creates a membrane torque. We also show that helix elongation impedes fission, hemifission is reached via a small transient pore, and coat disassembly assists fission. Our results have several testable structural consequences and help to reconcile mutual conflicting aspects between the two main present models of dynamin fission—the two-stage and the constrictase model.

When cells take up material from their surroundings, they must first transport this cargo across their outer membrane, a flexible sheet of tightly organized fat molecules that act as a barrier to the environment. Cells can achieve this by letting their membrane surround the object, pulling it inwards until it is contained in a pouch that bulges into the cell. This bag is then corded up so it splits off from the outer membrane. The ‘cord’ is a protein called dynamin, which is thought to form a tight spiral around the bag’s neck, closing it over and pinching it away. The structure of dynamin is fairly well known, and yet several theories compete to explain how it may snap the bag off the outer membrane.

Here, Pannuzzo et al. have created a computer simulation that faithfully replicates the geometry and the elasticity of the membrane and of dynamin, and used it to test different ways the protein could work. The first test featured simple constriction, where the dynamin spiral contracts around the membrane to pinch it; this only separated the bag from the membrane after implausibly tight constriction. The second test added elongation, with the spiral lengthening as well as reducing its diameter, but this further reduced the ability for the protein to snap off the membrane. The final test combined constriction and rotation, whereby dynamin ‘twirls’ as it presses on the neck of the bag: this succeeded in efficiently severing the membrane once the dynamin spiral disassembled. Indeed, the simulations suggested that dynamin might start to dismantle while it constricts, without compromising its role. In fact, getting rid of excess length as the protein contracts helps to dissolve any remnants of a membrane connection.

Defects in dynamin are associated with conditions such as centronuclear myopathy and Charcot‐Marie‐Tooth peripheral neuropathy. Recent research also indicates that the protein is involved in a much wider range of neurological disorders that include Alzheimer's, Parkinson's, Huntington's, and amyotrophic lateral sclerosis. The models created by Pannuzzo et al. are useful tools to understand how dynamin and similar proteins work and sometimes fail.

Lipid-dependence of target membrane stability during influenza viral fusion

Although influenza kills about a half million people each year, even after excluding pandemics, there is only one set of antiviral drugs: neuraminidase inhibitors. By using a new approach utilizing giant unilamellar vesicles and infectious X-31 influenza virus, and testing for the newly identified pore intermediate of membrane fusion, we observed ∼30-87% poration, depending upon lipid composition. Testing the hypothesis that spontaneous curvature (SC) of the lipid monolayer controls membrane poration, our Poisson model and Boltzmann energetic considerations suggest a transition from a leaky to a non-leaky fusion pathway depending on the SC of the target membrane. When the target membrane SC is below approximately -0.20 nm^-1 fusion between influenza virus and target membrane is predominantly non-leaky while above that fusion is predominantly leaky, suggesting that influenza hemagglutinin (HA)-catalyzed topological conversion of target membranes during fusion is associated with a loss of membrane integrity.

Cholesterol suppresses membrane leakage by decreasing water penetrability

Membrane fusion is a fundamental biological process that lies at the heart of enveloped virus infection, synaptic signaling, intracellular vesicle trafficking, gamete fertilization, and cell-cell fusion. Membrane fusion is initiated as two apposed membranes merge to a single bilayer called a hemifusion diaphragm. It is believed that the contents of the two fusing membranes are released through a fusion pore formed at the hemifusion diaphragm, and yet another possible pathway has been proposed in which an undefined pore may form outside the hemifusion diaphragm at the apposed membranes, leading to the so-called leaky fusion. Here, we performed all-atom molecular dynamics simulations to study the evolution of the hemifusion diaphragm structure with various lipid compositions. We found that the lipid cholesterol decreased water penetrability to inhibit leakage pore formation. Biochemical leakage experiments support these simulation results. This study may shed light on the underlying mechanism of the evolution pathways of the hemifusion structure, especially the understanding of content leakage during membrane fusion.

Rupturing the hemi-fission intermediate in membrane fission under tension: Reaction coordinates, kinetic pathways, and free-energy barriers

Membrane fission is a fundamental process in cells, involved inter alia in endocytosis, intracellular trafficking, and virus infection. Its underlying molecular mechanism, however, is only incompletely understood. Recently, experiments and computer simulation studies have revealed that dynamin-mediated membrane fission is a two-step process that proceeds via a metastable hemi-fission intermediate (or wormlike micelle) formed by dynamin's constriction. Importantly, this hemi-fission intermediate is remarkably metastable, i.e., its subsequent rupture that completes the fission process does not occur spontaneously but requires additional, external effects, e.g., dynamin's (unknown) conformational changes or membrane tension. Using simulations of a coarse-grained, implicit-solvent model of lipid membranes, we investigate the molecular mechanism of rupturing the hemi-fission intermediate, such as its pathway, the concomitant transition states, and barriers, as well as the role of membrane tension. The membrane tension is controlled by the chemical potential of the lipids, and the free-energy landscape as a function of two reaction coordinates is obtained by grand canonical Wang-Landau sampling. Our results show that, in the course of rupturing, the hemi-fission intermediate undergoes a "thinning → local pinching → rupture/fission" pathway, with a bottle-neck-shaped cylindrical micelle as a transition state. Although an increase of membrane tension facilitates the fission process by reducing the corresponding free-energy barrier, for biologically relevant tensions, the free-energy barriers still significantly exceed the thermal energy scale k_BT.

Metastable Prepores in Tension-Free Lipid Bilayers

The formation and closure of aqueous pores in lipid bilayers is a key step in various biophysical processes. Large pores are well described by classical nucleation theory, but the free-energy landscape of small, biologically relevant pores has remained largely unexplored. The existence of small and metastable "prepores" was hypothesized decades ago from electroporation experiments, but resolving metastable prepores from theoretical models remained challenging. Using two complementary methods-atomistic simulations and self-consistent field theory of a minimal lipid model-we determine the parameters for which metastable prepores occur in lipid membranes. Both methods consistently suggest that pore metastability depends on the relative volume ratio between the lipid head group and lipid tails: lipids with a larger head-group volume fraction (or shorter saturated tails) form metastable prepores, whereas lipids with a smaller head-group volume fraction (or longer unsaturated tails) form unstable prepores.

Ultraviolet-responsive self-assembled metallomicelles for photocontrollable catalysis of asymmetric sulfoxidation in water

Self-assembled metallomicelles with ultraviolet (UV)-controlled morphologies were constructed from a synthesized azobenzene-containing amphiphilic chiral salen Ti^IV catalyst.

Spontaneous onion shape vesicle formation and fusion of comb-like block copolymers studied by dissipative particle dynamics

The formation of an onion shape vesicle.

High Transmembrane Voltage Raised by Close Contact Initiates Fusion Pore

Membrane fusion lies at the heart of neuronal communication but the detailed mechanism of a critical step, fusion pore initiation, remains poorly understood. Here, through atomistic molecular dynamics simulations, a transient pore formation induced by a close contact of two apposed bilayers is firstly reported. Such a close contact gives rise to a high local transmembrane voltage that induces the transient pore formation. Through simulations on two apposed bilayers fixed at a series of given distances, the process in which two bilayers approaching to each other under the pulling force from fusion proteins for membrane fusion was mimicked. Of note, this close contact induced fusion pore formation is contrasted with previous reported electroporation under ad hoc applied external electric field or ionic charge in-balance. We show that the transmembrane voltage increases with the decrease of the distance between the bilayers. Below a critical distance, depending on the lipid composition, the local transmembrane voltage can be sufficiently high to induce the transient pores. The size of these pores is approximately 1~2 nm in diameter, which is large enough to allow passing of neurotransmitters. A resealing of the membrane pores resulting from the neutralization of the transmembrane voltage by ions through the pores was then observed. We also found that the membrane tension can either prolong the lifetime of transient pores or cause them to dilate for full collapse. This result provides a possible mechanism for fusion pore formation and regulation of pathway of fusion process.

Size-dependent protein segregation at membrane interfaces

Membrane interfaces formed at cell-cell junctions are associated with characteristic patterns of membrane protein organization, such as E-cadherin enrichment in epithelial junctional complexes and CD45 exclusion from the signaling foci of immunological synapses. To isolate the role of protein size in these processes, we reconstituted membrane interfaces in vitro using giant unilamellar vesicles decorated with synthetic binding and non-binding proteins. We show that size differences between binding and non-binding proteins can dramatically alter their organization at membrane interfaces in the absence of active contributions from the cytoskeleton, with as little as a ~5 nm increase in non-binding protein size driving its exclusion from the interface. Combining in vitro measurements with Monte Carlo simulations, we find that non-binding protein exclusion is also influenced by lateral crowding, binding protein affinity, and thermally-driven membrane height fluctuations that transiently limit access to the interface. This simple, sensitive, and highly effective means of passively segregating proteins has implications for signaling at cell-cell junctions and protein sorting at intracellular contact points between membrane-bound organelles.

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.

pH Alters PEG-mediated fusion of phosphatidylethanolamine-containing vesicles

Here, we examine the different mechanisms of poly(ethylene glycol)-mediated fusion of small unilamellar vesicles composed of dioleoylphosphatidylcholine/dioleoylphosphatidylethanolamine (DOPE)/sphingomyelin/cholesterol in a molar ratio of 35:30:15:20 at pH 7.4 versus pH 5. In doing so, we test the hypothesis that fusion of this lipid mixture should be influenced by differences in hydration of DOPE at these two pH values. An examination of the literature reveals that DOPE should be less hydrated at pH 5 (where influenza virus particles fuse with endosome membranes) than at pH 7.4 (where synaptic vesicles or HIV virus particles fuse with plasma membrane). Ensemble kinetic experiments revealed substantial differences in fusion of this plasma membrane mimetic system at these two pH values. The most dramatic difference was the observation of two intermediates at pH 5 but loss of one of these fusion intermediates at pH 7.4. Analysis of data collected at several temperatures also revealed that formation of the initial fusion intermediate (stalk) was favored at pH 7.4 due to increased activation entropy. Our observations support the hypothesis that the different negative intrinsic curvature of DOPE can account for different fusion paths and activation thermodynamics in steps of the fusion process at these two pH values. Finally, the effects of 2 mol % hexadecane on fusion at both pH values seemed to have similar origins for step 1 (promotion of acyl chain or hydrocarbon excursion into interbilayer space) and step 3 (reduction of interstice energy leading to expansion to a critical stalk radius). Different hexadecane effects on activation thermodynamics at these two pH values can also be related to altered DOPE hydration. The results support our kinetic model for fusion and offer insight into the critical role of phosphatidylethanolamine in fusion.

Wild-type and mutant hemagglutinin fusion peptides alter bilayer structure as well as kinetics and activation thermodynamics of stalk and pore formation differently: mechanistic implications

Viral fusion peptides are short N-terminal regions of type-1 viral fusion proteins that are critical for virus entry. Although the importance of viral fusion peptides in virus-cell membrane fusion is established, little is known about how they function. We report the effects of wild-type (WT) hemagglutinin (HA) fusion peptide and its G1S, G1V, and W14A mutants on the kinetics of poly(ethylene glycol)(PEG)-mediated fusion of small unilamellar vesicles composed of dioleoylphosphatidylcholine, dioleoylphosphatidylethanolamine, sphingomyelin, and cholesterol (molar ratio of 35:30:15:20). Time courses of lipid mixing, content mixing, and content leakage were obtained using fluorescence assays at multiple temperatures and analyzed globally using either a two-step or three-step sequential ensemble model of the fusion process to obtain the rate constant and activation thermodynamics of each step. We also monitored the influence of peptides on bilayer interfacial order, acyl chain order, bilayer free volume, and water penetration. All these data were considered in terms of a recently published mechanistic model for the thermodynamic transition states for each step of the fusion process. We propose that WT peptide catalyzes Step 1 by occupying bilayer regions vacated by acyl chains that protrude into interbilayer space to form the Step 1 transition state. It also uniquely contributes a positive intrinsic curvature to hemi-fused leaflets to eliminate Step 2 and catalyzes Step 3 by destabilizing the highly stressed edges of the hemi-fused microstructures that dominate the ensemble of the intermediate state directly preceding fusion pore formation. Similar arguments explain the catalytic and inhibitory properties of the mutant peptides and support the hypothesis that the membrane-contacting fusion peptide of HA fusion protein is key to its catalytic activity.

Expansion of the fusion stalk and its implication for biological membrane fusion

Over the past 20 years, it has been widely accepted that membrane fusion proceeds via a hemifusion step before opening of the productive fusion pore. An initial hourglass-shaped lipid structure, the fusion stalk, is formed between the adjacent membrane leaflets (cis leaflets). It remains controversial if and how fusion proteins drive the subsequent transition (expansion) of the stalk into a fusion pore. Here, we propose a comprehensive and consistent thermodynamic understanding in terms of the underlying free-energy landscape of stalk expansion. We illustrate how the underlying free energy landscape of stalk expansion and the concomitant pathway is altered by subtle differences in membrane environment, such as leaflet composition, asymmetry, and flexibility. Nonleaky stalk expansion (stalk widening) requires the formation of a critical trans-leaflet contact. The fusion machinery can mechanically enforce trans-leaflet contact formation either by directly enforcing the trans-leaflets in close proximity, or by (electrostatically) condensing the area of the cis leaflets. The rate of these fast fusion reactions may not be primarily limited by the energetics but by the forces that the fusion proteins are able to exert.

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.

The fusion mechanism of small polymersomes formed by rod-coil diblock copolymers

The fusion mechanism of polymersomes self-assembled by rod-coil copolymers is intrinsically different from that of liposomes due to the effect of chain topology on conformational entropy and molecular packing. The influences of membrane tension, coil-block length, rod-block length, mutual compatibility between the solvent and the rod-coil block, and π-π interaction strength on the fusion pathway are explored by dissipative particle dynamics. The fusion process of spontaneously formed polymersomes generally consists of four stages. In the kissing stage, hopping of rod-blocks forms a connection between two vesicles of a one-legged rod-coil copolymer. In the adhesion stage, a stalk is developed by a few link-up rods and then a stretched diaphragm with rods lying parallel to the stretching direction is formed in the hemi-fusion stage. Eventually, a pore is developed and expanded in the fusion stage. If the membrane tension (τ) is adjusted by deflation/inflation of the polymersomes, the hemi-fusion diaphragm disappears. As τ is reduced, multiple stalks take shape and lead to the formation of inverted micelles, which is the rate-determining step and raises the fusion time substantially. As τ is elevated, a neck is developed after the stalk formation. The fusion process is significantly accelerated. τ of spontaneously formed vesicles varies with the coil-block length, rod-block length, solvent quality, and π-π interaction strength. There exists a critical value of τ below which the fusion process cannot be completed and a hemi-fused polymersome is formed. In addition to τ, the anisotropic steric interactions within the rod layers also resist hopping of longer rod-blocks. The coil layers develop a barrier impeding fusion between vesicles with longer coil-blocks. Consequently, lowering the solvent quality for the coil-block or rod-block facilities the fusion process due to the formation of a thinner coil layer.

A novel assay for detecting fusion pore formation: implications for the fusion mechanism

Membrane fusion is broadly envisioned as a two- or three-step process proceeding from contacting bilayers through one or two semistable, nonlamellar lipidic intermediate structures to a fusion pore. A true fusion event requires mixing of contents between compartments and is monitored by the movement of soluble molecules between trapped compartments. We have used poly(ethylene glycol) (PEG) to rapidly generate an ensemble aggregated state A that proceeds sequentially through intermediates (I₁ and/or I₂) to a final fusion pore state (FP) with rate constants k₁, k₂, and k₃. Movement of moderately sized solutes (e.g., Tb³⁺/dipicolinic acid) has been used to detect pores assigned to intermediate states as well as to the final state (FP). Analysis of ensemble kinetic data has required that mixing of contents occurs with defined probabilities (αi) in each ensemble state, although it is unclear whether pores that form in different states are different. We introduce here a simple new assay that employs fluorescence resonance energy transfer (FRET) between a 6-carboxyfluorescein (donor) and tetramethylrhodamine (acceptor), which are covalently attached to complementary sequences of 10 bp oligonucleotides. Complementary sequences of fluorophore-labeled oligonucleotides were incorporated in vesicles separately, and the level of FRET increased in a simple exponential fashion during PEG-mediated fusion. The resulting rate constant corresponded closely to the slow rate constant of FP formation (k₃) derived from small molecule assays. Additionally, the total extent of oligonucleotide mixing corresponded to the fraction of content mixing that occurred in state FP in the small molecule assay. The results show that both large "final pores" and small (presumably transient) pores can form between vesicles throughout the fusion process. The implications of this result for the mechanism of membrane fusion are discussed.

Molecular dynamics simulation analysis of membrane defects and pore propensity of hemifusion diaphragms

Membrane fusion often exhibits slow dynamics in electrophysiological experiments, involving prespike foot and fusion pore-flickering, but the structural basis of such phenomena remains unclear. Hemifusion intermediates have been implicated in the early phase of membrane fusion. To elucidate the dynamics of formation of membrane defects and pores within the hemifusion diaphragm (HD), atomistic and coarse-grained models of hemifusion intermediates were constructed using dipalmitoylphosphatidylcholine or dioleoylphosphatidylcholine membranes. The work necessary to displace a lipid molecule to the hydrophobic core of the bilayer was measured. For a lipid within the HD with radius of 4 nm, the work was ∼80 kJ/mol, similar to that in a planar bilayer. The work was much less (∼40 kJ/mol) when the HD was surrounded by a steep stalk, i.e., stalk wings forming a large angle at the junction of three bilayers. In the latter case, the lipid displacement engendered formation of a pore contacting the HD rim. The work was similarly small (40 kJ/mol) for a small HD of 1.5 nm radius, where a pore formed and grew rapidly, quickly generating a toroidal structure (<40 ns). Combining the steep stalk and the small HD decreased the work further, although quantitative analysis was difficult because the latter system was not in a stable equilibrium state. Results suggest that fine tuning of fusion dynamics requires strict control of the HD size and the angle between the expanded stalk and HD. In additional free simulations, the steep stalk facilitated widening of a preformed pore contacting the HD rim.

Vesicles and vesicle fusion: coarse-grained simulations

Biological cells are highly dynamic, and continually move material around their own volume and between their interior and exterior. Much of this transport encapsulates the material inside phospholipid vesicles that shuttle to and from, fusing with, and budding from, other membranes. A feature of vesicles that is crucial for this transport is their ability to fuse to target membranes and release their contents to the distal side. In industry, some personal care products contain vesicles to help transport reagents across the skin, and research on drug formulation shows that packaging active compounds inside vesicles delays their clearance from the blood stream. In this chapter, we survey the biological role and physicochemical properties of phospholipids, and describe progress in coarse-grained simulations of vesicles and vesicle fusion. Because coarse-grained simulations retain only those molecular details that are thought to influence the large-scale processes of interest, they act as a model embodying our current understanding. Comparing the predictions of these models with experiments reveals the importance of the retained microscopic details and also the deficiencies that can suggest missing details, thereby furthering our understanding of the complex dynamic world of vesicles.

Nanoparticle-induced permeability of lipid membranes

Monte Carlo simulations using the bond fluctuation method with explicit solvent reveal the mechanism of enhanced permeability of lipid bilayers induced by the adsorption of nanoparticles with controlled hydrophobicity. Simulation results indicate an adsorption transition of nanoparticles on the bilayer in a certain range of relative degree of hydrophobicity. In this range the nanoparticles can translocate through the bilayer, reversibly destabilizing the structure of the bilayer and inducing enhanced permeability for water and small solutes. This transition is broader for amphiphilic nanoparticles.

Distinct Initial SNARE Configurations Underlying the Diversity of Exocytosis

The dynamics of exocytosis are diverse and have been optimized for the functions of synapses and a wide variety of cell types. For example, the kinetics of exocytosis varies by more than five orders of magnitude between ultrafast exocytosis in synaptic vesicles and slow exocytosis in large dense-core vesicles. However, in all cases, exocytosis is mediated by the same fundamental mechanism, i.e., the assembly of soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins. It is often assumed that vesicles need to be docked at the plasma membrane and SNARE proteins must be preassembled before exocytosis is triggered. However, this model cannot account for the dynamics of exocytosis recently reported in synapses and other cells. For example, vesicles undergo exocytosis without prestimulus docking during tonic exocytosis of synaptic vesicles in the active zone. In addition, epithelial and hematopoietic cells utilize cAMP and kinases to trigger slow exocytosis of nondocked vesicles. In this review, we summarize the manner in which the diversity of exocytosis reflects the initial configurations of SNARE assembly, including trans-SNARE, binary-SNARE, unitary-SNARE, and cis-SNARE configurations. The initial SNARE configurations depend on the particular SNARE subtype (syntaxin, SNAP25, or VAMP), priming proteins (Munc18, Munc13, CAPS, complexin, or snapin), triggering proteins (synaptotagmins, Doc2, and various protein kinases), and the submembraneous cytomatrix, and they are the key to determining the kinetics of subsequent exocytosis. These distinct initial configurations will help us clarify the common SNARE assembly processes underlying exocytosis and membrane trafficking in eukaryotic cells.

Cell entry of enveloped viruses

Infection of cells by enveloped viruses requires merger of the viral envelope membrane with target cell membranes, resulting in the formation of fusion pores through which the viral genome is released. Since lipid membranes do not mix spontaneously, the fusion process is energy-dependent and mediated by viral envelope glycoprotein complexes. Based on their structural and mechanistic properties, three distinct classes of viral fusion proteins have been identified to date. Despite their diversity, basic principles of viral membrane fusion, simultaneous engagement of both donor and target membrane and refolding into hairpin-like structures, have emerged as universally conserved. This article provides an overview of the basic principles of viral membrane fusion common to all enveloped viruses and discusses the specific structural and functional features of the different fusion protein classes by example of the paramyxovirus, flavivirus and rhabdovirus families.

Line-tension controlled mechanism for influenza fusion

Our molecular simulations reveal that wild-type influenza fusion peptides are able to stabilize a highly fusogenic pre-fusion structure, i.e. a peptide bundle formed by four or more trans-membrane arranged fusion peptides. We rationalize that the lipid rim around such bundle has a non-vanishing rim energy (line-tension), which is essential to (i) stabilize the initial contact point between the fusing bilayers, i.e. the stalk, and (ii) drive its subsequent evolution. Such line-tension controlled fusion event does not proceed along the hypothesized standard stalk-hemifusion pathway. In modeled influenza fusion, single point mutations in the influenza fusion peptide either completely inhibit fusion (mutants G1V and W14A) or, intriguingly, specifically arrest fusion at a hemifusion state (mutant G1S). Our simulations demonstrate that, within a line-tension controlled fusion mechanism, these known point mutations either completely inhibit fusion by impairing the peptide's ability to stabilize the required peptide bundle (G1V and W14A) or stabilize a persistent bundle that leads to a kinetically trapped hemifusion state (G1S). In addition, our results further suggest that the recently discovered leaky fusion mutant G13A, which is known to facilitate a pronounced leakage of the target membrane prior to lipid mixing, reduces the membrane integrity by forming a 'super' bundle. Our simulations offer a new interpretation for a number of experimentally observed features of the fusion reaction mediated by the prototypical fusion protein, influenza hemagglutinin, and might bring new insights into mechanisms of other viral fusion reactions.

Activation thermodynamics of poly(ethylene glycol)-mediated model membrane fusion support mechanistic models of stalk and pore formation

Membrane fusion, essential to eukaryotic life, is broadly envisioned as a three-step process proceeding from contacting bilayers through two semistable, nonlamellar lipidic intermediate states to a fusion pore. Here, we introduced a new, to our knowledge, experimental approach to gain insight into the nature of the transition states between initial, intermediate, and final states. Recorded time courses of lipid-mixing, content-mixing, and content-leakage associated with fusion of 23 nm vesicles in the presence of poly(ethylene glycol) at multiple temperatures were fitted globally to a three-step sequential model to yield rate constants and thereby activation thermodynamics for each step of the process, as well as probabilities of occurrence of lipid-mixing, content-mixing, or content-leakage in each state. Experiments with membranes containing hexadecane, known to reduce interstice energy in nonlamellar structures, provided additional insight into the nature of fusion intermediates and transition states. The results support a hypothesis for the mechanism of stalk formation (step-1) that involves acyl chain protrusions into the interbilayer contact region, a hypothesis for a step-2 mechanism involving continuous interconversion of semistable nonlamellar intermediates, and a hypothesis for step-3 (pore formation) mechanism involving correlated movement of whole lipid molecules into hydrophobic spaces created by geometry mismatch between intermediate structures.

Mitofusins and the mitochondrial permeability transition: the potential downside of mitochondrial fusion

Mitofusins (Mfn-1 and Mfn-2) are transmembrane proteins that bind and hydrolyze guanosine 5'-triphosphate to bring about the merging of adjacent mitochondrial membranes. This event is necessary for mitochondrial fusion, a biological process that is critical for organelle function. The broad effects of mitochondrial fusion on cell bioenergetics have been extensively studied, whereas the local effects of mitofusin activity on the structure and integrity of the fusing mitochondrial membranes have received relatively little attention. From the study of fusogenic proteins, theoretical models, and simulations, it has been noted that the fusion of biological membranes is associated with local perturbations on the integrity of the membrane that present in the form of lipidic holes which open on the opposing bilayers. These lipidic holes represent obligate intermediates that make the fusion process thermodynamically more favorable and at the same time induce leakage to the fusing membranes. In this perspectives article we present the relevant evidence selected from a spectrum of membrane fusion/leakage models and attempt to couple this information with observations conducted with cardiac myocytes or mitochondria deficient in Mfn-1 and Mfn-2. More specifically, we argue in favor of a situation whereby mitochondrial fusion in cardiac myocytes is coupled with outer mitochondrial membrane destabilization that is opportunistically employed during the process of mitochondrial permeability transition. We hope that these insights will initiate research on this new hypothesis of mitochondrial permeability transition regulation, a poorly understood mitochondrial function with significant consequences on myocyte survival.

Energetics of stalk intermediates in membrane fusion are controlled by lipid composition

We have used X-ray diffraction on the rhombohedral phospholipid phase to reconstruct stalk structures in different pure lipids and lipid mixtures with unprecedented resolution, enabling a quantitative analysis of geometry, as well as curvature and hydration energies. Electron density isosurfaces are used to study shape and curvature properties of the bent lipid monolayers. We observe that the stalk structure is highly universal in different lipid systems. The associated curvatures change in a subtle, but systematic fashion upon changes in lipid composition. In addition, we have studied the hydration interaction prior to the transition from the lamellar to the stalk phase. The results indicate that facilitating dehydration is the key to promote stalk formation, which becomes favorable at an approximately constant interbilayer separation of 9.0 ± 0.5 Å for the investigated lipid compositions.

How SNARE molecules mediate membrane fusion: recent insights from molecular simulations

SNARE molecules are the core constituents of the protein machinery that facilitate fusion of synaptic vesicles with the presynaptic plasma membrane, resulting in the release of neurotransmitter. On a molecular level, SNARE complexes seem to play a quite versatile and involved role during all stages of fusion. In addition to merely triggering fusion by forcing the opposing membranes into close proximity, SNARE complexes are now seen to also overcome subsequent fusion barriers and to actively guide the fusion reaction up to the expansion of the fusion pore. Here, we review recent advances in the understanding of SNARE-mediated membrane fusion by molecular simulations.

VISUALIZING SOFT MATTER: MESOSCOPIC SIMULATIONS OF MEMBRANES, VESICLES AND NANOPARTICLES

Biological membranes have properties and behavior that emerge from the propagation of the molecular characteristics of their components across many scales. Artificial smart materials, such as drug delivery vehicles and nanoparticles, often rely on modifying naturally-occurring soft matter, such as polymers and lipid vesicles, so that they possess useful behavior. Mesoscopic simulations allow in silico experiments to be easily and cheaply performed on complex, soft materials requiring as input only the molecular structure of the constituents at a coarse-grained level. They can therefore act as a guide to experimenters prior to performing costly assays. Additionally, mesoscopic simulations provide the only currently feasible window on the length and time scales relevant to important biophysical processes such as vesicle fusion. We describe here recent work using Dissipative Particle Dynamics simulations to explore the structure and behavior of amphiphilic membranes, the fusion of vesicles, and the interactions between rigid nanoparticles and soft surfaces.

Cardiomyocyte deletion of mitofusin-1 leads to mitochondrial fragmentation and improves tolerance to ROS-induced mitochondrial dysfunction and cell death

Molecular studies examining the impact of mitochondrial morphology on the mammalian heart have previously focused on dynamin related protein-1 (Drp-1) and mitofusin-2 (Mfn-2), while the role of the other mitofusin isoform, Mfn-1, has remained largely unexplored. In the present study, we report the generation and initial characterization of cardiomyocyte-specific Mfn-1 knockout (Mfn-1 KO) mice. Using electron microscopic analysis, we detect a greater prevalence of small, spherical mitochondria in Mfn-1 KO hearts, indicating that the absence of Mfn-1 causes a profound shift in the mitochondrial fusion/fission balance. Nevertheless, Mfn-1 KO mice exhibit normal left-ventricular function, and isolated Mfn-1 KO heart mitochondria display a normal respiratory repertoire. Mfn-1 KO myocytes are protected from mitochondrial depolarization and exhibit improved viability when challenged with reactive oxygen species (ROS) in the form of hydrogen peroxide (H(2)O(2)). Furthermore, in vitro studies detect a blunted response of KO mitochondria to undergo peroxide-induced mitochondrial permeability transition pore opening. These data suggest that Mfn-1 deletion confers protection against ROS-induced mitochondrial dysfunction. Collectively, we suggest that mitochondrial fragmentation in myocytes is not sufficient to induce heart dysfunction or trigger cardiomyocyte death. Additionally, our data suggest that endogenous levels of Mfn-1 can attenuate myocyte viability in the face of an imminent ROS overload, an effect that could be associated with the ability of Mfn-1 to remodel the outer mitochondrial membrane.

General hydrophobic interaction potential for surfactant/lipid bilayers from direct force measurements between light-modulated bilayers

We establish and quantify correlations among the molecular structures, interaction forces, and physical processes associated with light-responsive self-assembled surfactant monolayers or bilayers at interfaces. Using the surface forces apparatus (SFA), the interaction forces between adsorbed monolayers and bilayers of an azobenzene-functionalized surfactant can be drastically and controllably altered by light-induced conversion of trans and cis molecular conformations. These reversible conformation changes affect significantly the shape of the molecules, especially in the hydrophobic region, which induces dramatic transformations of molecular packing in self-assembled structures, causing corresponding modulation of electrostatic double layer, steric hydration, and hydrophobic interactions. For bilayers, the isomerization from trans to cis exposes more hydrophobic groups, making the cis bilayers more hydrophobic, which lowers the activation energy barrier for (hemi)fusion. A quantitative and general model is derived for the interaction potential of charged bilayers that includes the electrostatic double-layer force of the Derjaguin-Landau-Verwey-Overbeek theory, attractive hydrophobic interactions, and repulsive steric-hydration forces. The model quantitatively accounts for the elastic strains, deformations, long-range forces, energy maxima, adhesion minima, as well as the instability (when it exists) as two bilayers breakthrough and (hemi)fuse. These results have several important implications, including quantitative and qualitative understanding of the hydrophobic interaction, which is furthermore shown to be a nonadditive interaction.

Diffuse interface model of multicomponent vesicle adhesion and fusion

Prefusion and postfusion states of the biological fusion process between lipid bilayer vesicle membranes are studied in this paper. Based on the Helfrich-type continuum theory, a diffuse interface model is developed which describes the phase changes on the deformable vesicles via a scalar phase field function, and incorporates the adhesion effect between the different phases of the vesicles through a nonlocal interaction potential. Various equilibrium configurations in the prefusion and postfusion states are examined. The effects of spontaneous curvatures, bending, and Gaussian rigidities on the fusion process are discussed. Instead of considering only the regions in close contact as in many previous studies, the present approach allows us to include the energetic contributions from all parts of the vesicles. By carrying out simulations based on the gradient flow of the associated energy functional, we are also able to elucidate the dynamic transitions between the prefusion and postfusion states.