//Created by libAntimony v2.4 model soboleva_model_2005_version01__environment(time_) // Variable initializations: time_ = ; end model soboleva_model_2005_version01__R(R, f1, time_) // Assignment Rules: logR := log(R); // Rate Rules: R' = f1 - delta_R * R; // Variable initializations: R = 1.0; delta_R = 1.905; f1 = ; time_ = ; end model soboleva_model_2005_version01__U(U, alpha0, alpha1, d01, d12, B1, B2, P, f2, time_) // Assignment Rules: logU := log(U); // Rate Rules: U' = (f2 + d01 * B1 + d12 * B2) - (delta_U * U + alpha0 * P * U + alpha1 * B1 * U); // Variable initializations: U = 1.0; delta_U = 7.5; alpha0 = ; alpha1 = ; d01 = ; d12 = ; B1 = ; B2 = ; P = ; f2 = ; time_ = ; end model soboleva_model_2005_version01__B1(B1, alpha0, alpha1, d01, d12, U, B2, P, time_) // Assignment Rules: logB1 := log(B1); // Rate Rules: B1' = (alpha0 * P * U + d12 * B2) - (d01 * B1 + delta_b1 * B1 + alpha1 * B1 * U); // Variable initializations: B1 = 1.0; delta_b1 = 7.5; alpha0 = ; alpha1 = ; d01 = ; d12 = ; U = ; B2 = ; P = ; time_ = ; end model soboleva_model_2005_version01__B2(B2, alpha1, d12, U, B1, time_) // Assignment Rules: logB2 := log(B2); // Rate Rules: B2' = alpha1 * B1 * U - (d12 * B2 + delta_b2 * B2); // Variable initializations: B2 = 1.0; delta_b2 = 50.0; alpha1 = ; d12 = ; U = ; B1 = ; time_ = ; end model soboleva_model_2005_version01__P(P, P0, time_) // Assignment Rules: P := P0 + PE * q; q := piecewise( 0.0 , time_ < ti , 1.0 - exp(- G * (time_ - ti)) , ( time_ < T + ti) && (geq(time_, ti) ), exp(G * T - 1.0) * exp(G * (time_ - ti)) ); // Variable initializations: PE = 20.55; G = 35.6; ti = 1.0; T = 9.0; P0 = ; time_ = ; end model soboleva_model_2005_version01__f1(f1, x, B2) // Assignment Rules: f1 := a0 * (1.0 - B2 / (b + B2)) + (x * (1.0 - exp(- x / a)) * B2) / (b0 + B2); // Variable initializations: a = 3.58E5; a0 = 2.33E4; b = 100.0; b0 = 263.0; x = ; B2 = ; end model soboleva_model_2005_version01__x(x, time_, P, P0) // Rate Rules: x' = (s * (P - P0)) / (P + P1) - delta_x * x; // Variable initializations: x = 0.0; P1 = 20.0; s = 3.71E5; delta_x = 0.207; time_ = ; P = ; P0 = ; end model soboleva_model_2005_version01__f2(f2, R) // Assignment Rules: f2 := (kappa * R) / (kappa0 + R); // Variable initializations: kappa = 2.4E5; kappa0 = 6.55E3; R = ; end model soboleva_model_2005_version01__model_parameters(alpha0, alpha1, d01, d12, P0) // Variable initializations: alpha0 = 4.15; alpha1 = 3.02; d01 = 13.0; d12 = 4.72E3; P0 = 0.25; end model *soboleva_model_2005_version01____main() // Sub-modules, and any changes to those submodules: environment: soboleva_model_2005_version01__environment(time_); R: soboleva_model_2005_version01__R(R0, f10, time_); U: soboleva_model_2005_version01__U(U0, alpha0, alpha1, d01, d12, B10, B20, P0, f20, time_); B1: soboleva_model_2005_version01__B1(B10, alpha0, alpha1, d01, d12, U0, B20, P0, time_); B2: soboleva_model_2005_version01__B2(B20, alpha1, d12, U0, B10, time_); P: soboleva_model_2005_version01__P(P0, P00, time_); f1: soboleva_model_2005_version01__f1(f10, x0, B20); x: soboleva_model_2005_version01__x(x0, time_, P0, P00); f2: soboleva_model_2005_version01__f2(f20, R0); model_parameters: soboleva_model_2005_version01__model_parameters(alpha0, alpha1, d01, d12, P00); end