Next Article in Journal
Afatinib Overcomes Pemetrexed-Acquired Resistance in Non-Small Cell Lung Cancer Cells Harboring an EML4-ALK Rearrangement
Next Article in Special Issue
N-Cadherin Is Critical for the Survival of Germ Cells, the Formation of Steroidogenic Cells, and the Architecture of Developing Mouse Gonads
Previous Article in Journal
Cardiac Progenitor Cells from Stem Cells: Learning from Genetics and Biomaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Model Explaining the Role of CDC6 in the Diauxic Growth of CDK1 Activity during the M-Phase of the Cell Cycle

by
Mateusz Dębowski
1,
Zuzanna Szymańska
2,3,
Jacek Z. Kubiak
4,5,* and
Mirosław Lachowicz
1,*
1
Faculty of Mathematics, Informatics and Mechanics, Institute of Applied Mathematics and Mechanics, University of Warsaw, ul. Banacha 2, 02-097 Warsaw, Poland
2
Institute of Mathematics, Polish Academy of Sciences, ul. Śniadeckich 8, 00-656 Warsaw, Poland
3
ICM, University of Warsaw, ul. Tyniecka 15/17, 02-630 Warsaw, Poland
4
CNRS, Institute of Genetics and Development of Rennes, Univ Rennes, UMR 6290, Cell Cycle Group, Faculty of Medicine, F-35000 Rennes, France
5
Laboratory of Regenerative Medicine and Cell Biology, Military Institute of Hygiene and Epidemiology (WIHE), ul. Kozielska 4, 01-163 Warsaw, Poland
*
Authors to whom correspondence should be addressed.
Submission received: 6 November 2019 / Revised: 22 November 2019 / Accepted: 25 November 2019 / Published: 28 November 2019
(This article belongs to the Special Issue Manufacturing a Female Gamete: An Oocyte Story)

Abstract

:
In this paper we propose a role for the CDC 6 protein in the entry of cells into mitosis. This has not been considered in the literature so far. Recent experiments suggest that CDC 6 , upon entry into mitosis, inhibits the appearance of active CDK 1 and cyclin B complexes. This paper proposes a mathematical model which incorporates the dynamics of kinase CDK 1 , its regulatory protein cyclin B, the regulatory phosphatase CDC 25 and the inhibitor CDC 6 known to be involved in the regulation of active CDK 1 and cyclin B complexes. The experimental data lead us to formulate a new hypothesis that CDC 6 slows down the activation of inactive complexes of CDK 1 and cyclin B upon mitotic entry. Our mathematical model, based on mass action kinetics, provides a possible explanation for the experimental data. We claim that the dynamics of active complexes CDK 1 and cyclin B have a similar nature to diauxic dynamics introduced by Monod in 1949. In mathematical terms we state it as the existence of more than one inflection point of the curve defining the dynamics of the complexes.

1. Introduction

The mitotic cell cycle is an ordered sequence of events, grouped into four phases: G 1 , S , G 2 and M , during which the eukaryotic cell doubles its content, and divides into two daughter cells. The classical models for studies of cell cycle molecular machinery are oocytes and early embryos. These have the distinguishing property of being transcriptionally silent. This implies that the molecular machinery governing oocyte maturation and early embryo development is based on the maternal information accumulated during oocyte growth. While many different proteins regulate the progression of the cell cycle and the transitions between cell cycle phases, the major regulatory mechanisms are based on similar processes in all phases. Two main classes of proteins involved in cell cycle control are cyclins and enzymes called cyclin dependent kinases— CDKs . During individual phases a specific cyclin accumulates in the cell, associates with an appropriate kinase and with the help of other enzymes activates the kinase/cyclin complex. The appropriate level of an active complex triggers the transition to the next phase of the cell cycle.
A protein complex responsible for the transition from G 2 to M is formed from CDK 1 kinase and cyclin B, the latter being abbreviated as CYC B [3,4]. The complex of CDK 1 and CYC B , denoted by CDK 1 / CYC B , can be in one of two states, inactive or active, abbreviated as CDK 1 / CYC B N and CDK 1 / CYC B A , respectively [3]. CDK 1 / CYC B N is dephosphorylated by active CDC 25 phosphatase, denoted by CDC 25 A (inactive phosphatase CDC 25 is denoted by CDC 25 N ). CDC 25 A has the ability to pull away the phosphoryl group from two amino acids Tyr15 and Thr14 of the CDK 1 / CYC B N complex, making it active. This dephosphorylation results in the creation of CDK 1 / CYC B A . CDK 1 / CYC B A induces a cascade of phosphorylation of numerous substrates that change the character of cellular proteins from interphase to mitotic. These changes—necessary for mitotic progression—modify structures such as the cytoskeleton, membranes and DNA (condensation). Moreover, activation of CDC 25 phosphatase occurs due to its interaction with active CDK 1 / CYC B A complexes resulting in very powerful positive feedback between CDK 1 and CDC 25 that governs the CDK 1 activation upon the entry into M -phase.
Summarising, a subtle equilibrium between CDC 25 and CDK 1 enzymes is maintained at the beginning of the G 2 to M -phase transition process. The association of CDK 1 molecules with constantly synthesised CYC B results in the formation and accumulation of CDK 1 / CYC B N . This complex is inactive because it is phosphorylated on Tyr15/Thr14 since CDC 25 remains inactive. It was believed for a long time that, at some point, a spontaneous activation of the first molecules of CDK 1 / CYC B triggers a positive feedback between CDK 1 and CDC 25 . Active molecules of CDK 1 / CYC B start CDC 25 activation. Active CDC 25 , denoted as CDC 25 A , activates new molecules of CDK 1 / CYC B N . This in turn triggers a dramatic acceleration of biochemical events leading to full activation of the whole pool of CDK 1 / CYC B complexes present in the cell. Moreover, recently Vigneron et al. [5] have shown that another complex containing CDK 1 , namely, CDK 1 / CYC A , triggers the activation of CDK 1 / CYC B .
A number of mathematical models have been proposed to describe and understand mitotic cell cycle progression. For instance, some have used the stochastic approach to capture non-deterministic aspects of the process and emphasised that noise is an important factor influencing cell dynamics [6,7,8,9,10,11,12]. The common approach, however, is the deterministic one that usually is based on systems of ordinary differential equations. Some authors considered large systems of that type in order to investigate transitions between cell cycle phases. Such systems describing the dynamics of different proteins and enzymes were usually analysed only numerically, e.g., [13,14,15,16,17,18,19,20]. Some models have considered the activation of CDK 1 , e.g., [21,22,23]. There are cell cycle models based not only on standard ordinary differential equations, but also on delay differential equations [24,25]. For example, Busenberg et al. [24] contains a rigorous analysis of the model. Interesting analytic results can be found in other papers [26,27,28,29,30]. We emphasise the analytic results obtained by Ferrell et al. who described the cell cycle using ordinary differential equations and proved the existence of oscillatory dynamics of the cell cycle, switch-like behaviour of activity of CDK 1 and a bistability system [31,32].
The purpose of our work is to deepen the understanding of the cell cycle process. We are particularly focused on the G 2 to M -phase transition process and we aim at investigating the role of the CDC 6 protein in entering into the M -phase of the cell cycle. CDC 6 , up to now, was known as an essential ATPase active in the S phase and responsible for the initiation of DNA replication.
Recent experiments made in Xenopus laevis one-cell embryo cell-free extract suggest that CDC 6 has an important role in the delay of G 2 to M -phase transition [1]. This experimental system, however, simplified, as no nuclear DNA is present and all genetic information, both proteins and RNA, are of maternal origin and were accumulated in oocytes before the embryo development was triggered, allowed us to analyse the core molecular mechanisms of the cell cycle. We formulate a new hypothesis that explains this delay in terms of diauxic-like activation of CDK 1 / CYC B . In a different context, diauxic behaviour was introduced by Monod [2] in 1949. We propose a new mathematical model that captures this new hypothesis. We present the analysis and numerical simulations of this new model, suggesting how CDC 6 regulates the dynamics of CDK 1 / CYC B activation upon M -phase entry.
The structure of the paper is as follows. In Section 2 we present the material and methods used to obtain experimental results that inspired our subsequent research. In Section 3, at first, we introduce the biochemical model describing the basic events occurring when a cell enters mitosis. Then, we present the experimental results on the CDC 6 protein that motivate our work. We show both the data reprinted from El Dika et al. [1] in Figure 1 and original results in Figure 2. We formulate a new hypothesis that captures the role of CDC 6 in the process and the mathematical model corresponding to the biochemical one. Next, we present the numerical simulations of the proposed model and finally, Section 4 provides conclusions and directions for further research. In Appendix A we present the mathematical analysis of the presented model.

2. Material and Methods

2.1. Egg Collection and Activation

Spawn Xenopus laevis eggs were dejellied with 2% l-cysteine pH 7.81 in XB buffer (100 mM KCl, 1 mM MgCl2, 50 mM CaCl2, 10 mM HEPES and 50 mM sucrose pH 7.6). Next, they were washed in XB buffer, activated with 0.5 mg/mL calcium ionophore A23187 and extensively washed in XB.

2.2. Cell Free Extracts

Cytoplasmic extracts from calcium ionophore-activated one-cell embryos before the first embryonic mitosis were prepared according to El Dika et al. [1]. In short, embryos were cultured at 21 °C in XB buffer for 60–70 min postactivation, transferred into 5 mL ultraclearTM centrifuge tubes (Beckman Coulter, Roissy, France) in 0.5 mL of XB buffer containing 0.1 mM AEBSF, a protease inhibitor, at 4°. They were subjected to three consecutive centrifugations: The first short spin to remove XB excess and pack the embryos, the second 10,000× g spin at 4 °C for 10 min to separate the cell-free fractions, and the final 10,000× g clarification spin of the supernatant at 4 °C for 10 min. The supernatant was then incubated at 21 °C. Aliquots were taken out every 4 min and stored at −80 °C.

2.3. C D K 1 Activity Measurements

Samples of cell-free extracts were diluted in MPF buffer supplemented with: 0.5 mM sodium orthovanadate, 5 μg/μL of leupeptin, aprotinin, pepstatin and chymostatin, 0.4 mg/mL H1 histone (type III-S), 1 μ Ci [ γ 32P] ATP (specific activity: 3000 Ci/mmol; Amersham Biosciences, UK) and 0.8 mM ATP. After incubation at 30 °C for 30 min, phosphorylation reactions were stopped by adding Laemmli sample buffer and heated at 85 °C for 5 min. Histone H1 was separated by SDS-PAGE and incorporated radioactivity was measured by autoradiography of the gel using a STORM phosphorimager (Amersham Biosciences, Buckinghamshire, UK) followed by data analysis with ImageQuant 5.2 software.

2.4. C D C 6 Immunodepletion

Immunodepletion of CDC 6 from egg extracts was carried out using AffiPrep Protein A beads (Sigma, USA) conjugated with the anti- CDC 6 or with the preimmune serum overnight in 4 °C; 200 mL of beads were washed four times with XB buffer (pH 7.6) and incubated with 400 mL of extracts. After 30 min of incubation at 4 °C, extracts were centrifuged, beads were removed and supernatant was recovered. Two consecutive runs of immunodepletion were required to remove 90% of CDC 6 , as shown in El Dika et al. [1].

3. Results

3.1. Biochemical Model and the New Hypothesis

