#fernandez1

FIX: Empty

#reactions
v57:
 {1.0}Empty = {1.0}Ca
  Spine * k57
v58:
 {1.0}Ca = {1.0}Empty
  Spine * Ca * k58
vcat1:
 {1.0}D_CDK5 = {1.0}D75 + {1.0}CDK5
  Spine * D_CDK5 * kcat1
vcat10:
 {1.0}D75_PP2AP = {1.0}D + {1.0}PP2AP
  Spine * D75_PP2AP * kcat10
vcat11:
 {1.0}D137_CDK5 = {1.0}D75_137 + {1.0}CDK5
  Spine * D137_CDK5 * kcat11
vcat12:
 {1.0}D137_PKA = {1.0}D34_137 + {1.0}PKA
  Spine * D137_PKA * kcat12
vcat13:
 {1.0}D137_PP2C = {1.0}D + {1.0}PP2C
  Spine * D137_PP2C * kcat13
vcat14:
 {1.0}D34_75_CK1 = {1.0}D34_75_137 + {1.0}CK1
  Spine * D34_75_CK1 * kcat14
vcat15:
 {1.0}D34_75_PP2A = {1.0}D34 + {1.0}PP2A
  Spine * D34_75_PP2A * kcat15
vcat16:
 {1.0}D34_75_PP2AP = {1.0}D34 + {1.0}PP2AP
  Spine * D34_75_PP2AP * kcat16
vcat17:
 {1.0}D34_75_PP2B = {1.0}D75 + {1.0}PP2B
  Spine * D34_75_PP2B * kcat17
vcat18:
 {1.0}D34_137_CDK5 = {1.0}D34_75_137 + {1.0}CDK5
  Spine * D34_137_CDK5 * kcat18
vcat19:
 {1.0}D34_137_PP2B = {1.0}D137 + {1.0}PP2B
  Spine * D34_137_PP2B * kcat19
vcat2:
 {1.0}D_CK1 = {1.0}D137 + {1.0}CK1
  Spine * kcat2 * D_CK1
vcat20:
 {1.0}D34_137_PP2C = {1.0}D34 + {1.0}PP2C
  Spine * D34_137_PP2C * kcat20
vcat21:
 {1.0}D75_137_PKA = {1.0}D34_75_137 + {1.0}PKA
  Spine * D75_137_PKA * kcat21
vcat22:
 {1.0}D75_137_PP2A = {1.0}D137 + {1.0}PP2A
  Spine * D75_137_PP2A * kcat22
vcat23:
 {1.0}D75_137_PP2AP = {1.0}D137 + {1.0}PP2AP
  Spine * D75_137_PP2AP * kcat23
vcat24:
 {1.0}D75_137_PP2C = {1.0}D75 + {1.0}PP2C
  Spine * D75_137_PP2C * kcat24
vcat25:
 {1.0}D34_75_134_PP2A = {1.0}D34_137 + {1.0}PP2A
  Spine * D34_75_134_PP2A * kcat25
vcat26:
 {1.0}D34_75_137_PP2AP = {1.0}D34_137 + {1.0}PP2AP
  Spine * D34_75_137_PP2AP * kcat26
vcat27:
 {1.0}D34_75_137_PP2B = {1.0}D75_137 + {1.0}PP2B
  Spine * D34_75_137_PP2B * kcat27
vcat28:
 {1.0}D34_75_137_PP2C = {1.0}D34_75 + {1.0}PP2C
  Spine * D34_75_137_PP2C * kcat28
vcat29:
 {1.0}CK1P_PP2B = {1.0}CK1 + {1.0}PP2B
  Spine * CK1P_PP2B * kcat29
vcat3:
 {1.0}D_PKA = {1.0}D34 + {1.0}PKA
  Spine * D_PKA * kcat3
vcat30:
 {1.0}CK1 = {1.0}CK1P
  Spine * CK1 * kcat30
vcat31:
 {1.0}PDE_PKA = {1.0}PDEP + {1.0}PKA
  Spine * PDE_PKA * kcat31
vcat32:
 {1.0}PDEP = {1.0}PDE
  Spine * PDEP * kcat32
