Membrane Mechanics of Endocytosis in Cells with Turgor

Endocytosis is an essential process by which cells internalize a piece of plasma membrane and material from the outside. In cells with turgor, pressure opposes membrane deformations, and increases the amount of force that has to be generated by the endocytic machinery. To determine this force, and calculate the shape of the membrane, we used physical theory to model an elastic surface under pressure. Accurate fits of experimental profiles are obtained assuming that the coated membrane is highly rigid and preferentially curved at the endocytic site. The forces required from the actin machinery peaks at the onset of deformation, indicating that once invagination has been initiated, endocytosis is unlikely to stall before completion. Coat proteins do not lower the initiation force but may affect the process by the curvature they induce. In the presence of isotropic curvature inducers, pulling the tip of the invagination can trigger the formation of a neck at the base of the invagination. Hence direct neck constriction by actin may not be required, while its pulling role is essential. Finally, the theory shows that anisotropic curvature effectors stabilize membrane invaginations, and the loss of crescent-shaped BAR domain proteins such as Rvs167 could therefore trigger membrane scission.


Introduction
Endocytosis enables cells to internalize extracellular material and to recycle membrane components [1]. During this process, the plasma membrane is deformed into an invagination progressing inwards, which is severed and eventually released in the cytoplasm as a vesicle. Key endocytic components have been identified in several systems. This process usually involves membrane coating proteins (such as clathrin), their adaptors (such as epsins) and actin microfilaments together with associated factors [2]. We focused on the yeast model system, in which endocytosis is well characterized experimentally. The abundance and localization of the principal proteinaceous components have been measured as a function of time both in S. pombe [3,4] and S. cerevisiae [5,6]. To understand how these components work together to deform the membrane, it is necessary to consider the physical constraints under which the task is performed in vivo. While in animal cells, invaginations are opposed mostly by membrane tension and elasticity [7], turgor pressure strongly opposes invaginations in plants and fungi. In those cells, the difference of osmolarity with the outside causes a large pressure pushing membrane against the cell wall.
The magnitude of the pressure has been measured in different walled cells using various methods. For plant cells, studies converge to the range 0.2-1 MPa (see [9] for a careful review of the subject). In the yeast S. pombe, an effective pressure of 0.85 ± 0.15 MPa was derived from studying the buckling of the rod-like cells in micro fabricated chambers [10]. Based on the variation of volumes upon changes in osmolarity, the pressure in S. cerevisiae was recently estimated to be 0.6 ± 0.2 MPa [11], while an older study concluded 0.2 MPa for stationary phase cells [12]. It has been suggested that the pressure could be decreased locally by releasing osmolytes at the endocytic patch [13]. This would however practically not induce any local variation of pressure, because hydrostatic pressure gradients equilibrate at the speed of sound (possibly 1500 m/s). Given the size of yeast, this is considerably faster than an endocytic event (* 5s).
The pressure, thus of the order of 1 MPa * 1 pN/nm 2 , pushes the plasma membrane outward uniformly. This effect is balanced by an equal force from the cell wall, wherever the membrane is in contact with the cell wall (Fig 1, left). The pressure strongly opposes endocytic membrane invagination, as membrane and cell wall must come apart. The large required force is produced in yeast by an actin machinery, which forms a crosslinked network of rigid filaments around the invagination (see [8] and Fig 1,left). Actin polymerizes close to the basal plasma membrane [6]. The newly inserted F-actin at the bottom of the network is though to lift the entire network away from the cell wall (upward on Fig 1, left) [13]. For this force to be productive, the actin network should be attached to the tip of the invagination, and this is the role of the protein sla2, which is required for endocytosis-while several other proteins, including clathrin, are dispensable [5,14]. This essential role for actin is further supported by the fact that actin assembly precedes or coincides with membrane deformation [6,8]. The interplay between actin and pressure was nicely demonstrated by showing that impairing the arp2/3 complex (an actin nucleator) delayed the invagination, while decreasing the hydrostatic pressure (by adding sorbitol to the media) had the opposite effect [15].
Coat and associated proteins play a role in endocytosis, notably by recruiting the endocytic machinery, and keeping the actin network physically connected to the membrane. Proteins that bind to the lipid bilayer can also directly induce the membrane to curve [16,17]. Two familiar examples in yeast endocytosis are clathrin [1] and Rvs167 [18]. The induced curvature is expected to be qualitatively different for these two proteins, because clathrin proteins form regular triskelion [19] while Rvs167 is shaped as a crescent [20]. Moreover, clathrin adaptors can bind to membrane, cargo, actin and/or clathrin and have essential functions in endocytosis [21]. Additionally, some adaptor proteins can induce membrane curvature [22].
The shape of the membrane can be predicted by minimizing an effective deformation energy [23]. This approach has been used successfully in different systems [24][25][26], but the case of yeast endocytosis where the membrane detaches from the cell wall despite a large hydrostatic pressure has not been analyzed theoretically to our knowledge. Most previous work focused on the endocytosis of viruses in animal cells, which occurs by wrapping of the membrane around viral particles in the absence of hydrostatic pressure [7,27,28]; the constriction of vesicle necks was also studied recently [29], but without comparison to physiological results. Yeast endocytosis was modelled before, mostly without turgor [30,31]. One study considered a difference of pressure across the plasma membrane [32], but adopted a value of the pressure estimated in yeast spheroplasts [33]. The pressure in these cells, which lack a cell wall, is at least two orders of magnitude lower than under more physiological conditions [10], and the forces in this study are consequently severely underestimated.
Here, we integrate the contributions of pressure, coat proteins and membrane properties during the endocytic invagination, providing an estimate of the force that actin must exert to induce the invagination. We study the effect of coat proteins such as clathrin that induce isotropic curvature, and contrast it with crescent-shape proteins such as BAR-domain that induce curvature only in one direction. Finally, we discuss whether the combination of actin-mediated pulling together with the removal of Rvs167 is sufficient to lead to vesicle internalization.

