Email updates

Keep up to date with the latest news and content from BMC Biophysics and BioMed Central.

Open Access Research article

Using default constraints of the spindle assembly checkpoint to estimate the associated chemical rates

Khanh Dao Duc and David Holcman*

Author Affiliations

Institute for Biology (IBENS), Group of Computational Biology and Applied Mathematics, Ecole Normale Supérieure, 46 rue d'Ulm 75005 Paris, France

For all author emails, please log on.

BMC Biophysics 2012, 5:1  doi:10.1186/2046-1682-5-1

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/2046-1682/5/1


Received:11 November 2011
Accepted:19 January 2012
Published:19 January 2012

© 2012 Duc and Holcman; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background

Default activation of the spindle assembly checkpoint provides severe constraints on the underlying biochemical activation rates: on one hand, the cell cannot divide before all chromosomes are aligned, but on the other hand, when they are ready, the separation is quite fast, lasting a few minutes. Our purpose is to use these opposed constraints to estimate the associated chemical rates.

Results

To analyze the above constraints, we develop a markovian model to describe the dynamics of Cdc20 molecules. We compute the probability for no APC/C activation before time t, the distribution of Cdc20 at equilibrium and the mean time to complete APC/C activation after all chromosomes are attached.

Conclusions

By studying Cdc20 inhibition and the activation time, we obtain a range for the main chemical reaction rates regulating the spindle assembly checkpoint and transition to anaphase.

I. Background

A fundamental step in cell division consists in the alignment of each pair of chromosomes. This process occurs during metaphase, where centrosome nucleated microtubules interact with the chromosomes kinetochores to build the mitotic spindle. Only after all chromosomes have become aligned at the metaphase plate and when every kinetochore is properly attached to a bundle of microtubules, the cell enters anaphase. To prevent premature progression to anaphase, even if all-but-one of the kinetochores have been attached and the chromosomes are aligned, unattached or improperly attached kinetochores generate a signal inhibiting the anaphase activators. This process is called the spindle assembly checkpoint (SAC).

Although the exact mechanisms of the SAC and anaphase processes are still unclear, several key steps have been identified. Sister chromatids are initially bound by proteins such as cohesin. During anaphase onset, separase protein cleaves cohesin, thus allowing the sister chromatids to separate [1]. Usually, separase is prevented from cleaving cohesin through its association with another protein called securin. Securin can be ubiqitylated by the activated Cdc20-anaphase promoting complex/cyclosome (Cdc20-APC/C) at the kinetochores [2]. However, when the kinetochores are not all properly attached, the SAC enables the ubiquitylation and inhibition of Cdc20 binding with APC/C [2,3]. The mechanisms leading to Cdc20 ubiquitylation involve several proteins such as Mad2, BubR1, Bub3 [4]. Current models [5-10] of the checkpoint propose that Mad2 protein has a crucial role either to sequester Cdc20, or acts in conjunction with the BubR1 and Bub3 proteins to form an inhibitor called the Mitotic Checkpoint Complex (MCC). Mad2 changes conformation to bind Cdc20 tightly via a safety belt mechanism [5,7,8], while another inhibitory complex comprised of BubR1 and Bub3 has also been identified [9,10].

The SAC has been modeled at a molecular level, however the parameters used [11-14] may not necessarily reflect in vivo dynamics [15]. For example, these modelings do not take into account the finite number of binding sites for Cdc20. In addition, the constant flux assumption [11] made for molecules reaching a kinetochore impacts the APC/C activation, leading to an overestima-tion for the catalytic activity. We shall revisit here some of these major assumptions in the construction of our model. An improvement of these models was recently achieved [16,17] by taking into account the finite number of binding sites at kinetochores [16], leading also to an estimate for the MCC chemical rates associated with SAC.

In the present article, our purpose is to study the inhibition followed by its fast activation of Cdc20, which is the key activator of the anaphase promoting complex. As the number of kinetochores implied in the SAC is small, the forward binding rate of a chemical reaction as it is classically computed in the continuously concentrated limit cannot be applied. To adequately describe chemical reactions in microdomains [18,19], where targets such as kinetochores have to be reached by the anaphase activators, we use a stochastic approach. Using Markovian equations [18,20] to account for the binding dynamic associated with a finite number of molecules, we compute the time dependent probability that the spindle is not initiated before time t (formula 23) and then the mean time to induce anaphase (implicit formula 29). We apply our analysis to PTK2 cells and thus, we obtain some quantitative constraints on the Cdc20 production rate and the MCC concentration to guarantee strong inhibition of Cdc20 by the SAC. Using different parameter values (cell size, number of chromosomes...), our method can be extended to other cell types and organisms, providing a general framework to study the dynamics of activators during the spindle checkpoint and the anaphase transition.

