This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Shemesh, T. Articles by Kozlov, M. M. Search for Related Content PubMed PubMed Citation Articles by Shemesh, T. Articles by Kozlov, M. M.

Prefission Constriction of Golgi Tubular Carriers Driven by Local Lipid Metabolism: A Theoretical Model

Tom Shemesh ^*, Alberto Luini , Vivek Malhotra , Koert N. J. Burger  and Michael M. Kozlov ^*

^* Department of Physiology and Pharmacology, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel; Department of Cell Biology and Oncology, Istituto di Ricerche Farmacologiche Mario Negri, Santa Maria Imbaro (Chieti) Italy; Cell and Developmental Biology Department, Division of Biology, University of California at San Diego, La Jolla, California; and Department of Molecular Cell Biology, Institute of Biomembranes, Utrecht University, Utrecht, The Netherlands

Correspondence: Address reprint requests to Michael M. Kozlov, Department of Physiology and Pharmacology, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv 69978, Israel. Tel.: 972-3-640-7863; Fax: 972-3-640-9113; E-mail: michk{at}post.tau.ac.il.

	ABSTRACT

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
Membrane transport within mammalian cells is mediated by small vesicular as well as large pleiomorphic transport carriers (TCs). A major step in the formation of TCs is the creation and subsequent narrowing of a membrane neck connecting the emerging carrier with the initial membrane. In the case of small vesicular TCs, neck formation may be directly induced by the coat proteins that cover the emerging vesicle. However, the mechanism underlying the creation and narrowing of a membrane neck in the generation of large TCs remains unknown. We present a theoretical model for neck formation based on the elastic model of membranes. Our calculations suggest a lipid-driven mechanism with a central role for diacylglycerol (DAG). The model is applied to a well-characterized in vitro system that reconstitutes TC formation from the Golgi complex, namely the pearling and fission of Golgi tubules induced by CtBP/BARS, a protein that catalyzes the conversion of lysophosphatidic acid into phosphatidic acid. In view of the importance of a PA-DAG cycle in the formation of Golgi TCs, we assume that the newly formed phosphatidic acid undergoes rapid dephosphorylation into DAG. DAG possesses a unique molecular shape characterized by an extremely large negative spontaneous curvature, and it redistributes rapidly between the membrane monolayers and along the membrane surface. Coupling between local membrane curvature and local lipid composition results, by mutual enhancement, in constrictions of the tubule into membrane necks, and a related inhomogeneous lateral partitioning of DAG. Our theoretical model predicts the exact dimensions of the constrictions observed in the pearling Golgi tubules. Moreover, the model is able to explain membrane neck formation by physiologically relevant mole fractions of DAG.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
Transport of proteins between intracellular organelles occurs by means of membrane containers, also called transports carriers (TC) (Bonifacino and Lippincott-Schwartz, 2003; Rothman and Wieland, 1996). Schematically, one can distinguish between two structural types of TCs. One is represented by the small round vesicles of 50–100 nm in diameter initially covered by protein coats. Carriers of the second type are the large pleiomorphic, mostly tubular-saccular containers which mediate a substantial part of ER-Golgi and Golgi-plasma membrane transport in mammalian cells (Polishchuk et al., 2000). These are of variable size but often very large (up to 10 µm in length), travel along microtubules, and change their shape elastically during movement (Hirschberg et al., 1998; Polishchuk et al., 2000). For convenience, we will refer to the former type of transport carriers as vesicular TCs and to the second type as large tubular TCs.

A crucial step in TC generation is shaping of the membrane into a bud connected to the initial membrane by a narrow neck. The subsequent fission of the neck results in a separate membrane container. Neck formation is an energy-consuming process because it requires strong membrane bending which is opposed by the bilayer bending rigidity (Helfrich, 1990). What is the cell machinery supplying the energy for neck formation and constriction? Obviously, the molecular workers responsible for this membrane remodeling events are specific proteins and/or lipids.

It is commonly agreed that the major driving force for the formation of vesicular TCs comes from the coat-forming protein complexes—such as the clathrin-adaptor complex assisted by dynamin and its partners, and the coatomer complexes COPI and COPII (Kirchhausen, 2000)—which have been identified, characterized, and demonstrated to bend the lipid bilayers and generate small vesicles in vitro (Matsuoka et al., 1998; Spang et al., 1998; Takei et al., 1998).

In contrast and despite the crucial role large TCs play in intracellular membrane traffic, the mechanism of their formation remains largely unknown. Lipids such as polyphosphoinositides, phosphatidic acid (PA), and diacylglycerol (DAG) have long been recognized as crucial factors in the regulation of Golgi-plasma membrane transport mediated by large TCs (Baron and Malhotra, 2002; De Matteis et al., 2002; Kearns et al., 1998, 1997). However, they are commonly believed to play, primarily, a signaling role, whereas the origin of the force deforming the membranes and generating the carriers remained to be identified.

The goal of the present study is to address theoretically the mechanism by which membrane necks form in the Golgi tubules involved in Golgi-plasma membrane transport. Our analysis indicates that in this particular case the lipids can serve as the molecular workers, performing the mechanical work necessary for membrane remodeling.

