Length regulation of multiple flagella that self-assemble from a shared pool of components

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Appendix 3
Data availability
References
Article and author information
Metrics

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.

Introduction

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).

Figure 1

Download asset Open asset

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.

Figure 2

Download asset Open asset

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 $J_{i}$ , ballistic flux $J_{b}$ , and diffusive flux $J_{d}$ 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 $M_{f}$ , the injection flux $J_{i}$ of IFT particles into the flagellum satisfies

J_{i} = k_{on} M_{f} .

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 $M$ satisfies $M = M_{f} + M_{b} + M_{d}$ , where $M_{b}$ is the number of motors moving ballistically on the flagellum in IFT particles and $M_{d}$ is the number of motors moving diffusively. Therefore we may rewrite Equation 1 as

J = k_{on} (M - M_{b} - M_{d}) .

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 $J_{i}$ , diffusive flux $J_{b}$ , and ballistic flux $J_{d}$ are balanced, that is

J_{i} = J_{b} = J_{d} = : J .

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

J = \frac{k_{on} M}{1 + k_{on} L / v + k_{on} L^{2} / 2 D},

where $k_{on}$ is the rate constant of motor injection, $M$ is the total number of motors, $v$ is the ballistic motor speed in IFT, and $D$ 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 $x / y z$ are shorthand for $x / (y z)$ , for example $k_{on} L^{2} / 2 D$ is to be read as $k_{on} L^{2} / (2 D)$ .)

Table 1

Parameter values and definitions.

Symbol	Definition	Value	Units	References
Parameters
$L_{s s}$	Steady-state length	10–12	µm	Marshall and Rosenbaum, 2001; Rosenbaum et al., 1969
$T / N$	Tubulin pool per flagellum	38–47	µm	Marshall et al., 2005
$d$	Disassembly speed	0.5	µm/min	Marshall and Rosenbaum, 2001; Ludington et al., 2012
$v$	IFT speed	2.5–3	µm/s	Kozminski et al., 1993; Buisson et al., 2013
$D$	Diffusion coefficient	1.7	µm²/s	Chien et al., 2017
$γ k_{on} M / N$	Assembly rate per tubulin	2.3 × 10^-2–3.6 × 10^-2	min⁻¹	Fit
$k_{on}$	Injection rate constant	0.8–4	min⁻¹	Fit
$γ$	Prefactor in Equation 5	2.5 × 10^-4		Estimate (Appendix 2)
Variables
$N$	Number of flagella
$T_{f}$	Free tubulin		µm
$M$	Total motors
$M_{f}$	Free motors
$M_{b}$	Ballistic motors
$M_{d}$	Diffusing motors
$J$	Flux		min⁻¹
$c_{d} (x)$	Motor concentration		µm⁻¹
$\bar{c_{d}}$	Average concentration		µm⁻¹

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 $D ≫ L v$ , the $1 / L$ scaling of Marshall and Rosenbaum (2001) is recovered whereas in the limit $D ≪ L v$ , we recover the $1 / L^{2}$ 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 $L = 0$ .

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 $J$ of IFT particles times the amount of free tubulin $T_{f}$ , so that the growth rate is given by the following ODE:

\frac{d L}{d t} = γ J T_{f} - d,

where $d$ is the disassembly speed (assumed constant for now) and $γ$ is a constant. The total amount of tubulin $T$ is assumed to be conserved for now so that $T = T_{f} + L$ , 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

\frac{d L}{d t} = \frac{γ k_{on} M}{1 + k_{on} L / v + k_{on} L^{2} / 2 D} (T - L) - d .

As shown in Materials and methods, this equation yields a stable steady-state length $L_{s s}$ 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 $L_{1} (t)$ and $L_{2} (t)$ is given by the following ODEs:

\frac{d L_{1}}{d t} = γ J_{1} T_{f, 1} - d,

\frac{d L_{2}}{d t} = γ J_{2} T_{f, 2} - d,

where $J_{i}$ and $T_{f, i}$ for $i = 1, 2$ 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 $J_{i}$ and $T_{f, i}$ 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 $J_{i}$ and $T_{f, i}$ within our model. We investigate their consequences for length control by focusing on the solutions to the steady-state equations

0 = γ J_{1} T_{f, 1} - d,

0 = γ J_{2} T_{f, 2} - d .

Length control implies that these steady-state equations must yield a unique steady-state solution for $L_{1}$ and $L_{2}$ . 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.

Figure 3 with 6 supplements see all

Download asset Open asset

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 $M_{f, 1}$ and $M_{f, 2}$ in the following form:

\frac{d L_{1}}{d t} = γ k_{on} M_{f, 1} T_{f, 1} - d,

\frac{d L_{2}}{d t} = γ k_{on} M_{f, 2} T_{f, 2} - d .

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 $L_{1}$ and a number of motors $M_{1, b} + M_{1, d}$ 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 $L_{1} \to 0$ , $T_{1} \to (T_{1} - L_{1})$ , and $M_{1} \to (M_{1} - M_{1, b} - M_{1, d})$ . Because severing does not change any of the protein levels in the basal pools, $M_{f, i}$ and $T_{f, i}$ for $i = 1, 2$ 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 $J_{1} = J_{2} = : J$ and it can be shown by a straightforward generalization of Equation 4 (see Materials and methods) that the flux satisfies

J = \frac{k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} .

Initially the tubulin pools are equal, as are the flagellar lengths. However, after the loss of material due to severing, $T_{1} \neq T_{2}$ 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 $T_{f, 1} = T_{f, 2} = : T_{f}$ , with

T_{f} = T - L_{1} - L_{2} .

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

0 = γ J (T - L_{1} - L_{2}) - d,

0 = γ J (T - L_{1} - L_{2}) - d .

Because the two steady-state equations are identical, that is $d L_{1} / d t = d L_{2} / d t$ , 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

\frac{d (L_{1} - L_{2})}{d t} = 0,

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 $T$ and motor pool $M$ 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 $T$ and $M$ of tubulin and motors are replenished over a timescale $τ_{r}$ as the cell synthesizes new protein to achieve target protein levels ${\bar{T}}_{f}$ and ${\bar{M}}_{f}$ in the basal pool:

τ_{r} \frac{d T}{d t} = {\bar{T}}_{f} - T_{f},

τ_{r} \frac{d M}{d t} = {\bar{M}}_{f} - M_{f} .

In steady-state, the above equations become $T_{f} = {\bar{T}}_{f}$ and $M_{f} = {\bar{M}}_{f}$ . The steady-state equation for length then becomes $0 = γ k_{on} {\bar{M}}_{f} {\bar{T}}_{f} - d$ , 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 $T_{f}$ and $M_{f}$ are driven to their target values ${\bar{T}}_{f}$ and ${\bar{M}}_{f}$ 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 $d_{0} + d_{1} c_{d} (L)$ , where $c_{d} (L)$ is the concentration of diffusing motors at the tip of the flagellum and $d_{1} > 0$ . 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 $c_{d} (L) = J L / D$ (Equation 31 in Materials and methods and Figure 2b (inset)), so that we may rewrite the governing equations as

