// This model generated automatically from SBML
// WARNING: SBML events are not currently supported: Start
// WARNING: SBML events are not currently supported: Division

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// G2K_Creation: <=> G2K 
// G1K_Creation: <=> G1K 
// Cdc2Phos: G2K => PG2 
// G2R_Creation: G2K R => G2R 
// PG2R_Creation: PG2 R => PG2R 
// Rum1DegInG2R: G2R <=> G2K 
// Rum1Deg: R 
// Rum1DegInPG2R: PG2R <=> PG2 
// RumDegInG1R: G1R <=> G1K 
// G2K_dissoc: G2K 
// PG2_dissoc: PG2 
// G1K_Dissociation: G1K 
// PG2R_Dissociation: PG2R <=> R 
// G2R_Dissociation: G2R <=> R 
// G1R_Dissociation: G1R <=> R 
// G1R_Binding: G1K R => G1R 
// G2R_Dissociation_UbE: G2R <=> R 
// PG2R_Dissociation_UbE: PG2R <=> R 
// Rum1_Production: <=> R 
// Rum1_Deg_SPF: R 
// IE_Reaction: IEB => IE 
// UbE_Reaction: UbEB => UbE 
// UbE2_Reaction: UbE2B => UbE2 
// Wee1_Reaction: Wee1B => Wee1 
// Cdc25_Reaction: Cdc25B => Cdc25 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real Cell = 1 L;
  real mu = .00495;
  real Mass(time);
  real k1 = .015;
  real k2prime = .05;
  real k3 = .09375;
  real k4 = .1875;
  real k5 = .00175;
  real k6prime = 0;
  real k7 = 100;
  real k7r = .1;
  real k8 = 10;
  real k8r = .1;
  real kc = 1;
  real kcr = .25;
  real ki = .4;
  real kir = .1;
  real kp = 3.25;
  real ku = .2;
  real kur = .1;
  real ku2 = 1;
  real kur2 = .3;
  real kw = 1;
  real kwr = .25;
  real V2 = .25;
  real V2prime = .0075;
  real V6 = 7.5;
  real V6prime = .0375;
  real V25 = .5;
  real V25prime = .025;
  real Vw = .35;
  real Vwprime = .035;
  real Kmc = .1;
  real Kmcr = .1;
  real Kmi = .01;
  real Kmir = .01;
  real Kmp = .001;
  real Kmu = .01;
  real Kmur = .01;
  real Kmu2 = .05;
  real Kmur2 = .05;
  real Kmw = .1;
  real Kmwr = .1;
  real alpha = .25;
  real beta = .05;
  real Cig1 = 0;
  real k2(time);
  real k6(time);
  real kwee(time);
  real k25(time);
  real UbE(time) mol;
  real UbE2(time) mol;
  real Wee1(time) mol;
  real Cdc25(time) mol;
  real G2K(time) mol;
  real R(time) mol;
  real G1K(time) mol;
  real IE(time) mol;
  real PG2(time) mol;
  real G1R(time) mol;
  real G2R(time) mol;
  real PG2R(time) mol;
  real SPF(time) mol;
  real MPF(time) mol;
  real IEB(time) mol;
  real UbEB(time) mol;
  real UbE2B(time) mol;
  real Wee1B(time) mol;
  real Cdc25B(time) mol;
  real Rum1Total(time) mol;
  real Cdc13Total(time) mol;
  real Cig2Total(time) mol;
  real G2K_Creation(time) mole/min;
  real G1K_Creation(time) mole/min;
  real Cdc2Phos(time) mole/min;
  real G2R_Creation(time) mole/min;
  real PG2R_Creation(time) mole/min;
  real Rum1DegInG2R(time) mole/min;
  real Rum1Deg(time) mole/min;
  real Rum1DegInPG2R(time) mole/min;
  real RumDegInG1R(time) mole/min;
  real G2K_dissoc(time) mole/min;
  real PG2_dissoc(time) mole/min;
  real G1K_Dissociation(time) mole/min;
  real PG2R_Dissociation(time) mole/min;
  real G2R_Dissociation(time) mole/min;
  real G1R_Dissociation(time) mole/min;
  real G1R_Binding(time) mole/min;
  real G2R_Dissociation_UbE(time) mole/min;
  real PG2R_Dissociation_UbE(time) mole/min;
  real Rum1_Production(time) mole/min;
  real Rum1_Deg_SPF(time) mole/min;
  real IE_Reaction(time) mole/min;
  real UbE_Reaction(time) mole/min;
  real UbE2_Reaction(time) mole/min;
  real Wee1_Reaction(time) mole/min;
  real Cdc25_Reaction(time) mole/min;

  // equations
  when (time=time.min) Mass = .49;
  Mass:time = Mass*mu;
  k2 = UbE*V2+(1-UbE)*V2prime;
  k6 = UbE2*V6+(1-UbE2)*V6prime;
  kwee = Vwprime*(1-Wee1)+Vw*Wee1;
  k25 = Cdc25*V25+(1-Cdc25)*V25prime;
  when (time=time.min) UbE = .11;
  UbE:time = UbE_Reaction;
  when (time=time.min) UbE2 = 0;
  UbE2:time = UbE2_Reaction;
  when (time=time.min) Wee1 = 0;
  Wee1:time = Wee1_Reaction;
  when (time=time.min) Cdc25 = 0;
  Cdc25:time = Cdc25_Reaction;
  when (time=time.min) G2K = 0;
  G2K:time = G2K_Creation + -1*Cdc2Phos + -1*G2R_Creation + Rum1DegInG2R + -1*G2K_dissoc;
  when (time=time.min) R = .4;
  R:time = -1*G2R_Creation + -1*PG2R_Creation + -1*Rum1Deg + PG2R_Dissociation + G2R_Dissociation + G1R_Dissociation + -1*G1R_Binding + G2R_Dissociation_UbE + PG2R_Dissociation_UbE + Rum1_Production + -1*Rum1_Deg_SPF;
  when (time=time.min) G1K = 0;
  G1K:time = G1K_Creation + RumDegInG1R + -1*G1K_Dissociation + -1*G1R_Binding;
  when (time=time.min) IE = 0;
  IE:time = IE_Reaction;
  when (time=time.min) PG2 = 0;
  PG2:time = Cdc2Phos + -1*PG2R_Creation + Rum1DegInPG2R + -1*PG2_dissoc;
  when (time=time.min) G1R = 0;
  G1R:time = -1*RumDegInG1R + -1*G1R_Dissociation + G1R_Binding;
  when (time=time.min) G2R = 0;
  G2R:time = G2R_Creation + -1*Rum1DegInG2R + -1*G2R_Dissociation + -1*G2R_Dissociation_UbE;
  when (time=time.min) PG2R = 0;
  PG2R:time = PG2R_Creation + -1*Rum1DegInPG2R + -1*PG2R_Dissociation + -1*PG2R_Dissociation_UbE;
  SPF = Cig1+alpha*G1K+MPF;
  MPF = G2K+beta*PG2;
  IEB = 1-IE;
  UbEB = 1-UbE;
  UbE2B = 1-UbE2;
  Wee1B = 1-Wee1;
  Cdc25B = 1-Cdc25;
  Rum1Total = G1R+G2R+PG2R+R;
  Cdc13Total = G2K+G2R+PG2+PG2R;
  Cig2Total = G1K+G1R;
  G2K_Creation = k1;
  G1K_Creation = k5;
  Cdc2Phos = G2K*kwee-k25*PG2;
  G2R_Creation = G2K*k7*R-G2R*k7r;
  PG2R_Creation = k7*PG2*R-k7r*PG2R;
  Rum1DegInG2R = G2R*k4;
  Rum1Deg = k4*R;
  Rum1DegInPG2R = k4*PG2R;
  RumDegInG1R = G1R*k4;
  G2K_dissoc = G2K*k2;
  PG2_dissoc = k2*PG2;
  G1K_Dissociation = G1K*k6;
  PG2R_Dissociation = k2prime*PG2R;
  G2R_Dissociation = G2R*k2prime;
  G1R_Dissociation = G1R*k6prime;
  G1R_Binding = G1K*k8*R-G1R*k8r;
  G2R_Dissociation_UbE = G2R*k2;
  PG2R_Dissociation_UbE = k2*PG2R;
  Rum1_Production = k3;
  Rum1_Deg_SPF = kp*Mass*R*SPF/(Kmp+R);
  IE_Reaction = IEB*ki*MPF/(IEB+Kmi)-IE*kir/(IE+Kmir);
  UbE_Reaction = IE*ku*UbEB/(Kmu+UbEB)-kur*UbE/(Kmur+UbE);
  UbE2_Reaction = ku2*MPF*UbE2B/(Kmu2+UbE2B)-kur2*UbE2/(Kmur2+UbE2);
  Wee1_Reaction = kwr*Wee1B/(Kmwr+Wee1B)-kw*MPF*Wee1/(Kmw+Wee1);
  Cdc25_Reaction = Cdc25B*kc*MPF/(Cdc25B+Kmc)-Cdc25*kcr/(Cdc25+Kmcr);

  // variable properties
  Cell.sbmlRole="compartment";
  mu.sbmlRole="parameter";
  Mass.sbmlRole="parameter";
  k1.sbmlRole="parameter";
  k2prime.sbmlRole="parameter";
  k2prime.sbmlName="k2'";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6prime.sbmlRole="parameter";
  k6prime.sbmlName="k6'";
  k7.sbmlRole="parameter";
  k7r.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k8r.sbmlRole="parameter";
  kc.sbmlRole="parameter";
  kcr.sbmlRole="parameter";
  ki.sbmlRole="parameter";
  kir.sbmlRole="parameter";
  kp.sbmlRole="parameter";
  ku.sbmlRole="parameter";
  kur.sbmlRole="parameter";
  ku2.sbmlRole="parameter";
  kur2.sbmlRole="parameter";
  kw.sbmlRole="parameter";
  kwr.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  V2prime.sbmlRole="parameter";
  V2prime.sbmlName="V2'";
  V6.sbmlRole="parameter";
  V6prime.sbmlRole="parameter";
  V6prime.sbmlName="V6'";
  V25.sbmlRole="parameter";
  V25prime.sbmlRole="parameter";
  V25prime.sbmlName="V25'";
  Vw.sbmlRole="parameter";
  Vwprime.sbmlRole="parameter";
  Vwprime.sbmlName="Vw'";
  Kmc.sbmlRole="parameter";
  Kmcr.sbmlRole="parameter";
  Kmi.sbmlRole="parameter";
  Kmir.sbmlRole="parameter";
  Kmp.sbmlRole="parameter";
  Kmu.sbmlRole="parameter";
  Kmur.sbmlRole="parameter";
  Kmu2.sbmlRole="parameter";
  Kmur2.sbmlRole="parameter";
  Kmw.sbmlRole="parameter";
  Kmwr.sbmlRole="parameter";
  alpha.sbmlRole="parameter";
  beta.sbmlRole="parameter";
  Cig1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  kwee.sbmlRole="parameter";
  k25.sbmlRole="parameter";
  UbE.sbmlRole="species";
  UbE.sbmlCompartment="Cell";
  UbE2.sbmlRole="species";
  UbE2.sbmlCompartment="Cell";
  Wee1.sbmlRole="species";
  Wee1.sbmlCompartment="Cell";
  Cdc25.sbmlRole="species";
  Cdc25.sbmlCompartment="Cell";
  G2K.sbmlRole="species";
  G2K.sbmlCompartment="Cell";
  R.sbmlRole="species";
  R.sbmlCompartment="Cell";
  G1K.sbmlRole="species";
  G1K.sbmlCompartment="Cell";
  IE.sbmlRole="species";
  IE.sbmlCompartment="Cell";
  PG2.sbmlRole="species";
  PG2.sbmlCompartment="Cell";
  G1R.sbmlRole="species";
  G1R.sbmlCompartment="Cell";
  G2R.sbmlRole="species";
  G2R.sbmlCompartment="Cell";
  PG2R.sbmlRole="species";
  PG2R.sbmlCompartment="Cell";
  SPF.sbmlRole="species";
  SPF.sbmlCompartment="Cell";
  MPF.sbmlRole="species";
  MPF.sbmlCompartment="Cell";
  IEB.sbmlRole="species";
  IEB.sbmlCompartment="Cell";
  UbEB.sbmlRole="species";
  UbEB.sbmlCompartment="Cell";
  UbE2B.sbmlRole="species";
  UbE2B.sbmlCompartment="Cell";
  Wee1B.sbmlRole="species";
  Wee1B.sbmlCompartment="Cell";
  Cdc25B.sbmlRole="species";
  Cdc25B.sbmlCompartment="Cell";
  Rum1Total.sbmlRole="species";
  Rum1Total.sbmlCompartment="Cell";
  Cdc13Total.sbmlRole="species";
  Cdc13Total.sbmlCompartment="Cell";
  Cig2Total.sbmlRole="species";
  Cig2Total.sbmlCompartment="Cell";
  G2K_Creation.sbmlRole="rate";
  G1K_Creation.sbmlRole="rate";
  Cdc2Phos.sbmlRole="rate";
  G2R_Creation.sbmlRole="rate";
  PG2R_Creation.sbmlRole="rate";
  Rum1DegInG2R.sbmlRole="rate";
  Rum1Deg.sbmlRole="rate";
  Rum1DegInPG2R.sbmlRole="rate";
  RumDegInG1R.sbmlRole="rate";
  G2K_dissoc.sbmlRole="rate";
  PG2_dissoc.sbmlRole="rate";
  G1K_Dissociation.sbmlRole="rate";
  PG2R_Dissociation.sbmlRole="rate";
  G2R_Dissociation.sbmlRole="rate";
  G1R_Dissociation.sbmlRole="rate";
  G1R_Binding.sbmlRole="rate";
  G2R_Dissociation_UbE.sbmlRole="rate";
  PG2R_Dissociation_UbE.sbmlRole="rate";
  Rum1_Production.sbmlRole="rate";
  Rum1_Deg_SPF.sbmlRole="rate";
  IE_Reaction.sbmlRole="rate";
  UbE_Reaction.sbmlRole="rate";
  UbE2_Reaction.sbmlRole="rate";
  Wee1_Reaction.sbmlRole="rate";
  Cdc25_Reaction.sbmlRole="rate";
}