Useful experimental insight into this phenomenon came recently from observations on Golgi membrane constriction and fission made in isolated Golgi fractions incubated with CtBP/BARS (Weigert et al., 1999). CtBP/BARS plays an essential role in the formation of large TCs that carry proteins and lipids from the Golgi to the plasma membrane (A. Luini, unpublished data). In vitro, in the presence of long-chain acyl-CoAs, CtBP/BARS induces local constrictions in Golgi tubules resulting in narrow necks. Occasionally, the narrow necks were separated by regular distances along the tubule resulting in a shape similar to a pearl chain; in analogy to Bar-Ziv et al. (1999), this will be referred to as membrane pearling. In some cases pearling was followed by fission of the necks and membrane separation into TCs. The pearling phenomenon has been related to the lipid metabolic activity of CtBP/BARS. CtBP/BARS is a lipid acyltransferase, which, in the presence of acyl-CoAs, acylates lysophosphatidic acid (LPA) converting it into phosphatidic acid (PA) in the cytoplasmic leaflet of the Golgi tubule. The generation of membrane necks and the pearling of Golgi tubules observed in vitro is likely to be driven by the same mechanism as membrane neck formation and constriction in the course of detachment of large transport carriers from the Golgi complex. Therefore, elaboration of the pearling mechanism should have important implications for TC formation in live cells.

We propose a novel model for membrane pearling based on the assumption that the PA generated by CtBP/BARS undergoes enzymatic conversion into DAG, a transformation known to proceed rapidly in vivo (Nanjundan and Possmayer, 2003). DAG molecules are characterized by a unique molecular shape (Szule et al., 2002), which can strongly influence the shape of the whole membrane. The essence of our model lies in the coupling between local shape of the tubule and local concentration of the newly synthesized DAG. We account quantitatively for the documented features of pearling of the Golgi tubules (Weigert et al., 1999) and provide experimentally testable predictions on the partitioning of the newly formed DAG into the emerging membrane necks. Our model strongly suggests that physiological levels of DAG may drive the constriction of Golgi membrane tubules, a key step in the formation of transport carriers mediating Golgi-plasma membrane transport.

	PHENOMENOLOGICAL BACKGROUND

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
A lipid membrane is commonly described as a surface, whose shape is determined by the total, J, and Gaussian, K, curvatures (Appendix A). For the sake of brevity, by referring to J we will skip the word "total" and call it simply the curvature. The membrane bilayer is represented by its midplane with the curvature, J_mid, whereas the shapes of the outer and inner membrane monolayers are attributed to their neutral planes (Kozlov and Winterhalter, 1991) whose curvatures are denoted as J_out and J_in, respectively. Because of the opposite orientation of the two monolayers, their total curvatures, J_out, and J_in, have opposite signs. The exact definitions of the curvatures and the relationships between them are presented in Appendix A.

If a membrane monolayer is not subject to any external forces, it adopts a spontaneous shape, which is determined by the monolayer intrinsic structure and characterized by the spontaneous curvature, J_s (Helfrich, 1973). These concepts can be extended to the description of lipid molecules, whose influence on the membrane form is attributed to their effective molecular shape (Israelachvili et al., 1976) and the effective molecular spontaneous curvature, (Kozlov and Helfrich, 1992). A more detailed discussion of these notions is presented in Appendix A. The lipid molecules tend to reside in the membrane regions where the curvature, J, is as close as possible to the molecular spontaneous curvature, . This tendency is opposed by the effects of entropy of lipid mixing in the monolayers (see below).

Geometry of a Golgi tubule
Initially, a typical tubule has a shape of a nearly smooth cylinder connected at least at one end to a disklike cistern which should allow for a free exchange of the inner aqueous volume of the tubule with that of the much larger cistern (Weigert et al., 1999).

Shape transformation of a tubule upon the action of CtBP/BARS results in the emergence of local constriction sites separated by bulges (see Weigert et al., 1999; and this article, Fig. 1, b–d). The average tubular diameter in the constrictions is 11 nm, whereas in the bulges it reaches 48 nm. If present on the same tubule, the constriction sites are separated by regular intervals of 85 nm, which results in a pearl-chain-like shape.

View larger version (76K): [in this window] [in a new window]   FIGURE 1  Membrane shapes. (a) A family of Delaunay surfaces, ranging from a flat cylinder to a deformed state, all having the same total curvature. (b) Superposition of a Delaunay surface outline with an image of a deformed state of the Golgi tubule. c and d represent, respectively, the intermediate and the final states of pearling, where the sum of Rmin and Rmax have the same value. (e) The red line represents the theoretically derived contour of the constriction site. b–d have been published in Weigert et al. (1999); the scale bar in c and d represents 40 nm.

