// This model generated automatically from SBML
// WARNING: Renamed variable source to source_ due to JSim namespace restrictions

// unit definitions
import nsrunit;
unit conversion off;

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

// SBML reactions
// v1: source => CYCD1cdk4 
// v2: CYCD1cdk4 => sink 
// v3: source => CYCEcdk2 
// v4: CYCEcdk2 => sink 
// v5: source => Rbhypo unphosRb 
// v6: Rbhypo freeE2F => RbE2F 
// v7: RbE2F unphosRb => freeE2F 
// v8: source => MPF 
// v9: source => MPF 
// v10: source => MPF 
// v11: MPF => sink 

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

  // variable definitions
  real compartment = 1 L;
  real aD = .1;
  real aE = .038;
  real aM = .02;
  real h = .01;
  real qD = .6;
  real qE = .6;
  real qM = .8;
  real f = .05;
  real g = .132;
  real pS = .15;
  real pD = .12;
  real pE = .024;
  real pM = .12;
  real dD = .1;
  real dE = .05;
  real dM = .26;
  real ET = 1.5;
  real alpha1 = .2;
  real alpha2 = 11;
  real beta1 = 2;
  real beta2 = 44;
  real RT = 2;
  real GF = 1;
  real source_ = 1 M;
  real CYCD1cdk4(time) M;
  real sink = 1 M;
  real CYCEcdk2(time) M;
  real Rbhypo(time) M;
  real unphosRb(time) M;
  real freeE2F(time) M;
  real RbE2F(time) M;
  real MPF(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;

  // equations
  when (time=time.min) CYCD1cdk4 = .5;
  CYCD1cdk4:time = (v1 + -1*v2)/compartment;
  when (time=time.min) CYCEcdk2 = .3;
  CYCEcdk2:time = (v3 + -1*v4)/compartment;
  when (time=time.min) Rbhypo = .2;
  Rbhypo:time = (v5 + -1*v6)/compartment;
  when (time=time.min) unphosRb = 1;
  unphosRb:time = (v5 + -1*v7)/compartment;
  when (time=time.min) freeE2F = .7;
  freeE2F:time = (-1*v6 + v7)/compartment;
  when (time=time.min) RbE2F = .8;
  RbE2F:time = (v6 + -1*v7)/compartment;
  when (time=time.min) MPF = .7;
  MPF:time = (v8 + v9 + v10 + -1*v11)/compartment;
  v1 = aD*(GF/(h+GF));
  v2 = dD*CYCEcdk2*CYCD1cdk4;
  v3 = aE*(1+alpha1*GF/(beta1+GF))*(1+alpha2*RT/(beta2+RT))*(ET-RbE2F);
  v4 = dE*MPF*CYCEcdk2;
  v5 = pM*(RT-RbE2F-Rbhypo)*MPF/(qM+(RT-RbE2F-Rbhypo)+MPF);
  v6 = pS*(ET-RbE2F)*Rbhypo;
  v7 = pD*RbE2F*CYCD1cdk4/(qD+RbE2F+CYCD1cdk4)+pE*RbE2F*CYCEcdk2/(qE+RbE2F+CYCEcdk2);
  v8 = aM*CYCEcdk2;
  v9 = f*(ET-RbE2F);
  v10 = g*MPF*MPF*CYCEcdk2;
  v11 = dM*MPF;

  // variable properties
  compartment.sbmlRole="compartment";
  aD.sbmlRole="parameter";
  aE.sbmlRole="parameter";
  aM.sbmlRole="parameter";
  h.sbmlRole="parameter";
  qD.sbmlRole="parameter";
  qE.sbmlRole="parameter";
  qM.sbmlRole="parameter";
  f.sbmlRole="parameter";
  g.sbmlRole="parameter";
  pS.sbmlRole="parameter";
  pD.sbmlRole="parameter";
  pE.sbmlRole="parameter";
  pM.sbmlRole="parameter";
  dD.sbmlRole="parameter";
  dE.sbmlRole="parameter";
  dM.sbmlRole="parameter";
  ET.sbmlRole="parameter";
  alpha1.sbmlRole="parameter";
  alpha2.sbmlRole="parameter";
  beta1.sbmlRole="parameter";
  beta2.sbmlRole="parameter";
  RT.sbmlRole="parameter";
  GF.sbmlRole="parameter";
  source_.sbmlRole="species";
  source_.sbmlCompartment="compartment";
  CYCD1cdk4.sbmlRole="species";
  CYCD1cdk4.sbmlCompartment="compartment";
  sink.sbmlRole="species";
  sink.sbmlCompartment="compartment";
  CYCEcdk2.sbmlRole="species";
  CYCEcdk2.sbmlCompartment="compartment";
  Rbhypo.sbmlRole="species";
  Rbhypo.sbmlCompartment="compartment";
  unphosRb.sbmlRole="species";
  unphosRb.sbmlCompartment="compartment";
  freeE2F.sbmlRole="species";
  freeE2F.sbmlCompartment="compartment";
  RbE2F.sbmlRole="species";
  RbE2F.sbmlCompartment="compartment";
  MPF.sbmlRole="species";
  MPF.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
}