 Research article
 Open Access
 Published:
Analytical model for macromolecular partitioning during yeast cell division
BMC Biophysics volume 7, Article number: 10 (2014)
Abstract
Background
Asymmetric cell division, whereby a parent cell generates two sibling cells with unequal content and thereby distinct fates, is central to cell differentiation, organism development and ageing. Unequal partitioning of the macromolecular content of the parent cell which includes proteins, DNA, RNA, large proteinaceous assemblies and organelles can be achieved by both passive (e.g. diffusion, localized retention sites) and active (e.g. motordriven transport) processes operating in the presence of external polarity cues, internal asymmetries, spontaneous symmetry breaking, or stochastic effects. However, the quantitative contribution of different processes to the partitioning of macromolecular content is difficult to evaluate.
Results
Here we developed an analytical model that allows rapid quantitative assessment of partitioning as a function of various parameters in the budding yeast Saccharomyces cerevisiae. This model exposes quantitative degeneracies among the physical parameters that govern macromolecular partitioning, and reveals regions of the solution space where diffusion is sufficient to drive asymmetric partitioning and regions where asymmetric partitioning can only be achieved through additional processes such as motordriven transport. Application of the model to different macromolecular assemblies suggests that partitioning of protein aggregates and episomes, but not prions, is diffusionlimited in yeast, consistent with previous reports.
Conclusions
In contrast to computationally intensive stochastic simulations of particular scenarios, our analytical model provides an efficient and comprehensive overview of partitioning as a function of global and macromoleculespecific parameters. Identification of quantitative degeneracies among these parameters highlights the importance of their careful measurement for a given macromolecular species in order to understand the dominant processes responsible for its observed partitioning.
Background
Reproduction of the yeast Saccharomyces cerevisiae by budding is a classical example of asymmetric cell division. The future daughter yeast cell grows as a bud on the surface of the mother cell, receiving macromolecules and organelles that escape from the mother compartment through the narrow opening of the bud neck (Figure1). Movement of macromolecular content from the mother to the bud has been heavily investigated in the last decades. These studies identified various pathways that use motor proteins to transport and regulate the partitioning of organelles such as the nucleus, vacuoles, endoplasmic reticulum, peroxisomes, mitochondria and secretory vesicles (for reviews, see [1][3]).
Cells also contain other types of large structures that are close in size to membranebound organelles. Some are typically composed of both RNA and proteins and therefore referred to as ribonucleoprotein (RNP) bodies or RNA granules [4][6]. Protein granules can be composed of metabolic enzymes, protein degradation machinery or misfolded/damaged proteins in complex with chaperones [7][10]. Some of these assemblies exist constitutively, whereas others form in response to changes in the environment. Moreover, the nucleus contains extrachromosomal circular DNA (episomes), of which the extrachromosomal ribosomal DNA circles (ERCs), formed via recombination in the rDNA locus, are a prominent example. Notably, some of these large macromolecular assemblies, namely aggregates of misfolded/damaged proteins and episomes, are largely excluded from the bud during yeast cell division [11][13]. The mechanisms responsible for asymmetric partitioning of protein aggregates and episomes are a subject of ongoing discussion.
Asymmetric partitioning of oxidatively damaged protein aggregates is linked to yeast replicative aging [13],[14]. This partitioning appears to involve the actin cytoskeleton [13],[15],[16], possibly even in retrograde movement of aggregates from the bud into the mother [15],[17]. An alternative model, not involving active transport, suggests that slow diffusion of protein aggregates and yeast cell geometry may already largely account for aggregate retention in the mother cell [18]. Additionally, localization of protein aggregates to specialized compartments the juxtanuclear quality control compartment (JUNQ) and the insoluble protein deposit (IPOD) inclusions that undergo asymmetric segregation during cell division appears to provide a further contribution to aggregate exclusion from the bud [19][20].
In early studies, the binding of episomes to structures retained in the mother cell was hypothesized to explain their exclusion from the bud [11], thus contributing to the rejuvenation of the newly born cell [12]. This idea was later revived by a report that nuclear pores — large protein complexes embedded in the nuclear envelope — are retained in the mother compartment during cell division and bind ERCs [21]. However, several studies showed that nuclear pores are not retained in the mother cell [22][26] and that tethering to nuclear pores actually promotes segregation of episomes into the bud [27]. Moreover, using live cell tracking of episomes and 3D simulations, Gehlen et al. demonstrated that slow passive diffusion of episomes and the geometry of the yeast nucleus are largely sufficient to explain episome exclusion from the bud [28].
Here we argue that passive diffusion should always be considered when analyzing the partitioning of macromolecular assemblies during cell division. To facilitate the analysis of the role of diffusion, we developed an analytical model of diffusiondriven transfer of macromolecules between the mother and bud compartments of budding yeast. The model is based on a description of diffusive transfer according to the “narrow escape” approximation, which gives the mean time required for a diffusing particle to escape from a compartment through a narrow opening [29][32]. Application of this approximation allows conversion of a 3D diffusion problem in a complicated and dynamic geometry to a simple compartmental model for the transfer that is expressible in terms of a single nonlinear ordinary differential equation. The accuracy of this model was confirmed through comparison with the results of previous stochastic simulations [18],[28]. The phase space of all possible solutions of this simple model permits immediate examination of regions for which diffusion is sufficient to drive asymmetric partitioning and those where it is not (therefore implying the presence of other processes like active motordriven transfer to account for asymmetric partitioning). Finally, important degeneracies among the physical parameters that control the segregation of macromolecular content are revealed by the model.
Methods
Bidirectional motherbud transfer assuming a static geometry
According to the “narrow escape” approximation, the escape of particles with diffusion constant D from a spherical volume V (radius R) through a narrow circular opening of diameter d can be asymptotically described by an exponential decay with a mean escape time of:
where K _{ e } corresponds to the firstorder rate constant for escape [29][32]. This approximation is exact in the limit of d < R[29], with the mean escape time accurate to the percent level even for d ≈ R[31]. In budding yeast, the cytoplasmic content of the mother cell (of diameter ~4–6 μm) does not flow automatically into the bud compartment. Instead, macromolecules and organelles “escape” the mother compartment through the narrow constriction of the bud neck (of diameter ~1 μm), which can be mediated by either passive or active transfer processes. For transfer between the two compartments of the mother and bud, the rate constants for successful transfer between the two compartments are a factor of two lower than the rates of first escape from either compartment due to additional recrossing events at the bud neck [33],[34]. Applying the “narrow escape” approximation in a bidirectional sense to yeast budding leads to the following differential equation governing bidirectional motherbud transfer via passive diffusion:
with N _{ m } and N _{ b } the respective particle numbers in mother and bud compartments; K _{ f } = dD _{ m }/V _{ m } referring to forward transfer (mother to bud) and K _{ r } = dD _{ b }/V _{ b } to reverse transfer; D _{ m } and D _{ b } the compartmentspecific diffusion constants (reflecting possible differences in compartmental viscosities); and d the bud neck diameter. Total particle conservation, N _{ T } = N _{ m } + N _{ b }, implies dN _{ b }/dt = − dN _{ m }/dt. Note that the steady state solution (for equal compartment viscosities), N _{ b }/V _{ b } = N _{ m }/V _{ m }, satisfies the homogeneous concentration expected for diffusion in an enclosed geometry in the absence of particle sources or sinks. Rewriting Eq. 2 in terms of concentrations, yields:
Note the final equality of the rate constant coefficients for forward and reverse transfer in the case of equal compartmental viscosities, D _{ m } = D _{ b } = D. Generalization of Eq. 3 to
allows unified treatment of

