The Novak–Tyson Model is a non-linear dynamics framework developed in the context of cell-cycle control by Bela Novak and John J. Tyson. It is a prevalent theoretical model that describes a hysteretic, bistable bifurcation of which many biological systems have been shown to express.
Historical background
editBela Novak and John Tyson came from the Department of Biology at the Virginia Polytechnic Institute and State University in Blacksburg, Virginia, when this model was first published in the Journal of Cell Science in 1993.[1]
In 1990, two key papers were published that identified and characterized important dynamic relationships between cyclin and MPF in interphase-arrested frog egg extracts. The first was Solomon's 1990 Cell paper, titled "Cyclin activation of p34cdc2" and the second was Felix's 1990 Nature paper, titled "Triggering of cyclin degradation in interphase extracts of amphibian eggs of cdc2 kinase".[2][3] Solomon's paper showed a distinct cyclin concentration threshold for the activation of MPF.[3] Felix's paper looked at cyclin B degradation in these extracts and found that MPF degrades cyclin B in a concentration dependent and time-delayed manner.[2]
In response to these observations, three competing models were published in the next year, 1991, by Norel and Agur, Goldbeter, and Tyson.[4][5][6] These competing theories all attempted to model the experimental observations seen in the 1990 papers regarding the cyclin-MPF network.
The Norel and Agur model
editNorel and Agur's model proposes a mechanism where cyclin catalytically drives the production of MPF, which in turn autocatalyzes.[4] This model assumes that MPF activates cyclin degradation via APC activation, and it decouples cyclin degradation from MPF destruction.[4] However, this model is unable to recreate the observed cyclin dependent MPF activity relationship seen in Solomon's 1990 paper, as it shows no upper steady-state level of MPF activity.[7]
Goldbeter model
editGoldbeter proposed a model where cyclin also catalytically activates MPF, but without an autocatalytic, positive feedback loop.[5] The model describes a two-step process, where MPF first activates the APC, and then the APC drives cyclin degradation.[5] When graphing the MPF activity with respect to cyclin concentration, the model shows a sigmoidal shape, with a hypersensitive, threshold region similar to what was observed by Solomon.[7] However, this model depicts an effectively asymptotic plateau behavior at cyclin concentrations above the threshold, whereas the observed curve shows a steady increase in MPF activity at cyclin concentrations above the threshold.[7]
Tyson model
editIn Tyson's 1991 model, cyclin is a stoichiometric activator of Cdc2, as cyclin binds with phosphorylated Cdc2 to form preMPF, which is activated by Cdc25 to generate MPF.[6] Because Cdc25 itself is also activated by MPF, the conversion of preMPF to active MPF is a self-amplifying process in this model.[6] Tyson neglected the role of MPF in activating the APC, assuming that only a phosphorylated form of cyclin was rapidly degraded.[7] Tyson's model predicts an S-shaped curve, which is phenotypically consistent with Solomon's experimental results. However, this model generates additional lower turning point behavior on the S-curve that implies hysteresis when interpreted as a threshold.[7]
The Novak–Tyson model, first published in the paper titled "Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos", builds on the Goldbeter and Tyson 1991 models in order to generate a unifying theory, encapsulating the observed dynamics of the cyclin-MPF relationship.[1]
Model
editThe model proposes a complex set of feedback relationships that are mathematically defined by a series of rate constants and ordinary differential equations. It employs concepts seen in the previous models such as stoichiometric binding of Cdc2 and cyclin B, positive feedback loops through Cdc25 and Wee1, and delayed activation by MPF of the APC, but includes additional reactions such as that of Wee1 and Cdc25.[7] The result is a non-linear dynamic system with a similar S-shaped curve from Tyson's 1991 model.[7] In the process, this model makes four key predictions.
Discontinuous bistable hysteresis
editAccording to the Novak–Tyson model, rather than describing Solomon's observations as a sigmoidal switch as seen in the Goldbeter model, the threshold behavior of cyclin concentration dependent MPF activity is instead, a discontinuity of a bistable system.[1] Moreover, due to the S-shape dynamics, the Novak–Tyson model additionally predicts that the cyclin concentration threshold for activation is higher than the cyclin concentration threshold for inactivation; that is, this model predicts a dynamically hysteretic behavior.[1]
Critical slowing down
editSince the Novak–Tyson model predicts that the observed threshold is actually a discontinuity in the system dynamics, it additionally predicts a critical slowing down behavior near the threshold, which is a characteristic behavior of discontinuous bistable systems.[1]
Biochemical regulation
editSince the model predicts that MPF activation at the interphase-to-mitosis transition is governed by the turning point of an S-shaped curve, Novak and Tyson suggest that transition-delaying checkpoint signals biochemically move the turning point to larger values of cyclin B concentration.[1]
Phosphatase regulation
editNovak and Tyson predict that unreplicated DNA interferes with M-phase initiation by activating the phosphatases that oppose MPF in the positive feedback loops.[1] This prediction suggests a possible role for regulated serine/threonine protein phosphatases in cell cycle control.[7]
Model validation
editAt the time of publishing, the predictions from the paper were all experimentally untested and were based only off the signal pathways and mathematical modeling proposed by Novak and Tyson. However, since then, two papers have experimentally validated three of the four predictions listed above, namely the discontinuous bistable hysteresis, critical slowing down, and biochemical regulation predictions.[8][9]
Weaknesses
editAccording to the Novak and Tyson, this model, as with any biologically detailed, mathematically driven model, is heavily reliant on parameter estimation, especially given the mathematical complexity for this particular model. Ultimately these parameters are fit to experimental data, which is inherently susceptible to the compounded reliability of various experiments measuring various parameters.
References
edit- ^ a b c d e f g Novak, B.; Tyson, J. J. (December 1993). "Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos". Journal of Cell Science. 106 ( Pt 4) (4): 1153–1168. doi:10.1242/jcs.106.4.1153. ISSN 0021-9533. PMID 8126097.
- ^ a b Félix, Marie-Anne; Labbé, Jean-Claude; Dorée, Marcel; Hunt, Tim; Karsenti, Eric (July 1990). "Triggering of cyclin degradation in interphase extracts of amphibian eggs by cdc2 kinase". Nature. 346 (6282): 379–382. Bibcode:1990Natur.346..379F. doi:10.1038/346379a0. ISSN 0028-0836. PMID 2142754. S2CID 59267.
- ^ a b Solomon, M. J.; Glotzer, M.; Lee, T. H.; Philippe, M.; Kirschner, M. W. (1990-11-30). "Cyclin activation of p34cdc2". Cell. 63 (5): 1013–1024. doi:10.1016/0092-8674(90)90504-8. ISSN 0092-8674. PMID 2147872. S2CID 26120437.
- ^ a b c Norel, R.; Agur, Z. (1991-03-01). "A model for the adjustment of the mitotic clock by cyclin and MPF levels". Science. 251 (4997): 1076–1078. Bibcode:1991Sci...251.1076N. doi:10.1126/science.1825521. ISSN 0036-8075. PMID 1825521.
- ^ a b c Goldbeter, A. (1991-10-15). "A minimal cascade model for the mitotic oscillator involving cyclin and cdc2 kinase". Proceedings of the National Academy of Sciences of the United States of America. 88 (20): 9107–9111. Bibcode:1991PNAS...88.9107G. doi:10.1073/pnas.88.20.9107. ISSN 0027-8424. PMC 52661. PMID 1833774.
- ^ a b c Tyson, J. J. (1991-08-15). "Modeling the cell division cycle: cdc2 and cyclin interactions". Proceedings of the National Academy of Sciences of the United States of America. 88 (16): 7328–7332. Bibcode:1991PNAS...88.7328T. doi:10.1073/pnas.88.16.7328. ISSN 0027-8424. PMC 52288. PMID 1831270.
- ^ a b c d e f g h Tyson, John J.; Novák, Béla (December 2015). "Models in biology: lessons from modeling regulation of the eukaryotic cell cycle". BMC Biology. 13 (1): 46. doi:10.1186/s12915-015-0158-9. ISSN 1741-7007. PMC 4486427. PMID 26129844.
- ^ Sha, Wei; Moore, Jonathan; Chen, Katherine; Lassaletta, Antonio D.; Yi, Chung-Seon; Tyson, John J.; Sible, Jill C. (2003-02-04). "Hysteresis drives cell-cycle transitions in Xenopus laevis egg extracts". Proceedings of the National Academy of Sciences of the United States of America. 100 (3): 975–980. doi:10.1073/pnas.0235349100. ISSN 0027-8424. PMC 298711. PMID 12509509.
- ^ Pomerening, Joseph R.; Sontag, Eduardo D.; Ferrell, James E. (April 2003). "Building a cell cycle oscillator: hysteresis and bistability in the activation of Cdc2". Nature Cell Biology. 5 (4): 346–351. doi:10.1038/ncb954. ISSN 1465-7392. PMID 12629549. S2CID 11047458.