-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathEGFR_no_Shp.ant
329 lines (276 loc) · 12.3 KB
/
EGFR_no_Shp.ant
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
// Edited by HSW 5-14-24
// --------------------------------------------------------------------
// Created by libAntimony v2.8.1
model EGFR_combined()
// Compartments:
compartment cell, extra, intra;
// --------------------------------------------------------------------
//Species
// EGFR
species Lig in extra
species Egfr in cell
species iEgfr in intra
species Egfr_Lig in cell
species iEgfr_Lig in intra
species pEgfr_Lig in cell
species ipEgfr_Lig in intra
// Shp2
species Shp2 in cell
species pShp2 in cell
species aR_Shp2 in cell
species aR_pShp2 in cell
// GAREM
species Garem in cell
species Grb2 in cell
species Shp2 in cell
species Grb2_Garem in cell
species Grb2_pGarem in cell
species aR_Grb2_Garem in cell
species aR_Grb2_pGarem in cell
species Grb2_pGarem_pShp2 in cell
// SHC
species Shc1 in cell
species aR_Shc1 in cell
species aR_pShc1 in cell
species pShc1 in cell
// GAB
species Gab1 in cell
species Grb2_Gab1 in cell
species Grb2_pGab1 in cell
species aR_Grb2_Gab1 in cell
species aR_Grb2_pGab1 in cell
species Grb2_pGab1_pShp2 in cell
// SOS
species SOS1 in cell
species Grb2_SOS1 in cell
//
species Ad1_Grb2_SOS1 in cell
species iAd1_Grb2_SOS1 in cell
species iGrb2_pGarem_pShp2 in cell
//
species Ad2_Grb2_SOS1 in cell
species iAd2_Grb2_SOS1 in cell
species iGrb2_pGab1_pShp2 in cell
//
species Ad3_Grb2_SOS1 in cell
species iAd3_Grb2_SOS1 in cell
species ipShc1 in cell
// RAS
species Ras in cell
species aRas1 in cell
species aRas2 in cell
species aRas3 in cell
species Rasa1 in cell
// Total Species
species pEGFRtot in cell
species aRtot in cell
species pGAREMtot in cell
species pGAB1tot in cell
species pSHC1tot in cell
species aRAStot in cell
species pShp2tot in cell
species total_SOS in cell
// --------------------------------------------------------------------
// Reactions:
//EGFR
J1: => Egfr; Vr; //New Receptors
J2: Egfr => iEgfr; kt * Egfr; //Receptor constitutive internalization
J3: iEgfr => ; kx * iEgfr; //Internalized receptor Degradation
J4: Egfr -> Egfr_Lig; kf * Egfr * Lig - kr * Egfr_Lig; //Ligand_binding
J5: Egfr_Lig => iEgfr_Lig; ke * Egfr_Lig; //Receptor ligand endocytosis
J6: iEgfr_Lig => ; kh * iEgfr_Lig; //Internalized receptor Complex Degradation
J7: Egfr_Lig -> pEgfr_Lig; kp1 * Egfr_Lig - k_p1 * pEgfr_Lig; //Receptor Phosphorylation
J8: pEgfr_Lig => ipEgfr_Lig; ke * pEgfr_Lig; //Phosphorylated Receptor Endocytosis
J9: ipEgfr_Lig => ; kh * ipEgfr_Lig; //Phosphorylated internalized Receptor Complex Degradation
//Shp2
J10: pEgfr_Lig + Shp2 -> aR_Shp2; k7 * pEgfr_Lig * Shp2 - k_7 * aR_Shp2; //Shp2 receptor association
J11: aR_Shp2 => aR_pShp2; kp5 * aR_Shp2; //Shp2 phosphorylation
J12: aR_pShp2 => pEgfr_Lig + pShp2; k_7 * aR_pShp2; //EGFR-Shp2 dissociation
J13: pShp2 => Shp2; k_p5 * pShp2; //pShp2 dephosphorylation
//GAREM
J14: Garem + Grb2 -> Grb2_Garem; k1 * Garem * Grb2 - k_1 * Grb2_Garem; //Grb2_Garem complex formation
J15: pEgfr_Lig + Grb2_Garem -> aR_Grb2_Garem; k2 * pEgfr_Lig * Grb2_Garem - k_2 * aR_Grb2_Garem; //Grb2_Garem receptor association
J16: aR_Grb2_Garem => aR_Grb2_pGarem; kp2 * aR_Grb2_Garem; //Grb2_Garem phosphorylation
J17: aR_Grb2_pGarem => pEgfr_Lig + Grb2_pGarem; k_2 * aR_Grb2_pGarem; //Grb2_pGarem dissociation
J18: Grb2_pGarem => Grb2_Garem; k_p2 * Grb2_pGarem; //Grb2_pGarem dephosphorylation
J19: Grb2_pGarem + pShp2 -> Grb2_pGarem_pShp2; k3 * Grb2_pGarem * pShp2 - k_3 * Grb2_pGarem_pShp2; //Grb2_pGarem_pShp2 complex formation
//GAB
J20: Gab1 + Grb2 -> Grb2_Gab1; k4 * Gab1 * Grb2 - k_4 * Grb2_Gab1; //Grb2_Gab1 complex formation
J21: pEgfr_Lig + Grb2_Gab1 -> aR_Grb2_Gab1; k5 * pEgfr_Lig * Grb2_Gab1 - k_5 * aR_Grb2_Gab1; //Grb2_Gab1 receptor association
J22: aR_Grb2_Gab1 => aR_Grb2_pGab1; kp3 * aR_Grb2_Gab1; //Grb2_Gab1 phosphorylation
J23: aR_Grb2_pGab1 => pEgfr_Lig + Grb2_pGab1; k_5 * aR_Grb2_pGab1; //Grb2_pGab1 dissociation
J24: Grb2_pGab1 => Grb2_Gab1; k_p3 * Grb2_pGab1; //Grb2_pGab1 dephosphorylation
J25: Grb2_pGab1 + pShp2 -> Grb2_pGab1_pShp2; k3 * Grb2_pGab1 * pShp2 - k_3 * Grb2_pGab1_pShp2; //Grb2_pGab1_pShp2 complex formation
//SHC
J26: pEgfr_Lig + Shc1 -> aR_Shc1; k6 * pEgfr_Lig * Shc1 - k_6 * aR_Shc1; //SHC1 -pEGFR association
J27: aR_Shc1 => aR_pShc1; kp4 * aR_Shc1; //SHC1 phosphorylation
J28: aR_pShc1 => pEgfr_Lig + pShc1; k_6 * aR_pShc1; //pSHC1 dissociation
J29: pShc1 => Shc1; k_p4 * pShc1; //pSHC1 dephosphorylation
//SOS
J30: Grb2 + SOS1 -> Grb2_SOS1; k8 * Grb2 * SOS1 - k_8 * Grb2_SOS1; //Grb2_SOS1 binding
J31: Grb2_pGarem_pShp2 + Grb2_SOS1 -> Ad1_Grb2_SOS1; k11 * Grb2_pGarem_pShp2 * Grb2_SOS1 - k_11 * Ad1_Grb2_SOS1; //Grb2_SOS1 Garem binding
J32: Ad1_Grb2_SOS1 => iAd1_Grb2_SOS1; ke * Ad1_Grb2_SOS1; //Inhibition by endocytosis
J33: iAd1_Grb2_SOS1 => Grb2_SOS1 + iGrb2_pGarem_pShp2; kr1a * iAd1_Grb2_SOS1; //Recycling of the internalized complexes - first step
J34: iGrb2_pGarem_pShp2 => Grb2_pGarem + pShp2; kr2a * iGrb2_pGarem_pShp2; //Recycling of the internalized complexes - second step
J35: Grb2_pGab1_pShp2 + Grb2_SOS1 -> Ad2_Grb2_SOS1; k12 * Grb2_pGab1_pShp2 * Grb2_SOS1 - k_12 * Ad2_Grb2_SOS1; //Grb2_SOS1 Gab1 binding
J36: Ad2_Grb2_SOS1 => iAd2_Grb2_SOS1; ke * Ad2_Grb2_SOS1; //Inhibition by endocytosis
J37: iAd2_Grb2_SOS1 => Grb2_SOS1 + iGrb2_pGab1_pShp2; kr1b * iAd2_Grb2_SOS1; //Recycling of the internalized complexes - first step
J38: iGrb2_pGab1_pShp2 => Grb2_pGab1 + pShp2; kr2b * iGrb2_pGab1_pShp2; //Recycling of the internalized complexes - second step
J39: pShc1 + Grb2_SOS1 -> Ad3_Grb2_SOS1; k13 * pShc1 * Grb2_SOS1 - k_13 * Ad3_Grb2_SOS1; //Grb2_SOS1 Shc1 binding
J40: Ad3_Grb2_SOS1 => iAd3_Grb2_SOS1; ke * Ad3_Grb2_SOS1; //Inhibition by endocytosis
J41: iAd3_Grb2_SOS1 => Grb2_SOS1 + ipShc1; kr1c * iAd3_Grb2_SOS1; //Recycling of the internalized complexes - first step
J42: ipShc1 => pShc1; kr2c * ipShc1; //Recycling of the internalized complexes - second step
//RAS
J44: Ad1_Grb2_SOS1 + Ras => Ad1_Grb2_SOS1 + aRas1; k14 * Ad1_Grb2_SOS1 * Ras; //Activation by Grb2_SOS1 Garem
J45: Ad2_Grb2_SOS1 + Ras => Ad2_Grb2_SOS1 + aRas2; k14 * Ad2_Grb2_SOS1 * Ras; //Activation by Grb2_SOS1 Gab1
J46: Ad3_Grb2_SOS1 + Ras => Ad3_Grb2_SOS1 + aRas3; k14 * Ad3_Grb2_SOS1 * Ras; //Activation by Grb2_SOS1 Shc1
J47: aRas1 + Rasa1 => Ras + Rasa1; k15 * aRas1 * Rasa1; //Deactivation of Ras1
J48: aRas2 + Rasa1 => Ras + Rasa1; k15 * aRas2 * Rasa1; //Deactivation of Ras2
J49: aRas3 + Rasa1 => Ras + Rasa1; k15 * aRas3 * Rasa1; //Deactivation of Ras3
// Reporter species
pEGFRtot := pEgfr_Lig + ipEgfr_Lig + aR_Grb2_Garem + aR_Grb2_pGarem + aR_Grb2_Gab1 + aR_Grb2_pGab1 + aR_Shp2 + aR_pShp2 + aR_Shc1 + aR_pShc1
pShp2tot := aR_pShp2 + pShp2 + Grb2_pGarem_pShp2 + Ad1_Grb2_SOS1 + iAd1_Grb2_SOS1 + Grb2_pGab1_pShp2 + Ad2_Grb2_SOS1 + iAd2_Grb2_SOS1+ iGrb2_pGarem_pShp2 + iGrb2_pGab1_pShp2
pGAREMtot := aR_Grb2_pGarem + Grb2_pGarem + Grb2_pGarem_pShp2 + Ad1_Grb2_SOS1 +iAd1_Grb2_SOS1 + iGrb2_pGarem_pShp2
pGAB1tot := aR_Grb2_pGab1 + Grb2_pGab1 + Grb2_pGab1_pShp2 + Ad2_Grb2_SOS1 +iAd2_Grb2_SOS1 + iGrb2_pGab1_pShp2
pSHC1tot := aR_pShc1 + pShc1 + Ad3_Grb2_SOS1 +iAd3_Grb2_SOS1 + ipShc1
aRAStot := aRas1 + aRas2 + aRas3
total_SOS := SOS1 + Grb2_SOS1 + Ad1_Grb2_SOS1 + iAd1_Grb2_SOS1 + Ad2_Grb2_SOS1 + iAd2_Grb2_SOS1+ Ad3_Grb2_SOS1 + iAd3_Grb2_SOS1
// --------------------------------------------------------------------
// Initial Values:
// all species has nM units
// EGFR
// fix Lig
var Lig = 0.0;
Lig has nM;
Egfr = 135;
iEgfr = 0;
Egfr_Lig = .97;
iEgfr_Lig = 0.00;
pEgfr_Lig = .085;
ipEgfr_Lig =0;
aRtot := pEgfr_Lig + ipEgfr_Lig
// GAREM
Garem = 1.7
Grb2 = 39
Grb2_Garem = 2
Grb2_pGarem = 0.1
aR_Grb2_Garem = 0.163
aR_Grb2_pGarem = 0.001
// Grb2_pGarem_pShp2 = 0.002
Grb2_pGarem_pShp2 = 0.0
// Shp2
// Shp2 = 27.5;
Shp2 = 0.0;
// pShp2 = 0.27;
pShp2 = 0.0;
// aR_Shp2 = 0.23;
aR_Shp2 = 0.0;
// aR_pShp2 = 0.012;
aR_pShp2 = 0.0;
// SHC
Shc1 = 26.7
aR_Shc1 = 0.046
aR_pShc1 = 0.026
pShc1 = 0.31
/// GAB
Gab1 = 2.1
Grb2_Gab1 = 2.7
aR_Grb2_Gab1 = 0.08
aR_Grb2_pGab1 = 0
// Grb2_pGab1_pShp2 = 0.002
Grb2_pGab1_pShp2 = 0.0
Grb2_pGab1 = 0.08
// SOS
SOS1 = 0.13
Grb2_SOS1 = 3.1
// Ad1_Grb2_SOS1 = 0.007
Ad1_Grb2_SOS1 = 0.0
// iAd1_Grb2_SOS1 = 0.002
iAd1_Grb2_SOS1 = 0.0
// iGrb2_pGarem_pShp2 = 0.002
iGrb2_pGarem_pShp2 = 0.0
// Ad2_Grb2_SOS1 = 0.008
Ad2_Grb2_SOS1 = 0.0
// iAd2_Grb2_SOS1 = 0.002
iAd2_Grb2_SOS1 = 0.0
//iGrb2_pGab1_pShp2 = 0.002
iGrb2_pGab1_pShp2 = 0.0
Ad3_Grb2_SOS1 = 0.31
iAd3_Grb2_SOS1 = 0.08
ipShc1 = 0.08
Ras = 98.0
aRas1 = 0.08
aRas2 = 0.08
aRas3 = 0.67
Rasa1 = 98.0
// --------------------------------------------------------------------
// Compartment initializations:
cell = 1;
extra = 5000;
intra = 1;
// --------------------------------------------------------------------
// Parameters:
Vr = 0.176; // Rate of new EGFR appearing at cell surface
kt = 0.0013; // Turnover rate constant for empty EGFR
kx = 1.00E+00; // Internalized empty receptor degradation rate constant
kf = 1.59E-02; // Forward rate constant for Lig-EGFR binding
kr = 6.85E-03; // Reverse rate constant for Lig-EGFR dissociation
ke = 3.70E-03; // Endocytic rate constant for Lig-EGFR
kh = 5.60E-04; // Degradation rate constant for internalized Lig-EGFR
kp1 = 0.004; // Phosphorylation rate constant for active EGFR
kp2 = 0.0125; // Phosphorylation rate constant for EGFR substrate Garem
kp3 = 0.0072; // Phosphorylation rate constant for EGFR substrate Gab1
kp4 = 0.06; // Phosphorylation rate constant for EGFR substrate Shc1
kp5 = 0.0538; // Phosphorylation rate constant for EGFR substrate PTPN11
k_p1 = 7.00E-02; // Dephosphorylation rate constant for PY EGFR
k_p2 = 0.02; // Dephosphorylation rate constant for PY Garem
k_p3 = 7.20E-03; // Dephosphorylation rate constant for PY Gab1
k_p4 = 0.09; // Dephosphorylation rate constant for PY Shc1
k_p5 = 0.0461; // Dephosphorylation rate constant for PY PTPN11
k1 = 0.0005; // Forward rate constant for Garem-Grb2 binding
k_1 = 0.016; // Reverse rate constant for Garem-Grb2 dissociation
k2 = 1.7; // Forward rate constant for Garem-Grb2 binding to phosphorylated EGFR
k_2 = 1.8; // Reverse rate constant for Garem-Grb2 dissociating from phosphorylated EGFR
k3 = 0.0214; // Forward rate constant for Shp2 binding to doubly phosphorylated Garem or Gab1
k_3 = 1.5; // Reverse rate constant for Shp2 dissociating from doubly phosphorylated Garem or Gab1
k4 = 0.001; // Forward rate constant for Gab1-Grb2 binding
k_4 = .03; // Reverse rate constant for Gab1-Grb2 dissociation
k5 = 0.5; // Forward rate constant for Gab1-Grb2 binding to phosphorylated EGFR
k_5 = 1.46; // Reverse rate constant for Gab1-Grb2 dissociating from phosphorylated EGFR
k6 = 0.23; // Forward rate constant for Shc1 binding to phosphorylated EGFR
k_6 = 1.07; // Reverse rate constant for Shc1 dissociating from phosphorylated EGFR
k7 = 0.11; // Forward rate constant for EGFR-Shp2 binding
k_7 = 1.07; // Reverse rate constant for EGFR-Shp2 dissociation
k8 = 9.00E-05; // Forward rate constant for Grb2 - SOS1 binding
k_8 = 1.40E-04; // Reverse rate constant for Grb2_SOS1 dissociation
k11 = .09; // Forward rate constant for Grb2_SOS1 binding to Garem-Shp2 complex
k_11 = 0.011; // Reverse rate constant for Grb2_SOS1 dissociation from Garem-Shp2 complex
k12 = 0.105; // Forward rate constant for Grb2_SOS1 binding to Gab1-Shp2 complex
k_12 = 0.011; // Reverse rate constant for Grb2_SOS1 dissociation from Gab1-Shp2 complex
k13 = 0.02; // Forward rate constant for Grb2_SOS1 binding to pEGFR_pShc1 complex
k_13 = 0.043; // Reverse rate constant for Grb2_SOS1 dissociation from pEGFR_pShc1 complex
k14 = 0.18; // Activation of Ras by SOS-adaptor complexes
k15 = 0.0164; // Deactivation rate constant of active Ras by RASA1
kr1a = 1.50E-02; // Recycling of the internalized Garem complexes - first step
kr2a = 1.50E-02; // Recycling of the internalized Garem complexes - second step
kr1b = 1.50E-02; // Recycling of the internalized Gab1 complexes - first step
kr2b = 1.50E-02; // Recycling of the internalized Gab1 complexes - second step
kr1c = 1.50E-02; // Recycling of the internalized Shc1 complexes - first step
kr2c = 1.50E-02; // Recycling of the internalized Shc1 complexes - second step
// --------------------------------------------------------------------
// Unit definitions:
unit substance = 1e-9 mole;
unit time_unit = second;
unit nM = 1e-9 mole / litre;
unit per_nM = litre / 1e-9 mole;
unit per_sec = 1 / second;
unit nM_per_s = 1e-9 mole / (litre * second);
unit per_nMs = litre / (1e-9 mole * second);
// Display Names:
substance is "nmole";
time_unit is "s";
per_sec is "s-1";
nM_per_s is "nM/s";
per_nMs is "nM-1 s-1";
extra is "extracellular space";
end