#
# This file is automatically generated with 
# the System Biology Format Converter (http://sbfc.sourceforge.net/)
# from an SBML file.
#

#
# Model name = Vilar2006_TGFbeta
#
# is urn:miriam:biomodels.db:BIOMD0000000101
# is urn:miriam:biomodels.db:MODEL4023382414
# isDescribedBy urn:miriam:pubmed:16446785
#

# some function definitions that are allowed in SBML but not valid in xpp
ceil(x)=flr(1+x)

@delay=50


# Compartment: id = PM, name = Plasma membrane, constant
par PM=1.0

# Compartment: id = Endosome, name = Endosome, constant
par Endosome=1.0

# Parameter:   id =  ka, name = ka, constant
par ka=1.0

# Parameter:   id =  ligand, name = ligand
par ligand=3.0E-5

# Parameter:   id =  kcd, name = kcd, constant
par kcd=0.0277777778

# Parameter:   id =  klid, name = klid, constant
par klid=0.25

# Parameter:   id =  ki, name = ki, constant
par ki=0.3333333333333

# Parameter:   id =  pRI, name = pRI, constant
par pRI=8.0

# Parameter:   id =  kr, name = kr, constant
par kr=0.0333333333333333

# Parameter:   id =  alpha, name = alpha, constant
par alpha=1.0

# Parameter:   id =  pRII, name = pRII, constant
par pRII=4.0

# Reaction: id = v1, name = Ligand receptor complex formation

v1=ka*ligand*RI*RII

# Reaction: id = v2, name = Ligand receptor complex constitutive degradation

v2=kcd*lRIRII

# Reaction: id = v3, name = Ligand independent complex degradation

v3=klid*lRIRII

# Reaction: id = v4, name = Ligand receptor complex internalization

v4=ki*lRIRII

# Reaction: id = v5, name = RI synthesis

v5=pRI

# Reaction: id = v6, name = RI constitutive degradation

v6=kcd*RI

# Reaction: id = v7, name = RI internalization

v7=ki*RI

# Reaction: id = v8, name = RI recycling

v8=kr*RI_endo

# Reaction: id = v9, name = Ligand Receptor complex recycling

v9=kr*lRIRII_en

# Reaction: id = v10, name = RII synthesis

v10=pRII

# Reaction: id = v11, name = RII constitutive degradation

v11=kcd*RII

# Reaction: id = v12, name = RII internalization

v12=ki*RII

# Reaction: id = v13, name = RII recycling

v13=kr*RII_endo

# Species:   id = RI, name = Receptor 1, affected by kineticLaw

init RI=20.0
dRI/dt=(-1.0 * v1) + ( 1.0 * v5) + (-1.0 * v6) + (-1.0 * v7) + ( 1.0 * v8) + ( 1.0 * v9)

# Species:   id = RII, name = Receptor 2, affected by kineticLaw

init RII=20.0
dRII/dt=(-1.0 * v1) + ( 1.0 * v9) + ( 1.0 * v10) + (-1.0 * v11) + (-1.0 * v12) + ( 1.0 * v13)

# Species:   id = lRIRII, name = ligand receptor complex-plasma membrane, affected by kineticLaw

init lRIRII=0.0
dlRIRII/dt=( 1.0 * v1) + (-1.0 * v2) + (-1.0 * v3) + (-1.0 * v4)

# Species:   id = lRIRII_endo, name = ligand receptor complex-endosome, affected by kineticLaw

par lRIRII_en=40.0
aux lRIRII_en=lRIRII_en
dlRIRII_en/dt=( 1.0 * v4) + (-1.0 * v9)

# Species:   id = RI_endo, name = Receptor 1-endosome, affected by kineticLaw

init RI_endo=0.0
dRI_endo/dt=( 1.0 * v7) + (-1.0 * v8)

# Species:   id = RII_endo, name = Receptor 2 endosome, affected by kineticLaw

init RII_endo=0.0
dRII_endo/dt=( 1.0 * v12) + (-1.0 * v13)

# event : event_0000001
event_000=if (t >= 2500) then (1.5) else (0.2)
global 1 {event_000 - 1.1} {ligand=0.01}

@ meth=cvode, tol=1e-6, atol=1e-8
# @ maxstor=1e6
@ bound=40000, total=200
done