// 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
// R01975: space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA => Acetoacetyl_minus_CoA 
// R03026: Crotonoyl_minus_CoA <=> space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA 
// R04313: Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine <=> N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine 
// R00238: Acetyl_minus_CoA <=> Acetoacetyl_minus_CoA 

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 E3_space_4_space_minus_space_minus_ = 0 M;
  real E4_space_2_space_1_space_17 = 0 M;
  real E1_space_1_space_1_space_35 = 0 M;
  real E2_space_3_space_1_space_9 = 0 M;
  real Crotonoyl_minus_CoA(time) M;
  real space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA(time) M;
  real Acetoacetyl_minus_CoA(time) M;
  real Acetyl_minus_CoA(time) M;
  real Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine(time) M;
  real N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine(time) M;
  real R01975(time) katal;
  extern real R01975.kf;
  extern real R01975.kb;
  real R03026(time) katal;
  extern real R03026.kf;
  extern real R03026.kb;
  real R04313(time) katal;
  extern real R04313.kf;
  extern real R04313.kb;
  real R00238(time) katal;
  extern real R00238.kf;
  extern real R00238.kb;

  // equations
  when (time=time.min) Crotonoyl_minus_CoA = 0;
  Crotonoyl_minus_CoA:time = (-1*R03026)/uVol;
  when (time=time.min) space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA = 0;
  space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA:time = (-1*R01975 + R03026)/uVol;
  when (time=time.min) Acetoacetyl_minus_CoA = 0;
  Acetoacetyl_minus_CoA:time = (R01975 + R00238)/uVol;
  when (time=time.min) Acetyl_minus_CoA = 0;
  Acetyl_minus_CoA:time = (-1*R00238)/uVol;
  when (time=time.min) Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine = 0;
  Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine:time = (-1*R04313)/uVol;
  when (time=time.min) N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine = 0;
  N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine:time = (R04313)/uVol;
  R01975 = R01975.kf*space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA - R01975.kb*Acetoacetyl_minus_CoA;
  R03026 = R03026.kf*Crotonoyl_minus_CoA - R03026.kb*space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA;
  R04313 = R04313.kf*Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine - R04313.kb*N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine;
  R00238 = R00238.kf*Acetyl_minus_CoA - R00238.kb*Acetoacetyl_minus_CoA;

  // variable properties
  default.sbmlRole="compartment";
  uVol.sbmlRole="compartment";
  E3_space_4_space_minus_space_minus_.sbmlRole="species";
  E3_space_4_space_minus_space_minus_.sbmlCompartment="uVol";
  E4_space_2_space_1_space_17.sbmlRole="species";
  E4_space_2_space_1_space_17.sbmlCompartment="uVol";
  E1_space_1_space_1_space_35.sbmlRole="species";
  E1_space_1_space_1_space_35.sbmlCompartment="uVol";
  E2_space_3_space_1_space_9.sbmlRole="species";
  E2_space_3_space_1_space_9.sbmlCompartment="uVol";
  Crotonoyl_minus_CoA.sbmlRole="species";
  Crotonoyl_minus_CoA.sbmlCompartment="uVol";
  space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA.sbmlRole="species";
  space_S_space_minus_3_minus_Hydroxybutanoyl_minus_CoA.sbmlCompartment="uVol";
  Acetoacetyl_minus_CoA.sbmlRole="species";
  Acetoacetyl_minus_CoA.sbmlCompartment="uVol";
  Acetyl_minus_CoA.sbmlRole="species";
  Acetyl_minus_CoA.sbmlCompartment="uVol";
  Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine.sbmlRole="species";
  Protein_space_N6_space_N6_space_N6_minus_trimethyl_minus_L_minus_lysine.sbmlCompartment="uVol";
  N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine.sbmlRole="species";
  N6_space_N6_space_N6_minus_Trimethyl_minus_L_minus_lysine.sbmlCompartment="uVol";
  R01975.sbmlRole="rate";
  R01975.kf.sbmlRole="rate";
  R01975.kb.sbmlRole="rate";
  R03026.sbmlRole="rate";
  R03026.kf.sbmlRole="rate";
  R03026.kb.sbmlRole="rate";
  R04313.sbmlRole="rate";
  R04313.kf.sbmlRole="rate";
  R04313.kb.sbmlRole="rate";
  R00238.sbmlRole="rate";
  R00238.kf.sbmlRole="rate";
  R00238.kb.sbmlRole="rate";
}