II. Methods

Markovian modeling of APC/C activation and Cdc20 inhibition

We describe here the time evolution of the joint probability distribution of Cdc20 molecules and of APC/C complex activation, the later being responsible for the chromosome separation. In this model, APC/C is located on the chromosomes (figure 1), and is a target of the Cdc20 molecule, although there are some conflicting evidences that APC/C is located on the kinetochores [21]. This assumption can affect the binding rate, but does not impact the construction of our model. For a cell containing N chromosomes, the targets of the Cdc20 molecules are the N associated kinetochores, containing the APC/C complexes. When a Cdc20 molecule reaches a kinetochore, it activates the APC/C complex and this can trigger a cascade of reactions (detailed in the background section) leading to the separation of the sister chromatids. The goal of the spindle assembly checkpoint signal is to prevent this activation of APC/C by Cdc20, when at least one of the chromosomes is not properly attached to the microtubules responsible for the chromatids migration. The spindle assembly checkpoint signal consists in the production of proteins such as Mad2, BubR1 and Mad3 [4], generated by unattached kinetochores. These proteins diffuse in large quantity in the cell to inhibit the APC/C binding by Cdc20 molecules. Indeed, these proteins form with Cdc20 a complex called mitotic checkpoint complex (MCC). Similar complexes can be found in yeast, in which a BubR1-related Mad3 protein might inhibit Cdc20 as a pseudo-substrate [22-24]. The formation of this complex results in Cdc20 ubiquitylation, which prevents APC/C activation. In our model, MCC will represent the complex of inhibitory proteins before it binds Cdc20, in contrast with the usual terminology where MCC includes Cdc20. Cdc20 molecules are produced from the dissociation of a complex [25], which could be a subcomplex of Mad2 and Cdc20 resulting from MCC:Cdc20 disassembly, promoted by p31 [26]. In summary, the above chemical reactions can be summarized as

thumbnailFigure 1. A schematic view of the spindle assembly checkpoint and anaphase. A: Before all chromosomes are all attached, the mitotic checkpoint complex inhibits the Cdc20 molecules binding with APC/C to prevent premature separation of sister chromatids. This signal ubiquitylates Cdc20. B: When all chromosomes are properly attached, the inhibiting signal is shut down. C and D: Activation of APC/C triggers the separation of the chromatids and ultimately the anaphase.

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M1">View MathML</a>

(1)

where the rate λ measures the production of Cdc20 and k-1 the degradation, while μ is the arrival rate for a CdC20 to an APC/C site. We shall compute in the next paragraphs the joint probability

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M2">View MathML</a>

(2)

that the APC/C is not activated by any free Cdc20 molecule at time t. To compute this probability, we first derive a Markov equation. The difficulty is that this joint probability contains a discrete variable counting the number of Cdc20 molecules and a binary one, which monitors whether or not an activation of APC/C by Cdc20 has occurred before time t. The state space of this Markov process is completed by adding a state describing that activation occurred before time t. It is modeled as an absorbing state of probability

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M3">View MathML</a>

(3)

which accounts for all the activations which have happened before time t from all the states k. Starting with k active Cdc20 molecules and no activation, there are three possible transitions (figure 2): 1) one Cdc20 molecule is inhibited, so that k-1 active molecules are left 2) one Cdc20 molecule activates the APC/C 3) one Cdc20 molecule is generated, leading to the transition from k to k+1 active molecules. Thus, the probabilities pk satisfy the chemical master equations [18,20]

thumbnailFigure 2. Markov diagram for the probability of number of Cdc20 molecules. Cdc20 at state k is generated and destroyed at rate λ(S -- k) and k-1k respectively. Diffusing Cdc20 molecules bind to the APC/C complex to trigger the separation of sister chromatids.

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M4">View MathML</a>

(4)