Methods
We predict membrane shapes by minimizing a deformation energy, using a Helfrich-type Hamiltonian [24,34], in which membrane deformations are penalized by a bending rigidity κ and tension σ. We write P the difference of hydrostatic pressure between the inside and the outside of the cell. In addition, we assume a point force f a pulling the apex of the invagination, to represent the driving force generated by the actin cytoskeleton [6,13]. We thus implicitly Left: schematic representation of endocytosis. The membrane is connected mechanically to the actin network through sla2. Actin is polymerizing close to the basal membrane, thus generating an upwards (pulling) force f a . This force overcomes the rigidity κ, tension σ and the turgor pressure Π that all oppose invaginations. Right: four experimental membrane profiles at different stages, measured in an electron-microscopy study [8]. Scale bar: 20nm. assume that the forces produced by actin polymerization at the base of the invagination, and possibly by myosin motors or other processes, are transmitted to the tip of the invagination over the actin network (see Fig 1). We note κ the rigidity of the membrane together with its coat of proteins [16,17]. Moreover, we consider that the coat proteins curve the membrane either by scaffolding or inserting themselves in the membrane. We first describe proteins such as clathrin, that induce an isotropic curvature C 0 , with the same radius of curvature in both directions.
To write the membrane deformation energy, we introduce C ¼ ð 1 R 1 þ 1 R 2 Þ the local curvature, V the volume inside the invagination, S the surface area of the membrane that is not in contact with the wall, and L the height of the invagination. The total energy F of an invagination reads [26]: The rigidity term will tend to make invaginations as large as possible to minimize their curvature, while both pressure and tension will tend to make invaginations smaller, to minimize their volume and surface, respectively. Therefore an invagination dominated by pressure and rigidity will have a typical width R P ¼ ffiffiffiffiffiffiffiffiffiffiffi ffi k=2P 3 p while an invagination dominated by tension and rigidity will have a typical width l $ ffiffiffiffiffiffiffiffiffiffi ffi k=2s p . From measured values of P, σ (Table 1), and using the rigidity of a naked membrane (κ * 40 k b T), we find R P * 3 nm and λ > 9 nm. Since R P < λ in this case, we expect that invaginations will have a typical radius R P and tension will not significantly affect the shape of the invagination. In this estimate, we have used a lower bound for κ given by the rigidity of a pure lipid bilayer, but the actual value of κ will be higher because it should also include the stiffness provided by coat proteins. However, the statement R P < λ will be all the more true for higher values of κ because R P will increase slower than λ as κ becomes larger. We have used an upper bound for the membrane tension σ, but in reality it could be significantly smaller, since yeast cell have membrane furrows [35], which act as membrane reservoirs and limit the surface tension. Moreover, we will see later how fitting the experimental membrane shapes confirms that the contribution of tension is negligible.
Since the competition of pressure with membrane rigidity should control the shape of the invagination, the forces due to pressure and rigidity will have the same order of magnitude and will dominate those due to membrane tension. The total force exerted by pressure over the invagination should be PS 0 , where S 0 is the surface area of the wall that is not in contact with the membrane. This area is difficult to delineate in the electron micrographs [8], but using the Table 1. Physical characteristics of endocytosis in yeast.