A DAG molecule has a very pronounced conelike shape (Appendix A). Its spontaneous curvatures varies from to -1 nm^-1 depending on the length and the degree of unsaturation of the hydrocarbon chain (Szule et al., 2002). In fact, DAG is the lipid with the most negative spontaneous curvature determined so far. The spontaneous curvature of PA has not been measured quantitatively but its effective shape has been determined by nuclear magnetic resonance studies (Kooijman et al., 2003). At physiological pH and in mixture with other lipids, PA exhibits a cone shape with moderately negative spontaneous curvature, ^PA < 0, which is close to that of dioleoyl-phosphatidylethanolamine (Fuller and Rand, 2001). Under the same conditions LPA has an effective shape of inverted cone as expected for a molecule having one hydrophobic chain. The spontaneous curvature of oleoyl-LPA is positive, ^LPA > 0, and close to that of oleoyl-lysophosphatidylcholine (LPC) (Kooijman et al., 2003). According to preliminary data it constitutes (E. Kooijman and R. P. Rand, unpublished data). Summarizing, the conversion of LPA into DAG via PA results in transformation of the effective shape of the molecule from inverted cone to conelike, and thus to a change of the molecular spontaneous curvature from positive to strongly negative values.

	QUALITATIVE MODEL

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
We suggest the following mechanism for pearling of a Golgi tubule. In the initial cylindrical tubule the outer and the inner lipid monolayers have, respectively, a positive, J_out > 0, and a negative, J_in < 0, curvature (Appendix A). The DAG molecules generated in the outer monolayer as a result of LPA PA DAG conversion have a strongly negative spontaneous curvature, ^DAG < 0, which conflicts with J_out. As a result, the DAG molecules are "frustrated" and the outer monolayer accumulates the elastic stresses, which the system tends to relax. There are three routes the membrane goes to relieve the stresses:

1. Repartitioning of DAG from the outer to the inner membrane monolayer whose curvature, J_in < 0, matches in sign ^DAG.
   
2. Deviation of the tubular shape from the smooth cylinder.
   
3. Redistribution of DAG along the surfaces of the two monolayers.
   

Flipping of DAG to the inner monolayer decreases the stress, but does not relieve it entirely, for two reasons. First, a part of DAG molecules is retained in the outer monolayer by the entropic effects. Second, the curvature of the inner monolayer, -0.1 nm^-1 < J_in < 0, is not negative enough to completely accommodate the spontaneous curvature of DAG, -1 nm^-1.

To relax the stress further, the membrane develops periodic constrictions and bulges, thus adopting the pearl-chain-like shape (Fig. 1, b–d). In the constriction sites the bilayer has a necklike shape, where the monolayer curvatures J_in and J_out are shifted toward more negative values as compared to those in the bulges. Redistribution of DAG into the constriction sites stabilizes the necks and the overall energy decreases.

Thus, the morphological changes of the membrane tubules result from a coupling between local membrane composition and local shape (Leibler and Andelman, 1987; Markin, 1981).

In the following we present a detailed treatment of this scenario based on the elastic model of lipid monolayers.

	PHYSICAL MODEL

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
To analyze the pearling of the Golgi tubule, we have to compute the energy of different tubular membrane shapes and different distributions of the DAG molecules between the two monolayers and along the monolayer surfaces. For that we need to know the constraints imposed on the possible membrane configurations and have a model for the membrane free energy.

Constraints on the membrane shapes and lipid distributions
The initial cylindrical shape of the tubular membrane is characterized by a difference between the areas of the outer, A_out, and inner, A_in, membrane monolayers, A = A_out - A_in, which results from a global relationship between the monolayer areas of the whole complex including a Golgi cistern and the connected tubules. In the course of the pearling transition, A has to remain constant. This condition can be expressed by (Safran, 1994)

(1)

where the integration is performed over the bilayer midplane and the monolayer thickness (or, more accurately, the distance from the monolayer neutral surface—see Appendix A—to the bilayer midplane) is assumed to be equal for the two monolayers and constant along the membrane plane.

The local mole fraction of the newly synthesized molecules in a membrane monolayer will be denoted by . Although these molecules can repartition between the two membrane monolayers and along the monolayer planes, their average mole fraction in the membrane remains constant and equal to the _tot produced by the metabolic reaction. This condition leads to the second constraint expressed by

(2)

Monolayer free energy
For simplicity, we assume that the monolayer consists of two lipid components. The major one is characterized by the spontaneous curvature ₀ and represents the average of the lipids constituting the membrane before the onset of the lipid transformation. We will refer to ₀ as the background spontaneous curvature. The second minor component has the spontaneous curvature and corresponds to the newly formed molecules.

The free energy of a monolayer consists of the elastic energy of bending and the contribution of the entropy of lipid mixing.

The bending energy per unit area of a monolayer, f_B, is given by the Helfrich model (Helfrich, 1973) extended in Mitov (1978) as presented in Appendix B. For our study, only two energy contributions will be relevant,

(3)

where J and K are, respectively, the total and the Gaussian curvature of the monolayer, and J_s is its spontaneous curvature (Appendix A). The monolayer bending modulus, , has a value of 4 x 10^-20 J (Niggemann et al., 1995). The quadratic Gaussian modulus, , has not been directly measured. Its estimation, based on the model of trans-membrane lateral stress profile, shows that is negative and has an order of magnitude of , where ₀ 50 mN/m (Goetz and Helfrich, 1996; and Appendix B). Assuming 1.2 nm, the estimation gives .

We will assume the elastic moduli and to be independent of the membrane composition, as supported for by the experimental results (Leikin et al., 1996). In contrast, the monolayer spontaneous curvature, J_s, is often linear with respect to the mole fractions of membrane components (Kozlov and Helfrich, 1992; Leikin et al., 1996), and is given by

(4)

where = - ₀.

The contribution to the monolayer free energy from the entropy of lipid mixing can be presented in the mean field approximation as in Kozlov and Helfrich (1992) and Andelman et al. (1994),

