Biophysical Journal 87:999-1012 (2004)
© 2004 The Biophysical Society
Orientation and Interaction of Oblique Cylindrical Inclusions Embedded in a Lipid Monolayer: A Theoretical Model for Viral Fusion Peptides
Yonathan Kozlovsky *,
Joshua Zimmerberg 
and
Michael M. Kozlov *
* Department of Physiology and Pharmacology, Sackler Faculty of Medicine, Tel Aviv University, 69978 Tel Aviv, Israel; and Laboratory of Cellular and Molecular Biophysics, National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland USA
Correspondence: Address reprint requests to Michael M. Kozlov, E-mail: michk{at}post.tau.ac.il.
 |
ABSTRACT
|
|---|
We consider the elastic behavior of flat lipid monolayer embedding cylindrical inclusions oriented obliquely with respect to the monolayer plane. An oblique inclusion models a fusion peptide, a part of a specialized protein capable of inducing merger of biological membranes in the course of fundamental cellular processes. Although the crucial importance of the fusion peptides for membrane merger is well established, the molecular mechanism of their action remains unknown. This analysis is aimed at revealing mechanical deformations and stresses of lipid monolayers induced by the fusion peptides, which, potentially, can destabilize the monolayer structure and enhance membrane fusion. We calculate the deformation of a monolayer embedding a single oblique inclusion and subject to a lateral tension. We analyze the membrane-mediated interactions between two inclusions, taking into account bending of the monolayer and tilt of the hydrocarbon chains with respect to the surface normal. In contrast to a straightforward prediction that the oblique inclusions should induce tilt of the lipid chains, our analysis shows that the monolayer accommodates the oblique inclusion solely by bending. We find that the interaction between two inclusions varies nonmonotonically with the interinclusion distance and decays at large separations as square of the distance, similar to the electrostatic interaction between two electric dipoles in two dimensions. This long-range interaction is predicted to dominate the other interactions previously considered in the literature.
|
INTRODUCTION
|
|---|
A biological membrane binds multiple proteins, which play a fundamental biochemical and physiological role. A part of these proteins is totally or partially hydrophobic and, consequently, is embedded into the lipid matrix. Besides specific physiological and biochemical functions (see, e.g., Sackmann, 1995
), the integral proteins affect the mechanical properties of the membrane generating deformations and, possibly, rearrangements of the lipid matrix (Epand, 1998). In turn, the membrane deformations can mediate effective interactions between the integral proteins in the membrane plane and influence the protein conformation (Huang, 1986).
The important group of integral fusion proteins is represented by the so-called fusion peptides, which constitute fragments of large membrane proteins mediating biological membrane fusion (see, e.g., Epand, 1998; Skehel and Wiley, 2000). In the best-investigated cases of the proteins responsible for fusion of influenza and HIV viruses with the cell membranes (Skehel and Wiley, 2000), the fusion peptide is a short (
20 residues) hydrophobic N-terminal fragment of the extracellular part of the proteins. The ability of the whole protein to induce fusion is critically dependent upon the structure of the fusion peptide and its proper insertion into the target membrane (see for review Cohen and Melikyan, 2001). Moreover, isolated fusion peptides (separated from the body of the native protein) have been shown to mediate fusion of synthetic lipid bilayers (liposomes) (Davies et al., 1998; Epand et al., 1994; Epand and Epand, 1994).
Although it is generally accepted that the fusion peptides can mediate membrane fusion by inducing local structural changes in the lipid monolayer structure, the mechanism of this phenomenon remains unknown. Recently, a considerable effort has been undertaken to understand physics of lipid monolayers containing proteins (Aranda-Espinoza et al., 1996; Bohinc et al., 2003; Dan et al., 1993; Fattal and Benshaul, 1993, 1995; Fournier, 1999; Goulian, 1996; Goulian et al., 1993; Helfrich and Weikl, 2001; Netz and Pincus, 1995; Park and Lubensky, 1996; Weikl et al., 1998). In all these works the proteins were modeled as axisymmetric inclusions having a shape of cylindrical rods whose length is different from the thickness of the bilayer (Aranda-Espinoza et al., 1996; Dan et al., 1993; Fattal and Benshaul, 1993; Huang, 1986), conical (Weikl et al., 1998), or barrel-like (Fournier, 1999) molecules. At the same time, the fusion peptides insert into the membrane in conformations characterized by asymmetric distribution of the hydrophobic and hydrophilic residues along the effective molecular surface. As a result, they can adopt complicated intramembrane shapes (Han et al., 2001) and oblique orientations (Brasseur et al., 1997; Peuvot et al., 1999) characterized by the angle of 30–60° between the molecular axis and the normal to the monolayer surface.
The aim of this work is to analyze by means of the elastic model of the lipid monolayer the behavior of fusion peptides embedded obliquely into the membrane matrix, to address the type of membrane deformations induced by such inclusions, and to compute the membrane-mediated interaction between the inclusions. We show that although a priori one expects the oblique inclusion to induce a complicated monolayer deformation including perturbation of the orientation of the hydrocarbon chains, the model predicts that the deformation is, practically, limited by bending of a fragment of the monolayer surrounding the inclusion and the related rotation of the inclusion with respect to the membrane plane. The changes of chain orientation may become essential if the oblique inclusion is embedded into one monolayer of a lipid bilayer, whereas the second monolayer provides a considerable resistance to membrane bending. Analysis of this case is outside the scope of this study. The interinclusion interaction is predicted to have a character similar to the interaction of two-dimensional electric dipoles and qualitatively different from the interaction between the axisymmetric inclusions analyzed previously.
Statement of the problem
Our goal is to investigate the effects of the hydrophobic oblique inclusions on monolayer structure. The monolayer is subject to lateral tension, 
, and in the initial state preceding insertion of an inclusion it has a flat shape. The inclusion is modeled as an oblique circular cylinder, whose top base is hydrophilic and, therefore, located at the monolayer surface, whereas the side surface is hydrophobic and has to be inserted into the hydrocarbon moiety of the monolayer. The axis of the cylinder, represented by the unit vector 
is inclined by the angle 
with respect to the normal to the cylinder base 
(Fig. 1). The embedded inclusion can induce bending of the monolayer and tilt of the hydrocarbon chains of lipid molecules with respect to the normal of the monolayer surface (Hamm and Kozlov, 2000).
We first determine the conformation of the monolayer that embeds one oblique inclusion. Based on this analysis we will consider the structure of a monolayer containing two oblique inclusions, and analyze the monolayer-mediated interaction between them. Our analysis will be based on the tilt and bending theory of monolayer deformations (Hamm and Kozlov, 2000).
The embedded inclusion can exert several effects on the lipid monolayer. One can distinguish between effects that depend on the shape of the inclusion and those that are general and not shape specific. One such general effect is related to the large thermal fluctuations in the shape of lipid membranes (Helfrich, 1990). The embedded inclusions modify these fluctuations and change the free energy of the membrane, resulting in interaction between the inclusions (Netz and Pincus, 1995; Goulian et al., 1993). Another general effect arises from the lipid chains in direct contact with the inclusion. The lipid chains cannot penetrate the rigid inclusion, and, therefore, many of their chain conformations are prohibited and the entropy of the chains decreases (Fattal and Benshaul, 1993). In this work we do not consider the general effects above and focus on effects resulting from the specific oblique shape of the inclusion. To this end we find the difference between the energy of a monolayer embedding oblique circular cylinders and the energy of the reference state, which consists of a monolayer embedding right circular cylinders.
|
THEORETICAL MODEL
|
|---|
Shape and orientation of the inclusion
The inclusion is modeled as an oblique circular cylinder (Fig. 1 a). The radius of the cylinder base is a, and its height is taken to be equal to the monolayer thickness. The inclusion shape is determined by the unit vector normal to the top base and the unit vector 
directed along the cylinder axis. The angle between the two vectors, denoted by , will be called in this work the skew angle.
In the initial state the monolayer is flat (Fig. 1 b) but can be bent in the course of deformations, resulting in a change of the inclusion orientation (Fig. 1 c). To describe the system, we use a system of Cartesian coordinates, 
whose origin coincides with the center of the inclusion top base in the initial flat state (Fig. 1 b) of the membrane. The unit vectors of the Cartesian axes are denoted by 
In the initial state (Fig. 1 b) the z axis points along the initial direction of whereas the plane containing the 
and 
vectors forms the angle 
with the x axis. After the deformation (Fig. 1 c), orientation of the inclusion is described by two more angles: the polar, 
and azimuthal, ß, angles determining the new orientation of the top base normal with respect to the chosen Cartesian axes.
In the initial state, the cylinder orientation is given by
 | (1) |