(I)
passive diffusion with equal compartmental viscosities (D _{ m } = D _{ b } = D), for which k _{ r } = k _{ f } = K _{ f } = dD/V _{ m };

(II)
passive diffusion with different compartmental viscosities (D _{ m } ≠ D _{ b }), for which k _{ f } = K _{ f } = dD _{ m }/V _{ m } and k _{ r } = K _{ f } D _{ b }/D _{ m } = dD _{ b }/V _{ m };

(III)
active transport, for which k _{ f } and k _{ r } (as defined with respect to the mother) are determined by directional trafficking processes and can differ in magnitude.
Anomalous diffusion, ubiquitously observed in biological systems [35],[36], can also be considered in our model. In this case, the narrow escape approximation should still apply, requiring only appropriate modification of the relationship given in Eq. 1 between the mean escape time t _{ e } and the now anomalous diffusion constant, with this relationship dependent on the exact mathematical definition of anomalous diffusion.
Transport by convection (bulk flow or cytoplasmic streaming) can also influence macromolecular partitioning [37],[38]. This phenomenon is apparent in large cells such as oocytes or plant cells [39],[40]. However, to our knowledge, the role of bulk flow in yeast cell division has not been considered. In principle, cytoplasmic streaming between mother and bud compartments could result from bud growth during the division process or directed transport of vesicles and organelles into the bud. Our model provides a good testing ground for the contribution of such processes to macromolecular partitioning, especially in cases where diffusive transfer is insufficient to account for the experimentally observed partitioning.
Bidirectional motherbud transfer with a growing bud
In the above, we have implicitly considered a static motherbud geometry. Assuming a linear growth of the bud to a final fraction σ of the mother volume at the cell division time T (${V}_{b}\left(t\right)=\sigma \phantom{\rule{0.22em}{0ex}}{V}_{m}{\scriptscriptstyle \phantom{\rule{0.22em}{0ex}}\frac{t}{T}}$), gives for Eq. 2:
where we have implicitly assumed that bud growth is slow compared to the diffusive homogenization timescales within each compartment. Defining x = N _{ m }/N _{ T }, y = N _{ b }/N _{ T } = 1 − x, τ = t/T, α = k _{ f } T, and β = k _{ r } T/σ leads to the dimensionless form:
with final partitioning determined solely by the forward transfer efficiency α and the reverse transfer efficiency β. For passive diffusion, α = dD _{ m } T/V _{ m } and β = dD _{ b } T/(σ V _{ m }).
Influence of binding sites on bidirectional motherbud transfer
Binding sites in the bud (target sites) or in the mother (retention sites) can bias macromolecular partitioning. In the simplest scenario, binding sites are assumed far from saturation with on/off binding fast compared to the transfer timescales. For retention sites, the unbound fraction is then:
with B _{ m } the binding site concentration and K _{ d } the dissociation constant. Similarly, for target sites, the unbound fraction is:
Only this free fraction can undergo transfer (either passive or active). Such binding sites are accounted for in our model through inclusion of the free fractions in the forward and reverse transfer efficiencies: α = ϕ _{ m } k _{ f } T and β = ϕ _{ b } k _{ r } T/σ. In addition to binding sites, transient trapping of macromolecules on membranes or within organelles can similarly be accounted for through the reduction in their compartmental free fractions.
Model validation
The accuracy of our model is revealed by comparison to previously published 3D Monte Carlo simulations. For their 3D simulation of passively diffusing protein aggregates, Zhou et al. [18] obtained a bud partitioning of approximately 10% (determined from Figure 7(A) in Zhou et al. [18]), serving as a useful benchmark for our model. Using their parameters (σ = 0.61, R = 2.5 μm, d = 1.25 μm, D = 0.001 μm^{2}s^{‐ 1}, T = 90 min), we obtain 8.4%. However, Zhou et al. [18] did not include bud growth in their simulation. Neglect of bud growth, within the “narrow escape” approximation, requires direct integration of Eq. 2 for a static K _{ f } = dD/V _{ m } and K _{ r } = dD/V _{ b }. This leads to a simple exponential decay of the mother fraction:
Using their values, we obtain a final bud fraction of 9.0%, which is in good agreement with their value of approximately 10% (with both fractions far from the maximal value of 38% expected for full volumebased equilibration between the two compartments).
For their 3D simulation of nuclear episomes, which explicitly includes linear growth of the bud nuclear volume and is in all respects identical with the assumptions of our model, Gehlen et al. [28] obtained a bud nucleus partitioning of 25%. Using their parameters (σ = 0.64, mother nuclear radius R = 0.975 μm, d = 0.5 μm, D = 0.004 μm^{2} s^{‐ 1}, mitosis duration T = 9 min), we obtain a 17% bud partitioning, which is in reasonable agreement with the 25% of Gehlen et al. (with both values showing a significant deviation from the maximal result of 39% expected for full volumebased equilibration between the two compartments).
As further confirmation of the underlying assumptions of our model, we note that the intercompartmental gradients at the time of division in the displayed contour maps of both Zhou et al. [18] and Gehlen et al. [28] are confined to a narrow boundary layer at the bud neck with an otherwise flat concentration spanning each individual compartment; this importantly validates the “narrow escape”based assumption of efficient compartmental homogenization [31].
While these confirmatory results and observations are encouraging, more extensive comparisons of our analytical model with results from 3D simulations are necessary to determine its level of accuracy for different parameters and division geometries. Our model assumes the asymptotically valid “narrow escape” approximation of a uniform exponential distribution for particle transfer between the two compartmental volumes. This allowed us to convert a complicated 3D diffusion problem to a much simpler system of firstorder transfer between the two compartments. While the “narrow escape” approximation has been shown to be highly accurate for prediction of the mean escape time even for sizeable escape windows [31], the accuracy of the assumption of a uniform exponential escape (or, in our case, transfer) distribution still remains to be rigorously demonstrated for escape windows of different sizes [41].
Noisy partitioning for low numbers of particles
Reverting from fractions to numbers of molecules in Eq. 6 allows for convenient examination of noise in inheritance:
The dynamics of a single particle that can be present at any given time either in the mother (Q _{ m } = 1,Q _{ b } = 0) or in the bud (Q _{ m } = 0,Q _{ b } = 1) is given by:
The summed dynamics of N _{ T } independent particles is therefore:
As Eq. 12 is formally equivalent to Eq. 10, consideration of multiple particle dynamics in our model requires only knowledge of the probability of individual particle transfer. Due to this particle independence, the final motherbud distribution probabilities for N _{ T } particles are given by the corresponding terms in the binomial expansion of the single particle retention/transfer probabilities, ${\left(x+y\right)}^{{N}_{T}}$. The first term in this series, ${x}^{{N}_{T}}$, corresponds to the frequency of null transfer f with which all particles are retained in the mother compartment upon completion of cell division. For models in which mother/bud particle concentrations affect transfer rates (particles are not independent), simulations using the Gillespie algorithm [42] would be required to obtain transfer distributions for each N _{ T }.
Results and discussion
Application of the “narrow escape” approximation [29][32] in a bidirectional sense to the dividing yeast cell yields a simple compartmental model for diffusive particle transfer between the mother compartment (m) and the growing bud (b) (Figure 1 and Methods). In this model, the final partitioning of macromolecular content is determined by two dimensionless parameters — the forward (mother to bud) transfer efficiency α and the reverse (bud to mother) transfer efficiency β — that depend on the forward and reverse transfer rate constants (k _{ f } and k _{ r }), the duration of cell division (T) and the final budtomother volume ratio (σ = V _{ b }/V _{ m }). The effect of binding to nonpartitioned structures (e.g. the cell cortex) is accommodated in the model through the concomitant reduction of free fractions in each compartment (ϕ _{ m } and ϕ _{ b }) (Figure 1 and Methods). For passive diffusion, the transfer rate constants are defined, according to the narrow escape approximation, by the compartmental volumes (V _{ m } and V _{ b }), the bud neck diameter d and the compartmentspecific diffusion constants (D _{ m } and D _{ b }), yielding for the transfer efficiencies: α = ϕ _{ m } dD _{ m } T/V _{ m } and β = ϕ _{ b } dD _{ b } T/(σ V _{ m }). This formulation achieves a reduction of eight parameters to the two dimensionless transfer efficiencies α and β, which implies a high level of parameter degeneracy and shows that the same partitioning can be achieved with different parameter combinations. Interestingly, the parameters can be categorized into two groups: parameters that globally influence the partitioning of any molecular content (d,T,V _{ m },σ) and contentspecific parameters (D _{ m },D _{ b },ϕ _{ m },ϕ _{ b }).
The model allows calculation of partitioning trajectories over the course of cell division for any given set of parameters, including values of α and β for which partitioning does not reach equilibrium. Examples of temporal trajectories for different values of α and β are displayed in Figure 2a. Partitioning equilibration is not achieved for values of α and β roughly less than 1, defining a critical diffusion constant D _{ c } ≈ V _{ m }/(dT). Significant retention of molecular species i in the mother cell can be expected for D _{ i } < D _{ c } (Figure 2b,c).
To explore the model, we calculated the bud fraction of macromolecules at the end of the division process y = N _{ b }/N _{ T } for all possible sets of transfer efficiencies (Figure 3a). Note that the ratio ρ = α/β defines the asymptotic limit $y={\scriptscriptstyle \frac{\rho}{1+\rho}}$ along each 45° diagonal. The asymptotic limit can be reached by sufficiently slowing down cell division (increasing T) such that bidirectional transfer is efficiently equilibrated at each time point. Figure 3b shows the extent of compartmental equilibration, defined as the final bud concentration divided by the fully equilibrated concentration obtained in the asymptotic limit (along the corresponding 45° diagonal), attained at each α and α in the phase plot.
The blue diagonal in these plots (Figure 3) gives the set of solutions for particles undergoing unbiased passive diffusion in the case of final bud volume equal to the mother volume (corresponding to y _{eq} = 0.5 upon complete equilibration). To illustrate the scope of our model, consider the shift from the solution at α = 1, β = 1 (blue circle) to that at α = 0.1, β = 0.1 (blue square). This displacement could be achieved in several different ways: a tenfold increase in the mother volume; or a tenfold reduction in bud neck diameter, diffusion constant, or division time. Tenfold reductions of both compartmental free fractions due to efficient retention/target site binding would have the same effect. Our model therefore reveals the exact quantitative degeneracy among these various parameters, obviating much of the need for laborious stochastic simulations of different scenarios (e.g. Zhou et al. [18] and Gehlen et al. [28]). Note that particle enrichment in the bud (y > y _{eq}, corresponding to a higher concentration of particles in the bud) is, of course, not obtainable with unbiased passive diffusion. Solutions above or below the diagonal, the latter including the possibility of bud enrichment, are obtained with different compartmental diffusion constants (D _{ m } ≠ D _{ b }), asymmetric distribution of binding sites (leading to ϕ _{ m } ≠ ϕ _{ b }) or, more generally, with different transfer rate constants (k _{ f } ≠ k _{ r }) for the case of active transport between the compartments.
The offset red diagonal in Figure 3 gives the solution set more appropriate for S. cerevisiae, in which the final bud volume is 2/3 the mother volume (here equilibration corresponds to y _{eq} = 0.4). Again, only solutions below the red diagonal cover the possibility of particle enrichment in the bud (y > y _{eq}). Due to parameter degeneracy, an increase in the bud neck opening, an extension of the cell division duration, or a corresponding decrease in the cell size will move the partitioning solutions globally along the diagonal towards equilibration. Moleculespecific equilibration can, for example, be achieved by decreasing the number of binding sites in both mother and bud compartments.
For each solution set specified by the final budtomother volume ratio σ, fast diffusing particles are partitioned between mother and bud proportional to compartmental volumes, whereas particles with lower diffusion constants are progressively more retained in the mother (Figure 3b). Assembly of rapidly diffusing molecules into larger, less mobile structures could therefore be an important factor in cell fate determination. Misfolded or damaged proteins [10] as well as functional cytosolic content (e.g. metabolic enzymes [7],[8]) can form large supramolecular assemblies. But how slow should diffusion be to limit partitioning in yeast? Assuming typical parameters of S. cerevisiae cell division (mother cell radius R = 2.5 μm, bud neck diameter d = 1.25 μm, and division time T = 90 min) [18] yields a critical diffusion constant D _{ C } of 0.01 μm^{2}s^{1}. For comparison, the diffusion constant of the green fluorescent protein GFP is approximately 10 μm^{2}s^{1} in the yeast cytoplasm [43]. If macromolecular diffusion followed the StokesEinstein relation of D ≠ r ^{− 1}, where r is particle radius, only particles 1000fold larger than GFP (with radii >3 μm and therefore comparable to or larger than a typical haploid yeast cell) would diffuse slower than 0.01 μm^{2}s^{1}. However, macromolecular crowding within the cytoplasm generally leads to a much stronger dependence of diffusion constant on particle size than predicted by the StokesEinstein relation [44][46]. In fact, the diffusion constants of two factors associated with yeast replicative aging, protein aggregates marked with the Hsp104 chaperone and extrachromosomal rDNA circles (ERCs), are close to or significantly lower than the critical diffusion constant D _{ c } of 0.01 μm^{2}s^{1} (0.001 μm^{2}s^{1} for protein aggregates [18], 0.004 μm^{2}s^{1} for ERCs [28]), consistent with the hypothesis that passive diffusion may be sufficient to account for their observed asymmetric partitioning [18],[28] (Figure 3, Methods). Precise measurements of cell geometry, compartmentspecific diffusion constants, and interactions with subcellular compartments such as the endoplasmic reticulum [47], nucleus and vacuole [20] should help clarify the contribution of slow diffusion relative to mechanisms such as active budtomother transport [17],[18] and binding to retention sites in the mother.
Prions present another interesting case where diffusion could in principle limit partitioning. These infectious proteins can misfold and aggregate into βsheetrich amyloid fibers that catalyze the conversion of newly made proteins into the prion form. Prion propagation by selftemplated fiber growth, chaperonemediated fiber fragmentation and transmission of infectious aggregates to the bud during cell division can generate heritable phenotypes [48],[49] when at least one infectious aggregate is inherited by the daughter cell [50]. Experimental studies indicate that the size of prion aggregates affects their transmission from the mother to the bud, but the mechanisms underlying this size selection remain unclear [51]. According to our model, the probability of the bud receiving no aggregates (frequency of null transmission) depends inversely on the diffusion constant of the aggregates, strongly increasing for diffusion constants below the critical value of 0.01 μm^{2}s^{1} (Figure 4, Methods). However, measurements with fluorescence correlation spectroscopy indicate that prion aggregates typically diffuse ~50 fold faster than this critical diffusion constant [52]. Thus, diffusion alone appears insufficient to limit aggregate partitioning. Instead, sizedependent interactions of prion aggregates with protein quality control compartments retained in the mother [19],[20] (decreased ϕ _{ m }) may explain their exclusion from the bud. Further characterization of aggregate dynamics should help to test this hypothesis and to investigate the alternative possibility of active transport mechanisms.
As diffusion can limit the partitioning of macromolecular content in the context of narrow escape, how does the dividing yeast cell ensure equal partitioning of structures such as large membranous organelles? Decreasing the reverse transfer efficiency (e.g. by increasing the number of binding sites in the bud) can allow enrichment of fast diffusing particles in the bud but has little effect on the partitioning of slow diffusing particles (Figure 3a, solutions below the red diagonal and to the lower left). For the latter, increased forward transfer due to active transport (such that k _{ f } > k _{ r }) can obviously lead to bud enrichment. Indeed, most organelles in S. cerevisiae are actively transported into the bud by myosin motors acting along actin cables [2],[3]. However, directed transport requires precise monitoring to ensure that some copies of each organelle remain in the mother cell [1]. Although organelle transport appears to be directional (from the mother to the bud), it is noteworthy to consider that directionally is not strictly required. Active transport with frequent random directionality changes would also promote the partitioning equilibration of large organelles, and could be accounted for in our model as an increase in their apparent diffusion constants or, more precisely, as a form of anomalous diffusion (see Methods).
Conclusions
Our theoretical analysis provides a quantitative overview of the different factors that affect the segregation of cellular content during asymmetric cell division in yeast. In contrast to computationally intensive simulations of particular scenarios, the analytical model derived above provides immediate quantitative assessment of partitioning for arbitrary scenarios defined by their underlying physical parameters. Moreover, our analytical model directly shows how these different physical parameters are related to each other. Analysis of the results of this model suggests that distinct strategies are required to achieve a comparable amount of partitioning of macromolecules and organelles with different mobilities (i.e. different diffusion constants and/or free fractions). Consequently, precise measurements of all parameters are clearly needed to understand how the partitioning of each macromolecular species is accomplished, and to determine whether passive diffusion is sufficient to account for experimentally observed partitioning or whether other processes such as active transport or convection should also be considered. Using fluorescent protein fusions, the mobility of labeled macromolecular species can be investigated with photoconversion, photobleaching or fluorescence correlation spectroscopy, which permit investigation of average mobilities on different time scales, whereas single particle tracking can provide access to the mobilities of individual molecules [53]. Disturbances, including mutant strains, are frequently used to experimentally probe the partitioning process. Since most of the parameters governing partitioning are degenerate, following the partitioning of distinct macromolecular species experimentally upon particular disturbances could be useful for determining whether global parameters or particlespecific parameters are impacted by the disturbance.
References
 1.
