Length regulation of multiple flagella that self-assemble from a shared pool of components
Abstract
The single-celled green algae Chlamydomonas reinhardtii with its two flagella—microtubule-based structures of equal and constant lengths—is the canonical model organism for studying size control of organelles. Experiments have identified motor-driven transport of tubulin to the flagella tips as a key component of their length control. Here we consider a class of models whose key assumption is that proteins responsible for the intraflagellar transport (IFT) of tubulin are present in limiting amounts. We show that the limiting-pool assumption is insufficient to describe the results of severing experiments, in which a flagellum is regenerated after it has been severed. Next, we consider an extension of the limiting-pool model that incorporates proteins that depolymerize microtubules. We show that this ‘active disassembly’ model of flagellar length control explains in quantitative detail the results of severing experiments and use it to make predictions that can be tested in experiments.
https://doi.org/10.7554/eLife.42599.001Introduction
The size regulation of cellular organelles is a fundamental problem in biology (Marshall, 2016; Milo and Phillips, 2015). For example, nuclear size is tightly coupled with cell size across a wide range of species (Hara and Merten, 2015) and the loss of this coupling mechanism is implicated in various types of cancer (Zink et al., 2004).
A striking example of organelle size control in eukaryotes is the single-celled algae Chlamydomonas reinhardtii (Figure 1), which uses two flagella to move through its aqueous environment. The backbone of each flagellum is an assembly known as the axoneme that consists of nine microtubule doublets arranged in a ring around a central pair of microtubules (Fawcett and Porter, 1954; Witman et al., 1972). Unlike the dynamic instability of cytoplasmic microtubules, which can alternate between rapidly shortening ‘catastrophe’ and stable ‘rescue’ states depending on whether or not the tip is bound to GTP, microtubules in the axoneme exist in a highly stable state (Behnke and Forer, 1967; Orbach and Howard, 2019). This stability reflects a tight control over flagellar lengths, the loss of which has dramatic physiological consequences; mutants with longer flagella have decreased swimming velocities and beat frequencies (Khona et al., 2013) compared to wild type cells and mutants with unequal flagellar lengths are observed to spin around in circles (Tam et al., 2003).
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-fig1-v2.tif/full/617,/0/default.jpg)
Experimental background.
(a) Electron microscopy images of the biflagellate green algae Chlamydomonas and its flagella captured by Elisa Vannuccini and Pietro Lupetti (University of Siena, Italy) and reproduced from Morga and Bastin (2013) under the Creative Commons Attribution License CC BY 2.0 [http://creativecommons.org/licenses/by/2.0]). The inset shows the whole organism (scale bar 5 µm) and the close-up shows the flagellar basal body (BB), transition zone (TZ), and cell wall (CW) (scale bar 1 µm). (b) Severing experiments: after one flagellum is severed, the two flagella equalize at a shorter length and then grow together to the original steady-state length. This is demonstrated in experimental data of 20 severing experiments from Ludington et al. (2012) provided by the authors. The green and blue shaded regions show the mean plus or minus one standard deviation.
A key process contributing to the assembly of flagella is the continual transport of proteins from the flagellar base to tip and back. The original evidence for this intraflagellar transport (IFT) was provided around 25 years ago by experimental observations in Chlamydomonas of particles moving processively along the flagellum at constant speed (Kozminski et al., 1993). In the time since, a significant body of work has revealed the many proteins and biochemical pathways that coordinate this complex process in Chlamydomonas and other organisms such as C. elegans, as described in several review articles (Prevo et al., 2017; Scholey, 2003; Scholey, 2008; Cole, 2003; Rosenbaum et al., 1999; Rosenbaum and Witman, 2002; Rosenbaum et al., 1999).
Already at the time of its discovery, IFT was hypothesized to play a role in flagellar length control by transporting building blocks to the tip of the flagellum (Kozminski et al., 1993). IFT particles containing tubulin are transported along the flagellum by two different motor proteins: kinesin-2 transports IFT particles from the flagellar base to tip (the anterograde direction) whereas dynein carries IFT particles from the tip to the base (the retrograde direction). As shown in subsequent work (Song and Dentler, 2001; Marshall and Rosenbaum, 2001; Buisson et al., 2013), these flagellar proteins are continually exchanged with those localized to the basal body of the flagellum, represented schematically in Figure 2a as a basal pool.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-fig2-v2.tif/full/617,/0/default.jpg)
IFT in Chlamydomonas with diffusive return of kinesin-2 to the base.
(a) Kinesin-2, dynein, and tubulin combine to form a complex with other IFT components in the basal pool and are injected into the flagellum. The kinesin-2 motors move toward the tip of the flagellum, where the complex eventually breaks down into its constituent parts. Dynein motors move in the retrograde direction carrying back some, but not all IFT components. Specifically, kinesin motors are not part of the retrograde IFT complex and they diffuse back to the base. In steady-state, the injection flux , ballistic flux , and diffusive flux are all equal. (b) IFT flux vs. length, comparison between the quasi-steady state approximation (Equation 4) and stochastic simulations (Appendix 1). The shaded red area represents the mean and standard deviation computed from 10 realizations of the stochastic simulations. (inset) Kinesin-2 concentration vs. position along the flagellum (steady-state approximation [Equation 31]).
These observations motivated the development of mathematical models of flagellar length dynamics such as the balance point model (Marshall and Rosenbaum, 2001; Marshall et al., 2005). In the balance point model, there is a continual competition between assembly and disassembly, and either the rate of assembly, the rate of disassembly, or both, may be length-dependent. The steady-state length is determined by the point at which the assembly and disassembly processes come into balance.
Up to now, the balance point model has been considered primarily for a single flagellum and the focus has been on assembly being the length dependent process that leads to length control (Marshall and Rosenbaum, 2001). However, other formulations of the balance point model, such as the case in which the disassembly rate is also length-dependent, have received considerably less attention. In this work, we revisit the balance point model and particularly the assumption of a constant rate of disassembly in the context of simultaneous length control of multiple flagella assembled from a shared pool of biomolecules.
We limit our theoretical exploration to the space of models defined by the following processes: IFT particle assembly and injection at the flagellar base, motion of IFT proteins along the flagellum, and tubulin polymerization and depolymerization at the flagellar tip (see schematic Figure 2a). We further assume that IFT particle injection satisfies first-order chemical kinetics and allow for a control mechanism that regulates protein levels in the basal pool. Note that this model space does not include all possibilities. In particular, it does not include the time-of-flight model (Wren et al., 2013), in which additional reactions affect protein state inside the flagellum (e.g. proteins enter in an activated state and deactivate at some rate).
Nevertheless, our work retains a high degree of generality. For example, it allows for different modes of coupling between proteins in basal pools (shared or separate), different modes of IFT motion (ballistic or diffusive and processive or non-processive), and general depolymerizing activity. Within the model space outlined, our main results hold independent of these details. Notably, we find that length-independent disassembly of microtubules cannot account for the experimental results, whereas incorporating length-dependent disassembly (e.g. through the ballistic-to-diffusive motion of a depolymerizing protein) leads to reasonable agreement with the experiments.
We analyze these models at two levels of detail, both using detailed agent-based stochastic simulations and a reduced description in terms of ordinary differential equations (ODEs) and compare the results of simulations to experimental data. It is known from experiments on Chlamydomonas that (i) its two flagella reach steady-state lengths of about 10 µm, and that (ii) the flagellar lengths are correlated since, if one flagellum is severed, the remaining flagellum shortens until it reaches the length of the growing, previously severed flagellum (Rosenbaum et al., 1969). This latter protocol is typically referred to as a ‘severing’ or ‘long-zero’ experiment.
As shown in Figure 1b, the severing experiments contain two separate timescales. There is an initial fast timescale of 10–20 min over which the flagellar lengths first equalize, followed by a slower timescale of 60–90 min over which the lengths increase simultaneously to their original lengths. With respect to the short time behavior, we find that none of the candidate models assuming constant disassembly capture the rapid length equalization observed experimentally. Moreover, if an external control mechanism is added that replenishes proteins lost due to severing over a slower timescale, the constant disassembly models within our model space fail to control lengths at all.
Motivated by experimental observations of microtubule-depolymerizing proteins within the flagellum (Piao et al., 2009; Luo et al., 2011; Hilton et al., 2013), we subsequently consider a model in which the disassembly rate is dependent on the local concentration of a depolymerizing protein at the flagellar tips. This gives rise to a length-dependent disassembly rate, which is a departure from existing models that assume constant disassembly. We find that this model is consistent with the rapid length equalization observed in experiments. Further, upon adding a control mechanism on protein levels in the basal protein pool, the model is able to capture the slow return to the original steady-state lengths observed experimentally. In the Discussion, we examine this model in light of the currently available experimental data and discuss possible candidates for the depolymerizing protein.
Results
Limiting motor pool gives rise to length-dependent injection rate
It is known from experiments that the injection rate of IFT particles in Chlamydomonas is length-dependent (Dentler, 2005). In this section, we show how length-dependent injection is a consequence of mass action kinetics of proteins available in limiting amounts.
We first consider the flagellar length dynamics of a single flagellum. Our approach initially closely follows Hendel et al. (2018), although important differences will appear later on in the case of two flagella. The key biochemical variables are the tubulin dimers that make up the flagellar axoneme and the molecular motors kinesin-2 and dynein that transport IFT particles from the base to the tip and back. IFT particles combine in the basal pool with kinesin-2, tubulin and dynein to form a complex that is injected into the flagellum (Cole et al., 1998); see the schematic in Figure 2a.
Here we consider the case that the rate-limiting molecule is kinesin-2, although similar results would apply to any other protein being rate-limiting for IFT assembly. Denoting the number of free molecular motors in the pool at the flagellar base by , the injection flux of IFT particles into the flagellum satisfies
Note that we consider the cell volume to be fixed, in which case biomolecule numbers and concentrations may be used interchangeably.
First we consider the case in which the total number of motors is conserved. Later on we will also consider the case in which motor concentrations in basal pools are regulated by an external control mechanism, motivated by the severing experiments described previously in which length recovery indicates replenishment of protein levels.
Kinesin-2 has been shown to undergo ballistic transport in the anterograde direction and diffusive motion in the retrograde direction (Chien et al., 2017). The total number of motors satisfies , where is the number of motors moving ballistically on the flagellum in IFT particles and is the number of motors moving diffusively. Therefore we may rewrite Equation 1 as
The intraflagellar dynamics are fast compared to changes in length. The timescale of flagellar length dynamics, for example the recovery time after severing, is of the order of 10 min, whereas the molecular motors involved in IFT take at most tens of seconds to traverse the length of the flagellum by either ballistic or diffusive motion. Based on this separation of timescales, we treat IFT as a quasi-steady state process in which the injection flux , diffusive flux , and ballistic flux are balanced, that is
As shown in Materials and methods, the quasi-steady state assumption of flux balance together with mass action kinetics of injection expressed in Equation 1 imply that
where is the rate constant of motor injection, is the total number of motors, is the ballistic motor speed in IFT, and is the diffusion coefficient of motors in the flagellum.
We have validated the quasi-steady state assumption used to derive Equation 4 by comparing the IFT particle flux and the concentration of diffusing motors in the flagellum to the results of stochastic simulations (see Appendix 1), in which the dynamics of motors as well as the dynamics of microtubule assembly are taken into account explicitly; see Figure 2b. Parameter values used in simulations are provided in Table 1. (Note that in Equation 4 and throughout the manuscript, terms written in the form are shorthand for , for example is to be read as .)
Parameter values and definitions.
https://doi.org/10.7554/eLife.42599.004Symbol | Definition | Value | Units | References |
---|---|---|---|---|
Parameters | ||||
Steady-state length | 10–12 | µm | Marshall and Rosenbaum, 2001; Rosenbaum et al., 1969 | |
Tubulin pool per flagellum | 38–47 | µm | Marshall et al., 2005 | |
Disassembly speed | 0.5 | µm/min | Marshall and Rosenbaum, 2001; Ludington et al., 2012 | |
IFT speed | 2.5–3 | µm/s | Kozminski et al., 1993; Buisson et al., 2013 | |
Diffusion coefficient | 1.7 | µm2/s | Chien et al., 2017 | |
Assembly rate per tubulin | 2.3 × 10-2 –3.6 × 10-2 | min−1 | Fit | |
Injection rate constant | 0.8–4 | min−1 | Fit | |
Prefactor in Equation 5 | 2.5 × 10-4 | Estimate (Appendix 2) | ||
Variables | ||||
Number of flagella | ||||
Free tubulin | µm | |||
Total motors | ||||
Free motors | ||||
Ballistic motors | ||||
Diffusing motors | ||||
Flux | min−1 | |||
Motor concentration | µm−1 | |||
Average concentration | µm−1 |
IFT particle injection arising from a finite number of motors shared between the flagellum and the basal pool therefore leads to a length-dependent flux. This result holds true regardless of the identity of the rate-limiting IFT protein. However, the scaling of formula Equation 4 with length depends on whether ballistic or diffusive transport dominates. In the limit , the scaling of Marshall and Rosenbaum (2001) is recovered whereas in the limit , we recover the scaling of Hendel et al. (2018). Note that a distinction between our model and these previous works is the presence of a constant term in the denominator of Equation 4, which implies that in our formulation the flux does not blow up at .
In our model the assembly rate is determined by the rate of tubulin transport to the flagellar tip, as in Marshall and Rosenbaum (2001) and Hendel et al. (2018). Given mass-action kinetics of IFT particle assembly in the basal pool, this is simply the flux of IFT particles times the amount of free tubulin , so that the growth rate is given by the following ODE:
where is the disassembly speed (assumed constant for now) and is a constant. The total amount of tubulin is assumed to be conserved for now so that , where tubulin is measured in units of corresponding flagellar length. As previously mentioned, later on we will consider the case in which protein levels in basal pools are not conserved and are instead monitored and regulated by an external control mechanism.
Substituting the expression Equation 4 for the flux into the above growth rate results in
As shown in Materials and methods, this equation yields a stable steady-state length for the single flagellum.
Limiting-pool mechanisms alone cannot account for the rapid length equalization observed in severing experiments
Whereas for a single flagellum the limiting-pool mechanism leads to a stable steady-state length, we show next that this mechanism is insufficient for the simultaneous length control of two flagella. A general limiting-pool model for the dynamics of two flagella having lengths and is given by the following ODEs:
where and for denote the fluxes of IFT particles into the two flagella and the amounts of free tubulin in their basal pools, respectively. As we will show next, the particular forms of and depend on how the pools are coupled. In particular, the fluxes and free tubulin will be equal for the two flagella if the basal proteins are held in a common shared pool.
Severing experiments illustrate that the flagella are coupled. We consider various modes of coupling, that is shared or separate motor pools and shared or separate tubulin pools, as depicted in Figure 3a, that give rise to different forms of and within our model. We investigate their consequences for length control by focusing on the solutions to the steady-state equations
Length control implies that these steady-state equations must yield a unique steady-state solution for and . We consider the steady-state lengths before and after severing. It is known from experiments (Figure 1b) that initially the two flagella have equal lengths, and that after severing, which is accompanied by a loss of material (e.g. tubulin and motors lost from the severed flagellum), there is a rapid equalization of flagellar lengths. This initial, fast equalization of lengths leads to flagella that are shorter than they were before the severing. These experimental observations may be used to reject candidate models. Here we focus on the short-time dynamics after severing and therefore do not account for protein replenishment; however, as we show in the next section, incorporating protein replenishment is incompatible with the constant disassembly models considered here.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-fig3-v2.tif/full/617,/0/default.jpg)
Length dynamics of two flagella assembling from shared pools of building blocks.
(a) Flagellar assembly in Chlamydomonas reinhardtii and modes of coupling between basal proteins pools. (b) Simulations of the severing experiment using different modes of coupling between basal pools of proteins. In all cases severing occurs at time zero. (i) In the case of separate pools of both motors and tubulin (Equations 44 and 45), the unsevered flagellum does not decrease in length. (ii) When only motors are shared (Equations 49 and 50), the flagellar lengths do not equalize after severing. (iii) When only tubulin is shared (Equations 56 and 57), the flagellar lengths do not equalize after severing. (iv) When both tubulin and motors are shared (Equations 58 and 59), there is no unique steady-state. The difference between the two flagellar lengths undergoes a random walk, as shown by the mean square displacement (inset, average over 10 runs).
As we shall see, the constant disassembly models within our model space do not yield the rapid length equalization observed experimentally. As shown in Figure 3b and considered next on a case-by-case basis, if the lengths are in steady-state prior to severing, they are also in steady-state immediately after severing as well so that no length equalization occurs. This is because in these models the assembly and disassembly rates depend only on free protein levels in the basal protein pools, which are unaffected by severing. To be more precise, we may rewrite the governing equations Equations 7 and 8 in terms of the number of free motors and in the following form:
In order to simulate a severing experiment, we first allow the two flagella to reach steady-state. Prior to severing, there is an amount of tubulin and a number of motors loaded on the first flagellum corresponding to the numbers of motors undergoing ballistic or diffusive motion (see Materials and methods). Severing is modeled by setting , , and . Because severing does not change any of the protein levels in the basal pools, and for are the same before and after severing. Moreover, the linear concentration profile of diffusing motors remains linear in the truncated flagellum; no redistribution of motors is required to maintain quasi-steady state. As a consequence the constant disassembly models predict no length recovery of the severed flagellum, contrary to what is observed.
Tubulin separate, motors separate
The case of separate tubulin pools and individual motor pools leads to two uncoupled instances of (6). This is fundamentally inconsistent with the coupling observed in severing experiments, in particular the significant decrease in the length of the unsevered flagellum (Figure 3b(i)). In this model, the unsevered flagellum does not change length after severing, and therefore it may be ruled out.
Tubulin separate, motors shared
When motors are shared through a common pool, IFT particles are injected into either flagellum with equal probability. Therefore and it can be shown by a straightforward generalization of Equation 4 (see Materials and methods) that the flux satisfies
Initially the tubulin pools are equal, as are the flagellar lengths. However, after the loss of material due to severing, and the lengths do not equalize (see Figure 3b(ii) and Figure 3—video 1). Because of the separate tubulin pools, severing leads to asymmetrical tubulin depletion and unequal steady-state lengths after severing. Therefore, we may rule out the model.
Tubulin shared, motors separate
This case is analogous to the previous model in which only motors are shared, but this time only the tubulin pools are shared so that , with
As shown in Figure 3b(iii) and Figure 3—video 2, this model does not capture the length equalization observed in severing experiments and therefore may be ruled out. See Materials and methods for details.
Tubulin shared, motors shared
The flux resulting from the shared motor assumption is the same as in Equation 13, and we are left with the steady-state equations
Because the two steady-state equations are identical, that is , this model does not account for the simultaneous positive and negative growth rates for the two flagella observed in the severing experiment.
More strikingly, subtracting the steady-state equations yields
so that the difference in lengths is not controlled at all. Note that this result is independent of the parameters. Indeed, a similar conclusion was reached in the context of actin filaments (Mohapatra et al., 2017), in which it was observed that sharing all biomolecules between filaments does not yield simultaneous length control. In the context of the full stochastic simulations, this degeneracy is manifested by the difference in lengths undergoing a random walk (inset to Figure 3b(iv) and Figure 3—video 3).
The above analysis shows that, regardless of the manner in which tubulin and motors are shared between flagella, the constant disassembly models we have considered are unable to explain the results of severing experiments. Although sharing either tubulin or motors, but not both, yields a unique steady-state, these models do not agree with the rapid length equalization observed in severing experiments. This motivates us to extend our study beyond the models considered thus far.
Controlling protein levels in the basal pool is incompatible with the constant disassembly models considered
So far, we have assumed that the tubulin pool and motor pool are fixed throughout the simulation. While this assumption is reasonable for the fast initial phase of the severing experiment in which the flagellar lengths rapidly equalize, the slower second phase of recovery to the original steady state lengths requires replenishment of proteins back to their original levels. This was shown experimentally by using cycloheximide at the time of severing to block the synthesis of new proteins (Rosenbaum et al., 1969), resulting in shorter flagella that did not recover to their original lengths.
It may seem plausible that adding an external control mechanism that replenishes protein levels would lead to length equalization, thus resolving the issue of unequal steady-state lengths after severing. However, as we next show, adding such a control mechanism on free proteins levels in the basal pool does not lead to length equalization. Instead, in this case the constant disassembly models we have considered completely fail to control lengths.
To incorporate control on protein levels into our model, we assume that the total levels and of tubulin and motors are replenished over a timescale as the cell synthesizes new protein to achieve target protein levels and in the basal pool:
In steady-state, the above equations become and . The steady-state equation for length then becomes , that is length drops out of the assembly term completely! Therefore, this external control mechanism actually destabilizes the flagellar lengths (see Figure 3—videos 4—6 to observe this destabilization for various modes of coupling). Note that although we have used the simplest case of linear feedback Equations 18 and 19 to illustrate the point, this argument is general and does not depend on the details of the control mechanism. The only requirement is that the free protein levels and are driven to their target values and in steady state. This argument gives another compelling reason to look beyond the constant disassembly models considered thus far.
Tubulin shared, motors shared and concentration-dependent disassembly
We next consider a model that allows for full exchange of IFT components between basal protein pools and replaces the constant disassembly assumption with a concentration-dependent disassembly rate. The assumption of a constant disassembly rate was based on experiments on mutants in which IFT was disabled (Marshall et al., 2005). However, subsequent experiments in organisms with intact IFT led to 50-fold greater disassembly rates than those measured in the absence of IFT (Ludington et al., 2012).
Experimental observations that some kinesin species (e.g. kinesin-13) participate in microtubule disassembly (Piao et al., 2009) provide a potential biochemical basis for IFT-dependent disassembly. In what follows we take the disassembly rate to depend on the concentration of a depolymerizing protein. This is a reasonable model for a depolymerizer that is non-processive in its depolymerization activity in that it removes at most a few tubulin subunits before falling off into a deactivated state. (The case of processive depolymerizers is investigated in Appendix 3).
We will assume that the depolymerizer has the same motion as kinesin-2—uninterrupted ballistic motion to the tip followed by diffusive motion to the base—resulting in a linear concentration profile. This would be the case for any non-motile protein that is transported ballistically to the flagellar tip as IFT cargo and diffuses back to the flagellar base. Note however that the ballistic-to-diffusive assumption is not essential for the model; so long as there is a gradient in the depolymerizer concentration, the conclusion of simultaneous length control holds. The essential ingredient is the non-constant concentration along the length of the flagellum, which in this case is achieved by ballistic anterograde motion and diffusive return.
Given that the formulas for the steady-state flux apply for any rate-limiting IFT protein undergoing ballistic-to-diffusive motion along the flagella, for convenience we assume in what follows that the depolymerizer is the rate-limiting protein. However, this assumption is made only for convenience; in the more general case that the depolymerizer and the rate-limiting IFT protein are different, the same results are obtained with suitably rescaled parameters, as shown in Appendix 3.
We replace the assumption of constant disassembly by a disassembly speed of the form , where is the concentration of diffusing motors at the tip of the flagellum and . In general the disassembly rate may be an arbitrary function of concentration, in which case this model may be viewed as a first-order Taylor series expansion valid near steady-state. The flux and concentration at the flagellar tip are related by (Equation 31 in Materials and methods and Figure 2b (inset)), so that we may rewrite the governing equations as
where as before in the case of shared motors
This model yields simultaneous length control and length equalization after severing (Figure 4 and Figure 4—video 1). Subtracting Equation 21 from Equation 20, it follows immediately that , and solving for the steady-state length results in
In Appendix 3 we show that this solution is stable using linear stability analysis. Therefore, concentration-dependent disassembly yields simultaneous length control when all biomolecules are shared between flagella and is consistent with the rapid length equalization observed after severing Figure 4b. The presence of a concentration gradient is a critical ingredient in this model and here it is achieved by ballistic transport to the flagellar tip with diffusive return. The concentration gradient makes the disassembly rates length-dependent and yields independent equations for the steady-state lengths.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-fig4-v2.tif/full/617,/0/default.jpg)
Concentration-dependent disassembly model: simultaneous length control is achieved using shared tubulin and shared depolymerizers.
(a) The depolymerizer moves ballistically to the flagellar tip and diffuses back, (b) The model Equations 20 and 21 captures rapid length equalization, (c) Protein replenishment with timescale min is included through Equations 18 and 19 and the model is fit to severing data from Ludington et al. (2012) (and to Rosenbaum et al., 1969, see inset). For the fit to Ludington et al. (2012), we use fitting parameters of , , , , and ; all other parameters are as in Table 1. For the fit to Rosenbaum et al., 1969, we use the same fitting parameters except for , , and .
Unlike the constant disassembly models we have considered, for which a limiting-pool mechanism is essential for length control, concentration-dependent disassembly yields length control under mild assumptions including the case that all biomolecules are in excess. Nevertheless, limiting-pools of biomolecules are needed to capture the depletion effects observed in severing experiments on Chlamydomonas, for example the shortening of the unsevered flagellum and the previously-mentioned absence of length recovery after cyclohexamide treatment (Rosenbaum et al., 1969). In Appendix 3 we explain that limiting pools of IFT motors are necessary for agreement with data whereas tubulin may either be limited or in excess. The presence of a limiting-pool once again raises the question, now in the context of concentration-dependent disassembly, of how biomolecules may be shared between flagella. In Appendix 3 we show that all relevant biomolecule pools must be shared for the concentration-dependent disassembly model to capture the rapid length equalization observed.
Another feature of the concentration-dependent disassembly model is that it allows for an external control mechanism on protein levels in the basal pool, unlike the constant disassembly models we have considered. As shown in Figure 4c and Figure 4—video 2, upon including protein replenishment via Equations 18 and 19 on a timescale of mins, the recovery of the flagella back to their original lengths is in reasonable agreement with experimental data.
Generalization to flagella
The concentration-dependent disassembly model may be generalized to arbitrary flagellar number , and here we demonstrate simultaneous length control in the case of flagella (Figure 5).
Because motors are shared, the injection fluxes are equal and for all , with satisfying
and length dynamics given by
where we have applied the boundary condition as before. Taking any pairwise difference between the th and th equations at steady-state yields immediately , so that the steady-state lengths are equal to
for all . Stability follows from analyzing the linearized equations, as shown in Appendix 3.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-fig5-v2.tif/full/617,/0/default.jpg)
The concentration-dependent disassembly model generalizes to arbitrary flagellar number .
We solve (Equations 25) with flagella and the larger shared pool , , and ; otherwise all parameters are as in Table 1.
Discussion
In this work, we capture the aspects of IFT essential for length control, that is the motor-driven transport of tubulin across the flagellum, to explore models of flagellar length dynamics. Our theoretical framework makes it possible to investigate the consequences of biomolecule exchange between flagella on length control. In our initial exploration, in which we take the disassembly rate to be constant, we find that sharing both tubulin and motors leads to an indeterminate system of equations regardless of the details of the model, whereas sharing either tubulin or motors, but not both, results in simultaneous length control of both flagella. However, by examining the steady-state lengths immediately before and after severing and accounting for depletion in both the tubulin and motor pools, we observe that none of the constant disassembly models we have considered are able to capture the length equalization observed in experiments.
Given that the constant disassembly models we have considered are unable to explain the experiments, we have proposed a model in which disassembly depends on the local concentration of a depolymerizing protein at the tip of the flagellum (Figure 4). In this ‘active disassembly’ model, a length-dependent concentration at the tip is achieved by assuming that the depolymerizer undergoes ballistic motion to the tip and returns diffusively. This model agrees with the results of severing experiments. As suggested by the title of Hendel et al. (2018), the linear concentration gradient generated by diffusion acts as a ruler. However, in our model diffusion must be combined with a mechanism such as concentration-dependent disassembly; the constant disassembly models we have considered are inconsistent with severing experiments regardless of whether the motor dynamics are ballistic, diffusive, or some combination.
We remark on the differences between our model and Hendel et al. (2018), in which simultaneous length control was obtained using a balance point model with constant disassembly. As we describe in Materials and methods, the formulation of Hendel et al. (2018) is very similar to Equations 56 and 57 derived in the case of shared tubulin pools and separate motor pools with constant disassembly. Whereas these equations do not lead to length equalization after severing in the case of no protein replenishment (Figure 3b(iii)), the model of Hendel et al. (2018) achieves length equalization by replenishing the total number of motors on the flagellum, that is through an additional control mechanism that instantaneously adds the motors lost through severing back into the basal pool. The importance of replenishment for the model appears to be inconsistent with severing experiments that show length equalization occurs even when protein synthesis is blocked using cyclohexamide (Rosenbaum et al., 1969). In contrast, the concentration-dependent disassembly model achieves length equalization with or without replenishment (Figure 4b and Figure 4c). Note further that controlling motor number in the flagellum is not equivalent to controlling protein concentrations in the basal pool. As shown in Results, when basal pool concentrations are controlled according to Equations 18 and 19, there is a breakdown of length control for the constant disassembly models.
Our model predicts that the tubulin and depolymerizer pools must both be shared for the concentration-dependent disassembly model to capture the rapid length equalization observed (Appendix 3). This illustrates the dramatic consequences in behavior that can occur when biomolecules are shared between compartments, and highlights the importance of knowing which proteins are exchanged between the basal pools in the context of flagellar length control. In particular, having a protein that is not exchanged can provide simultaneous length control by a limiting-pool mechanism, but it introduces asymmetries that contradict the length equalization observed in severing experiments.
In addition to our claim that the depolymerization rate is non-constant and dependent on length, our model leads to testable predictions that may be useful in identifying candidate depolymerizers. For example, according to our model the depolymerizer active in length control is not uniformly distributed along the flagellum; its concentration increases toward the flagellar tip. This could be tested experimentally by fluorescently labeling candidate depolymerizers and studying their concentration profiles along the flagellum as recently done to characterize the concentration profile of kinesin-13 in Giardia (McInally et al., 2019).
In our model the essential ingredient that leads to simultaneous length control is the presence of a depolymerizing protein with a concentration gradient along the flagellum. How proteins can develop and maintain such concentration profiles is therefore one of the key questions raised by the model. Indeed, concentration gradients (unassociated with depolymerizing activity) were observed already in Hendel et al. (2018) as the result of ballistic-to-diffusive motion. Although here as proof of principle this concentration gradient is achieved by the mechanism of ballistic-to-diffusive motion, our main results are independent of the detailed form of this concentration gradient and how it is generated. As shown in Appendix 3—figure 3b, our model allows for depolymerizers with nonlinear concentration profiles and different patterns of motion, for example exponential concentration distributions such as those recently observed in Giardia (McInally et al., 2019) and those generated by motile proteins that bind and unbind to cytoskeletal filaments, as theorized in the context of actin-myosin systems (Naoz et al., 2008; Orly et al., 2014; Pinkoviezky and Gov, 2014; Pinkoviezky and Gov, 2017; Yochelis et al., 2015).
Our main results do not rely on many of the details of the model. Although here we have taken depolymerase activity to depend linearly on concentration, the model generalizes to the non-linear case in a straightforward manner. For a local depolymerization rate that is an arbitrary function of concentration, a Taylor series expansion may be performed as in Klein et al. (2005) to obtain the corresponding linearized system discussed in Appendix 3. Further, whereas here we have explored the case of non-motile depolymerizers transported to the flagellar tip by IFT, motile proteins could in principle aggregate at the flagellar tip independent of IFT. Although we are not aware of any such examples in Chlamydomonas, interestingly in budding yeast the motile protein Kip3p has been shown to depolymerize microtubules in a length-dependent manner (Varga et al., 2009). The aggregation of motile depolymerizing proteins has also been demonstrated in previous theoretical studies of microtubule length control (Klein et al., 2005; Johann et al., 2012; Reese et al., 2014); note however that these previous works differed from our model of flagellar IFT in that they considered isolated microtubules surrounded by a constant concentration bath and/or significant steric interactions between motile proteins.
Finally, although in principle concentration-dependent disassembly yields length control even the case that all biomolecules are in excess, in Chlamydomonas severing experiments have shown that depletion effects are important (e.g. cyclohexamide treatment yields short flagella that do not return to their original steady-state lengths [Rosenbaum et al., 1969]). Within our model such depletion effects arise through limiting-pools of proteins, and in Appendix 3 we explain that limiting pools of IFT motors are necessary for agreement with data whereas tubulin may either be limited or in excess.
Although our results suggest that having a depolymerizer—one which is ballistically transported to the tip and then diffuses back—provides an appealing model for simultaneous length control, such a depolymeriser has yet to be identified in Chlamydomonas. Further experiments such as the single molecule turnaround experiments pioneered recently in C. elegans (Mijalkovic et al., 2018) are needed to establish the identity of the hypothesized depolymerizer. While recent experiments have shown kinesin-13 to be involved in length control in Giardia (McInally et al., 2019), the observation that only negligible amounts of flagellar kinesin-13 are present at steady-state (Wang et al., 2013) appears to preclude it from being the candidate depolymerizer of our model. Other candidates include aurora-like kinase CALK, which has been shown to influence disassembly through its state of phosphorylation (Luo et al., 2011; Cao et al., 2013), and CNK2, a NIMA-related protein kinase known to localize to flagella (Bradley and Quarmby, 2005) whose absence yields Chlamydomonas with abnormally long flagella and decreased disassembly rates (Hilton et al., 2013).
On the side of theory, a promising avenue to further test the active disassembly model and discriminate between different possibilities is to study length fluctuations about steady-state. This could be done using agent-based stochastic simulations, stochastic differential equations, or a combination of the two. The fluctuation spectra of each model provides a signature that can be used to assess both the general model framework and to test the autocorrelation timescales predicted by each model (Amir and Balaban, 2018).
Materials and methods
Single flagellum
We first consider a single flagellum with time-dependent length . Given the highly regular structure of the axoneme revealed by cryo-electron microscopy (Bui et al., 2008; Barber et al., 2012), we assume a constant cross-sectional area in which case the flagellar geometry is fully described by its length. In our model the flagellar assembly rate is proportional to the flux of IFT particles times the tubulin carried per particle. The tubulin carried is proportional to the amount of free tubulin , assuming mass action kinetics in the basal pool (i.e. constant probability per time of tubulin binding to an IFT particle). For now we take the disassembly speed to be equal to a constant . This yields the growth rate
As described in Results, we assume that motors are conserved having total number , where is the number of motors freely available in the basal protein pool, is the number of motors moving ballistically on the flagellum in IFT particles, and is the number of motors moving diffusively.
We assume a limiting pool of tubulin in addition to the limiting pool of motors, that is tubulin is conserved with total amount . As the flagellum grows it incorporates more tubulin and the size of the free tubulin pool decreases. (In reality , where is the amount of tubulin undergoing IFT, but this correction is negligible; the amount of tubulin moving ballistically in IFT satisfies , and consequently for the parameters contained in Table 1).
Flux balance
Request a detailed protocolIn our model the flux, or injection rate, is proportional to the number of free molecular motors so that
according to mass action kinetics with first-order rate constant . By mass action and conservation of motors, the flux of motors may be expressed as
The ballistic flux is related to concentration in a simple manner. It satisfies , where is the average concentration of motors moving in the anterograde direction and is the anterograde velocity. (As mentioned previously, it follows from the quasi-steady state assumption that the injection flux, anterograde flux, and retrograde flux are equal so that there is a single flux .) Therefore
In quasi-steady state, the diffusive flux must equal the injection rate . By Fick’s law, for constant the concentration profile of the diffusing motors is linear, that is
in which is the average concentration along the flagellum and is the diffusion constant (Figure 2b (inset)). We treat the flagellar base as a diffusive sink by fixing the boundary condition , which assumes that motors in the basal pool cannot leak diffusively into the flagellum; instead they attach to the microtubules in the axoneme and move directionally toward the tip (More general boundary conditions are discussed in Appendix 2.) This implies that
Therefore
for the diffusively-moving motors. Substituting the expressions Equation 30 and Equation 33 for and into Equation 29 results in
The denominator is a quadratic function in length, and it is interesting to note that the flux has a similar functional form to the familiar substrate production rate in Michaelis-Menten enzyme kinetics (Fall, 2002); this is because of the separation of timescales assumption invoked in both derivations.
Using the above expression for the flux in the growth rate Equation 27 together with the relation gives
Solving Equation 35 for the steady-state results in a quadratic equation for . One root is always negative, leaving the solution
which is positive provided that . (We remind the reader that terms such as are to be interpreted as .) This inequality provides a theoretical lower limit on the product of total motors and tubulin needed to obtain a positive steady-state length. When the inequality is not satisfied, that is , the disassembly term dominates and the length shrinks to zero.
We evaluate the stability of this solution by linearizing about . Expanding to first order in , we find that the steady-state is stable, that is with a positive constant given by
Based on the parameters estimated in Appendix 2, the associated timescale is approximately 15 min, which is consistent with experiment. This timescale is long compared to the few tens of seconds needed for molecular motors to traverse the flagellum in IFT, which justifies a posteriori our approximation of IFT as a quasi-steady state process.
We next consider the parameter space associated with the length dynamics. Introducing the nondimensional length and nondimensional time , we may rewrite Equation 35 in terms of the dimensionless parameters , , , and as
We interpret these parameters as follows: is the ratio of disassembly and assembly rates, is the fraction of the tubulin pool taken up by the flagellum at steady-state, is the ratio of the ballistic timescale of IFT transport to the injection timescale , and is an analogous ratio of the diffusive timescale to the injection timescale. (We could equivalently think of and as ratios of lengthscales related to the same physical processes.)
In terms of the experimentally measured parameters and those estimated in Appendix 2, we find , , , and . The relatively small values of and lead us to consider the limit (i.e. no tubulin depletion) and (i.e. instantaneous ballistic motion). In this limit, we have
nearly recovering the model of Hendel et al. (2018) with the distinction that, as mentioned above, in our model there is an additional constant term in the denominator. Note however that the essential difference between our model and Hendel et al. (2018) lies in their effective control mechanism on the number of motors loaded on the flagella, which is not captured by any differences in these formulas (see Discussion).
Two flagella
In the case of two flagella with lengths and the length dynamics are given by
where and for denote the fluxes and free amounts of tubulin for the two flagella, which may be equal when the biomolecule pools are shared. We consider various modes of coupling between the flagella giving rise to different forms of and and their consequences for length control. To assess whether a model achieves simultaneous length control we analyze the stability of solutions to the following steady-state equations:
Here, we focus on short-time behavior, that is whether a candidate model yields rapid length equalization, and do not include the protein replenishment that takes place over a longer timescale.
Tubulin separate, motors separate
Request a detailed protocolThe presence of separate tubulin pools and separate motor pools leads to two uncoupled instances of the single flagellum dynamics, that is
Setting and leads to steady state lengths given by Equation 36.
Tubulin separate, motors shared
Request a detailed protocolIn the case of separate tubulin pools, we have and . The flux may be calculated according to
where the factor of one-half comes from assuming equal injection probability into either flagellum. Further, and
so that
This yields the flagellar length dynamics
The steady-state equations are given by
When it follows from subtracting the above equations that . The steady-state equations are identical to the corresponding steady-state Equation 35 for a single flagellum with replaced by . Therefore the steady-state length satisfies Equation 36 upon rescaling :
The steady state lengths are only equal if , which is not the case after asymmetrical depletion of tubulin pools by severing (see Figure 3bii).
Tubulin shared, motors separate
Request a detailed protocolWe next consider the case in which tubulin is shared but the motor pools are separate. The separate motor pools yield decoupled fluxes identical to Equations 44 and 45:
which leads to the systems of equations
This system of equations is similar to the case of no sharing given by Equations 44-45, with the notable exception that the equations are coupled through the shared tubulin pool term . The resulting steady-state equations are identical to those of Equation 35 for a single flagellum, with replaced by and replaced by . Therefore the steady-state lengths satisfy Equation 36 upon rescaling and . This model yields simultaneous length control, and the resulting steady-state lengths satisfy only if . It follows that the steady-state lengths are unequal after severing because one of the two motor pools is depleted.
The model equations have a similar form to existing models (Marshall et al., 2005; Hendel et al., 2018), in which the assembly rates involve a factor of and either a or -dependence in the denominator, for as discussed earlier in the context of a single growing flagellum. Although the equations are similar, the absence of length equalization in our model (Figure 3biii) contrasts with the length equalization achieved in Hendel et al. (2018) by an additional control mechanism that instantaneously replenishes the number of motors on the flagellum after severing. As noted in the Discussion, the importance of protein replenishment for the model appears to be inconsistent with experimental results (Rosenbaum et al., 1969), which show that length equalization occurs even in the absence of new protein synthesis.
Tubulin shared, motors shared
Request a detailed protocolWe finally consider the case in which both tubulin and motors are shared through a common pool. By the shared tubulin pool assumption . Further, by the shared motor pool assumption the injection rates satisfy . Therefore
in which satisfies Equation 48. We are left with the steady-state equations
These equations are identical, so that there is only a single equation for the two unknowns and and the steady-state lengths are indeterminate provided that the disassembly rate is constant. Note that this conclusion holds regardless of the particular form of the flux.
We next use linear stability analysis to demonstrate this breakdown of simultaneous length in greater detail. Let and denote any one of the infinitely-many possible solutions to Equations 60 and 61. Letting and be the deviations from steady-state such that and , linearizing about any one of these solutions yields a matrix equation of the form
for . The 2 × 2 matrix above has an vanishing eigenvalue, as we now show. Diagonalizing in terms of the sum and difference gives
There is a vanishing eigenvalue associated to the difference of lengths, that is perturbations from steady-state in the difference of lengths do not decay on a finite timescale. This is consistent with previous results from stochastic simulations that the model with constant disassembly in which all biomolecules are shared does not yield simultaneous length control (Mohapatra et al., 2017). Noise must be included to observe this result; if fluctuations are not included, as in the deterministic ODE, any initial state with the correct sum in lengths appears stable. This is because the zero eigenvalue causes such states to be marginally stable.
Appendix 1
Stochastic simulation
The simulations are done using an agent-based model in which individual motors are tracked at each point of time, similar to that used in Hendel et al. (2018).
The simulations are run for a fixed amount of time discretized into small time steps of size . Parameters such as the diffusion constant , velocity , disassembly rate , and injection rate are fixed. Four different variables are tracked for each motor:
A variable indicating whether motors are located in the basal pool (value 0) or the flagellum (value 1). In the case of multiple flagella with shared motors, each flagellum is assigned a different positive integer.
A Boolean variable indicating whether motors are undergoing ballistic or diffusive motion.
The position of motor for .
The amount of tubulin carried by the motor ( at the time of injection).
The following processes take place at each time step:
The length of each flagellum for is decreased by a constant amount (in the case of constant disassembly) or by an amount (in the case of concentration-dependent disassembly), where the concentration is computed by counting the number of diffusing motors within 1 µm of the flagellar tip.
Motors are injected into the flagella with a probability , where is the number of motors in the basal pool.
Motors undergoing ballistic motion are advanced along the flagellum by a constant amount .
Motors undergoing diffusive motion are advanced by a random amount drawn from a normal distribution with mean 0 and variance .
If a motor undergoing ballistic motion reaches the tip of the flagellum, the flagellar length is increased by the amount of tubulin carried by the motor. The motion of the motor is changed to diffusive and its position is reset to , that is the location of the flagellar tip.
If the position of a motor moving diffusively exceeds the length of flagella, that is , it is reflected according to .
If a diffusive motor reaches the base of the flagellum, it is taken up into the basal pool and stops diffusing. This enforces the absorbing boundary condition .
By advancing the above processes in time, we obtain the length profiles .
To simulate the severing experiments, the length of one flagellum is shortened at a particular time and the motors and tubulin within the severed part of the flagellum are lost.
For those simulations with no protein replenishment (Figure 3—videos 1, 2, 3 and Figure 4—video 1) we fix the total amount of tubulin and the total number of motors . For those simulations in which we incorporate the replenishment of proteins in the basal pool (Figure 3—videos 4, 5, 6 and Figure 4—video 2), the target protein numbers and are fixed along with the timescale of replenishment . At each time step, there is a possibility of adding or removing a protein to or from the basal pool. The probability of addition or removal is proportional to the difference between the target number and the number of proteins currently in the basal pool.
The videos were made using the following parameters: , , , , , , , . In Figure 3—videos 4–6, the replenishment timescale is .
Appendix 2
Parameter estimation
To estimate the model parameters, we fit to experimentally-measured data. First, we extrapolate from Marshall et al. (2005) to estimate that the growth rate upon severing a flagellum is approximately 0.4 µm/min. This implies
We next estimate the product . To do so, we use the measured disassembly speed (Marshall and Rosenbaum, 2001) and an estimated initial tubulin pool of . This yields –. The estimate for the tubulin pool of a single flagellum is based on the reported amount 76–94 μm for the total tubulin shared by two flagella (Marshall et al., 2005), allowing for some tubulin loss after severing. The estimate for is not expected to be very precise, particularly given that the total tubulin pool size from Marshall et al. (2005) was itself obtained by fitting the parameters of a related model to data on mutants with extra flagella.
We subsequently estimate by fitting it to the steady-state length observed in the severing experiment. Returning to Equation 35, we may express in terms of quantities that have already been measured or estimated:
Plugging in the values for ,, and above as well as the measured values (Kozminski et al., 1993; Buisson et al., 2013), (Chien et al., 2017), and measured in Ludington et al. (2012), we obtain . See the parameter values and definitions in Table 1.
Boundary conditions
Here we discuss potential boundary conditions besides the diffusive sink described in Materials and methods. Repeating the previous calculation for a general prescribed concentration at the flagellar base yields
Therefore, to first order in the effect of nonzero is equivalent to rescaling in Equation 4.
We also consider the scenario in which the basal pool concentration depends on the number of free molecular motors through , where and are geometric parameters representing the flagellar cross sectional area and the basal pool volume , respectively. Mass action kinetics of injection may be used to rewrite the basal concentration as . Substituting this expression into Equation 66 results in
Therefore, to first order in the effect of the boundary condition in this case is equivalent to rescaling in Equation 4.
Appendix 3
Concentration-dependent disassembly model
Separate pools
We test various modes of coupling within the concentration-dependent depolymerization model and find that only the case of fully shared biomolecule pools is consistent with data (Appendix 3—figure 1). Unlike the constant disassembly models we have considered, there is some degree of length recovery after severing when only one biomolecule is shared. However, there is still no length equalization after severing because of the asymmetical depletion caused by separate pools.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-app3-fig1-v2.tif/full/617,/0/default.jpg)
Concentration-dependent disassembly model: simulations of severing experiment with different modes of coupling between basal pools (and no replenishment of protein levels).
We conclude that biomolecule pools are fully shared after ruling out all models that disagree with the rapid length equalization that occurs in severing experiments. (i) In the case of separate pools of both tubulin and depolymerizers, the unsevered flagellum does not decrease in length. (ii) When only depolymerizers are shared, the flagellar lengths do not equalize after severing. (iii) When only tubulin is shared, the flagellar lengths do not equalize after severing. (iv) When both tubulin and depolymerizers are shared, the flagella lengths rapidly equalize after severing, as observed in experiments.
Biomolecules in excess
Flagellar growth in Chlamydomonas is limited by the supply of proteins. As discussed in Results, the importance of depletion effects is evidenced by the shortening of the long flagellum in severing experiments and by the cyclohexamide experiments of Rosenbaum et al., 1969 that blocked new protein synthesis and resulted in shorter flagella. Within the proposed active disassembly model this limiting-pool mechanism may arise in three ways: (i) both tubulin and motors are limited, (ii) tubulin is limited and motors are in excess, or (iii) tubulin is in excess and motors are limited. In Figure 4 we considered case (i) given by Equations 20 and 21 in which both biomolecules are limited. Here for completeness we consider the other two cases as well. Case (ii) may immediately be ruled out: if the rate-limiting IFT protein were in excess, their injection rate would not depend on flagellar length as observed in experiments (Dentler, 2005). A limiting pool of IFT proteins is therefore necessary for agreement with the severing data. We further consider the following model in which tubulin is in excess:
where is constant since tubulin is in excess and the flux is given as before by
Simulations show that the model Equations 68 and 69 is consistent with severing data (Appendix 3—figure 2). Therefore, a limiting pool of IFT proteins is necessary for agreement with data, whereas tubulin may either be limited or in excess.
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-app3-fig2-v2.tif/full/617,/0/default.jpg)
Concentration-dependent disassembly model Equation 68 and 69 with tubulin in excess yields rapid length equalization consistent with severing experiments, here with replenishment timescale .
https://doi.org/10.7554/eLife.42599.021Linearization
Stability is established by linearizing about steady-state. We perform the linearization for the concentration-dependent disassembly model and discuss how the result generalizes to simultaneous length control with arbitrary flagellar number.
Let and be the deviations from steady-state such that and . From (13), the flux is given by
Keeping only terms up to first order in and in the denominator,
so that the first-order Taylor series expansion yields
Substituting this expression for the flux into the dynamical equations Equations 20 and 21 results in
where we have retained terms up to first order in and . Defining
we may write this system in the matrix form
The matrix above is the sum of a rank-one matrix and a diagonal perturbation. Roughly speaking, corresponds to the shared quantities whereas corresponds to the independent quantities. This system may be diagonalized in terms of the sum and difference to yield
Note that the eigenvalues
are both negative, so that the steady-state is stable, since is clearly positive and implies the positivity of . The fact that and are distinct is noteworthy as it provides a possible means to extract two independent parameters from experiment.
Arbitrary flagellar number
This rank-one plus diagonal matrix structure also applies to the case of arbitrary flagellar number discussed in Results. Let denote the deviation of the flagellum from steady-state for . The linearized equations satisfy
where and the matrix is of the form
with the rank-one matrix satisfying for all and the identity matrix, for example for N = 3
is straightforward to diagonalize. The vector corresponding to the sum of all lengths is an eigenvector of :
Further, for any vector such that , we have and
therefore it is an eigenvector with eigenvalue . Note that the pairwise differences
for form a convenient basis for the space of such vectors whose components sum to zero. In terms of the basis consisting of the sum and differences in lengths, the evolution equations diagonalize with eigenvalues
for the sum and differences, respectively, so that is an eigenvalue of multiplicity 1 and is an eigenvalue of multiplicity . Note that in addition to the explicit dependence of and on there is an implicit number-dependence through the steady-state length (and potentially the size of the pool ). Since , the dynamics involve two distinct timescales and corresponding to these two eigenvalues.
Different depolymerizer and rate-limiting IFT protein
In the description of the concentration-dependent disassembly model in Results, for convenience we made the assumption that the depolymerizer is the rate-limiting IFT protein. Here, we consider the more general case that the depolymerizer is not the rate-limiting IFT protein, but rather a different protein that is carried to the flagellar tip by IFT and present in excess in the basal pool. The formula derived in the manuscript for the flux represents the flux of IFT particles, and in this case the flux of depolymerizer is , where the capacity of depolymerizers per IFT particle is assumed constant. Assuming the depolymerizer diffuses back from the flagellar tip with diffusivity , the concentration of diffusing depolymerizer proteins satisfies
so that the concentration at position along the flagellum satisfies under the boundary condition . (Primes are used to distinguish the parameters for the depolymerizer from those of the rate-limiting IFT protein.) In this case the dynamical equations for length become
Therefore Equations 20 and 21 derived in Results remain valid in this more general case upon making the identification and . We have verified these results through agent-based simulations (Appendix 3—figure 3a).
![](https://iiif.elifesciences.org/lax/42599%2Felife-42599-app3-fig3-v2.tif/full/617,/0/default.jpg)
Simultaneous length control is achieved by different versions of the concentration-dependent disassembly model.
(a) Results from agent-based simulations with different populations of depolymerizer and rate-limiting IFT protein, using capacity , (inset) Concentrations along the flagellum of diffusing IFT proteins and depolymerizers, (b) An exponential profile is generated by motile proteins that bind and unbind from the flagellum, move ballistically in the anterograde direction when bound, and undergo diffusion when unbound, (inset, fit to Equation 89), (c) Agent-based simulations of processive depolymerization, in which depolymerizers diffusing near the flagellar tip have some probability of binding and removing tubulin at a fixed rate. We perform simulations of severing over a range of mean depolymerization times from 10 s to 600 s.
Exponential concentration gradient
To show that our model allows for non-linear concentration gradients, we set up agent based simulations which would lead to an exponential concentration distribution, similar to Naoz et al. (2008).
In our simulations, after injection motors travel ballistically to the tip of the flagellum where they begin to diffuse back toward the base as discussed in Appendix 1. In addition motors may detach with rate to undergo diffusive motion and reattach with rate to switch from diffusive to ballistic motion. This leads to an exponential distribution with the concentration profile of diffusing motors being
where as before = injection rate = and is given by
The concentration profile of diffusing motors with a fit to the above equation is shown along with the length regulation of flagella in Appendix 3—figure 3b for rates s−1 and s−1. The results are found to be qualitatively similar to the case of a linear concentration gradient, supporting our claim that concentration-dependent disassembly is able to control lengths independent of the precise form of the concentration gradient.
Processive depolymerizers
By taking the disassembly rate to depend on local depolymerizer concentration, we have implicitly assumed that the depolymerizer acts non-processively. However, this is not a fundamental restriction. We have used our agent-based model to explore a mechanism of processive depolymerization, in which depolymerizers diffusing within 1 micron of the flagellar tip bind at rate 9.4 × 10-2 s−1 and depolymerize at a rate of 9.7 × 10-4 µm/s. The depolymerization duration is drawn from an exponential distribution with prescribed mean. We have verified that this model achieves simultaneous length control over a range of mean depolymerization times. In Appendix 3—figure 3c we show the results of simulations using mean depolymerization times from 10 s (corresponding to 0.01 µm) up to 10 min (corresponding to 0.6 µm). For mean depolymerization times up to 100 s, we find the results to be qualitatively similar to the non-processive case, whereas for longer mean depolymerization times the steady-state lengths exhibit oscillations about the steady-state. The emergence of oscillations is not entirely surprising since processive depolymerization effectively introduces into the equations a time delay, which is known to give rise to oscillations in many contexts (Richard, 2003).
Data availability
All data analyzed during this study are contained in the published studies cited in the references. Source code of the simulations used in our work can be found here: https://github.com/pkar96/Agent-based-simulation (copy archived at https://github.com/elifesciences-publications/Agent-based-simulation).
References
-
Learning from noise: how observing stochasticity may aid microbiologyTrends in Microbiology 26:376–385.https://doi.org/10.1016/j.tim.2018.02.003
-
Evidence for four classes of microtubules in individual cellsJournal of Cell Science 2:169–192.
-
A NIMA-related kinase, Cnk2p, regulates both flagellar length and cell size in ChlamydomonasJournal of Cell Science 118:3317–3326.https://doi.org/10.1242/jcs.02455
-
Molecular architecture of inner dynein arms in situ in Chlamydomonas reinhardtii flagellaThe Journal of Cell Biology 183:923–932.https://doi.org/10.1083/jcb.200808050
-
Intraflagellar transport proteins cycle between the flagellum and its baseJournal of Cell Science 126:327–338.https://doi.org/10.1242/jcs.117069
-
Intraflagellar transport (IFT) during assembly and disassembly of Chlamydomonas flagellaThe Journal of Cell Biology 170:649–659.https://doi.org/10.1083/jcb.200412021
-
BookComputational Cell Biology Interdisciplinary Applied Mathematics, 20New York: Springer-Verlag.https://doi.org/10.1007/b97701
-
A study of the fine structure of ciliated epitheliaJournal of Morphology 94:221–281.https://doi.org/10.1002/jmor.1050940202
-
Length regulation of active biopolymers by molecular motorsPhysical Review Letters 108:258103.https://doi.org/10.1103/PhysRevLett.108.258103
-
Anomalies in the motion dynamics of long-flagella mutants of Chlamydomonas reinhardtiiJournal of Biological Physics 39:1–14.https://doi.org/10.1007/s10867-012-9282-8
-
Filament depolymerization by motor moleculesPhysical Review Letters 94:108102.https://doi.org/10.1103/PhysRevLett.94.108102
-
Organelle size equalization by a constitutive processCurrent Biology 22:2173–2179.https://doi.org/10.1016/j.cub.2012.09.040
-
Flagellar length control system: testing a simple model based on intraflagellar transport and turnoverMolecular Biology of the Cell 16:270–278.https://doi.org/10.1091/mbc.e04-07-0586
-
Cell geometry: how cells count and measure sizeAnnual Review of Biophysics 45:49–64.https://doi.org/10.1146/annurev-biophys-062215-010905
-
Physical model for the geometry of actin-based cellular protrusionsBiophysical Journal 107:576–587.https://doi.org/10.1016/j.bpj.2014.05.040
-
Intraflagellar transport: mechanisms of motor action, cooperation, and cargo deliveryThe FEBS Journal 284:2905–2931.https://doi.org/10.1111/febs.14068
-
Intraflagellar transport: the eyes have itThe Journal of Cell Biology 144:385–388.https://doi.org/10.1083/jcb.144.3.385
-
Intraflagellar transportNature Reviews Molecular Cell Biology 3:813–825.https://doi.org/10.1038/nrm952
-
Intraflagellar TransportAnnual Review of Cell and Developmental Biology 19:423–443.https://doi.org/10.1146/annurev.cellbio.19.111401.091318
-
Intraflagellar transport motors in cilia: moving along the cell's antennaThe Journal of Cell Biology 180:23–29.https://doi.org/10.1083/jcb.200709133
-
Flagellar protein dynamics in ChlamydomonasJournal of Biological Chemistry 276:763.https://doi.org/10.1074/jbc.M103184200
-
Defective flagellar assembly and length regulation in LF3 null mutants in ChlamydomonasThe Journal of Cell Biology 163:597–607.https://doi.org/10.1083/jcb.200307143
-
Flagellar regeneration requires cytoplasmic microtubule depolymerization and kinesin-13Journal of Cell Science 126:1531–1540.https://doi.org/10.1242/jcs.124255
-
A differential cargo-loading model of ciliary length regulation by IFTCurrent Biology 23:2463–2471.https://doi.org/10.1016/j.cub.2013.10.044
Article and author information
Author details
Funding
National Science Foundation (CAREER 1752024)
- Ariel Amir
National Science Foundation (DMS-1502851)
- Thomas G Fai
National Science Foundation (DMR-1610737)
- Jane Kondev
National Science Foundation (MRSEC-1420382)
- Jane Kondev
National Science Foundation (PHY-1806818)
- Prathitha Kar
Alfred P. Sloan Foundation
- Ariel Amir
Kavli Foundation
- Ariel Amir
Simons Foundation
- Lishibanya Mohapatra
- Jane Kondev
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We acknowledge useful discussions with Wallace Marshall, William Ludington, Shashank Shekhar, and Shane Mcinally. We thank Jie Lin and Ethan Levien for reading a draft of this manuscript and offering helpful feedback. We acknowledge support from National Science Foundation grants DMS-1502851 (TGF), grant DMR-1610737 (JK), MRSEC-1420382 (JK), the Simons Foundation (JK, LM), the A P Sloan Foundation (AA), NSF CAREER award 1752024 (AA), and the Kavli Institute (AA).
Copyright
© 2019, Fai et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,357
- views
-
- 351
- downloads
-
- 28
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution and dynamics. We present HILARy (High-precision Inference of Lineages in Antibody Repertoires), an efficient, fast and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and 𝑑𝑁∕𝑑𝑆 ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.
-
- Physics of Living Systems
During locomotion, soft-bodied terrestrial animals solve complex control problems at substrate interfaces, but our understanding of how they achieve this without rigid components remains incomplete. Here, we develop new all-optical methods based on optical interference in a deformable substrate to measure ground reaction forces (GRFs) with micrometre and nanonewton precision in behaving Drosophila larvae. Combining this with a kinematic analysis of substrate-interfacing features, we shed new light onto the biomechanical control of larval locomotion. Crawling in larvae measuring ~1 mm in length involves an intricate pattern of cuticle sequestration and planting, producing GRFs of 1–7 µN. We show that larvae insert and expand denticulated, feet-like structures into substrates as they move, a process not previously observed in soft-bodied animals. These ‘protopodia’ form dynamic anchors to compensate counteracting forces. Our work provides a framework for future biomechanics research in soft-bodied animals and promises to inspire improved soft-robot design.