After the deformation the top base normal is represented by
 | (2) |
The orientation of cylinder axis, 
in the deformed state is given by a rather complicated relationship, which can be simplified if deviates from by a small angle, 
In this case, the monolayer deformations and the related change of the inclusion orientation expressed by the polar angle-
have to be small as well, 
and the cylinder axis is given, approximately, by the equation
 | (3) |
where
 | (4) |
and that takes into account the contributions up to the second order in and 
Elastic energy of the monolayer
The first contribution to the monolayer elastic energy is related to the deformations of tilt and splay of the hydrocarbon chains of lipid molecules. To account for this energy of the monolayer we use the elastic model presented in Hamm and Kozlov (2000) and further developed in Kozlovsky and Kozlov (2002). The monolayer is described by the shape of its neutral surface lying at the interface between the polar heads and the hydrocarbon chains (Leikin et al., 1996). It is represented by the height, 
of the monolayer neutral surface above the 
plane (Fig. 1 c). The normal to the surface is given by
 | (5) |
The average orientation of the hydrocarbon chains at each point of the monolayer surface is expressed by a unit vector 
Three deformations contribute to the monolayer elastic energy. The first is tilt 
of the chain orientation 
with respect to the normal to the monolayer surface 
determined by
 | (6) |
The second and third are splay, 
and saddle splay, 
of the hydrocarbon chains. The latter deformations include additive contributions from the monolayer bending and tilt variation along the monolayer surface (Hamm and Kozlov, 2000). They can be expressed as the first and second order invariants of the tensor 
which is a covariant gradient of the chain director calculated along the monolayer surface (Appendix C). The splay is the covariant divergence of the chain director, 
whereas the saddle splay is the determinant of the director gradient, 
(Hamm and Kozlov, 2000; Kozlovsky and Kozlov, 2002). In the case of a bent monolayer with vanishing tilt, 
the splay and saddle splay reduce to the total, J, and Gaussian, K, curvatures of the monolayer surface (Hamm and Kozlov, 2000).
The structure of the monolayer is characterized by its spontaneous curvature, 
and by its saddle splay modulus, 
(Helfrich, 1973). The resistance of the monolayer to deformation is accounted by the monolayer bending, 
, and tilt, 
moduli (Hamm and Kozlov, 1998, 2000). The elastic energy per monolayer unit area related to the reference state of a flat monolayer with vanishing tilt is given by
 | (7) |
