begin name
fernandez1
end name

begin reactions
v57 {1.0}$Empty = {1.0}Ca
v58 {1.0}Ca = {1.0}$Empty
vcat1 {1.0}DCDK5 = {1.0}D75 + {1.0}CDK5
vcat10 {1.0}D75PP2AP = {1.0}D + {1.0}PP2AP
vcat11 {1.0}D137CDK5 = {1.0}D75137 + {1.0}CDK5
vcat12 {1.0}D137PKA = {1.0}D34137 + {1.0}PKA
vcat13 {1.0}D137PP2C = {1.0}D + {1.0}PP2C
vcat14 {1.0}D3475CK1 = {1.0}D3475137 + {1.0}CK1
vcat15 {1.0}D3475PP2A = {1.0}D34 + {1.0}PP2A
vcat16 {1.0}D3475PP2AP = {1.0}D34 + {1.0}PP2AP
vcat17 {1.0}D3475PP2B = {1.0}D75 + {1.0}PP2B
vcat18 {1.0}D34137CDK5 = {1.0}D3475137 + {1.0}CDK5
vcat19 {1.0}D34137PP2B = {1.0}D137 + {1.0}PP2B
vcat2 {1.0}DCK1 = {1.0}D137 + {1.0}CK1
vcat20 {1.0}D34137PP2C = {1.0}D34 + {1.0}PP2C
vcat21 {1.0}D75137PKA = {1.0}D3475137 + {1.0}PKA
vcat22 {1.0}D75137PP2A = {1.0}D137 + {1.0}PP2A
vcat23 {1.0}D75137PP2AP = {1.0}D137 + {1.0}PP2AP
vcat24 {1.0}D75137PP2C = {1.0}D75 + {1.0}PP2C
vcat25 {1.0}D3475134PP2A = {1.0}D34137 + {1.0}PP2A
vcat26 {1.0}D3475137PP2AP = {1.0}D34137 + {1.0}PP2AP
vcat27 {1.0}D3475137PP2B = {1.0}D75137 + {1.0}PP2B
vcat28 {1.0}D3475137PP2C = {1.0}D3475 + {1.0}PP2C
vcat29 {1.0}CK1PPP2B = {1.0}CK1 + {1.0}PP2B
vcat3 {1.0}DPKA = {1.0}D34 + {1.0}PKA
vcat30 {1.0}CK1 = {1.0}CK1P
vcat31 {1.0}PDEPKA = {1.0}PDEP + {1.0}PKA
vcat32 {1.0}PDEP = {1.0}PDE
vcat33 {1.0}PP2APKA = {1.0}PP2AP + {1.0}PKA
vcat34 {1.0}PP2AP = {1.0}PP2A
vcat4 {1.0}D34CDK5 = {1.0}D3475 + {1.0}CDK5
vcat44 {1.0}cAMPPDE = {1.0}AMP + {1.0}PDE
vcat45 {1.0}cAMPPDEP = {1.0}AMP + {1.0}PDEP
vcat5 {1.0}D34CK1 = {1.0}D34137 + {1.0}CK1
vcat6 {1.0}D34PP2B = {1.0}D + {1.0}PP2B
vcat7 {1.0}D75CK1 = {1.0}CK1 + {1.0}D75137
vcat8 {1.0}D75PKA = {1.0}D3475 + {1.0}PKA
vcat9 {1.0}D75PP2A = {1.0}D + {1.0}PP2A
voff1 {1.0}DCDK5 = {1.0}D + {1.0}CDK5
voff10 {1.0}D75PP2AP = {1.0}D75 + {1.0}PP2AP
voff11 {1.0}D137CDK5 = {1.0}D137 + {1.0}CDK5
voff12 {1.0}D137PKA = {1.0}D137 + {1.0}PKA
voff13 {1.0}D137PP2C = {1.0}D137 + {1.0}PP2C
voff14 {1.0}D3475CK1 = {1.0}D3475 + {1.0}CK1
voff15 {1.0}D3475PP2A = {1.0}D3475 + {1.0}PP2A
voff16 {1.0}D3475PP2AP = {1.0}D3475 + {1.0}PP2AP
voff17 {1.0}D3475PP2B = {1.0}D3475 + {1.0}PP2B
voff18 {1.0}D34137CDK5 = {1.0}D34137 + {1.0}CDK5
voff19 {1.0}D34137PP2B = {1.0}D34137 + {1.0}PP2B
voff2 {1.0}DCK1 = {1.0}D + {1.0}CK1
voff20 {1.0}D34137PP2C = {1.0}D34137 + {1.0}PP2C
voff21 {1.0}D75137PKA = {1.0}D75137 + {1.0}PKA
voff22 {1.0}D75137PP2A = {1.0}D75137 + {1.0}PP2A
voff23 {1.0}D75137PP2AP = {1.0}D75137 + {1.0}PP2AP
voff24 {1.0}D75137PP2C = {1.0}D75137 + {1.0}PP2C
voff25 {1.0}D3475134PP2A = {1.0}D3475137 + {1.0}PP2A
voff26 {1.0}D3475137PP2AP = {1.0}D3475137 + {1.0}PP2AP
voff27 {1.0}D3475137PP2B = {1.0}D3475137 + {1.0}PP2B
voff28 {1.0}D3475137PP2C = {1.0}D3475137 + {1.0}PP2C
voff29 {1.0}CK1PPP2B = {1.0}CK1P + {1.0}PP2B
voff3 {1.0}DPKA = {1.0}D + {1.0}PKA
voff31 {1.0}PDEPKA = {1.0}PDE + {1.0}PKA
voff33 {1.0}PP2APKA = {1.0}PP2A + {1.0}PKA
voff35 {1.0}PP2BinactiveCa2 = {1.0}PP2Binactive + {2.0}Ca
voff36 {1.0}PP2B = {1.0}PP2BinactiveCa2 + {2.0}Ca
voff37 {1.0}cAMPR2C2 = {1.0}R2C2 + {1.0}cAMP
voff38 {1.0}cAMP2R2C2 = {1.0}cAMPR2C2 + {1.0}cAMP
voff39 {1.0}cAMP3R2C2 = {1.0}cAMP2R2C2 + {1.0}cAMP
voff4 {1.0}D34CDK5 = {1.0}D34 + {1.0}CDK5
voff40 {1.0}cAMP4R2C2 = {1.0}cAMP3R2C2 + {1.0}cAMP
voff41 {1.0}cAMP4R2C2 = {1.0}cAMP4R2C + {1.0}PKA
voff44 {1.0}cAMPPDE = {1.0}cAMP + {1.0}PDE
voff45 {1.0}cAMPPDEP = {1.0}cAMP + {1.0}PDEP
voff5 {1.0}D34CK1 = {1.0}D34 + {1.0}CK1
voff6 {1.0}D34PP2B = {1.0}D34 + {1.0}PP2B
voff7 {1.0}D75CK1 = {1.0}D75 + {1.0}CK1
voff8 {1.0}D75PKA = {1.0}D75 + {1.0}PKA
voff9 {1.0}D75PP2A = {1.0}D75 + {1.0}PP2A
von1 {1.0}D + {1.0}CDK5 = {1.0}DCDK5
von10 {1.0}D75 + {1.0}PP2AP = {1.0}D75PP2AP
von11 {1.0}D137 + {1.0}CDK5 = {1.0}D137CDK5
von12 {1.0}D137 + {1.0}PKA = {1.0}D137PKA
von13 {1.0}D137 + {1.0}PP2C = {1.0}D137PP2C
von14 {1.0}D3475 + {1.0}CK1 = {1.0}D3475CK1
von15 {1.0}D3475 + {1.0}PP2A = {1.0}D3475PP2A
von16 {1.0}D3475 + {1.0}PP2AP = {1.0}D3475PP2AP
von17 {1.0}D3475 + {1.0}PP2B = {1.0}D3475PP2B
von18 {1.0}D34137 + {1.0}CDK5 = {1.0}D34137CDK5
von19 {1.0}D34137 + {1.0}PP2B = {1.0}D34137PP2B
von2 {1.0}D + {1.0}CK1 = {1.0}DCK1
von20 {1.0}D34137 + {1.0}PP2C = {1.0}D34137PP2C
von21 {1.0}D75137 + {1.0}PKA = {1.0}D75137PKA
von22 {1.0}D75137 + {1.0}PP2A = {1.0}D75137PP2A
von23 {1.0}D75137 + {1.0}PP2AP = {1.0}D75137PP2AP
von24 {1.0}D75137 + {1.0}PP2C = {1.0}D75137PP2C
von25 {1.0}D3475137 + {1.0}PP2A = {1.0}D3475134PP2A
von26 {1.0}D3475137 + {1.0}PP2AP = {1.0}D3475137PP2AP
von27 {1.0}D3475137 + {1.0}PP2B = {1.0}D3475137PP2B
von28 {1.0}D3475137 + {1.0}PP2C = {1.0}D3475137PP2C
von29 {1.0}CK1P + {1.0}PP2B = {1.0}CK1PPP2B
von3 {1.0}D + {1.0}PKA = {1.0}DPKA
von31 {1.0}PDE + {1.0}PKA = {1.0}PDEPKA
von33 {1.0}PP2A + {1.0}PKA = {1.0}PP2APKA
von35 {1.0}PP2Binactive + {2.0}Ca = {1.0}PP2BinactiveCa2
von36 {1.0}PP2BinactiveCa2 + {2.0}Ca = {1.0}PP2B
von37 {1.0}R2C2 + {1.0}cAMP = {1.0}cAMPR2C2
von38 {1.0}cAMPR2C2 + {1.0}cAMP = {1.0}cAMP2R2C2
von39 {1.0}cAMP2R2C2 + {1.0}cAMP = {1.0}cAMP3R2C2
von4 {1.0}D34 + {1.0}CDK5 = {1.0}D34CDK5
von40 {1.0}cAMP3R2C2 + {1.0}cAMP = {1.0}cAMP4R2C2
von41 {1.0}cAMP4R2C + {1.0}PKA = {1.0}cAMP4R2C2
von42 {1.0}cAMP4R2 + {1.0}PKA = {1.0}cAMP4R2C
von43 {1.0}cAMP4R2C = {1.0}cAMP4R2 + {1.0}PKA
von44 {1.0}cAMP + {1.0}PDE = {1.0}cAMPPDE
von45 {1.0}cAMP + {1.0}PDEP = {1.0}cAMPPDEP
von5 {1.0}D34 + {1.0}CK1 = {1.0}D34CK1
von6 {1.0}D34 + {1.0}PP2B = {1.0}D34PP2B
von7 {1.0}D75 + {1.0}CK1 = {1.0}D75CK1
von8 {1.0}D75 + {1.0}PKA = {1.0}D75PKA
von9 {1.0}D75 + {1.0}PP2A = {1.0}D75PP2A
end reactions