vcat33:
 {1.0}PP2A_PKA = {1.0}PP2AP + {1.0}PKA
  Spine * PP2A_PKA * kcat33
vcat34:
 {1.0}PP2AP = {1.0}PP2A
  Spine * PP2AP * kcat34
vcat4:
 {1.0}D34_CDK5 = {1.0}D34_75 + {1.0}CDK5
  Spine * D34_CDK5 * kcat4
vcat44:
 {1.0}cAMP_PDE = {1.0}AMP + {1.0}PDE
  Spine * cAMP_PDE * kcat44
vcat45:
 {1.0}cAMP_PDEP = {1.0}AMP + {1.0}PDEP
  Spine * cAMP_PDEP * kcat45
vcat5:
 {1.0}D34_CK1 = {1.0}D34_137 + {1.0}CK1
  Spine * D34_CK1 * kcat5
vcat6:
 {1.0}D34_PP2B = {1.0}D + {1.0}PP2B
  Spine * D34_PP2B * kcat6
vcat7:
 {1.0}D75CK1 = {1.0}CK1 + {1.0}D75_137
  Spine * D75CK1 * kcat7
vcat8:
 {1.0}D75_PKA = {1.0}D34_75 + {1.0}PKA
  Spine * D75_PKA * kcat8
vcat9:
 {1.0}D75_PP2A = {1.0}D + {1.0}PP2A
  Spine * D75_PP2A * kcat9
voff1:
 {1.0}D_CDK5 = {1.0}D + {1.0}CDK5
  Spine * D_CDK5 * koff1
voff10:
 {1.0}D75_PP2AP = {1.0}D75 + {1.0}PP2AP
  Spine * D75_PP2AP * koff10
voff11:
 {1.0}D137_CDK5 = {1.0}D137 + {1.0}CDK5
  Spine * D137_CDK5 * koff11
voff12:
 {1.0}D137_PKA = {1.0}D137 + {1.0}PKA
  Spine * D137_PKA * koff12
voff13:
 {1.0}D137_PP2C = {1.0}D137 + {1.0}PP2C
  Spine * D137_PP2C * koff13
voff14:
 {1.0}D34_75_CK1 = {1.0}D34_75 + {1.0}CK1
  Spine * D34_75_CK1 * koff14
voff15:
 {1.0}D34_75_PP2A = {1.0}D34_75 + {1.0}PP2A
  Spine * D34_75_PP2A * koff15
voff16:
 {1.0}D34_75_PP2AP = {1.0}D34_75 + {1.0}PP2AP
  Spine * D34_75_PP2AP * koff16
voff17:
 {1.0}D34_75_PP2B = {1.0}D34_75 + {1.0}PP2B
  Spine * D34_75_PP2B * koff17
voff18:
 {1.0}D34_137_CDK5 = {1.0}D34_137 + {1.0}CDK5
  Spine * D34_137_CDK5 * koff18
voff19:
 {1.0}D34_137_PP2B = {1.0}D34_137 + {1.0}PP2B
  Spine * D34_137_PP2B * koff19
voff2:
 {1.0}D_CK1 = {1.0}D + {1.0}CK1
  Spine * koff2 * D_CK1
voff20:
 {1.0}D34_137_PP2C = {1.0}D34_137 + {1.0}PP2C
  Spine * D34_137_PP2C * koff20
voff21:
 {1.0}D75_137_PKA = {1.0}D75_137 + {1.0}PKA
  Spine * D75_137_PKA * koff21
voff22:
 {1.0}D75_137_PP2A = {1.0}D75_137 + {1.0}PP2A
  Spine * D75_137_PP2A * koff22
voff23:
 {1.0}D75_137_PP2AP = {1.0}D75_137 + {1.0}PP2AP
  Spine * D75_137_PP2AP * koff23
voff24:
 {1.0}D75_137_PP2C = {1.0}D75_137 + {1.0}PP2C
  Spine * D75_137_PP2C * koff24
voff25:
 {1.0}D34_75_134_PP2A = {1.0}D34_75_137 + {1.0}PP2A
  Spine * D34_75_134_PP2A * koff25