Further contribution to the elastic energy is related to the monolayer lateral tension 
Because the monolayer is undergoing bending deformation, its area, S, is larger than the projection area, A, onto the 
plane, along which acts the tension The area element, 
and the corresponding element of the projection area, 
are related by
 | (8) |
where 
is the two-dimensional gradient of the height calculated along the x-y plane.
The energy of tension related to monolayer bending is given by
 | (9) |
The total elastic energy related to reference state, which is characterized by vanishing deformations: 
is expressed for the case of small deformations by
 | (10) |
The energy of small deformations
We consider oblique inclusions whose shape deviates only slightly from a right circular cylinder, as determined by the condition The deformations of tilt, splay, and saddle splay of the lipid monolayer produced by such inclusions will be small as expressed by
 | (11) |
where 
is the monolayer thickness. These conditions are prerequisites for the validity of the elastic model of tilt and splay deformations we are using (Hamm and Kozlov, 2000). Because the monolayer deforms only slightly with respect to the flat state, we also assume that the slope of the monolayer surface is small everywhere, 
Taking into account the contributions up to the second order in 
and 
the average orientation of the lipid chains, as derived from Eqs. 5 and 6, is presented by
 | (12) |
the splay of the lipid chains is given by
 | (13) |
where 
and the elastic energy (Eq. 10) has the form
 | (14) |