CYC B concentration gradually increases during the G 2 phase (cf. Equation (5)). CYC B pairs with protein kinase CDK 1 creates an inactive (phosphorylated) complex— CDK 1 / CYC B N (cf. Equation (1)). Inactive complex CDK 1 / CYC B N upon its interaction with phosphatase CDC 25 A becomes activated, thus the concentration of active complexes CDK 1 / CYC B A increases (cf. Equation (2)). Conversely, complex CDK 1 / CYC B A activates phosphatase CDC 25 N causing the appearance of more CDC 25 A (cf. Equation (3)). Summarising, CDC 25 A and CDK 1 / CYC B A form a positive feedback loop. The M phase begins when the concentration of active CDK 1 / CYC B A exceeds the threshold value.
Recent experimental studies provoke intriguing questions about the role of the CDC 6 protein in slowing down the activation of CDK 1 / CYC B N complexes. Figure 1 shows the concentration of CDK 1 / CYC B A (from a biochemical point of view, simply the CDK 1 activity) obtained on the basis of molecular experiments in two cases: (a) With CDC 6 ; and (b) without CDC 6 (after removal of CDC 6 from the experimental system) [1]. In the experimental setting with CDC 6 , one can notice a slower increase in the concentration of CDK 1 / CYC B A . Therefore, in the experimental system with CDC 6 the entry into mitosis is delayed. Our main goal is to explain the role of CDC 6 in the observed phenomenon.
The diauxic growth of CDK 1 activity was clearly noticed in previous studies of Xenopus laevis one-cell embryo cell-free extracts ([33]: Figure 1A bottom, Figure 2A right, Figure 3A right, [34]: Figure 2A bottom and [28]: Figure 1V). Furthermore, in our own research, we always observed the same type of behaviour of CDK 1 ([35]: Figures 1A, 2A, 3A and 6A, [36]: Figures 2A, 3A, 6A and 7A,B and [37]: Figures 6A and 7B). Moreover, the diauxic growth of CDK 1 activity is not an artefact due to the cell-free system because it was also observed in individual Xenopus laevis one-cell embryos ([38]: Figure 1A); however, it is more clear in the vegetal hemisphere where the CDK 1 activation is delayed and proceeds with lower dynamics ([38]: Figure 1B). The precise dynamics of the diauxic growth of CDK 1 , i.e., inflection times and slope of the curve, varies form one experiment to another. For this reason, the average curves showing the dynamics of CDK 1 activation upon the M -phase entry in Xenopus laevis one-cell embryos do not preserve the diauxic character ([35]: Figure 4B), where the average curve of 16 independent experiments does not show any inflection points. In Figure 2 in the current paper we show two examples of the fast and slow growth of CDK 1 activity in two independent experiments illustrating this problem well. The average curve of these two experiments also does not show the inflection points clearly visible in each experimental curve (data not shown).
We hypothesise that CDC 6 binds to CDK 1 / CYC B N and creates a new CDK 1 / CYC B / CDC 6 complex, preventing CDK 1 / CYC B N from being activated by CDC 25 A phosphatase (cf. Equation (4)). The resulting CDK 1 / CYC B / CDC 6 complexes constantly break down into CDC 6 and CDK 1 / CYC B N that constantly associate again. The more CDK 1 / CYC B N accumulate in the cell, the more CDK 1 / CYC B N complexes are activated by residual CDC 25 A . The formation of CDK 1 / CYC B / CDC 6 prolongs a very slow increase in the appearance of CDK 1 / CYC B A complexes. A slowdown in CDK 1 / CYC B A increase is visible as the flattening of the curve in Figure 1a. The experimental data leads to a new hypothesis on the mutual interaction between CDC 6 and CDK 1 / CYC B N , which determines the dynamics of CDK 1 / CYC B A upon mitotic entry. Our mathematical model, based on the law of mass action, bolsters this hypothesis. We suggest that the dynamics of CDK 1 / CYC B A are similar to diauxic dynamics introduced by Monod [2]. In mathematical terms, we state it as the existence of more than one inflection point on the curve defining the dynamics of the complexes, cf. Figure 3. Indeed, in the present model, we observe three or four inflexion points.
The second part of our hypothesis is that the reaction speed of CDK 1 / CYC B N and CDC 6 binding depends on active CDK 1 / CYC B A in a switch like mode. This means that when the concentration of CDK 1 / CYC B A is less than the concentration value, then the reaction speed of CDK 1 / CYC B / CDC 6 formation is low. When the CDK 1 / CYC B A is higher than the threshold value the reaction speed becomes much faster resulting in a two-step CDK 1 activation visible in biological experiments as an inflection of the activation curve of CDK 1 . Experimental data show clearly that this activation depends on the presence of CDC 6 [1].
Summarising, we consider the biochemical model that takes into account eight species, the descriptions of which are provided in Table 1, whereas the scheme of their mutual interactions is provided in Figure 4.
We consider the following five reactions (i.e., Equations (1)–(5)).
CDK 1 + CYC B α 1 CDK 1 / CYC B N ,
CDK 1 / CYC B N + CDC 25 A α 2 CDK 1 / CYC B A + CDC 25 A ,
CDK 1 / CYC B A + CDC 25 N α 3 CDK 1 / CYC B A + CDC 25 A ,
hyp : CDK 1 / CYC B N + CDC 6 δ α 4 · f ( CDK 1 / CYC B A ) CDK 1 / CYC B / CDC 6 ,
β β · K c y c B CYC B
We want to emphasise that Equations (1)–(3) and (5) correspond to the current state of knowledge. Equation (4) reflects the new hypothesis and is our contribution to understanding the phenomenon. In summary, taking Equation (4) into account is the first part of the new hypothesis. The speed of Equation (4) described by the function f is the second part of the hypothesis.

3.2. Mathematical Model

Assuming mass action kinetics for Equations (1)–(5) we transform the biochemical model into the system of eight ordinary differential equations (ODEs) with the following notation
x = CDK 1 , x a = CDK 1 / CYC B A , x n = CDK 1 / CYC B N , y a = CDC 25 A , y n = CDC 25 N , z = CDC 6 , w = CDK 1 / CYC B / CDC 6 , c = CYC B .
We have
x ˙ = α 1 x c , x ˙ a = α 2 x n y a , x ˙ n = α 1 x c α 2 x n y a α 4 f ( x a ) x n z + δ w , y ˙ a = α 3 x a y n , y ˙ n = α 3 x a y n , z ˙ = α 4 f ( x a ) x n z + δ w , w ˙ = α 4 f ( x a ) x n z δ w , c ˙ = α 1 x c + β ( K CYC B c ) ,
where α 1 , α 2 , α 3 , α 4 , β , δ are positive parameters and f ( x ) = ω + ν x k v t h k + x k . The function f is the Hill function that describes switch-like behaviour, where ν is a positive coefficient, k is a Hill coefficient, v t h is the threshold value of the switch and ω is the basic rate when x a = 0 . The function f describes the reaction rate of CDK 1 / CYC B N associated with CDC 6 resulting in the formation of CDK 1 / CYC B / CDC 6 complexes (cf. Equation (4)). In the system of ordinary differential equations, Equation (6) appears in the 7th equation describing the dynamics of CDK 1 / CYC B / CDC 6 and, due to the law of mass action, in the 3rd and 6th equations describing the dynamics of CDK 1 / CYC B N and CDC 6 , respectively. The process of CDK 1 / CYC B / CDC 6 formation seems to be highly nonlinear and we assume its rate to be CDK 1 / CYC B A dependent. There exists a similar mechanism governing interactions between CDK 1 and CDC 6 in S -phase. If CDK 1 / CYC B A is low then the majority of CDC 6 is not phosphorylated. However, with an increase of CDK 1 / CYC B A more phosphorylated CDC 6 appears in the cell [39]. The function f is bounded as it plays the role of a rate coefficient. The typical way of modelling such a nonlinear dependence is based on the Hill function, see, e.g., [40].
Taking into consideration the biological constraints, we propose the following initial data
x ( 0 ) = K CDK 1 ε x a ε x n ε w > 0 , x a ( 0 ) = ε x a K CDK 1 , x n ( 0 ) = ε x n K CDK 1 , y a ( 0 ) = ε y a K CDC 25 , y n ( 0 ) = K CDC 25 ε y a , z ( 0 ) = K CDC 6 ε w , w ( 0 ) = ε w min { K CDK 1 , K CDC 6 } , c ( 0 ) = 0 .
Equation (6) have the following conservation laws
x a + x n + x + w = K CDK 1 , y a + y n = K CDC 25 , z + w = K CDC 6 ,
where K CDK 1 , K CDC 25 , K CDC 6 denote constants given at the initial time. In Appendix A we provide the mathematical analysis of the model.
We provide the standard non-dimensionalisation of Equation (6). In other words we relate all considered variables to their characteristic values. With the substitution
x * = x K CDK 1 , x a * = x a K CDK 1 , x n * = x n K CDK 1 , w * = w K CDK 1 , y a * = y a K CDC 25 , y n * = y n K CDC 25 , z * = z K CDK 1 , c * = c K CYC B , t * = β t , γ = K CDC 6 K CDK 1 , ν t h * = ν t h K C D K 1 , δ * = δ β , α 1 * = α 1 K CYC B β , α 2 * = α 2 K CDC 25 β , α 3 * = α 3 K CDK 1 β , α 4 * = α 4 K CDK 1 β , ε x a * = ε x a K CDK 1 , ε x n * = ε x n K CDK 1 , ε y a * = ε y a K CDC 25 , ε w * = ε w K CDK 1 ,
and omitting the stars for simplicity, we obtain
x ˙ = α 1 x c , x ˙ a = α 2 x n y a , x ˙ n = α 1 x c α 2 x n y a α 4 f ( x a ) x n z + δ w , y ˙ a = α 3 x a y n , y ˙ n = α 3 x a y n , z ˙ = α 4 f ( x a ) x n z + δ w , w ˙ = α 4 f ( x a ) x n z δ w , c ˙ = α 1 x c + ( 1 c ) .
By Equation (8) it follows
x a + x n + x + w = 1 , y a + y n = 1 , z + w = γ .
By Equation (7) we obtain
x ( 0 ) = 1 ε x a ε x n ε w > 0 , x a ( 0 ) = ε x a 1 , x n ( 0 ) = ε x n 1 , y a ( 0 ) = ε y a 1 , y n ( 0 ) = 1 ε y a , z ( 0 ) = γ ε w , w ( 0 ) = ε w min { 1 , γ } , c ( 0 ) = 0 .
From the mathematical analysis presented in Appendix A we deduce that if the system contains even a small amount of CDK 1 / CYC B A or CDC 25 A then CDK 1 / CYC B A and CDC 25 A converge to full activation. This result is consistent with biological observations, because if the initial concentration of CDK 1 / CYC B A or CDC 25 A is positive then the positive feedback loop starts and the biological system tends to its equilibrium state (called S 2 ) defined by the maximal concentrations of CDK 1 / CYC B A and CDC 25 A . If the initial concentrations of CDK 1 / CYC B A or CDC 25 A are equal to zero, then the positive feedback loop does not start and the biological system tends to another equilibrium state (called S 1 ) defined by the concentrations of CDK 1 / CYC B A and CDC 25 A equal to 0. Small perturbations of the initial concentrations from zero to positive values change the equilibrium points, and this is the biological reason for S 1 being unstable and S 2 being asymptotically stable.
We note that a further simplification of the reduced model, Equation (A3), considered in Appendix A is reasonable. For example taking x = 0 and c = 1 we may reduce this system to a system of three equations
x ˙ a = α 2 ( 1 γ + z x a ) y a , y ˙ a = α 3 x a ( 1 y a ) , z ˙ = α 4 f ( x a ) ( 1 γ + z x a ) z + δ ( γ z ) .