begin rate equations
v57 = k57*Spine
v58 = Spine*v58k58*Ca[t]
vcat1 = Spine*vcat1kcat1*DCDK5[t]
vcat10 = Spine*vcat10kcat10*D75PP2AP[t]
vcat11 = Spine*vcat11kcat11*D137CDK5[t]
vcat12 = Spine*vcat12kcat12*D137PKA[t]
vcat13 = Spine*vcat13kcat13*D137PP2C[t]
vcat14 = Spine*vcat14kcat14*D3475CK1[t]
vcat15 = Spine*vcat15kcat15*D3475PP2A[t]
vcat16 = Spine*vcat16kcat16*D3475PP2AP[t]
vcat17 = Spine*vcat17kcat17*D3475PP2B[t]
vcat18 = Spine*vcat18kcat18*D34137CDK5[t]
vcat19 = Spine*vcat19kcat19*D34137PP2B[t]
vcat2 = Spine*vcat2kcat2*DCK1[t]
vcat20 = Spine*vcat20kcat20*D34137PP2C[t]
vcat21 = Spine*vcat21kcat21*D75137PKA[t]
vcat22 = Spine*vcat22kcat22*D75137PP2A[t]
vcat23 = Spine*vcat23kcat23*D75137PP2AP[t]
vcat24 = Spine*vcat24kcat24*D75137PP2C[t]
vcat25 = Spine*vcat25kcat25*D3475134PP2A[t]
vcat26 = Spine*vcat26kcat26*D3475137PP2AP[t]
vcat27 = Spine*vcat27kcat27*D3475137PP2B[t]
vcat28 = Spine*vcat28kcat28*D3475137PP2C[t]
vcat29 = Spine*vcat29kcat29*CK1PPP2B[t]
vcat3 = Spine*vcat3kcat3*DPKA[t]
vcat30 = Spine*vcat30kcat30*CK1[t]
vcat31 = Spine*vcat31kcat31*PDEPKA[t]
vcat32 = Spine*vcat32kcat32*PDEP[t]
vcat33 = Spine*vcat33kcat33*PP2APKA[t]
vcat34 = Spine*vcat34kcat34*PP2AP[t]
vcat4 = Spine*vcat4kcat4*D34CDK5[t]
vcat44 = Spine*vcat44kcat44*cAMPPDE[t]
vcat45 = Spine*vcat45kcat45*cAMPPDEP[t]
vcat5 = Spine*vcat5kcat5*D34CK1[t]
vcat6 = Spine*vcat6kcat6*D34PP2B[t]
vcat7 = Spine*vcat7kcat7*D75CK1[t]
vcat8 = Spine*vcat8kcat8*D75PKA[t]
vcat9 = Spine*vcat9kcat9*D75PP2A[t]
voff1 = Spine*voff1koff1*DCDK5[t]
voff10 = Spine*voff10koff10*D75PP2AP[t]
voff11 = Spine*voff11koff11*D137CDK5[t]
voff12 = Spine*voff12koff12*D137PKA[t]
voff13 = Spine*voff13koff13*D137PP2C[t]
voff14 = Spine*voff14koff14*D3475CK1[t]
voff15 = Spine*voff15koff15*D3475PP2A[t]
voff16 = Spine*voff16koff16*D3475PP2AP[t]
voff17 = Spine*voff17koff17*D3475PP2B[t]
voff18 = Spine*voff18koff18*D34137CDK5[t]
voff19 = Spine*voff19koff19*D34137PP2B[t]
voff2 = Spine*voff2koff2*DCK1[t]
voff20 = Spine*voff20koff20*D34137PP2C[t]
voff21 = Spine*voff21koff21*D75137PKA[t]
voff22 = Spine*voff22koff22*D75137PP2A[t]
voff23 = Spine*voff23koff23*D75137PP2AP[t]
voff24 = Spine*voff24koff24*D75137PP2C[t]
voff25 = Spine*voff25koff25*D3475134PP2A[t]
voff26 = Spine*voff26koff26*D3475137PP2AP[t]
voff27 = Spine*voff27koff27*D3475137PP2B[t]
voff28 = Spine*voff28koff28*D3475137PP2C[t]
voff29 = Spine*voff29koff29*CK1PPP2B[t]
voff3 = Spine*voff3koff3*DPKA[t]
voff31 = Spine*voff31koff31*PDEPKA[t]
voff33 = Spine*voff33koff33*PP2APKA[t]
voff35 = Spine*voff35koff35*PP2BinactiveCa2[t]
voff36 = Spine*voff36koff36*PP2B[t]
voff37 = Spine*voff37koff37*cAMPR2C2[t]
voff38 = Spine*voff38koff38*cAMP2R2C2[t]
voff39 = Spine*voff39koff39*cAMP3R2C2[t]
voff4 = Spine*voff4koff4*D34CDK5[t]
voff40 = Spine*voff40koff40*cAMP4R2C2[t]
voff41 = Spine*voff41koff41*cAMP4R2C2[t]
voff44 = Spine*voff44koff44*cAMPPDE[t]
voff45 = Spine*voff45koff45*cAMPPDEP[t]
voff5 = Spine*voff5koff5*D34CK1[t]
voff6 = Spine*voff6koff6*D34PP2B[t]
voff7 = Spine*voff7koff7*D75CK1[t]
voff8 = Spine*voff8koff8*D75PKA[t]
voff9 = Spine*voff9koff9*D75PP2A[t]
von1 = Spine*von1kon1*CDK5[t]*D[t]
von10 = Spine*von10kon10*D75[t]*PP2AP[t]
von11 = Spine*von11kon11*CDK5[t]*D137[t]
von12 = Spine*von12kon12*D137[t]*PKA[t]
von13 = Spine*von13kon13*D137[t]*PP2C[t]
von14 = Spine*von14kon14*CK1[t]*D3475[t]
von15 = Spine*von15kon15*D3475[t]*PP2A[t]
von16 = Spine*von16kon16*D3475[t]*PP2AP[t]
von17 = Spine*von17kon17*D3475[t]*PP2B[t]
von18 = Spine*von18kon18*CDK5[t]*D34137[t]
von19 = Spine*von19kon19*D34137[t]*PP2B[t]
von2 = Spine*von2kon2*CK1[t]*D[t]
von20 = Spine*von20kon20*D34137[t]*PP2C[t]
von21 = Spine*von21kon21*D75137[t]*PKA[t]
von22 = Spine*von22kon22*D75137[t]*PP2A[t]
von23 = Spine*von23kon23*D75137[t]*PP2AP[t]
von24 = Spine*von24kon24*D75137[t]*PP2C[t]
von25 = Spine*von25kon25*D3475137[t]*PP2A[t]
von26 = Spine*von26kon26*D3475137[t]*PP2AP[t]
von27 = Spine*von27kon27*D3475137[t]*PP2B[t]
von28 = Spine*von28kon28*D3475137[t]*PP2C[t]
von29 = Spine*von29kon29*CK1P[t]*PP2B[t]
von3 = Spine*von3kon3*D[t]*PKA[t]
von31 = Spine*von31kon31*PDE[t]*PKA[t]
von33 = Spine*von33kon33*PKA[t]*PP2A[t]
von35 = Spine*von35kon35*Ca[t]^2*PP2Binactive[t]
von36 = Spine*von36kon36*Ca[t]^2*PP2BinactiveCa2[t]
von37 = Spine*von37kon37*cAMP[t]*R2C2[t]
von38 = Spine*von38kon38*cAMP[t]*cAMPR2C2[t]
von39 = Spine*von39kon39*cAMP[t]*cAMP2R2C2[t]
von4 = Spine*von4kon4*CDK5[t]*D34[t]
von40 = Spine*von40kon40*cAMP[t]*cAMP3R2C2[t]
von41 = Spine*von41kon41*cAMP4R2C[t]*PKA[t]
von42 = Spine*von42kon42*cAMP4R2[t]*PKA[t]
von43 = Spine*von43kon43*cAMP4R2C[t]
von44 = Spine*von44kon44*cAMP[t]*PDE[t]
von45 = Spine*von45kon45*cAMP[t]*PDEP[t]
von5 = Spine*von5kon5*CK1[t]*D34[t]
von6 = Spine*von6kon6*D34[t]*PP2B[t]
von7 = Spine*von7kon7*CK1[t]*D75[t]
von8 = Spine*von8kon8*D75[t]*PKA[t]
von9 = Spine*von9kon9*D75[t]*PP2A[t]
end rate equations