where the integration is performed over the projection of the monolayer surface onto the plane, 
Note that deviation of tilt from the x-y plane contributes to higher than quadratic order terms to the energy and, hence, has to be neglected in Eq. 14. The energy of saddle splay deformation, is given in Eq. 14 by the terms proportional to 
as derived in Appendix C. It can be shown (Appendix C) that integral of this term over the membrane area vanishes so that it does not contribute either to the total energy of membrane deformation or to the energy of interaction between the inclusions (Appendix C). Hence, we omit the saddle splay contribution in the calculations below.
Boundary conditions at the inclusion boundary
The inclusion inserted into the monolayer imposes several conditions on the lipid chains adjacent to its boundary. Because the inclusion is a circular cylinder, it is convenient to use polar coordinates, 
and the corresponding unit vectors, 
We call 
the azimuthal angle to be consistent with previous sections. The origin of the coordinate system, 
is set at the center of the top base of the oblique circular cylinder. If the top base is inclined by the angle (Fig. 1 c; Eq. 2), its projection onto the plane will be an ellipse, described by
 | (15) |
but, because within the accuracy up to the linear order in the boundary between the monolayer and the inclusion in the x-y plane is a circle, 
The monolayer must match the inclusion boundary. The height of the monolayer surface at the boundary of the inclusion that changes with the orientation of the inclusion base, (Eq. 2), is
 | (16) |
The parameter 
is the vertical displacement of the center of the inclusion top base. If the top base of the inclusion is not inclined, 
all points of the top base have the same height 
The values of the parameters 
and ß will be determined from energy minimization.
The boundary condition for lipid chains adjacent to the inclusion is that the chain orientation, 
must be tangent to the inclusion surface. This condition can be expressed by
 | (17) |
Otherwise, a void is created between the inclusion surface and the lipid chains or alternatively, the lipid chains penetrate the inclusion. The explicit form of this boundary condition taking into account Eqs. 3 and 12 is given by
 | (18) |
where 
is the radial component of the tilt vector, 
Both boundary conditions (Eqs. 16 and 18) are accurate up to quadratic order in and
Euler-Lagrange equations
We have to determine the distribution of height 
and tilt 
along the monolayer surface, for which the elastic energy F adopts its minimal value. The energy variation presented in Appendix A results in the following Lagrange equations for and
 | (19a) |
 | (19b) |
whose solutions have to satisfy the conditions
 | (20a) |
 | (20b) |
 | (20c) |
The height function that satisfies these equations can be written as a sum of two parts,
 | (21) |
where 
solves Laplace equation, 
and 
solves Helmholtz equation, 
It is convenient to introduce the characteristic length of the system, 
defined by
 | (22) |
The resulting solution of Euler-Lagrange equations for the height function, whose slope decays for large 
is
 | (23) |
The terms including the modified Bessel functions 
and multiplied by the arbitrary amplitudes 
correspond to 
whereas the terms, which are multiplied by the arbitrary amplitudes 
correspond to 
The arbitrary constant angles are denoted by 
and 
and 
is a constant length. These constants along with the amplitudes will be determined from the boundary conditions. The tilt function is derived from the component of the height function,
 | (24) |
For functions that minimize the energy and satisfy the boundary conditions (Eqs. 16, 18, and 20), the expression for the energy (Eq. 14) can be simplified (Appendix A) and becomes
 | (25) |
Although the equations above are given for a single inclusion, they also apply for a monolayer containing several oblique inclusions. In this case, the boundary conditions (Eqs. 16, 18, and 20) apply independently to each inclusion and the energy (Eq. 25) consists of a sum of integrals over the boundaries of all inclusion.
|
RESULTS
|
|---|
A single inclusion: monolayer shape and the energy
In the case of one oblique inclusion embedded into the initially flat monolayer, the symmetry of the system requires that the solution is proportional to odd powers of 
Accordingly, the height (Eq. 23) and the tilt (Eq. 24) functions have the form
 | (26) |
Determination of the unknown constants 
and based on Eqs. 20a, 20b, and 20c and the boundary conditions (Eqs. 16 and 18), result in
 | (27) |
where the dimensionless parameter, 
is defined as
 | (28) |
decreases monotonically to zero as a function of growing 
or 
and approaches its maximum value, 
when 
and 
The height of the monolayer embedding the oblique inclusion is determined by
 | (29) |
The base of the inclusion is inclined by the angle 
(Eqs. 2 and 27). A profile of a monolayer embedding an oblique inclusion is shown in Fig. 2.
The tilt of the hydrocarbon chains is given by
 | (30) |
The energy of a monolayer embedding a single oblique-cylindrical inclusion relative to the energy of embedding a right-cylindrical inclusion (Eq. 25) is given in the quadratic approximation in the skew angle by
 | (31) |