\frac{d L_{1}}{d t} = γ J (T - L_{1} - L_{2}) - d_{0} - d_{1} \frac{J L_{1}}{D},

\frac{d L_{2}}{d t} = γ J (T - L_{1} - L_{2}) - d_{0} - d_{1} \frac{J L_{2}}{D},

where as before in the case of shared motors

J = \frac{k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} .

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 $L_{1, s s} = L_{2, s s} = : L_{s s}$ , and solving for the steady-state length results in

L_{s s} = (\frac{D}{v} + \frac{γ D M}{d_{0}} + \frac{M d_{1}}{2 d_{0}}) ⟮ - 1 + \sqrt{1 + (\frac{d_{0} T}{γ M D}) \frac{1 - 2 d_{0} / γ k_{on} M T}{{(1 + d_{0} / γ M v + d_{1} / 2 γ D)}^{2}}} ⟯ .

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.

Figure 4 with 2 supplements see all

Download asset Open asset

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 $τ_{r} = 5$ 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 $L_{s s} = 12.5 μ m$ , $T = 28 μ m$ , $d_{0} = 2 μ m / m i n$ , $d_{1} = 2.7 \times 10^{- 3} μ m^{2} / m i n$ , and $γ k_{on} M = 2.5 {min}^{- 1}$ ; all other parameters are as in Table 1. For the fit to Rosenbaum et al. (1969), we use the same fitting parameters except for $τ_{r} = 8 m i n$ , $L_{s s} = 10.2 μ m$ , and $d_{1} = 8.4 \times 10^{- 3} μ m^{2} / m i n$ .

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 $τ_{r} = 10$ mins, the recovery of the flagella back to their original lengths is in reasonable agreement with experimental data.

Generalization to $N > 2$ flagella

The concentration-dependent disassembly model may be generalized to arbitrary flagellar number $N$ , and here we demonstrate simultaneous length control in the case of $N = 8$ flagella (Figure 5).

Because motors are shared, the injection fluxes are equal and $J_{i} = J$ for all $i = 1 \dots N$ , with $J$ satisfying

J = \frac{k_{on} M / N}{1 + k_{on} (\sum_{i = 1}^{N} L_{i}) / N v + k_{on} (\sum_{i = 1}^{N} L_{i}^{2}) / 2 N D},

and length dynamics given by

\frac{d L_{i}}{d t} = γ J (T - \sum_{j = 1}^{N} L_{j}) - d_{0} - d_{1} \frac{J L_{i}}{D}, i = 1 \dots N,

where we have applied the boundary condition $c_{d, i} (0) = 0$ as before. Taking any pairwise difference between the $i$ ^th and $j$ ^th equations at steady-state yields immediately $L_{i, s s} = L_{j, s s}$ , so that the steady-state lengths are equal to

L_{i, s s} = (\frac{D}{v} + \frac{γ D M}{d_{0}} + \frac{M d_{1}}{N d_{0}}) (- 1 + \sqrt{1 + (\frac{2 d_{0} T}{N γ M D}) \frac{1 - N d_{0} / γ k_{on} M T}{{(1 + d_{0} / γ M v + d_{1} / N γ D)}^{2}}}),

for all $i = 1, \dots, N$ . Stability follows from analyzing the linearized equations, as shown in Appendix 3.

Figure 5

Download asset Open asset

The concentration-dependent disassembly model generalizes to arbitrary flagellar number $N$ .

We solve (Equations 25) with $N = 8$ flagella and the larger shared pool $T = 336 μ m$ , $d_{0} = 0.1 μ m / m i n$ , and $d_{1} = 120 μ m^{2} / m i n$ ; 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 $L (t)$ . 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 $J$ of IFT particles times the tubulin carried per particle. The tubulin carried is proportional to the amount of free tubulin $T_{f}$ , 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 $d$ . This yields the growth rate

\frac{d L}{d t} = γ J T_{f} - d .

As described in Results, we assume that motors are conserved having total number $M = M_{f} + M_{b} + M_{d}$ , where $M_{f}$ is the number of motors freely available in the basal protein pool, $M_{b}$ is the number of motors moving ballistically on the flagellum in IFT particles, and $M_{d}$ 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 $T = T_{f} + L$ . As the flagellum grows it incorporates more tubulin and the size of the free tubulin pool decreases. (In reality $T = T_{b} + T_{f} + L$ , where $T_{b}$ is the amount of tubulin undergoing IFT, but this correction is negligible; the amount of tubulin moving ballistically in IFT satisfies $T_{b} / T_{f} < 2 γ k_{on} M L_{s s} / v$ , and consequently $T_{b} / T_{f} < 2.6 \times 10^{- 3}$ for the parameters contained in Table 1).

Flux balance

Request a detailed protocol

In our model the flux, or injection rate, is proportional to the number of free molecular motors $M_{f}$ so that

J = k_{on} M_{f},

according to mass action kinetics with first-order rate constant $k_{on}$ . By mass action and conservation of motors, the flux of motors may be expressed as

J = k_{on} (M - M_{b} - M_{d}) .

The ballistic flux is related to concentration in a simple manner. It satisfies $J = \bar{c_{a}} v$ , where $\bar{c_{a}}$ is the average concentration of motors moving in the anterograde direction and $v$ 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 $J$ .) Therefore

M_{b} = L \bar{c_{a}} = L J / v .

In quasi-steady state, the diffusive flux $D \partial c_{d} / \partial x$ must equal the injection rate $J$ . By Fick’s law, for constant $J$ the concentration profile $c_{d} (x)$ of the diffusing motors is linear, that is

c_{d} (x) = \bar{c_{d}} + \frac{J}{D} (x - \frac{L}{2}),

in which $\bar{c_{d}}$ is the average concentration along the flagellum and $D$ is the diffusion constant (Figure 2b (inset)). We treat the flagellar base as a diffusive sink by fixing the boundary condition $c_{d} (0) = 0$ , 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

\bar{c_{d}} = \frac{J L}{2 D} .

Therefore

M_{d} = \bar{c_{d}} L = (\frac{J L}{2 D}) L = \frac{J L^{2}}{2 D},

for the diffusively-moving motors. Substituting the expressions Equation 30 and Equation 33 for $M_{b}$ and $M_{d}$ into Equation 29 results in

J = \frac{k_{on} M}{1 + k_{on} L / v + k_{on} L^{2} / 2 D} .

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 $T_{f} = T - L$ gives

\frac{d L}{d t} = \frac{γ k_{on} M}{1 + k_{on} L / v + k_{on} L^{2} / 2 D} (T - L) - d .

