#bray1



#reactions
complex_r1:
 {1.0}TTW = {1.0}TT + {1.0}W
  cell * (k1 * TTW - k2 * TT * W)
complex_r10:
 {1.0}TTWW = {1.0}TTW + {1.0}W
  cell * (k1 * TTWW - k2 * TTW * W)
complex_r11:
 {1.0}WWAA = {1.0}W + {1.0}WAA
  cell * (k1 * WWAA - k2 * W * WAA)
complex_r12:
 {1.0}TTWAA = {1.0}TTAA + {1.0}W
  cell * (k1 * TTWAA - k2 * TTAA * W)
complex_r2:
 {1.0}WAA = {1.0}W + {1.0}AA
  cell * (k1 * WAA - k2 * W * AA)
complex_r3:
 {1.0}TTWAA = {1.0}TT + {1.0}WAA
  cell * (k1 * TTWAA - k2 * TT * WAA)
complex_r4:
 {1.0}TTWWAA = {1.0}TTW + {1.0}WAA
  cell * (k1 * TTWWAA - k2 * TTW * WAA)
complex_r5:
 {1.0}TTWWAA = {1.0}TTWW + {1.0}AA
  cell * (k1 * TTWWAA - k2 * TTWW * AA)
complex_r6:
 {1.0}TTWWAA = {1.0}TT + {1.0}WWAA
  cell * (k1 * TTWWAA - k2 * TT * WWAA)
complex_r7:
 {1.0}TTAA = {1.0}TT + {1.0}AA
  cell * (k1 * TTAA - k2 * TT * AA)
complex_r8:
 {1.0}TTWAA = {1.0}TTW + {1.0}AA
  cell * (k1 * TTWAA - k2 * TTW * AA)
complex_r9:
 {1.0}TTWWAA = {1.0}TTWAA + {1.0}W
  cell * (k1 * TTWWAA - k2 * TTWAA * W)
phosphorylation_r1:
 {1.0}TTWWAA = {1.0}TTWWAAp
  cell * k1 * TTWWAA
phosphorylation_r10:
 {1.0}Bp = {1.0}B
  cell * k1 * Bp
phosphorylation_r2:
 {1.0}AA = {1.0}AAp
  cell * k1 * AA
phosphorylation_r3:
 {1.0}WAA = {1.0}WAAp
  cell * k1 * WAA
phosphorylation_r4:
 {1.0}WWAA = {1.0}WWAAp
  cell * k1 * WWAA
phosphorylation_r5:
 {1.0}TTAA = {1.0}TTAAp
  cell * k1 * TTAA
phosphorylation_r6:
 {1.0}TTWAA = {1.0}TTWAAp
  cell * k1 * TTWAA
phosphorylation_r7:
 {1.0}Y = {1.0}Yp
  cell * k1 * Y
phosphorylation_r8:
 {1.0}Yp = {1.0}Y
  cell * k1 * Yp
phosphorylation_r9:
 {1.0}Yp + {1.0}Z = {1.0}Y + {1.0}Z
  cell * k1 * Yp * Z
phosphotransfer_r1:
 {1.0}B + {1.0}AAp = {1.0}Bp + {1.0}AA
  cell * k1 * B * AAp
phosphotransfer_r10:
 {1.0}Y + {1.0}TTAAp = {1.0}Yp + {1.0}TTAA
  cell * k1 * Y * TTAAp
phosphotransfer_r11:
 {1.0}Y + {1.0}TTWAAp = {1.0}Yp + {1.0}TTWAA
  cell * k1 * Y * TTWAAp
phosphotransfer_r12:
 {1.0}Y + {1.0}TTWWAAp = {1.0}Yp + {1.0}TTWWAA
  cell * k1 * Y * TTWWAAp
phosphotransfer_r2:
 {1.0}B + {1.0}WAAp = {1.0}Bp + {1.0}WAA
  cell * k1 * B * WAAp
phosphotransfer_r3:
 {1.0}B + {1.0}WWAAp = {1.0}Bp + {1.0}WWAA
  cell * k1 * B * WWAAp