The energy must be en even function of because changing 
is equivalent to rotating the inclusion orientation by 
(Eq. 3). Therefore, the next order term of F will be proportional to 
We now consider the realistic values of the parameters determining the monolayer shape, tilt, and energy (Eqs. 29–31). The tilt modulus was estimated to be 
(Hamm and Kozlov, 1998). The maximal tension a lipid bilayer can sustain before it ruptures is about 
(the exact value depends on the lipid species and the time of experiment) (Evans and Needham, 1987) meaning that the maximal monolayer tension is about 
Therefore, in realistic lipid monolayers 
so that the characteristic decay length of membrane deformations (Eq. 22) becomes 
A typical cross section radius of a rod representing such membrane inclusion as -helical peptide is about 
(Macosko et al., 1997). Using this value we obtain that 
meaning that 
Therefore, we have that both 
and 
so that can be taken at its maximal value 
It follows from the above estimations that the slope of the surface is small everywhere, 
consistent with our approximation. The tilt function is smaller by the factor 
than the height gradient, 
(Eqs. 26 and 31). Because, as shown above, this factor is small, 
the deformations of a monolayer produced by an oblique inclusion practically do not include tilt of the lipid chains. This is illustrated in Fig. 2 where we show the calculated height profile of such a monolayer embedding an oblique inclusion. Indeed, according to the figure, the monolayer contacts the inclusion boundary at a right angle, and, therefore, the lipid chains are not required to tilt. Hence, the monolayer bending is the major deformation resulting from embedding an oblique inclusion.
Interaction between two inclusions
We consider two oblique inclusions separated by a distance 
which is much larger than the inclusion radius, 
The skew angles of the inclusions are 
and 
and their orientations in the membrane plane are characterized by the azimuthal angles 
and 
which are determined with respect to the axis x connecting the centers of the inclusion bases (Fig. 3). Our goal is to calculate the energy of the membrane-mediated interaction between the inclusions, 
as a function of their separation and orientation. The calculation presented in (Appendix C) and accounting for contributions up to third order in 
and results in
 | (32) |
Simplified expressions for the interaction energy can be obtained for the two extreme interaction regimes determined by the ratio 
For large separations, 
the modified Bessel functions 
decay exponentially, so that the energy becomes
 | (33) |
In the opposite regime of 
the argument of the modified Bessel functions is small. Using the approximate form of the modified Bessel functions for small argument
 | (34) |
the interaction energy becomes
 | (35) |
The character of the interaction between the inclusions depends on the distance L between them and on their orientation with respect to the axis x connecting their centers as well as their mutual orientation, as determined by the sum of the azimuthal angles, 
It follows from Eqs. 33 and 35 that the interaction changes from repulsive to attractive or vice versa when the inclusions approach each other from long to short distances. The interaction is attractive at large and repulsive at small distances if 
In case of orientations with 
the interaction changes from repulsive at large separation to attractive at small one.
The energy change in the whole range of the interinclusion distances L is illustrated in Fig. 4 for four representative inclusion orientations. For convenience, the energy is represented in dimensionless form, 
where 
is a factor accounting for all parameters of the system, which are independent on the inclusion separation and orientation,
 | (36) |