Solving Equation 35 for the steady-state results in a quadratic equation for $L_{s s}$ . One root is always negative, leaving the solution

L_{s s} = (\frac{D}{v} + \frac{γ D M}{d}) (- 1 + \sqrt{1 + (\frac{2 d T}{γ M D}) \frac{1 - d / γ k_{on} M T}{{(1 + d / γ M v)}^{2}}}),

which is positive provided that $T > d / γ k M$ . (We remind the reader that terms such as $d / γ k_{on} M$ are to be interpreted as $d / (γ k_{on} M)$ .) 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 $T \leq d / γ k_{on} M$ , the disassembly term dominates and the length shrinks to zero.

We evaluate the stability of this solution by linearizing about $L_{s s}$ . Expanding to first order in $Δ L := L - L_{s s}$ , we find that the steady-state is stable, that is $d (Δ L) / d t = - λ (Δ L)$ with $λ$ a positive constant given by

λ = (\frac{γ k_{on} M}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}) (1 + \frac{k_{on} / v + k_{on} L_{s s} / D}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (T - L_{s s})) .

Based on the parameters estimated in Appendix 2, the associated timescale $τ := 1 / λ$ 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 $\tilde{L} = L / L_{s s}$ and nondimensional time $\tilde{t} = t γ k_{on} M T / L_{s s}$ , we may rewrite Equation 35 in terms of the dimensionless parameters $π_{1} = d / γ k_{on} M T$ , $π_{2} = L_{s s} / T$ , $π_{3} = k_{on} L_{s s} / v$ , and $π_{4} = k_{on} L_{s s}^{2} / 2 D$ as

\frac{d \tilde{L}}{d \tilde{t}} = \frac{1 - π_{2} \tilde{L}}{1 + π_{3} \tilde{L} + π_{4} {\tilde{L}}^{2}} - π_{1} .

