// 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: Hser => Phser 
// vCys: Phser Cys => Cystathionine Phi 
// vThr: Phser => Thr Phi 

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

  // variable definitions
  real compartment = 1 L;
  real V1 = 1;
  real kcat2 = 30;
  real KmCYS = 460;
  real KmPHSER = 2500;
  real Ki2 = 2E3;
  real Ki3 = 1E3;
  real CGS = .7;
  real TS = 5;
  real AdoMet = 20;
  real Hser = 0 M;
  real Phser(time) M;
  real Cys = 15 M;
  real Cystathionine = 0 M;
  real Phi = 1E4 M;
  real Thr = 0 M;
  real v1(time) katal;
  real vCys(time) katal;
  real vThr(time) katal;

  // equations
  when (time=time.min) Phser = 0;
  Phser:time = (v1 + -1*vCys + -1*vThr)/compartment;
  v1 = V1;
  vCys = CGS*(kcat2/(1+KmCYS/Cys))*Phser/(Phser+KmPHSER*(1+Phi/Ki2)/(1+KmCYS/Cys));
  vThr = TS*((.42+3.5*AdoMet^2/73)/(1+AdoMet^2/73))*Phser/((1+Phi/Ki3)*(250*(1+AdoMet/.5)/(1+AdoMet/1.09))/(1+AdoMet^2/142)+Phser);

  // variable properties
  compartment.sbmlRole="compartment";
  V1.sbmlRole="parameter";
  kcat2.sbmlRole="parameter";
  KmCYS.sbmlRole="parameter";
  KmPHSER.sbmlRole="parameter";
  Ki2.sbmlRole="parameter";
  Ki3.sbmlRole="parameter";
  CGS.sbmlRole="parameter";
  TS.sbmlRole="parameter";
  AdoMet.sbmlRole="parameter";
  Hser.sbmlRole="species";
  Hser.sbmlCompartment="compartment";
  Phser.sbmlRole="species";
  Phser.sbmlCompartment="compartment";
  Cys.sbmlRole="species";
  Cys.sbmlCompartment="compartment";
  Cystathionine.sbmlRole="species";
  Cystathionine.sbmlCompartment="compartment";
  Phi.sbmlRole="species";
  Phi.sbmlCompartment="compartment";
  Thr.sbmlRole="species";
  Thr.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  vCys.sbmlRole="rate";
  vThr.sbmlRole="rate";
}