Fagarasanu A, Mast FD, Knoblach B, Rachubinski RA: Molecular mechanisms of organelle inheritance: lessons from peroxisomes in yeast. Nat Rev Mol Cell Biol. 2010, 11: 644654. 10.1038/nrm2960.
 2.
Fagarasanu A, Rachubinski RA: Orchestrating organelle inheritance in Saccharomyces cerevisiae. Curr Opin Microbiol. 2007, 10: 528538. 10.1016/j.mib.2007.10.002.
 3.
Lowe M, Barr FA: Inheritance and biogenesis of organelles in the secretory pathway. Nat Rev Mol Cell Biol. 2007, 8: 429439. 10.1038/nrm2179.
 4.
Anderson P, Kedersha N: RNA granules. J Cell Biol. 2006, 172: 803808. 10.1083/jcb.200512082.
 5.
Erickson SL, LykkeAndersen J: Cytoplasmic mRNP granules at a glance. J Cell Sci. 2011, 124: 293297. 10.1242/jcs.072140.
 6.
Mao YS, Zhang B, Spector DL: Biogenesis and function of nuclear bodies. Trends Genet. 2011, 27: 295306. 10.1016/j.tig.2011.05.006.
 7.
Narayanaswamy R, Levy M, Tsechansky M, Stovall GM, O'Connell JD, Mirrielees J, Ellington AD, Marcotte EM: Widespread reorganization of metabolic enzymes into reversible assemblies upon nutrient starvation. Proc Natl Acad Sci USA. 2009, 106: 1014710152. 10.1073/pnas.0812771106.
 8.
