# 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] Shujing Shi, Jicai Huang, Yang Kuang. Global dynamics in a tumor-immune model with an immune checkpoint inhibitor. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1149-1170. doi: 10.3934/dcdsb.2020157 [2] Jean-Paul Chehab. Damping, stabilization, and numerical filtering for the modeling and the simulation of time dependent PDEs. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021002 [3] Linhao Xu, Marya Claire Zdechlik, Melissa C. Smith, Min B. Rayamajhi, Don L. DeAngelis, Bo Zhang. Simulation of post-hurricane impact on invasive species with biological control management. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 4059-4071. doi: 10.3934/dcds.2020038 [4] Michal Beneš, Pavel Eichler, Jakub Klinkovský, Miroslav Kolář, Jakub Solovský, Pavel Strachota, Alexandr Žák. Numerical simulation of fluidization for application in oxyfuel combustion. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 769-783. doi: 10.3934/dcdss.2020232 [5] Xin Guo, Lexin Li, Qiang Wu. Modeling interactive components by coordinate kernel polynomial models. Mathematical Foundations of Computing, 2020, 3 (4) : 263-277. doi: 10.3934/mfc.2020010 [6] Min Chen, Olivier Goubet, Shenghao Li. Mathematical analysis of bump to bucket problem. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5567-5580. doi: 10.3934/cpaa.2020251 [7] Philippe G. Ciarlet, Liliana Gratie, Cristinel Mardare. Intrinsic methods in elasticity: a mathematical survey. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 133-164. doi: 10.3934/dcds.2009.23.133 [8] M. Dambrine, B. Puig, G. Vallet. A mathematical model for marine dinoflagellates blooms. Discrete & Continuous Dynamical Systems - S, 2021, 14 (2) : 615-633. doi: 10.3934/dcdss.2020424 [9] Riadh Chteoui, Abdulrahman F. Aljohani, Anouar Ben Mabrouk. Classification and simulation of chaotic behaviour of the solutions of a mixed nonlinear Schrödinger system. Electronic Research Archive, , () : -. doi: 10.3934/era.2021002 [10] Chun Liu, Huan Sun. On energetic variational approaches in modeling the nematic liquid crystal flows. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 455-475. doi: 10.3934/dcds.2009.23.455 [11] Vieri Benci, Sunra Mosconi, Marco Squassina. Preface: Applications of mathematical analysis to problems in theoretical physics. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020446 [12] Urszula Ledzewicz, Heinz Schättler. On the role of pharmacometrics in mathematical models for cancer treatments. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 483-499. doi: 10.3934/dcdsb.2020213 [13] Jakub Kantner, Michal Beneš. Mathematical model of signal propagation in excitable media. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 935-951. doi: 10.3934/dcdss.2020382 [14] Hong Niu, Zhijiang Feng, Qijin Xiao, Yajun Zhang. A PID control method based on optimal control strategy. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 117-126. doi: 10.3934/naco.2020019 [15] Omid Nikan, Seyedeh Mahboubeh Molavi-Arabshai, Hossein Jafari. Numerical simulation of the nonlinear fractional regularized long-wave model arising in ion acoustic plasma waves. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020466 [16] Jun Zhou. Lifespan of solutions to a fourth order parabolic PDE involving the Hessian modeling epitaxial growth. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5581-5590. doi: 10.3934/cpaa.2020252 [17] Niklas Kolbe, Nikolaos Sfakianakis, Christian Stinner, Christina Surulescu, Jonas Lenz. Modeling multiple taxis: Tumor invasion with phenotypic heterogeneity, haptotaxis, and unilateral interspecies repellence. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 443-481. doi: 10.3934/dcdsb.2020284 [18] Huijuan Song, Bei Hu, Zejia Wang. Stationary solutions of a free boundary problem modeling the growth of vascular tumors with a necrotic core. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 667-691. doi: 10.3934/dcdsb.2020084 [19] Yining Cao, Chuck Jia, Roger Temam, Joseph Tribbia. Mathematical analysis of a cloud resolving model including the ice microphysics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 131-167. doi: 10.3934/dcds.2020219 [20] Martin Kalousek, Joshua Kortum, Anja Schlömerkemper. Mathematical analysis of weak and strong solutions to an evolutionary model for magnetoviscoelasticity. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 17-39. doi: 10.3934/dcdss.2020331

2018 Impact Factor: 1.313