Volume
V * 4 × 10 4 nm 3 [8] Surface of invagination S * 6 × 10 3 nm 2 [8] Height L up to 120 nm [8] Tip radius R t * 12 nm dimension at the tip of the invagination R t (Table 1) to obtain an ersatz pR 2 t , we estimate the magnitude of the force to * 300 pN.
Given this force, and a membrane viscosity η m * 10 −8 kg.s −1 [39], we can calculate the time needed to pull a membrane tube over a distance L * 100 nm [8]. Because the resulting time scale Lηm/f * 10 −5 s is considerably shorter than the duration of an endocytic event (a few seconds [8]), the membrane has plenty of time to reach a static equilibrium, at any stage of its evolution. In this quasi-static regime, we can derive the membrane shape equations from minimizing the energy in Eq (1), and following earlier work, we also assume that they are axisymmetric for simplicity (see supplementary information S.I. 1.1 in S1 Text).
To compare our predictions to experimental membrane profiles, the measured shapes ( Fig  1, right) were projected to an axisymmetric profile (Fig 2). For a given set of C 0 , σ, R P , only one value of f a allows the theoretical profile to match the invagination shape, and therefore three parameters could be varied to fit experimental membrane profiles (see supplementary information S.I. 1.5 in S1 Text for fitting procedure). Generally, we found that measured profiles could be fitted precisely (Fig 2) except for long invaginations (L > 80nm, see comment in the discussion). A given experimental profile can usually be fitted by a range of parameters rather than a unique set, but the values of the fitting parameters are remarkably consistent. The tension is negligible in all cases as expected, and the values of R P and κ can be interpreted and used to calculate the magnitude of the apical driving force f a . We will usually express lengths as multiples of R P and forces in terms of f P ¼ 4pPR 2 P , since this makes our results independent of the values P and κ, which are known only within an interval, without any loss of generality.
At this stage, we emphasize that the values of κ and C 0 cannot be known a priori because these parameters represent properties of the coated membrane which is a complex assembly of interacting lipids and proteins, and not a simple lipid bilayer. In the following, we determine the rigidity κ and curvature C 0 from the fit, assuming them to be constant over the membrane. In reality the coated membrane is expected to be inhomogeneous, but allowing these parameters to vary spatially would drastically increase the parameter space, such that fitting would be Predicted membrane profiles (solid lines) and rectified experimental profiles (light blue fill) measured on electron tomograms [8]. Three parameters were varied to obtain the best fit: the scale R Π * 15 -25nm, the surface tension σ and the spontaneous curvature C 0 . Usually a range of parameter values (including σ > 0) is suitable to fit each profile. The experimental membrane profiles (Fig 1, right) have been here straightened to an axisymmetric shape. Details of the method are found in the supplementary information (S.I. 1.5 in S1 Text). both underdetermined and computationally expensive. Moreover, the overall quality of the fits obtained by assuming κ and C 0 to be homogeneous does not warrant an increase in the complexity of the theory. We discuss the implications of heterogeneity in the membrane later.