Noree C, Sato BK, Broyer RM, Wilhelm JE: Identification of novel filamentforming proteins in Saccharomyces cerevisiae and Drosophila melanogaster. J Cell Biol. 2010, 190: 541551. 10.1083/jcb.201003001.
 9.
Laporte D, Salin B, DaignanFornier B, Sagot I: Reversible cytoplasmic localization of the proteasome in quiescent yeast cells. J Cell Biol. 2008, 181: 737745. 10.1083/jcb.200711154.
 10.
Tyedmers J, Mogk A, Bukau B: Cellular strategies for controlling protein aggregation. Nat Rev Mol Cell Biol. 2010, 11: 777788. 10.1038/nrm2993.
 11.
Murray AW, Szostak JW: Pedigree analysis of plasmid segregation in yeast. Cell. 1983, 34: 961970. 10.1016/00928674(83)905536.
 12.
Sinclair DA, Guarente L: Extrachromosomal rDNA circlesa cause of aging in yeast. Cell. 1997, 91: 10331042. 10.1016/S00928674(00)804936.
 13.
Aguilaniu H, Gustafsson L, Rigoulet M, Nyström T: Asymmetric inheritance of oxidatively damaged proteins during cytokinesis. Science. 2003, 299: 17511753. 10.1126/science.1080418.
 14.