phosphotransfer_r4:
 {1.0}B + {1.0}TTAAp = {1.0}Bp + {1.0}TTAA
  cell * k1 * B * TTAAp
phosphotransfer_r5:
 {1.0}B + {1.0}TTWAAp = {1.0}Bp + {1.0}TTWAA
  cell * k1 * B * TTWAAp
phosphotransfer_r6:
 {1.0}B + {1.0}TTWWAAp = {1.0}Bp + {1.0}TTWWAA
  cell * k1 * B * TTWWAAp
phosphotransfer_r7:
 {1.0}Y + {1.0}AAp = {1.0}Yp + {1.0}AA
  cell * k1 * Y * AAp
phosphotransfer_r8:
 {1.0}Y + {1.0}WAAp = {1.0}Yp + {1.0}WAA
  cell * k1 * Y * WAAp
phosphotransfer_r9:
 {1.0}Y + {1.0}WWAAp = {1.0}Yp + {1.0}WWAA
  cell * k1 * Y * WWAAp


#parameters
Hill = 4.0
complex_r1_k1 = 0.00365
complex_r1_k2 = 1000000.0
complex_r2_k1 = 0.00894
complex_r2_k2 = 1000000.0
complex_r3_k1 = 297.0
complex_r3_k2 = 1000000.0
complex_r4_k1 = 0.64
complex_r4_k2 = 1000000.0
complex_r5_k1 = 0.112
complex_r5_k2 = 1000000.0
complex_r6_k1 = 0.0229
complex_r6_k2 = 1000000.0
complex_r7_k1 = 39.3
complex_r7_k2 = 1000000.0
complex_r8_k1 = 727.0
complex_r8_k2 = 1000000.0
complex_r9_k1 = 7.87e-06
complex_r9_k2 = 1000000.0
complex_r10_k1 = 0.0511
complex_r10_k2 = 1000000.0
complex_r11_k1 = 0.102
complex_r11_k2 = 1000000.0
complex_r12_k1 = 0.0676
complex_r12_k2 = 1000000.0
phosphorylation_r1_k1 = 15.5
phosphorylation_r2_k1 = 0.0227
phosphorylation_r3_k1 = 0.0227
phosphorylation_r4_k1 = 0.0227
phosphorylation_r5_k1 = 0.0227
phosphorylation_r6_k1 = 0.0227
phosphorylation_r7_k1 = 0.00124
phosphorylation_r8_k1 = 0.037
phosphorylation_r9_k1 = 500000.0
phosphorylation_r10_k1 = 0.35
phosphotransfer_r1_k1 = 6000000.0
phosphotransfer_r2_k1 = 6000000.0
phosphotransfer_r3_k1 = 6000000.0
phosphotransfer_r4_k1 = 6000000.0
phosphotransfer_r5_k1 = 6000000.0
phosphotransfer_r6_k1 = 6000000.0
phosphotransfer_r7_k1 = 30000000.0
phosphotransfer_r8_k1 = 30000000.0
phosphotransfer_r9_k1 = 30000000.0
phosphotransfer_r10_k1 = 30000000.0
phosphotransfer_r11_k1 = 30000000.0
phosphotransfer_r12_k1 = 30000000.0
cell = 1.41


#initial values
AA = 2.5e-06
AAp = 0.0
B = 2e-06
Bp = 0.0
SetYp = 1.63e-06
TT = 2.5e-06
TTAA = 0.0
TTAAp = 0.0
TTW = 0.0
TTWAA = 0.0
TTWAAp = 0.0
TTWW = 0.0
TTWWAA = 0.0
TTWWAAp = 0.0
W = 5e-06
WAA = 0.0
WAAp = 0.0
WWAA = 0.0
WWAAp = 0.0
Y = 1e-05
Yp = 0.0
Z = 2e-05


#assignment rules
!F Bias = 1.0 - pow(Yp, Hill) / (2.333 * pow(SetYp, Hill) + pow(Yp, Hill))