article: We present a Master Equation approach to calculating polymerization dynamics and force generation by branched actin networks at membranes. The method treats the time evolution of the F-actin distribution in three dimensions, with branching included as a directional spreading term. It is validated by comparison with stochastic simulations of force generation by actin polymerization at obstacles coated with actin “nucleation promoting factors” (NPFs). The method is then used to treat the dynamics of actin polymerization and force generation during endocytosis in yeast, using a model in which NPFs form a ring around the endocytic site, centered by a spot of molecules attaching the actin network strongly to the membrane. We find that a spontaneous actin filament nucleation mechanism is required for adequate forces to drive the process, that partial inhibition of branching and polymerization lead to different characteristic responses, and that a limited range of polymerization-rate values provide effective invagination and obtain correct predictions for the effects of mutations in the active regions of the NPFs.
  Forces exerted by polymerization of monomeric actin (G-actin) into filamentous actin (F-actin) are crucial for bending the cell membrane in many important cellular processes, including cytokinesis, cell migration, and, under some conditions, endocytosis [1]. Specifically, actin is required for yeast clathrin-mediated endocytosis (CME), a central mechanism that controls cellular signaling, nutrient uptake and membrane recycling [2]. CME is driven by a transient protein patch, in which different proteins appear in a well-defined sequence [3], including actin and its nucleators. The actin patch bends a small portion of the cell membrane into a highly curved invagination that encloses extracellular substances. The invagination is later severed and its contents, as well as lipids and membrane proteins, are released into the cytoplasm. For this membrane-bending process, the actin network needs to exert both pulling forces and pushing forces (see Fig 1). The required pushing forces are several pN per filament [4,5], mainly to overcome the large (∼0. 2MPa) osmotic pressure difference [6] (turgor pressure) across the membrane, because of the small number (∼102) [7] of actin filaments at each endocytic site. The machinery driving CME constitutes a coupled mechanochemical network [8]. Force regulates protein dynamics via processes such as the slowing of actin polymerization by opposing force; conversely polymerization of actin and assembly of curvature-generating proteins generate force. We are only beginning to understand the basic properties of this network. Protein dynamics during CME have been extensively studied via fluorescence imaging methods. Assembly of endocytic proteins (EPs), including F-actin, was first quantified [3,9] using relative fluorescence intensities. Later, Ref. [10] developed a method for measurement of the absolute counts of the EPs in fission yeast (Schizosaccharomyces pombe). In Ref. [7], and later in Ref. [11], absolute counts were measured in budding yeast (Saccharomyces cerevisiae). These studies have suggested a count of about 6000 polymerized-actin subunits at the endocytic site [2], with the counts of other proteins typically in the range of 50 to 300. Actin nucleators, or “NPFs”, precede actin polymerization, which proceeds over a period of about 15 seconds. These quantitative measurements have inspired several quantitative modeling studies of dynamics of the EPs [4,8, 12,13]. The mechanical aspects of CME are less well understood due to the difficulties of measuring forces on a scale of tens of nanometers in vivo. Balance of forces on the actin network requires that inward pulling forces at the center of the endocytic site are opposed by equal pushing forces from the outer regions of the actin network [4]. The mechanics of bending the cell membrane in CME were studied in detail in a recent model [5] based on the “Helfrich” free energy density [14]. The authors calculated the energy-minimizing shape functions of the membrane during endocytosis, using parameters fitted to electron microscopy tomography data [15]. However, the dynamics of the actin force were not included when obtaining the shape functions, so it was not possible to calculate a time-dependent shape nor to include the mechanochemical feedbacks driving the protein and shape dynamics. In Ref. [4], one of us treated actin as an actively growing gel simulated using a finite element method (FEM), and thus predicted endocytic invagination dynamics. However, the actin growth was modeled with a simple phenomenological description, which was not quantitatively compared to experiment. It is thus unclear how well the growing gel represents the actual dendritic structure of the network [16]. The mechanics of stationary CME membrane profiles were also treated in Refs. [17] and [18], but again neither treated the protein dynamics in the process. A mechanochemical model of CME was proposed in Ref. [8]. It contained several types of feedback interactions, which indirectly impacted actin polymerization. This model explained several traits of endocytic mutants. However, the treatment of actin polymerization was highly simplified, and the model used an extremely low value of the turgor pressure. A major challenge in developing a complete description of the mechanochemical network driving CME is to accurately model the actin network and its interaction with the cell membrane. In Ref. [12], we and others developed a stochastic model of actin polymerization during CME. The F-actin in the model was modeled via a stochastic-growth method that gave an explicit three-dimensional actin network, with parameters fitted to experimental data. The force opposing actin polymerization was assumed to “kick in” when the network reached a certain size. This work revealed some important feedback mechanisms between actin and its nucleators, required for CME. But the membrane mechanics were oversimplifed by using a step-function force opposing actin polymerization, and the membrane profile was not obtained explicitly. Additionally, this stochastic model required a large ensemble of repeated calculations to produce meaningful results. It was thus difficult to quantitatively fit the model, for use in other potential applications. Furthermore, stochastic models are difficult to use in studying possible oscillatory behaviors and bifurcations of the actin network that require the recognition of subtle changes in the F-actin count. Deterministic rate-equation approaches, as in Ref. [19] and parts of Ref. [12], would thus be more convenient for experimental fitting and capturing subtle effects. However, such rate equations cannot treat the mechanics and geometry of the actin network because they describe only the average F-actin count rather than the spatial distribution of F-actin. A number of deterministic reaction-diffusion approaches have improved on rate-equation approaches by modeling the spatial distribution of F-actin explicitly, using various assumptions about dynamics and the interaction of F-actin with the cell membrane [20–28]. However, such methods have not explicitly included the oblique branched geometry of the network, and their force-generation component has not been validated by comparison to stochastic-simulation results. Ref. [29] developed a two-dimensional treatment of the actin network from a spatially dependent rate equation, explicitly treating branching angles. The mesh size was small in comparison with the size of the cell, but still large enough to treat the coarse-grained density of F-actin. This method was shown to give promising results for global cell properties such as migration. However, it is not clear to what extent it can be applied to processes such as CME, which are fully three-dimensional and have significant structure at very small length scales. In this paper, we propose a Master Equation (ME) method to treat the reciprocal interactions of polymerizing actin and its nucleators with a bending membrane, and apply it to CME. The ME method describes the spatial distribution of F-actin using a single simulation for a given set of parameter values, while having nearly the realism of the stochastic-growth approach implemented with a large numbers of runs. The methodology explicitly includes the branching geometry, and is validated by comparison with stochastic simulations for the case of an actin network pushing an obstacle. It uses a mesh size smaller than the characteristic size of the actin network, to calculate a probability distribution function (pdf) of actin subunits at given points in time and space. It builds on the work of Ref. [29] by treating a three-dimensional geometry, using a smaller mesh size (about 2 nm vs. 100 nm) that allows better treatment of actin-based forces, and calculating the pdf rather than the coarse-grained actin density. It differs from the reaction-diffusion approaches above in its more complete description of both the orientation and length of new branches in three dimensions, in its more accurate treatment of force generation, and in the use of a pdf. We apply the ME method to a mechanochemical model of CME in budding yeast that treats the time courses of F-actin, its nucleator Las17, and the deformation of the membrane. The model integrates the chemical variables F-actin and Las17 (slow), and the membrane shape variables (fast) into one dynamically interacting system, and shows how the actin network bends the cell membrane in real time. The model accomplishes several important goals: 1) a theory of dendritic actin polymerization that is mechanistically realistic, numerically accurate and computationally efficient, 2) a mechanochemical model of the dynamics of the cell membrane driven by the actin network during CME, and 3) a more accurate model of CME that quantitatively determines several core parameters that were “floating” (not determined by experimental data) in the previous model [12]. The results of the ME model are consistent with experimental data [7,12] for protein dynamics and the effects of mutations. New predictions from the model include the following: i) a spontaneous nucleation mechanism is required in the central portion of the endocytic site, ii) controlled inhibition of branching and/or polymerization lead to characteristic behaviors of the peak counts of actin and its main nucleator, and iii) a certain range of polymerization rates is required for robust invagination, and correct prediction of the peak counts of actin and its main nucleator in mutants.
  Our model treats dendritic actin network growth in the presence of capping, with new filaments created as branches induced by a planar distribution of nucleation-promoting factors (NPFs), as on a membrane or hard substrate. The model describes polymerization in three dimensions, but we introduce it in two dimensions first, for the sake of clarity. The dynamics of the F-actin probability distribution function ρ in the network are treated by a master equation including branching, spontaneous nucleation and severing: ∂ ρ (x, y, t) ∂ t = ∫ 0 l (y, t) k b r (x, y + y ′) 2 N (t) [ ρ (x - y ′, y + y ′, t) + ρ (x + y ′, y + y ′, t) ] d y ′ ︸ branching + k n u c (x, y, t) ︸ nucleation - k s e v ρ (x, y, t) ︸ severing. (1) Here kbr (x, y) is the branching rate constant, which gives the rate of branching per unit length of F-actin, per molecule of NPF in the membrane. The rate knuc (r, y, t) describes the amount of actin generated per unit time by spontaneous nucleation, while the decay rate ksev is assumed to be controlled by cofilin-driven severing; l (y, t) is the projection onto the y-axis of the length of the filaments added at each time step. The coordinate x is in the plane of the membrane, y is perpendicular to the membrane, and t is the time variable. We use a frame of reference in which the existing actin filaments are stationary, so convective terms are not required. N (t) is the number of molecules of the NPF on the membrane. We assume that the NPFs are uniformly distributed over the membrane, either because they diffuse rapidly in the membrane, or because their initial distribution is uniform. We do not explicitly treat the assembly of curvature-generating proteins. Rather, they are included as in Ref. [5], as a contribution to the forces acting on the membrane. In the model, we take all the filaments to have either a 45° or −45° angle with respect to the y direction, in line with the oblique alignment generally found in dendritic networks. New filaments instantly polymerize to a final length, whose projection on the y-axis is l (y, t) (the filament length multiplied by cos45 ° = 1 / 2). The length is determined by the force-dependent polymerization rate and the capping rate. We assume that capping (and thus the growth of a filament to its final length) occurs on time scales faster than the evolution of the invagination. The validity of this approximation is discussed below. At a given time, only new filaments that branch from certain F-actin subunits can increase ρ (x, y, t). These subunits are included in the integral in Eq (1). The two branching directions correspond to the x ± y′ and y + y′ terms in the integral in Eq (1). In practice, considering Fig 2 as an example, we use the dimensionless length l ¯ (y, t) = l (y, t) / (a / 2) normalized by the length of the actin monomer (a) projected onto the y-axis; l ¯ (y, t) is thus the number of subunits in a given new filament. We discretize Eq 1 accordingly (see Eq. S4) to describe the discrete spatial distribution of F-actin in the network. Fig 2 illustrates the case l ¯ (y, t) = 4 from Eq. S4. This two-dimensional version of the model could be applied to actin networks growing on a long strip of NPF, as in Ref [30]. We extend Eq (1) into three dimensions using a cylindrical coordinate system (see S1 Fig) The spatial coordinates become r, θ and y, and Eq (1) becomes ∂ ρ (r, y, t) ∂ t = 1 2 π ∫ 0 l (y, t) d y ′ k b r (r, y + y ′) N (t) ∫ 0 2 π ρ [ R (r, θ), y + y ′, t ] d θ - k s e v ρ (r, y, t) + k n u c (r, y, t), (2) where R = y ′ 2 + r 2 - 2 r y ′ cos θ is the radial coordinate of the base of a branch having its end at radius r. Again, we discretize the time and spatial dependence, as described in Eq. S5. More details of the simulation procedure are described in the Supplementary Material. This azimuthally symmetric ME approach is directly applicable to templated-nucleation experiments of the type described in Ref. [31]. We assess the validity of the ME by examining how it treats a basic mechanochemical problem: an actin network pushing an obstacle that exerts a constant force opposing polymerization. The obstacle is coated with a ring-like NPF region, mimicking the Las17 ring in Ref. [12]; similar results are obtained for other NPF distributions, including a rectangle and a complete circle. In order to focus on the treatment of branched actin network growth, we ignore negative-feedback effects [12] of actin onto the NPFs. We compare ME results with stochastic simulation results. We start the calculation with 400 filaments, each of which has 50 subunits. For each force value, we run calculations for 200 seconds, and for the stochastic case run an ensemble of 100 simulation runs. The results, shown in Fig 3, show that the ME method agrees quantitatively with the stochastic simulations. We find that F-actin count depends linearly on the external force, and that the velocity is independent on the external force. Both findings are also consistent with the previous stochastic study in Ref. [32]. As mentioned above, key approximation of the ME is that filaments are assumed to grow instantaneously to their final lengths l (y, t), controlled by capping and force. This approximation is valid if the time scale of capping is much shorter than the characteristic time over which the actin count varies in the process of interest. The capping rate is on the order of 1s−1 in budding yeast [33], so the approximation should hold reasonably well for processes occurring on time scales of several seconds or more. This holds for the endocytosis system and model studied in Ref. [12], since the time scale of invagination is on the order of ten seconds. In addition, in our previous study, the assumption of instantaneous polymerization/capping gave results very similar to those of the explicit-polymerization methods when the same parameters were used, as shown by the “Four Variable” model in the Supplemental Material of [12]. Thus the range of validity of the ME approach includes the endocytosis model studied below, but otherwise will vary from case to case. We apply the ME method to CME by solving the actin network dynamics described by Eq (2), while simultaneously calculating membrane shape dynamics using the analysis of Ref. [5]. First, many possible mechanical equilibria corresponding to a range of force values are obtained as in Ref. [5]. These equilibria are given as a shape function for each value of the force. Then, the actin polymerization dynamics are calculated using the ME according to the shape function, and the process is repeated. We focus on branching induced by the NPF Las17, which is the strongest one in budding yeast. Thus N = Las17 count. Recent superresolution data indicate that the NPF Las17 accumulates in a ring-shaped region, while the protein Sla2, which links F-actin to the cell membrane, accumulates at a spot inside the ring [34]. Therefore we divide the membrane into a ring of pushing forces corresponding to Las17 and a spot of pulling force corresponding to Sla2. The pushing forces are generated by the growth of branched filaments in the network. The pulling forces act on spontaneously nucleated filaments, which are assumed to form a passive layer attached to both Sla2 and the branched network. Sla2, with the membrane attached to it, is thus pulled back with the retrograde flow of the network (see Fig 1), as suggested in Ref. [11]. Actin polymerization also occurs mainly near the membrane [35]. We thus use the following form for the branching rate constant: kbr (r, y) ={ kbrmaxexp[ − (y−yL) 22σbr2 ]ifyL<y<yL+ybrandrLin<r<rLout0otherwise (3) where k b r m a x is the maximum value of kbr, σbr is the width of the branching region (a precise definition is given in the Supplemental Material), ybr is a cutoff imposed for numerical convenience, and r L i n and r L o u t are the inner and outer radii of the Las17 ring. Eq. (S33) gives the corresponding formula for knuc. In practice, we use dimensionless rates k ¯ b r m a x and k ¯ n u c m a x, as in Eqs. S4 and S5. The spatial branching and nucleation functions are shown in Fig 4. Note that even though kbr cuts off sharply at y = yL, the dynamics of Eq 2 will result in a small component of ρ penetrating past yL. This portion of the ρ is used to calculate the pushing force of the actin onto the membrane below. In our previous work [12], we demonstrated a crucial negative-feedback effect of actin branching on the Las17 count, and here we treat the Las17 dynamics using a similar rate equation: dN (t) dt=N (t) 2[ Nfull−N (t) ]−αNFbr (t), (4) where Nfull is the maximum possible count of Las17 (from 2-d packing considerations), α is the probability that a branching event will cause Las17 to dissociate from the membrane (thus being inactivated), and F b r (t) = 2 π ∫ y L y L + y b r d y ∫ r L i n r L o u t (a 2) k b r (r, y) ρ (r, y, t) r d r, (5) is the number of new branches created per unit time per Las17 molecule. The probability α reflects the strength of binding of the Las17 to the membrane, and a = 2. 7 nm is the step size per added subunit. The factor of a / 2 is the projection of a onto the y-axis. The forces from actin polymerization deform the membrane from one mechanical equilibrium to another, as indicated in Fig 5. The use of equilibrium shape functions is justified, because the kinetics of the actin-membrane system are determined by the slowly varying actin network shape (timescale > 1s), but relaxation of the membrane occurs much faster, at the speed of sound (timescale ∼ 0. 001s). Because the net force on the actin network is exceedingly small [4], the pushing force from the actin network must balance the pulling force, as shown in S2 Fig. The pushing force is generated by network growth in the outer region comprising the Las17 ring, and is thus denoted fout. It is calculated by allowing the actin network to protrude slightly into the membrane, according to Eq 2, and imposing a linear repulsive force between this portion of the actin and the cell membrane (see Supplemental Material for details). The pulling force is exerted in the inner region corresponding to the Sla2 spot, and is thus denoted fin. The balancing forces fin and fout produce a deformation described by the membrane shape function ym[rm (s) ], where ym is membrane height and rm is the radial coordinate in the membrane. During each time step, actin first polymerizes according to Eq 2. The membrane deformation, or invagination depth yI is then determined from the extent of actin polymerization by a procedure implementing a “molecular clutch” based on the amount of F-actin. The possibility of such a mechanism is supported by findings [36] that a clutch transmits forces from the actin cytoskeleton to the extracellular matrix or other cells. A clutch should also be present in CME because there must be a transition in mechanical behavior with increasing F-actin count. When the F-actin count is small, there is very little actin material in the central region of the endocytic patch. Therefore there is almost nothing for the outer filaments, which are moving backwards in retrograde flow, to “grab” onto. This makes it impossible for the growing network to exert a pulling force. On the other hand, when the F-actin count is larger, there is enough material at the center to transmit the force generated by the growing filament to the endocytic coat proteins We implement the clutch as follows. Up to a certain minimum value of the F-actin count, Fmin, yI is taken to vanish; for larger values of F, the actin network is assumed to be completely rigid, and yI is driven by the difference in polymerization rates between the outside and the inside. Details are given in the Supporting Information. Given yI, the deformation profile is chosen from the pretabulated set (see S3 Fig) as the one with invagination y I p r e closest to yI. The force is chosen by linear interpolation between the force of the profile with y I p r e and that of another profile with y I p r e ± 0. 1 R Π, so that y I p r e - 0. 1 R Π < y I < y I p r e + 0. 1 R Π. The updated shape function and force are used in the next step to determine the branching region, spontaneous-nucleation region and dynamics of actin polymerization (see Fig 5). This approach should describe the dynamics of the membrane deformation well, since the differences between successive membrane shape functions are relatively small. The initial actin distribution is a ring of filaments represented by ρ (r, y, 0) = (1 / 2 π r) δ (r - r L i n) θ (y - y L) θ (y L + l m a x - y), where r L i n is the inner limit of the Las17 ring. Further, at the beginning of the simulation L (0) = 20, yS = yL = 0, and all forces vanish.
  We use the experimentally measured time courses of F-actin (F) and Las17 number (N) [12] as our fitting targets. We use the four quantities (k0, k ¯ b r m a x, ksev, and α) (see Table 1) as our fitting parameters, and regard the rest of parameters as “fixed” (see Table 1 and S1 Table). The values of the fitting parameters are obtained by minimizing the mean-square difference between the measured F and N time courses on one hand, and the model and the experimental data on the other hand: ϵ=1nN∑i=1nN[ Nmod (ti) −Nexp (ti) ]2[ max (Nexp) ]2+1nF∑i=1nF[ Fmod (ti) −Fexp (ti) ]2[ max (Fexp) ]2, (6) while keeping the “fixed” parameters unchanged. Here nN and nF are the number of experimental data points of Las17 and F-actin in the time courses. Fexp is obtained from measurements of the time course of Abp1, as in Ref. [12]. At each step of the the fitting process, we randomly vary (k0, k ¯ b r m a x, ksev, α) and calculate ϵ. The new values are accepted if ϵ is lower. The above computation is repeated until ϵ does not decrease despite a large number of attempts (typically about 300). Then the values of the fitting parameters are found for a given set of “fixed” parameters. This process requires a large number of trial calculations (about 300) for one set of the “fixed” parameters, which is nevertheless manageable within the ME method. In the stochastic simulations, one needs to repeat the calculation about 1000 times for the same fitting parameter values to obtain adequate statistics. Thus, a total in the range of 300,000 runs are needed for one set of “fixed” parameters, which is a very demanding computational load. In Ref. [12], we estimated the fitting parameters by first pre-fitting a simplified four-variable rate-equation model to the experimental time courses, and then fitting the stochastic model to experimental maxima and lifetimes starting with the pre-fitted parameter values. This process is less efficient and accurate than the automatic fitting process used here, which is difficult to incorporate into the stochastic model. There are two additional “fixed” parameters that are estimated via either other experimental data or physical constraints (see Table 1). The zero-force dimensionless filament length l ¯ m a x (the maximal number of subunits) in Table 1 is estimated as l ¯ m a x = k o n G / k c a p, where kon = 11. 6μM−1 s−1 [37] is the on-rate constant, G = 5. 3μM [33] is the free-actin concentration, and the capping rate kcap is taken to be 1 s−1 [33]. The spontaneous nucleation rate parameter k n u c m a x is fixed by a combination of two constraints: i) that adequate pulling forces can be generated, and ii) that sufficient invagination can be obtained. Additional parameters, including those describing the geometry, are given in S1 Table. Fig 6 shows how the actin network invaginates the membrane over time. The membrane forms an ‘Ω’ shaped invagination after 16 seconds into the simulation or about 5 seconds after actin polymerization starts, consistent with observations in electron micrographs [15]. The time courses of the F-actin count F and Las17 count N also reveal a good fit to the experimental data in Ref [12], shown in Fig 7A. Note that the actin distribution extends slightly below the plane of the membrane, and outside the Las17 ring. This occurs because of the nonlocal dynamics of Eq 2. The spreading outside the Las17 ring is physically expected because of the nonzero filament length. The portion below the membrane plane is an approximation used to calculate the pushing force generated by the actin, as described in the Supplemental Material. The three-dimensional distributions predicted here could be tested by superresolution microscopy methods with resolution on the scale of tens of nanometers. Such methods [34] have found that F-actin forms a hemispherical shape, and Las17 forms a ring, as in the present model. Electron microscopy data in the literature [15,38,39] are also consistent with an F-actin hemisphere. The “acidic” regions of Las17 and other yeast NPFs are believed to control their binding to Arp2/3 complex, and therefore their NPF activity. The mutant containing mutations of both the Las17 and the NPF Pan1, las17 pan1Δacidic (abbreviated as LPΔA) should have a strong reduction in the Las17 branching activity. We choose this mutant to avoid possible compensatory effects from the nucleation-promoting activity of Pan1. As in Ref. [12], we model the mutation via a 40% reduction in k b r m a x. Fig 7B shows that the model matches the measured LPΔA phenotype [7] well, with an accuracy comparable to that of our previous model [12]. Remarkably, the F-actin count is actually increased by reducing kbr. This counter-intuitive phenotype results from a competition between a direct effect and an indirect effect. The direct effect is the reduction in branching rate per molecule of Las17 caused by the mutation. The indirect effect is the resulting increase in Las17 caused by the reduced branching rate, due to the negative-feedback effect described in Eq (4). This increase will tend to increase the F-actin count. For the conditions considered here, the indirect effect outweighs the direct one. As indicated in Fig 1, actin filaments in the central region are required to exert pulling forces. We assume that these filaments arise from spontaneous nucleation (not requiring NPFs at the endocytic site). Possible sources of the spontaneously nucleated actin filaments can be severed filament fragments [40], or nucleation via Dip1, which is independent of NPFs [41]. A minimum value of k ¯ n u c m a x is required to exert adequate pulling force, since reducing k ¯ n u c m a x reduces the number of filaments in the central region (see Fig 8B–8D). This increases the pulling force per filament, and eventually causes them to detach from the membrane. The maximum pulling force that a membrane-attached filament can sustain is not known. But Ref. [42] gives a quantitative measurement (> 40pN) of the rupture force between a single actin filament and the crosslinking proteins filamin and α-actinin in vitro, at low loading rates. In red blood cells [43], the interaction force between the actin cytoskeleton and the membrane was found to be ∼10pN per filament in a model fitted to experimental data. However the cytoskeleton-membrane interactions in yeast could be very different from those in red blood cells. We thus based our estimate of the maximum actin filament pulling force on the measured rupture force. To estimate the pulling force per filament in our model, we divided the total maximum pulling force of 725 pN obtained from the membrane-energetics analysis (Supplementary Material) by the number of pulling filaments, estimated as the total F-actin count inside r L i n divided by 2 y n u c / a, a dimensionless measure of the height of the nucleation layer. In the default model, this procedure gave a pulling force of ∼30pN per filament, below the rupture force [42]. When spontaneous nucleation was suppressed by reducing k ¯ n u c m a x, the pulling force exceeded the rupture force, as shown in Fig 8A. Thus a minimum rate of spontaneous filament nucleation in the central region is required for the actin network to pull the membrane without rupture of the actin-membrane interactions. This prediction might be tested by deletion of the protein Dip1, which could participate in an NPF-independent actin polymerization pathway [41]. If Dip1 nucleates filaments in the central region, its deletion should have two main effects. First, the actin hemisphere should be sparser in the middle. as shown in Fig 8B. This could be verified by superresolution images of dip1Δ cells. Second, in the dip1Δ cells, reduced nucleation should lower the number of pulling filaments (see Fig 8A), causing rupture if the force per filament exceeds the rupture force. This would reduce the efficiency of invagination. Ref. [41] found that in dip1Δ cells, only 40% of the total patches were internalized, compared to ∼90% in WT cells. Internalization of 30% of the patches was also delayed delayed by over 20 seconds. Both of these effects could be due to the reduced number of filaments in the central region. On the other hand, we find that too large a magnitude of k ¯ n u c m a x also disrupts invagination. We increased k ¯ n u c m a x by a factor of 2 and refitted. In comparison with the default invagination y l m a x = 54 nm, we obtained a maximum invagination y l m a x = 30 nm, which we consider to be a failed event. Therefore, possible values of k ¯ n u c m a x are in a range limited by the constraints of i) adequate pulling force and ii) adequate invagination. We investigate the responses of the actin-membrane system to the drugs CK-666, which inhibits branching, and Latrunculin A (LatA), which inhibits polymerization and to some extent branching. S4 Fig summarizes the effects on Fmax and y I m a x of a broad range of combinations of parameter changes. Note than in S4 Fig frame D, breaking of pulling filaments bearing very large loads is not included. Inclusion of this effect will lead to failure of invagination at small k ¯ n u c m a x, as discussed above.
  The ME enables a “one-shot” approach to stochastic actin network dynamics by treating the pdf. It represents a statistical average of many ensembles of the same actin network calculated by the equivalent stochastic simulation, which is referred as “the infinite population limit” in Ref. [47]. We feel that the use of a statistical average is legitimate, because endocytosis in yeast is highly stereotypical. The behavior of the actin network typically displays a single mode value [7] (the most likely value) plus a range of fluctuations, instead of having multiple mode values. The “one-shot” nature of the ME simplifies the process of making modifications and refitting the model to experimental data. For both ME and stochastic methods, calculation of branching processes causes the largest computational load. The ME approach is most efficient when the branching region is small because it scans the physical space for possible branching events. The stochastic simulation, on the other hand, is most efficient when the total number of filaments is small, because it scans the subunits for possible branching sites. Thus, the ME is computationally most powerful for treating dense actin networks with relatively small regions of NPF, especially with spatial symmetry. The ME approach could be applied to many other problems involving branched network growth, as long as the branching is generated by a flat NPF region at a membrane or surface. For cell migration, the cytoskeleton can sometimes be simplified to a two dimensional network [48] or even a one dimensional network [49] with NPFs close to the cell membrane. The ME method is well suited to this type of problem. Notice that because the ME method is mechanochemical, it could be used to describe mechanical feedback effects on cell migration, if combined with a “Helfrich”-type calculation of the cell membrane forces. The ME method is also well-suited for treating filopodium and lamellipodium geometries where the membrane forces can be calculated straightforwardly from the membrane geometry. In order to develop a practical method for treating a problem as complex as endocytosis, we have made several simplifying assumptions. We have developed a computationally efficient master equation (ME) approach for calculating the spatial distribution of F-actin branched networks growing in the presence of mechanical forces. The approach was validated by comparison with stochastic-simulation results. It was then used to develop a mechanochemical model of clathrin-mediated endocytosis in yeast (CME), which treats both the actin network and the cell membrane realistically. The mechanochemical model was used to reveal the time evolution of the actin-membrane system during CME, to quantitatively estimate unknown parameter values, and to predict several important mechanisms in CME that are unseen or omitted in previous models. These predictions provide possible directions for experiments in CME, especially for superresolution microscopy and drug treatments. Beyond CME, the new ME approach provides possible applications to a wide range of problems involving the spatial distribution of branched actin polymerization.
