# American Institute of Mathematical Sciences

April  2018, 15(2): 429-440. doi: 10.3934/mbe.2018019

## Mathematical analysis and modeling of DNA segregation mechanisms

 Department of Mathematics and Computer Science, University of Jena, Ernst-Abbe-Platz 2,07743 Jena, Germany

* Corresponding author: Bashar Ibrahim.

Received  January 02, 2017 Accepted  May 04, 2017 Published  June 2017

The precise regulation of cell life division is indispensable to the reliable inheritance of genetic material, i.e. DNA, in successive generations of cells. This is governed by dedicated biochemical networks which ensure that all requirements are met before transition from one phase to the next. The Spindle Assembly Checkpoint (SAC) is an evolutionarily mechanism that delays mitotic progression until all chromosomes are properly linked to the mitotic spindle. During some asymmetric cell divisions, such as those observed in budding yeast, an additional mechanism, the Spindle Position Checkpoint (SPOC), is required to delay exit from mitosis until the mitotic spindle is correctly aligned. These checkpoints are complex and their elaborate spatiotemporal dynamics are challenging to understand intuitively. In this study, bistable mathematical models for both activation and silencing of mitotic checkpoints were constructed and analyzed. A one-parameter bifurcation was computed to show the realistic biochemical switches considering all signals. Numerical simulations involving systems of ODEs and PDEs were performed over various parameters, to investigate the effect of the diffusion coefficient. The results provide systems-level insights into mitotic transition and demonstrate that mathematical analysis constitutes a powerful tool for investigation of the dynamic properties of complex biomedical systems.

Citation: Bashar Ibrahim. Mathematical analysis and modeling of DNA segregation mechanisms. Mathematical Biosciences & Engineering, 2018, 15 (2) : 429-440. doi: 10.3934/mbe.2018019
##### References:

show all references

