// 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
// R03245: space_R_space_minus_5_minus_Phosphomevalonate <=> space_R_space_minus_5_minus_Diphosphomevalonate 
// R02082: space_R_space_minus_Mevalonate <=> space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA 
// R01121: space_R_space_minus_5_minus_Diphosphomevalonate <=> Isopentenyl_space_diphosphate 
// R01123: Isopentenyl_space_diphosphate => Dimethylallyl_space_diphosphate 
// R02245: space_R_space_minus_Mevalonate <=> space_R_space_minus_5_minus_Phosphomevalonate 

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

  // variable definitions
  real default = 1 L;
  real uVol = 1 L;
  real E5_space_3_space_3_space_2 = 0 M;
  real E4_space_1_space_1_space_33 = 0 M;
  real E2_space_7_space_4_space_2 = 0 M;
  real E2_space_7_space_1_space_36 = 0 M;
  real E1_space_1_space_1_space_34 = 0 M;
  real Dimethylallyl_space_diphosphate(time) M;
  real Isopentenyl_space_diphosphate(time) M;
  real space_R_space_minus_5_minus_Diphosphomevalonate(time) M;
  real space_R_space_minus_5_minus_Phosphomevalonate(time) M;
  real space_R_space_minus_Mevalonate(time) M;
  real space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA(time) M;
  real R03245(time) katal;
  extern real R03245.kf;
  extern real R03245.kb;
  real R02082(time) katal;
  extern real R02082.kf;
  extern real R02082.kb;
  real R01121(time) katal;
  extern real R01121.kf;
  extern real R01121.kb;
  real R01123(time) katal;
  extern real R01123.kf;
  extern real R01123.kb;
  real R02245(time) katal;
  extern real R02245.kf;
  extern real R02245.kb;

  // equations
  when (time=time.min) Dimethylallyl_space_diphosphate = 0;
  Dimethylallyl_space_diphosphate:time = (R01123)/uVol;
  when (time=time.min) Isopentenyl_space_diphosphate = 0;
  Isopentenyl_space_diphosphate:time = (R01121 + -1*R01123)/uVol;
  when (time=time.min) space_R_space_minus_5_minus_Diphosphomevalonate = 0;
  space_R_space_minus_5_minus_Diphosphomevalonate:time = (R03245 + -1*R01121)/uVol;
  when (time=time.min) space_R_space_minus_5_minus_Phosphomevalonate = 0;
  space_R_space_minus_5_minus_Phosphomevalonate:time = (-1*R03245 + R02245)/uVol;
  when (time=time.min) space_R_space_minus_Mevalonate = 0;
  space_R_space_minus_Mevalonate:time = (-1*R02082 + -1*R02245)/uVol;
  when (time=time.min) space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA = 0;
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA:time = (R02082)/uVol;
  R03245 = R03245.kf*space_R_space_minus_5_minus_Phosphomevalonate - R03245.kb*space_R_space_minus_5_minus_Diphosphomevalonate;
  R02082 = R02082.kf*space_R_space_minus_Mevalonate - R02082.kb*space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA;
  R01121 = R01121.kf*space_R_space_minus_5_minus_Diphosphomevalonate - R01121.kb*Isopentenyl_space_diphosphate;
  R01123 = R01123.kf*Isopentenyl_space_diphosphate - R01123.kb*Dimethylallyl_space_diphosphate;
  R02245 = R02245.kf*space_R_space_minus_Mevalonate - R02245.kb*space_R_space_minus_5_minus_Phosphomevalonate;

  // variable properties
  default.sbmlRole="compartment";
  uVol.sbmlRole="compartment";
  E5_space_3_space_3_space_2.sbmlRole="species";
  E5_space_3_space_3_space_2.sbmlCompartment="uVol";
  E4_space_1_space_1_space_33.sbmlRole="species";
  E4_space_1_space_1_space_33.sbmlCompartment="uVol";
  E2_space_7_space_4_space_2.sbmlRole="species";
  E2_space_7_space_4_space_2.sbmlCompartment="uVol";
  E2_space_7_space_1_space_36.sbmlRole="species";
  E2_space_7_space_1_space_36.sbmlCompartment="uVol";
  E1_space_1_space_1_space_34.sbmlRole="species";
  E1_space_1_space_1_space_34.sbmlCompartment="uVol";
  Dimethylallyl_space_diphosphate.sbmlRole="species";
  Dimethylallyl_space_diphosphate.sbmlCompartment="uVol";
  Isopentenyl_space_diphosphate.sbmlRole="species";
  Isopentenyl_space_diphosphate.sbmlCompartment="uVol";
  space_R_space_minus_5_minus_Diphosphomevalonate.sbmlRole="species";
  space_R_space_minus_5_minus_Diphosphomevalonate.sbmlCompartment="uVol";
  space_R_space_minus_5_minus_Phosphomevalonate.sbmlRole="species";
  space_R_space_minus_5_minus_Phosphomevalonate.sbmlCompartment="uVol";
  space_R_space_minus_Mevalonate.sbmlRole="species";
  space_R_space_minus_Mevalonate.sbmlCompartment="uVol";
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA.sbmlRole="species";
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA.sbmlCompartment="uVol";
  R03245.sbmlRole="rate";
  R03245.kf.sbmlRole="rate";
  R03245.kb.sbmlRole="rate";
  R02082.sbmlRole="rate";
  R02082.kf.sbmlRole="rate";
  R02082.kb.sbmlRole="rate";
  R01121.sbmlRole="rate";
  R01121.kf.sbmlRole="rate";
  R01121.kb.sbmlRole="rate";
  R01123.sbmlRole="rate";
  R01123.kf.sbmlRole="rate";
  R01123.kb.sbmlRole="rate";
  R02245.sbmlRole="rate";
  R02245.kf.sbmlRole="rate";
  R02245.kb.sbmlRole="rate";
}