Erjavec N, Larsson L, Grantham J, Nyström T: Accelerated aging and failure to segregate damaged proteins in Sir2 mutants can be suppressed by overproducing the protein aggregationremodeling factor Hsp104p. Genes Dev. 2007, 21: 24102421. 10.1101/gad.439307.
 15.
Liu B, Larsson L, Caballero A, Hao X, Oling D, Grantham J, Nyström T: The polarisome is required for segregation and retrograde transport of protein aggregates. Cell. 2010, 140: 257267. 10.1016/j.cell.2009.12.031.
 16.
Tessarz P, Schwarz M, Mogk A, Bukau B: The yeast AAA + chaperone Hsp104 is part of a network that links the actin cytoskeleton with the inheritance of damaged proteins. Mol Cell Biol. 2009, 29: 37383745. 10.1128/MCB.0020109.
 17.
Liu B, Larsson L, Franssens V, Hao X, Hill SM, Andersson V, Höglund D, Song J, Yang X, Oling D, Grantham J, Winderickx J, Nyström T: Segregation of protein aggregates involves actin and the polarity machinery. Cell. 2011, 147: 959961. 10.1016/j.cell.2011.11.018.
 18.
Zhou C, Slaughter BD, Unruh JR, Eldakak A, Rubinstein B, Li R: Motility and segregation of Hsp104associated protein aggregates in budding yeast. Cell. 2011, 147: 11861196. 10.1016/j.cell.2011.11.002.
 19.
