Skip to content
Snippets Groups Projects
Commit 73496c54 authored by herrTilda's avatar herrTilda
Browse files

cleaned M1

parent 14196bba
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 32 additions and 37 deletions
scripts/Models/M1.txt
+ 32
− 37
View file @ 73496c54
......@@ -31,12 +31,11 @@ d/dt(GLUT4) = -v7f+v7b %%24
d/dt(S6K) = v9b1-v9f1 %%26
d/dt(S6Kp) = v9f1-v9b1 %%27
d/dt(S6) = v9b2-v9f2 %%28
d/dt(S6p) = v9f2-v9b2 %%29
d/dt(S6p) = v9f2-v9b2 %%29
d/dt(G_p) = EGP+Ra-E-U_ii-k_1*G_p+k_2*G_t
d/dt(G_t) = (-U_id)+k_1*G_p-k_2*G_t
d/dt(I_p) = (-m_2*I_p)-m_4*I_p+m_1*I_l
d/dt(G_t) = (-U_id)+k_1*G_p-k_2*G_t
d/dt(I_l) = (-m_1*I_l)-m_3*I_l+m_2*I_p+S
d/dt(I_p) = (-m_2*I_p)-m_4*I_p+m_1*I_l
d/dt(Q_sto1) = -k_gri*Q_sto1
d/dt(Q_sto2) = (-k_empt*Q_sto2)+k_gri*Q_sto1
d/dt(Q_gut) = (-k_abs*Q_gut)+k_empt*Q_sto2
......@@ -74,8 +73,7 @@ GLUT4(0) = 100
S6K(0) = 100
S6Kp(0) = 0
S6(0) = 100
S6p(0) = 0
S6p(0) = 0
G_p(0) = 178
G_t(0) = 130
I_l(0) = 4.5
......@@ -137,7 +135,7 @@ k9b2 = 30.9967 %%42
km9 = 5872.68 %%43
n9 = 0.985466 %%44
kbf = 1e+06 %%45
nC = 2.1e-06 %%46
scaleModel = 2.1e-06 %%46
V_G = 1.88 %%47
k_1 = 0.065 %%48
......@@ -166,38 +164,33 @@ k_p3 = 0.009 %%70
k_p4 = 0.0618 %%71
k_i = 0.0079 %%72
V_m0 = 2.5 %%73
V_mX = 0.047 %%74
K_m0 = 225.59 %%75
V_f0 = 2.5 %%76
V_fX = 0.047 %%77
K_f0 = 225.59 %%78
p_2U = 0.0331 %%79
part = 0.2 %%80
K = 2.3 %%81
alpha = 0.05 %%82
beta = 0.11 %%83
gamma = 0.5 %%84
k_e1 = 0.0005 %%85
k_e2 = 339 %%86
D = 78000 %%87
V_m0 = 2.5 %%74
V_mX = 0.047 %%75
K_m0 = 225.59 %%76
V_f0 = 2.5 %%77
V_fX = 0.047 %%78
K_f0 = 225.59 %%79
p_2U = 0.0331 %%80
part = 0.2 %%81
K = 2.3 %%82
alpha = 0.05 %%83
beta = 0.11 %%84
gamma = 0.5 %%85
k_e1 = 0.0005 %%86
k_e2 = 339 %%87
D = 78000 %%88
K_l0 = 225.59 %%88
V_l0 = 2.5 %%89
V_lX = 0.047 %%90
pf = 34
U_ii = 0.8
part_m = 0.5
part_f = 0.5
part_l = 1.7
K_l0 = 225.59
V_l0 = 2.5
V_lX = 0.047
********** MODEL VARIABLES
aa = 5/2/(1-b)/D
cc = 5/2/d/D
EGP = k_p1-k_p2*G_p-k_p3*I_d-k_p4*I_po
V_mmax = (1-part)*(V_m0+V_mX*INS)
V_fmax = part*(V_f0+V_fX*INS)
V_fmax = part*(V_f0+V_fX*INS)
V_lmax = (1-part)*(V_l0+V_lX*G_t)
Ee = k_e1*(G_p - k_e2)
S = gamma*I_po
I = I_p/V_I
......@@ -207,11 +200,13 @@ m_3 = HE*m_1/(1-HE)
Q_sto = Q_sto1+Q_sto2
Ra = f*k_abs*Q_gut/BW
k_empt = k_min+(k_max-k_min)/2*(tanh(aa*(Q_sto-b*D))-tanh(cc*(Q_sto-d*D))+2)
U_idm = part_m*V_mmax*G_t/(K_m0+G_t)
U_idf = part_f*nC*(k8*GLUT4m*G_t/pf + glut1*G_t/pf + kbf*(INS+5)*1e-3)
U_idm = 0.5*V_mmax*G_t/(K_m0+G_t)
U_idf = 0.5*scaleModel*(k8*GLUT4m*G_t/34 + glut1*G_t/34 + kbf*(INS+5)*1e-3)
U = U_ii+U_id
S_po = Y+K*(EGP+Ra-E-U_ii-k_1*G_p+k_2*G_t)/V_G+S_b
U_idl = part_l*U_idm
U_id = U_idm+U_idf+U_idl
U_liver = 45*U_idm/27
U_ii = 23*U_liver/45
U_id = U_idm+U_idf+U_liver
********** MODEL REACTIONS
v1a = IR*k1a*(INS+5)*1e-3
......@@ -259,4 +254,4 @@ glucoseuptake = k8*GLUT4m*(G_t*5/170) + glut1*(G_t*5/170) + kbf*((INS+5)*1e-3)
********** MODEL EVENTS
event1 = gt(G_p,339),E,Ee
event2 = lt(G_p,339),E,0
********** MODEL MATLAB FUNCTIONS
\ No newline at end of file
********** MODEL MATLAB FUNCTIONS
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment