#file arabiposis_branch.ode
#from Curien et al., 2003, Eur. J. Biochemistry 270, p. 4615-4627
#Problem 5.6.7

#declare parameters
par kcatCGS=30 
par KmCGScys=460
par KmCGSPhser=15000
 
par JPhser=0.3
 
par Cys=250
par CGS=0.7
par ADOMET=20
par TS=5

#declare dynamics

Phser'= Jphser - TS*(4.9*0.000001+ 5.6*0.0001*ADOMET^2.9/(32^2.9+ADOMET^2.9))*Phser - (kcatCGS/(1+KmCGScys/Cys))*CGS*Phser/(((KmCGSPhser/(1+KmCGScys/Cys)))+Phser)

#reaction fluxes:

aux vTS=TS*(4.9*0.000001+ 5.6*0.0001*ADOMET^2.9/(32^2.9+ADOMET^2.9))*Phser
aux vCGS=(kcatCGS/(1+KmCGScys/Cys))*CGS*Phser/(((KmCGSPhser/(1+KmCGScys/Cys)))+Phser)

#declare initial conditions
init Phser=0

# run the simulation in the GUI by selecting (I)nitialconds|(G)o

@ maxstor=1000000
@ total=2000
@ bound=1000

done