#ferreira1



#reactions
v1a:
 {1.0}Lysine + {1.0}Glucose = {1.0}Schiff
  compartment * p1 * k1a * Glucose * Lysine
v1b:
 {1.0}Schiff = {1.0}Lysine + {1.0}Glucose
  compartment * k1b * Schiff
v2a:
 {1.0}Schiff = {1.0}Amadori
  compartment * p2 * k2a * Schiff
v2b:
 {1.0}Amadori = {1.0}Schiff
  compartment * p2 * k2b * Amadori
v3:
 {1.0}Glucose = {1.0}Glyoxal
  compartment * ox * p3 * k3 * pow(Glucose / 0.25, 0.36)
v4:
 {1.0}Amadori = {1.0}CML
  compartment * ox * p4 * k4 * Amadori
v5:
 {1.0}Lysine + {1.0}Glyoxal = {1.0}CML
  compartment * ox * p5 * k5 * Glyoxal * Lysine
v5b:
 {1.0}Glyoxal = $pool
  compartment * k5b * Glyoxal
v6:
 {1.0}Schiff = {1.0}CML
  compartment * ox * p6 * k3 * pow(Schiff / 0.25, 0.36)
v7a:
 $pool = {1.0}Lysine
  compartment * 0.05 * ox * p7 * k3 * pow(Schiff / 0.25, 0.36)
v7b:
 $pool = {1.0}Glyoxal
  compartment * 0.005 * ox * p7 * k3 * pow(Schiff / 0.25, 0.36)
v7c:
 {1.0}Schiff = $pool
  compartment * ox * p7 * k3 * pow(Schiff / 0.25, 0.36)


#parameters
v1a_p1 = 0.115
v1a_k1a = 0.09
v1b_k1b = 0.36
v2a_p2 = 0.75
v2a_k2a = 0.033
v2b_p2 = 0.75
v2b_k2b = 0.0012
v3_ox = 1.0
v3_p3 = 1.0
v3_k3 = 7.92e-07
v4_ox = 1.0
v4_p4 = 1.0
v4_k4 = 8.6e-05
v5_ox = 1.0
v5_p5 = 1.0
v5_k5 = 0.019
v5b_k5b = 0.0017
v6_ox = 1.0
v6_p6 = 2.7
v6_k3 = 7.92e-07
v7a_ox = 1.0
v7a_p7 = 60.0
v7a_k3 = 7.92e-07
v7b_ox = 1.0
v7b_p7 = 60.0
v7b_k3 = 7.92e-07
v7c_ox = 1.0
v7c_p7 = 60.0
v7c_k3 = 7.92e-07
compartment = 1.0


#initial values
Amadori = 0.0
CML = 0.0
Glucose = 0.25
Glyoxal = 0.0
Lysine = 0.0034
Schiff = 0.0


#assignment rules
!F Fraction_Amadori = Amadori / 0.0034
!F Fraction_CML = CML / 0.0034