voff26:
 {1.0}D34_75_137_PP2AP = {1.0}D34_75_137 + {1.0}PP2AP
  Spine * D34_75_137_PP2AP * koff26
voff27:
 {1.0}D34_75_137_PP2B = {1.0}D34_75_137 + {1.0}PP2B
  Spine * D34_75_137_PP2B * koff27
voff28:
 {1.0}D34_75_137_PP2C = {1.0}D34_75_137 + {1.0}PP2C
  Spine * D34_75_137_PP2C * koff28
voff29:
 {1.0}CK1P_PP2B = {1.0}CK1P + {1.0}PP2B
  Spine * CK1P_PP2B * koff29
voff3:
 {1.0}D_PKA = {1.0}D + {1.0}PKA
  Spine * D_PKA * koff3
voff31:
 {1.0}PDE_PKA = {1.0}PDE + {1.0}PKA
  Spine * PDE_PKA * koff31
voff33:
 {1.0}PP2A_PKA = {1.0}PP2A + {1.0}PKA
  Spine * PP2A_PKA * koff33
voff35:
 {1.0}PP2BinactiveCa2 = {1.0}PP2Binactive + {2.0}Ca
  Spine * PP2BinactiveCa2 * koff35
voff36:
 {1.0}PP2B = {1.0}PP2BinactiveCa2 + {2.0}Ca
  Spine * PP2B * koff36
voff37:
 {1.0}cAMP_R2C2 = {1.0}R2C2 + {1.0}cAMP
  Spine * cAMP_R2C2 * koff37
voff38:
 {1.0}cAMP2_R2C2 = {1.0}cAMP_R2C2 + {1.0}cAMP
  Spine * cAMP2_R2C2 * koff38
voff39:
 {1.0}cAMP3_R2C2 = {1.0}cAMP2_R2C2 + {1.0}cAMP
  Spine * cAMP3_R2C2 * koff39
voff4:
 {1.0}D34_CDK5 = {1.0}D34 + {1.0}CDK5
  Spine * D34_CDK5 * koff4
voff40:
 {1.0}cAMP4_R2C2 = {1.0}cAMP3_R2C2 + {1.0}cAMP
  Spine * cAMP4_R2C2 * koff40
voff41:
 {1.0}cAMP4_R2C2 = {1.0}cAMP4_R2C + {1.0}PKA
  Spine * cAMP4_R2C2 * koff41
voff44:
 {1.0}cAMP_PDE = {1.0}cAMP + {1.0}PDE
  Spine * cAMP_PDE * koff44
voff45:
 {1.0}cAMP_PDEP = {1.0}cAMP + {1.0}PDEP
  Spine * cAMP_PDEP * koff45
voff5:
 {1.0}D34_CK1 = {1.0}D34 + {1.0}CK1
  Spine * D34_CK1 * koff5
voff6:
 {1.0}D34_PP2B = {1.0}D34 + {1.0}PP2B
  Spine * D34_PP2B * koff6
voff7:
 {1.0}D75CK1 = {1.0}D75 + {1.0}CK1
  Spine * D75CK1 * koff7
voff8:
 {1.0}D75_PKA = {1.0}D75 + {1.0}PKA
  Spine * D75_PKA * koff8
voff9:
 {1.0}D75_PP2A = {1.0}D75 + {1.0}PP2A
  Spine * D75_PP2A * koff9
von1:
 {1.0}D + {1.0}CDK5 = {1.0}D_CDK5
  Spine * kon1 * D * CDK5
von10:
 {1.0}D75 + {1.0}PP2AP = {1.0}D75_PP2AP
  Spine * D75 * PP2AP * kon10
von11:
 {1.0}D137 + {1.0}CDK5 = {1.0}D137_CDK5
  Spine * D137 * CDK5 * kon11
von12:
 {1.0}D137 + {1.0}PKA = {1.0}D137_PKA
  Spine * D137 * PKA * kon12
von13:
 {1.0}D137 + {1.0}PP2C = {1.0}D137_PP2C
  Spine * D137 * PP2C * kon13
von14:
 {1.0}D34_75 + {1.0}CK1 = {1.0}D34_75_CK1
  Spine * D34_75 * CK1 * kon14
