TITLE Synapsin/actin/vesicle interactions

PARAMETER {

k1on = 1.0e4 (/M/sec) k1off = 1.0e-3 (/sec)

k2on = 1.0e3 (/M/sec) k2off = 1.0e-2 (/sec)

k3on = 1.0e3 (/M/sec) k3off = 1.0e-2 (/sec)

k4on = 1.0e2 (/M/sec) k4off = 1.0e-1 (/sec)

kcat1 = 1.0e3 kcat2 = 1.0e-3

S0 = 1.0e-5 (mol/l)

V0 = 1.0e-6 (mol/l)

A0 = 1.0e-6 (mol/l)

}

INDEPENDENT {

time FROM 0 TO 300 WITH 300 (sec) }

PLOT BoundV VS time

ASSIGNED {

BoundV FreeV

}

STATE {

S V A SV SA SVA SP SPV SPA SPVA

}

INITIAL {

SV=0

SA=0

SVA=0

SP=0

SPV=0

SPA=0

SPVA=0

}

KINETIC StoSVA {

~S + V <-> SV (k1on,k1off)

~S + A <-> SA (k2on,k2off)

~SV + A <-> SVA (k2on,k2off)

~SA + V <-> SVA (k1on,k1off)

~S <-> SP (kcat1,kcat2)

~SP + V <-> SPV (k3on,k3off)

~SP + A <-> SPA (k4on,k4off)

~SPV + A <-> SPVA (k4on,k4off)

~SPA + V <-> SPVA (k3on,k3off)

~SA <-> SPA (kcat1,kcat2)

~SVA <-> SPVA (kcat1,kcat2)

}

BREAKPOINT {

SOLVE StoSVA

BoundV = (SVA + SPVA)*1e7

FreeV = (V + SV +SPV)*1E7

}