// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

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

// SBML reactions
// v1: S => x1 
// v2: x1 => x2 
// v3: x2 => x3 
// v4: x3 => P 
// vtrscsyn: mRsource => mR 
// vtrscdeg: mR => mRsink 
// vtrnlsyn: esource => e 
// vtrnldeg: e => esink 

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

  // variable definitions
  real compartment = 1 L;
  real k1 = 100;
  real K1s = .05;
  real L = .8;
  real cs = 1E-4;
  real cp = 1E-4;
  real K1x3 = 1;
  real K1x1 = 1E3;
  real n = 3;
  real Keq1 = 100;
  real k2 = 50;
  real K2x1 = 1;
  real K2x2 = 1;
  real Keq2 = 10;
  real k3 = 75;
  real K3x2 = 1;
  real K3x3 = 10;
  real Keq3 = 1;
  real V4 = 100;
  real K4x3 = .4;
  real K4P = 1;
  real Keq4 = 10;
  real Ktrx3 = .5;
  real ntr = 2;
  real ktrscdeg = .5;
  real ktrnlsyn = 1;
  real ktrnldeg = .5;
  real S = 1 M;
  real x1(time) M;
  real e(time) M;
  real x3(time) M;
  real x2(time) M;
  real P = 1 M;
  real mRsource = 1 M;
  real mR(time) M;
  real mRsink = 1 M;
  real esource = 1 M;
  real esink = 1 M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real vtrscsyn(time) katal;
  real vtrscdeg(time) katal;
  real vtrnlsyn(time) katal;
  real vtrnldeg(time) katal;

  // equations
  when (time=time.min) x1 = 1;
  x1:time = (v1 + -1*v2)/compartment;
  when (time=time.min) e = 1;
  e:time = (vtrnlsyn + -1*vtrnldeg)/compartment;
  when (time=time.min) x3 = 1;
  x3:time = (v3 + -1*v4)/compartment;
  when (time=time.min) x2 = 1;
  x2:time = (v2 + -1*v3)/compartment;
  when (time=time.min) mR = 1;
  mR:time = (vtrscsyn + -1*vtrscdeg)/compartment;
  v1 = e*k1*S*(1-x1/(Keq1*S))*(((1+S/K1s+x1/K1x1)/(1+cs*S/K1s+cp*x1/K1x1))^n+cs*L*(1+S/K1s+x1/K1x1)*(1+x3/K1x3)^n/(1+cs*S/K1s+cp*x1/K1x1))/(K1s*(1+S/K1s+x1/K1x1)*(((1+S/K1s+x1/K1x1)/(1+cs*S/K1s+cp*x1/K1x1))^n+L*(1+x3/K1x3)^n));
  v2 = e*k2*x1*(1-x2/(Keq2*x1))/(K2x1*(1+x1/K2x1+x2/K2x2));
  v3 = e*k3*x2*(1-x3/(Keq3*x2))/(K3x2*(1+x2/K3x2+x3/K3x3));
  v4 = V4*(1-P/(Keq4*x3))*x3/(K4x3*(1+P/K4P+x3/K4x3));
  vtrscsyn = (1+(x3/Ktrx3)^ntr)^(-1);
  vtrscdeg = ktrscdeg*mR;
  vtrnlsyn = ktrnlsyn*mR;
  vtrnldeg = e*ktrnldeg;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  K1s.sbmlRole="parameter";
  L.sbmlRole="parameter";
  cs.sbmlRole="parameter";
  cp.sbmlRole="parameter";
  K1x3.sbmlRole="parameter";
  K1x1.sbmlRole="parameter";
  n.sbmlRole="parameter";
  Keq1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  K2x1.sbmlRole="parameter";
  K2x2.sbmlRole="parameter";
  Keq2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  K3x2.sbmlRole="parameter";
  K3x3.sbmlRole="parameter";
  Keq3.sbmlRole="parameter";
  V4.sbmlRole="parameter";
  K4x3.sbmlRole="parameter";
  K4P.sbmlRole="parameter";
  Keq4.sbmlRole="parameter";
  Ktrx3.sbmlRole="parameter";
  ntr.sbmlRole="parameter";
  ktrscdeg.sbmlRole="parameter";
  ktrnlsyn.sbmlRole="parameter";
  ktrnldeg.sbmlRole="parameter";
  S.sbmlRole="species";
  S.sbmlCompartment="compartment";
  x1.sbmlRole="species";
  x1.sbmlCompartment="compartment";
  e.sbmlRole="species";
  e.sbmlCompartment="compartment";
  x3.sbmlRole="species";
  x3.sbmlCompartment="compartment";
  x2.sbmlRole="species";
  x2.sbmlCompartment="compartment";
  P.sbmlRole="species";
  P.sbmlCompartment="compartment";
  mRsource.sbmlRole="species";
  mRsource.sbmlCompartment="compartment";
  mR.sbmlRole="species";
  mR.sbmlCompartment="compartment";
  mRsink.sbmlRole="species";
  mRsink.sbmlCompartment="compartment";
  esource.sbmlRole="species";
  esource.sbmlCompartment="compartment";
  esink.sbmlRole="species";
  esink.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  vtrscsyn.sbmlRole="rate";
  vtrscdeg.sbmlRole="rate";
  vtrnlsyn.sbmlRole="rate";
  vtrnldeg.sbmlRole="rate";
}