We interpret these parameters as follows: $π_{1}$ is the ratio of disassembly and assembly rates, $π_{2}$ is the fraction of the tubulin pool taken up by the flagellum at steady-state, $π_{3} = τ_{b} / τ_{i}$ is the ratio of the ballistic timescale $τ_{b} := L_{s s} / v$ of IFT transport to the injection timescale $τ_{i} := k_{on}^{- 1}$ , and $π_{4} = τ_{d} / τ_{i}$ is an analogous ratio of the diffusive timescale $τ_{d} = L_{s s}^{2} / 2 D$ to the injection timescale. (We could equivalently think of $π_{3}$ and $π_{4}$ 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 $π_{1} \approx 0.4$ , $π_{2} \approx 0.2$ , $π_{3} \approx 0.1$ , and $π_{4} \approx 0.8$ . The relatively small values of $π_{2}$ and $π_{3}$ lead us to consider the limit $π_{2} \to 0$ (i.e. no tubulin depletion) and $π_{3} \to 0$ (i.e. instantaneous ballistic motion). In this limit, we have

\frac{d L}{d t} = \frac{γ k_{on} M T}{1 + k_{on} L^{2} / 2 D} - d,

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 $L_{1} (t)$ and $L_{2} (t)$ the length dynamics are given by

\frac{d L_{1}}{d t} = γ J_{1} T_{f, 1} - d,

\frac{d L_{2}}{d t} = γ J_{2} T_{f, 2} - d,

where $J_{i}$ and $T_{f, i}$ for $i = 1, 2$ 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 $J_{i}$ and $T_{f, i}$ 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:

0 = γ J_{1} T_{f, 1} - d,

0 = γ J_{2} T_{f, 2} - d .

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 protocol

The presence of separate tubulin pools and separate motor pools leads to two uncoupled instances of the single flagellum dynamics, that is

\frac{d L_{1}}{d t} = \frac{γ k_{on} M_{1}}{1 + k_{on} L_{1} / v + k_{on} L_{1}^{2} / 2 D} (T_{1} - L_{1}) - d,

\frac{d L_{2}}{d t} = \frac{γ k_{on} M_{2}}{1 + k_{on} L_{2} / v + k_{on} L_{2}^{2} / 2 D} (T_{2} - L_{2}) - d .

Setting $M_{1} = M_{2} = M$ and $T_{1} = T_{2} = T$ leads to steady state lengths given by Equation 36.

Tubulin separate, motors shared

Request a detailed protocol

In the case of separate tubulin pools, we have $T_{f, 1} = T_{1} - L_{1}$ and $T_{f, 2} = T_{2} - L_{2}$ . The flux may be calculated according to

J = \frac{1}{2} k_{on} M_{f} = \frac{1}{2} k_{on} (M - M_{b} - M_{d}),

where the factor of one-half comes from assuming equal injection probability into either flagellum. Further, $M_{b} = J (L_{1} + L_{2}) / v$ and

M_{d} = \bar{c_{d, 1}} L_{1} + \bar{c_{d, 2}} L_{2} = J (\frac{L_{1}^{2}}{2 D} + \frac{L_{2}^{2}}{2 D}),

so that

J = \frac{k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} .

This yields the flagellar length dynamics

\frac{d L_{1}}{d t} = \frac{γ k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} (T_{1} - L_{1}) - d,

\frac{d L_{2}}{d t} = \frac{γ k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} (T_{2} - L_{2}) - d .

The steady-state equations are given by

0 = γ J (T_{1} - L_{1, s s}) - d,

0 = γ J (T_{2} - L_{2, s s}) - d .

When $T_{1} = T_{2}$ it follows from subtracting the above equations that $L_{1, s s} = L_{2, s s}$ . The steady-state equations are identical to the corresponding steady-state Equation 35 for a single flagellum with $M$ replaced by $M / 2$ . Therefore the steady-state length satisfies Equation 36 upon rescaling $M \to M / 2$ :

L_{1, s s} = L_{2, s s} = (\frac{D}{v} + \frac{γ D M}{2 d}) (- 1 + \sqrt{1 + (\frac{4 d T}{γ M D}) \frac{1 - 2 d / γ k_{on} M T}{{(1 + 2 d / γ M v)}^{2}}}) .

The steady state lengths are only equal if $T_{1} = T_{2}$ , which is not the case after asymmetrical depletion of tubulin pools by severing (see Figure 3b(ii)).

Tubulin shared, motors separate

Request a detailed protocol

We 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:

J_{1} = \frac{k_{on} M_{1}}{1 + k_{on} L_{1} / v + k_{on} L_{1}^{2} / 2 D},

J_{2} = \frac{k_{on} M_{2}}{1 + k_{on} L_{2} / v + k_{on} L_{2}^{2} / 2 D},

which leads to the systems of equations

\frac{d L_{1}}{d t} = \frac{γ k_{on} M_{1}}{1 + k_{on} L_{1} / v + k_{on} L_{1}^{2} / 2 D} (T - L_{1} - L_{2}) - d,

\frac{d L_{2}}{d t} = \frac{γ k_{on} M_{2}}{1 + k_{on} L_{2} / v + k_{on} L_{2}^{2} / 2 D} (T - L_{1} - L_{2}) - d .

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 $T - L_{1} - L_{2}$ . The resulting steady-state equations are identical to those of Equation 35 for a single flagellum, with $T$ replaced by $T / 2$ and $γ$ replaced by $2 γ$ . Therefore the steady-state lengths satisfy Equation 36 upon rescaling $T \to T / 2$ and $γ \to 2 γ$ . This model yields simultaneous length control, and the resulting steady-state lengths satisfy $L_{1, s s} = L_{2, s s}$ only if $M_{1} = M_{2} = M$ . 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 $T - L_{1} - L_{2}$ and either a $1 / L_{i}$ or $1 / L_{i}^{2}$ -dependence in the denominator, for $i = 1, 2$ 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 protocol

We finally consider the case in which both tubulin and motors are shared through a common pool. By the shared tubulin pool assumption $T_{f, 1} = T_{f, 2} = T - L_{1} - L_{2}$ . Further, by the shared motor pool assumption the injection rates satisfy $J_{1} = J_{2} \equiv J$ . Therefore

\frac{d L_{1}}{d t} = γ J (T - L_{1} - L_{2}) - d,

\frac{d L_{2}}{d t} = γ J (T - L_{1} - L_{2}) - d,

in which $J$ satisfies Equation 48. We are left with the steady-state equations

0 = γ J (T - L_{1, s s} - L_{2, s s}) - d,

0 = γ J (T - L_{1, s s} - L_{2, s s}) - d .

These equations are identical, so that there is only a single equation for the two unknowns $L_{1, s s}$ and $L_{2, s s}$ 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 $L_{1, s s}$ and $L_{2, s s}$ denote any one of the infinitely-many possible solutions to Equations 60 and 61. Letting $Δ L_{1}$ and $Δ L_{2}$ be the deviations from steady-state such that $L_{1} = L_{1, s s} + Δ L_{1}$ and $L_{2} = L_{2, s s} + Δ L_{2}$ , linearizing about any one of these solutions yields a matrix equation of the form

\frac{d}{d t} (\begin{matrix} Δ L_{1} \\ Δ L_{2} \end{matrix}) = - (\begin{matrix} a & a \\ a & a \end{matrix}) (\begin{matrix} Δ L_{1} \\ Δ L_{2} \end{matrix}),

for $a > 0$ . The 2 × 2 matrix above has an vanishing eigenvalue, as we now show. Diagonalizing in terms of the sum $Σ = Δ L_{1} + Δ L_{2}$ and difference $Γ = Δ L_{1} - Δ L_{2}$ gives

\frac{d}{d t} (\begin{matrix} Σ \\ Γ \end{matrix}) = - (\begin{matrix} 2 a & 0 \\ 0 & 0 \end{matrix}) (\begin{matrix} Σ \\ Γ \end{matrix}) .

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 $Δ t$ . Parameters such as the diffusion constant $D$ , velocity $v$ , disassembly rate $d$ , and injection rate $k_{o n}$ 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 $x_{j} (t)$ of motor $j$ for $j = 1, \dots, M$ .
The amount of tubulin $t_{j} (t)$ carried by the motor ( $t_{j} = γ T_{f}$ at the time of injection).

The following processes take place at each time step:

The length of each flagellum $L_{i} (t)$ for $i = 1, 2$ is decreased by a constant amount $d Δ t$ (in the case of constant disassembly) or by an amount $(d_{0} + d_{1} c_{d} (L)) Δ t$ (in the case of concentration-dependent disassembly), where the concentration $c_{d} (L)$ 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 $k_{on} M_{f} Δ t$ , where $M_{f}$ is the number of motors in the basal pool.
Motors undergoing ballistic motion are advanced along the flagellum by a constant amount $v Δ t$ .
Motors undergoing diffusive motion are advanced by a random amount drawn from a normal distribution with mean 0 and variance $2 D Δ t$ .
If a motor undergoing ballistic motion reaches the tip of the flagellum, the flagellar length is increased by the amount of tubulin $t_{j}$ carried by the motor. The motion of the motor is changed to diffusive and its position is reset to $L (t)$ , that is the location of the flagellar tip.
If the position of a motor moving diffusively exceeds the length of flagella, that is $x_{j} (t) > L (t)$ , it is reflected according to $x_{j} (t) \to L (t) - (x_{j} (t) - L (t))$ .
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 $c_{d} (0) = 0$ .

By advancing the above processes in time, we obtain the length profiles $L_{i} (t)$ .

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 $T$ and the total number of motors $M$ . 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 ${\bar{T}}_{f}$ and ${\bar{M}}_{f}$ are fixed along with the timescale of replenishment $τ_{r}$ . 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: $D = 1.7 μ m^{2} / s$ , $k_{on} = 0.075 s^{- 1}$ , $T = 38 μ m / f l a g e l l u m$ , $M = 200 / f l a g e l l u m$ , $γ = 2.5 \times 10^{- 4}$ , $d_{0} = 0.01 μ m / s$ , $d_{1} = 1.7 \times 10^{- 3} μ m^{2} / s$ , $v = 2.5 μ m / s$ . In Figure 3—videos 4–6, the replenishment timescale is $τ_{r} = 300 s$ .

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

γ k_{on} M T - d \approx 0.4 .

We next estimate the product $γ k_{on} M$ . To do so, we use the measured disassembly speed $d = 0.5 μ m / m i n$ (Marshall and Rosenbaum, 2001) and an estimated initial tubulin pool of $T = 25 - 40 μ m$ . This yields $γ k_{on} M = 2.3 \times 10^{- 2}$ – $3.6 \times 10^{- 2} {min}^{- 1}$ . The estimate $T$ 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 $γ k_{on} M$ 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 $k_{on}$ by fitting it to the steady-state length observed in the severing experiment. Returning to Equation 35, we may express $k_{on}$ in terms of quantities that have already been measured or estimated:

k_{on} \approx {(\frac{L_{s s}}{v} + \frac{L_{s s}^{2}}{2 D})}^{- 1} (\frac{γ k M}{d} (T - L_{s s}) - 1) .

Plugging in the values for $γ k M$ , $d$ , and $T$ above as well as the measured values $v = 150 μ m / m i n$ (Kozminski et al., 1993; Buisson et al., 2013), $D = 102 μ m^{2} / m i n$ (Chien et al., 2017), and $L_{s s} = 10 μ m$ measured in Ludington et al. (2012), we obtain $k_{on} = 0.8 - 4.5 m i n^{- 1}$ . See the parameter values and definitions in Table 1.

Boundary conditions

Here we discuss potential boundary conditions besides the diffusive sink $c_{d} (0) = 0$ described in Materials and methods. Repeating the previous calculation for a general prescribed concentration $c_{d} (0) = c_{0}$ at the flagellar base yields

J = \frac{k_{on} (M - c_{0} L)}{1 + k_{on} L / v + k_{on} L^{2} / 2 D} .

Therefore, to first order in $Δ L = L - L_{s s}$ the effect of nonzero $c_{0}$ is equivalent to rescaling $M$ in Equation 4.

We also consider the scenario in which the basal pool concentration depends on the number of free molecular motors $M_{f}$ through $c_{d} (0) = A M_{f} / V$ , where $A$ and $V$ are geometric parameters representing the flagellar cross sectional area $A$ and the basal pool volume $V$ , respectively. Mass action kinetics of injection $J = k_{on} M_{f}$ may be used to rewrite the basal concentration as $c_{d} (0) = A J / (k_{on} V)$ . Substituting this expression into Equation 66 results in

J = \frac{k_{on} M}{1 + (k_{on} / v + A / V) L + k_{on} L^{2} / 2 D} .

Therefore, to first order in $Δ L$ the effect of the boundary condition in this case is equivalent to rescaling $v \to v / (1 + (A / V) (v / k_{on}))$ 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.

Appendix 3—figure 1

Download asset Open asset

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:

\frac{d L_{1}}{d t} = γ J T - d_{0} - d_{1} \frac{J L_{1}}{D},

\frac{d L_{2}}{d t} = γ J T - d_{0} - d_{1} \frac{J L_{2}}{D},

where $T$ is constant since tubulin is in excess and the flux $J$ is given as before by

J = \frac{k_{on} M / 2}{1 + k_{on} (L_{1} + L_{2}) / 2 v + k_{on} (L_{1}^{2} + L_{2}^{2}) / 4 D} .

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.

Appendix 3—figure 2

Download asset Open asset

Concentration-dependent disassembly model Equation 68 and 69 with tubulin in excess yields rapid length equalization consistent with severing experiments, here with replenishment timescale $τ_{r} = 5 m i n$ .

Linearization

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 $Δ L_{1}$ and $Δ L_{2}$ be the deviations from steady-state such that $L_{1} = L_{s s} + Δ L_{1}$ and $L_{2} = L_{s s} + Δ L_{2}$ . From (13), the flux $J$ is given by

J = \frac{k_{on} M / 2}{\begin{array}{ll} 1 + k_{on} ((L_{s s} + Δ L_{1}) + (L_{s s} + Δ L_{2})) / 2 v + k_{on} ((L_{s s} + Δ L_{1})^{2} + (L_{s s} + Δ L_{2})^{2}) / 4 D \end{array}},

Keeping only terms up to first order in $Δ L_{1}$ and $Δ L_{2}$ in the denominator,

J = \frac{k_{on} M / 2}{\begin{array}{ll} 1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D + k_{on} (Δ L_{1} + Δ L_{2}) / 2 v + k_{on} L_{s s} (Δ L_{1} + Δ L_{2}) / 2 D \\ + 𝒪 (Δ L_{1}^{2}) + 𝒪 (Δ L_{1} \cdot Δ L_{2}) + 𝒪 (Δ L_{2}^{2}), \end{array}}

so that the first-order Taylor series expansion yields

\begin{array}{ll} J = \frac{k_{on} M / 2}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (1 - \frac{k_{on} (Δ L_{1} + Δ L_{2}) / 2 v + k_{on} L_{s s} (Δ L_{1} + Δ L_{2}) / 2 D}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}) \\ + 𝒪 (Δ L_{1} \cdot Δ L_{2}) + 𝒪 (Δ L_{2}^{2}), \end{array}