Kaganovich D, Kopito R, Frydman J: Misfolded proteins partition between two distinct quality control compartments. Nature. 2008, 454: 10881095. 10.1038/nature07195.
 20.
Spokoini R, Moldavski O, Nahmias Y, England JL, Schuldiner M, Kaganovich D: Confinement to organelleassociated inclusion structures mediates asymmetric inheritance of aggregated protein in budding yeast. Cell Rep. 2012, 2: 738747. 10.1016/j.celrep.2012.08.024.
 21.
Shcheprova Z, Baldi S, Frei SB, Gonnet G, Barral Y: A mechanism for asymmetric segregation of age during yeast budding. Nature. 2008, 454: 728734.
 22.
Khmelinskii A, Keller PJ, Lorenz H, Schiebel E, Knop M: Segregation of yeast nuclear pores. Nature. 2010, 466: E110.1038/nature09255.
 23.
Makio T, Lapetina DL, Wozniak RW: Inheritance of yeast nuclear pore complexes requires the Nsp1p subcomplex. J Cell Biol. 2013, 203: 187196. 10.1083/jcb.201304047.
 24.
Colombi P, Webster BM, Fröhlich F, Lusk CP: The transmission of nuclear pore complexes to daughter cells requires a cytoplasmic pool of Nsp1. J Cell Biol. 2013, 203: 215232. 10.1083/jcb.201305115.
 25.