(5)

where a is the area per lipid molecule, assumed to be equal for the different lipids constituting the monolayer. In the following, the subscripts in and out indicate the values describing the inner and the outer membrane monolayers, respectively, whereas the subscript mid indicating the bilayer midsurface is omitted.

The main equations and the outline of analysis
The total free energy of the membrane is the sum of the elastic and entropic contributions of the two monolayers,

(6)

To proceed, we insert Eqs. 3–5 in Eq. 6, and express all curvatures and the area elements through those of the midsurface of the bilayer using the relationships presented in Appendix A. The resulting final form of Eq. 6 is

(7)

where the integration is performed over the area of the midsurface.

We will minimize the energy (Eq. 7) with respect to the distributions of curvatures, J and K, and the DAG mole fractions, _in and _out, along the membrane surface. The minimization will take into account the constraints (Eqs. 1–2). The total DAG mole fraction, _tot, will serve as a major parameter controlling the shape of the tubule.

	RESULTS

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
Shapes of the tubules: Delaunay surfaces
The tubule shapes have to satisfy the constraint (Eq. 1), which requires a constant difference between the monolayer areas, A. This condition restricts only the integral of the total curvature, whereas the local values of J may vary along the membrane surface. Despite this possibility, we further restrict ourselves by considering the shapes of constant J along the membrane surface.

The axial-symmetric shapes of constant total curvature, J = const, are called Delaunay surfaces (Delaunay, 1841). A family of such shapes, all having the same J, is illustrated in Fig. 1 a. The surfaces have a periodic shape and each of them is characterized by two radii, R_max and R_min, in the widest and the narrowest cross-section, respectively. The constant total curvature of a Delaunay surface is

(8)

The length of periodicity of a Delaunay surface, which can be measured as a distance between two closest constrictions, is

(9)

A cylinder is a Delaunay surface with R_min = R_max = R_cyl. According to Eq. 9, all the Delaunay surfaces having the same curvature J are characterized by the same sum of the maximal and minimal cross-section radii,

(10)

The shapes of the Golgi tubules observed in Weigert et al. (1999) are well described by the Delaunay surfaces. Indeed, the average radii measured for the tubular constrictions, R_min 6 nm, and bulges, R_max 24 nm, require, according to Eq. 10, the periodicity length of L 85 nm. This is in excellent agreement with the experimental value. We also directly superimpose an outline of a Delaunay surface, determined by the parameters above, with a representative experimental image (Fig. 1 e). Bearing in mind that the image is a two-dimensional projection of a three-dimensional shape, which may not be oriented exactly in the imaging plane, the theoretical prediction is in accord with the observed shape. Finally, we verify whether the relationship (Eq. 10) is satisfied for a typical intermediate pearled shape (Fig. 1 c) and one that corresponds to the almost vanishing radii of constrictions (Fig. 1 d). In both cases, R_min + R_max has the same value of 30 nm.

The agreement between the Delaunay surfaces and the shapes of the Golgi tubules confirms the correctness of our assumption of constancy of the total curvature J. To further verify this assumption we have analyzed numerically deviations of the tubular shape from the Delaunay surface and the related variations of J along the membrane. Minimization of the bending energy, which changes in this case in the second order of J, confirmed that the shape of the lowest bending energy is that of constant J.

Although the midplane of the bilayer described by a Delaunay surface has a constant curvature, J, its Gaussian curvature, K, changes from positive values, K > 0, in the bulges to negative values, K < 0, in the constriction sites. This results in variation along the membrane surface of the total curvatures of the two monolayers, J_out and J_in, which differ from the midplane curvature J (Appendix A). In the constriction sites they become more negative than in the bulge, which stimulates repartitioning of the DAG molecules into the constriction sites.

It has to be emphasized that consideration of the shapes of membrane tubules belonging to the Delaunay family of surfaces requires accounting for the energy contributions up to the fourth order in the principal curvatures. Indeed, it can be readily seen that, due to constancy of J, all lower order energy contributions either remain constant or their integrals vanish according to Gauss-Bonnett theorem.

Partitioning of DAG
Finding the optimal distribution of DAG in the two monolayers, subject to the constraint of Eq. 3, has to be performed numerically. However, it is instructive to consider first a simplified case where the DAG mole fractions in the two monolayers are set equal, _in = _out = . For a cylindrical shape, DAG is distributed homogeneously over the membrane monolayers, = _tot. Development of constrictions is accompanied by the emergence of the Gaussian curvature, K, of the tubular surface. The analytical solution can be found for expansion of in orders of K assuming that the total DAG mole fraction is small, _tot << 1. The details of the computation are presented in Appendix C. The resulting distribution of DAG accounting for the contributions up to the second order in K is

(11)

where K² is the surface average of the Gaussian curvature squared. According to Eq. 11 the DAG molecules, whose spontaneous curvature is more negative than that of the background, < 0, are concentrated in the constrictions where K < 0, as illustrated in Fig. 2.

View larger version (18K): [in this window] [in a new window]   FIGURE 2  The optimal distribution of the DAG mole fraction along the length of the membrane tubule. (a) The numeric results for (1) the inner monolayer, (2) the outer monolayer, (3) initial even distribution in the inner monolayer, and (4) initial even distribution in the outer monolayer. (b) The outline of the Delaunay surface and the schematic illustration of DAG repartitioning. Note that this illustration describes the excess DAG in the constricted region.

