//Created by libAntimony v2.4 model swat_2004__environment1(time_) // Variable initializations: time_ = ; end model swat_2004__Core_Module(time_, pRB, E2F1, pRB_p, CycD_a, k_16, k_61) // Rate Rules: pRB' = ((k_1 * (E2F1 / (Km_1 + E2F1)) * (J_11 / (J_11 + pRB)) * (J_61 / (J_61 + pRB_p)) - k_16 * pRB * CycD_a) + k_61 * pRB_p) - theta_pRB * pRB; E2F1' = (k_p + k_2 * ((power(a, 2) + power(E2F1, 2)) / (power(Km_2, 2) + power(E2F1, 2))) * (J_12 / (J_12 + pRB)) * (J_62 / (J_62 + pRB_p))) - theta_E2F1 * E2F1; // Variable initializations: time_ = ; pRB = 0.1; E2F1 = 0.1; pRB_p = ; CycD_a = ; k_1 = 1; Km_1 = 0.5; J_11 = 0.5; J_61 = 5; k_16 = 0.4; k_61 = 0.3; theta_pRB = 0.005; theta_E2F1 = 0.1; k_p = 0.05; k_2 = 1.6; a = 0.04; Km_2 = 4; J_12 = 5; J_62 = 8; end model swat_2004__first_extension(time_, pRB, E2F1, pRB_p, CycD_a, AP_1, CycD_i, pRB_pp, CycE_a, Fm, k_16, k_61, k_67, k_76) // Rate Rules: pRB_p' = (((k_16 * pRB * CycD_a - k_61 * pRB_p) - k_67 * pRB_p * CycE_a) + k_76 * pRB_pp) - theta_pRB_p * pRB_p; CycD_a' = (k_34 * CycD_i * (CycD_a / (Km_4 + CycD_a)) - k_43 * CycD_a) - theta_CycD_a * CycD_a; AP_1' = (Fm + k_25 * E2F1 * (J_15 / (J_15 + pRB)) * (J_65 / (J_65 + pRB_p))) - theta_AP_1 * AP_1; CycD_i' = ((k_3 * AP_1 + k_23 * E2F1 * (J_13 / (J_13 + pRB)) * (J_63 / (J_63 + pRB_p)) + k_43 * CycD_a) - k_34 * CycD_i * (CycD_a / (Km_4 + CycD_a))) - theta_CycD_i * CycD_i; // Variable initializations: time_ = ; pRB = ; E2F1 = ; pRB_p = 0.1; CycD_a = 0.1; AP_1 = 0.1; CycD_i = 0.1; pRB_pp = ; CycE_a = ; Fm = 0.0057; k_3 = 0.05; k_23 = 0.3; J_13 = 0.002; J_63 = 2; k_43 = 0.01; k_34 = 0.04; Km_4 = 0.3; theta_CycD_i = 0.023; theta_CycD_a = 0.03; theta_AP_1 = 0.01; k_25 = 0.9; J_15 = 0.001; J_65 = 6; k_16 = ; k_61 = ; k_67 = 0.7; k_76 = 0.1; theta_pRB_p = 0.06; end model swat_2004__Second_extension(time_, pRB_p, CycE_a, E2F1, pRB, pRB_pp, CycE_i, k_67, k_76) // Rate Rules: CycE_a' = (k_89 * CycE_i * (CycE_a / (Km_9 + CycE_a)) - k_98 * CycE_a) - theta_CycE_a * CycE_a; pRB_pp' = (k_67 * pRB_p * CycE_a - k_76 * pRB_pp) - theta_pRB_pp * pRB_pp; CycE_i' = ((k_28 * E2F1 * (J_18 / (J_18 + pRB)) * (J_68 / (J_68 + pRB_p)) + k_98 * CycE_a) - k_89 * CycE_i * (CycE_a / (Km_9 + CycE_a))) - theta_CycE_i * CycE_i; // Variable initializations: time_ = ; pRB_p = ; CycE_a = 0.1; E2F1 = ; pRB = ; pRB_pp = 0.1; CycE_i = 0.1; k_67 = ; k_76 = ; theta_pRB_pp = 0.04; k_28 = 0.06; J_18 = 0.6; J_68 = 7; k_98 = 0.01; k_89 = 0.07; Km_9 = 0.005; theta_CycE_i = 0.06; theta_CycE_a = 0.05; end model *swat_2004____main() // Sub-modules, and any changes to those submodules: environment1: swat_2004__environment1(time_); Core_Module: swat_2004__Core_Module(time_, pRB, E2F1, pRB_p, CycD_a, k_16, k_61); first_extension: swat_2004__first_extension(time_, pRB, E2F1, pRB_p, CycD_a, AP_1, CycD_i, pRB_pp, CycE_a, Fm, k_16, k_61, k_67, k_76); Second_extension: swat_2004__Second_extension(time_, pRB_p, CycE_a, E2F1, pRB, pRB_pp, CycE_i, k_67, k_76); end