Winey M, Yarar D, Giddings TH, Mastronarde DN: Nuclear pore complex number and distribution throughout the Saccharomyces cerevisiae cell cycle by threedimensional reconstruction from electron micrographs of nuclear envelopes. Mol Biol Cell. 1997, 8: 21192132. 10.1091/mbc.8.11.2119.
 26.
MenendezBenito V, van Deventer SJ, JimenezGarcia V, RoyLuzarraga M, van Leeuwen F, Neefjes J: Spatiotemporal analysis of organelle and macromolecular complex inheritance. Proc Natl Acad Sci USA. 2013, 110: 175180. 10.1073/pnas.1207424110.
 27.
Khmelinskii A, Meurer M, Knop M, Schiebel E: Artificial tethering to nuclear pores promotes partitioning of extrachromosomal DNA during yeast asymmetric cell division. Curr Biol. 2011, 21: R17R18. 10.1016/j.cub.2010.11.034.
 28.
Gehlen LR, Nagai S, Shimada K, Meister P, Taddei A, Gasser SM: Nuclear geometry and rapid mitosis ensure asymmetric episome segregation in yeast. Curr Biol. 2011, 21: 2533. 10.1016/j.cub.2010.12.016.
 29.
Schuss Z, Singer A, Holcman D: The narrow escape problem for diffusion in cellular microdomains. Proc Natl Acad Sci USA. 2007, 104: 1609816103. 10.1073/pnas.0706599104.
 30.
Bénichou O, Voituriez R: Narrowescape time problem: time needed for a particle to exit a confining domain through a small window. Phys Rev Lett. 2008, 100: 16810510.1103/PhysRevLett.100.168105.
 31.
Cheviakov AF, Reimer AS, Ward MJ: Mathematical modeling and numerical computation of narrow escape problems. Phys Rev E Stat Nonlin Soft Matter Phys. 2012, 85: 02113110.1103/PhysRevE.85.021131.
 32.
Rupprecht JF, Bénichou O, Grebenkov DS, Voituriez R: Exact mean exit time for surfacemediated diffusion. Phys Rev E Stat Nonlin Soft Matter Phys. 2012, 86: 04113510.1103/PhysRevE.86.041135.
 33.
Schuss Z: Equilibrium and recrossings of the transition state: what can be learned from diffusion?. J Phys Chem C. 2010, 114: 2032020334. 10.1021/jp101674y.
 34.
Holcman D, Schuss Z: Control of flux by narrow passages and hidden targets in cellular biology. Rep Prog Phys. 2013, 76: 07460110.1088/00344885/76/7/074601.
 35.
Barkai E, Garini Y, Metzler R: Strange kinetics of single molecules in living cells. Physics Today. 2012, 65: 2935. 10.1063/PT.3.1677.
 36.
Höfling F, Franosch T: Anomalous transport in the crowded world of biological cells. Rep Prog Phys. 2013, 76: 04660210.1088/00344885/76/4/046602.
 37.