As a result of the molecular rearrangement within each monolayer, the DAG level in the constriction sites is 30–50% (depending on the monolayer) higher than the constant absolute value of the mole fraction in the initial cylinder, although throughout the remainder of the shape the mole fraction is very close to the initial constant value (Fig. 2). The latter is because the area of the constriction sites is small compared to the overall surface area so that even when the high DAG mole fraction at those parts is taken into account, the number of DAG molecules there remains small and has little effect on the mole fraction elsewhere in the membrane.

Criterion for pearling
The total energy of the tubule, F_tot, is obtained by integration of Eq. 7 over the membrane area while accounting for Eq. 11. In general, for the self-consistency of the elastic model, a term proportional to the Gaussian curvature and several terms of higher order in the curvatures have to be added to Eq. 3, as presented in the Appendix B. However, because of the constancy of the total curvature J of the shapes we consider, and due to the Gauss-Bonnet theorem (see Helfrich, 1990), all these terms do not contribute to the energy of the tubule constriction. For the case of symmetric trans-monolayer distribution of DAG we obtain for the energy

(12)

where A_tot is the total area of the tube, and

(13)

The parameter characterizing the membrane shape and distinguishing between different members of the Delaunay family of surfaces is the average of the Gaussian curvature squared K². For a cylinder the Gaussian curvature vanishes, K² = 0. The tighter the constrictions of a Delaunay surface are, the larger is K².

The energetically preferable shape of the tubule is determined by the sign of the coefficient B (Eq. 13), which depends on the total mole fraction of DAG, _tot. If B > 0, the growth of K² results in an increase of F_tot. In this case, to have a minimal energy, the system keeps the cylindrical shape with vanishing Gaussian curvature, K = 0, and, according to Eq. 11, the DAG molecules are distributed homogeneously along the membrane surface.

In the opposite case where B < 0, the emerging Gaussian curvature, K² > 0, results in decrease of the energy, F_tot < 0. This means that the tubule deviates from the cylindrical shape and undergoes pearling. This is accompanied by repartitioning of the DAG molecules into the constrictions of the pearled tubule (Eq. 11 and Fig. 2).

Pearling starts when the DAG mole fraction _tot exceeds a critical value,

(14)

Based on Eq. 13, this critical value is

(15)

Note that Eq. 15 is meaningful only for the parameter range where adopts positive values. Based on Eqs. 12 and 13, if the values of the parameters result, formally, according to Eq. 15, in zero or negative values of , the cylindrical shape is predicted to be unstable intrinsically, i.e., to undergo pearling also in a case of homogeneous membrane.

All parameters determining the critical DAG mole fraction (Eq. 15) are known from the experiments, except for the quadratic Gaussian modulus . For the latter we have only an order-of-magnitude theoretical estimation (Goetz and Helfrich, 1996; and Appendix B). However, we can fit the value of to the experimental value of . According to the assumption of our model, DAG results from dephosphorylation of PA. Therefore, the critical DAG mole fraction, , must have a value similar to the experimentally found PA mole fraction of 0.25% needed to trigger constriction and pearling of the Golgi tubules (Weigert et al., 1999). Inserting , a = 0.6 nm², = 1.2 nm, and = 4 x 10^-20 J into Eq. 15, we find that = 0.0025 if = -1.12 x 10^-37 J x m². This value is very close to the previous estimation, -10^-37 J x m² (Goetz and Helfrich 1996), meaning that the model is able to explain triggering of the pearling of the Golgi tubules by the experimentally predicted mole fractions of DAG.

To determine the critical DAG concentration, , for an unrestricted DAG redistribution between the monolayers and for arbitrary values of _tot, we have performed numeric analysis, the results of which are illustrated in Fig. 3. The numeric and analytical results for the dependence of the coefficient B on the total DAG mole fraction, _tot, nearly coincide in the range of _tot << 1. The obtained numeric results for , which correspond to B(_tot) = 0, as indicated in Fig. 3 b, are very close to those given by Eq. 15 for small _tot.

View larger version (18K): [in this window] [in a new window]   FIGURE 3  A plot of the coefficient B (13) determining the pearling onset as a function of DAG mole fraction. (a) Numerically computed value for B for the case of free redistribution of DAG between the monolayers (solid line) and analytically derived B for assumed equal partitioning of DAG between the monolayers (dashed line). (b) The same plot, centered and magnified around the critical point where B changes sign. The critical concentration beyond which the flat cylinder loses stability and undergoes pearling (B < 0) is indicated by an asterisk.

	DISCUSSION

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
We present a mechanism by which DAG molecules resulting from LPA acylation to PA and subsequent dephosphorylation induce pearling of a Golgi tubule, which is a development of periodic constrictions and bulges in the initially cylindrical shape. Within a constriction site the membrane has the shape of a narrow neck, identical to that expected to form in vivo at an intermediate stage of transport carrier formation. Fission of such a neck results in carrier detachment. We suggest that this mechanism underlies the generation of the large pleiomorphic carriers that mediate transport from the Golgi complex to the plasma membrane.