3.3. Numerical Simulations

To carry out the numerical simulations we use the Runge-Kutta 4th order method provided by Matlab. Parameters values used to carry out the numerical simulations are given in Table 2. Figure 5 shows the concentrations of CDK 1 , CDK 1 / CYC B A , CDK 1 / CYC B N , CDK 1 / CYC B / CDC 6 . The most interesting curve is CDK 1 / CYC B A , where we observe three inflection points. The concentrations of species containing CDK 1 are shown in Figure 6 with the concentration of CDC 6 set to zero. Figure 7 shows the difference in activation: The timing and dynamics of the activation of CDK 1 / CYC B A in the presence and absence of CDC 6 . When CDC 6 is present the activation has more than one inflection point, and mitosis starts later, whereas when CDC 6 is absent, the activation is fast. In Figure 5 and Figure 7 we observe diauxic-type behaviour for the curve of CDK 1 / CYC B A . According to our hypothesis, this is related to the mutual interaction between CDC 6 and CDK 1 / CYC B N . We link this kind of behaviour with the existence of multiple (three or four in this case) inflection points in the curve of CDK 1 / CYC B A . The rigorous investigation of this fact leads to the analysis of behaviour of the second derivative of CDK 1 / CYC B A and more precisely its number of zeros. Figure 8 presents the graphs of the second derivative of x a obtained for the reduced system of three equations, Equation (13).
We may note that the numerical result given in Figure 8 has a rigorous nature as is visible by a careful estimation of the error of Matlab approximation. For example, considering Equation (13), we may provide (following the idea of [41]) the detailed analysis of the error
E n = y n y ( t n ) ,
where y ( t n ) is the value of the true solution at point t n and y n is the approximation of the solution at point t n , showing that
| E n | < 4.5 · 10 10 · h 4 .
Taking h = 1 1000 the error is smaller than the variation in the second derivative. Moreover, the rounding error can be neglected because the machine epsilon (see [42]) is sufficiently small compared with h. The details of the estimation are not reported here. This leads to the conclusion that the result stating the number of zeros of the second derivative has a rigorous nature and there the number of inflection points of the variable x a is either three or four, which give diauxic behaviour.

4. Discussion

The proposed model captures the most important characteristics of the diauxic growth of CDK 1 activation observed in biochemical experiments. Based on our previous experimental results [1] we claim that CDC 6 is the most important factor which causes the inflection of the CDK 1 activation curve. We have shown for the first time that CDC 6 is an inhibitory protein acting on CDK 1 during M -phase. Making use of our modelling setting, we hypothesise that CDC 6 binds to CDK 1 / CYC B N forming CDK 1 / CYC B / CDC 6 . CDK 1 / CYC B / CDC 6 formation results in slower activation of CDK 1 and consequently a delayed entry into M -phase. From a biochemical perspective our results, both experimental and modelling, are particularly interesting because the inhibitory effect of CDC 6 on CDK 1 activation during M -phase was not shown previously.
The second part of our hypothesis stands for the switch-like dependence of the reaction rate of CDK 1 / CYC B N binding to CDC 6 resulting in the formation of CDK 1 / CYC B / CDC 6 (Equation (4)). Our assumption that the mentioned reaction rate depends on CDK 1 / CYC B A provides a good qualitative explanation of the observed diauxic dynamics of CDK 1 activation. Further biological research is needed to investigate what molecular modification is necessary for this switch-like pattern of CDK 1 activation. We can postulate that CDK 1 , at some threshold, phosphorylates CDC 6 triggering the abrupt increase in CDC 6 affinity to CDK 1 / CYC B N . From a more general perspective, the slow rate of CDK 1 activation is very likely important for the physiological course of mitotic processes such as chromatin condensation or spindle formation in such a large cell as the Xenopus laevis one-cell embryo.
One of the main goals of the paper was to describe mathematically the diauxic behaviour of CDK 1 activation in the presence of CDC 6 protein. However, this kind of approach has a wide spectrum of use and may be applied to a large variety of problems. Usually such complex biological systems are very difficult to treat rigorously from a mathematical point of view. The analysis of the error gives a chance for a rigorous statement based on the numerical simulations. We leave the details of this approach for a forthcoming paper.
Our results may affect the understanding of the process of cancerogenesis since CDC 6 and its interactions with CDK 1 play an important role in mitotic regulation and in cancer etiology [43,44]. The CDC 6 role in M -phase regulation is not limited only to the mitotic cell cycle as shown in the current paper and in El Dika et al. [1], but also to the meiotic regulation in oocytes [45,46,47,48]. The requirement of CDC 6 for the meiotic spindle formation in mice and Xenopus laevis oocytes suggests that it can also be involved in mitotic spindle formation. The diauxic growth of CDK 1 activity determined by CDC 6 may be in relation to the proper dynamics of spindle assembly not only through the fine tuning of microtubule dynamics, but also by the proper coordination with other players like actin filaments [49,50].