Howard J, Grill SW, Bois JS: Turing's next steps: the mechanochemical basis of morphogenesis. Nat Rev Mol Cell Biol. 2011, 12: 392398. 10.1038/nrm3120.
 38.
Goehring NW, Grill SW: Cell polarity: mechanochemical patterning. Trends Cell Biol. 2013, 23: 7280. 10.1016/j.tcb.2012.10.009.
 39.
Gutzeit H, Koppa R: Timelapse film analysis of cytoplasmic streaming during late oogenesis of Drosophila. J Embryol Exp Morphol. 1982, 67: 101111.
 40.
Goldstein RE, Tuval I, van de Meent JW: Microfluidics of cytoplasmic streaming and its implications for intracellular transport. Proc Natl Acad Sci USA. 2008, 105: 36633667. 10.1073/pnas.0707223105.
 41.
Mattos TG, MejíaMonasterio C, Metzler R, Oshanin G: First passages in bounded domains: when is the mean first passage time meaningful?. Phys Rev E Stat Nonlin Soft Matter Phys. 2012, 86: 03114310.1103/PhysRevE.86.031143.
 42.
Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81: 23402361. 10.1021/j100540a008.
 43.
Slaughter BD, Schwartz JW, Li R: Mapping dynamic protein interactions in MAP kinase signaling using livecell fluorescence fluctuation spectroscopy and imaging. Proc Natl Acad Sci USA. 2007, 104: 2032020325. 10.1073/pnas.0710336105.
 44.
LubyPhelps K: Probing the structure of cytoplasm. J Cell Biol. 1986, 102: 20152022. 10.1083/jcb.102.6.2015.
 45.
LubyPhelps K, Castle PE, Taylor DL, Lanni F: Hindered diffusion of inert tracer particles in the cytoplasm of mouse 3 T3 cells. Proc Natl Acad Sci USA. 1987, 84: 49104913. 10.1073/pnas.84.14.4910.
 46.
Lukacs GL, Haggie P, Seksek O, Lechardeur D, Freedman N, Verkman AS: Sizedependent DNA mobility in cytoplasm and nucleus. J Biol Chem. 2000, 275: 16251629. 10.1074/jbc.275.3.1625.
 47.
EscusaToret S, Vonk WIM, Frydman J: Spatial sequestration of misfolded proteins by a dynamic chaperone pathway enhances cellular fitness during stress. Nat Cell Biol. 2013, 15: 12311243. 10.1038/ncb2838.
 48.
Crow ET, Li L: Newly identified prions in budding yeast, and their possible functions. Semin Cell Dev Biol. 2011, 22: 452459. 10.1016/j.semcdb.2011.03.003.
 49.
Winkler J, Tyedmers J, Bukau B, Mogk A: Chaperone networks in protein disaggregation and prion propagation. J Struct Biol. 2012, 179: 152160. 10.1016/j.jsb.2012.05.002.
 50.
Verges KJ, Smith MH, Toyama BH, Weissman JS: Strain conformation, primary structure and the propagation of the yeast prion [PSI+]. Nat Struct Mol Biol. 2011, 18: 493499. 10.1038/nsmb.2030.
 51.
Derdowski A, Sindi SS, Klaips CL, DiSalvo S, Serio TR: A size threshold limits prion transmission and establishes phenotypic diversity. Science. 2010, 330: 680683. 10.1126/science.1197785.
 52.
KawaiNoma S, Pack CG, Kojidani T, Asakawa H, Hiraoka Y, Kinjo M, Haraguchi T, Taguchi H, Hirata A: In vivo evidence for the fibrillar structures of Sup35 prions in yeast cells. J Cell Biol. 2010, 190: 223231. 10.1083/jcb.201002149.
 53.
Miyawaki A: Proteins on the move: insights gained from fluorescent protein technologies. Nat Rev Mol Cell Biol. 2011, 12: 656668. 10.1038/nrm3199.
Acknowledgements
We thank Gilles Charvin, Joseph D. Barry and Elmar Schiebel for discussions and comments on the manuscript. We thank the reviewers for helpful comments and suggestions. This work was supported in parts by a grant from the German Research Foundation (DFG), Sonderforschungsbereich (SFB) 1036 (TP10).
Author information
Affiliations
Corresponding authors
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
All authors conceived the research. AKi developed the model with input from AKh and MK. All authors discussed the results and wrote the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Kinkhabwala, A., Khmelinskii, A. & Knop, M. Analytical model for macromolecular partitioning during yeast cell division. BMC Biophys 7, 10 (2014). https://doi.org/10.1186/s1362801400106
Received:
Accepted:
Published:
Keywords
 Asymmetric cell division
 Yeast mitosis
 Cellular transport processes
 Organelle segregation
 Organelle inheritance
 Diffusion
 Narrow escape
 Bud neck
 Protein aggregate
 Prion
 Extrachromosomal rDNA circle
 ERC
 Nuclear pore complex
 NPC