The force driving neck formation is generated as a consequence of local lipid metabolism. It is an elastic force, which results from transformation of LPA molecules characterized by a positive spontaneous curvature into DAG molecules which have a strongly negative spontaneous curvature. The newly formed DAG redistributes between the two membrane monolayers and tends to partition into the emerging membrane necks, thus enhancing formation of the latter.

We show that pearling has a phase-transition character. It occurs abruptly when the total mole fraction of the LPA molecules converted into DAG exceeds a critical value, , estimated to be a couple of tenths of mole percent of total lipid.

Assumptions and limitations of the model
We want to emphasize that the major quantitative result of the model, the critical mole fraction of the newly synthesized phospholipid necessary to trigger pearling, , is very sensitive to the value of the quadratic Gaussian modulus, . At the same time, this modulus has never been measured and its value can be only roughly estimated based on the intermembrane stress profile (see Goetz and Helfrich, 1996, and Appendix B). Changes of , which may be related to variations of lipid compositions between membrane tubules or even along the same tubule, can change considerably the critical molar ratio resulting in pearling. Furthermore, according to Eq. 15, the numerical results for also depend explicitly on the background spontaneous curvature, ₀, which can only be estimated. Furthermore, the background lipids do not all have the same spontaneous curvature and can undergo repartitioning, which is, de facto, insufficient to drive pearling, but can contribute to the effect of the newly synthesized DAG. Based on all these limitations, the results of our model have a qualitative rather than quantitative character, and although providing a physical mechanism of the pearling phenomenon, they cannot give an exact prediction of the critical mole fraction . Due to the uncertainty in the values of and ₀, we expect that Golgi tubule constriction in vivo may occur at higher or even lower mole fractions of DAG than that estimated in the present study.

Our model implies that the newly synthesized molecules are evenly distributed over the tubular membrane, and that the mechanical properties of the tubule, such as the elastic moduli of its membrane and deformability of its lumen, are homogeneous. However, a Golgi tubule in vivo is certainly inhomogeneous in terms of both its lipid composition and elastic properties. The metabolic reactions produce, most probably, locally elevated concentrations of the corresponding lipids; the tubular lumen may be partially filled by proteins limiting the ability of the tubule to undergo constriction; and finally, the membrane elastic moduli can change along the membrane being influenced by proteins bound to the membrane surface. In such a realistic case, the results of the model can be applied locally to a limited region of the tubule where DAG is synthesized and the lumen properties allow for membrane deformation into a neck. Accordingly, formation of just one or a few constrictions rather than pearling of the whole tubule would be expected if the local concentration of DAG reaches the critical DAG mole fraction .

Is the LPA-DAG a unique metabolic reaction resulting in pearling?
We have discussed a specific case of membrane pearling, which results from LPA PA DAG transformation catalyzed by CtBP/BARS and a PA-phosphatase. However, the model and the resulting criteria, Eqs. 14 and 15, can be generalized and account for membrane neck formation driven by any metabolic lipid transformation. The major parameter determining the ability of a specific lipid conversion reaction to induce membrane neck formation is the difference between the spontaneous curvature of the newly formed lipid, , and the background spontaneous curvature, ₀.

The value of ₀ can be estimated based on the average lipid composition of the Golgi complex (van Meer 1998) and data available on the spontaneous curvatures of individual lipids (Chen and Rand, 1997; Leikin et al., 1996; Szule et al., 2002, and references therein). The major lipid components of the Golgi complex are phosphatidylcholines (PC) and phosphatidylethanolamines (PE). The spontaneous curvature of unsaturated PC such as dioleoyl-phosphatidylcholine (DOPC) is close to _DOPC -0.1 nm^-1, whereas of saturated PCs is less negative (Sjolund et al., 1987). The spontaneous curvature of unsaturated PE such as dioleoyl-phosphatidylethanolamine (DOPE) is _DOPE -0.3 nm^-1 (Kozlov et al., 1994; Rand and Fuller, 1994). The spontaneous curvature of the saturated PEs has so far not been measured, but is expected to have values more negative than those of PCs due to a smaller size of the hydrated headgroup (Rand and Parsegian, 1989). For estimations we assume the background spontaneous curvature to be ₀ = -0.1 nm^-1.

The molar fraction of the newly synthesized lipid, , necessary to induce pearling is inversely proportional to ()² (see Eq. 15). This means that the larger the deviation of the spontaneous curvature of the newly synthesized lipid from the background value ₀, the more effective this lipid is in inducing pearling.

Another consequence of this relationship is that not only lipids such as DAG, whose spontaneous curvature is more negative, < 0, but also lipids with more positive than the background value, > 0, may induce neck formation. An example of the latter reaction may be de-acylation of lipids into lysolipids such as lysophosphatidylcholine (LPC) characterized by a positive spontaneous curvature of _LPC 0.26 nm^-1 (oleoyl-LPC); see Fuller and Rand (2001). The difference is, however, that the newly formed molecules with < 0 concentrate in the neck region where the Gaussian curvature is negative, K < 0, whereas molecules with > 0 tend to repartition into the bulges with K > 0.