von15:
 {1.0}D34_75 + {1.0}PP2A = {1.0}D34_75_PP2A
  Spine * D34_75 * PP2A * kon15
von16:
 {1.0}D34_75 + {1.0}PP2AP = {1.0}D34_75_PP2AP
  Spine * D34_75 * PP2AP * kon16
von17:
 {1.0}D34_75 + {1.0}PP2B = {1.0}D34_75_PP2B
  Spine * D34_75 * PP2B * kon17
von18:
 {1.0}D34_137 + {1.0}CDK5 = {1.0}D34_137_CDK5
  Spine * D34_137 * CDK5 * kon18
von19:
 {1.0}D34_137 + {1.0}PP2B = {1.0}D34_137_PP2B
  Spine * D34_137 * PP2B * kon19
von2:
 {1.0}D + {1.0}CK1 = {1.0}D_CK1
  Spine * D * CK1 * kon2
von20:
 {1.0}D34_137 + {1.0}PP2C = {1.0}D34_137_PP2C
  Spine * D34_137 * PP2C * kon20
von21:
 {1.0}D75_137 + {1.0}PKA = {1.0}D75_137_PKA
  Spine * D75_137 * PKA * kon21
von22:
 {1.0}D75_137 + {1.0}PP2A = {1.0}D75_137_PP2A
  Spine * D75_137 * PP2A * kon22
von23:
 {1.0}D75_137 + {1.0}PP2AP = {1.0}D75_137_PP2AP
  Spine * D75_137 * PP2AP * kon23
von24:
 {1.0}D75_137 + {1.0}PP2C = {1.0}D75_137_PP2C
  Spine * D75_137 * PP2C * kon24
von25:
 {1.0}D34_75_137 + {1.0}PP2A = {1.0}D34_75_134_PP2A
  Spine * D34_75_137 * PP2A * kon25
von26:
 {1.0}D34_75_137 + {1.0}PP2AP = {1.0}D34_75_137_PP2AP
  Spine * D34_75_137 * PP2AP * kon26
von27:
 {1.0}D34_75_137 + {1.0}PP2B = {1.0}D34_75_137_PP2B
  Spine * D34_75_137 * PP2B * kon27
von28:
 {1.0}D34_75_137 + {1.0}PP2C = {1.0}D34_75_137_PP2C
  Spine * D34_75_137 * PP2C * kon28
von29:
 {1.0}CK1P + {1.0}PP2B = {1.0}CK1P_PP2B
  Spine * CK1P * PP2B * kon29
von3:
 {1.0}D + {1.0}PKA = {1.0}D_PKA
  Spine * D * PKA * kon3
von31:
 {1.0}PDE + {1.0}PKA = {1.0}PDE_PKA
  Spine * PDE * PKA * kon31
von33:
 {1.0}PP2A + {1.0}PKA = {1.0}PP2A_PKA
  Spine * PP2A * PKA * kon33
von35:
 {1.0}PP2Binactive + {2.0}Ca = {1.0}PP2BinactiveCa2
  Spine * PP2Binactive * Ca * Ca * kon35
von36:
 {1.0}PP2BinactiveCa2 + {2.0}Ca = {1.0}PP2B
  Spine * PP2BinactiveCa2 * Ca * Ca * kon36
von37:
 {1.0}R2C2 + {1.0}cAMP = {1.0}cAMP_R2C2
  Spine * R2C2 * cAMP * kon37
von38:
 {1.0}cAMP_R2C2 + {1.0}cAMP = {1.0}cAMP2_R2C2
  Spine * cAMP_R2C2 * cAMP * kon38
von39:
 {1.0}cAMP2_R2C2 + {1.0}cAMP = {1.0}cAMP3_R2C2
  Spine * cAMP2_R2C2 * cAMP * kon39
von4:
 {1.0}D34 + {1.0}CDK5 = {1.0}D34_CDK5
  Spine * D34 * CDK5 * kon4
von40:
 {1.0}cAMP3_R2C2 + {1.0}cAMP = {1.0}cAMP4_R2C2
  Spine * cAMP3_R2C2 * cAMP * kon40