and the distance is normalized by the characteristic length,
Consider first the inclusions, which are directed oppositely but along the x axis direction. This corresponds to the cases where their azimuthal angles are 
and 
what can be presented schematically by two arrows 
; or 
and 
(). These two cases are illustrated in (Fig. 3, b and c), respectively. All other orientations of the particles mentioned below can be presented analogously. The interaction energy shown in Fig. 4 (1) represents attraction at large L, repulsion at small L, and a shallow minimum of 
at 
The same character of interaction is predicted for two inclusions, which are directed parallel to each other but perpendicularly to the x axis (either 
or 
), what corresponds to the 
The energy profile illustrated in Fig. 4 (2), shows a minimum 
which is deeper than in the previous case, and is reached at smaller distances 
A different behavior characterized by repulsion at large separation and attraction at small ones is predicted for the cases where the inclusions are oriented parallel to each other and to the x axis, as determined 
() or 
() and illustrated in (Fig. 4 (3)); or the inclusion are antiparallel and directed perpendicularly to the x axis, as described by 
(either 
or ) and represented in (Fig. 4 (4)). In all these configurations, transition from the repulsion to attraction requires overcoming an energy barrier, which is larger in the case of antiparallel than parallel inclusions. Note that the two latter curves (Fig. 4 (3 and 4)) are mirror images of the two former ones (Fig. 4 (1 and 2)).
Summarizing the results presented in Fig. 4, the system reaches the state of the lowest energy when the two inclusions adopt the antiparallel orientation perpendicular to the x axis and approach each other to the closest distance (Fig. 4 (4)).
|
DISCUSSION
|
|---|
We addressed the elastic effects produced by insertion into a lipid monolayer of an amphiphilic inclusion modeling a fusion peptide. Our major goal was to investigate the effects resulting from obliqueness of the protein orientation within the membrane (Brasseur et al., 1997; Peuvot et al., 1999). Therefore, we modeled the peptide shape as a slanted cylinder and did not account for potentially more complicated forms such as a boomerang-like conformation, which has been observed for fusion peptide of a specific family of influenza virus hemagglutinins (Han et al., 2001). We have addressed the monolayer deformations resulting from accommodation of such inclusion in a way that its hydrophilic top base lies in the plane of the lipid polar heads, whereas the hydrophobic body resides in the hydrocarbon moiety of the monolayer. Based on this analysis, we have calculated the interaction between two oblique inclusions mediated by the monolayer deformations.
Because of a skewed shape of the inclusion, one could expect a priori that it generates deformation of tilt of the adjacent lipid hydrocarbon chains with respect to the monolayer surface and the related splay of the chains, the latter resulting from the tilt variation along the surface and from a possible monolayer bending. Our analysis based on the elastic model for tilt and splay deformation (Hamm and Kozlov, 2000) has demonstrated that for the realistic values of the monolayer elastic moduli an oblique inclusion is accommodated in the expense of the "softest mode" of deformation, which is the monolayer bending, whereas the chain tilt, practically, does not come into play. This result may have consequences for interpretation of the experimental results on orientation of inclusions, such as fusion peptides, inserted into lipid bilayers (Brasseur et al., 1997; Peuvot et al., 1999). It has to be taken into account that an apparent tilting of the peptides with respect to the membrane surface may result from local bending of the latter, while the axis of the insertion remains perpendicular to the membrane plane, as illustrated in Fig. 2.
It is important to note that our analysis does not include explicitly the effects of the second membrane monolayer, which is not penetrated by the oblique inclusion, but must undergo a bending deformation coupled to that of the inclusion containing monolayer. Bending resistance of the second monolayer increases the overall energy of the membrane curvature and may, therefore, favor some extent of the tilt deformation. According to our estimations, if the two monolayers have the value of bending rigidity of 
this effect is small and membrane bending remains the essential deformation. In a rare case where the second monolayer has a large bending rigidity or for an experimental setup where the oblique inclusion is inserted into a lipid monolayer attached to a rigid support, bending may be suppressed and the tilt of the hydrocarbon chains becomes the leading deformation.
Analysis of the interaction has shown that, depending on the orientation of the inclusions in the membrane plane, the interaction energy can have a minimum corresponding to an equilibrium interinclusion separation, or, alternatively, the inclusions tend to approach each other to zero distance after having overcome an energy barrier. The minimal final energy of the system corresponds to the latter case of the mutual approach to vanishing distance, which is accompanied by adopting the inclusions of an antiparallel orientation perpendicular to the axis connecting the centers of their bases.
An interesting feature of the interinclusion interaction is that at large separations it decays as 
(Eq. 33). This scaling of the interaction energy differs from what has been found for other membrane-mediated interactions such as that originating from the "hydrophobic mismatch", which decays exponentially with the separation L (Aranda-Espinoza et al., 1996; Dan et al., 1993), or the interaction originating from the thermal shape fluctuations of the membrane decaying as 
(Goulian et al., 1993). This unusual 
decay can be understood by considering the formal analogy between the height of the monolayer surface, and the electrical potential generated by two-dimensional dipoles. The height is the solution of the Euler-Lagrange equations (Eq. 19) and has two components (Eq. 23): satisfying the Helmholtz equation, 
and decaying exponentially at large distances, and 
which is the solution of the Laplace equation, and has a power law dependence. Hence, at large distances the component will dominate. The boundary conditions at the inclusion boundary (Eqs. 16 and 18) have polar symmetry. According to both the Laplace equation and the symmetry, determination of is similar, mathematically, to the electrostatic dipole problem. Indeed, it can be easily shown that the interaction energy of two-dimensional dipoles is given by
 | (37) |
