Tyson1991 - Cell Cycle 2 var

  public model
Model Identifier
BIOMD0000000006
Short description
Tyson1991 - Cell Cycle 2 var

Mathematical model of the interactions of cdc2 and cyclin.

Description taken from the original Cellerator version of the model ( Tyson (1991, 2 variables) at http://www.cellerator.org ).

This model is described in the article:

Tyson JJ.
Proc. Natl. Acad. Sci. U.S.A. 1991; 88(16); 7328-32

Abstract:

The proteins cdc2 and cyclin form a heterodimer (maturation promoting factor) that controls the major events of the cell cycle. A mathematical model for the interactions of cdc2 and cyclin is constructed. Simulation and analysis of the model show that the control system can operate in three modes: as a steady state with high maturation promoting factor activity, as a spontaneous oscillator, or as an excitable switch. We associate the steady state with metaphase arrest in unfertilized eggs, the spontaneous oscillations with rapid division cycles in early embryos, and the excitable switch with growth-controlled division cycles typical of nonembryonic cells.

This is a two variable reduction of the larger 6-variable model published in the same paper. The equations are:

u'= k4(v-u)(alpha+u^2)-k6*u
v'=kappa-k6*u
z= v-u
with kappa = k1[aa]/[CT]

In the present implementation, an additional variable z is introduced with z = v-u is made, so that the different variables be interpreted as follows:

u=[activeMPF]/[CT]
v=([cyclin]+[preMPF]+[activeMPF])/[CT]
z=([ cyclin]+[preMPF])/[CT]
with [CT]=[CDC2]+{CDC2P]+[preMPF]+[aMPF].

The reactions included are only to show the flows between z and u, and do not influence the species, as they all are set to boundaryCondition=True , meaning, that they are only determined by the rate rules (explicit differential equations) and assignment rules.

If you set boundaryCondition=False and remove the rate rules for v, u and the the assignment rule for z, you get the more symmetrical, but equivalent, version from the Cellerator repository:

u'= k4*z*(alpha+u^2)-k6*u
z'=kappa-z*(alpha+u^2)

To the extent possible under law, all copyright and related or neighbouring rights to this encoded model have been dedicated to the public domain worldwide. Please refer to CC0 Public Domain Dedication for more information.

Format
SBML (L2V4)
Related Publication
  • Modeling the cell division cycle: cdc2 and cyclin interactions.
  • Tyson JJ
  • Proceedings of the National Academy of Sciences of the United States of America , 8/ 1991 , Volume 88 , pages: 7328-7332 , PubMed ID: 1831270
  • Department of Biology, Virginia Polytechnic Institute and State University, Blacksburg 24061.
  • The proteins cdc2 and cyclin form a heterodimer (maturation promoting factor) that controls the major events of the cell cycle. A mathematical model for the interactions of cdc2 and cyclin is constructed. Simulation and analysis of the model show that the control system can operate in three modes: as a steady state with high maturation promoting factor activity, as a spontaneous oscillator, or as an excitable switch. We associate the steady state with metaphase arrest in unfertilized eggs, the spontaneous oscillations with rapid division cycles in early embryos, and the excitable switch with growth-controlled division cycles typical of nonembryonic cells.
Contributors
Submitter of the first revision: Nicolas Le Novère
Submitter of this revision: Nicolas Le Novère
Modellers: Nicolas Le Novère

Metadata information

isDescribedBy
PubMed 1831270
hasTaxon
Taxonomy Opisthokonta
isVersionOf
Gene Ontology mitotic cell cycle
hasVersion

Curation status
Curated

Tags

Connected external resources

SBGN view in Newt Editor

Name Description Size Actions

Model files

BIOMD0000000006_url.xml SBML L2V4 representation of Tyson1991 - Cell Cycle 2 var 15.10 KB Preview | Download

Additional files

BIOMD0000000006-biopax2.owl Auto-generated BioPAX (Level 2) 8.27 KB Preview | Download
BIOMD0000000006-biopax3.owl Auto-generated BioPAX (Level 3) 11.00 KB Preview | Download
BIOMD0000000006.m Auto-generated Octave file 3.45 KB Preview | Download
BIOMD0000000006.pdf Auto-generated PDF file 153.66 KB Preview | Download
BIOMD0000000006.png Auto-generated Reaction graph (PNG) 4.90 KB Preview | Download
BIOMD0000000006.sci Auto-generated Scilab file 1.04 KB Preview | Download
BIOMD0000000006.svg Auto-generated Reaction graph (SVG) 9.38 KB Preview | Download
BIOMD0000000006.vcml Auto-generated VCML file 900.00 Bytes Preview | Download
BIOMD0000000006.xpp Auto-generated XPP file 1.65 KB Preview | Download
BIOMD0000000006_manual.png Manually generated Reaction graph (PNG) 4.90 KB Preview | Download
BIOMD0000000006_manual.svg Manually generated Reaction graph (SVG) 9.38 KB Preview | Download
BIOMD0000000006_urn.xml Auto-generated SBML file with URNs 15.14 KB Preview | Download

  • Model originally submitted by : Nicolas Le Novère
  • Submitted: Sep 13, 2005 1:32:12 PM
  • Last Modified: May 16, 2013 3:38:56 PM
Revisions
  • Version: 2 public model Download this version
    • Submitted on: May 16, 2013 3:38:56 PM
    • Submitted by: Nicolas Le Novère
    • With comment: Current version of Tyson1991 - Cell Cycle 2 var
  • Version: 1 public model Download this version
    • Submitted on: Sep 13, 2005 1:32:12 PM
    • Submitted by: Nicolas Le Novère
    • With comment: Original import of Tyson1991_CellCycle_2var

(*) You might be seeing discontinuous revisions as only public revisions are displayed here. Any private revisions unpublished model revision of this model will only be shown to the submitter and their collaborators.

Legends
: Variable used inside SBML models


Species
Reactions
Reactions Rate Parameters
z = v-u [] []
u => EmptySet k6*u k6 = 1.0
z => u k4*z*(k4prime/k4+u^2) k4 = 180.0; k4prime = 0.018
v = kappa-k6*u kappa-k6*u kappa = 0.015; k6 = 1.0
EmptySet => z kappa kappa = 0.015
u = k4*(v-u)*(alpha+u^2)-k6*u k4*(v-u)*(alpha+u^2)-k6*u k4 = 180.0; alpha = NaN; k6 = 1.0
Curator's comment:
(added: 08 Apr 2011, 02:56:21, updated: 08 Apr 2011, 02:56:21)
Nullclines (left) and time course trajectories (right) of the 2 variable system for various values of kappa in the u-v plane. The values of kappa were chosen to give similar results to figure 4 in the original publication, but as the exact values used are not known, they only show approximately the same results. The right graph shows time courses trajectories (200 time units) for the various values of kappa, and the u-nullcline (orange). the high and low values of kappa (kappa = 0.13, blue and kappa = 0.0075, green)lead to stable steady states, the intermediate one (kappa=0.018, yellow) approaches a limit cycle. The results were calculated and plotted using XPPaut and the Biomodels.net sbml2xpp converter.