The production rate of Cdc20 molecules is proportional to the number of remaining available complex molecules given by λ(S - k). Indeed the Cdc20 molecules are produced during metaphase by dissociation from a pool of complexes, which limits the level of Cdc20 to a maximum of S molecules [25]. In addition, we shall emphasize that there is another interpretation of equation (4): indeed, dissociation of the MCC:Cdc20 complex, producing Cdc20 molecule, leads also to equation (4). When none of the N target kinetochores have been activated, the arrival rate for a Cdc20 molecule to an APC/C is μN, where <a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M5">View MathML</a> and τ is the mean time for a Cdc20 molecule to reach the APC/C site. This mean time can be approximated by [27-29]

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M6">View MathML</a>

(5)

where D is the diffusion coefficient of a Cdc20 molecule, a is the radius of the APC/C complex, R the radius of the cell. As the mechanisms underlying the production and the regulation of MCC are still unclear, we consider that the MCC concentration is homogeneous and remains constant over time to guarantee a robust inhibition of Cdc20. Thus, k-1 is given by the Smoluchovski formula for the binding rate of a Brownian particle

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M7">View MathML</a>

(6)

where [MCC] is the concentration of MCC, uniform over the cell and b is the radius of the Cdc20 binding site. When the SAC starts, no free Cdc20 molecules are present in the cell, thus we choose for the initial conditions pk(0) = δk,0. Because there can only be S Cdc20 molecules, we have for all time t, pk(t) = 0 and k > S.

The probability for no activation

To quantify the inhibition capacity of the SAC, we estimate the probability P(t) that at time t, no APC/C has been activated, so that no chromosomal migration could have been initiated. This probability is given by

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M8">View MathML</a>

(7)

We shall compute P(t) using the generating function

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M9">View MathML</a>

(8)

Using equation (4), f satisfies a first order PDE

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M10">View MathML</a>

(9)

Using the characteristics method, we look for a solution of

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M11">View MathML</a>

(10)

which is of Riccati type. From the classical substitution <a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M12">View MathML</a> we obtain the linear second order differential equation

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M13">View MathML</a>

(11)

Thus, the solution for characteristics is

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M14">View MathML</a>

(12)

where C is a constant and r1 and r2 are the two roots of the quadratic polynomial associated with (11)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M15">View MathML</a>

(13)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M16">View MathML</a>

(14)

Along one of these characteristics, f satisfies the linear first order ODE

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M17">View MathML</a>

(15)

The general solution for equation (15) is

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M18">View MathML</a>

(16)

At time t = 0, we have pk(0) = δk0, initial conditions set K = 1. To find p(t), we shall select the characteristic for which at time t, xC(t) = 1. Solving this, yields to

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M19">View MathML</a>

(17)

and we obtain the characteristic

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M20">View MathML</a>

(18)

Finally, the probability P(t) of no activation is

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M21">View MathML</a>

(19)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M22">View MathML</a>

(20)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M23">View MathML</a>

(21)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M24">View MathML</a>

(22)

Finally,

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M25">View MathML</a>

(23)

P is a decreasing function of time, and remains constant for λ = 0 and μ = 0. In figure 3, we plot P as a function of time for different values of λ and k-1, and as a function of λ and k-1 at a given time. It is a decreasing function of λ (increasing the Cdc20 production rate decreases the probability of activation) and a decreasing function of k-1 (increasing the inhibition of Cdc20 increases the probability for no activation).

thumbnailFigure 3. The probability P for no activation during the SAC is represented as a function of the time and the rates λ and k-1. A :We plot P(t) as a function of time for different values of λ and k-1 = 1. B: as a function of time for different values of k-1 and λ = 0.1. C: at time t = 1200s as a function of λ and k-1. The parameters are given in table 1.

The distribution of Cdc20 at equilibrium

After the last chromosome attached and thus all kinetochores are properly positioned, the inhibition of APC/C:Cdc20 binding is suppressed and anaphase can start. The initial condition for the number of Cdc20 molecules for this new phase is the one obtained at equilibrium from the previous phase, in which Cdc20 is produced and destroyed by the SAC. When there are k Cdc20 molecules, the production rate is given by λ(S-k) and the destruction rate k-1(k+1). Thus, the probability pk(t|NA) that k Cdc20 molecules are inside the cell, conditioned that no activation has occurred, satisfies the Master equations

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M26">View MathML</a>

(24)

The equilibrium probabilities pk(∞) and the mean number <a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M27">View MathML</a> for such a system is [30]

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M28">View MathML</a>

(25)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M29">View MathML</a>

(26)

When the SAC is suppressed, Cdc20 is no longer inhibited and can activate APC/C to trigger anaphase. Using the distribution computed here we compute in the next section the mean time for complete separation of sister chromatids during anaphase.

Activation of APC/C

When all kinetochores are properly attached, the SAC is shut down and the activation of APC/C:Cdc20 complex triggers a cascade of reactions leading to cohesin ubiquitylation at the chromosome sites [1]. To study the time for such activation, we consider that production and degradation of MCC are fast enough so that the MCC concentration decreases rapidly once all the kinetochores are attached. In that case, we can neglect the transient time for the rate k-1 to decay to 0 and thus we take the equilibrium Cdc20 concentrations as the initial conditions for the activation of APC/C. For our analysis, we further consider that the time for all kinetochores to be attached is not too short compared to the degradation and production time scale, so that the CdC20 concentration is close to equilibrium. When there are k Cdc20 present in the cell and m of them are bound to APC/C, it results that the association rate is μ(N - m)(k - m) [18]. We shall now estimate the joint probabilities that there are k Cdc20 molecules, and m activated APC/C by Cdc20

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M43">View MathML</a>

From the state (k,m), the transition rate to activation of APC/C located on another kinetochore is then μ(k - m)(N - m). Thus, we get the following Markov chain (represented in figure 4)

thumbnailFigure 4. Schematic representation of the Markov Chain associated with the joint probability pk,m(t) to have k bounds APC/C and m free Cdc20 molecules.

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M30">View MathML</a>

(27)

with the initial condition

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M31">View MathML</a>

(28)

computed in equation (25). In that case, the mean time τ that all APC/C are activated is obtained by analyzing a continuous markov process that reaches a given threshold [20]. Using formula 12 of [20], the mean time to threshold is expressed as a sum

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M32">View MathML</a>

(29)

where <a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M33">View MathML</a>. Integrating the system of equation (28) from 0 to + ∞ with the initial conditions

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M34">View MathML</a>

(30)

leads to

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M35">View MathML</a>

(31)

In practice, we solve this linear system of equations numerically and in figure 5, we plot the mean τ as a function of the parameters k-1 and λ. We find that τ is a decreasing function of λ (the faster Cdc20 is produced, the faster the threshold of bindings is reached) and an increasing function of the rate k-1 (inhibition decreases the number of Cdc20 at equilibrium and thus the time to reach the threshold, after the source of inhibition is terminated). These variations go in the opposite direction compared to the probability of no activation during the SAC. Thus we expect that using the probability P and this mean time τ will lead to limit the range of the parameters λ and k-1 as we will describe now.

thumbnailFigure 5. The time τ is plotted as a function of the parameters λ and k-1. The parameter valuers are given in table 1.

III. Results

Quantitative constraints on the rates λ and k-1

We now apply our previous modeling to determine the rates of production λ and the backward binding rate k-1. Indeed, during SAC, a strong inhibition signal imposes that the probability for no activation remains very high and thus, the degradation rate k-1 has to be high enough compared to the production rate λ. In contrast, a fast activation during anaphase forces the mean time to activate all the kinetochores to be short, thus the production rate λ has to be high. These opposite constraints allow us to determine a range for the parameters λ and k-1. We use the following quantitative constraints.

1. First, the probability P of no activation remains high enough during the time τ1 where all chro mosomes get properly attached in the metaphase plate. It has been estimated that τ1 ≈ 20 min [31]. Thus by fixing a threshold of 0.95 for the probabil ity P(τ1) that no activation occurred before time τ1

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M36">View MathML</a>

(32)

2. Second, during the anaphase onset, the time 〈τs〉 for all chromosomes to get separated is short. Since APC/C activation triggers the chromosome sep aration, we can consider that τS is the time for all APC/C to get activated. Indeed, biophysical data [32] suggest that τs should be limited in time τ' ≈ 10 min. Thus,

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M37">View MathML</a>

(33)

Using formula (23) for the probability P(τ1) and integrating numerically the time 〈τS〉 from the matrix equation (31), we determine a range of validity for these parameters by a geometrical domain Ω represented in figure 6, as the intersection Ω = Ω1 ∩ Ω2, where

thumbnailFigure 6. Representation of the domain Ω (red). A: The probability for no activation at time τ = 10min as a function of parameters λ and k-1. Other parameters are those of table I B : The mean time to threshold as a as a function of parameters λ and k-1. Other parameters are those of table I. C : The curve 1 is given as the level line associated τ = 10 min in figure B. The curve 2 is the level line associated with the probability P given by P(τ) = 0.95 in figure A. These two curves determine the domain Ω satisfying conditions C1 and C2.

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M38">View MathML</a>

(34)

and

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M39">View MathML</a>

(35)

We tested the prediction of our model on PTK2 cells, originating from kangaroo rat kidney, used in studies on mitosis because there are only a few large chromosomes and the cells remain flattened during mitosis. For these cells, the concentration of bound complex from which Cdc20 is produced is approximately 50 nM [25] during interphase. This concentration is of the same order as the one reported in [13], and is equivalent to 3000 molecules restricted in a volume of 100 μm3 (for a flat cell of size 10 μm and of height 1 μm, leading to a volume of 10 × 10 × 1 μm3). For larger cells, the number of molecules can be multiplied by 10 or 100. Thus, during the SAC, it is tempting to think that the system escapes the stochastic limit. However, because the number of Cdc20 is small at early metaphase and limited by the inhibition of MCC, the stochastic regime is still controlling the behavior of the system and in addition, the inhibition is strong enough to maintain a low level of Cdc20. In figure 6, we represented the two domains Ω1 and Ω2, the first is on the left of curve 1, while the second is on the right side of curve 2. Our analysis can be generalized by changing the two conditions C1 and C2 for specific cell types. The other parameters are summarized in table 1. Surprisingly, Ω is not bounded, but it provides an interesting and new range for the parameters.

IV. Discussion and Conclusion

Based on the two main constraints C1 and C2, we presented here a Markovian analysis to estimate two fundamental rates regulating the spindle assembly checkpoint. This idea of using physical and timing/inhibition constraints was used before [11] to infer SAC characteristics and compare different models. In a different mathematical framework, our study is based on describing precisely the role of the SAC, which is to prevent a premature separation of chromosomes. Such an event is stochastic, determined by an accidental binding which leads to activate APC/C. Estimating parameters such as chemical rates in the context of stochastic systems is challenging, as new methods and tools have to be used and developed, notably in statistical inference of Markov processes [19,33]. In our case, we use a Markovian approach to relate the chemical rates to the characteristics of the SAC observed at the cellular level. To guarantee an inhibition strong enough to prevent accidental binding, the probability that no activation occurred has to be high enough and plays a crucial role in determining the validity of the production rate parameters. In contrast to other quantitative studies of the SAC, we also study the premises of separation of sister chromatids during anaphase. As the SAC determines the amount of Cdc20 when anaphase starts, the time for activation of all APC/C located on kinetochores should not be neglected. Finally, in Figure 6, we obtain a range for the rates λS (which is the production rate at the beginning of SAC) and k-1. Actually, this range approximatively depends on the ratio <a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M40">View MathML</a>, which satisfies

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M41">View MathML</a>

(36)

The constraint ρ ≥ 50 gives the minimum value required to produce enough Cdc20 molecules to activate the APC/C before 10 minutes. However, this ratio should not be too high, because an overproduction of Cdc20 could trigger a premature anaphase and thus ρ cannot also be too large, limited to 90. To close the domain Ω, a third constraint can be added by providing an upper bound for k-1. Because k-1 is given by the Smolu-chovski formula 2πbD[MCC], limiting the concentration [MCC] would precisely limit the chemical rate k-1. For example, when the number of MCC is in the range of 10000 (which corresponds to the Mad2 concentration of 200nM found in [34], and used in [16,17]), we obtain that k-1 ≈ 24 and in that case, we approximatively get for the production rate

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/1/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/1/mathml/M42">View MathML</a>

(37)

For example, fixing the value λS = 0.3, we find that anaphase is triggered after a mean time of 239 s (4 minutes), while the probability for no activation at time t = 20min is P = 0.96, which satisfies the biophysical constraints C1 and C2 described above. To close the domain Ω more tightly, several considerations can be suggested. The fact that once free/active, Cdc20 needs to find a securin-separase complex and bring it to the APC/C (alternatively find the APC/C first and not get ubiquitylated) might constrain the time somewhat further. We also did not take into account the effective time to get to equilibrium for Cdc20 or the time to clear inhibition of MCC. If the number of MCC is too large, this time cannot be neglected and would provide an explanation for overexpressing inhibitors that prevent anaphase [35]. It would be interesting to account for the dynamics of MCC [26,36], and how it can influence the transition phase between SAC and anaphase onset [2,35]. Finally, the present study can be extended to various cell geometry with different size and with different number of chromosomes.

Authors' contributions

KDD and DH designed research and wrote the paper. KDD and DH performed research. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank J. Pines, Z. Xu and D. Peric for fruitful discussions. D. H research is supported by an ERC-Starting Grant.

References

  1. Uhlmann F, Lottspeich F, Nasmyth K: Nature.

    Sister-chromatid separation at anaphase onset is promoted by cleavage of the cohesin subunit Scc1 1999, 400(6739):19-21. OpenURL

  2. Nilsson J, Yekezare M, Minshull J, Pines J: The APC/C maintains the spindle assembly checkpoint by targeting Cdc20 for destruction.

    Nature Cell Biology 2008, 10:1411-1420. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Passmore L, Booth R, Vnien-Bryan C, Ludtke S, Fioretto C, Johnson L, Chiu W, Barford D: Structural Analysis of the Anaphase-Promoting Complex Reveals Multiple Active Sites and Insights into Polyubiquitylation.

    Molecular Cell 2005, 20:6. OpenURL

  4. Sear RP, Howard M: Modeling dual pathways for the metazoan spindle assembly checkpoint.

    Proc Natl Acad Sci USA 2006, 103(45):16758-63. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Sironi L, et al.: Crystal structure of the tetrameric Mad1-Mad2 core complex: implications of a safety belt binding mechanism for the spindle checkpoint.

    The EMBO journal 2002, 21:24962506. OpenURL

  6. Vigneron S, Prieto S, Bernis C, Labbe JC, Castro A, Lorca T: Kinetochore Localization of Spindle Checkpoint Proteins: Who Controls Whom?

    Molecular Biology of the Cell 2004, 15:45844596. OpenURL

  7. Mapelli M, Massimiliano L, Santaguida S, Musacchio A: The mad2 conformational dimer: structure and implications for the spindle assembly checkpoint.

    Cell 2007, 131:730743. OpenURL

  8. Luo X, Tang Z, Rizo J, Yu H: The Mad2 spindle checkpoint protein undergoes similar major conformational changes upon binding to either Mad1 or Cdc20.

    Molecular cell 2002, 9:5971. OpenURL

  9. Tang Z, Bharadwaj R, Li B, Yu H: Mad2-Independent inhibition of APCCdc20 by the mitotic checkpoint protein BubR1.

    Developmental cell 2001, 1:227237. OpenURL

  10. Fang G: Checkpoint protein BubR1 acts synergisti-cally with Mad2 to inhibit anaphase-promoting complex.

    Molecular biology of the cell 2002, 13:755766. OpenURL

  11. Doncic A, Ben-Jacob E, Barkai : Evaluating putative mechanisms of the mitotic spindle checkpoint.

    PNAS 2005, 102(18):6332-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Doncic A, Ben-Jacob E, Einav S, Barkai1N : Reverse Engineering of the Spindle Assembly Checkpoint.

    PLoS One 2009., 4(8) OpenURL

  13. Ibrahim B, Diekmann S, Schmitt E, Dittrich P: In-Silico Modeling of the Mitotic Spindle Assembly Checkpoint.

    PLoS ONE 2008., 3(2) OpenURL

  14. Mistry H, MacCallum D, Jackson R, Chaplain M, Davidson F: Modeling the temporal evolution of the spindle assembly checkpoint and role of Aurora B kinase.

    PNAS 2008., 105(51) OpenURL

  15. Ciliberto A, Shah JV: A quantitative systems view of the spindle assembly checkpoint.

    The EMBO Journal 2009, 28:2162-2173. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Lohel M, Ibrahim B, Diekmann S, Dittrich P: The role of localization in the operation of the mitotic spindle assembly checkpoint.

    Cell Cycle 2009, 8(16):2650-60.

    Epub 2009 Aug 29

    PubMed Abstract | Publisher Full Text OpenURL

  17. Simonetta M, Manzoni R, Mosca R, Mapelli M, Mas-similiano L, Vink M, Novak B, Musacchio A, Ciliberto A: The influence of catalysis on mad2 activation dynamics.

    PLoS Biol 2009, 7(1):e10. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Holcman D, Schuss Z: Stochastic chemical reactions in microdomains.

    J Chem Phys 2005, 122(11):114710. PubMed Abstract | Publisher Full Text OpenURL

  19. Wilkinson DJ: Stochastic modelling for quantitative description of heterogeneous biological systems.

    Nature Reviews Genetics 2009, 10(2):122-133. PubMed Abstract | Publisher Full Text OpenURL

  20. Dao Duc K, Holcman D: Threshold activation for stochastic chemical reactions.

    Phys Rev E 2010, 81:041107. OpenURL

  21. Howell BJ, Moree B, Farrar EM, Stewart S, Fang G, Salmon ED: Spindle checkpoint protein dynamics at kinetochores in living cells.

    Curr Biol 2004, 14(11):953-64.

    2004

    PubMed Abstract | Publisher Full Text OpenURL

  22. Pan J, Chen R-H: Spindle checkpoint regulates Cdc20p stability in Saccharomyces cerevisiae.

    Genes & Dev 2004, 18:1439-1451. PubMed Abstract | Publisher Full Text OpenURL

  23. Sczaniecka M, Feoktistova A, May KM, Chen JS, Blyth J, Gould KL, Hardwick KG: The Spindle Checkpoint Functions of Mad3 and Mad2 depend on a Mad3 KEN Box-mediated Interaction with Cdc20-Anaphase-promoting Complex (APC/C).

    J Biol Chem 2008, 283(34):23039-23047. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Burton JL, Solomon MJ: Mad3p, a pseudosubstrate inhibitor of APCCdc20 in the spindle assembly checkpoint.

    Genes Dev 2007, 21(6):655-67. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Wang Z, Shah JV, Berns MW, Cleveland DW: In vivo quantitative studies of dynamic intracellu-lar processes using fluorescence correlation spectroscopy.

    Biophysical Journal 2006, 91. OpenURL

  26. Teichner A, Eytan E, Sitry-Shevah D, Miniowitz-Shemtov S, Dumin E, Gromis J, Hershko A: p31comet Promotes disassembly of the mitotic checkpoint complex in an ATP-dependent process.

    Proc Natl Acad Sci USA 2011, 108(8):3187-92.

    Epub 2011 Feb 7

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Holcman D, Schuss Z: "Escape through a small opening: receptor trafficking in a synaptic membrane".

    J of Statistical Physics 2004, 117(5/6):191-230. OpenURL

  28. Singer A, Schuss Z, Holcman D, Eisenberg B: Narrow Escape I.

    J Stat Phys 2006, 122(3):437-463. Publisher Full Text OpenURL

  29. Schuss Z, Singer A, Holcman D: The narrow escape problem for diffusion in cellular microdomains.

    Proc Natl Acad Sci USA 2007, 104:1609816103. OpenURL

  30. Kleinrock L: Queueing systems. Wiley; 1976.

  31. Rieder CL, Schultz A, Cole R, Sluder G: Anaphase onset in vertebrate somatic cells is controlled by a checkpoint that monitors sister kinetochore attachment to the spindle.

    J Cell Biol 1994, 127(5):13011310. OpenURL

  32. Meraldi P, Draviam VM, Sorger PK: Timing and Checkpoints in the Regulation of Mitotic Progression.

    Developmental Cell 2004, 7:4560. OpenURL

  33. Opper M, Sanguinetti G: Variational inference for Markov jump processes. In Advances in Neural Information Processing Systems. MIT Press; 2008:1105-1112. OpenURL

  34. Shah JV, Botvinick E, Bonday Z, Furnari F, Berns M, Cleveland DW: Dynamics of centromere and kinetochore proteins; implications for checkpoint signaling and silencing.

    Curr Biol 2004, 14:942952. OpenURL

  35. Wassmann K, Niault T, Maro B: Metaphase I arrest upon activation of the Mad2-dependent spindle checkpoint in mouse ovocytes.

    Current Biology 2003, 13(18):1596-1608. PubMed Abstract | Publisher Full Text OpenURL

  36. Ibrahim B, Schmidt E, Dittrich P, Diekmann S: MCC Assembly is not Combined With Full Cdc20 Sequestering.

    Biosystems 2009, 95(1):35-50. PubMed Abstract | Publisher Full Text OpenURL

  37. Herzog F, Primorac I, Dube P, Lenart P, Sander B, Mechtler K, Stark H, Peters JM: Structure of the anaphase-promoting complex/cyclosome interacting with a mitotic checkpoint complex.

    Science 2009, 323(5920):1477-81. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL