Upload new file

CRA_Updated_Version/CRA_Toolbox_bash/Examples/EGFR_IGF1R_model.m

+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];
\ No newline at end of file