Substituting this expression for the flux into the dynamical equations Equations 20 and 21 results in

\begin{array}{ll} \frac{d (Δ L_{1})}{d t} & = \frac{d_{1} k_{on} M / 2 D}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (- Δ L_{1}) \\ + \frac{γ k_{on} M / 2}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (1 + \frac{(k_{on} / 2 v + k_{on} L_{s s} / 2 D) (T - (2 + d_{1} / γ D) L_{s s})}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}) (- Δ L_{1} - Δ L_{2}), \end{array}

\begin{array}{ll} \frac{d (Δ L_{2})}{d t} & = \frac{d_{1} k_{on} M / 2 D}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (- Δ L_{2}) \\ + \frac{γ k_{on} M / 2}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (1 + \frac{(k_{on} / 2 v + k_{on} L_{s s} / 2 D) (T - (2 + d_{1} / γ D) L_{s s})}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}) (- Δ L_{1} - Δ L_{2}), \end{array}

where we have retained terms up to first order in $Δ L_{1}$ and $Δ L_{2}$ . Defining

\begin{array}{ll} a & = \frac{γ k_{on} M / 2}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D} (1 + \frac{(k_{on} / 2 v + k_{on} L_{s s} / 2 D) (T - (2 + d_{1} / γ D) L_{s s})}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}), \\ b & = \frac{d_{1} k_{on} M / 2 D}{1 + k_{on} L_{s s} / v + k_{on} L_{s s}^{2} / 2 D}, \end{array}

we may write this system in the matrix form

\frac{d}{d t} (\begin{matrix} Δ L_{1} \\ Δ L_{2} \end{matrix}) = - (\begin{matrix} a + b & a \\ a & a + b \end{matrix}) (\begin{matrix} Δ L_{1} \\ Δ L_{2} \end{matrix}) .

The matrix above is the sum of a rank-one matrix and a diagonal perturbation. Roughly speaking, $a$ corresponds to the shared quantities whereas $b$ corresponds to the independent quantities. This system may be diagonalized in terms of the sum $Σ = Δ L_{1} + Δ L_{2}$ and difference $Γ = Δ L_{1} - Δ L_{2}$ to yield

\frac{d}{d t} (\begin{matrix} Σ \\ Γ \end{matrix}) = - (\begin{matrix} 2 a + b & 0 \\ 0 & b \end{matrix}) (\begin{matrix} Σ \\ Γ \end{matrix}) .

Note that the eigenvalues

λ_{Σ} = - (2 a + b),

λ_{Γ} = - b,

are both negative, so that the steady-state is stable, since $b$ is clearly positive and $T > 2 L_{s s}$ implies the positivity of $a$ . 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 $N$ discussed in Results. Let $Δ L_{i} := L_{i} - L_{s s}$ denote the deviation of the $i th$ flagellum from steady-state for $i = 1, \dots, N$ . The linearized equations satisfy

Δ 𝐋 = M Δ 𝐋,

where $Δ 𝐋 = (Δ L_{1}, \dots, Δ L_{N})$ and the matrix $M$ is of the form

M = - a R - b I,

with $R = {𝟏𝟏}^{T}$ the rank-one matrix satisfying $R_{i j} = 1$ for all $i, j$ and $I$ the $N \times N$ identity matrix, for example for N = 3

M = - a (\begin{matrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{matrix}) - b (\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}) = (\begin{matrix} - a - b & - a & - a \\ - a & - a - b & - a \\ - a & - a & - a - b \end{matrix}) .