summary: Endocytosis is a dynamic process by which cells internalize substances from outside the cell. Especially in yeast, endocytosis is mechanically demanding due to the high pressure difference across the cell membrane, or turgor pressure. Polymerization of a branched actin network is the major process providing the mechanical force to overcome the turgor pressure. Understanding the kinetics of the actin network, and the mechanical interaction between the actin network and the cell membrane, is thus crucial for the study of endocytosis. We develop an efficient mathematical framework for actin dynamics that can realistically incorporate these two features, thus providing a practical method for quantitatively modeling actin dynamics during endocytosis. The resulting model mechanistically reveals that spontaneous nucleation at the center of the endocytic site is required for successful endocytosis, distinguishes the roles of branching and polymerization, and predicts several other experimentally testable outcomes. The accuracy and efficiency of the method, in describing both mechanics and chemistry, render it applicable to a broad field of membrane-bending processes.
keywords: cell motility
  actin filaments
  classical mechanics
  cell processes
  condensed matter physics
  cellular structures and organelles
  contractile proteins
  polymer chemistry
  actins
  actin polymerization
  proteins
  chemistry
  polymerization
  cell membranes
  physics
  endocytosis
  biochemistry
  cytoskeletal proteins
  cell biology
  secretory pathway
  biology and life sciences
  chemical reactions
  physical sciences
  nucleation
title: A master equation approach to actin polymerization applied to endocytosis in yeast