##### References:
Schematics illustrating the intracellular signaling of mitotic transition control mechanisms. (A) Dependency of spindle assembly checkpoint signaling on microtubule attachment. Kinetochores which are not attached to the spindle apparatus generate a 'wait'-signal. Chemically SAC precludes the mitotic progression by inhibiting the activation of APC/C, presumably through sequestration of the APC/C-activator Cdc20. Even a single unattached or misattached kinetochore can maintain the spindle checkpoint. If all kinetochores are correctly attached from opposite poles to the mitotic spindle, SAC is turned off and APC/C:Cdc20 formation is turned on. Paired chromosomes are held together by protein complexes called cohesin, which is depredated by active APC/C:Cdc20 complexes. $n$ refers to the number of chromosomes, as SAC is conserved from yeast to human ($n=46$ in human cells). (B) Dependency of spindle position checkpoint signaling on the correct alignment of spindle pole bodies. If the cell progresses into anaphase with a misaligned spindle, SPOC delays mitotic exit to provide the cell with extra time to make the correction. The roles of Mad2 and Cdc20 in the SAC are similar to the respective roles of Bfa1 and Tem1 in the SPOC. SPOC prevents exit from mitosis through inhibition of the MEN-activator Tem1 by Bfa1 until the spindle is properly aligned. Tem1 'inactive' indicates the GDP-bound form; 'active' refers to the GTP-bound form.
Spindle Assembly Checkpoint (SAC) model set-up and analysis of all 92 kinetochores. (A) Schematic representation of the biochemical reaction network of the SAC mechanism. SAC proteins/complexes are shown as nodes, and the interactions between them as edges. Abbreviations are APC/C, anaphase promoting complex/cyclosome; Cdc20, cell division cycle 20 homolog and MCC, Mitotic Checkpoint Complex. Unattached kinetochores enhance the production of the mitotic checkpoint complex (MCC) consisting of BubR1, Bub3, Mad2, and Cdc20. Eventually, MCC binds tightly to and inhibits the APC/C in a manner preventing Cdc20 from interacting with mitotic APC/C. Immediately after the last kinetochore attachment to microtubules, the inhibitors dissolve, eventually resulting in active APC/C. This reactivation process is known as SAC silencing, where APC/C plays a role in its feedback loop. (B) Single parameter bifurcation curve. Shown is the number of attached kinetochores versus the total concentration of the MCC inhibitor. Stable node points and stable steady states are indicated by solid lines, while unstable saddle points are shown by dashed lines. Stable and unstable steady states meet at saddle-node bifurcation points, indicated by solid circles. When the number of attached kinetochore raises above approximately 91.98, the SAC checkpoint disengages and APC/C is activated. As the cell enters anaphase, MCC lowers back to zero. The black dashed line explains how the switch flips from the SAC-active state to the SAC-inactive state as the number of attached kinetochores increase. The number of attached kinetochores in the cell determines whether MCC activity is high (SAC active state) or low (SAC silence state). The bifurcation curve can be shifted to left or right depending on the values of $k_3$ and $k_2$. The earliest shift (to the left) can take place at high values of both $k_3$ as well as $k_1$ (1.1 and 0.5 respectively). (C) Numerical simulations using ODEs, SAC model. Dynamical behavior of core SAC component concentration is plotted versus time. All parameter settings are according to Table 1 (see text for more details).
Spindle Position Checkpoint model with both SPBs signals. (A) The biochemical reactions governing the SPOC mechanism. Shown are SPOC proteins/complexes and the interactions between them. Abbreviations are (see for details [2,3]): BUB2, budding uninhibited by benomyl; BFA1, byrfour-alike-1; SPB, spindle pole body; and Tem1, Ras-like GTPase. (B) SPOC bifurcation curve, parameterized as the number of aligned SPBs versus the total concentration of the Bfa1:Bub2:Tem1. Unstable saddle points are shown by dashed lines, and stable node point steady states are indicated by solid lines. As the number of aligned SPBs increases above approximately 1.99 (about to be both correctly aligned), the SPOC checkpoint switches off. As the cell finalizes anaphase with both SPBs aligned, Tem1 activity becomes high and Bfa1:Bub2 becomes inactive. The black dashed line represents how the switch flips from the SPOC-active state to the SPOC-inactive state as the number of aligned SPBs increase. The bifurcation curve is sensitive to some parameters, and it can be very slightly shifted to left or right. (C) Numerical simulations of SPOC model. Dynamical behavior is very similar qualitatively but not quantitatively to that of the SAC model. The concentrations are drawn versus time. Tem1 concentration is kept low, as long as the SPBs are not aligned. After about 7 minutes, SPOC switches off and Tem1 is rapidly activated. Meanwhile, Bfa1 curve displays a steep downwards trend. All parameter settings are according to Table 1.
Numerical simulation of SAC Reaction-Diffusion system with spherical symmetry. (A) Simulation of two PDEs, one per components (i.e. MCC and APC/C). (B) Simulation of three PDEs, one per component (MCC, APC/C, and MCC:APC/C). Both panels were generated for various diffusion constants and showed no qualitative changes; therefore, shown are typical curves for MCC:APC/C (to the left) and APC/C (to the right).
Kinetic Parameters of the SAC and the SPOC Models
 $\mathbf{Symbol}$ $\mathbf{SAC value}$ $\mathbf{SPOC value}$ $\mathbf{Remark}$ Initial amount APC/C 0.09 $\mu M$ [37] MCC 0.15 $\mu M$ [15] Tem1 0.06 $\mu M$ [4] Bfa1 0.04 $\mu M$ [4] Bub2 0.04 $\mu M$ [4] Diffusion constants MCC 1-20 ${\mu m^2}{s^{-1}}$ This study APC/C 1.8 ${\mu m^2}{s^{-1}}$ [40] Cdc20 19.5 ${\mu m^2}{s^{-1}}$ [40] Mad2 5 ${\mu m^2}{s^{-1}}$ [13] Environment Radius of the kinetochore 0.1$\mu m$ 0.01$\mu m$ [5] Radius of the cell 10 $\mu m$ 4$\mu m$ [21] Rate constants kinetochores or SPBs 0-92 0-2 [16] $k_1$ $1-100 s^{-1}$ $1-50 s^{-1}$ This study $k_2$ $50-100 \mu\text{M}^{-1}s^{-1}$ $10-50 s^{-1}$ This study $k_{-2}$ $0.008-0.08 s^{-1}$ $0.001-0.08 s^{-1}$ This study $k_{3}$ $0.005-0.5 s^{-1}$ $0.001-0.5 s^{-1}$ This study
 $\mathbf{Symbol}$ $\mathbf{SAC value}$ $\mathbf{SPOC value}$ $\mathbf{Remark}$ Initial amount APC/C 0.09 $\mu M$ [37] MCC 0.15 $\mu M$ [15] Tem1 0.06 $\mu M$ [4] Bfa1 0.04 $\mu M$ [4] Bub2 0.04 $\mu M$ [4] Diffusion constants MCC 1-20 ${\mu m^2}{s^{-1}}$ This study APC/C 1.8 ${\mu m^2}{s^{-1}}$ [40] Cdc20 19.5 ${\mu m^2}{s^{-1}}$ [40] Mad2 5 ${\mu m^2}{s^{-1}}$ [13] Environment Radius of the kinetochore 0.1$\mu m$ 0.01$\mu m$ [5] Radius of the cell 10 $\mu m$ 4$\mu m$ [21] Rate constants kinetochores or SPBs 0-92 0-2 [16] $k_1$ $1-100 s^{-1}$ $1-50 s^{-1}$ This study $k_2$ $50-100 \mu\text{M}^{-1}s^{-1}$ $10-50 s^{-1}$ This study $k_{-2}$ $0.008-0.08 s^{-1}$ $0.001-0.08 s^{-1}$ This study $k_{3}$ $0.005-0.5 s^{-1}$ $0.001-0.5 s^{-1}$ This study
 [1] James W. Webber, Sean Holman. Microlocal analysis of a spindle transform. Inverse Problems & Imaging, 2019, 13 (2) : 231-261. doi: 10.3934/ipi.2019013 [2] Qingwen Hu, Huan Zhang. Stabilization of turning processes using spindle feedback with state-dependent delay. Discrete & Continuous Dynamical Systems - B, 2018, 23 (10) : 4329-4360. doi: 10.3934/dcdsb.2018167 [3] Sung-Hwa Lim, Se Won Lee, Byoung-Hoon Lee, Seongil Lee, Ho Woo Lee. Stochastic method for power-aware checkpoint intervals in wireless environments: Theory and application. Journal of Industrial & Management Optimization, 2012, 8 (4) : 969-986. doi: 10.3934/jimo.2012.8.969 [4] Eugene Kashdan, Dominique Duncan, Andrew Parnell, Heinz Schättler. Mathematical methods in systems biology. Mathematical Biosciences & Engineering, 2016, 13 (6) : i-ii. doi: 10.3934/mbe.201606i [5] Avner Friedman. Conservation laws in mathematical biology. Discrete & Continuous Dynamical Systems - A, 2012, 32 (9) : 3081-3097. doi: 10.3934/dcds.2012.32.3081 [6] Avner Friedman. PDE problems arising in mathematical biology. Networks & Heterogeneous Media, 2012, 7 (4) : 691-703. doi: 10.3934/nhm.2012.7.691 [7] Cristina Anton, Jian Deng, Yau Shu Wong, Yile Zhang, Weiping Zhang, Stephan Gabos, Dorothy Yu Huang, Can Jin. Modeling and simulation for toxicity assessment. Mathematical Biosciences & Engineering, 2017, 14 (3) : 581-606. doi: 10.3934/mbe.2017034 [8] Loïc Louison, Abdennebi Omrane, Harry Ozier-Lafontaine, Delphine Picart. Modeling plant nutrient uptake: Mathematical analysis and optimal control. Evolution Equations & Control Theory, 2015, 4 (2) : 193-203. doi: 10.3934/eect.2015.4.193 [9] Urszula Ledzewicz, Shuo Wang, Heinz Schättler, Nicolas André, Marie Amélie Heng, Eddy Pasquier. On drug resistance and metronomic chemotherapy: A mathematical modeling and optimal control approach. Mathematical Biosciences & Engineering, 2017, 14 (1) : 217-235. doi: 10.3934/mbe.2017014 [10] Monique Chyba, Benedetto Piccoli. Special issue on mathematical methods in systems biology. Networks & Heterogeneous Media, 2019, 14 (1) : ⅰ-ⅱ. doi: 10.3934/nhm.20191i [11] Erika T. Camacho, Christopher M. Kribs-Zaleta, Stephen Wirkus. The mathematical and theoretical biology institute - a model of mentorship through research. Mathematical Biosciences & Engineering, 2013, 10 (5&6) : 1351-1363. doi: 10.3934/mbe.2013.10.1351 [12] Avner Friedman, Wenrui Hao. Mathematical modeling of liver fibrosis. Mathematical Biosciences & Engineering, 2017, 14 (1) : 143-164. doi: 10.3934/mbe.2017010 [13] Ana I. Muñoz, José Ignacio Tello. Mathematical analysis and numerical simulation of a model of morphogenesis. Mathematical Biosciences & Engineering, 2011, 8 (4) : 1035-1059. doi: 10.3934/mbe.2011.8.1035 [14] Jinhae Park, Feng Chen, Jie Shen. Modeling and simulation of switchings in ferroelectric liquid crystals. Discrete & Continuous Dynamical Systems - A, 2010, 26 (4) : 1419-1440. doi: 10.3934/dcds.2010.26.1419 [15] Xiaoli Yang, Jin Liang, Bei Hu. Minimization of carbon abatement cost: Modeling, analysis and simulation. Discrete & Continuous Dynamical Systems - B, 2017, 22 (7) : 2939-2969. doi: 10.3934/dcdsb.2017158 [16] Michael Herty. Modeling, simulation and optimization of gas networks with compressors. Networks & Heterogeneous Media, 2007, 2 (1) : 81-97. doi: 10.3934/nhm.2007.2.81 [17] Gang Bao. Mathematical modeling of nonlinear diffracvtive optics. Conference Publications, 1998, 1998 (Special) : 89-99. doi: 10.3934/proc.1998.1998.89 [18] Edward J. Allen. Derivation and computation of discrete-delay and continuous-delay SDEs in mathematical biology. Mathematical Biosciences & Engineering, 2014, 11 (3) : 403-425. doi: 10.3934/mbe.2014.11.403 [19] Heiko Enderling, Alexander R.A. Anderson, Mark A.J. Chaplain, Glenn W.A. Rowe. Visualisation of the numerical solution of partial differential equation systems in three space dimensions and its importance for mathematical models in biology. Mathematical Biosciences & Engineering, 2006, 3 (4) : 571-582. doi: 10.3934/mbe.2006.3.571 [20] John D. Nagy. The Ecology and Evolutionary Biology of Cancer: A Review of Mathematical Models of Necrosis and Tumor Cell Diversity. Mathematical Biosciences & Engineering, 2005, 2 (2) : 381-418. doi: 10.3934/mbe.2005.2.381

2018 Impact Factor: 1.313