$M$ is straightforward to diagonalize. The vector $𝐯_{1} := 𝟏 = {(1, \dots, 1)}^{T}$ corresponding to the sum of all lengths is an eigenvector of $M$ :

\begin{array}{ll} M 1 & = - a R 1 - b I 1 \\ = - (a N + b) 1 . \end{array}

Further, for any vector $𝐱 = (x_{1}, \dots, x_{N})$ such that $\sum_{i = 1}^{N} x_{i} = 0$ , we have $R 𝐱 = 0$ and

M 𝐱 = - b 𝐱,

therefore it is an eigenvector with eigenvalue $λ = - b$ . Note that the pairwise differences

\begin{array}{ll} v_{k} = (0, . . .0, \underset{(k - 1)^{s t} e n t r y}{\underset{⏟}{1}}, \underset{k^{t h} e n t r y}{\underset{⏟}{- 1}}, 0, . . ., 0) \end{array}

for $k = 2, \dots, N$ form a convenient basis for the space of such vectors $𝐱$ whose components sum to zero. In terms of the basis ${𝐯_{1}, 𝐯_{2}, \dots, 𝐯_{N}}$ consisting of the sum and differences in lengths, the evolution equations diagonalize with eigenvalues

λ_{Σ} (N) = - a (N + b),

λ_{Γ} (N) = - b,

for the sum and differences, respectively, so that $λ_{Σ} (N)$ is an eigenvalue of multiplicity 1 and $λ_{Γ} (N)$ is an eigenvalue of multiplicity $N - 1$ . Note that in addition to the explicit dependence of $λ_{Σ} (N)$ and $λ_{Γ} (N)$ on $N$ there is an implicit number-dependence through the steady-state length (and potentially the size of the pool $T$ ). Since $λ_{Σ} (N) \neq λ_{Γ} (N)$ , the dynamics involve two distinct timescales $τ_{Σ} (N) = λ_{Σ} {(N)}^{- 1}$ and $τ_{Γ} (N) = λ_{Γ} {(N)}^{- 1}$ 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 $J$ represents the flux of IFT particles, and in this case the flux of depolymerizer is $K J$ , where the capacity $K$ of depolymerizers per IFT particle is assumed constant. Assuming the depolymerizer diffuses back from the flagellar tip with diffusivity $D^{'}$ , the concentration $c_{d}^{'}$ of diffusing depolymerizer proteins satisfies

D^{'} \frac{\partial c_{d}^{'}}{\partial x} = K J,

so that the concentration at position $x$ along the flagellum satisfies $c_{d}^{'} (x) = c_{0}^{'} + \frac{K J}{D^{'}} x$ under the boundary condition $c_{d}^{'} (0) = c_{0}^{'}$ . (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

\frac{d L_{1}}{d t} = γ J (T - L_{1} - L_{2}) - d_{0} - d_{1} (c_{0}^{'} + \frac{K J L_{1}}{D^{'}}),

\frac{d L_{2}}{d t} = γ J (T - L_{1} - L_{2}) - d_{0} - d_{1} (c_{0}^{'} + \frac{K J L_{2}}{D^{'}}),

Therefore Equations 20 and 21 derived in Results remain valid in this more general case upon making the identification $d_{0} \to d_{0} + d_{1} c_{0}^{'}$ and $d_{1} \to \frac{K D}{D^{'}} d_{1}$ . We have verified these results through agent-based simulations (Appendix 3—figure 3a).

Appendix 3—figure 3