begin parameters
cAMPCadelay = 50.0
cAMPdelay = 400.0
k57 = "2.5e-08"
konhigh = "6.6e-06"
konlow = "2.5e-08"
spikeduration = 2.0
spikeinterval = 2.0
Empty = 0.0
von1kon1 = 5600000.0
voff1koff1 = 12.0
vcat1kcat1 = 3.0
von2kon2 = 4400000.0
voff2koff2 = 12.0
vcat2kcat2 = 3.0
von3kon3 = 5600000.0
voff3koff3 = 10.8
vcat3kcat3 = 2.7
von4kon4 = 5600000.0
von5kon5 = 4400000.0
von6kon6 = 10000000.0
voff4koff4 = 12.0
vcat4kcat4 = 3.0
voff5koff5 = 12.0
vcat5kcat5 = 3.0
vcat6kcat6 = 4.0
voff6koff6 = 16.0
von7kon7 = 4400000.0
von8kon8 = 5600000.0
von9kon9 = 3800000.0
von10kon10 = 17000000.0
voff7koff7 = 12.0
vcat7kcat7 = 3.0
vcat8kcat8 = 0.0
voff8koff8 = 10.8
vcat9kcat9 = 10.0
voff9koff9 = 24.0
vcat10kcat10 = 24.0
voff10koff10 = 40.0
von11kon11 = 5600000.0
von12kon12 = 5600000.0
von13kon13 = 7500000.0
voff11koff11 = 12.0
vcat11kcat11 = 3.0
voff12koff12 = 10.8
vcat12kcat12 = 2.7
vcat13kcat13 = 3.0
voff13koff13 = 12.0
von14kon14 = 4400000.0
von18kon18 = 5600000.0
voff18koff18 = 12.0
voff14koff14 = 12.0
vcat14kcat14 = 3.0
vcat18kcat18 = 3.0
von21kon21 = 5600000.0
vcat21kcat21 = 0.0
voff21koff21 = 10.8
von17kon17 = 10000000.0
voff17koff17 = 1600.0
vcat17kcat17 = 4.0
von19kon19 = 75000.0
vcat19kcat19 = 0.03
voff19koff19 = 0.12
von27kon27 = 75000.0
voff27koff27 = 120.0
vcat27kcat27 = 0.03
von15kon15 = 3800000.0
vcat15kcat15 = 10.0
voff15koff15 = 24.0
von22kon22 = 3800000.0
vcat22kcat22 = 10.0
voff22koff22 = 24.0
von25kon25 = 3800000.0
vcat25kcat25 = 10.0
voff25koff25 = 24.0
von16kon16 = 17000000.0
vcat16kcat16 = 24.0
voff16koff16 = 40.0
von23kon23 = 17000000.0
vcat23kcat23 = 24.0
voff23koff23 = 40.0
vcat26kcat26 = 24.0
von26kon26 = 17000000.0
voff26koff26 = 40.0
von20kon20 = 7500000.0
vcat20kcat20 = 3.0
voff20koff20 = 12.0
von24kon24 = 7500000.0
vcat24kcat24 = 3.0
voff24koff24 = 12.0
von28kon28 = 7500000.0
vcat28kcat28 = 3.0
voff28koff28 = 12.0
von29kon29 = 30000000.0
voff29koff29 = 24.0
vcat29kcat29 = 6.0
vcat30kcat30 = 1.0
von31kon31 = 6000000.0
vcat31kcat31 = 9.0
voff31koff31 = 36.0
vcat32kcat32 = 0.1
von33kon33 = 10000000.0
voff33koff33 = 16.0
vcat33kcat33 = 4.0
vcat34kcat34 = 5.0
von35kon35 = 1000000000000000.0
von36kon36 = 3000000000000000.0
voff35koff35 = 1.0
voff36koff36 = 1.0
von37kon37 = 54000000.0
von38kon38 = 54000000.0
von39kon39 = 75000000.0
von40kon40 = 75000000.0
voff37koff37 = 33.0
voff38koff38 = 33.0
voff39koff39 = 110.0
voff40koff40 = 32.5
von41kon41 = 18000000.0
voff41koff41 = 60.0
von42kon42 = 18000000.0
von43kon43 = 60.0
von44kon44 = 2520000.0
voff44koff44 = 40.0
vcat44kcat44 = 10.0
von45kon45 = 5040000.0
voff45koff45 = 80.0
vcat45kcat45 = 20.0
v58k58 = 1.7
Spine = "1e-15"
end parameters

