# Kholod06BistSwitch_a.ode.ode
#
# Kholodenko (2006) Nat Rev Mol Cell Bio 7, p 165
# Box 2 Example (a)
#

Mp'=vkin(Mp,Mtot-Mp)-vphos(Mp,Mtot-Mp)

vkin(Mp,M)=((k_kincat*Ekin*M)/(Km1+M))*((1+A*(Mp/Ka))/(1+(Mp/Ka)))
vphos(Mp,M)=k_phoskin*Ephos*Mp/(Km2+Mp)

par k_kincat=1,A=100,Ka=500,Km1=500,Ekin=80
par k_phoskin=1,Km2=10,Ephos=200
par Mtot=300

@ total=20,xp=S,yp=C,dt=0.01,xlo=0,xhi=20,ylo=0,yhi=300,bounds=1000
done