Results
The coated membrane is rigid and preferentially curved The fit always provided a well determined value for R P in the range 15-25 nm, and since , this allowed us to estimate the ratio between the bending rigidity and the pressure. Assuming P * 1 MPa one gets κ * 2000 k B T, while a conservative value of the pressure (0.2 MPa) would yield κ * 400 k B T. The coated membrane is therefore much more rigid than a pure lipid bilayer (κ * 40 k B T), leading to wider invaginations. The lower estimate corresponds to the expected rigidity of clathrin coated membranes in vitro (300 k B T, [40]), while the higher value is also realistic, because other proteins besides clathrin are concentrated at the endocytic site [14]. Adaptor proteins could drastically stiffen the coat by interacting with both clathrin and the membrane [22]. By increasing the effective rigidity, coat proteins enlarge the invaginations, but as a consequence, the overall resisting force is also increased. Importantly, our estimate of the minimal force needed to sustain the invagination, f a * 3000 pN is not based on the parameters that are poorly determined by the fit (such as σ).
We were also able to determine the value of the spontaneous curvature C 0 * 0.4/R P . While this parameter is less well determined, it corresponds to a radius of curvature 2/C 0 * 100 nm, which is consistent with the known characteristics of the coat proteins. In particular, purified clathrin form in vitro spherical cages with an average radius of 35 nm [41]. The conditions in vivo are likely to be different however, since more proteins coat the membrane, and the membrane itself can also contribute to C 0 . Note that we have assumed that curvature applied everywhere to the membrane, rather than in a restricted region as can be expected in reality, but this point will be discussed later.

Large forces are required to initiate endocytosis
By solving the membrane shape equations for a range of invagination heights L, we could compute the required apical force during different stages of the invagination. While the magnitude of the force is set by f P ¼ 4pPR 2 P , its variations as a function of L are very instructive. The force has a non-zero value for L = 0, for the determined values of the spontaneous curvature C 0 * 0.4/R P (Fig 3, blue curve). f a increases and reaches a maximum for L * R P , and then decreases to a minimum value for L * 3R P . For longer invaginations a plateau is reached. This is very different from a membrane without spontaneous curvature (Fig 3, black line) or for a tension-dominated membrane [26], for which the force is minimum for L = 0. In the presence of spontaneous curvature, the initiation force at L = 0 is high, the force peaks for smaller invaginations, and the pulling force is generally lowered for longer tubes (compare black and blue lines in Fig 3).
The fact that an important initial force f 0 = f a (L = 0) is required to start the invagination is biologically meaningful. We estimated f 0 * 0.8f P , corresponding to 60% of the maximal force. Because the membrane is nearly flat if L ! 0, we could obtain analytically p (see supplementary information S.I. 2.2 in S1 Text). With σ negligible, this expression reduces to f 0 = 4πκC 0 , which summarizes how coat proteins can affect the initial force, either by changing the spontaneous curvature, or the effective rigidity of the coated membrane. By itself, spontaneous curvature is not able to lift the membrane off the wall, but rather increases the force needed to start of the invagination.
Because the force rises sharply with the height of the invagination, actin should be needed to lift the membrane already at the earliest steps of endocytosis, as recently observed [6,8]. This idea is consistent with the observation that the delay before significant membrane deformation is observed depends on the competition between actin and pressure [15]. It has been suggested that when the membrane is still flat, actin could pull at one site of the membrane while simultaneously pushing on a ring-shaped zone surrounding this site (Fig 1, left). This scenario was supported by the recent observation that sla1, an actin organizer, forms rings at endocytic sites on flat membranes, possibly indicating where pushing forces are applied [6]. Because the membrane is still juxtaposed to the cell wall at this stage, the pulling and pushing forces generated by the actin network have to be balanced.
Finally, a key feature of the force profile is that the most demanding part of the invagination occurs early during endocytosis (L < R P ). In mechanical terms, this can lead to a snap-through transition [32], in which once a critical threshold is reached, subsequent stages spontaneously follow because they do not require additional efforts. This may explain that while the duration of the endocytic early phase is highly sensitive to the competition between actin nucleation and hydrostatic pressure, later stages are largely insensitive to pressure [15]. This snap-through force profile also elegantly explains why so few retraction events were observed experimentally (less than 1% of endocytic events fail to complete [5]), despite the inherent stochasticity of a biological machine such as endocytosis, that only contains a few hundred actin filaments [3].