Author Contributions

Conceptualization, M.D., Z.S., J.Z.K. and M.L.; formal analysis, M.D.; investigation, M.D., Z.S., J.Z.K. and M.L.; methodology, M.D., Z.S., J.Z.K. and M.L.; supervision, J.Z.K. and M.L.; writing–original draft, M.D. and Z.S.; writing–review and editing, M.D., Z.S., J.Z.K. and M.L.

Funding

M.D. was supported by a PhD scholarship of National Science Centre, Poland Grant 2017/25/B/ST1/00051. Z.S. was supported by the National Science Centre, Poland Grant 2017/26/M/ST1/00783. J.Z.K. was supported by the grant from Polish Ministry of National Defence “Kościuszko” 571/2016/DA. M.L. was supported by the National Science Centre, Poland Grant 2017/25/B/ST1/00051.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Mathematical Analysis of the Model

The global existence, uniqueness, positiveness and boundedness of the solutions follow directly from the form of Equation (10). Keeping in mind Equation (12), we obtain
x ( t ) , x a ( t ) , x n ( t ) , y a ( t ) , y n ( t ) , c ( t ) [ 0 , 1 ] , z ( t ) [ 0 , γ ] , w ( t ) [ 0 , min ( 1 , γ ) ] ,
for all t > 0 . From Equation (11) we obtain
x n = 1 x a x w , y n = 1 y a , w = γ z .
Using Equation (A2) we reduce Equation (10) to
x ˙ = α 1 x c , x ˙ a = α 2 ( 1 γ + z x a x ) y a , y ˙ a = α 3 x a ( 1 y a ) , z ˙ = α 4 f ( x a ) ( 1 γ + z x a x ) z + δ ( γ z ) c ˙ = α 1 x c + ( 1 c ) .
Referring to equilibrium points we formulate the following proposition.
Proposition A1.
In the set [ 0 , 1 ] 3 × [ 0 , γ ] × [ 0 , 1 ] there exist only two equilibrium points of Equation (A3):
(1) 
S 1 = 0 , 0 , 0 , 1 2 γ 1 δ α 4 ω + 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω , 1 ,
(2) 
S 2 = 0 , 1 , 1 , γ , 1 .
Proof. 
Let the equilibrium point of Equation (A3) be denoted by ( x ¯ , x ¯ a , y ¯ a , z ¯ , c ¯ ) . We have
α 1 x ¯ c ¯ = 0 , α 2 ( 1 γ + z ¯ x ¯ a x ¯ ) y ¯ a = 0 , α 3 x ¯ a ( 1 y ¯ a ) = 0 , α 4 f ( x ¯ a ) ( 1 γ + z ¯ x ¯ a x ¯ ) z ¯ + δ ( γ z ¯ ) = 0 , α 1 x ¯ c ¯ + ( 1 c ¯ ) = 0 .
From the fifth equation of Equation (A4) we obtain c ¯ 0 , and from the first and fifth equations of Equation (A4) we obtain x ¯ = 0 , c ¯ = 1 . Analysing the second equation of Equation (A4) we consider two cases
  • y ¯ a = 0 . Then from the third equation of Equation (A4) we get x ¯ a = 0 . Putting these results into the fourth equation of Equation (A4), with f ( 0 ) = ω , we obtain
    α 4 ω ( 1 γ + z ¯ ) z ¯ + δ ( γ z ¯ ) = 0 .
    Equation (A5) is a quadratic equation and has two solutions
    z ¯ 1 = 1 2 γ 1 δ α 4 ω 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω , z ¯ 2 = 1 2 γ 1 δ α 4 ω + 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω .
    Solution z ¯ 1 is negative because
    γ 1 δ α 4 ω 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω < < γ 1 δ α 4 ω | γ 1 δ α 4 ω | 0 .
    Solution z ¯ 2 is positive because
    γ 1 δ α 4 ω + 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω > > γ 1 δ α 4 ω + | γ 1 δ α 4 ω | 0 .
    We may note that z ¯ 2 γ . Taking into consideration Equation (A1) we obtain the equilibrium point
    x ¯ , x ¯ a , y ¯ a , z ¯ , c ¯ = 0 , 0 , 0 , 1 2 γ 1 δ α 4 ω + 1 γ + δ α 4 ω 2 + 4 δ γ α 4 ω , 1 .
  • y ¯ a 0 . From the second equation of Equation (A4) we get 1 γ + z ¯ x ¯ a x ¯ = 0 . Substituting this to the fourth equation of Equation (A4) we obtain z ¯ = γ . Now by x ¯ = 0 we have x ¯ a = 1 . Then from the third equation of Equation (A4) we obtain y ¯ a = 1 . We have the following equilibrium point
    x ¯ , x ¯ a , y ¯ a , z ¯ , c ¯ = 0 , 1 , 1 , γ , 1 .
