Field theoretic study of bilayer membrane fusion III: membranes with leaves of different composition

Kinetic pathway and micromechanics of fusion/fission for polyelectrolyte vesicles

Despite the wide existence of vesicles in living cells as well as their important applications like drug delivery, the underlying mechanism of vesicle fusion/fission remains under debate. Classical models cannot fully explain recent observations in experiments and simulations. Here, we develop a constrained self-consistent field theory that allows tracking the shape evolution and free energy as a function of center-of-mass separation distance. Fusion and fission are described in a unified framework. Both the kinetic pathway and the mechanical response can be simultaneously captured. By taking vesicles formed by polyelectrolytes as a model system, we predict discontinuous transitions between the three morphologies: parent vesicle with a single cavity, hemifission/hemifusion, and two separated child vesicles, as a result of breaking topological isomorphism. With the increase in inter-vesicle repulsion, we observe a great reduction in the cleavage energy, indicating that vesicle fission can be achieved without hemifission, in good agreement with simulation results. The force-extension relationship elucidates typical plasticity for separating two vesicles. The super extensibility in the mechanical response of vesicle is in stark contrast to soft particles with other morphologies, such as cylinder and sphere. Our work elucidates the fundamental physical chemistry based on intrinsic topological features of vesicle fusion/fission, which provides insights into various phenomena observed in experiments and simulations.

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.

Direct observation of intermediate states in model membrane fusion

We introduce a novel assay for membrane fusion of solid supported membranes on silica beads and on coverslips. Fusion of the lipid bilayers is induced by bringing an optically trapped bead in contact with the coverslip surface while observing the bead's thermal motion with microsecond temporal and nanometer spatial resolution using a three-dimensional position detector. The probability of fusion is controlled by the membrane tension on the particle. We show that the progression of fusion can be monitored by changes in the three-dimensional position histograms of the bead and in its rate of diffusion. We were able to observe all fusion intermediates including transient fusion, formation of a stalk, hemifusion and the completion of a fusion pore. Fusion intermediates are characterized by axial but not lateral confinement of the motion of the bead and independently by the change of its rate of diffusion due to the additional drag from the stalk-like connection between the two membranes. The detailed information provided by this assay makes it ideally suited for studies of early events in pure lipid bilayer fusion or fusion assisted by fusogenic molecules.

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.

Minimum membrane bending energies of fusion pores

Membranes fuse by forming highly curved intermediates, culminating in structures described as fusion pores. These hourglass-like figures that join two fusing membranes have high bending energies, which can be estimated using continuum elasticity models. Fusion pore bending energies depend strongly on shape, and the present study developed a method for determining the shape that minimizes bending energy. This was first applied to a fusion pore modeled as a single surface and then extended to a more realistic model treating a bilayer as two monolayers. For the two-monolayer model, fusion pores were found to have metastable states with energy minima at particular values of the pore diameter and bilayer separation. Fusion pore energies were relatively insensitive to membrane thickness but highly sensitive to spontaneous curvature and membrane asymmetry. With symmetrical bilayers and monolayer spontaneous curvatures of -0.1 nm(-1) (a typical value) separated by 6 nm (closest distance determined by repulsive hydration forces), fusion pore formation required 43-65 kT. The pore radius of approximately 2.25 nm fell within the range estimated from conductance measurements. With bilayer separation >6 nm, fusion pore formation required less energy, suggesting that protein scaffolds can promote fusion by bending membranes toward one another. With nonzero spontaneous monolayer curvature, the shape that minimized the energy change during fusion pore formation differed from the shape that minimized its energy after it formed. Thus, a nascent fusion pore will relax spontaneously to a new shape, consistent with the experimentally observed expansion of nascent fusion pores during viral fusion.

Model of SNARE-mediated membrane adhesion kinetics