Download asset Open asset

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 $K = 2$ , (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 $k_{d}$ to undergo diffusive motion and reattach with rate $k_{a}$ to switch from diffusive to ballistic motion. This leads to an exponential distribution with the concentration profile of diffusing motors being

C (x) = \frac{J v}{\sqrt{d^{2} k_{d}^{2} + 4 k_{a} v^{2} d}} (\exp (λ_{+} x) - \exp (λ_{-} x)),

where as before $J$ = injection rate = $k_{o n} M_{f} / 2$ and $λ_{\pm}$ is given by

λ_{\pm} = \frac{- 1 \pm \sqrt{1 + (4 k_{a} / D) {(v / k_{d})}^{2}}}{2 v / k_{d}} .

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 $k_{a} = 6.3 \times 10^{- 2}$ s⁻¹ and $k_{d} = 3.13 \times 10^{- 4}$ s⁻¹. 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⁻¹ 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

1. Amir A
2. Balaban NQ
(2018) Learning from noise: how observing stochasticity may aid microbiology
Trends in Microbiology 26:376–385.

https://doi.org/10.1016/j.tim.2018.02.003
- PubMed
- Google Scholar
(2012) Three-dimensional structure of the radial spokes reveals heterogeneity and interactions with dyneins in Chlamydomonas flagella
Molecular Biology of the Cell 23:111–120.

https://doi.org/10.1091/mbc.e11-08-0692
- PubMed
- Google Scholar
1. Behnke O
2. Forer A
(1967)
Evidence for four classes of microtubules in individual cells

Journal of Cell Science 2:169–192.
- PubMed
- Google Scholar
1. Bradley BA
2. Quarmby LM
(2005) A NIMA-related kinase, Cnk2p, regulates both flagellar length and cell size in Chlamydomonas
Journal of Cell Science 118:3317–3326.

https://doi.org/10.1242/jcs.02455
- PubMed
- Google Scholar
1. Bui KH
2. Sakakibara H
3. Movassagh T
4. Oiwa K
5. Ishikawa T
(2008) Molecular architecture of inner dynein arms in situ in Chlamydomonas reinhardtii flagella
The Journal of Cell Biology 183:923–932.

https://doi.org/10.1083/jcb.200808050
- PubMed
- Google Scholar
(2013) Intraflagellar transport proteins cycle between the flagellum and its base
Journal of Cell Science 126:327–338.

https://doi.org/10.1242/jcs.117069
- PubMed
- Google Scholar
1. Cao M
2. Meng D
3. Wang L
4. Bei S
5. Snell WJ
6. Pan J
(2013) Activation loop phosphorylation of a protein kinase is a molecular marker of organelle size that dynamically reports flagellar length
PNAS 110:12337–12342.

https://doi.org/10.1073/pnas.1302364110
- PubMed
- Google Scholar
1. Chien A
2. Shih SM
3. Bower R
4. Tritschler D
5. Porter ME
6. Yildiz A
(2017) Dynamics of the IFT machinery at the ciliary tip
eLife 6:e28606.

https://doi.org/10.7554/eLife.28606
- PubMed
- Google Scholar
(1998) Chlamydomonas Kinesin-II–dependent Intraflagellar Transport (IFT): IFT Particles Contain Proteins Required for Ciliary Assembly in Caenorhabditis elegans Sensory Neurons
The Journal of Cell Biology 141:993–1008.

https://doi.org/10.1083/jcb.141.4.993
- Google Scholar
1. Cole DG
(2003) The intraflagellar transport machinery of Chlamydomonas reinhardtii
Traffic 4:435–442.

https://doi.org/10.1034/j.1600-0854.2003.t01-1-00103.x
- PubMed
- Google Scholar
1. Dentler W
(2005) Intraflagellar transport (IFT) during assembly and disassembly of Chlamydomonas flagella
The Journal of Cell Biology 170:649–659.

https://doi.org/10.1083/jcb.200412021
- PubMed
- Google Scholar
Book
1. Fall CP
(2002) Computational Cell Biology Interdisciplinary Applied Mathematics, 20
New York: Springer-Verlag.

https://doi.org/10.1007/b97701
- Google Scholar
1. Fawcett DW
2. Porter KR
(1954) A study of the fine structure of ciliated epithelia
Journal of Morphology 94:221–281.

https://doi.org/10.1002/jmor.1050940202
- Google Scholar
1. Hara Y
2. Merten CA
(2015) Dynein-Based accumulation of membranes regulates nuclear expansion in xenopus laevis egg extracts
Developmental Cell 33:562–575.

https://doi.org/10.1016/j.devcel.2015.04.016
- PubMed
- Google Scholar
(2018) Diffusion as a ruler: modeling kinesin diffusion as a Length Sensor for Intraflagellar Transport
Biophysical Journal 114:663–674.

https://doi.org/10.1016/j.bpj.2017.11.3784
- PubMed
- Google Scholar
(2013) The kinases LF4 and CNK2 control ciliary length by feedback regulation of assembly and disassembly rates
Current Biology 23:2208–2214.

https://doi.org/10.1016/j.cub.2013.09.038
- PubMed
- Google Scholar
(2012) Length regulation of active biopolymers by molecular motors
Physical Review Letters 108:258103.

https://doi.org/10.1103/PhysRevLett.108.258103
- PubMed
- Google Scholar
1. Khona DK
2. Rao VG
3. Motiwalla MJ
4. Varma PC
5. Kashyap AR
6. Das K
7. Shirolikar SM
8. Borde L
9. Dharmadhikari JA
10. Dharmadhikari AK
11. Mukhopadhyay S
12. Mathur D
13. D'Souza JS
(2013) Anomalies in the motion dynamics of long-flagella mutants of Chlamydomonas reinhardtii
Journal of Biological Physics 39:1–14.

https://doi.org/10.1007/s10867-012-9282-8
- PubMed
- Google Scholar
(2005) Filament depolymerization by motor molecules
Physical Review Letters 94:108102.

https://doi.org/10.1103/PhysRevLett.94.108102
- PubMed
- Google Scholar
(1993) A motility in the eukaryotic flagellum unrelated to flagellar beating
PNAS 90:5519–5523.

https://doi.org/10.1073/pnas.90.12.5519
- PubMed
- Google Scholar
1. Ludington WB
2. Shi LZ
3. Zhu Q
4. Berns MW
5. Marshall WF
(2012) Organelle size equalization by a constitutive process
Current Biology 22:2173–2179.

https://doi.org/10.1016/j.cub.2012.09.040
- PubMed
- Google Scholar
1. Luo M
2. Cao M
3. Kan Y
4. Li G
5. Snell W
6. Pan J
(2011) The phosphorylation state of an aurora-like kinase marks the length of growing flagella in Chlamydomonas
Current Biology 21:586–591.

https://doi.org/10.1016/j.cub.2011.02.046
- PubMed
- Google Scholar
(2005) Flagellar length control system: testing a simple model based on intraflagellar transport and turnover
Molecular Biology of the Cell 16:270–278.

https://doi.org/10.1091/mbc.e04-07-0586
- PubMed
- Google Scholar
1. Marshall WF
(2016) Cell geometry: how cells count and measure size
Annual Review of Biophysics 45:49–64.

https://doi.org/10.1146/annurev-biophys-062215-010905
- PubMed
- Google Scholar
1. Marshall WF
2. Rosenbaum JL
(2001) Intraflagellar transport balances continuous turnover of outer doublet microtubules: implications for flagellar length control
The Journal of Cell Biology 155:405–414.

https://doi.org/10.1083/jcb.200106141
- PubMed
- Google Scholar
Preprint
(2019) Length-dependent disassembly maintains four different flagellar lengths in Giardia
bioRxiv.

https://doi.org/10.1101/647115
- Google Scholar
(2018) Single-Molecule turnarounds of Intraflagellar transport at the C. elegans Ciliary Tip
Cell Reports 25:1701–1707.

https://doi.org/10.1016/j.celrep.2018.10.050
- PubMed
- Google Scholar
Book
1. Milo R
2. Phillips R
(2015)
Cell Biology by the Numbers

Garland Science.
- Google Scholar
(2017) The Limiting-Pool mechanism fails to control the size of multiple organelles
Cell Systems 4:559–567.

https://doi.org/10.1016/j.cels.2017.04.011
- PubMed
- Google Scholar
1. Morga B
2. Bastin P
(2013) Getting to the heart of intraflagellar transport using Trypanosoma and Chlamydomonas models: the strength is in their differences
Cilia 2:16.

https://doi.org/10.1186/2046-2530-2-16
- PubMed
- Google Scholar
1. Naoz M
2. Manor U
3. Sakaguchi H
4. Kachar B
5. Gov NS
(2008) Protein localization by actin treadmilling and molecular motors regulates stereocilia shape and treadmilling rate
Biophysical Journal 95:5706–5718.

https://doi.org/10.1529/biophysj.108.143453
- PubMed
- Google Scholar
1. Orbach R
2. Howard J
(2019) The dynamic and structural properties of axonemal tubulins support the high length stability of cilia
Nature Communications 10:1838.

https://doi.org/10.1038/s41467-019-09779-6
- PubMed
- Google Scholar
1. Orly G
2. Naoz M
3. Gov NS
(2014) Physical model for the geometry of actin-based cellular protrusions
Biophysical Journal 107:576–587.

https://doi.org/10.1016/j.bpj.2014.05.040
- PubMed
- Google Scholar
1. Piao T
2. Luo M
3. Wang L
4. Guo Y
5. Li D
6. Li P
7. Snell WJ
8. Pan J
(2009) A microtubule depolymerizing kinesin functions during both flagellar disassembly and flagellar assembly in Chlamydomonas
PNAS 106:4713–4718.

https://doi.org/10.1073/pnas.0808671106
- PubMed
- Google Scholar
1. Pinkoviezky I
2. Gov NS
(2014) Traffic jams and shocks of molecular motors inside cellular protrusions
Physical Review E 89:052703.

https://doi.org/10.1103/PhysRevE.89.052703
- Google Scholar
1. Pinkoviezky I
2. Gov NS
(2017) Exclusion and hierarchy of time scales lead to spatial segregation of molecular motors in cellular protrusions
Physical Review Letters 118:018102.

https://doi.org/10.1103/PhysRevLett.118.018102
- PubMed
- Google Scholar
(2017) Intraflagellar transport: mechanisms of motor action, cooperation, and cargo delivery
The FEBS Journal 284:2905–2931.

https://doi.org/10.1111/febs.14068
- PubMed
- Google Scholar
(2014) Molecular mechanisms for microtubule length regulation by kinesin-8 and XMAP215 proteins
Interface Focus 4:20140031.

https://doi.org/10.1098/rsfs.2014.0031
- PubMed
- Google Scholar
1. Richard J-P
(2003) Time-delay systems: an overview of some recent advances and open problems
Automatica 39:1667–1694.

https://doi.org/10.1016/S0005-1098(03)00167-5
- Google Scholar
(1969) Flagellar elongation and shortening in Chlamydomonas. The use of cycloheximide and colchicine to study the synthesis and assembly of flagellar proteins
The Journal of Cell Biology 41:600–619.

https://doi.org/10.1083/jcb.41.2.600
- PubMed
- Google Scholar
(1999) Intraflagellar transport: the eyes have it
The Journal of Cell Biology 144:385–388.

https://doi.org/10.1083/jcb.144.3.385
- PubMed
- Google Scholar
1. Rosenbaum JL
2. Witman GB
(2002) Intraflagellar transport
Nature Reviews Molecular Cell Biology 3:813–825.

https://doi.org/10.1038/nrm952
- PubMed
- Google Scholar
1. Scholey JM
(2003) Intraflagellar Transport
Annual Review of Cell and Developmental Biology 19:423–443.

https://doi.org/10.1146/annurev.cellbio.19.111401.091318
- Google Scholar
1. Scholey JM
(2008) Intraflagellar transport motors in cilia: moving along the cell's antenna
The Journal of Cell Biology 180:23–29.

https://doi.org/10.1083/jcb.200709133
- PubMed
- Google Scholar
1. Song L
2. Dentler WL
(2001) Flagellar protein dynamics in Chlamydomonas
Journal of Biological Chemistry 276:763.

https://doi.org/10.1074/jbc.M103184200
- Google Scholar
(2003) Defective flagellar assembly and length regulation in LF3 null mutants in Chlamydomonas
The Journal of Cell Biology 163:597–607.

https://doi.org/10.1083/jcb.200307143
- PubMed
- Google Scholar
1. Varga V
2. Leduc C
3. Bormuth V
4. Diez S
5. Howard J
(2009) Kinesin-8 motors act cooperatively to mediate length-dependent microtubule depolymerization
Cell 138:1174–1183.

https://doi.org/10.1016/j.cell.2009.07.032
- PubMed
- Google Scholar
1. Wang L
2. Piao T
3. Cao M
4. Qin T
5. Huang L
6. Deng H
7. Mao T
8. Pan J
(2013) Flagellar regeneration requires cytoplasmic microtubule depolymerization and kinesin-13
Journal of Cell Science 126:1531–1540.

https://doi.org/10.1242/jcs.124255
- PubMed
- Google Scholar
(1972) Chlamydomonas flagella. I. isolation and electrophoretic analysis of microtubules, matrix, membranes, and mastigonemes
The Journal of Cell Biology 54:507–539.

https://doi.org/10.1083/jcb.54.3.507
- PubMed
- Google Scholar
1. Wren KN
2. Craft JM
3. Tritschler D
4. Schauer A
5. Patel DK
6. Smith EF
7. Porter ME
8. Kner P
9. Lechtreck KF
(2013) A differential cargo-loading model of ciliary length regulation by IFT
Current Biology 23:2463–2471.

https://doi.org/10.1016/j.cub.2013.10.044
- PubMed
- Google Scholar
1. Yochelis A
2. Ebrahim S
3. Millis B
4. Cui R
5. Kachar B
6. Naoz M
7. Gov NS
(2015) Self-organization of waves and pulse trains by molecular motors in cellular protrusions
Scientific Reports 5:13521.

https://doi.org/10.1038/srep13521
- PubMed
- Google Scholar
(2004) Nuclear structure in Cancer cells
Nature Reviews Cancer 4:677–687.

https://doi.org/10.1038/nrc1430
- PubMed
- Google Scholar

Article and author information

Author details

Thomas G Fai

Department of Mathematics, Brandeis University, Waltham, United States

Contribution
Conceptualization, Formal analysis, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0383-5217
Lishibanya Mohapatra

Department of Physics, Brandeis University, Waltham, United States

Contribution
Conceptualization, Formal analysis, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared
Prathitha Kar

Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, United States

Contribution
Formal analysis, Validation, Investigation, Visualization

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4091-6860
Jane Kondev

Department of Physics, Brandeis University, Waltham, United States

Contribution
Conceptualization, Formal analysis, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7522-7144
Ariel Amir

Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, United States

Contribution
Conceptualization, Formal analysis, Writing—original draft, Writing—review and editing

For correspondence
arielamir@seas.harvard.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-2611-0139

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 PHY-1806818 (PK), 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

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,386

views
352

downloads
28

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Thomas G Fai
Lishibanya Mohapatra
Prathitha Kar
Jane Kondev
Ariel Amir

(2019)

Length regulation of multiple flagella that self-assemble from a shared pool of components

eLife 8:e42599.

https://doi.org/10.7554/eLife.42599

Categories and tags

Research organism

Chlamydomonas reinhardtii

Share this article

Cite this article

Experimental background.

IFT in Chlamydomonas with diffusive return of kinesin-2 to the base.

Parameter values and definitions.

Length dynamics of two flagella assembling from shared pools of building blocks.

Concentration-dependent disassembly model: simultaneous length control is achieved using shared tubulin and shared depolymerizers.

The concentration-dependent disassembly model generalizes to arbitrary flagellar number N.

Concentration-dependent disassembly model: simulations of severing experiment with different modes of coupling between basal pools (and no replenishment of protein levels).

Concentration-dependent disassembly model Equation 68 and 69 with tubulin in excess yields rapid length equalization consistent with severing experiments, here with replenishment timescale τr=5min.

Simultaneous length control is achieved by different versions of the concentration-dependent disassembly model.

Author details

Thomas G Fai

Contribution

Competing interests

Lishibanya Mohapatra

Contribution

Competing interests

Prathitha Kar

Contribution

Competing interests

Jane Kondev

Contribution

Competing interests

Ariel Amir

Contribution

For correspondence

Competing interests

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)

Categories and tags

Research organism

Further reading

The concentration-dependent disassembly model generalizes to arbitrary flagellar number $N$ .

Concentration-dependent disassembly model Equation 68 and 69 with tubulin in excess yields rapid length equalization consistent with severing experiments, here with replenishment timescale $τ_{r} = 5 m i n$ .