| Initiation and Dynamics of Hemifusion in Lipid Bilayers Biophysical Journal, Volume 85, Issue 1, 1 July 2003, Pages 381-389 Guy Hed and S.A. Safran Abstract One approach to the understanding of fusion in cells and model membranes involves stalk formation and expansion of the hemifusion diaphragm. We predict theoretically the initiation of hemifusion by stalk expansion and the dynamics of mesoscopic hemifusion diaphragm expansion in the light of recent experiments and theory that suggested that hemifusion is driven by intramembrane tension far from the fusion zone. Our predictions include a square-root scaling of the hemifusion zone size on time as well as an estimate of the minimal tension for initiation of hemifusion. Whereas a minimal amount of pressure is evidently needed for stalk formation, it is not necessarily required for stalk expansion. The energy required for tension-induced fusion is much smaller than that required for pressure-driven fusion. Abstract | Full Text | PDF (361 kb) |
| Calculation of Free Energy Barriers to the Fusion of Small Vesicles Biophysical Journal, Volume 94, Issue 5, 1 March 2008, Pages 1699-1706 J.Y. Lee and M. Schick Abstract 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 for biological membranes, fusion involving small vesicles should proceed with little or no further input of energy. Abstract | Full Text | PDF (200 kb) |
| Field Theoretic Study of Bilayer Membrane Fusion: II. Mechanism of a Stalk-Hole Complex Biophysical Journal, Volume 90, Issue 3, 1 February 2006, Pages 915-926 K. Katsov, M. Müller and M. Schick Abstract We use self-consistent field theory to determine structural and energetic properties of intermediates and transition states involved in bilayer membrane fusion. In particular, we extend our original calculations from those of the standard hemifusion mechanism, which was studied in detail in the first article of this series, to consider a possible alternative to it. This mechanism involves non-axial stalk expansion, in contrast to the axially symmetric evolution postulated in the classical mechanism. Elongation of the initial stalk facilitates the nucleation of holes and leads to destabilization of the fusing membranes via the formation of a stalk-hole complex. We study properties of this complex in detail, and show how transient leakage during fusion, previously predicted and recently observed in experiment, should vary with lipid architecture and tension. We also show that the barrier to fusion in the alternative mechanism is lower than that of the standard mechanism by a few over most of the relevant region of system parameters, so that this alternative mechanism is a viable alternative to the standard pathway. We emphasize that any mechanism, such as this alternative one, which affects, even modestly, the line tension of a hole in a membrane, affects greatly the ability of that membrane to undergo fusion. Abstract | Full Text | PDF (436 kb) |
Copyright © 2007 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 92, Issue 8, L61-L63, 15 April 2007
doi:10.1529/biophysj.106.103341
Biophysical Letters
Avishay Efrat*, †, Leonid V. Chernomordik‡ and Michael M. Kozlov*,
, 
* Department of Physiology and Pharmacology, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
† Faculty of Sciences, Holon Institute of Technology, Holon, Israel
‡ Section on Membrane Biology, Laboratory of Cellular and Molecular Biophysics, National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland
Address reprint requests and inquiries to Michael M. Kozlov.Fusion of biological membranes boils down to merger of their bilayer lipid matrices 1,2. Currently, it is largely accepted that merger of protein-free lipid bilayers as well as fusion between viral and cell membranes, intracellular and developmental membrane fusion reactions proceed according to a universal scenario called the fusion-through-hemifusion pathway 3,4. Within this model, fusion begins with formation of a minimal lipidic bridge between the contacting monolayers of the apposing membranes referred to as the fusion stalk. Further expansion of the fusion stalk into a hemifusion diaphragm and its decay into a growing fusion pore completes the fusion process 1.
Recently, an extensive analysis of the energies of the fusion intermediates was performed aimed at revealing the rate-limiting step of the fusion reaction (see for review Chernomordik and Kozlov 1). It was shown that for lipid compositions typical for an average biological membrane, the stalk energy is only slightly larger than 40kT (where kT is the product of the Boltzmann constant and the absolute temperature), and can be overcome in the biologically relevant timescale 5. In most cases, the fusion stages downstream the stalk formation and, especially, the fusion pore growth are more energy consuming 1.
Despite its current recognition and successful applications, the fusion-through-hemifusion model leaves open important questions about the structures and energies of lipid intermediates preceding the fusion stalk formation and a possibility that these intermediates limit the rate of the whole fusion reaction. Indeed, to form a stalk, the membranes have to establish, at least locally, a dehydrated contact allowing them to perturb the continuity of their surfaces and merge, without exposure of the hydrophobic moieties of lipids to the aqueous surrounding. Such membrane contact requires overcoming the resistance of the powerful short-range repulsion forces referred to as the hydration forces 6. Straightforward estimations based on the representation of the prestalk intermediates as conventional wave-like membrane bulges approaching each other and creating extended dehydrated contacts predict the energies of hundreds of kT. Such intermediates are unfeasible and, hence, a more delicate structure must emerge on the way from two separated flat membranes to the fusion stalk.
Here we suggest an energetically feasible structure of prestalk intermediate, which will be referred to as the point-like protrusion (PLP), and is illustrated in Fig. 1.
A characteristic feature of the PLP is its sharp tip allowing for establishment of a point-like rather than extended dehydrated contact between the membranes. Obviously, this shape minimizes the energy of the hydration repulsion between the membranes. At the same time, a sharp bend of the membrane monolayers, may result in a large elastic energy making the structure unfeasible. Although aware of the latter issue, we were encouraged by the results of our previous modeling of the fusion stalk showing that interplay between the splay of the lipid hydrocarbon chains and their tilt with respect to the monolayer surface results in modest elastic energy costs of similar kinks of the monolayer profiles 7.
Below, we analyze the overall energy of the PLP (Fig. 1) using the elastic tilt-splay model for the membrane deformations 8 and the hydration force model for the intermembrane repulsion 6. The results below have to be considered as semiquantitative.
We consider a point-like protrusion formed on one membrane, while the second membrane remains flat (Fig. 1). The PLP configuration is determined by the tip angle φ, the distance dw between the flat parts of the two membranes and the requirement that no empty void forms inside the structure. The PLP monolayers undergo deformations of splay of the lipid hydrocarbon chains and tilt of the chains with respect to the monolayer surface (Fig. 1). Because we are using the essentially continuous description of a system consisting of discrete lipid molecules, the meaning of a point-like membrane contact (Fig. 1) is that of a contact with an area smaller than the area per lipid molecule a≤1nm2.
The elastic theory of the tilt and splay deformations was presented in Hamm and Kozlov 8 and applied to calculation of the energy of fusion intermediates 7. In brief, the tilt-splay energy per unit area of lipid monolayer is given by
![]() | (1) |
where the splay
and tilt
are related to the unit vector
of the average direction of the hydrocarbon chains and the unit vector
normal to the monolayer surface by
(∇ being the two-dimensional divergence operator at the monolayer surface) and
and
are the monolayer splay (bending) and tilt elastic moduli, respectively, while Js is the monolayer spontaneous curvature (spontaneous splay). The total elastic energy is given by integration of Eq. (1) over the surfaces of the two monolayers constituting the PLP.
To calculate the hydration repulsion energy of PLP formation, we have to account for the discreteness of distribution, on the membrane surface, of the effective centers generating hydration of the lipid-water interface. The characteristic distance between the hydration centers has to be of the order of dimension of a lipid polar head δ*≈1nm. Hence, the discreteness effects should be irrelevant and can be neglected if extended membrane contacts of tens of nanometers and larger are produced by smooth membrane bulges. At the same time, these effects must become significant for the point-like membrane contacts considered in this work, which are of the order of the discreteness length.
The energy of hydration repulsion per unit area of the monolayer surface is given in Rand and Parsegian 6
![]() | (2) |
To find the optimal conformation of the PLP, we calculate the sum of its elastic and hydration energies and minimize it, numerically, with respect to the PLP shape and distribution along the PLP surface of the tilt and splay of the hydrocarbon chains of the two PLP monolayers. Full details of the numerical technique we use for the energy calculation are presented in Kozlovsky and Kozlov 7.
We explored whether the sharpness of the tip of the PLP indeed decreases its overall energy, how the PLP energy depends on the initial transmembrane distance dw, and to what extent the intrinsic properties of the membrane such as its monolayer spontaneous curvature Js and the discreteness of the hydration centers δ* influence the probability of the PLP formation.
The PLP energy dependence on the tip angle φ, which determines the PLP sharpness, is presented in Fig. 2. We found that the change of φ from almost zero (smooth bulge) to 45° decreases the energy by nearly an order of magnitude (the angle range between 0° and 20° is not shown). Requirement of small tilt implied by the model (Eq. (1)) does not allow us to extend the results beyond φ=45° but, as indicated by Fig. 2, the energy continues to decrease with growing sharpness.
The PLP energy increases considerably with the distance dw (Figure 3a), is rather sensitive to the discreteness parameter δ*, and depends weakly on the monolayer spontaneous curvature Js (Figure 3b). The PLP energy depends on dw stronger, and on Js weaker than the energy of the stalk (red line). The PLP energy is smaller than that of the stalk for the experimentally relevant values of Js∼−0.1nm−1 and dw≥3nm provided that the distance between the hydration centers has a feasible value of δ*≈1nm.
We found that, in contrast to the straightforward expectations, the elastic energy of the monolayer deformation in PLP adopts modest values of the order of 10kT (not shown). Further, we obtained that, due to the shape of PLP, along with the discreteness of the hydration centers, also the hydration repulsion energy of PLP formation remains in the range of a few tens of kT.
Altogether we showed that for the relevant lipid compositions of membrane monolayers, the overall energy of PLP is lower than that of the fusion stalk, which may exceed the values presented in Fig. 3 due to contribution of the saddle-splay elasticity 9. Consequently, PLP does not limit the rate of hemifusion and the lipid dependency of hemifusion is determined by the stalk and the hemifusion diaphragm. Moreover, a point-like dehydrated contact between PLP and the target membrane facilitates stalk formation. Summarizing, PLP is a feasible prestalk intermediate completing the model pathway of fusion reaction. Experiments and/or molecular dynamics simulations are needed to confirm formation of PLP at the very initial steps of membrane fusion.
The financial support for M.M.K. by the Israel Science Foundation, the Binational United States-Israel Science Foundation, the Marie Curie Network “Flippases” and the Intramural Research Program of the National Institute of Child Health and Human Development, National Institutes of Health, is gratefully acknowledged.
1. (2003). Protein-lipid interplay in fusion and fission of biological membranes. Annu. Rev. Biochem. 72, 175–207. CrossRef | PubMed
2. (1999). Membrane fusion and exocytosis. Annu. Rev. Biochem. 68, 863–911. CrossRef | PubMed
3. (2006). Membranes of the world uniteǃ. J. Cell Biol. 175, 201–207. CrossRef | PubMed
4. (1998). Secretory and viral fusion may share mechanistic events with fusion between curved lipid bilayers. Proc. Natl. Acad. Sci. USA 95, 9274–9279. CrossRef | PubMed
5. (2001). A quantitative model for membrane fusion based on low-energy intermediates. Proc. Natl. Acad. Sci. USA 98, 7235–7240. CrossRef | PubMed
6. (1989). Hydration forces between phospholipid bilayers. Biochim. Biophys. Acta 988, 351–376. PubMed
7. (2002). Stalk model of membrane fusion: solution of energy crisis. Biophys. J. 88, 882–895. PubMed
8. (1998). Tilt model of inverted amphiphilic mesophases. Eur. Phys. J. B 6, 519–528. PubMed
9. (2004). Stalk phase formation: effects of dehydration and saddle splay modulus. Biophys. J. 87, 2508–2521. Abstract | Full Text | PDF (225 kb) | CrossRef | PubMed