//Created by libAntimony v2.4 model plant_1981__environment(time_) // Variable initializations: time_ = ; end model plant_1981__membrane(V, Vs, time_, c, sI, yI, xT, xK) // Assignment Rules: Vs := 1 * g_I * V + 1 * g_T; g_I := 1 * ((V_H_Na - V_H_K) / (V_I - V_K)); g_T := 1 * ((V_H_Na * V_K - V_I * V_H_K) / (V_I - V_K)); // Rate Rules: V' = (g_I * power(sI, 3.0) * yI + g_T * xT) * (V_I - V) + (g_K * power(xK, 4.0) + g_P * c * power(Kp + c, -1)) * (V_K - V) + g_L * (V_L - V); // Variable initializations: V = -54; V_I = 30.0; V_K = -75.0; V_L = -40.0; V_H_Na = 115.0; V_H_K = -12.0; g_K = 0.3; g_L = 0.003; g_P = 0.03; Kp = 0.5; time_ = ; c = ; sI = ; yI = ; xT = ; xK = ; end model plant_1981__sI_gate(sI, Vs) // Assignment Rules: sI := alpha_m / (alpha_m + beta_m); alpha_m := (0.1 * (50.0 - Vs)) / -exp((50.0 - Vs) / 10.0); beta_m := 4.0 * exp((25.0 - Vs) / 18.0); // Variable initializations: Vs = ; end model plant_1981__yI_gate(yI, Vs, time_) // Assignment Rules: ZI := alpha_h / (alpha_h + beta_h); alpha_h := 0.07 * exp((25.0 - Vs) / 20.0); beta_h := 1.0 / (exp((55.0 - Vs) / 10.0) + 1.0); tau_yI := 12.5 / (alpha_h + beta_h); // Rate Rules: yI' = (ZI - yI) / tau_yI; // Variable initializations: yI = 0.1; Vs = ; time_ = ; end model plant_1981__xT_gate(xT, Vs, time_) // Assignment Rules: sT := 1.0 / (exp(0.15 * (-50.0 - Vs)) + 1.0); // Rate Rules: xT' = (sT - xT) / tau_xT; // Variable initializations: xT = 0.1; tau_xT = 235.0; Vs = ; time_ = ; end model plant_1981__calcium_concentration(c, time_, V, xT) // Rate Rules: c' = rho * (K_c * xT * (V_Ca - V) - c); // Variable initializations: c = 0.1; V_Ca = 140.0; rho = 0.0003; K_c = 0.0085; time_ = ; V = ; xT = ; end model plant_1981__xK_gate(xK, Vs, time_) // Assignment Rules: alpha_n := (0.01 * (55.0 - Vs)) / (exp((55.0 - Vs) / 10.0) - 1.0); beta_n := 0.125 * exp((45.0 - Vs) / 80.0); tau_xK := 12.5 / (alpha_n + beta_n); sK := alpha_n / (alpha_n + beta_n); // Rate Rules: xK' = (sK - xK) / tau_xK; // Variable initializations: xK = 0.1; Vs = ; time_ = ; end model *plant_1981____main() // Sub-modules, and any changes to those submodules: environment: plant_1981__environment(time_); membrane: plant_1981__membrane(V, Vs, time_, c, sI, yI, xT, xK); sI_gate: plant_1981__sI_gate(sI, Vs); yI_gate: plant_1981__yI_gate(yI, Vs, time_); xT_gate: plant_1981__xT_gate(xT, Vs, time_); calcium_concentration: plant_1981__calcium_concentration(c, time_, V, xT); xK_gate: plant_1981__xK_gate(xK, Vs, time_); end