Spontaneous curvature leads to neck constriction
Above a certain height, invagination shapes predicted by theory exhibited a neck, even though there is no constriction force in the model (Fig 2). The only external force is a vertical lifting force applied on the apex of the invagination, and the neck appears as a consequence of the spontaneous curvature (Fig 2, last profile). Higher spontaneous curvatures will produce smaller necks, up to a point where the theory predicts a neck of zero radius, corresponding to a shape instability.
The apparition of a localized neck as a result of global spontaneous curvature has been described theoretically before [29]. We have more precisely characterized the nature of the shape instability (see supplementary information S.I. 3 in S1 Text). We can understand its biological consequences in two limit cases: with a given value of C 0 with increasing length L, or with a given length L and increasing C 0 .
Firstly, for small values of the spontaneous curvature, the membrane shape evolves smoothly upon increasing the length L of the invagination, and f a (L) is continuous. Above a threshold C Ã 0 however, we can distinguish two branches: one corresponding to short, quasitubular invaginations, and another one corresponding to larger spheroidal invaginations (Fig  4a-4c, C 0 = 0.45/R P ). The two branches overlap, i.e. for a range of lengths, there are two possible membrane shapes for a given L, as illustrated on the example profiles (Fig 4a). Above a second threshold C ÃÃ 0 , the two branches cease to overlap, such that there is a range of L where there is no equilibrium membrane shape (Fig 4d-4f, C 0 = 1/R P ). The shapes corresponding to the edges of these regions are illustrated in (Fig 4d). The transition from the short invagination to the tall spheroidal shape corresponds to a decrease in total membrane deformation energy (Fig 4f). We could compute the values of C Ã 0 and C ÃÃ 0 as a function of tension ( Supplementary  Fig. 4). The same phenomena can be observed for a given length L by increasing C 0 . There is a threshold C þ 0 ðLÞ at which the predicted radius is zero, see Fig 5. For larger spontaneous curvatures C 0 > C þ 0 ðLÞ, no stable membrane shape exist. This is in agreement with the existence of a zero-radius neck at a critical value of C 0 [29].
This shape instability is not a pearling instability (a shape instability in membrane tubes with tension [42]), because it is inhibited by membrane tension (see supplementary information S.I. 3 in S1 Text). Rather, it stems from the energetic cost of shape defects. Indeed, invaginations cannot be perfectly tubular, and the tip and base of the invagination can be viewed as defects that increase the membrane conformation energy. When the invagination takes an almost spherical shape, the tip defect is eliminated, and the cost of the base defect is minimized by making the neck infinitesimal. Indeed, we can see that the transition from tubular to spheroidal shapes reduces the total energy as the reduction in bending energy exceeds the increase of pressure-associated energy (Fig 4f). This membrane shape instability can facilitate membrane scission, a crucial step of endocytosis.

Membrane coating heterogeneity can lead to instability
Membrane coating is not homogeneous: clathrin, sla2 locate at the tip of invagination, while Rvs167 is usually located at the neck [6]. Many other proteins can associate to membranes, including the F-Bar protein Bzz1 [43], epsins such as Ent1 [44], and others [5,45]. While we do not know the exact localization, the mechanical properties and the interactions of all these proteins, it is usually believed that the clathrin-coated tip is more rigid than the base of the invagination [32]. This could also cause a higher spontaneous curvature at the tip, and additional spontaneous curvature could also stem from lipid asymmetry in the membrane leaflets.
We modeled this possible heterogeneity in rigidity as a region with higher rigidity κ (the tip), a region with a lower rigidity κ min (the base), and a small transition zone (see supplementary information S.I. 1.4 in S1 Text). Even in the absence of spontaneous curvature, there is an instability if κ min is too small compared to κ, above a threshold that depends on the surface area having higher rigidity. This instability resembles the curvature instability described above, and also stem from the destabilization of the base of the invagination (see supplementary information S.I. 1.4 in S1 Text). Overall, a membrane that is more flexible at the base is more prone to shape instability. The results that were derived for a homogeneous membrane should thus remain qualitatively valid for a heterogeneous membrane as well.