where 
and 
are the dipole moments, and 
are the angles between the dipole moments and the axis connecting them, and L is the distance between the dipoles. The electrostatic energy (Eq. 37) has the same scaling with the distance as the energy between the oblique inclusions. Interestingly, it can be shown that this analogy is valid also for small separations, (Eq. 35).
Tension-free monolayers
In this work we have considered the monolayer subject to a lateral tension, which, in addition to determining the elastic energy, suppresses the thermal undulation of the membrane. As a result, we did not have to address the issue of the undulation-mediated interaction between the inclusions, which plays a leading role in the case of stress-free membrane (Goulian et al., 1993). Although the real lipid monolayers are in most cases subject to some tension, it is important to mention shortly the predictions of our model for the case of 
Our analysis shows that in the case of vanishing lateral tension, insertion of an oblique inclusion does not induce tilt deformation and the distribution of height of the monolayer surface generated by a single inclusion is given by
 | (38) |
where r is the distance form the center of the inclusion top base and all other notations are the same as defined above.
The interaction energy between two oblique inclusions reduces to
 | (39) |
The corresponding Euler-Lagrange equations have been derived and solved in Weikl et al. (1998). Calculation of the interaction energy, by the method described above leads to
 | (40) |
The interaction is always repulsive and decays with separation L much faster than the interaction in the presence of lateral tension (Eq. 33). Because of its fast decay 
the interaction (Eq. 40) can be neglected in comparison to the interaction mediated by the thermal membrane undulation, which has been calculated for two right circular cylinders (Goulian et al., 1993; Helfrich and Weikl, 2001; Park and Lubensky, 1996) and shown to scale as 
|
CONCLUSION
|
|---|
It is commonly believed, that fusion proteins promote the fusion reaction by inserting into the target membrane and reducing the energy of fusion intermediates such as membrane stalks, hemifusion diaphragms, and fusion pores. Indeed, such action of the fusion proteins can be envisaged if one assumes that, in accord to their orientation inside the membranes, a fusion protein has an effective shape of oblique cylinder that promotes tilt of the hydrocarbon chains of lipid molecules. In this case, insertion of the fusion peptides into the fusion stalks could relieve the tilt of the chains, and fill the packing defect in the middle of the stalk and in the rim of the hemifusion diaphragm (Kozlovsky et al., 2002; Kozlovsky and Kozlov, 2002). However, the analysis above shows that an oblique cylindrical inclusion inserted into a lipid monolayer results in an anisotropic bending of the monolayer surface around the inclusion, but does not lead to tilt of hydrocarbon chains. According to our estimates, this monolayer bending cannot be the sole factor determining the membrane rearrangements. Therefore, the fusion peptides must have another property allowing them to drive the fusion reaction. This conclusion is consistent with data on mutations in the fusion peptide region that indicate a more specific role in fusion of the exact chemical side group of each amino acid in the fusion peptide conservative substitutions leading to changes in fusion phenotypes (Cross et al., 2001; Qiao et al., 1999). Possibly, the role of the fusion peptide is to create connections between the fusion proteins leading to formation of protein coats, which promote transformation of the early fusion intermediates into fusion pores (Chernomordik and Kozlov, 2003; Kozlov and Chernomordik, 2002). In addition, the membrane shape generated by a fusion peptide (Fig. 2) is reminiscent of the border between liquid ordered and disordered lipids domains. Perhaps the fusion peptide is designed as an oblique inclusion to reduce the line tension of lipid microdomains asymmetrically (in the outer cellular and viral leaflets only), and, hence, to stabilize the small size of microdomain needed for fusion.
|
Appendix A
|
|---|
In this section we derive the Euler-Lagrange equations and the related boundary conditions, solve these equations, and find a simple expression for the elastic energy of a monolayer.
Euler-Lagrange equations
The energy of a nearly flat monolayer with bending and tilt deformations is Eq. 14:
 | (A1) |
To derive the Euler-Lagrange equations, we consider small variations in the height of the monolayer surface, 
and the tilt of the hydrocarbon chains, 
The variation in the energy resulting from the height variation, 
is given by:
We manipulate the first term of the integrand
and obtain
The second term of the integrand can be presented in the form:
leading to
The first term of the integrand can be transformed into a boundary integral by the divergence theorem.
The variation in the energy resulting from the tilt variation, 
is given by:
We present the first term of the integrand in the form
and obtain
The first integral can be transformed into a boundary integral by the divergence theorem. The variation in the energy to linear order when both the surface height and the tilt are varied is the sum of the two variations:
 | (A2) |