Referring to the stability of the equilibrium points.
Proposition A2.
The equilibrium points S 1 and S 2 are unstable and asymptotically stable, respectively.
Proof. 
The Jacobi matrix J of Equation (A3) has the form
J = α 1 c 0 0 0 α 1 x α 2 y a α 2 y a α 2 ( 1 γ + z x a x ) α 2 y a 0 0 α 3 ( 1 y a ) α 3 x a 0 0 α 4 f ( x a ) z A 1 0 A 2 0 α 1 c 0 0 0 1 α 1 x ,
where
A 1 = α 4 f ( x a ) 1 γ + z x a x z + α 4 f ( x a ) z , A 2 = α 4 f ( x a ) 1 γ + z x a x α 4 f ( x a ) z δ .
We study the stability of the equilibrium point S 1 . For convenience we denote S 1 = ( 0 , 0 , 0 , z ¯ , 1 ) , where z ¯ = 1 2 γ 1 δ α 4 ω + 1 γ + δ α 4 ω 2 4 δ γ α 4 ω . The corresponding Jacobi matrix is
J ( S 1 ) = α 1 0 0 0 0 0 0 α 2 ( 1 γ + z ¯ ) 0 0 0 α 3 0 0 0 a 4 ω z ¯ α 4 ω z ¯ 0 α 4 ω ( 1 γ + 2 z ¯ ) δ 0 α 1 0 0 0 1 , ,
because f ( 0 ) = ω , f ( 0 ) = 0 . J ( S 1 ) is a block matrix, so we have
det J S 1 λ I = det P 1 λ I · det P 2 λ I · det P 3 λ I ,
where
P 1 = α 1 ,
P 2 = 0 α 1 ( 1 γ + z ¯ ) 0 α 3 0 0 α 4 ω z ¯ 0 α 4 ω 1 γ + 2 z ¯ δ ,
P 3 = 1 .
The characteristic polynomial w P 2 ( λ ) of matrix P 2 is
w P 2 ( λ ) = λ 2 α 4 ω ( 1 γ + 2 z ¯ ) δ λ + + α 4 ω ( 1 γ + 2 z ¯ ) + δ + λ α 1 α 3 ( 1 γ + z ¯ ) , w P 2 ( 0 ) = α 1 α 3 ( 1 γ + z ¯ ) α 4 ω ( 1 γ + 2 z ¯ ) + δ .
Now we examine the sign of 1 γ + z ¯ .
1 γ + z ¯ = 1 2 1 γ δ α 4 ω + γ 1 + δ α 4 ω 2 + 4 δ α 4 ω > > 1 2 1 γ δ α 4 ω + | γ 1 + δ α 4 ω | 0 ,
It follows that 1 γ + 2 z ¯ > 1 γ + z ¯ > 0 , then w P 2 ( 0 ) > 0 and w P 2 ( ) = . Therefore, there exits λ 0 > 0 such that w P 2 ( λ 0 ) = 0 . The matrix J ( S 1 ) has a positive eigenvalue; thus, the equilibrium point S 1 is unstable.
Next we study the stability of the point S 2 . The corresponding Jacobi matrix is
J ( S 2 ) = α 1 0 0 0 0 α 2 α 2 0 α 2 0 0 0 α 3 0 0 α 4 f ( 1 ) γ α 4 f ( 1 ) γ 0 α 4 f ( 1 ) γ δ 0 α 1 0 0 0 1 .
J ( S 2 ) is a block matrix; thus, we have
det J S 1 λ I = det R 1 λ I det R 2 λ I det R 3 λ I ,
where
R 1 = α 1 ,
R 2 = α 2 0 α 2 0 α 3 0 α 4 f ( 1 ) γ 0 α 4 f ( 1 ) γ δ ,
R 3 = 1 .
The characteristic polynomial w R 2 of matrix R 2 is:
w R 2 ( λ ) = λ 3 + v 2 λ 2 + v 1 λ α 2 α 3 δ ,
where
v 1 = ( α 3 α 4 f ( 1 ) γ + α 2 δ + α 3 δ + α 2 α 3 ) < 0 , v 2 = ( α 4 f ( 1 ) γ + δ + α 2 + α 3 ) < 0 .
Each coefficient of polynomial w R 2 is negative, so each root has a negative real part. Matrices R 1 and R 3 have eigenvalues α 1 and 1 , respectively, which are negative. Each eigenvalue of matrix J ( S 2 ) has a negative real part, so the stationary state S 2 is asymptotically stable. □
We may note that from Equation (A3) we easily obtain the following remark.
Remark A1.
(a)
If x a ( 0 ) = 0 = y a ( 0 ) , then x ( t ) , x a ( t ) , y a ( t ) , z ( t ) , c ( t ) S 1 as t .
(b)
If x a ( 0 ) > 0 or y a ( 0 ) > 0 , then x ( t ) , x a ( t ) , y a ( t ) , z ( t ) , c ( t ) S 2 as t .