begin initial conditions
AMP[0] = AMPi
CDK5[0] = CDK5i
CK1[0] = CK1i
CK1P[0] = CK1Pi
CK1PPP2B[0] = CK1PPP2Bi
Ca[0] = Cai
D[0] = Di
D137[0] = D137i
D137CDK5[0] = D137CDK5i
D137PKA[0] = D137PKAi
D137PP2C[0] = D137PP2Ci
D34[0] = D34i
D34137[0] = D34137i
D34137CDK5[0] = D34137CDK5i
D34137PP2B[0] = D34137PP2Bi
D34137PP2C[0] = D34137PP2Ci
D3475[0] = D3475i
D3475134PP2A[0] = D3475134PP2Ai
D3475137[0] = D3475137i
D3475137PP2AP[0] = D3475137PP2APi
D3475137PP2B[0] = D3475137PP2Bi
D3475137PP2C[0] = D3475137PP2Ci
D3475CK1[0] = D3475CK1i
D3475PP2A[0] = D3475PP2Ai
D3475PP2AP[0] = D3475PP2APi
D3475PP2B[0] = D3475PP2Bi
D34CDK5[0] = D34CDK5i
D34CK1[0] = D34CK1i
D34PP2B[0] = D34PP2Bi
D75[0] = D75i
D75CK1[0] = D75CK1i
D75137[0] = D75137i
D75137PKA[0] = D75137PKAi
D75137PP2A[0] = D75137PP2Ai
D75137PP2AP[0] = D75137PP2APi
D75137PP2C[0] = D75137PP2Ci
D75PKA[0] = D75PKAi
D75PP2A[0] = D75PP2Ai
D75PP2AP[0] = D75PP2APi
DCDK5[0] = DCDK5i
DCK1[0] = DCK1i
DPKA[0] = DPKAi
PDE[0] = PDEi
PDEP[0] = PDEPi
PDEPKA[0] = PDEPKAi
PKA[0] = PKAi
PP2A[0] = PP2Ai
PP2AP[0] = PP2APi
PP2APKA[0] = PP2APKAi
PP2B[0] = PP2Bi
PP2Binactive[0] = PP2Binactivei
PP2BinactiveCa2[0] = PP2BinactiveCa2i
PP2C[0] = PP2Ci
R2C2[0] = R2C2i
cAMP[0] = cAMPi
cAMP2R2C2[0] = cAMP2R2C2i
cAMP3R2C2[0] = cAMP3R2C2i
cAMP4R2[0] = cAMP4R2i
cAMP4R2C[0] = cAMP4R2Ci
cAMP4R2C2[0] = cAMP4R2C2i
cAMPPDE[0] = cAMPPDEi
cAMPPDEP[0] = cAMPPDEPi
cAMPR2C2[0] = cAMPR2C2i
end initial conditions

