EGFR_IGF1R_model.m
2.85 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
function dx= EGFR_IGF1R_model(t,x,p,u,xT)
% Parameters
gamma_EGFR=p(1);
gamma_IGF1R=p(2);
kd_P90Rsk=p(3);
kd_PIK3_active=p(4);
kd_Akt=p(5);
k_SOS_EGFR=p(6);
KM_SOS_EGFR=p(7);
k_SOS_IGF1R=p(8);
KM_SOS_IGF1R=p(9);
k_DSOS_p90Rsk=p(10);
KM_DSOS_p90Rsk=p(11);
k_Ras_SOS=p(12);
KM_Ras_SOS=p(13);
k_Ras_RasGab=p(14);
KM_Ras_RasGab=p(15);
k_Raf_Ras=p(16);
KM_Raf_Ras=p(17);
k_Raf_RafPP=p(18);
KM_Raf_RafPP=p(19);
k_Raf_Akt=p(20);
KM_Raf_Akt=p(21);
k_Mek_Raf=p(22);
KM_Mek_Raf=p(23);
k_MEK_PP2A=p(24);
KM_MEK_PP2A=p(25);
k_ERK_MEK=p(26);
KM_ERK_MEK=p(27);
k_Erk_PP2A=p(28);
KM_Erk_PP2A=p(29);
k_p90Rsk_Erk=p(30);
KM_p90Rsk_Erk=p(31);
k_PIK3_Ras=p(32);
KM_PIK3_Ras=p(33);
k_PIK3_IGF1R=p(34);
KM_PIK3_IGF1R=p(35);
k_PIK3_EGFR=p(36);
KM_PIK3_EGFR=p(37);
k_Akt_PIK3=p(38);
KM_Akt_PIK3=p(39);
% Input values
RafPP=u(1);
PP2A=u(2);
RasGapActive=u(3);
% State Variables
EGFR_active=x(1);
IGF1R_active=x(2);
SOS=x(3);
%DSOS=x(4);
Ras_active=x(4);
Raf_active=x(5);
Mek_active=x(6);
Erk_active=x(7);
P90Rsk_active=x(8);
PIK3_active=x(9);
Akt_active=x(10);
% Total concentrations
DSOS=xT(1);
Ras_tot=xT(2);
Raf_tot=xT(3);
Mek_tot=xT(4);
Erk_tot=xT(5);
P90Rsk_tot=xT(6);
PIK3_tot=xT(7);
Akt_tot=xT(8);
% ODEs
dEGFR_active=-gamma_EGFR*EGFR_active;
dIGF1R_active=-gamma_IGF1R*IGF1R_active;
dSOS=k_SOS_EGFR*EGFR_active*(DSOS-SOS)/(KM_SOS_EGFR+(DSOS-SOS))+k_SOS_IGF1R*IGF1R_active*(DSOS-SOS)/(KM_SOS_IGF1R+(DSOS-SOS))-k_DSOS_p90Rsk*P90Rsk_active*SOS/(KM_DSOS_p90Rsk+SOS);
%dDSOS=-k_SOS_EGFR*EGFR_active*DSOS/(KM_SOS_EGFR+DSOS)-k_SOS_IGF1R*IGF1R_active*DSOS/(KM_SOS_IGF1R+DSOS) +k_DSOS_p90Rsk*P90Rsk_active*SOS/(KM_DSOS_p90Rsk+SOS);
dRas_active=k_Ras_SOS*SOS*(Ras_tot-Ras_active)/(KM_Ras_SOS+(Ras_tot-Ras_active))-k_Ras_RasGab*RasGapActive*Ras_active/(KM_Ras_RasGab+Ras_active);
dRaf_active=k_Raf_Ras*Ras_active*(Raf_tot-Raf_active)/(KM_Raf_Ras+(Raf_tot-Raf_active))-k_Raf_RafPP*RafPP*Raf_active/(KM_Raf_RafPP+Raf_active)-k_Raf_Akt*Akt_active*Raf_active/(KM_Raf_Akt+Raf_active);
dMek_active=k_Mek_Raf*Raf_active*(Mek_tot-Mek_active)/(KM_Mek_Raf+(Mek_tot-Mek_active))-k_MEK_PP2A*PP2A*Mek_active/(KM_MEK_PP2A+Mek_active);
dErk_active=k_ERK_MEK*Mek_active*(Erk_tot-Erk_active)/(KM_ERK_MEK+(Erk_tot-Erk_active))-k_Erk_PP2A*PP2A*Erk_active/(KM_Erk_PP2A+Erk_active);
dP90Rsk_active=k_p90Rsk_Erk*Erk_active*(P90Rsk_tot-P90Rsk_active)/(KM_p90Rsk_Erk+(P90Rsk_tot-P90Rsk_active))-kd_P90Rsk*P90Rsk_active;
dPIK3_active=k_PIK3_IGF1R*IGF1R_active*(PIK3_tot-PIK3_active)/(KM_PIK3_IGF1R+(PIK3_tot-PIK3_active))+k_PIK3_EGFR*EGFR_active*(PIK3_tot-PIK3_active)/(KM_PIK3_EGFR+(PIK3_tot-PIK3_active))+k_PIK3_Ras*Ras_active*(PIK3_tot-PIK3_active)/(KM_PIK3_Ras+(PIK3_tot-PIK3_active))-kd_PIK3_active*PIK3_active;
dAkt_active=k_Akt_PIK3*PIK3_active*(Akt_tot-Akt_active)/(KM_Akt_PIK3+(Akt_tot-Akt_active))-kd_Akt*Akt_active;
dx=[dEGFR_active;dIGF1R_active;dSOS;dRas_active;dRaf_active;dMek_active;dErk_active;dP90Rsk_active;dPIK3_active;dAkt_active];