von41:
 {1.0}cAMP4_R2C + {1.0}PKA = {1.0}cAMP4_R2C2
  Spine * cAMP4_R2C * PKA * kon41
von42:
 {1.0}cAMP4_R2 + {1.0}PKA = {1.0}cAMP4_R2C
  Spine * cAMP4_R2 * PKA * kon42
von43:
 {1.0}cAMP4_R2C = {1.0}cAMP4_R2 + {1.0}PKA
  Spine * cAMP4_R2C * kon43
von44:
 {1.0}cAMP + {1.0}PDE = {1.0}cAMP_PDE
  Spine * cAMP * PDE * kon44
von45:
 {1.0}cAMP + {1.0}PDEP = {1.0}cAMP_PDEP
  Spine * cAMP * PDEP * kon45
von5:
 {1.0}D34 + {1.0}CK1 = {1.0}D34_CK1
  Spine * D34 * CK1 * kon5
von6:
 {1.0}D34 + {1.0}PP2B = {1.0}D34_PP2B
  Spine * D34 * PP2B * kon6
von7:
 {1.0}D75 + {1.0}CK1 = {1.0}D75CK1
  Spine * D75 * CK1 * kon7
von8:
 {1.0}D75 + {1.0}PKA = {1.0}D75_PKA
  Spine * D75 * PKA * kon8
von9:
 {1.0}D75 + {1.0}PP2A = {1.0}D75_PP2A
  Spine * D75 * PP2A * kon9


#parameters
cAMP_Ca_delay = 50.0
cAMP_delay = 400.0
k57 = 2.5e-08
kon_high = 6.6e-06
kon_low = 2.5e-08
spike_duration = 2.0
spike_interval = 2.0
Empty = 0.0
von1_kon1 = 5600000.0
voff1_koff1 = 12.0
vcat1_kcat1 = 3.0
von2_kon2 = 4400000.0
voff2_koff2 = 12.0
vcat2_kcat2 = 3.0
von3_kon3 = 5600000.0
voff3_koff3 = 10.8
vcat3_kcat3 = 2.7
von4_kon4 = 5600000.0
von5_kon5 = 4400000.0
von6_kon6 = 10000000.0
voff4_koff4 = 12.0
vcat4_kcat4 = 3.0
voff5_koff5 = 12.0
vcat5_kcat5 = 3.0
vcat6_kcat6 = 4.0
voff6_koff6 = 16.0
von7_kon7 = 4400000.0
von8_kon8 = 5600000.0
von9_kon9 = 3800000.0
von10_kon10 = 17000000.0
voff7_koff7 = 12.0
vcat7_kcat7 = 3.0
vcat8_kcat8 = 0.0
voff8_koff8 = 10.8
vcat9_kcat9 = 10.0
voff9_koff9 = 24.0
vcat10_kcat10 = 24.0
voff10_koff10 = 40.0
von11_kon11 = 5600000.0
von12_kon12 = 5600000.0
von13_kon13 = 7500000.0
voff11_koff11 = 12.0
vcat11_kcat11 = 3.0
voff12_koff12 = 10.8
vcat12_kcat12 = 2.7
vcat13_kcat13 = 3.0
voff13_koff13 = 12.0
von14_kon14 = 4400000.0
von18_kon18 = 5600000.0
voff18_koff18 = 12.0
voff14_koff14 = 12.0
vcat14_kcat14 = 3.0
vcat18_kcat18 = 3.0
von21_kon21 = 5600000.0
vcat21_kcat21 = 0.0
voff21_koff21 = 10.8
von17_kon17 = 10000000.0
voff17_koff17 = 1600.0
vcat17_kcat17 = 4.0
von19_kon19 = 75000.0
vcat19_kcat19 = 0.03
voff19_koff19 = 0.12
von27_kon27 = 75000.0
voff27_koff27 = 120.0
vcat27_kcat27 = 0.03
von15_kon15 = 3800000.0
vcat15_kcat15 = 10.0
voff15_koff15 = 24.0
von22_kon22 = 3800000.0
vcat22_kcat22 = 10.0
voff22_koff22 = 24.0
von25_kon25 = 3800000.0
vcat25_kcat25 = 10.0
voff25_koff25 = 24.0
von16_kon16 = 17000000.0
vcat16_kcat16 = 24.0
voff16_koff16 = 40.0
von23_kon23 = 17000000.0
vcat23_kcat23 = 24.0
voff23_koff23 = 40.0
vcat26_kcat26 = 24.0
von26_kon26 = 17000000.0
voff26_koff26 = 40.0
von20_kon20 = 7500000.0
vcat20_kcat20 = 3.0
voff20_koff20 = 12.0
von24_kon24 = 7500000.0
vcat24_kcat24 = 3.0
voff24_koff24 = 12.0
von28_kon28 = 7500000.0
vcat28_kcat28 = 3.0
voff28_koff28 = 12.0
von29_kon29 = 30000000.0
voff29_koff29 = 24.0
vcat29_kcat29 = 6.0
vcat30_kcat30 = 1.0
von31_kon31 = 6000000.0
vcat31_kcat31 = 9.0
voff31_koff31 = 36.0
vcat32_kcat32 = 0.1
von33_kon33 = 10000000.0
voff33_koff33 = 16.0
vcat33_kcat33 = 4.0
vcat34_kcat34 = 5.0
von35_kon35 = 1000000000000000.0
von36_kon36 = 3000000000000000.0
voff35_koff35 = 1.0
voff36_koff36 = 1.0
von37_kon37 = 54000000.0
von38_kon38 = 54000000.0
von39_kon39 = 75000000.0
von40_kon40 = 75000000.0
voff37_koff37 = 33.0
voff38_koff38 = 33.0
voff39_koff39 = 110.0
voff40_koff40 = 32.5
von41_kon41 = 18000000.0
voff41_koff41 = 60.0
von42_kon42 = 18000000.0
von43_kon43 = 60.0
von44_kon44 = 2520000.0
voff44_koff44 = 40.0
vcat44_kcat44 = 10.0
von45_kon45 = 5040000.0
voff45_koff45 = 80.0
vcat45_kcat45 = 20.0
v58_k58 = 1.7
Spine = 1e-15