begin initial values
AMPi = 0.0
CDK5i = 2e-07
CK1i = 1.66e-07
CK1Pi = 0.0
CK1PPP2Bi = 0.0
Cai = 0.0
Di = 4.98e-06
D137i = 0.0
D137CDK5i = 0.0
D137PKAi = 0.0
D137PP2Ci = 0.0
D34i = 0.0
D34137i = 0.0
D34137CDK5i = 0.0
D34137PP2Bi = 0.0
D34137PP2Ci = 0.0
D3475i = 0.0
D3475134PP2Ai = 0.0
D3475137i = 0.0
D3475137PP2APi = 0.0
D3475137PP2Bi = 0.0
D3475137PP2Ci = 0.0
D3475CK1i = 0.0
D3475PP2Ai = 0.0
D3475PP2APi = 0.0
D3475PP2Bi = 0.0
D34CDK5i = 0.0
D34CK1i = 0.0
D34PP2Bi = 0.0
D75i = 0.0
D75CK1i = 0.0
D75137i = 0.0
D75137PKAi = 0.0
D75137PP2Ai = 0.0
D75137PP2APi = 0.0
D75137PP2Ci = 0.0
D75PKAi = 0.0
D75PP2Ai = 0.0
D75PP2APi = 0.0
DCDK5i = 0.0
DCK1i = 0.0
DPKAi = 0.0
PDEi = 2e-06
PDEPi = 0.0
PDEPKAi = 0.0
PKAi = 0.0
PP2Ai = 2e-07
PP2APi = 0.0
PP2APKAi = 0.0
PP2Bi = 0.0
PP2Binactivei = 3.32e-07
PP2BinactiveCa2i = 0.0
PP2Ci = 1.33e-07
R2C2i = 6.64e-06
cAMPi = 0.0
cAMP2R2C2i = 0.0
cAMP3R2C2i = 0.0
cAMP4R2i = 0.0
cAMP4R2Ci = 0.0
cAMP4R2C2i = 0.0
cAMPPDEi = 0.0
cAMPPDEPi = 0.0
cAMPR2C2i = 0.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
ca_on2 = (t >= cAMPCadelay + cAMPdelay + spikeduration + spikeinterval, {k57 -> konhigh}, 0)
ca_on9 = (t >= cAMPCadelay + cAMPdelay + 8*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_off3 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 2*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_off4 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 3*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_on5 = (t >= cAMPCadelay + cAMPdelay + 4*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_on10 = (t >= cAMPCadelay + cAMPdelay + 9*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_off2 = (t >= cAMPCadelay + cAMPdelay + 2*spikeduration + spikeinterval, {k57 -> konlow}, 0)
cAMP_pulse = (t >= cAMPdelay, {cAMP[t] -> 6.5999999999999995*^-6}, 0)
ca_off1 = (t >= cAMPCadelay + cAMPdelay + spikeduration, {k57 -> konlow}, 0)
Ca_on3 = (t >= cAMPCadelay + cAMPdelay + 2*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_on6 = (t >= cAMPCadelay + cAMPdelay + 5*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_on4 = (t >= cAMPCadelay + cAMPdelay + 3*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_on1 = (t >= cAMPCadelay + cAMPdelay, {k57 -> konhigh}, 0)
ca_on8 = (t >= cAMPCadelay + cAMPdelay + 7*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_off6 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 5*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_off5 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 4*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_off10 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 9*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_on7 = (t >= cAMPCadelay + cAMPdelay + 6*(spikeduration + spikeinterval), {k57 -> konhigh}, 0)
ca_off9 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 8*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_off8 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 7*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
ca_off7 = (t >= cAMPCadelay + cAMPdelay + spikeduration + 6*(spikeduration + spikeinterval), {k57 -> konlow}, 0)
end events

begin process annotations
end process annotations

begin species annotations
AMP http://identifiers.org/chebi/CHEBI:16027
CDK5 http://identifiers.org/uniprot/Q03114
Ca http://identifiers.org/chebi/CHEBI:29108
PP2C http://identifiers.org/interpro/IPR015655
cAMP http://identifiers.org/chebi/CHEBI:17489
end species annotations

begin units
end units