The relationship between and can be expressed in the form of a phase diagram computed according to Eq. 15 and shown in Fig. 4 for = 4.14 x 10^-20 J, a = 0.6 nm², = 1.2 nm, and = -1.18 x 10^-37 J x m². For reasons of clarity, we present it in terms of , where is the absolute value of the difference in the spontaneous curvatures. The bell-shaped line in the (, R_c)-plane, called the phase boundary, separates ranges of the parameters that result in a smooth cylindrical shape (outside the "bell" ), from those that result in the pearled (inside the "bell") shapes. The dotted line (a) in Fig. 4 clearly illustrates that if the product of a metabolic reaction results in a lipid whose spontaneous curvature is not sufficiently different from the background level (), such a reaction cannot induce pearling, not even at a very high mole fraction of the synthesized lipid.

View larger version (13K): [in this window] [in a new window]   FIGURE 4  Phase diagram of the tubules. The y axis represents the mole fraction of the nonbilayer lipid and the x axis represents , where is the difference between the spontaneous curvature of this lipid and that of the background. Region I corresponds to the smooth cylindrical shape. Region II corresponds to the pearled shape. The numerically computed phase boundary is indicated by the solid line. The vertical lines represent (a) a lipid unable to induce pearling, and (b) a lipid capable of pearling the tubule.

The critical mole fractions, , for several lipids, based on the assumption that the background spontaneous curvature is ₀ = -0.1 nm^-1, are presented in Table 1, which also includes the predictions of the model for the case where the synthesized lipid remains in the outer monolayer of the Golgi tubule (i.e., does not undergo transbilayer movement). It can be seen that only DAG, whose spontaneous curvature deviates most from the background value, can induce pearling at feasible mole fractions. The for other lipids are larger than those expected in cell membranes under physiological conditions. Specifically, assuming that the spontaneous curvature of PA is close to that of DOPE (Kooijman et al., 2003), its membrane mole fraction necessary for pearling to occur has to reach 3%, which is >10x higher than the average PA amount in isolated Golgi membranes incubated with CtBP/BARS and acylCoA (Weigert et al., 1999). In addition, although the transbilayer movement of PA has not been studied in Golgi membranes, it is likely to be slower than that of DAG. If the newly synthesized PA remains in the outer monolayer its mole fraction necessary for pearling is predicted to be as high as 9%. This result confirms our assumption that lipid-induced pearling requires conversion of PA into DAG.

Sources and roles for DAG in the Golgi complex: a blue-collar worker?
The conclusion that DAG is the most probable lipid candidate for driving neck formation in a Golgi tubule is supported by the general phenomenology collected on this lipid in studies of the Golgi complex.

DAG is constitutively synthesized in the Golgi complex. Formation of PA from LPA through CtBP/BARS is only one of the possible metabolic routes leading to DAG. Another source of PA in the Golgi is the removal of the headgroup from glycerolipids, especially PC, by PLD (Exton, 2002); PA can then be converted very quickly to DAG by PA phosphatases (Nanjundan and Possmayer, 2003). Furthermore, a potentially efficient DAG generating machinery independent of PA is the enzyme sphingomyelin synthase, which transfers the headgroup of PC to ceramide, producing sphingomyelin and DAG in the lumenal leaflet of the Golgi (Fang et al., 1998). This is probably the main DAG synthetic route in yeast (Fang et al., 1998). Finally, DAG can be generated by PLC acting on PIs, which can be synthesized, and are abundant, in the Golgi (De Matteis et al., 2002). DAG can also be consumed by metabolic activities present in the Golgi. The important ones are the CDP choline pathway which uses DAG to make PC in yeast (Fang et al., 1998) and the PC-synthase pathway in mammalian cells (Henneberry et al., 2002). In addition, DAG can be very rapidly converted into PA by DAG kinases (Nanjundan and Possmayer, 2003), of which there are numerous isoforms in mammalian cells; or it can be deacylated by DAG lipases (Brindley et al., 2002). Thus, although the Golgi DAG levels resulting from this complex metabolic network are difficult to predict, there is little doubt that local DAG synthesis may result in DAG concentrations predicted to induce tubule narrowing and pearling (this article).

Not only is DAG present in the Golgi, but it is also required for protein transport from the trans-Golgi network to the plasma membrane, both in yeast, as indicated by extensive genetic analysis (Kearns et al., 1997), and in mammalian cells, as shown more recently by biochemical experiments (Baron and Malhotra, 2002). The possible models of DAG action have been recently reviewed (Huijbregts et al., 2000; Kearns et al., 1998). There are two main non-exclusive possibilities. One, so far favored in the literature, is based on the signaling function of DAG. Indeed, DAG has been shown to bind and recruit to the trans-Golgi network protein kinase D (Baron and Malhotra, 2002), which regulates the fission of the transport carriers destined to the plasma membrane (Liljedahl et al., 2001; Van Lint et al., 2002), and the recruitment and activation of ARF-GAP1, an important regulator of vesicle formation and fission (Yang et al., 2002). These experiments also indicate that the likely execution point of DAG in the formation of carriers is tubule fission. Finally, a signaling role of DAG is suggested by genetic evidence in yeast (Yanagisawa et al., 2002).

The alternative model suggests that DAG molecules, having a small and electrically neutral polar head, can redistribute into the internal membrane monolayer and modify the physical properties of membrane, stimulating in this way membrane budding (Kearns et al., 1998). In this case, the DAG molecules do not serve as messengers but rather as "blue-collar workers," directly carrying out the mechanical work. The model we present belongs to the latter kind of mechanisms and suggests a specific mode of DAG-mediated membrane shaping.