#initial values
AMP = 0.0
CDK5 = 2e-07
CK1 = 1.66e-07
CK1P = 0.0
CK1P_PP2B = 0.0
Ca = 0.0
D = 4.98e-06
D137 = 0.0
D137_CDK5 = 0.0
D137_PKA = 0.0
D137_PP2C = 0.0
D34 = 0.0
D34_137 = 0.0
D34_137_CDK5 = 0.0
D34_137_PP2B = 0.0
D34_137_PP2C = 0.0
D34_75 = 0.0
D34_75_134_PP2A = 0.0
D34_75_137 = 0.0
D34_75_137_PP2AP = 0.0
D34_75_137_PP2B = 0.0
D34_75_137_PP2C = 0.0
D34_75_CK1 = 0.0
D34_75_PP2A = 0.0
D34_75_PP2AP = 0.0
D34_75_PP2B = 0.0
D34_CDK5 = 0.0
D34_CK1 = 0.0
D34_PP2B = 0.0
D75 = 0.0
D75CK1 = 0.0
D75_137 = 0.0
D75_137_PKA = 0.0
D75_137_PP2A = 0.0
D75_137_PP2AP = 0.0
D75_137_PP2C = 0.0
D75_PKA = 0.0
D75_PP2A = 0.0
D75_PP2AP = 0.0
D_CDK5 = 0.0
D_CK1 = 0.0
D_PKA = 0.0
PDE = 2e-06
PDEP = 0.0
PDE_PKA = 0.0
PKA = 0.0
PP2A = 2e-07
PP2AP = 0.0
PP2A_PKA = 0.0
PP2B = 0.0
PP2Binactive = 3.32e-07
PP2BinactiveCa2 = 0.0
PP2C = 1.33e-07
R2C2 = 6.64e-06
cAMP = 0.0
cAMP2_R2C2 = 0.0
cAMP3_R2C2 = 0.0
cAMP4_R2 = 0.0
cAMP4_R2C = 0.0
cAMP4_R2C2 = 0.0
cAMP_PDE = 0.0
cAMP_PDEP = 0.0
cAMP_R2C2 = 0.0


#assignment rules