SNARE proteins are conserved components of the core fusion machinery driving diverse membrane adhesion and fusion processes in the cell. In many cases micron-sized membranes adhere over large areas before fusion. Reconstituted in vitro assays have helped isolate SNARE mechanisms in small membrane adhesion-fusion and are emerging as powerful tools to study large membrane systems by use of giant unilamellar vesicles (GUVs). Here we model SNARE-mediated adhesion kinetics in SNARE-reconstituted GUV-GUV or GUV-supported bilayer experiments. Adhesion involves many SNAREs whose complexation pulls apposing membranes into contact. The contact region is a tightly bound rapidly expanding patch whose growth velocity v(patch) increases with SNARE density Gamma(snare). We find three patch expansion regimes: slow, intermediate, fast. Typical experiments belong to the fast regime where v(patch) ~ (Gamma(snare)(2/3) depends on SNARE diffusivities and complexation binding constant. The model predicts growth velocities ~10 - 300 microm/s. The patch may provide a close contact region where SNAREs can trigger fusion. Extending the model to a simple description of fusion, a broad distribution of fusion times is predicted. Increasing SNARE density accelerates fusion by boosting the patch growth velocity, thereby providing more complexes to participate in fusion. This quantifies the notion of SNAREs as dual adhesion-fusion agents.

Mechanics of membrane fusion

Diverse membrane fusion reactions in biology involve close contact between two lipid bilayers, followed by the local distortion of the individual bilayers and reformation into a single, merged membrane. We consider the structures and energies of the fusion intermediates identified in experimental and theoretical work on protein-free lipid bilayers. On the basis of this analysis, we then discuss the conserved fusion-through-hemifusion pathway of merger between biological membranes and propose that the entire progression, from the close juxtaposition of membrane bilayers to the expansion of a fusion pore, is controlled by protein-generated membrane stresses.

Calculation of free energy barriers to the fusion of small vesicles

The fusion of small vesicles, either with a planar bilayer or with one another, is studied using a microscopic model in which the bilayers are composed of hexagonal- and lamellar-forming amphiphiles. The free energy of the system is obtained within the self-consistent field approximation. We find that the free energy barrier to form the initial stalk is hardly affected by the radius of the vesicle, but that the barrier to expand the hemifusion diaphragm and form a fusion pore decreases rapidly as the radius decreases. As a consequence, once the initial barrier to stalk formation is overcome, one which we estimate at 13 k(B)T for biological membranes, fusion involving small vesicles should proceed with little or no further input of energy.

Dependence of the energies of fusion on the intermembrane separation: optimal and constrained

We calculate the characteristic energies of fusion between planar bilayers as a function of the distance between them, measured from the hydrophobic/hydrophilic interface of one of the two nearest, cis, leaves to the other. The two leaves of each bilayer are of equal composition: 0.6 volume fraction of a lamellar-forming amphiphile, such as dioleoylphosphatidylcholine, and 0.4 volume fraction of a hexagonal-forming amphiphile, such as dioleoylphosphatidylethanolamine. Self-consistent field theory is employed to solve the model. We find that the largest barrier to fusion is that to create the metastable stalk. This barrier is the smallest, about 14.6k(B)T, when the bilayers are at a distance about 20% greater than the thickness of a single leaf, a distance which would correspond to between 2 and 3 nm for typical bilayers. The very size of the protein machinery which brings the membranes together can prevent them from reaching this optimum separation. For even modestly larger separations, we find a linear rate of increase of the free energy with distance between bilayers for the metastable stalk itself and for the barrier to the creation of this stalk. We estimate these rates for biological membranes to be about 7.1k(B)Tnm and 16.7 k(B)Tnm, respectively. The major contribution to this rate comes from the increased packing energy associated with the hydrophobic tails. From this we estimate, for the case of hemagglutinin, a free energy of 38k(B)T for the metastable stalk itself and a barrier to create it of 73 k(B)T. Such a large barrier would require that more than a single hemagglutinin molecule be involved in the fusion process, as is observed.

Control of membrane fusion mechanism by lipid composition: predictions from ensemble molecular dynamics

Membrane fusion is critical to biological processes such as viral infection, endocrine hormone secretion, and neurotransmission, yet the precise mechanistic details of the fusion process remain unknown. Current experimental and computational model systems approximate the complex physiological membrane environment for fusion using one or a few protein and lipid species. Here, we report results of a computational model system for fusion in which the ratio of lipid components was systematically varied, using thousands of simulations of up to a microsecond in length to predict the effects of lipid composition on both fusion kinetics and mechanism. In our simulations, increased phosphatidylcholine content in vesicles causes increased activation energies for formation of the initial stalk-like intermediate for fusion and of hemifusion intermediates, in accordance with previous continuum-mechanics theoretical treatments. We also use our large simulation dataset to quantitatively compare the mechanism by which vesicles fuse at different lipid compositions, showing a significant difference in fusion kinetics and mechanism at different compositions simulated. As physiological membranes have different compositions in the inner and outer leaflets, we examine the effect of such asymmetry, as well as the effect of membrane curvature on fusion. These predicted effects of lipid composition on fusion mechanism both underscore the way in which experimental model system construction may affect the observed mechanism of fusion and illustrate a potential mechanism for cellular regulation of the fusion process by altering membrane composition.

Membrane fusion is the transport process used for neurotransmitter release, insulin secretion, and infection by enveloped viruses. The precise mechanism of fusion is not yet understood, nor is the means by which membrane properties such as composition and curvature affect the fusion process. Here, we use molecular-dynamics simulations of lipid vesicle fusion under different lipid compositions to generate a more detailed explanation for how composition controls membrane fusion. We predict that lipid composition affects both the initial process of forming a contact “stalk” between two vesicles and the formation of a metastable “hemifused” intermediate. These two roles act in concert to change both the rate of fusion and the level of detectable fusion intermediates. We also present initial results on fusion of vesicles at different membrane curvatures. Recent experimental results suggest that the creation of highly curved membranes is important to fusion of synaptic vesicles. Our simulations cover a curvature regime similar to these experimental systems. In combination with previous results, we predict that the effect of lipid composition on fusion is general across different membrane curvatures, but that the rate of fusion is controlled by both composition and curvature.