where the line element of the boundary is 
and the normal to the boundary is 
The condition for obtaining the energy minimum is that the first variation vanishes, 
The variations in the functions, 
and 
are arbitrary and independent. From the vanishing of the surface integral we obtain the Euler-Lagrange equations:
|
| (A3a) |
 | (A3b) |
We must also demand that the boundary integral vanishes (Eq. A2). We now specifically consider an oblique inclusion in a monolayer. The monolayer has one boundary around the inclusion and a second boundary far from the inclusion (at infinity). Because the deformations caused by the oblique inclusion are localized, at the far boundary they have to decay to zero so that the boundary integral vanishes. We are left with the boundary integral around the inclusion. The inclusion boundary is a circle. The boundary element is 
and the normal to the boundary is 
The boundary integral is set to be zero
 | (A4) |
The variations at the boundary, 
and 
are not independent because they must satisfy the boundary conditions (Eqs. 16 and 18).
Solution of the Euler-Lagrange equations
The divergence of Eq. A3b can be presented as
its insertion into Eq. A3a results in
and elimination 
from Eq. A3a leads to
 | (A5) |
We define the characteristic length of the system
so that the Eq. A5 adopts the form
 | (A6) |
One solution of the equation satisfies
This is Laplace equation in two dimensions, whose solution in polar coordinates is
where and 
are arbitrary constant angles. However, because we are interested only in surfaces with decaying slope, described by the asymptotic boundary condition 
we retain only the corresponding terms
where is a constant. To obtain the second solution of Eq. A6, we denote 
and obtain
Solution of this equation can be presented as the sum of two parts,
where is the solution of Laplace equation, and is the decaying solution of Helmholtz equation, 
having the form
where are the modified Bessel functions. The tilt is derived from h to be
Concluding, the solution of Euler-Lagrange equations is
For h and satisfying the Euler-Lagrange equations, the energy integral can be simplified. We present the energy (Eq. A1) as
and perform the same derivation we used to calculate the first variation of the energy. As a result we obtain
where the surface integral vanishes and only the boundary integral contributes to the energy. The latter can be further simplified to
 | (A7) |
Boundary conditions
The condition for vanishing of the boundary integral part of the variation (Eq. A4) can be written in a similar form as Eq. A7:
 | (A8) |
The boundary condition (Eq. 16) contains three free parameters 
and 
They are found from the requirement that given an independent arbitrary variation of theses parameters, 
and 
the energy variation (Eq. A7) vanishes. The possible variation of the boundary condition (Eq. 16) is
Instead, we can present it as
where 
and 
are linear combinations of 
and 
The possible variation of the second boundary condition (Eq. 18) is
Inserting the two last expressions into the condition (Eq. A7), we obtain,
Because the three variations 
and 
are independent, each of the three integrals must vanish, and we obtain the three conditions
Including these conditions and the boundary conditions (Eqs. 16 and 18), the energy (Eq. A7) can be further simplified
|
Appendix B
|
|---|
Interaction between two oblique inclusions
We first define for each inclusion a coordinate system whose origin is at the center of its top base. The polar coordinates are denoted 
and 
for the first and second inclusion, respectively. In both coordinate systems, the azimuthal angle is measured with respect to the direction from the first inclusion toward the second one. Thus, this direction is given by 
and also by 
The two oblique inclusions are slanted by the angles and called the skew angles, in the orientation of the azimuthal angles 
and 
respectively (Fig. 3).
We denote by 
the solution for the height of the monolayer surface, h, determined by a single isolated inclusion (Eq. 29). The height resulting from two distant inclusions can be presented as a sum of the solutions for two isolated inclusions, 
with a small correction, 
Each correction 
is given by the full series (Eq. 23) and the total height is:
 | (B1) |
The Laplace operator is invariant upon translation, 
TOP
ABSTRACT
INTRODUCTION
In this figure the orientation of the oblique cylinder is along the x axis, which means that 
(c) The monolayer undergoes bending deformation. 
(ß is not shown).
and 
(3) A side view of two oppositely oriented inclusions with 
for several orientations of the inclusions: (1) 
or 
; (2) 
; (3) 