Alternative models for pearling
Besides the suggested model, two more conventional mechanisms of pearling need to be discussed briefly.

An external pulling force applied to the tubule and producing tension in its membrane may result in pearling (Bar-Ziv and Moses, 1994; Markin et al., 1999). However, there is no source for such a pulling force in the Golgi in vitro. Indeed, the tubules, as seen in the images of Weigert et al. (1999), are not stretched and, consequently, cannot be under tension. Moreover, in contrast to the Golgi tubules (Fig. 1), the experimentally observed shapes resulting from the tension-induced pearling differ considerably from the Delaunay surfaces (Goldstein et al., 1996) providing additional evidence that pearling of the Golgi tubules is not driven by membrane tension.

Another possible mechanism is based on the bilayer-couple effects or increase of the overall bilayer spontaneous curvature at a constant volume-to-area ratio of the tubule (Tsafrir et al., 2001). Pearling may be driven by an increase in the monolayer area difference, A, resulting from acylation of the LPA molecules in the outer monolayer of the tubule and the related increase of their molecular area. This mechanism requires that the internal volume of the tubule does not change. This condition, most probably, does not hold because the tubules are connected to the Golgi cisternae. Nevertheless, let us estimate the change in A necessary to produce the degree of pearling observed experimentally. To do so we consider a flat cylinder as the initial state and a row of connected spheres as the final state of the membrane. Assuming that both the membrane area and the enclosed volume are conserved we relate the cylinder radius, R_cyl, and the sphere radius, R_sphere, by

(18)

which is equivalent to the relationship between the curvatures

(19)

Inserting Eq. 19 in Eq. 2, and accounting for our more stringent requirement that the curvature remains locally constant, J_mid = const, we obtain that the difference of the monolayer areas in the spheres is related to that in the cylinder by A_sphere 1.3 x A_cyl. Hence, pearling of the tubule driven by the bilayer-couple effect requires an increase of the monolayer area difference of 30%. How large should the mole fraction of LPA converted to PA be, to provide such increase in A? The molecular area of PA should be 45% larger than that of LPA (Demel et al., 1992). We then obtain that LPA has to constitute 90% of the lipids of the initial membrane. This is much larger than the mole fraction of LPA in the Golgi membranes, which typically does not exceed 0.6% (Weigert et al., 1999). Hence, the bilayer-couple mechanism cannot explain pearling of Golgi tubules induced by CtBP/BARS.

	PHYSIOLOGICAL SIGNIFICANCE AND CONCLUSIONS

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
Formation of transport carriers involves membrane budding and subsequent fission of the membrane neck connecting the forming carrier with the initial membrane. A necessary step in this process is narrowing of the neck connecting the bud with the initial membrane. Our model describes a lipid-driven mechanism for this process taking, as a paradigm, the pearling of Golgi tubules induced by the combined action of CtBP/BARS generating PA, and a PA-phosphatase, converting the PA into DAG.

Is the DAG-based mechanism of tubule pearling described here a valid model for in vivo transformations of Golgi tubules? Tubule constrictions are technically difficult to observe in thin sections of cells prepared for EM but, when suitable techniques (e.g., electron tomography and scanning EM) have been employed, they have been seen (Lindsey and Ellisman, 1985). DAG, as noted, is both present in the Golgi and needed for carrier fission in vivo (Baron and Malhotra, 2002). Thus, together, these observations indicate that the formation of constrictions regularly distributed along Golgi tubules as observed in cells may well be mediated by a local rise in DAG level.

Finally, what is the relationship between the membrane neck formation we describe and the next step, fission of the membrane neck (Kozlovsky and Kozlov, 2003)? Do the DAG molecules also promote, in addition to pearling, severing of the membrane necks, or, alternatively, are further lipid transformations or an involvement of proteins such as PKD (Baron and Malhotra, 2002) necessary to complete fission and carrier formation? These questions cannot be answered by the present model and motivate future studies.

	APPENDIX A

TOP ABSTRACT INTRODUCTION PHENOMENOLOGICAL BACKGROUND QUALITATIVE MODEL PHYSICAL MODEL RESULTS DISCUSSION PHYSIOLOGICAL SIGNIFICANCE AND... APPENDIX A APPENDIX B APPENDIX C ACKNOWLEDGEMENTS REFERENCES

 
Geometrical description of membranes
Membrane curvature and spontaneous curvature
Generally, surface geometry is determined at each point by two principal curvatures, c₁ and c₂ (Fig. A1). For description of membranes (Helfrich, 1973), one uses the combinations, J = c₁ + c₂, called the total curvature, and K = c₁ x c₂, referred to as the Gaussian curvature, the two having a profound geometrical meaning (Spivak, 1970).

View larger version (10K): [in this window] [in a new window]   FIGURE A1  An illustration of an area element with the principal curvatures c1 and c2 indicated.

View larger version (26K): [in this window] [in a new window]   FIGURE A2  An illustration of a lipid bilayer structure and individual lipid shapes. (a) A representation of a bilayer midsurface and the neutral surfaces of its constituting monolayers. (b–d) The effective shapes of lipids with zero, positive, and negative spontaneous curvature, respectively.