begin name clancy1 end name begin reactions reaction0000001 {1.0}c3 = {1.0}c2 reaction0000002 {1.0}c2 = {1.0}c1 reaction0000003 {1.0}c1 = {1.0}o reaction0000004 {1.0}o = {1.0}i reaction0000005 {1.0}c1 = {1.0}i end reactions begin rate equations reaction0000001 = cell*(-(b*c2[t]) + a*c3[t]) reaction0000002 = cell*(-(bin*c1[t]) + ain*c2[t]) reaction0000003 = cell*(aa*c1[t] - bb*o[t]) reaction0000004 = cell*(-(ai*i[t]) + bi*o[t]) reaction0000005 = cell*(aa*c1[t] - u*i[t]) end rate equations begin parameters F = 96.485 R = 8.314 Temp = 310.0 ain = 2.172 bin = 1.077 ki = 140.0 ko = 5.4 vhold = -40.0 vtest = 0.0 cell = 1.0 end parameters begin initial conditions c1[0] = c1i c2[0] = c2i c3[0] = c3i i[0] = ii o[0] = oi end initial conditions begin initial values c1i = 0.0 c2i = 0.0 c3i = 1.0 ii = 0.0 oi = 0.06 end initial values begin assignment rules aa := 0.0655*E^(0.05547153*(-36 + v)) ik := Gk*(v - vk)*o[t] a := 0.0555*E^(0.05547153*(-12 + v)) b := 0.002357/E^(0.036588*v) Gk := 0.0135*ko^0.59 bi := 1.03007202707417*E^(0.000942*v)*(ko^(-1))^0.3 u := (ai*bb)/bi bb := 0.0029357/E^(0.02158*v) vk := (R*Temp*Log[ko/ki])/F ai := 1.9755/(E^(0.02352*(25 + v))*ko) v := Piecewise[{{vtest, 50 < t && t <= 300}}, vhold] end assignment rules begin function definitions end function definitions begin events end events begin process annotations end process annotations begin species annotations end species annotations begin units end units