References

  1. El Dika, M.; Laskowska-Kaszub, K.; Koryto, M.; Dudka, D.; Prigent, C.; Tassan, J.P.; Kloc, M.; Polanski, Z.; Borsuk, E.; Kubiak, J.Z. CDC6 controls dynamics of the first embryonic M-phase entry and progression via CDK1 inhibition. Dev. Biol. 2014, 396, 67–80. [Google Scholar] [CrossRef] [PubMed]
  2. Monod, J. The Growth of Bacterial Cultures. Annu. Rev. Microbiol. 1949, 3, 371–394. [Google Scholar] [CrossRef]
  3. Alberts, B.; Bray, D.; Hopkin, K.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Essential Cell Biology, 4th ed.; Garland Science/Taylor & Francis Group: New York, NY, USA, 2014. [Google Scholar]
  4. Morgan, D.O. The Cell Cycle: Principles of Control; New Science Press: London, UK, 2007. [Google Scholar]
  5. Vigneron, S.; Sundermann, L.; Labbé, J.C.; Pintard, L.; Radulescu, O.; Castro, A.; Lorca, T. Cyclin A-Cdk1-dependent phosphorylation of Bora is the triggering factor promoting mitotic entry. Dev. Cell 2018, 45, 637–650. [Google Scholar] [CrossRef] [PubMed]
  6. Ball, D.A.; Adames, N.R.; Reischmann, N.; Barik, D.; Franck, C.T.; Tyson, J.J.; Peccoud, J. Measurement and modeling of transcriptional noise in the cell cycle regulatory network. Cell Cycle 2013, 12, 3203–3218. [Google Scholar] [CrossRef] [PubMed]
  7. Braunewell, S.; Bornholdt, S. Superstability of the yeast cell-cycle dynamics: Ensuring causality in the presence of biochemical stochasticity. J. Theor. Biol. 2007, 245, 638–643. [Google Scholar] [CrossRef]
  8. Ge, H.; Qian, H.; Qian, M. Synchronized dynamics and nonequilibrium steady states in a stochastic yeast cell-cycle network. Math. Biosci. 2008, 211, 132–152. [Google Scholar] [CrossRef]
  9. Kar, S.; Baumann, W.T.; Paul, M.R.; Tyson, J.J. Exploring the roles of noise in the eukaryotic cell cycle. Proc. Natl. Acad. Sci. USA 2009, 106, 6471–6476. [Google Scholar] [CrossRef]
  10. Mura, I.; Csikász-Nagy, A. Stochastic Petri Net extension of a yeast cell cycle model. J. Theor. Biol. 2008, 254, 850–860. [Google Scholar] [CrossRef]
  11. Okabe, Y.; Sasai, M. Stable stochastic dynamics in yeast cell cycle. Biophys. J. 2007, 93, 3451–3459. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Qian, M.; Ouyang, Q.; Deng, M.; Li, F.; Tang, C. Stochastic model of the yeast cell-cycle network. Physica D 2006, 219, 35–39. [Google Scholar] [CrossRef]
  13. Aguda, B.D.; Tang, Y. The kinetics origins of the restriction point in the mammalian cell cycle. Cell Prolif. 1999, 32, 321–335. [Google Scholar] [CrossRef] [PubMed]
  14. Borisuk, M.T.; Tyson, J.J. Bifurcation analysis of a model of mitotic control in frog eggs. J. Theor. Biol. 1998, 195, 69–85. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, K.C.; Calzone, L.; Csikász-Nagy, A.; Cross, F.R.; Novak, B.; Tyson, J.J. Integrative analysis of cell cycle control in budding yeast. Mol. Biol. Cell 2004, 15, 3841–3862. [Google Scholar] [CrossRef] [PubMed]
  16. Ciliberto, A.; Novak, B.; Tyson, J.J. Mathematical model of the morphogenesis checkpoint in budding yeast. J. Cell Biol. 2003, 163, 1243–1254. [Google Scholar] [CrossRef]
  17. Li, B.; Shao, B.; Yu, C.; Ouyang, Q.; Wang, H. A mathematical model for cell size control in fission yeast. J. Theor. Biol. 2010, 264, 771–781. [Google Scholar] [CrossRef]
  18. Novak, B.; Tyson, J.J. Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos. J. Cell Sci. 1993, 106, 1153–1168. [Google Scholar]
  19. Novak, B.; Tyson, J.J. Modeling the cell division cycle: M-phase trigger, oscillations, and size control. J. Theor. Biol. 1993, 165, 101–134. [Google Scholar] [CrossRef]
  20. Novak, B.; Tyson, J.J. Modeling the control of DNA replication in fission yeast. Proc. Natl. Acad. Sci. USA 1997, 94, 9147–9152. [Google Scholar] [CrossRef]
  21. Charvin, G.; Oikonomou, C.; Siggia, E.D.; Cross, F.R. Origin of irreversibility of cell cycle start in budding yeast. PLoS Biol. 2010, 8, e1000284. [Google Scholar] [CrossRef]
  22. Goldbeter, A. A minimal cascade model for the mitotic oscillator involving cyclin and cdc2 kinase. Proc. Natl. Acad. Sci. USA 1991, 88, 9107–9111. [Google Scholar] [CrossRef]
  23. Goldbeter, A.; Guilmot, J.M. Arresting the mitotic oscillator and the control of cell proliferation: Insights from a cascade model for cdc2 kinase activation. Experientia 1996, 52, 212–216. [Google Scholar] [CrossRef] [PubMed]
  24. Busenberg, S.; Tang, B. Mathematical models of the early embryonic cell cycle: The role of MPF activation and cyclic degradation. J. Math. Biol. 1994, 32, 573–596. [Google Scholar] [CrossRef] [PubMed]
  25. Srividhya, J.; Gopinathan, M.S. A simple time delay model for eukaryotic cell cycle. J. Theor. Biol. 2006, 241, 617–627. [Google Scholar] [CrossRef] [PubMed]
  26. Angeli, D.; Ferrell, J.E., Jr.; Sontag, E.D. Detection of multistability, bifurcation, and hysteresis in the large class of biological positive-feedback systems. Proc. Natl. Acad. Sci. USA 2004, 101, 1822–1827. [Google Scholar] [CrossRef] [PubMed]
  27. Chang, J.B.; Ferrell, J.E., Jr. Mitotic trigger waves and the spatial coordination of the Xenopus cell cycle. Nature 2013, 500, 603–607. [Google Scholar] [CrossRef]
  28. Pomerening, J.R.; Kim, S.Y.; Ferrell, J.E., Jr. System-level dissection of the cell cycle oscillator: Bypassing positive feedback produces damped oscillations. Cell 2005, 122, 565–578. [Google Scholar] [CrossRef]
  29. Pomerening, J.R.; Sontag, E.D.; Ferrell, J.E., Jr. Building a cell cycle oscillator: hysteresis and bistability in the activation of Cdc2. Nat. Cell Biol. 2003, 5, 346–351. [Google Scholar] [CrossRef]
  30. Trunnell, N.B.; Poon, A.C.; Kim, S.Y.; Ferrell, J.E., Jr. Ultrasensitivity in the regulation of Cdc25C by Cdk1. Mol. Cell 2011, 41, 263–274. [Google Scholar] [CrossRef]
  31. Ferrell, J.E., Jr.; Pomerening, J.R.; Kim, S.Y.; Trunnell, N.B.; Xiong, W.; Huang, C.Y.; Machleder, E.M. Simple, realistic models of complex biological processes: Positive feedback and bistability in a cell fate switch and a cell cycle oscillator. FEBS Lett. 2009, 583, 3999–4005. [Google Scholar] [CrossRef]
  32. Ferrell, J.E., Jr.; Tsai, T.Y.; Yang, Q. Modeling the cell cycle: Why do certain circuits oscillate? Cell 2011, 144, 874–885. [Google Scholar] [CrossRef]
  33. Walter, S.A.; Guadagno, T.M.; Ferrell, J.E., Jr. Induction of a G2-phase arrest in Xenopus egg extracts by activation of p42 mitogen-activated protein kinase. Mol. Biol. Cell 1997, 8, 2157–2169. [Google Scholar] [CrossRef] [PubMed]
  34. Bitangcol, J.C.; Chau, A.S.; Stadnick, E.; Lohka, M.J.; Dicken, B.; Shibuya, E.K. Activation of the p42 mitogen-activated protein kinase pathway inhibits Cdc2 activation and entry into M-phase in cycling Xenopus egg extracts. Mol. Biol. Cell 1998, 9, 451–467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Chesnel, F.; Vignaux, F.; Richard-Parpaillon, L.; Huguet, A.; Kubiak, J.Z. Differences in regulation of the first two M-phases in Xenopus laevis embryo cell-free extracts. Dev. Biol. 2005, 285, 358–375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Chesnel, F.; Bazile, F.; Pascal, A.; Kubiak, J.Z. Cyclin B dissociation from CDK1 precedes its degradation upon MPF inactivation in mitotic extracts of Xenopus laevis embryos. Cell Cycle 2006, 5, 1687–1698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Bazile, F.; Pascal, A.; Karaiskou, A.; Chesnel, F.; Kubiak, J.Z. Absence of reciprocal feedback between MPF and ERK2 MAP kinase in mitotic Xenopus laevis embryo cell-free extract. Cell Cycle 2007, 6, 489–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Pérez-Mongiovi, D.; Chang, P.; Houliston, E. A propagated wave of MPF activation accompanies surface contraction waves at first mitosis in Xenopus. J. Cell Sci. 1998, 111, 385–393. [Google Scholar]
  39. Al-Zain, A.; Schroeder, L.; Sheglov, A.; Ikui, A.E. Cdc6 degradation requires phosphodegron created by GSK-3 and Cdk1 for SCFCdc4 recognition in Saccharomyces cerevisiae. Mol. Biol. Cell 2015, 26, 2609–2619. [Google Scholar] [CrossRef]
  40. Monod, J.; Wyman, J.; Changeux, J.P. On the nature of allosteric transition: A plausible model. J. Mol. Biol. 1965, 12, 88–118. [Google Scholar] [CrossRef]
  41. Lotkin, M. On the accuracy of Runge Kutta’s method. Math. Tables Other Aids Comput. 1951, 5, 128–133. [Google Scholar] [CrossRef]
  42. Atkinson, K.; Han, W.; Stewart, D. Numerical Solution of Ordinary Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  43. Kalfalah, F.M.; Berg, E.; Christensen, M.O.; Linka, R.M.; Dirks, W.G.; Boege, F.; Mielke, C. Spatio-temporal regulation of the human licensing factor Cdc6 in replication and mitosis. Cell Cycle 2015, 14, 1704–1715. [Google Scholar] [CrossRef] [Green Version]
  44. Borlado, L.R.; Méndez, J. CDC6: From DNA replication to cell cycle checkpoints and oncogenesis. Carcinogenesis 2008, 29, 237–243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Anger, M.; Stein, P.; Schultz, R.M. CDC6 requirement for spindle formation during maturation of mouse oocytes. Biol. Reprod. 2005, 72, 188–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Narasimhachar, Y.; Webster, D.R.; Gard, D.L.; Coué, M. Cdc6 is required for meiotic spindle assembly in Xenopus oocytes. Cell Cycle 2012, 11, 524–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Daldello, E.M.; Le, T.; Poulhe, R.; Jessus, C.; Haccard, O.; Dupré, A. Control of Cdc6 accumulation by Cdk1 and MAPK is essential for completion of oocyte meiotic divisions in Xenopus. J. Cell Sci. 2015, 128, 2482–2496. [Google Scholar] [CrossRef] [Green Version]
  48. Daldello, E.M.; Le, T.; Poulhe, R.; Jessus, C.; Haccard, O.; Dupré, A. Correction: Control of Cdc6 accumulation by Cdk1 and MAPK is essential for completion of oocyte meiotic divisions in Xenopus. J. Cell Sci. 2018, 131, 2482–2496. [Google Scholar] [CrossRef] [Green Version]
  49. Roeles, J.; Tsiavaliaris, G. Actin-microtubule interplay coordinates spindle assembly in human oocytes. Nat. Commun. 2019, 10, 1–10. [Google Scholar] [CrossRef] [Green Version]
  50. Mogessie, B.; Scheffler, K.; Schuh, M. Assembly and positioning of the oocyte meiotic spindle. Annu. Rev. Cell Dev. Biol. 2018, 34, 381–403. [Google Scholar] [CrossRef] [Green Version]