BAR-domain proteins stabilize the invagination
BAR-domain proteins are elongated crescent-like objects [46]. It was shown recently that the number of Rvs167 molecules, a BAR-domain coat protein, increases rapidly in the late stages of the endocytic process [6]. They are thought to induce anisotropic curvature in the membrane. Following previous work, we assumed that they favor curvature only orthogonally to the symmetry axis (Z), i.e. they give a favorite radius of curvature R 0 [47], rather than a favorite total curvature. Noting Γ the rigidity of a membrane coated with BAR-domain proteins, the contribution to the energy Eq (1) reads: This implicitly assumes that the coverage is uniform on the membrane, whereas in reality Rvs167 proteins are localized at the neck of the invagination [6]. This simplification was necessary however to derive the membrane shape equations (see supplementary information S.I. 1.2 in S1 Text) and sufficient to understand in essence how anisotropic curvature can affect the invagination. It allowed us to predict the membrane shape in the presence of both isotropic curvature effectors (with parameters C 0 , κ) and anisotropic curvature effectors (with parameters R 0 , Γ).
We found that anisotropic curvature inhibits the membrane shape instability. The membrane is therefore stabilized and longer invaginations can grow continuously even for high values of C 0 (Fig 6, bottom, where we used R 0 = R P for simplicity). This is in agreement with in vivo observations that invaginations without Rvs167 cannot grow longer than about 60 nm. As a corollary, removing BAR-domain proteins is a possible mean of triggering membrane scission. In our theory, this corresponds to Γ ! 0, which indeed destabilizes the membrane (Fig 6,  top). This possibility is supported by the recent observation that membrane scission is synchronous with the disappearance of Rvs167 [6].

Discussion
Using a general model for membrane mechanics, we could accurately fit experimental profiles shorter than 80nm, even though we assumed the membrane to be homogeneous, i.e. with constant rigidity and spontaneous curvature over the surface of the deformed membrane. In combination with membrane rigidity, we can expect pressure to be the dominant factor opposing membrane invagination during yeast endocytosis, while membrane tension should be negligible. This statement is derived from the dimension of the invagination, and from the scale of pressure determined experimentally, and is thus independent of the details of the model. We estimated the force required to pull the invagination based on the value of turgor pressure, and on the value of the rigidity that was determined by fitting the experimental curves. While the exact value of the force also depends on other parameters and in particular on the height of the invagination (see Fig 3), its scale is primarily determined by pressure and the width of the invagination. For the measured range of pressure P * 0.2-1 MPa, the force scale is f P * 1000-5000 pN. This is significantly larger than previously estimated (1-1000 pN [13,48,49]). The corresponding range of value for the rigidity is κ * 400-2000 k B T. This is much stiffer than a pure lipid bilayer, because the membrane is heavily coated at the endocytic site. It was suggested that phase boundaries could play a role by generating a line tension, favorable to membrane budding. However the typical line tension (of the order of * 0.4 pN [50]) is much smaller than f P , and such phenomenon can therefore be discarded.
The fact that the scale of the invagination is determined by the ratio of rigidity and pressure (R P ¼ ffiffiffiffiffiffiffiffiffiffiffi ffi k=2P 3 p ) suggests a possible tradeoff. In the absence of any reinforcement, the naked membrane could only make small invaginations with a width of 3 nm. The coat enlarges the invagination by a factor * 10 by increasing the local rigidity, but at the same time also increases the required pulling force. The mechanical properties of the coat thus likely represent an optimal tradeoff between increasing the radius of the invagination (making individual endocytic events more productive) and limiting the driving force required from the actin machinery. The ratio κ/P is likely to have been adjusted in the course of evolution, and is not expected to be conserved across species.
In general, coat proteins can induce a negative tension that can favor tubulation [51], because of their adhesion energy with the membrane. This adhesion energy has been estimated to ω * 10 −4 N/m in the case of clathrin [52], while other coat proteins may have adhesion energies of ω * 10 −3 N/m [7]. Using the largest of the two values, we find a pulling force πωR * 30pN, which is insufficient to drive the invagination against the turgor pressure. Moreover, we found that spontaneous curvature actually increases the initial pulling force in these conditions, and consequently coat proteins do not help to initially lift the membrane. This counter- Top: Theoretical membrane profiles for a membrane coated both with isotropic curvature effectors (such as clathrin) with C 0 = 1/R Π and anisotropic membrane effectors (such as bar domain proteins) with R 0 = R Π . The first three shapes correspond to three steps of invagination process with Γ = κ. The last shape was obtained with Γ ! 0, i.e. by removing the BAR-domain proteins. The removal leads to an infinitesimal neck that promote scission. Bottom: Required force as a function of invagination height, with and without anisotropic spontaneous curvature (Γ = 0 or κ), using σ = 0, C 0 = 1/R Π and R 0 = R Π . In the presence of anisotropic curvature induced by BAR-domain proteins, there is no shape instability (the blue curve is continuous). Comparing with the black curve shows that the shape instability was suppressed. intuitive result is in agreement with observations that clathrin is not necessary to initiate curvature during endocytosis in yeast [8]. This result is specific to the situation where a membrane subjected to a large pressure has to be pulled away from its supporting wall. In this configuration, the invagination must have regions with positive (at the tip) and negative (at the rim) curvature, and the energetic cost of the latter are significant in the presence of positive spontaneous curvature.
The good quality of the fits (Fig 2) indicates that the experimental invagination profiles are consistent with a pulling point force, directed away from the cell wall. In the cell, the forces generated in the actin network are probably transferred to the tip of the invagination over a finite area corresponding to the sla2 proteins. In the model, the point force is effectively propagated over the tip of the invagination, due to the inferred rigidity of the coated membrane. As long as the coat structure remains sufficiently rigid, these two situations should be equivalent.
We find that under physiological conditions of pressure, direct constriction by actin in a mechanism similar to cytokinesis is not necessary to explain the shape of the invagination, or specifically the formation of a neck. This result is similar to what was reported in the absence of turgor [29], and stems from having an isotropic spontaneous curvature that is everywhere the same on the membrane.
With isotropic spontaneous curvature, the membrane can have a shape instability (Figs 4 and 5). A more complicated model that included some level of inhomogeneity (S.I. 1.4 in S1 Text), showed that this shape instability is facilitated if the coated membrane is less rigid at the base than at the tip of the invagination (S.I. 3.2 in S1 Text). We thus expect the shape instability to be promoted by the inhomogeneity of the system, in particular by the fact that the membrane tip is covered with clathrin while the base most likely is not. Qualitatively, pressure is pushing on the neck, thus constricting it, while membrane tension is pulling on it. Thus the shape instability is promoted by pressure and inhibited by membrane tension (S.I. 3.1 in S1 Text).
Crescent-shaped curvature effectors such as Rvs167 were modelled as inducing anisotropic curvature. The analysis showed that their presence stabilizes the invagination, suggesting that Rvs167 removal can cause membrane destabilization and scission. It is known that membrane fission is promoted by amphipathic helix insertion and inhibited by crescent BAR domains [53], and a theoretical argument was made comparing idealized tubes to spherical vesicles. Here we show that crescent-shaped proteins also stabilize intermediate shapes of endocytosis, in what appears to be a very generic phenomenon.
It was reported that Rvs167-deprived cells have shorter invaginations [8]. A possible explanation is that BAR domain proteins stabilize invaginations, as in our theory. An alternative possibility is that Rvs167 would be needed to lower the pulling force required at later stages of endocytosis. However, our calculations indicate that the highest forces are required during the initial stages of endocytosis (Fig 3). Indeed experimentally, the growth speed of the later stages seems independent from the pressure [15], which is the main force opposing endocytosis. This indicates that stabilization of the shape by BAR domains, rather than reduction of force, is the most likely explanation of the observed shortening of invaginations in Rvs167 mutants.
The homogeneous theory could not fit large membrane deformations (> 80nm). Indeed, Rvs167 appears at the neck for large invaginations [6], and probably induce inhomogeneous properties near the base. This effect was not included in the inhomogeneous model described in (S.I. 1.4 in S1 Text), because introducing several different regions makes it technically intractable.
Endocytosis is ubiquitous in nature and in other model systems such as animal cells, it appears to involve a similar set of molecules. The physical conditions can however vary greatly in different cells, and this should be reflected in some of the requirements that have constrained the evolution of the endocytic machinery. In the absence of significant hydrostatic pressure, membrane coating is sufficient to generate invaginations [7], while dissipative processes involving actin and/or dynamin machinery appear necessary to membrane severing [54,55]. For the yeast system that we examined, pulling the membrane away from the cell wall seems sufficient to induce a complete budding event. The next crucial step is to decipher the mechanism by which the actin cytoskeleton is able to exert the required force. Precise quantitative information is available to do so [4,6] and theoretical work is under way [13]. Ultimately, this will offer a unified model incorporating both actin cytoskeletal dynamics and membrane mechanics, based on the experimental observations. Supporting Information S1 Text. Additional theory. 1. Membrane shape. 2. Forces. 3. Shape instability (PDF)