1000 ug/ml is max

2 - ku * Plocal - ke * Pex rM - dM * M T_macro + T_senescence - T_local * R * kon_t + koff_t * RT - T * dT rmit - dmit * MIT rD - dD * D ke * Pex ku * Plocal

time>95 Pex_start

(time-96)/12

if(glucu_sulf_total>0 or oxidation_total > 0,(glucu_sulf_total/(glucu_sulf_total+oxidation_total))*100,0) if(glucu_sulf_total>0 or oxidation_total > 0,(oxidation_total/(glucu_sulf_total+oxidation_total))*100,0) if(efflux_total>0 or uptake_total > 0,(efflux_total/(efflux_total+uptake_total))*100,0) if(efflux_total>0 or uptake_total > 0,(uptake_total/(efflux_total+uptake_total))*100,0)

if(macrophages_total > 100,switch+1,switch+0) if(M_pp > alltime_M_pp_max, M_pp, alltime_M_pp_max) 1

1

ku * Plocal - Pin * kP450 * P450_cell - Pin * kg - ks * S * Pin + N * kN Pin * kP450 * P450_cell - N * kGSH * G_cell - N * kpsh - N * kN T_local * R * kon_t - koff_t * RT N * kpsh - C * dc ks_p53 + k_dp * p53p - k_p * p53 * C - kd_p53 * p53 - kd_p53_mdm2 * MDM2 * p53 k_p * p53 * C - k_dp * p53p - kd_p53p * p53p - kd_p53p_mdm2 * MDM2 * p53p ks_mdm2 + (ks_mdm2_p53p * p53p^4) / (Km_mdm2^4 + p53p^4) - kd_mdm2 * MDM2 ks_p21 + (ks_p21_p53p * p53p^4) / (Km_p21^4 + p53p^4) - kd_p21 * p21 + RT * k_stim bGSH - dGSH * G_cell - N * kGSH * G_cell bk + k_mitogen * MIT_local - k_inhibitor * p21 * K - dk * K - T_local * R * kon_t + koff_t * RT N * kGSH * G_cell Pin * ks * S - ks * S * Pin + bs - ds * S kg * Pin if(neigh_hepmacro > 0, nt_hepmacro + 1, 0) divisions + 1 7 rand_uni(0,1)<p_proliferation && K > 4.6863425 C > 4 && rand_uni(0,1) < p_necrosis 4

senescence == 1 && neigh_hepmacro > 0 && nt_hepmacro > 60 && rand_uni(0,1) < 0.45 K < 4.6863375 or senescence > 1 6 0.01 senescence + 1 0.0 sqrt(((cell.center.x-dimension)^2)+((cell.center.y-dimension)^2)) (distance/radius) (d_p450+(1/(1+exp(b*(distance_normal-a)))))/c_p450 (d_gsh+(1/(1+exp(-b*(distance_normal-a)))))/c_gsh N + C + NG dGSH*G_dist

macro_neigh > 0 && macro_neigh_time > limit_neigh && rand_uni(0,1) > 0.1

if(macro_neigh > 0, macro_neigh_time+1, macro_neigh_time+0) lifetime_necrotic + 1

if(lifetime_necrotic > 1,1,0) periportal_area if(M_pp>0.0001 && M_pp > alltime_M_pp_max*0.9 , p_macrophage,0)

((lifetime_macrophage > 200 and M_mean < 1 * 10^-5) && rand_uni(0,1) < 0.005) lifetime_macrophage + 1 if(nec_neigh > 0, nec_neigh_time+1, nec_neigh_time+0)

if(nec_neigh > 0, T_macro_on, 0) if(pv_neigh > 0, 600, if(nec_neigh > 0, 2,100))

if(NR_necrotic > 0, NR_time +5, 0)

if(D_local > 0.5, rM_true, 0) 2 3