Figure 1. CDK 1 / CYC B A activity in the control extract containing physiological amounts of CDC 6 (a) and in the extract immunodepleted of CDC 6 (b). Note a slow and diauxic growth of CDK 1 / CYC B A activity in the control extract (a) and the very rapid activation in the absence of CDC 6 (b). Curves reprinted from El Dika et al. [1].
Figure 1. CDK 1 / CYC B A activity in the control extract containing physiological amounts of CDC 6 (a) and in the extract immunodepleted of CDC 6 (b). Note a slow and diauxic growth of CDK 1 / CYC B A activity in the control extract (a) and the very rapid activation in the absence of CDC 6 (b). Curves reprinted from El Dika et al. [1].
Cells 08 01537 g001
Figure 2. Differences in dynamics of CDK 1 / CYC B A activation curves in control extracts containing physiological amounts of CDC 6 . Two extreme examples are shown: Rapid activation taking 16 min (a) and slow activation taking 28 min (b). Note that the inflection points of the curves appear at different moments in relation to the maximum activity.
Figure 2. Differences in dynamics of CDK 1 / CYC B A activation curves in control extracts containing physiological amounts of CDC 6 . Two extreme examples are shown: Rapid activation taking 16 min (a) and slow activation taking 28 min (b). Note that the inflection points of the curves appear at different moments in relation to the maximum activity.
Cells 08 01537 g002
Figure 3. The smoothed curve obtained on the basis of experimental data presented in Figure 1a. Red circles indicate approximate location of inflection points for the setting with CDC 6 upon M -phase entry.
Figure 3. The smoothed curve obtained on the basis of experimental data presented in Figure 1a. Red circles indicate approximate location of inflection points for the setting with CDC 6 upon M -phase entry.
Cells 08 01537 g003
Figure 4. The schematic diagram of the considered system. Colours of arrows and dots correspond to colours of Equations (1)–(5). For simplicity we do not consider the potential marginal separation of the complex CDK 1 / CYC B N into CDK 1 and CYC B . On the diagram we indicate this by “*”.
Figure 4. The schematic diagram of the considered system. Colours of arrows and dots correspond to colours of Equations (1)–(5). For simplicity we do not consider the potential marginal separation of the complex CDK 1 / CYC B N into CDK 1 and CYC B . On the diagram we indicate this by “*”.
Cells 08 01537 g004
Figure 5. Concentration of CDK 1 , CDK 1 / CYC B A , CDC 6 in the presence of CDC 6 .
Figure 5. Concentration of CDK 1 , CDK 1 / CYC B A , CDC 6 in the presence of CDC 6 .
Cells 08 01537 g005
Figure 6. Concentration of CDK 1 , CDK 1 / CYC B A , CDC 6 in the absence of CDC 6 .
Figure 6. Concentration of CDK 1 , CDK 1 / CYC B A , CDC 6 in the absence of CDC 6 .
Cells 08 01537 g006
Figure 7. Comparison between concentration of CDK 1 / CYC B A in the absence and presence of CDC 6 . Solid line—system with CDC 6 ; dotted line—system without CDC 6 .
Figure 7. Comparison between concentration of CDK 1 / CYC B A in the absence and presence of CDC 6 . Solid line—system with CDC 6 ; dotted line—system without CDC 6 .
Cells 08 01537 g007
Figure 8. Graphs presenting second derivatives of x a showing the number of zeros, which indicates the number of inflection points. (a) corresponds to the case with the second derivative starting from a positive value and having three zeros. (b) corresponds to the case with the second derivative starting from a negative value and having four zeros.
Figure 8. Graphs presenting second derivatives of x a showing the number of zeros, which indicates the number of inflection points. (a) corresponds to the case with the second derivative starting from a positive value and having three zeros. (b) corresponds to the case with the second derivative starting from a negative value and having four zeros.
Cells 08 01537 g008
Table 1. Species of the biochemical model involved in Equations (1)–(5).
Table 1. Species of the biochemical model involved in Equations (1)–(5).
SpeciesDescription
CDK 1 cyclin-dependent kinase 1
CYC B cyclin B
CDK 1 / CYC B A active complex of CDK and CYC B
CDK 1 / CYC B N inactive complex of CDK and CYC B
CDC 25 A active phosphatase CDC25
CDC 25 N inactive phosphatase CDC25
CDC 6 cell division cycle 6 ATPase
CDK 1 / CYC B / CDC 6 complex of CDK 1 / CYC B N and CDC 6
Table 2. Values of parameters used in simulations.
Table 2. Values of parameters used in simulations.
ParameterValueParameterValueParameterValue
α 1 1 δ 4 ε x a 0.001
α 2 30k20 ε x n 0.001
α 3 1 ν 8 ε y a 0
α 4 7 ν t h 0.25 ε w 0
γ 1 ω 0.6

Share and Cite

MDPI and ACS Style

Dębowski, M.; Szymańska, Z.; Kubiak, J.Z.; Lachowicz, M. Mathematical Model Explaining the Role of CDC6 in the Diauxic Growth of CDK1 Activity during the M-Phase of the Cell Cycle. Cells 2019, 8, 1537. https://0-doi-org.brum.beds.ac.uk/10.3390/cells8121537

AMA Style

Dębowski M, Szymańska Z, Kubiak JZ, Lachowicz M. Mathematical Model Explaining the Role of CDC6 in the Diauxic Growth of CDK1 Activity during the M-Phase of the Cell Cycle. Cells. 2019; 8(12):1537. https://0-doi-org.brum.beds.ac.uk/10.3390/cells8121537

Chicago/Turabian Style

Dębowski, Mateusz, Zuzanna Szymańska, Jacek Z. Kubiak, and Mirosław Lachowicz. 2019. "Mathematical Model Explaining the Role of CDC6 in the Diauxic Growth of CDK1 Activity during the M-Phase of the Cell Cycle" Cells 8, no. 12: 1537. https://0-doi-org.brum.beds.ac.uk/10.3390/cells8121537

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop