 Research article
 Open Access
Using default constraints of the spindle assembly checkpoint to estimate the associated chemical rates
 Khanh Dao Duc^{1} and
 David Holcman^{1}Email author
https://doi.org/10.1186/2046168251
© Duc and Holcman; licensee BioMed Central Ltd. 2012
 Received: 11 November 2011
 Accepted: 19 January 2012
 Published: 19 January 2012

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.
Keywords
 Sister Chromatid
 Spindle Assembly Checkpoint
 Binding Rate
 Chemical Master Equation
 PTK2 Cell
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 allbutone 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 Cdc20anaphase promoting complex/cyclosome (Cdc20APC/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 overestimation 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
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
The distribution of Cdc20 at equilibrium
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
III. Results
Quantitative constraints on the rates λ and k_{1}
 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}
 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,
.
Parameter  Description  Value 

V  cell volume  ≈ 100 μm ^{3} [25] 
D  Diffusion coefficient of Cdc20  D = 20 μm^{2}/sec [25] 
a  radius of APC/C complex  ≈ 10^{2} μm [3] 
b  radius of MCC binding site  ≈ 2nm [37] 
N  Number of chromosomes  13 [25] 
S  Initial number of complexes  3000 [25] 
μ  binding rate $\left(=\frac{4aD}{V}\right)$  ≈ 2.10^{4}s^{1} 
IV. Discussion and Conclusion
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 securinseparase 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.
Declarations
Acknowledgements
We would like to thank J. Pines, Z. Xu and D. Peric for fruitful discussions. D. H research is supported by an ERCStarting Grant.
Authors’ Affiliations
References
 Uhlmann F, Lottspeich F, Nasmyth K: Nature. Sisterchromatid separation at anaphase onset is promoted by cleavage of the cohesin subunit Scc1. 1999, 400 (6739): 1921.Google Scholar
 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: 14111420. 10.1038/ncb1799.View ArticleGoogle Scholar
 Passmore L, Booth R, VnienBryan C, Ludtke S, Fioretto C, Johnson L, Chiu W, Barford D: Structural Analysis of the AnaphasePromoting Complex Reveals Multiple Active Sites and Insights into Polyubiquitylation. Molecular Cell. 2005, 20: 6View ArticleGoogle Scholar
 Sear RP, Howard M: Modeling dual pathways for the metazoan spindle assembly checkpoint. Proc Natl Acad Sci USA. 2006, 103 (45): 1675863. 10.1073/pnas.0603174103.ADSView ArticleGoogle Scholar
 Sironi L, et al: Crystal structure of the tetrameric Mad1Mad2 core complex: implications of a safety belt binding mechanism for the spindle checkpoint. The EMBO journal. 2002, 21: 24962506View ArticleGoogle Scholar
 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: 45844596View ArticleGoogle Scholar
 Mapelli M, Massimiliano L, Santaguida S, Musacchio A: The mad2 conformational dimer: structure and implications for the spindle assembly checkpoint. Cell. 2007, 131: 730743View ArticleGoogle Scholar
 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: 5971View ArticleGoogle Scholar
 Tang Z, Bharadwaj R, Li B, Yu H: Mad2Independent inhibition of APCCdc20 by the mitotic checkpoint protein BubR1. Developmental cell. 2001, 1: 227237View ArticleGoogle Scholar
 Fang G: Checkpoint protein BubR1 acts synergistically with Mad2 to inhibit anaphasepromoting complex. Molecular biology of the cell. 2002, 13: 755766View ArticleGoogle Scholar
 Doncic A, BenJacob E, Barkai : Evaluating putative mechanisms of the mitotic spindle checkpoint. PNAS. 2005, 102 (18): 63327. 10.1073/pnas.0409142102.ADSView ArticleGoogle Scholar
 Doncic A, BenJacob E, Einav S, Barkai1N : Reverse Engineering of the Spindle Assembly Checkpoint. PLoS One. 2009, 4 (8):Google Scholar
 Ibrahim B, Diekmann S, Schmitt E, Dittrich P: InSilico Modeling of the Mitotic Spindle Assembly Checkpoint. PLoS ONE. 2008, 3 (2):Google Scholar
 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):Google Scholar
 Ciliberto A, Shah JV: A quantitative systems view of the spindle assembly checkpoint. The EMBO Journal. 2009, 28: 21622173. 10.1038/emboj.2009.186.View ArticleGoogle Scholar
 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): 265060. 10.4161/cc.8.16.9383. Epub 2009 Aug 29View ArticleGoogle Scholar
 Simonetta M, Manzoni R, Mosca R, Mapelli M, Massimiliano L, Vink M, Novak B, Musacchio A, Ciliberto A: The influence of catalysis on mad2 activation dynamics. PLoS Biol. 2009, 7 (1): e1010.1371/journal.pbio.1000010.View ArticleGoogle Scholar
 Holcman D, Schuss Z: Stochastic chemical reactions in microdomains. J Chem Phys. 2005, 122 (11): 11471010.1063/1.1849155.ADSView ArticleGoogle Scholar
 Wilkinson DJ: Stochastic modelling for quantitative description of heterogeneous biological systems. Nature Reviews Genetics. 2009, 10 (2): 122133. 10.1038/nrg2509.MathSciNetView ArticleGoogle Scholar
 Dao Duc K, Holcman D: Threshold activation for stochastic chemical reactions. Phys Rev E. 2010, 81: 041107ADSView ArticleGoogle Scholar
 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): 95364. 10.1016/j.cub.2004.05.053. 2004View ArticleGoogle Scholar
 Pan J, Chen RH: Spindle checkpoint regulates Cdc20p stability in Saccharomyces cerevisiae. Genes & Dev. 2004, 18: 14391451. 10.1101/gad.1184204.View ArticleGoogle Scholar
 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 Boxmediated Interaction with Cdc20Anaphasepromoting Complex (APC/C). J Biol Chem. 2008, 283 (34): 2303923047. 10.1074/jbc.M803594200.View ArticleGoogle Scholar
 Burton JL, Solomon MJ: Mad3p, a pseudosubstrate inhibitor of APCCdc20 in the spindle assembly checkpoint. Genes Dev. 2007, 21 (6): 65567. 10.1101/gad.1511107.View ArticleGoogle Scholar
 Wang Z, Shah JV, Berns MW, Cleveland DW: In vivo quantitative studies of dynamic intracellular processes using fluorescence correlation spectroscopy. Biophysical Journal. 2006, 91Google Scholar
 Teichner A, Eytan E, SitryShevah D, MiniowitzShemtov S, Dumin E, Gromis J, Hershko A: p31comet Promotes disassembly of the mitotic checkpoint complex in an ATPdependent process. Proc Natl Acad Sci USA. 2011, 108 (8): 318792. 10.1073/pnas.1100023108. Epub 2011 Feb 7ADSView ArticleGoogle Scholar
 Holcman D, Schuss Z: "Escape through a small opening: receptor trafficking in a synaptic membrane". J of Statistical Physics. 2004, 117 (5/6): 191230.MathSciNetMATHGoogle Scholar
 Singer A, Schuss Z, Holcman D, Eisenberg B: Narrow Escape I. J Stat Phys. 2006, 122 (3): 437463. 10.1007/s1095500580266.ADSMathSciNetView ArticleMATHGoogle Scholar
 Schuss Z, Singer A, Holcman D: The narrow escape problem for diffusion in cellular microdomains. Proc Natl Acad Sci USA. 2007, 104: 1609816103ADSView ArticleGoogle Scholar
 Kleinrock L: Queueing systems. 1976, WileyGoogle Scholar
 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): 13011310View ArticleGoogle Scholar
 Meraldi P, Draviam VM, Sorger PK: Timing and Checkpoints in the Regulation of Mitotic Progression. Developmental Cell. 2004, 7: 4560View ArticleGoogle Scholar
 Opper M, Sanguinetti G: Variational inference for Markov jump processes. Advances in Neural Information Processing Systems. 2008, MIT Press, 11051112.Google Scholar
 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: 942952Google Scholar
 Wassmann K, Niault T, Maro B: Metaphase I arrest upon activation of the Mad2dependent spindle checkpoint in mouse ovocytes. Current Biology. 2003, 13 (18): 15961608. 10.1016/j.cub.2003.08.052.View ArticleGoogle Scholar
 Ibrahim B, Schmidt E, Dittrich P, Diekmann S: MCC Assembly is not Combined With Full Cdc20 Sequestering. Biosystems. 2009, 95 (1): 3550. 10.1016/j.biosystems.2008.06.007.View ArticleGoogle Scholar
 Herzog F, Primorac I, Dube P, Lenart P, Sander B, Mechtler K, Stark H, Peters JM: Structure of the anaphasepromoting complex/cyclosome interacting with a mitotic checkpoint complex. Science. 2009, 323 (5920): 147781. 10.1126/science.1163300.ADSView ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.