forked from BhallaLab/FindSim
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathsynSynth5.g
More file actions
4375 lines (4369 loc) · 335 KB
/
synSynth5.g
File metadata and controls
4375 lines (4369 loc) · 335 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
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
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
//genesis
// kkit Version 11 flat dumpfile
// Saved on Sun Apr 1 17:17:23 2018
include kkit {argv 1}
FASTDT = 0.001
SIMDT = 0.001
CONTROLDT = 0.1
PLOTDT = 0.1
MAXTIME = 100
TRANSIENT_TIME = 2
VARIABLE_DT_FLAG = 0
DEFAULT_VOL = 1.049e-15
VERSION = 11.0
setfield /file/modpath value ~/scripts/modules
kparms
//genesis
initdump -version 3 -ignoreorphans 1
simobjdump doqcsinfo filename accessname accesstype transcriber developer \
citation species tissue cellcompartment methodology sources \
model_implementation model_validation x y z
simobjdump table input output alloced step_mode stepsize x y z
simobjdump xtree path script namemode sizescale
simobjdump xcoredraw xmin xmax ymin ymax
simobjdump xtext editable
simobjdump xgraph xmin xmax ymin ymax overlay
simobjdump xplot pixflags script fg ysquish do_slope wy
simobjdump group xtree_fg_req xtree_textfg_req plotfield expanded movealone \
link savename file version md5sum mod_save_flag x y z
simobjdump geometry size dim shape outside xtree_fg_req xtree_textfg_req x y \
z
simobjdump kpool DiffConst CoInit Co n nInit mwt nMin vol slave_enable \
geomname xtree_fg_req xtree_textfg_req x y z
simobjdump kreac kf kb notes xtree_fg_req xtree_textfg_req x y z
simobjdump kenz CoComplexInit CoComplex nComplexInit nComplex vol k1 k2 k3 \
keepconc usecomplex notes xtree_fg_req xtree_textfg_req link x y z
simobjdump stim level1 width1 delay1 level2 width2 delay2 baselevel trig_time \
trig_mode notes xtree_fg_req xtree_textfg_req is_running x y z
simobjdump xtab input output alloced step_mode stepsize notes editfunc \
xtree_fg_req xtree_textfg_req baselevel last_x last_y is_running x y z
simobjdump kchan perm gmax Vm is_active use_nernst notes xtree_fg_req \
xtree_textfg_req x y z
simobjdump transport input output alloced step_mode stepsize dt delay clock \
kf xtree_fg_req xtree_textfg_req x y z
simobjdump proto x y z
simobjdump text str
simundump geometry /kinetics/geometry 0 1.0529e-15 3 sphere "" white black 0 \
0 0
simundump text /kinetics/notes 0 ""
call /kinetics/notes LOAD \
""
simundump group /kinetics/PKC 0 0 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 0 0 0
simundump text /kinetics/PKC/notes 0 ""
call /kinetics/PKC/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_Ca 0 0.0 3.7208e-17 3.7208e-17 2.3506e-11 \
2.3506e-11 0 0 6.3175e+05 0 /kinetics/geometry red 0 0 1 0
simundump text /kinetics/PKC/PKC_Ca/notes 0 ""
call /kinetics/PKC/PKC_Ca/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_DAG_AA_p 0 0.0 4.9137e-18 4.9137e-18 \
3.1042e-12 3.1042e-12 0 0 6.3175e+05 0 /kinetics/geometry cyan 0 2 1 0
simundump text /kinetics/PKC/PKC_DAG_AA_p/notes 0 ""
call /kinetics/PKC/PKC_DAG_AA_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_Ca_AA_p 0 0.0 1.7501e-16 1.7501e-16 \
1.1056e-10 1.1056e-10 0 0 6.3175e+05 0 /kinetics/geometry orange 0 4 1 0
simundump text /kinetics/PKC/PKC_Ca_AA_p/notes 0 ""
call /kinetics/PKC/PKC_Ca_AA_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_Ca_memb_p 0 0.0 1.3896e-17 1.3896e-17 \
8.7788e-12 8.7788e-12 0 0 6.3175e+05 0 /kinetics/geometry pink 0 6 1 0
simundump text /kinetics/PKC/PKC_Ca_memb_p/notes 0 ""
call /kinetics/PKC/PKC_Ca_memb_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_DAG_memb_p 0 0.0 9.4352e-21 9.4352e-21 \
5.9607e-15 5.9607e-15 0 0 6.3175e+05 0 /kinetics/geometry yellow 0 8 1 0
simundump text /kinetics/PKC/PKC_DAG_memb_p/notes 0 ""
call /kinetics/PKC/PKC_DAG_memb_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_basal_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 0 10 1 0
simundump text /kinetics/PKC/PKC_basal_p/notes 0 ""
call /kinetics/PKC/PKC_basal_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_AA_p 0 0.0 1.8134e-17 1.8134e-17 1.1456e-11 \
1.1456e-11 0 0 6.3175e+05 0 /kinetics/geometry cyan 0 12 1 0
simundump text /kinetics/PKC/PKC_AA_p/notes 0 ""
call /kinetics/PKC/PKC_AA_p/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_Ca_DAG 0 0.0 8.4632e-23 8.4632e-23 \
5.3466e-17 5.3466e-17 0 0 6.3175e+05 0 /kinetics/geometry white 0 14 1 0
simundump text /kinetics/PKC/PKC_Ca_DAG/notes 0 ""
call /kinetics/PKC/PKC_Ca_DAG/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_DAG 0 0.0 1.161e-16 1.161e-16 7.3348e-11 \
7.3348e-11 0 0 6.3175e+05 0 /kinetics/geometry white 0 16 1 0
simundump text /kinetics/PKC/PKC_DAG/notes 0 "CoInit was .0624"
call /kinetics/PKC/PKC_DAG/notes LOAD \
"CoInit was .0624"
simundump kpool /kinetics/PKC/PKC_DAG_AA 0 0.0 2.5189e-19 2.5189e-19 \
1.5913e-13 1.5913e-13 0 0 6.3175e+05 0 /kinetics/geometry white 0 1 16 0
simundump text /kinetics/PKC/PKC_DAG_AA/notes 0 ""
call /kinetics/PKC/PKC_DAG_AA/notes LOAD \
""
simundump kpool /kinetics/PKC/PKC_cytosolic 0 0.0 1 1 6.3175e+05 6.3175e+05 0 \
0 6.3175e+05 0 /kinetics/geometry white 0 3 16 0
simundump text /kinetics/PKC/PKC_cytosolic/notes 0 \
"Marquez et al J. Immun 149,2560(92) est 1e6/cell for chromaffin cells We will use 1 uM as our initial concen"
call /kinetics/PKC/PKC_cytosolic/notes LOAD \
"Marquez et al J. Immun 149,2560(92) est 1e6/cell for chromaffin cells We will use 1 uM as our initial concen"
simundump kpool /kinetics/PKC/PKC_active 0 0.0 0 2.1196e-16 1.339e-10 0 0 0 \
6.3175e+05 0 /kinetics/geometry red 0 9 16 0
simundump text /kinetics/PKC/PKC_active/notes 0 ""
call /kinetics/PKC/PKC_active/notes LOAD \
""
simundump kenz /kinetics/PKC/PKC_active/PKC_act_raf 0 0 0 0 0 6.3175e+05 \
4.51e-07 16 4 0 0 "" yellow red "" 9 14 0
simundump text /kinetics/PKC/PKC_active/PKC_act_raf/notes 0 \
"Rate consts from Chen et al Biochem 32, 1032 (1993) k3 = k2 = 4 k1 = 9e-5 recalculated gives 1.666e-5, which is not very different. Looks like k3 is rate-limiting in this case: there is a huge amount of craf locked up in the enz complex. Let us assume a 10x higher Km, ie, lower affinity. k1 drops by 10x. Also changed k2 to 4x k3. Lowerd k1 to 1e-6 to balance 10X DAG sensitivity of PKC"
call /kinetics/PKC/PKC_active/PKC_act_raf/notes LOAD \
"Rate consts from Chen et al Biochem 32, 1032 (1993) k3 = k2 = 4 k1 = 9e-5 recalculated gives 1.666e-5, which is not very different. Looks like k3 is rate-limiting in this case: there is a huge amount of craf locked up in the enz complex. Let us assume a 10x higher Km, ie, lower affinity. k1 drops by 10x. Also changed k2 to 4x k3. Lowerd k1 to 1e-6 to balance 10X DAG sensitivity of PKC"
simundump kenz /kinetics/PKC/PKC_active/PKC_inact_GAP 0 0 0 0 0 6.3175e+05 \
9.02e-06 16 4 0 0 "" yellow red "" 9 13 0
simundump text /kinetics/PKC/PKC_active/PKC_inact_GAP/notes 0 \
"Rate consts copied from PCK-act-raf This reaction inactivates GAP. The idea is from the Boguski and McCormick review."
call /kinetics/PKC/PKC_active/PKC_inact_GAP/notes LOAD \
"Rate consts copied from PCK-act-raf This reaction inactivates GAP. The idea is from the Boguski and McCormick review."
simundump kenz /kinetics/PKC/PKC_active/PKC_act_GEF 0 0 0 0 0 6.3175e+05 \
9.02e-06 16 4 0 0 "" yellow red "" 9 12 0
simundump text /kinetics/PKC/PKC_active/PKC_act_GEF/notes 0 \
"Rate consts from PKC-act-raf. This reaction activates GEF. It can lead to at least 2X stim of ras, and a 2X stim of MAPK over and above that obtained via direct phosph of c-raf. Note that it is a push-pull reaction, and there is also a contribution through the phosphorylation and inactivation of GAPs. The original PKC-act-raf rate consts are too fast. We lower K1 by 10 X"
call /kinetics/PKC/PKC_active/PKC_act_GEF/notes LOAD \
"Rate consts from PKC-act-raf. This reaction activates GEF. It can lead to at least 2X stim of ras, and a 2X stim of MAPK over and above that obtained via direct phosph of c-raf. Note that it is a push-pull reaction, and there is also a contribution through the phosphorylation and inactivation of GAPs. The original PKC-act-raf rate consts are too fast. We lower K1 by 10 X"
simundump kenz /kinetics/PKC/PKC_active/PKC_phosph_neurogranin 0 0 0 0 0 \
6.3175e+05 1.5334e-07 2.34 0.58 0 0 "" red red "" 9 11 0
simundump text /kinetics/PKC/PKC_active/PKC_phosph_neurogranin/notes 0 \
"Rates from Huang et al ABB 305:2 570-580 1993"
call /kinetics/PKC/PKC_active/PKC_phosph_neurogranin/notes LOAD \
"Rates from Huang et al ABB 305:2 570-580 1993"
simundump kenz /kinetics/PKC/PKC_active/PKC_phosph_ng_CaM 0 0 0 0 0 \
6.3175e+05 9.2004e-08 1.4 0.35 0 0 "" red red "" 9 10 0
simundump text /kinetics/PKC/PKC_active/PKC_phosph_ng_CaM/notes 0 \
"Rates are 60% those of PKC-phosph-neurogranin. See Huang et al ABB 305:2 570-580 1993"
call /kinetics/PKC/PKC_active/PKC_phosph_ng_CaM/notes LOAD \
"Rates are 60% those of PKC-phosph-neurogranin. See Huang et al ABB 305:2 570-580 1993"
simundump kenz /kinetics/PKC/PKC_active/phosph_AC2 0 0 0 0 0 6.3175e+05 \
9.02e-07 16 4 0 0 "" red red "" 9 9 0
simundump text /kinetics/PKC/PKC_active/phosph_AC2/notes 0 \
"Phorbol esters have little effect on AC1 or on the Gs-stimulation of AC2. So in this model we are only dealing with the increase in basal activation of AC2 induced by PKC k1 = 1.66e-6 k2 = 16 k3 =4"
call /kinetics/PKC/PKC_active/phosph_AC2/notes LOAD \
"Phorbol esters have little effect on AC1 or on the Gs-stimulation of AC2. So in this model we are only dealing with the increase in basal activation of AC2 induced by PKC k1 = 1.66e-6 k2 = 16 k3 =4"
simundump kreac /kinetics/PKC/PKC_act_by_Ca 0 9.4974e-07 0.5 "" white 0 1 3 0
simundump text /kinetics/PKC/PKC_act_by_Ca/notes 0 \
"Need est of rate of assoc of Ca and PKC. Assume it is fast The original parameter-searched kf of 439.4 has been scaled by 1/6e8 to account for change of units to n. Kf now 8.16e-7, kb=.6085 Raised kf to 1e-6 to match Ca curve, kb to .5"
call /kinetics/PKC/PKC_act_by_Ca/notes LOAD \
"Need est of rate of assoc of Ca and PKC. Assume it is fast The original parameter-searched kf of 439.4 has been scaled by 1/6e8 to account for change of units to n. Kf now 8.16e-7, kb=.6085 Raised kf to 1e-6 to match Ca curve, kb to .5"
simundump kreac /kinetics/PKC/PKC_act_by_DAG 0 1.2663e-08 8.6348 "" white 0 3 \
3 0
simundump text /kinetics/PKC/PKC_act_by_DAG/notes 0 \
"Need est of rate. Assume it is fast Obtained from param search kf raised 10 X : see Shinomura et al PNAS 88 5149-5153 1991. kf changed from 3.865e-7 to 2.0e-7 in line with closer analysis of Shinomura data. 26 June 1996: Corrected DAG data: reduce kf 15x from 2e-7 to 1.333e-8"
call /kinetics/PKC/PKC_act_by_DAG/notes LOAD \
"Need est of rate. Assume it is fast Obtained from param search kf raised 10 X : see Shinomura et al PNAS 88 5149-5153 1991. kf changed from 3.865e-7 to 2.0e-7 in line with closer analysis of Shinomura data. 26 June 1996: Corrected DAG data: reduce kf 15x from 2e-7 to 1.333e-8"
simundump kreac /kinetics/PKC/PKC_Ca_to_memb 0 1.2705 3.5026 "" white 0 5 3 0
simundump text /kinetics/PKC/PKC_Ca_to_memb/notes 0 ""
call /kinetics/PKC/PKC_Ca_to_memb/notes LOAD \
""
simundump kreac /kinetics/PKC/PKC_DAG_to_memb 0 1 0.1 "" white 0 7 3 0
simundump text /kinetics/PKC/PKC_DAG_to_memb/notes 0 \
"Raise kb from .087 to 0.1 to match data from Shinomura et al. Lower kf from 1.155 to 1.0 to match data from Shinomura et al."
call /kinetics/PKC/PKC_DAG_to_memb/notes LOAD \
"Raise kb from .087 to 0.1 to match data from Shinomura et al. Lower kf from 1.155 to 1.0 to match data from Shinomura et al."
simundump kreac /kinetics/PKC/PKC_act_by_Ca_AA 0 1.8995e-09 0.1 "" white 0 9 \
3 0
simundump text /kinetics/PKC/PKC_act_by_Ca_AA/notes 0 \
"Schaechter and Benowitz We have to increase Kf for conc scaling Changed kf to 2e-9 on Sept 19, 94. This gives better match."
call /kinetics/PKC/PKC_act_by_Ca_AA/notes LOAD \
"Schaechter and Benowitz We have to increase Kf for conc scaling Changed kf to 2e-9 on Sept 19, 94. This gives better match."
simundump kreac /kinetics/PKC/PKC_act_by_DAG_AA 0 2 0.2 "" white 0 11 3 0
simundump text /kinetics/PKC/PKC_act_by_DAG_AA/notes 0 \
"Assume slowish too. Schaechter and Benowitz"
call /kinetics/PKC/PKC_act_by_DAG_AA/notes LOAD \
"Assume slowish too. Schaechter and Benowitz"
simundump kreac /kinetics/PKC/PKC_basal_act 0 1 50 "" white 0 13 3 0
simundump text /kinetics/PKC/PKC_basal_act/notes 0 \
"Initial basal levels were set by kf = 1, kb = 20. In model, though, the basal levels of PKC activity are higher."
call /kinetics/PKC/PKC_basal_act/notes LOAD \
"Initial basal levels were set by kf = 1, kb = 20. In model, though, the basal levels of PKC activity are higher."
simundump kreac /kinetics/PKC/PKC_act_by_AA 0 1.8995e-10 0.1 "" white 0 15 3 \
0
simundump text /kinetics/PKC/PKC_act_by_AA/notes 0 \
"Raise kf from 1.667e-10 to 2e-10 to get better match to data."
call /kinetics/PKC/PKC_act_by_AA/notes LOAD \
"Raise kf from 1.667e-10 to 2e-10 to get better match to data."
simundump kreac /kinetics/PKC/PKC_n_DAG 0 9.4974e-10 0.1 "" white 0 17 3 0
simundump text /kinetics/PKC/PKC_n_DAG/notes 0 \
"kf raised 10 X based on Shinomura et al PNAS 88 5149-5153 1991 closer analysis of Shinomura et al: kf now 1e-8 (was 1.66e-8). Further tweak. To get sufficient AA synergy, increase kf to 1.5e-8 26 June 1996: Corrected DAG levels: reduce kf by 15x from 1.5e-8 to 1e-9"
call /kinetics/PKC/PKC_n_DAG/notes LOAD \
"kf raised 10 X based on Shinomura et al PNAS 88 5149-5153 1991 closer analysis of Shinomura et al: kf now 1e-8 (was 1.66e-8). Further tweak. To get sufficient AA synergy, increase kf to 1.5e-8 26 June 1996: Corrected DAG levels: reduce kf by 15x from 1.5e-8 to 1e-9"
simundump kreac /kinetics/PKC/PKC_n_DAG_AA 0 2.8492e-08 2 "" white 0 2 14 0
simundump text /kinetics/PKC/PKC_n_DAG_AA/notes 0 \
"Reduced kf to 0.66X to match Shinomura et al data. Initial: kf = 3.3333e-9 New: 2e-9 Newer: 2e-8 kb was 0.2 now 2."
call /kinetics/PKC/PKC_n_DAG_AA/notes LOAD \
"Reduced kf to 0.66X to match Shinomura et al data. Initial: kf = 3.3333e-9 New: 2e-9 Newer: 2e-8 kb was 0.2 now 2."
simundump group /kinetics/PLA2 0 1 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 20 0 0
simundump text /kinetics/PLA2/notes 0 ""
call /kinetics/PLA2/notes LOAD \
""
simundump kpool /kinetics/PLA2/PLA2_cytosolic 0 0.0 0.4 0.4 2.527e+05 \
2.527e+05 0 0 6.3175e+05 0 /kinetics/geometry yellow 1 20 1 0
simundump text /kinetics/PLA2/PLA2_cytosolic/notes 0 \
"Calculated cytosolic was 20 nm from Wijkander and Sundler However, Leslie and Channon use about 400 nM. Need to confirm, but this is the value I use here. Another recalc of W&S gives 1uM"
call /kinetics/PLA2/PLA2_cytosolic/notes LOAD \
"Calculated cytosolic was 20 nm from Wijkander and Sundler However, Leslie and Channon use about 400 nM. Need to confirm, but this is the value I use here. Another recalc of W&S gives 1uM"
simundump kpool /kinetics/PLA2/PLA2_Ca_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry yellow 1 22 1 0
simundump text /kinetics/PLA2/PLA2_Ca_p/notes 0 ""
call /kinetics/PLA2/PLA2_Ca_p/notes LOAD \
""
simundump kenz /kinetics/PLA2/PLA2_Ca_p/kenz 0 0 0 0 0 6.3175e+05 2.0295e-06 \
21.6 5.4 0 0 "" yellow red "" 22 3 0
simundump text /kinetics/PLA2/PLA2_Ca_p/kenz/notes 0 \
"10 x raise oct22 12 x oct 24, set k2 = 4 * k3"
call /kinetics/PLA2/PLA2_Ca_p/kenz/notes LOAD \
"10 x raise oct22 12 x oct 24, set k2 = 4 * k3"
simundump kpool /kinetics/PLA2/PIP2_PLA2_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 1 24 1 0
simundump text /kinetics/PLA2/PIP2_PLA2_p/notes 0 ""
call /kinetics/PLA2/PIP2_PLA2_p/notes LOAD \
""
simundump kenz /kinetics/PLA2/PIP2_PLA2_p/kenz 0 0 0 0 0 6.3175e+05 \
4.1492e-06 44.16 11.04 0 0 "" cyan red "" 24 3 0
simundump text /kinetics/PLA2/PIP2_PLA2_p/kenz/notes 0 \
"10 X raise oct 22 12 X further raise oct 24 to allow for correct conc of enzyme"
call /kinetics/PLA2/PIP2_PLA2_p/kenz/notes LOAD \
"10 X raise oct 22 12 X further raise oct 24 to allow for correct conc of enzyme"
simundump kpool /kinetics/PLA2/PIP2_Ca_PLA2_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 1 26 1 0
simundump text /kinetics/PLA2/PIP2_Ca_PLA2_p/notes 0 ""
call /kinetics/PLA2/PIP2_Ca_PLA2_p/notes LOAD \
""
simundump kenz /kinetics/PLA2/PIP2_Ca_PLA2_p/kenz 0 0 0 0 0 6.3175e+05 \
1.353e-05 144 36 0 0 "" cyan red "" 26 3 0
simundump text /kinetics/PLA2/PIP2_Ca_PLA2_p/kenz/notes 0 \
"10 x raise oct 22 12 x and rescale for k2 = 4 * k3 convention oct 24 Increase further to get the match to expt, which was spoilt due to large accumulation of PLA2 in the enzyme complexed forms. Lets raise k3, leaving the others at k1 = 1.5e-5 and k2 = 144 since they are large already."
call /kinetics/PLA2/PIP2_Ca_PLA2_p/kenz/notes LOAD \
"10 x raise oct 22 12 x and rescale for k2 = 4 * k3 convention oct 24 Increase further to get the match to expt, which was spoilt due to large accumulation of PLA2 in the enzyme complexed forms. Lets raise k3, leaving the others at k1 = 1.5e-5 and k2 = 144 since they are large already."
simundump kpool /kinetics/PLA2/DAG_Ca_PLA2_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 1 28 1 0
simundump text /kinetics/PLA2/DAG_Ca_PLA2_p/notes 0 ""
call /kinetics/PLA2/DAG_Ca_PLA2_p/notes LOAD \
""
simundump kenz /kinetics/PLA2/DAG_Ca_PLA2_p/kenz 0 0 0 0 0 6.3175e+05 \
2.255e-05 240 60 0 0 "" pink red "" 28 3 0
simundump text /kinetics/PLA2/DAG_Ca_PLA2_p/kenz/notes 0 \
"10 X raise oct 22 12 X raise oct 24 + conversion to k2 =4 * k3"
call /kinetics/PLA2/DAG_Ca_PLA2_p/kenz/notes LOAD \
"10 X raise oct 22 12 X raise oct 24 + conversion to k2 =4 * k3"
simundump kpool /kinetics/PLA2/APC 0 0.0 30.001 30.001 1.8953e+07 1.8953e+07 \
0 0 6.3175e+05 4 /kinetics/geometry yellow 1 30 1 0
simundump text /kinetics/PLA2/APC/notes 0 \
"arachodonylphosphatidylcholine is the favoured substrate from Wijkander and Sundler, JBC 202 pp 873-880, 1991. Their assay used 30 uM substrate, which is what the kinetics in this model are based on. For the later model we should locate a more realistic value for APC."
call /kinetics/PLA2/APC/notes LOAD \
"arachodonylphosphatidylcholine is the favoured substrate from Wijkander and Sundler, JBC 202 pp 873-880, 1991. Their assay used 30 uM substrate, which is what the kinetics in this model are based on. For the later model we should locate a more realistic value for APC."
simundump kpool /kinetics/PLA2/PLA2_p_Ca 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 1 32 1 0
simundump text /kinetics/PLA2/PLA2_p_Ca/notes 0 \
"Phosphorylated form of PLA2. Still need to hook it up using kinases. PKA: Wightman et al JBC 257 pp6650 1982 PKC: Many refs, eg Gronich et al JBC 263 pp 16645, 1988 but see Lin etal MAPK: Lin et al, Cell 72 pp 269, 1993. Show 3x with MAPK but not PKC alone Do not know if there is a Ca requirement for active phosphorylated state."
call /kinetics/PLA2/PLA2_p_Ca/notes LOAD \
"Phosphorylated form of PLA2. Still need to hook it up using kinases. PKA: Wightman et al JBC 257 pp6650 1982 PKC: Many refs, eg Gronich et al JBC 263 pp 16645, 1988 but see Lin etal MAPK: Lin et al, Cell 72 pp 269, 1993. Show 3x with MAPK but not PKC alone Do not know if there is a Ca requirement for active phosphorylated state."
simundump kenz /kinetics/PLA2/PLA2_p_Ca/kenz 0 0 0 0 0 6.3175e+05 4.51e-05 \
480 120 0 0 "" orange red "" 32 3 0
simundump text /kinetics/PLA2/PLA2_p_Ca/kenz/notes 0 \
"This form should be 3 to 6 times as fast as the Ca-only form. I have scaled by 4x which seems to give a 5x rise. 10x raise Oct 22 12 x oct 24, changed k2 = 4 * k3"
call /kinetics/PLA2/PLA2_p_Ca/kenz/notes LOAD \
"This form should be 3 to 6 times as fast as the Ca-only form. I have scaled by 4x which seems to give a 5x rise. 10x raise Oct 22 12 x oct 24, changed k2 = 4 * k3"
simundump kpool /kinetics/PLA2/PLA2_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 1 34 1 0
simundump text /kinetics/PLA2/PLA2_p/notes 0 ""
call /kinetics/PLA2/PLA2_p/notes LOAD \
""
simundump kreac /kinetics/PLA2/PLA2_Ca_act 0 1.5829e-06 0.1 "" white 1 21 3 0
simundump text /kinetics/PLA2/PLA2_Ca_act/notes 0 \
"Leslie and Channon BBA 1045 (1990) 261-270 fig6 pp267."
call /kinetics/PLA2/PLA2_Ca_act/notes LOAD \
"Leslie and Channon BBA 1045 (1990) 261-270 fig6 pp267."
simundump kreac /kinetics/PLA2/PIP2_PLA2_act 0 1.8995e-09 0.5 "" white 1 23 3 \
0
simundump text /kinetics/PLA2/PIP2_PLA2_act/notes 0 ""
call /kinetics/PLA2/PIP2_PLA2_act/notes LOAD \
""
simundump kreac /kinetics/PLA2/PIP2_Ca_PLA2_act 0 1.8995e-08 0.1 "" white 1 \
25 3 0
simundump text /kinetics/PLA2/PIP2_Ca_PLA2_act/notes 0 ""
call /kinetics/PLA2/PIP2_Ca_PLA2_act/notes LOAD \
""
simundump kreac /kinetics/PLA2/DAG_Ca_PLA2_act 0 4.7487e-09 4 "" white 1 27 3 \
0
simundump text /kinetics/PLA2/DAG_Ca_PLA2_act/notes 0 \
"27 June 1996 Scaled kf down by 0.015 from 3.33e-7 to 5e-9 to fit with revised DAG estimates and use of mole-fraction to calculate eff on PLA2."
call /kinetics/PLA2/DAG_Ca_PLA2_act/notes LOAD \
"27 June 1996 Scaled kf down by 0.015 from 3.33e-7 to 5e-9 to fit with revised DAG estimates and use of mole-fraction to calculate eff on PLA2."
simundump kreac /kinetics/PLA2/Degrade_AA 0 0.4 0 "" white 1 29 3 0
simundump text /kinetics/PLA2/Degrade_AA/notes 0 \
"I need to check if the AA degradation pathway really leads back to APC. Anyway, it is a convenient buffered pool to dump it back into. For the purposes of the full model we use a rate of degradation of 0.2/sec Raised decay to 0.4 : see PLA35.g notes for Feb17"
call /kinetics/PLA2/Degrade_AA/notes LOAD \
"I need to check if the AA degradation pathway really leads back to APC. Anyway, it is a convenient buffered pool to dump it back into. For the purposes of the full model we use a rate of degradation of 0.2/sec Raised decay to 0.4 : see PLA35.g notes for Feb17"
simundump kreac /kinetics/PLA2/PLA2_p_Ca_act 0 9.4974e-06 0.1 "" white 1 31 3 \
0
simundump text /kinetics/PLA2/PLA2_p_Ca_act/notes 0 \
"To start off, same kinetics as the PLA2-Ca-act direct pathway. Oops ! Missed out the Ca input to this pathway first time round. Let's raise the forward rate about 3x to 5e-6. This will let us reduce the rather high rates we have used for the kenz on PLA2*-Ca. In fact, it may be that the rates are not that different, just that this pathway for getting the PLA2 to the memb is more efficien...."
call /kinetics/PLA2/PLA2_p_Ca_act/notes LOAD \
"To start off, same kinetics as the PLA2-Ca-act direct pathway. Oops ! Missed out the Ca input to this pathway first time round. Let's raise the forward rate about 3x to 5e-6. This will let us reduce the rather high rates we have used for the kenz on PLA2*-Ca. In fact, it may be that the rates are not that different, just that this pathway for getting the PLA2 to the memb is more efficien...."
simundump kreac /kinetics/PLA2/dephosphorylate_PLA2_p 0 0.17 0 "" white 1 33 \
3 0
simundump text /kinetics/PLA2/dephosphorylate_PLA2_p/notes 0 ""
call /kinetics/PLA2/dephosphorylate_PLA2_p/notes LOAD \
""
simundump kpool /kinetics/PLA2/Arachidonic_Acid 0 0.0 0 0 0 0 0 0 6.3175e+05 \
0 /kinetics/geometry darkgreen 1 21 16 0
simundump text /kinetics/PLA2/Arachidonic_Acid/notes 0 ""
call /kinetics/PLA2/Arachidonic_Acid/notes LOAD \
""
simundump group /kinetics/PLCbeta 0 2 black x 0 0 "" defaultfile \
defaultfile.g 0 0 0 40 0 0
simundump text /kinetics/PLCbeta/notes 0 ""
call /kinetics/PLCbeta/notes LOAD \
""
simundump kpool /kinetics/PLCbeta/PLC 0 0.0 0.8 0.8 5.054e+05 5.054e+05 0 0 \
6.3175e+05 0 /kinetics/geometry cyan 2 40 1 0
simundump text /kinetics/PLCbeta/PLC/notes 0 \
"Total PLC = 0.8 uM see Ryu et al JBC 262 (26) pp 12511 1987"
call /kinetics/PLCbeta/PLC/notes LOAD \
"Total PLC = 0.8 uM see Ryu et al JBC 262 (26) pp 12511 1987"
simundump kpool /kinetics/PLCbeta/Inositol 0 0.0 0 0 0 0 0 0 6.3175e+05 4 \
/kinetics/geometry green 2 42 1 0
simundump text /kinetics/PLCbeta/Inositol/notes 0 ""
call /kinetics/PLCbeta/Inositol/notes LOAD \
""
simundump kpool /kinetics/PLCbeta/PLC_Ca 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 2 44 1 0
simundump text /kinetics/PLCbeta/PLC_Ca/notes 0 ""
call /kinetics/PLCbeta/PLC_Ca/notes LOAD \
""
simundump kenz /kinetics/PLCbeta/PLC_Ca/PLC_Ca 0 0 0 0 0 6.3175e+05 \
3.7884e-06 40 10 0 0 "" cyan red "" 44 3 0
simundump text /kinetics/PLCbeta/PLC_Ca/PLC_Ca/notes 0 \
"From Sternweis et al Phil Trans R Soc Lond 1992, also matched by Homma et al. k1 = 1.5e-5, now 4.2e-6 k2 = 70/sec; now 40/sec k3 = 17.5/sec; now 10/sec Note that the wording in Sternweis et al is ambiguous re the Km."
call /kinetics/PLCbeta/PLC_Ca/PLC_Ca/notes LOAD \
"From Sternweis et al Phil Trans R Soc Lond 1992, also matched by Homma et al. k1 = 1.5e-5, now 4.2e-6 k2 = 70/sec; now 40/sec k3 = 17.5/sec; now 10/sec Note that the wording in Sternweis et al is ambiguous re the Km."
simundump kpool /kinetics/PLCbeta/PLC_Ca_Gq 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 2 46 1 0
simundump text /kinetics/PLCbeta/PLC_Ca_Gq/notes 0 \
"This should really be labelled PLC-G*GTP-Ca. This is the activated form of the enzyme. Mahama and Linderman assume that the IP3 precursors are not rate-limiting, but I will include those for completeness as they may be needed later."
call /kinetics/PLCbeta/PLC_Ca_Gq/notes LOAD \
"This should really be labelled PLC-G*GTP-Ca. This is the activated form of the enzyme. Mahama and Linderman assume that the IP3 precursors are not rate-limiting, but I will include those for completeness as they may be needed later."
simundump kenz /kinetics/PLCbeta/PLC_Ca_Gq/PLCb_Ca_Gq 0 0 0 0 0 6.3175e+05 \
7.216e-05 192 48 0 0 "" cyan red "" 46 3 0
simundump text /kinetics/PLCbeta/PLC_Ca_Gq/PLCb_Ca_Gq/notes 0 \
"From Sternweis et al, Phil Trans R Soc Lond 1992, and the values from other refs eg Homma et al JBC 263(14) pp6592 1988 match. k1 = 5e-5/sec k2 = 240/sec; now 120/sec k3 = 60/sec; now 30/sec Note that the wording in Sternweis et al is ambiguous wr. to the Km for Gq vs non-Gq states of PLC. K1 is still a bit too low. Raise to 7e-5 9 Jun 1996: k1 was 0.0002, changed to 5e-5"
call /kinetics/PLCbeta/PLC_Ca_Gq/PLCb_Ca_Gq/notes LOAD \
"From Sternweis et al, Phil Trans R Soc Lond 1992, and the values from other refs eg Homma et al JBC 263(14) pp6592 1988 match. k1 = 5e-5/sec k2 = 240/sec; now 120/sec k3 = 60/sec; now 30/sec Note that the wording in Sternweis et al is ambiguous wr. to the Km for Gq vs non-Gq states of PLC. K1 is still a bit too low. Raise to 7e-5 9 Jun 1996: k1 was 0.0002, changed to 5e-5"
simundump kpool /kinetics/PLCbeta/PLC_Gq 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 2 48 1 0
simundump text /kinetics/PLCbeta/PLC_Gq/notes 0 ""
call /kinetics/PLCbeta/PLC_Gq/notes LOAD \
""
simundump kreac /kinetics/PLCbeta/Act_PLC_Ca 0 4.7487e-06 1 "" white 2 41 3 0
simundump text /kinetics/PLCbeta/Act_PLC_Ca/notes 0 \
"Affinity for Ca = 1uM without AlF, 0.1 with: from Smrcka et al science 251 pp 804-807 1991 so [Ca].kf = kb so kb/kf = 1 * 6e5 = 1/1.66e-6 11 June 1996: Raised affinity to 5e-6 to maintain balance. See notes."
call /kinetics/PLCbeta/Act_PLC_Ca/notes LOAD \
"Affinity for Ca = 1uM without AlF, 0.1 with: from Smrcka et al science 251 pp 804-807 1991 so [Ca].kf = kb so kb/kf = 1 * 6e5 = 1/1.66e-6 11 June 1996: Raised affinity to 5e-6 to maintain balance. See notes."
simundump kreac /kinetics/PLCbeta/Degrade_IP3 0 2.5 0 "" white 2 43 3 0
simundump text /kinetics/PLCbeta/Degrade_IP3/notes 0 \
"The enzyme is IP3 5-phosphomonesterase. about 45K. Actual products are Ins(1,4)P2, and cIns(1:2,4,5)P3. review in Majerus et al Science 234 1519-1526, 1986. Meyer and Stryer 1988 PNAS 85:5051-5055 est decay of IP3 at 1-3/sec"
call /kinetics/PLCbeta/Degrade_IP3/notes LOAD \
"The enzyme is IP3 5-phosphomonesterase. about 45K. Actual products are Ins(1,4)P2, and cIns(1:2,4,5)P3. review in Majerus et al Science 234 1519-1526, 1986. Meyer and Stryer 1988 PNAS 85:5051-5055 est decay of IP3 at 1-3/sec"
simundump kreac /kinetics/PLCbeta/Degrade_DAG 0 0.15 0 "" white 2 45 3 0
simundump text /kinetics/PLCbeta/Degrade_DAG/notes 0 \
"These rates are the same as for degrading IP3, but I am sure that they could be improved. Lets double kf to 0.2, since the amount of DAG in the cell should be <= 1uM. Need to double it again, for the same reason. kf now 0.5 27 June 1996 kf is now 0.02 to get 50 sec time course 30 Aug 1997: Raised kf to 0.11 to accomodate PLC_gamma 27 Mar 1998: kf now 0.15 for PLC_gamma"
call /kinetics/PLCbeta/Degrade_DAG/notes LOAD \
"These rates are the same as for degrading IP3, but I am sure that they could be improved. Lets double kf to 0.2, since the amount of DAG in the cell should be <= 1uM. Need to double it again, for the same reason. kf now 0.5 27 June 1996 kf is now 0.02 to get 50 sec time course 30 Aug 1997: Raised kf to 0.11 to accomodate PLC_gamma 27 Mar 1998: kf now 0.15 for PLC_gamma"
simundump kreac /kinetics/PLCbeta/Act_PLC_by_Gq 0 3.9889e-05 1 "" white 2 47 \
3 0
simundump text /kinetics/PLCbeta/Act_PLC_by_Gq/notes 0 \
"Affinity for Gq is > 20 nM (Smrcka et al Science251 804-807 1991) so [Gq].kf = kb so 40nM * 6e5 = kb/kf = 24e3 so kf = 4.2e-5, kb =1"
call /kinetics/PLCbeta/Act_PLC_by_Gq/notes LOAD \
"Affinity for Gq is > 20 nM (Smrcka et al Science251 804-807 1991) so [Gq].kf = kb so 40nM * 6e5 = kb/kf = 24e3 so kf = 4.2e-5, kb =1"
simundump kreac /kinetics/PLCbeta/Inact_PLC_Gq 0 0.0133 0 "" white 2 49 3 0
simundump text /kinetics/PLCbeta/Inact_PLC_Gq/notes 0 \
"This process is assumed to be directly caused by the inactivation of the G*GTP to G*GDP. Hence, kf = .013 /sec = 0.8/min, same as the rate for Inact-G. kb = 0 since this is irreversible. We may be interested in studying the role of PLC as a GAP. If so, the kf would be faster here than in Inact-G"
call /kinetics/PLCbeta/Inact_PLC_Gq/notes LOAD \
"This process is assumed to be directly caused by the inactivation of the G*GTP to G*GDP. Hence, kf = .013 /sec = 0.8/min, same as the rate for Inact-G. kb = 0 since this is irreversible. We may be interested in studying the role of PLC as a GAP. If so, the kf would be faster here than in Inact-G"
simundump kreac /kinetics/PLCbeta/PLC_bind_Gq 0 3.9889e-06 1 "" white 2 51 3 \
0
simundump text /kinetics/PLCbeta/PLC_bind_Gq/notes 0 \
"this binding does not produce active PLC. This step was needed to implement the described (Smrcka et al) increase in affinity for Ca by PLC once Gq was bound. The kinetics are the same as the binding step for Ca-PLC to Gq. June 1996: Changed the kf to 4.2e-5 to 4.2e-6 to preserve balance around the reactions."
call /kinetics/PLCbeta/PLC_bind_Gq/notes LOAD \
"this binding does not produce active PLC. This step was needed to implement the described (Smrcka et al) increase in affinity for Ca by PLC once Gq was bound. The kinetics are the same as the binding step for Ca-PLC to Gq. June 1996: Changed the kf to 4.2e-5 to 4.2e-6 to preserve balance around the reactions."
simundump kreac /kinetics/PLCbeta/PLC_Gq_bind_Ca 0 4.7487e-05 1 "" white 2 53 \
3 0
simundump text /kinetics/PLCbeta/PLC_Gq_bind_Ca/notes 0 \
"this step has a high affinity for Ca, from Smrcka et al. 0.1uM so kf /kb = 1/6e4 = 1.666e-5:1. See the Act-PLC-by-Gq reac. 11 Jun 1996: Raised kf to 5e-5 based on match to conc-eff curves from Smrcka et al."
call /kinetics/PLCbeta/PLC_Gq_bind_Ca/notes LOAD \
"this step has a high affinity for Ca, from Smrcka et al. 0.1uM so kf /kb = 1/6e4 = 1.666e-5:1. See the Act-PLC-by-Gq reac. 11 Jun 1996: Raised kf to 5e-5 based on match to conc-eff curves from Smrcka et al."
simundump kpool /kinetics/PLCbeta/PC 0 0.0 0 0 0 0 0 0 6.3175e+05 4 \
/kinetics/geometry green 2 50 1 0
simundump text /kinetics/PLCbeta/PC/notes 0 \
"Phosphatidylcholine is the main (around 55%) metabolic product of DAG, follwed by PE (around 25%). Ref is Welsh and Cabot, JCB35:231-245(1987)"
call /kinetics/PLCbeta/PC/notes LOAD \
"Phosphatidylcholine is the main (around 55%) metabolic product of DAG, follwed by PE (around 25%). Ref is Welsh and Cabot, JCB35:231-245(1987)"
simundump kpool /kinetics/PLCbeta/DAG 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry green 2 44 16 0
simundump text /kinetics/PLCbeta/DAG/notes 0 ""
call /kinetics/PLCbeta/DAG/notes LOAD \
""
simundump kpool /kinetics/PLCbeta/IP3 0 0.0 0.73 0.73 4.6118e+05 4.6118e+05 0 \
0 6.3175e+05 0 /kinetics/geometry pink 2 46 16 0
simundump text /kinetics/PLCbeta/IP3/notes 0 \
"Peak IP3 is perhaps 15 uM, basal <= 0.2 uM."
call /kinetics/PLCbeta/IP3/notes LOAD \
"Peak IP3 is perhaps 15 uM, basal <= 0.2 uM."
simundump kpool /kinetics/PLCbeta/PIP2 0 0.0 10 10 6.3175e+06 6.3175e+06 0 0 \
6.3175e+05 4 /kinetics/geometry green 2 52 1 0
simundump text /kinetics/PLCbeta/PIP2/notes 0 ""
call /kinetics/PLCbeta/PIP2/notes LOAD \
""
simundump group /kinetics/mGluRAntag 0 3 black x 0 0 "" defaultfile \
defaultfile.g 0 0 0 60 0 0
simundump text /kinetics/mGluRAntag/notes 0 ""
call /kinetics/mGluRAntag/notes LOAD \
""
simundump kpool /kinetics/mGluRAntag/mGluRAntag 0 0.0 0 0 0 0 0 0 6.3175e+05 \
4 /kinetics/geometry seagreen 3 60 1 0
simundump text /kinetics/mGluRAntag/mGluRAntag/notes 0 \
"I am implementing this as acting only on the Rec-Gq complex, based on a more complete model PLC_Gq48.g which showed that the binding to the rec alone contributed only a small amount."
call /kinetics/mGluRAntag/mGluRAntag/notes LOAD \
"I am implementing this as acting only on the Rec-Gq complex, based on a more complete model PLC_Gq48.g which showed that the binding to the rec alone contributed only a small amount."
simundump kpool /kinetics/mGluRAntag/Blocked_Rec_Gq 0 0.0 0 0 0 0 0 0 \
6.3175e+05 0 /kinetics/geometry seagreen 3 62 1 0
simundump text /kinetics/mGluRAntag/Blocked_Rec_Gq/notes 0 ""
call /kinetics/mGluRAntag/Blocked_Rec_Gq/notes LOAD \
""
simundump kreac /kinetics/mGluRAntag/Antag_bind_Rec_Gq 0 9.4974e-05 0.01 "" \
white 3 61 3 0
simundump text /kinetics/mGluRAntag/Antag_bind_Rec_Gq/notes 0 \
"The rate consts give a total binding affinity of only"
call /kinetics/mGluRAntag/Antag_bind_Rec_Gq/notes LOAD \
"The rate consts give a total binding affinity of only"
simundump group /kinetics/MAPK 0 4 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 80 0 0
simundump text /kinetics/MAPK/notes 0 ""
call /kinetics/MAPK/notes LOAD \
""
simundump kpool /kinetics/MAPK/MAPK_p_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 4 80 1 0
simundump text /kinetics/MAPK/MAPK_p_p/notes 0 ""
call /kinetics/MAPK/MAPK_p_p/notes LOAD \
""
simundump kenz /kinetics/MAPK/MAPK_p_p/MAPK_p_p 0 0 0 0 0 6.3175e+05 \
5.863e-06 80 20 0 0 "" orange red "" 80 3 0
simundump text /kinetics/MAPK/MAPK_p_p/MAPK_p_p/notes 0 \
"Km = 25uM @ 50 uM ATP and 1mg/ml MBP (huge XS of substrate) Vmax = 4124 pmol/min/ml at a conc of 125 pmol/ml of enz, so: k3 = .5/sec (rate limiting) k1 = (k2 + k3)/Km = (.5 + 0)/(25*6e5) = 2e-8 (#/cell)^-1 #s from Sanghera et al JBC 265 pp 52 , 1990. From Nemenoff et al JBC 268(3):1960-1964 - using Sanghera's 1e-4 ratio of MAPK to protein, we get k3 = 7/sec from 1000 pmol/min/mg fig 5"
call /kinetics/MAPK/MAPK_p_p/MAPK_p_p/notes LOAD \
"Km = 25uM @ 50 uM ATP and 1mg/ml MBP (huge XS of substrate) Vmax = 4124 pmol/min/ml at a conc of 125 pmol/ml of enz, so: k3 = .5/sec (rate limiting) k1 = (k2 + k3)/Km = (.5 + 0)/(25*6e5) = 2e-8 (#/cell)^-1 #s from Sanghera et al JBC 265 pp 52 , 1990. From Nemenoff et al JBC 268(3):1960-1964 - using Sanghera's 1e-4 ratio of MAPK to protein, we get k3 = 7/sec from 1000 pmol/min/mg fig 5"
simundump kenz /kinetics/MAPK/MAPK_p_p/MAPK_p_p_feedback 0 0 0 0 0 6.3175e+05 \
2.9315e-06 40 10 0 0 "" orange red "" 80 4 0
simundump text /kinetics/MAPK/MAPK_p_p/MAPK_p_p_feedback/notes 0 \
"Ueki et al JBC 269(22):15756-15761 show the presence of this step, but not the rate consts, which are derived from Sanghera et al JBC 265(1):52-57, 1990, see the deriv in the MAPK* notes."
call /kinetics/MAPK/MAPK_p_p/MAPK_p_p_feedback/notes LOAD \
"Ueki et al JBC 269(22):15756-15761 show the presence of this step, but not the rate consts, which are derived from Sanghera et al JBC 265(1):52-57, 1990, see the deriv in the MAPK* notes."
simundump kenz /kinetics/MAPK/MAPK_p_p/phosph_Sos 0 0 0 0 0 6.3175e+05 \
2.9315e-05 40 10 0 0 "" orange red "" 80 5 0
simundump text /kinetics/MAPK/MAPK_p_p/phosph_Sos/notes 0 \
"See Porfiri and McCormick JBC 271:10 pp5871 1996 for the existence of this step. We'll take the rates from the ones used for the phosph of Raf by MAPK. Sep 17 1997: The transient activation curve matches better with k1 up by 10 x."
call /kinetics/MAPK/MAPK_p_p/phosph_Sos/notes LOAD \
"See Porfiri and McCormick JBC 271:10 pp5871 1996 for the existence of this step. We'll take the rates from the ones used for the phosph of Raf by MAPK. Sep 17 1997: The transient activation curve matches better with k1 up by 10 x."
simundump kenz /kinetics/MAPK/MAPK_p_p/cluster_phospho_S6K 0 0 0 0 0 \
6.3175e+05 4.4094e-07 4 1 0 0 "" red red "" 80 6 0
simundump text /kinetics/MAPK/MAPK_p_p/cluster_phospho_S6K/notes 0 ""
call /kinetics/MAPK/MAPK_p_p/cluster_phospho_S6K/notes LOAD \
""
simundump kenz /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_p_p 0 0 0 0 0 6.3175e+05 \
6.8186e-08 0.64 0.16 0 0 "" red orange "" 80 7 0
simundump text /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_p_p/notes 0 ""
call /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_p_p/notes LOAD \
""
simundump kenz /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_phospho 0 0 0 0 0 \
6.3175e+05 6.8186e-08 0.64 0.16 0 0 "" red 47 "" 80 8 0
simundump text /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_phospho/notes 0 ""
call /kinetics/MAPK/MAPK_p_p/MAPK_4E_BP_phospho/notes LOAD \
""
simundump kpool /kinetics/MAPK/craf_1 0 0.0 0.2 0.2 1.2635e+05 1.2635e+05 0 0 \
6.3175e+05 0 /kinetics/geometry pink 4 82 1 0
simundump text /kinetics/MAPK/craf_1/notes 0 \
"Couldn't find any ref to the actual conc of craf-1 but I should try Strom et al Oncogene 5 pp 345 In line with the other kinases in the cascade, I estimate the conc to be 0.2 uM. To init we use 0.15, which is close to equil"
call /kinetics/MAPK/craf_1/notes LOAD \
"Couldn't find any ref to the actual conc of craf-1 but I should try Strom et al Oncogene 5 pp 345 In line with the other kinases in the cascade, I estimate the conc to be 0.2 uM. To init we use 0.15, which is close to equil"
simundump kpool /kinetics/MAPK/craf_1_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 4 84 1 0
simundump text /kinetics/MAPK/craf_1_p/notes 0 ""
call /kinetics/MAPK/craf_1_p/notes LOAD \
""
simundump kenz /kinetics/MAPK/craf_1_p/MEK_phospho 0 0 0 0 0 6.3175e+05 \
4.7948e-06 0.42 0.105 0 0 "" 0 black "" 84 3 0
simundump text /kinetics/MAPK/craf_1_p/MEK_phospho/notes 0 ""
call /kinetics/MAPK/craf_1_p/MEK_phospho/notes LOAD \
""
simundump kenz /kinetics/MAPK/craf_1_p/MEKp_phospho 0 0 0 0 0 6.3175e+05 \
4.7948e-06 0.42 0.105 0 0 "" 63 black "" 84 4 0
simundump text /kinetics/MAPK/craf_1_p/MEKp_phospho/notes 0 ""
call /kinetics/MAPK/craf_1_p/MEKp_phospho/notes LOAD \
""
simundump kpool /kinetics/MAPK/MAPKK 0 0.0 0.18001 0.18001 1.1372e+05 \
1.1372e+05 0 0 6.3175e+05 0 /kinetics/geometry pink 4 86 1 0
simundump text /kinetics/MAPK/MAPKK/notes 0 \
"Conc is from Seger et al JBC 267:20 pp14373 (1992) mwt is 45/46 Kd We assume that phosphorylation on both ser and thr is needed for activiation. See Kyriakis et al Nature 358 417 1992 Init conc of total is 0.18"
call /kinetics/MAPK/MAPKK/notes LOAD \
"Conc is from Seger et al JBC 267:20 pp14373 (1992) mwt is 45/46 Kd We assume that phosphorylation on both ser and thr is needed for activiation. See Kyriakis et al Nature 358 417 1992 Init conc of total is 0.18"
simundump kpool /kinetics/MAPK/MAPK 0 0.0 0.36 0.36 2.2743e+05 2.2743e+05 0 0 \
6.3175e+05 0 /kinetics/geometry pink 4 88 1 0
simundump text /kinetics/MAPK/MAPK/notes 0 \
"conc is from Sanghera et al JBC 265 pp 52 (1990) A second calculation gives 3.1 uM, from same paper. They est MAPK is 1e-4x total protein, and protein is 15% of cell wt, so MAPK is 1.5e-5g/ml = 0.36uM. which is closer to our first estimate. Lets use this."
call /kinetics/MAPK/MAPK/notes LOAD \
"conc is from Sanghera et al JBC 265 pp 52 (1990) A second calculation gives 3.1 uM, from same paper. They est MAPK is 1e-4x total protein, and protein is 15% of cell wt, so MAPK is 1.5e-5g/ml = 0.36uM. which is closer to our first estimate. Lets use this."
simundump kpool /kinetics/MAPK/craf_1_p_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry hotpink 4 90 1 0
simundump text /kinetics/MAPK/craf_1_p_p/notes 0 \
"Negative feedback by MAPK* by hyperphosphorylating craf-1* gives rise to this pool. Ueki et al JBC 269(22):15756-15761, 1994"
call /kinetics/MAPK/craf_1_p_p/notes LOAD \
"Negative feedback by MAPK* by hyperphosphorylating craf-1* gives rise to this pool. Ueki et al JBC 269(22):15756-15761, 1994"
simundump kpool /kinetics/MAPK/MAPK_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 4 92 1 0
simundump text /kinetics/MAPK/MAPK_p/notes 0 \
"Haystead et al FEBS Lett. 306(1) pp 17-22 show that phosphorylation is strictly sequential, first tyr185 then thr183."
call /kinetics/MAPK/MAPK_p/notes LOAD \
"Haystead et al FEBS Lett. 306(1) pp 17-22 show that phosphorylation is strictly sequential, first tyr185 then thr183."
simundump kpool /kinetics/MAPK/MAPKK_p_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 4 94 1 0
simundump text /kinetics/MAPK/MAPKK_p_p/notes 0 \
"MAPKK phosphorylates MAPK on both the tyr and thr residues, first tyr then thr. Refs: Seger et al JBC267:20 pp 14373 1992 The MAPKK itself is phosphorylated on ser as well as thr residues. Let us assume that the ser goes first, and that the sequential phosphorylation is needed. See Kyriakis et al Nature 358 417-421 1992"
call /kinetics/MAPK/MAPKK_p_p/notes LOAD \
"MAPKK phosphorylates MAPK on both the tyr and thr residues, first tyr then thr. Refs: Seger et al JBC267:20 pp 14373 1992 The MAPKK itself is phosphorylated on ser as well as thr residues. Let us assume that the ser goes first, and that the sequential phosphorylation is needed. See Kyriakis et al Nature 358 417-421 1992"
simundump kenz /kinetics/MAPK/MAPKK_p_p/MAPKKtyr 0 0 0 0 0 6.3175e+05 \
2.4354e-05 0.6 0.15 0 0 "" pink red "" 94 3 0
simundump text /kinetics/MAPK/MAPKK_p_p/MAPKKtyr/notes 0 \
"The actual MAPKK is 2 forms from Seger et al JBC 267:20 14373(1992) Vmax = 150nmol/min/mg From Haystead et al FEBS 306(1):17-22 we get Km=46.6nM for at least one of the phosphs. Putting these together: k3=0.15/sec, scale to get k2=0.6. k1=0.75/46.6nM=2.7e-5"
call /kinetics/MAPK/MAPKK_p_p/MAPKKtyr/notes LOAD \
"The actual MAPKK is 2 forms from Seger et al JBC 267:20 14373(1992) Vmax = 150nmol/min/mg From Haystead et al FEBS 306(1):17-22 we get Km=46.6nM for at least one of the phosphs. Putting these together: k3=0.15/sec, scale to get k2=0.6. k1=0.75/46.6nM=2.7e-5"
simundump kenz /kinetics/MAPK/MAPKK_p_p/MAPKKthr 0 0 0 0 0 6.3175e+05 \
2.4354e-05 0.6 0.15 0 0 "" pink red "" 94 4 0
simundump text /kinetics/MAPK/MAPKK_p_p/MAPKKthr/notes 0 \
"Rate consts same as for MAPKKtyr."
call /kinetics/MAPK/MAPKK_p_p/MAPKKthr/notes LOAD \
"Rate consts same as for MAPKKtyr."
simundump kpool /kinetics/MAPK/MAPKK_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 4 96 1 0
simundump text /kinetics/MAPK/MAPKK_p/notes 0 \
"Intermediately phophorylated, assumed inactive, form of MAPKK"
call /kinetics/MAPK/MAPKK_p/notes LOAD \
"Intermediately phophorylated, assumed inactive, form of MAPKK"
simundump kpool /kinetics/MAPK/Raf_p_GTP_Ras 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 4 81 16 0
simundump text /kinetics/MAPK/Raf_p_GTP_Ras/notes 0 ""
call /kinetics/MAPK/Raf_p_GTP_Ras/notes LOAD \
""
simundump kenz /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.1 0 0 0 0 0 \
6.3175e+05 4.961e-06 0.42 0.105 0 0 "" red red "" 81 14 0
simundump text /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.1/notes 0 \
"Kinetics are the same as for the craf-1* activity, ie., k1=1.1e-6, k2=.42, k3 =0.105 These are based on Force et al PNAS USA 91 1270-1274 1994. These parms cannot reach the observed 4X stim of MAPK. So lets increase the affinity, ie, raise k1 10X to 1.1e-5 Lets take it back down to where it was. Back up to 5X: 5.5e-6"
call /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.1/notes LOAD \
"Kinetics are the same as for the craf-1* activity, ie., k1=1.1e-6, k2=.42, k3 =0.105 These are based on Force et al PNAS USA 91 1270-1274 1994. These parms cannot reach the observed 4X stim of MAPK. So lets increase the affinity, ie, raise k1 10X to 1.1e-5 Lets take it back down to where it was. Back up to 5X: 5.5e-6"
simundump kenz /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.2 0 0 0 0 0 \
6.3175e+05 4.961e-06 0.42 0.105 0 0 "" red red "" 81 13 0
simundump text /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.2/notes 0 \
"Same kinetics as other c-raf activated forms. See Force et al PNAS 91 1270-1274 1994. k1 = 1.1e-6, k2 = .42, k3 = 1.05 raise k1 to 5.5e-6"
call /kinetics/MAPK/Raf_p_GTP_Ras/Raf_p_GTP_Ras.2/notes LOAD \
"Same kinetics as other c-raf activated forms. See Force et al PNAS 91 1270-1274 1994. k1 = 1.1e-6, k2 = .42, k3 = 1.05 raise k1 to 5.5e-6"
simundump kpool /kinetics/MAPK/craf_1_p_ser259 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry 27 4 83 16 0
simundump text /kinetics/MAPK/craf_1_p_ser259/notes 0 ""
call /kinetics/MAPK/craf_1_p_ser259/notes LOAD \
""
simundump kreac /kinetics/MAPK/mGluR_barr2_Raf_scaffolding 0 2.3743e-09 0.001 \
"" white 4 81 3 0
simundump text /kinetics/MAPK/mGluR_barr2_Raf_scaffolding/notes 0 ""
call /kinetics/MAPK/mGluR_barr2_Raf_scaffolding/notes LOAD \
""
simundump group /kinetics/Ras 0 5 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 100 0 0
simundump text /kinetics/Ras/notes 0 ""
call /kinetics/Ras/notes LOAD \
""
simundump kpool /kinetics/Ras/GEF_Gprot_bg 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry hotpink 5 100 1 0
simundump text /kinetics/Ras/GEF_Gprot_bg/notes 0 \
"Guanine nucleotide exchange factor. This activates raf by exchanging bound GDP with GTP. I have left the GDP/GTP out of this reaction, it would be trivial to put them in. See Boguski & McCormick. Possible candidate molecules: RasGRF, smgGDS, Vav (in dispute). rasGRF: Kcat= 1.2/min Km = 680 nM smgGDS: Kcat: 0.37 /min, Km = 220 nM. vav: Turnover up over baseline by 10X,"
call /kinetics/Ras/GEF_Gprot_bg/notes LOAD \
"Guanine nucleotide exchange factor. This activates raf by exchanging bound GDP with GTP. I have left the GDP/GTP out of this reaction, it would be trivial to put them in. See Boguski & McCormick. Possible candidate molecules: RasGRF, smgGDS, Vav (in dispute). rasGRF: Kcat= 1.2/min Km = 680 nM smgGDS: Kcat: 0.37 /min, Km = 220 nM. vav: Turnover up over baseline by 10X,"
simundump kenz /kinetics/Ras/GEF_Gprot_bg/GEF_bg_act_Ras 0 0 0 0 0 6.3175e+05 \
2.9766e-07 0.08 0.02 0 0 "" hotpink red "" 100 3 0
simundump text /kinetics/Ras/GEF_Gprot_bg/GEF_bg_act_Ras/notes 0 \
"Kinetics based on the activation of Gq by the receptor complex in the Gq model (in turn based on the Mahama and Linderman model) k1 = 2e-5, k2 = 1e-10, k3 = 10 (I do not know why they even bother with k2). Lets put k1 at 2e-6 to get a reasonable equilibrium More specific values from, eg.g: Orita et al JBC 268(34) 25542-25546 from rasGRF and smgGDS: k1=3.3e-7; k2 = 0.08, k3 = 0.02"
call /kinetics/Ras/GEF_Gprot_bg/GEF_bg_act_Ras/notes LOAD \
"Kinetics based on the activation of Gq by the receptor complex in the Gq model (in turn based on the Mahama and Linderman model) k1 = 2e-5, k2 = 1e-10, k3 = 10 (I do not know why they even bother with k2). Lets put k1 at 2e-6 to get a reasonable equilibrium More specific values from, eg.g: Orita et al JBC 268(34) 25542-25546 from rasGRF and smgGDS: k1=3.3e-7; k2 = 0.08, k3 = 0.02"
simundump kpool /kinetics/Ras/inact_GEF 0 0.0 0.1 0.1 63175 63175 0 0 \
6.3175e+05 0 /kinetics/geometry hotpink 5 102 1 0
simundump text /kinetics/Ras/inact_GEF/notes 0 \
"Assume that SoS is present only at 50 nM. Revised to 100 nM to get equil to experimentally known levels."
call /kinetics/Ras/inact_GEF/notes LOAD \
"Assume that SoS is present only at 50 nM. Revised to 100 nM to get equil to experimentally known levels."
simundump kenz /kinetics/Ras/inact_GEF/basal_GEF_activity 0 0 0 0 0 \
6.3175e+05 1.5104e-08 0.08 0.02 0 0 "" red hotpink "" 102 3 0
simundump text /kinetics/Ras/inact_GEF/basal_GEF_activity/notes 0 ""
call /kinetics/Ras/inact_GEF/basal_GEF_activity/notes LOAD \
""
simundump kpool /kinetics/Ras/GEF_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry hotpink 5 104 1 0
simundump text /kinetics/Ras/GEF_p/notes 0 \
"phosphorylated and thereby activated form of GEF. See, e.g. Orita et al JBC 268:34 25542-25546 1993, Gulbins et al. It is not clear whether there is major specificity for tyr or ser/thr."
call /kinetics/Ras/GEF_p/notes LOAD \
"phosphorylated and thereby activated form of GEF. See, e.g. Orita et al JBC 268:34 25542-25546 1993, Gulbins et al. It is not clear whether there is major specificity for tyr or ser/thr."
simundump kenz /kinetics/Ras/GEF_p/GEF_p_act_Ras 0 0 0 0 0 6.3175e+05 \
2.9766e-07 0.08 0.02 0 0 "" hotpink red "" 104 3 0
simundump text /kinetics/Ras/GEF_p/GEF_p_act_Ras/notes 0 \
"Kinetics same as GEF-bg-act-ras"
call /kinetics/Ras/GEF_p/GEF_p_act_Ras/notes LOAD \
"Kinetics same as GEF-bg-act-ras"
simundump kpool /kinetics/Ras/GTP_Ras 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 5 106 1 0
simundump text /kinetics/Ras/GTP_Ras/notes 0 \
"Only a very small fraction (7% unstim, 15% stim) of ras is GTP-bound. Gibbs et al JBC 265(33) 20437"
call /kinetics/Ras/GTP_Ras/notes LOAD \
"Only a very small fraction (7% unstim, 15% stim) of ras is GTP-bound. Gibbs et al JBC 265(33) 20437"
simundump kpool /kinetics/Ras/GDP_Ras 0 0.0 0.2 0.2 1.2635e+05 1.2635e+05 0 0 \
6.3175e+05 0 /kinetics/geometry pink 5 108 1 0
simundump text /kinetics/Ras/GDP_Ras/notes 0 \
"GDP bound form. See Rosen et al Neuron 12 1207-1221 June 1994. the activation loop is based on Boguski and McCormick Nature 366 643-654 93 Assume Ras is present at about the same level as craf-1, 0.2 uM. Hallberg et al JBC 269:6 3913-3916 1994 estimate upto 5-10% of cellular Raf is assoc with Ras. Given that only 5-10% of Ras is GTP-bound, we need similar amounts of Ras as Raf."
call /kinetics/Ras/GDP_Ras/notes LOAD \
"GDP bound form. See Rosen et al Neuron 12 1207-1221 June 1994. the activation loop is based on Boguski and McCormick Nature 366 643-654 93 Assume Ras is present at about the same level as craf-1, 0.2 uM. Hallberg et al JBC 269:6 3913-3916 1994 estimate upto 5-10% of cellular Raf is assoc with Ras. Given that only 5-10% of Ras is GTP-bound, we need similar amounts of Ras as Raf."
simundump kpool /kinetics/Ras/GAP_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 5 110 1 0
simundump text /kinetics/Ras/GAP_p/notes 0 ""
call /kinetics/Ras/GAP_p/notes LOAD \
""
simundump kpool /kinetics/Ras/GAP 0 0.0 0.002 0.002 1263.5 1263.5 0 0 \
6.3175e+05 0 /kinetics/geometry red 5 112 1 0
simundump text /kinetics/Ras/GAP/notes 0 \
"GTPase-activating proteins. See Boguski and McCormick. Turn off Ras by helping to hydrolyze bound GTP. This one is probably NF1, ie., Neurofibromin as it is inhibited by AA and lipids, and expressed in neural cells. p120-GAP is also a possible candidate, but is less regulated. Both may exist at similar levels. See Eccleston et al JBC 268(36) pp27012-19 Level=.002"
call /kinetics/Ras/GAP/notes LOAD \
"GTPase-activating proteins. See Boguski and McCormick. Turn off Ras by helping to hydrolyze bound GTP. This one is probably NF1, ie., Neurofibromin as it is inhibited by AA and lipids, and expressed in neural cells. p120-GAP is also a possible candidate, but is less regulated. Both may exist at similar levels. See Eccleston et al JBC 268(36) pp27012-19 Level=.002"
simundump kenz /kinetics/Ras/GAP/GAP_inact_Ras 0 0 0 0 0 6.3175e+05 \
7.4394e-05 40 10 0 0 "" red red "" 112 3 0
simundump text /kinetics/Ras/GAP/GAP_inact_Ras/notes 0 \
"From Eccleston et al JBC 268(36)pp27012-19 get Kd < 2uM, kcat - 10/sec From Martin et al Cell 63 843-849 1990 get Kd ~ 250 nM, kcat = 20/min I will go with the Eccleston figures as there are good error bars (10%). In general the values are reasonably close. k1 = 1.666e-3/sec, k2 = 1000/sec, k3 = 10/sec (note k3 is rate-limiting) 5 Nov 2002: Changed ratio term to 4 from 100. Now we have k...."
call /kinetics/Ras/GAP/GAP_inact_Ras/notes LOAD \
"From Eccleston et al JBC 268(36)pp27012-19 get Kd < 2uM, kcat - 10/sec From Martin et al Cell 63 843-849 1990 get Kd ~ 250 nM, kcat = 20/min I will go with the Eccleston figures as there are good error bars (10%). In general the values are reasonably close. k1 = 1.666e-3/sec, k2 = 1000/sec, k3 = 10/sec (note k3 is rate-limiting) 5 Nov 2002: Changed ratio term to 4 from 100. Now we have k1=8.25e-5; k2=40, k3=10. k3 is still rate-limiting."
simundump kpool /kinetics/Ras/inact_GEF_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry hotpink 5 114 1 0
simundump text /kinetics/Ras/inact_GEF_p/notes 0 \
"Phosphorylation-inactivated form of GEF. See Hordijk et al JBC 269:5 3534-3538 1994 and Buregering et al EMBO J 12:11 4211-4220 1993"
call /kinetics/Ras/inact_GEF_p/notes LOAD \
"Phosphorylation-inactivated form of GEF. See Hordijk et al JBC 269:5 3534-3538 1994 and Buregering et al EMBO J 12:11 4211-4220 1993"
simundump kpool /kinetics/Ras/CaM_GEF 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 5 116 1 0
simundump text /kinetics/Ras/CaM_GEF/notes 0 \
"See Farnsworth et al Nature 376 524-527 1995"
call /kinetics/Ras/CaM_GEF/notes LOAD \
"See Farnsworth et al Nature 376 524-527 1995"
simundump kenz /kinetics/Ras/CaM_GEF/CaM_GEF_act_Ras 0 0 0 0 0 6.3175e+05 \
2.9766e-07 0.08 0.02 0 0 "" pink red "" 116 3 0
simundump text /kinetics/Ras/CaM_GEF/CaM_GEF_act_Ras/notes 0 \
"Kinetics same as GEF-bg_act-ras"
call /kinetics/Ras/CaM_GEF/CaM_GEF_act_Ras/notes LOAD \
"Kinetics same as GEF-bg_act-ras"
simundump kreac /kinetics/Ras/Ras_act_craf 0 3.799e-05 0.5 "" white 5 101 3 0
simundump text /kinetics/Ras/Ras_act_craf/notes 0 \
"Assume the binding is fast and limited only by the amount of Ras* available. So kf=kb/[craf-1] If kb is 1/sec, then kf = 1/0.2 uM = 1/(0.2 * 6e5) = 8.3e-6 Later: Raise it by 10 X to 4e-5 From Hallberg et al JBC 269:6 3913-3916 1994, 3% of cellular Raf is complexed with Ras. So we raise kb 4x to 4 This step needed to memb-anchor and activate Raf: Leevers et al Nature 369 411-414 (I don't....."
call /kinetics/Ras/Ras_act_craf/notes LOAD \
"Assume the binding is fast and limited only by the amount of Ras* available. So kf=kb/[craf-1] If kb is 1/sec, then kf = 1/0.2 uM = 1/(0.2 * 6e5) = 8.3e-6 Later: Raise it by 10 X to 4e-5 From Hallberg et al JBC 269:6 3913-3916 1994, 3% of cellular Raf is complexed with Ras. So we raise kb 4x to 4 This step needed to memb-anchor and activate Raf: Leevers et al Nature 369 411-414 (I don't...."
simundump kreac /kinetics/Ras/bg_act_GEF 0 9.4974e-06 1 "" white 5 103 3 0
simundump text /kinetics/Ras/bg_act_GEF/notes 0 \
"SoS/GEF is present at 50 nM ie 3e4/cell. BetaGamma maxes out at 9e4. Assume we have 1/3 of the GEF active when the BetaGamma is 1.5e4. so 1e4 * kb = 2e4 * 1.5e4 * kf, so kf/kb = 3e-5. The rate of this equil should be reasonably fast, say 1/sec"
call /kinetics/Ras/bg_act_GEF/notes LOAD \
"SoS/GEF is present at 50 nM ie 3e4/cell. BetaGamma maxes out at 9e4. Assume we have 1/3 of the GEF active when the BetaGamma is 1.5e4. so 1e4 * kb = 2e4 * 1.5e4 * kf, so kf/kb = 3e-5. The rate of this equil should be reasonably fast, say 1/sec"
simundump kreac /kinetics/Ras/dephosph_GEF 0 1 0 "" white 5 105 3 0
simundump text /kinetics/Ras/dephosph_GEF/notes 0 ""
call /kinetics/Ras/dephosph_GEF/notes LOAD \
""
simundump kreac /kinetics/Ras/Ras_intrinsic_GTPase 0 0.0001 0 "" white 5 107 \
3 0
simundump text /kinetics/Ras/Ras_intrinsic_GTPase/notes 0 \
"This is extremely slow (1e-4), but it is significant as so little GAP actually gets complexed with it that the total GTP turnover rises only by 2-3 X (see Gibbs et al, JBC 265(33) 20437-20422) and Eccleston et al JBC 268(36) 27012-27019 kf = 1e-4"
call /kinetics/Ras/Ras_intrinsic_GTPase/notes LOAD \
"This is extremely slow (1e-4), but it is significant as so little GAP actually gets complexed with it that the total GTP turnover rises only by 2-3 X (see Gibbs et al, JBC 265(33) 20437-20422) and Eccleston et al JBC 268(36) 27012-27019 kf = 1e-4"
simundump kreac /kinetics/Ras/dephosph_GAP 0 0.1 0 "" white 5 109 3 0
simundump text /kinetics/Ras/dephosph_GAP/notes 0 \
"Assume a reasonably good rate for dephosphorylating it, 1/sec"
call /kinetics/Ras/dephosph_GAP/notes LOAD \
"Assume a reasonably good rate for dephosphorylating it, 1/sec"
simundump kreac /kinetics/Ras/CaM_bind_GEF 0 9.4974e-05 1 "" white 5 111 3 0
simundump text /kinetics/Ras/CaM_bind_GEF/notes 0 \
"We have no numbers for this. It is probably between the two extremes represented by the CaMKII phosph states, and I have used guesses based on this. kf=1e-4 kb=1 The reaction is based on Farnsworth et al Nature 376 524-527 1995"
call /kinetics/Ras/CaM_bind_GEF/notes LOAD \
"We have no numbers for this. It is probably between the two extremes represented by the CaMKII phosph states, and I have used guesses based on this. kf=1e-4 kb=1 The reaction is based on Farnsworth et al Nature 376 524-527 1995"
simundump kreac /kinetics/Ras/dephosph_inact_GEF_p 0 1 0 "" white 5 113 3 0
simundump text /kinetics/Ras/dephosph_inact_GEF_p/notes 0 ""
call /kinetics/Ras/dephosph_inact_GEF_p/notes LOAD \
""
simundump group /kinetics/EGFR 0 6 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 120 0 0
simundump text /kinetics/EGFR/notes 0 ""
call /kinetics/EGFR/notes LOAD \
""
simundump kpool /kinetics/EGFR/EGFR 0 0.0 0.16666 0.16666 1.0529e+05 \
1.0529e+05 0 0 6.3175e+05 0 /kinetics/geometry red 6 120 1 0
simundump text /kinetics/EGFR/EGFR/notes 0 \
"Berkers et al JBC 266 say 22K hi aff recs. Sherrill and Kyte Biochemistry 35 use range 4-200 nM. These match, lets use them."
call /kinetics/EGFR/EGFR/notes LOAD \
"Berkers et al JBC 266 say 22K hi aff recs. Sherrill and Kyte Biochemistry 35 use range 4-200 nM. These match, lets use them."
simundump kpool /kinetics/EGFR/L.EGFR 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 6 122 1 0
simundump text /kinetics/EGFR/L.EGFR/notes 0 \
"This is terribly simplified: there are many interesting intermediate stages, including dimerization and assoc with adapter molecules like Shc, that contribute to the activation of the EGFR."
call /kinetics/EGFR/L.EGFR/notes LOAD \
"This is terribly simplified: there are many interesting intermediate stages, including dimerization and assoc with adapter molecules like Shc, that contribute to the activation of the EGFR."
simundump kenz /kinetics/EGFR/L.EGFR/Ca.PLC_g_phospho 0 0 0 0 0 6.3175e+05 \
4.51e-06 0.8 0.2 0 0 "" red red "" 122 3 0
simundump text /kinetics/EGFR/L.EGFR/Ca.PLC_g_phospho/notes 0 \
"Hsu et al JBC 266:1 603-608 1991 Km = 385 +- 100 uM, Vm = 5.1 +-1 pmol/min/ug for PLC-771. Other sites have similar range, but are not stim as much by EGF. k1 = 2.8e-2/385/6e5 = 1.2e-10. Phenomenally slow. But Sherrill and Kyte say turnover # for angiotensin II is 5/min for cell extt, and 2/min for placental. Also see Okada et al for Shc rates which are much faster."
call /kinetics/EGFR/L.EGFR/Ca.PLC_g_phospho/notes LOAD \
"Hsu et al JBC 266:1 603-608 1991 Km = 385 +- 100 uM, Vm = 5.1 +-1 pmol/min/ug for PLC-771. Other sites have similar range, but are not stim as much by EGF. k1 = 2.8e-2/385/6e5 = 1.2e-10. Phenomenally slow. But Sherrill and Kyte say turnover # for angiotensin II is 5/min for cell extt, and 2/min for placental. Also see Okada et al for Shc rates which are much faster."
simundump kenz /kinetics/EGFR/L.EGFR/SHC_phospho 0 0 0 0 0 6.3175e+05 \
1.804e-06 0.8 0.2 0 0 "" red red "" 122 4 0
simundump text /kinetics/EGFR/L.EGFR/SHC_phospho/notes 0 \
"Rates from Okada et al JBC 270:35 pp 20737 1995 Km = 0.70 to 0.85 uM, Vmax = 4.4 to 5.0 pmol/min. Unfortunately the amount of enzyme is not known, the prep is only partially purified. Time course of phosph is max within 30 sec, falls back within 20 min. Ref: Sasaoka et al JBC 269:51 32621 1994. Use k3 = 0.1 based on this tau."
call /kinetics/EGFR/L.EGFR/SHC_phospho/notes LOAD \
"Rates from Okada et al JBC 270:35 pp 20737 1995 Km = 0.70 to 0.85 uM, Vmax = 4.4 to 5.0 pmol/min. Unfortunately the amount of enzyme is not known, the prep is only partially purified. Time course of phosph is max within 30 sec, falls back within 20 min. Ref: Sasaoka et al JBC 269:51 32621 1994. Use k3 = 0.1 based on this tau."
simundump kpool /kinetics/EGFR/EGF 0 0.0 0 0 0 0 0 0 6.3175e+05 4 \
/kinetics/geometry red 6 124 1 0
simundump text /kinetics/EGFR/EGF/notes 0 ""
call /kinetics/EGFR/EGF/notes LOAD \
""
simundump kpool /kinetics/EGFR/Internal_L.EGFR 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 6 126 1 0
simundump text /kinetics/EGFR/Internal_L.EGFR/notes 0 ""
call /kinetics/EGFR/Internal_L.EGFR/notes LOAD \
""
simundump kreac /kinetics/EGFR/act_EGFR 0 6.6482e-06 0.25 "" white 6 121 3 0
simundump text /kinetics/EGFR/act_EGFR/notes 0 \
"Affinity of EGFR for EGF is complex: depends on [EGFR]. We'll assume fixed [EGFR] and use exptal affinity ~20 nM (see Sherrill and Kyte Biochem 1996 35 5705-5718, Berkers et al JBC 266:2 922-927 1991, Sorokin et al JBC 269:13 9752-9759 1994). Tau =~2 min (Davis et al JBC 263:11 5373-5379 1988) or Berkers Kass = 6.2e5/M/sec, Kdiss=3.5e-4/sec. Sherrill and Kyte have Hill Coeff=1.7"
call /kinetics/EGFR/act_EGFR/notes LOAD \
"Affinity of EGFR for EGF is complex: depends on [EGFR]. We'll assume fixed [EGFR] and use exptal affinity ~20 nM (see Sherrill and Kyte Biochem 1996 35 5705-5718, Berkers et al JBC 266:2 922-927 1991, Sorokin et al JBC 269:13 9752-9759 1994). Tau =~2 min (Davis et al JBC 263:11 5373-5379 1988) or Berkers Kass = 6.2e5/M/sec, Kdiss=3.5e-4/sec. Sherrill and Kyte have Hill Coeff=1.7"
simundump kreac /kinetics/EGFR/Internalize 0 0.002 0.00033 "" white 6 123 3 0
simundump text /kinetics/EGFR/Internalize/notes 0 \
"See Helin and Beguinot JBC 266:13 1991 pg 8363-8368. In Fig 3 they have internalization tau about 10 min, equil at about 20% EGF available. So kf = 4x kb, and 1/(kf + kb) = 600 sec so kb = 1/3K = 3.3e-4, and kf = 1.33e-3. This doesn't take into account the unbound receptor, so we need to push the kf up a bit, to 0.002"
call /kinetics/EGFR/Internalize/notes LOAD \
"See Helin and Beguinot JBC 266:13 1991 pg 8363-8368. In Fig 3 they have internalization tau about 10 min, equil at about 20% EGF available. So kf = 4x kb, and 1/(kf + kb) = 600 sec so kb = 1/3K = 3.3e-4, and kf = 1.33e-3. This doesn't take into account the unbound receptor, so we need to push the kf up a bit, to 0.002"
simundump group /kinetics/Sos 0 7 black x 0 0 "" defaultfile defaultfile.g 0 \
0 0 0 20 0
simundump text /kinetics/Sos/notes 0 ""
call /kinetics/Sos/notes LOAD \
""
simundump kpool /kinetics/Sos/SHC_p.Sos.Grb2 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry brown 7 0 21 0
simundump text /kinetics/Sos/SHC_p.Sos.Grb2/notes 0 ""
call /kinetics/Sos/SHC_p.Sos.Grb2/notes LOAD \
""
simundump kenz /kinetics/Sos/SHC_p.Sos.Grb2/Sos.Ras_GEF 0 0 0 0 0 6.3175e+05 \
2.9766e-07 0.08 0.02 0 0 "" brown red "" 0 23 0
simundump text /kinetics/Sos/SHC_p.Sos.Grb2/Sos.Ras_GEF/notes 0 ""
call /kinetics/Sos/SHC_p.Sos.Grb2/Sos.Ras_GEF/notes LOAD \
""
simundump kpool /kinetics/Sos/SHC 0 0.0 0.50001 0.50001 3.1588e+05 3.1588e+05 \
0 0 6.3175e+05 0 /kinetics/geometry orange 7 2 21 0
simundump text /kinetics/Sos/SHC/notes 0 \
"There are 2 isoforms: 52 KDa and 46 KDa (See Okada et al JBC 270:35 pp 20737 1995). They are acted up on by the EGFR in very similar ways, and apparently both bind Grb2 similarly, so we'll bundle them together here. Sasaoka et al JBC 269:51 pp 32621 1994 show immunoprecs where it looks like there is at least as much Shc as Grb2. So we'll tentatively say there is 1 uM of Shc."
call /kinetics/Sos/SHC/notes LOAD \
"There are 2 isoforms: 52 KDa and 46 KDa (See Okada et al JBC 270:35 pp 20737 1995). They are acted up on by the EGFR in very similar ways, and apparently both bind Grb2 similarly, so we'll bundle them together here. Sasaoka et al JBC 269:51 pp 32621 1994 show immunoprecs where it looks like there is at least as much Shc as Grb2. So we'll tentatively say there is 1 uM of Shc."
simundump kpool /kinetics/Sos/SHC_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 7 4 21 0
simundump text /kinetics/Sos/SHC_p/notes 0 ""
call /kinetics/Sos/SHC_p/notes LOAD \
""
simundump kpool /kinetics/Sos/Sos_p.Grb2 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 7 6 21 0
simundump text /kinetics/Sos/Sos_p.Grb2/notes 0 ""
call /kinetics/Sos/Sos_p.Grb2/notes LOAD \
""
simundump kpool /kinetics/Sos/Grb2 0 0.0 1 1 6.3175e+05 6.3175e+05 0 0 \
6.3175e+05 0 /kinetics/geometry orange 7 8 21 0
simundump text /kinetics/Sos/Grb2/notes 0 \
"There is probably a lot of it in the cell: it is also known as Ash (abundant src homology protein I think). Also Waters et al JBC 271:30 18224 1996 say that only a small fraction of cellular Grb is precipitated out when SoS is precipitated. As most of the Sos seems to be associated with Grb2, it would seem like there is a lot of the latter. Say 1 uM. I haven't been able to find a decent....."
call /kinetics/Sos/Grb2/notes LOAD \
"There is probably a lot of it in the cell: it is also known as Ash (abundant src homology protein I think). Also Waters et al JBC 271:30 18224 1996 say that only a small fraction of cellular Grb is precipitated out when SoS is precipitated. As most of the Sos seems to be associated with Grb2, it would seem like there is a lot of the latter. Say 1 uM. I haven't been able to find a decent...."
simundump kpool /kinetics/Sos/Sos.Grb2 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry orange 7 10 21 0
simundump text /kinetics/Sos/Sos.Grb2/notes 0 ""
call /kinetics/Sos/Sos.Grb2/notes LOAD \
""
simundump kpool /kinetics/Sos/Sos_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 7 12 21 0
simundump text /kinetics/Sos/Sos_p/notes 0 ""
call /kinetics/Sos/Sos_p/notes LOAD \
""
simundump kpool /kinetics/Sos/Sos 0 0.0 0.1 0.1 63175 63175 0 0 6.3175e+05 0 \
/kinetics/geometry red 7 14 21 0
simundump text /kinetics/Sos/Sos/notes 0 \
"I have tried using low (0.02 uM) initial concs, but these give a very flat response to EGF stim although the overall activation of Ras is not too bad. I am reverting to 0.1 because we expect a sharp initial response, followed by a decline. Sep 17 1997: The transient activation curve looks better with [Sos] = 0.05. Apr 26 1998: Some error there, it is better where it was at 0.1"
call /kinetics/Sos/Sos/notes LOAD \
"I have tried using low (0.02 uM) initial concs, but these give a very flat response to EGF stim although the overall activation of Ras is not too bad. I am reverting to 0.1 because we expect a sharp initial response, followed by a decline. Sep 17 1997: The transient activation curve looks better with [Sos] = 0.05. Apr 26 1998: Some error there, it is better where it was at 0.1"
simundump kpool /kinetics/Sos/SHC_p_Grb2_clx 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry 42 7 16 21 0
simundump text /kinetics/Sos/SHC_p_Grb2_clx/notes 0 ""
call /kinetics/Sos/SHC_p_Grb2_clx/notes LOAD \
""
simundump kreac /kinetics/Sos/SHC_p_dephospho 0 0.0016667 0 "" white 7 1 23 0
simundump text /kinetics/Sos/SHC_p_dephospho/notes 0 \
"Time course of decline of phosph is 20 min. Part of this is the turnoff time of the EGFR itself. Lets assume a tau of 10 min for this dephosph. It may be wildly off."
call /kinetics/Sos/SHC_p_dephospho/notes LOAD \
"Time course of decline of phosph is 20 min. Part of this is the turnoff time of the EGFR itself. Lets assume a tau of 10 min for this dephosph. It may be wildly off."
simundump kreac /kinetics/Sos/SHC_bind_Sos.Grb2 0 7.9142e-07 0.1 "" white 7 3 \
23 0
simundump text /kinetics/Sos/SHC_bind_Sos.Grb2/notes 0 \
"Sasaoka et al JBC 269:51 pp 32621 1994, table on pg 32623 indicates that this pathway accounts for about 50% of the GEF activation. (88% - 39%). Error is large, about 20%. Fig 1 is most useful in constraining rates. Chook et al JBC 271:48 pp 30472, 1996 say that the Kd is 0.2 uM for Shc binding to EGFR. The Kd for Grb direct binding is 0.7, so we'll ignore it."
call /kinetics/Sos/SHC_bind_Sos.Grb2/notes LOAD \
"Sasaoka et al JBC 269:51 pp 32621 1994, table on pg 32623 indicates that this pathway accounts for about 50% of the GEF activation. (88% - 39%). Error is large, about 20%. Fig 1 is most useful in constraining rates. Chook et al JBC 271:48 pp 30472, 1996 say that the Kd is 0.2 uM for Shc binding to EGFR. The Kd for Grb direct binding is 0.7, so we'll ignore it."
simundump kreac /kinetics/Sos/Grb2_bind_Sos_p 0 3.9573e-08 0.0168 "" white 7 \
5 23 0
simundump text /kinetics/Sos/Grb2_bind_Sos_p/notes 0 \
"Same rates as Grb2_bind_Sos: Porfiri and McCormick JBC 271:10 pp 5871 1996 show that the binding is not affected by the phosph."
call /kinetics/Sos/Grb2_bind_Sos_p/notes LOAD \
"Same rates as Grb2_bind_Sos: Porfiri and McCormick JBC 271:10 pp 5871 1996 show that the binding is not affected by the phosph."
simundump kreac /kinetics/Sos/dephosph_Sos 0 0.001 0 "" white 7 7 23 0
simundump text /kinetics/Sos/dephosph_Sos/notes 0 \
"The only clue I have to these rates is from the time courses of the EGF activation, which is around 1 to 5 min. The dephosph would be expected to be of the same order, perhaps a bit longer. Lets use 0.002 which is about 8 min. Sep 17: The transient activation curve matches better with kf = 0.001"
call /kinetics/Sos/dephosph_Sos/notes LOAD \
"The only clue I have to these rates is from the time courses of the EGF activation, which is around 1 to 5 min. The dephosph would be expected to be of the same order, perhaps a bit longer. Lets use 0.002 which is about 8 min. Sep 17: The transient activation curve matches better with kf = 0.001"
simundump kreac /kinetics/Sos/Grb2_bind_Sos 0 3.9573e-08 0.0168 "" white 7 9 \
23 0
simundump text /kinetics/Sos/Grb2_bind_Sos/notes 0 \
"As there are 2 SH3 domains, this reaction could be 2nd order. I have a Kd of 22 uM from peptide binding (Lemmon et al JBC 269:50 pg 31653). However, Chook et al JBC 271:48 pg30472 say it is 0.4uM with purified proteins, so we believe them. They say it is 1:1 binding."
call /kinetics/Sos/Grb2_bind_Sos/notes LOAD \
"As there are 2 SH3 domains, this reaction could be 2nd order. I have a Kd of 22 uM from peptide binding (Lemmon et al JBC 269:50 pg 31653). However, Chook et al JBC 271:48 pg30472 say it is 0.4uM with purified proteins, so we believe them. They say it is 1:1 binding."
simundump kreac /kinetics/Sos/Grb2_bind_SHC 0 1.5829e-06 1 "" white 7 11 23 0
simundump text /kinetics/Sos/Grb2_bind_SHC/notes 0 ""
call /kinetics/Sos/Grb2_bind_SHC/notes LOAD \
""
simundump group /kinetics/PLC_g 0 8 black x 0 0 "" defaultfile defaultfile.g \
0 0 0 20 20 0
simundump text /kinetics/PLC_g/notes 0 ""
call /kinetics/PLC_g/notes LOAD \
""
simundump kpool /kinetics/PLC_g/PLC_g 0 0.0 0.82001 0.82001 5.1804e+05 \
5.1804e+05 0 0 6.3175e+05 0 /kinetics/geometry pink 8 20 21 0
simundump text /kinetics/PLC_g/PLC_g/notes 0 \
"Amount from Homma et al JBC 263:14 6592-6598 1988"
call /kinetics/PLC_g/PLC_g/notes LOAD \
"Amount from Homma et al JBC 263:14 6592-6598 1988"
simundump kpool /kinetics/PLC_g/PLC_g_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 8 22 21 0
simundump text /kinetics/PLC_g/PLC_g_p/notes 0 ""
call /kinetics/PLC_g/PLC_g_p/notes LOAD \
""
simundump kpool /kinetics/PLC_g/Ca.PLC_g 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 8 24 21 0
simundump text /kinetics/PLC_g/Ca.PLC_g/notes 0 ""
call /kinetics/PLC_g/Ca.PLC_g/notes LOAD \
""
simundump kenz /kinetics/PLC_g/Ca.PLC_g/PIP2_hydrolysis 0 0 0 0 0 6.3175e+05 \
1.1397e-06 56 14 0 1 "" pink red "" 24 23 0
simundump text /kinetics/PLC_g/Ca.PLC_g/PIP2_hydrolysis/notes 0 \
"Mainly Homma et al JBC 263:14 1988 pp 6592, but these parms are the Ca-stimulated form. It is not clear whether the enzyme is activated by tyrosine phosphorylation at this point or not. Wahl et al JBC 267:15 10447-10456 1992 say that the Ca_stim and phosph form has 7X higher affinity for substrate than control. This is close to Wahl's figure 7, which I am using as reference."
call /kinetics/PLC_g/Ca.PLC_g/PIP2_hydrolysis/notes LOAD \
"Mainly Homma et al JBC 263:14 1988 pp 6592, but these parms are the Ca-stimulated form. It is not clear whether the enzyme is activated by tyrosine phosphorylation at this point or not. Wahl et al JBC 267:15 10447-10456 1992 say that the Ca_stim and phosph form has 7X higher affinity for substrate than control. This is close to Wahl's figure 7, which I am using as reference."
simundump kpool /kinetics/PLC_g/Ca.PLC_g_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry pink 8 26 21 0
simundump text /kinetics/PLC_g/Ca.PLC_g_p/notes 0 ""
call /kinetics/PLC_g/Ca.PLC_g_p/notes LOAD \
""
simundump kenz /kinetics/PLC_g/Ca.PLC_g_p/PIP2_hydrolysis 0 0 0 0 0 \
6.3175e+05 2.2794e-05 228 57 0 1 "" pink red "" 26 23 0
simundump text /kinetics/PLC_g/Ca.PLC_g_p/PIP2_hydrolysis/notes 0 \
"Mainly Homma et al JBC 263:14 1988 pp 6592, but these parms are the Ca-stimulated form. It is not clear whether the enzyme is activated by tyrosine phosphorylation at this point or not. Wahl et al JBC 267:15 10447-10456 1992 say that this has 7X higher affinity for substrate than control."
call /kinetics/PLC_g/Ca.PLC_g_p/PIP2_hydrolysis/notes LOAD \
"Mainly Homma et al JBC 263:14 1988 pp 6592, but these parms are the Ca-stimulated form. It is not clear whether the enzyme is activated by tyrosine phosphorylation at this point or not. Wahl et al JBC 267:15 10447-10456 1992 say that this has 7X higher affinity for substrate than control."
simundump kpool /kinetics/PLC_g/PLCg_basal 0 0.0 0.00070002 0.00070002 442.24 \
442.24 0 0 6.3175e+05 0 /kinetics/geometry 33 8 28 21 0
simundump text /kinetics/PLC_g/PLCg_basal/notes 0 ""
call /kinetics/PLC_g/PLCg_basal/notes LOAD \
""
simundump kenz /kinetics/PLC_g/PLCg_basal/PLC_g_phospho 0 0 0 0 0 6.3175e+05 \
1.2714e-05 2 0.5 0 0 "" red 33 "" 28 23 0
simundump text /kinetics/PLC_g/PLCg_basal/PLC_g_phospho/notes 0 ""
call /kinetics/PLC_g/PLCg_basal/PLC_g_phospho/notes LOAD \
""
simundump kreac /kinetics/PLC_g/Ca_act_PLC_g 0 0.00028492 10 "" white 8 21 23 \
0
simundump text /kinetics/PLC_g/Ca_act_PLC_g/notes 0 \
"Nice curves from Homma et al JBC 263:14 6592-6598 1988 Fig 5c. The activity falls above 10 uM, but that is too high to reach physiologically anyway, so we'll ignore the higher pts and match the lower ones only. Half-max at 1 uM. But Wahl et al JBC 267:15 10447-10456 1992 have half-max at 56 nM which is what I'll use."
call /kinetics/PLC_g/Ca_act_PLC_g/notes LOAD \
"Nice curves from Homma et al JBC 263:14 6592-6598 1988 Fig 5c. The activity falls above 10 uM, but that is too high to reach physiologically anyway, so we'll ignore the higher pts and match the lower ones only. Half-max at 1 uM. But Wahl et al JBC 267:15 10447-10456 1992 have half-max at 56 nM which is what I'll use."
simundump kreac /kinetics/PLC_g/Ca_act_PLC_g_p 0 1.8995e-05 10 "" white 8 23 \
23 0
simundump text /kinetics/PLC_g/Ca_act_PLC_g_p/notes 0 \
"Again, we refer to Homma et al and Wahl et al, for preference using Wahl. Half-Max of the phosph form is at 316 nM. Use kb of 10 as this is likely to be pretty fast. Did some curve comparisons, and instead of 316 nM giving a kf of 5.27e-5, we will use 8e-5 for kf. 16 Sep 97. As we are now phosphorylating the Ca-bound form, equils have shifted. kf should now be 2e-5 to match the curves."
call /kinetics/PLC_g/Ca_act_PLC_g_p/notes LOAD \
"Again, we refer to Homma et al and Wahl et al, for preference using Wahl. Half-Max of the phosph form is at 316 nM. Use kb of 10 as this is likely to be pretty fast. Did some curve comparisons, and instead of 316 nM giving a kf of 5.27e-5, we will use 8e-5 for kf. 16 Sep 97. As we are now phosphorylating the Ca-bound form, equils have shifted. kf should now be 2e-5 to match the curves."
simundump kreac /kinetics/PLC_g/dephosph_PLC_g 0 0.05 0 "" white 8 25 23 0
simundump text /kinetics/PLC_g/dephosph_PLC_g/notes 0 ""
call /kinetics/PLC_g/dephosph_PLC_g/notes LOAD \
""
simundump kreac /kinetics/PLC_g/PLC_g_p_dephospho 0 0.07 0 "" white 8 27 23 0
simundump text /kinetics/PLC_g/PLC_g_p_dephospho/notes 0 ""
call /kinetics/PLC_g/PLC_g_p_dephospho/notes LOAD \
""
simundump group /kinetics/CaMKII 0 9 black x 0 0 "" defaultfile defaultfile.g \
0 0 0 40 20 0
simundump text /kinetics/CaMKII/notes 0 ""
call /kinetics/CaMKII/notes LOAD \
""
simundump kpool /kinetics/CaMKII/CaMKII 0 0.0 70.001 70.001 4.4223e+07 \
4.4223e+07 0 0 6.3175e+05 0 /kinetics/geometry palegreen 9 40 21 0
simundump text /kinetics/CaMKII/CaMKII/notes 0 \
"Huge conc of CaMKII. In PSD it is 20-40% of protein, so we assume it is around 2.5% of protein in spine as a whole. This level is so high it is unlikely to matter much if we are off a bit. This comes to about 70 uM."
call /kinetics/CaMKII/CaMKII/notes LOAD \
"Huge conc of CaMKII. In PSD it is 20-40% of protein, so we assume it is around 2.5% of protein in spine as a whole. This level is so high it is unlikely to matter much if we are off a bit. This comes to about 70 uM."
simundump kpool /kinetics/CaMKII/CaMKII_CaM 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry palegreen 9 42 21 0
simundump text /kinetics/CaMKII/CaMKII_CaM/notes 0 ""
call /kinetics/CaMKII/CaMKII_CaM/notes LOAD \
""
simundump kpool /kinetics/CaMKII/CaMKII_thr286_p_CaM 0 0.0 0 0 0 0 0 0 \
6.3175e+05 0 /kinetics/geometry palegreen 9 44 21 0
simundump text /kinetics/CaMKII/CaMKII_thr286_p_CaM/notes 0 \
"From Hanson and Schulman, the thr286 is responsible for autonomous activation of CaMKII."
call /kinetics/CaMKII/CaMKII_thr286_p_CaM/notes LOAD \
"From Hanson and Schulman, the thr286 is responsible for autonomous activation of CaMKII."
simundump kpool /kinetics/CaMKII/CaMKII_p_p_p 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry cyan 9 46 21 0
simundump text /kinetics/CaMKII/CaMKII_p_p_p/notes 0 \
"From Hanson and Schulman, the CaMKII does a lot of autophosphorylation just after the CaM is released. This prevents further CaM binding and renders the enzyme quite independent of Ca."
call /kinetics/CaMKII/CaMKII_p_p_p/notes LOAD \
"From Hanson and Schulman, the CaMKII does a lot of autophosphorylation just after the CaM is released. This prevents further CaM binding and renders the enzyme quite independent of Ca."
simundump kpool /kinetics/CaMKII/CaMKII_thr286 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry red 9 48 21 0
simundump text /kinetics/CaMKII/CaMKII_thr286/notes 0 \
"I am not sure if we need to endow this one with a lot of enzs. It is likely to be a short-lived intermediate, since it will be phosphorylated further as soon as the CAM falls off."
call /kinetics/CaMKII/CaMKII_thr286/notes LOAD \
"I am not sure if we need to endow this one with a lot of enzs. It is likely to be a short-lived intermediate, since it will be phosphorylated further as soon as the CAM falls off."
simundump kpool /kinetics/CaMKII/CaMKII_thr306 0 0.0 0 0 0 0 0 0 6.3175e+05 0 \
/kinetics/geometry palegreen 9 50 21 0
simundump text /kinetics/CaMKII/CaMKII_thr306/notes 0 \
"This forms due to basal autophosphorylation, but I think it has to be considered as a pathway even if some CaM is floating around. In either case it will tend to block further binding of CaM, and will not display any enzyme activity. See Hanson and Schulman JBC 267:24 pp17216-17224 1992"
call /kinetics/CaMKII/CaMKII_thr306/notes LOAD \
"This forms due to basal autophosphorylation, but I think it has to be considered as a pathway even if some CaM is floating around. In either case it will tend to block further binding of CaM, and will not display any enzyme activity. See Hanson and Schulman JBC 267:24 pp17216-17224 1992"
simundump kpool /kinetics/CaMKII/tot_CaM_CaMKII 0 0.0 0 0 0 0 0 0 6.3175e+05 \
0 /kinetics/geometry green 9 43 31 0
simundump text /kinetics/CaMKII/tot_CaM_CaMKII/notes 0 ""
call /kinetics/CaMKII/tot_CaM_CaMKII/notes LOAD \
""
simundump kenz /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305 0 0 0 0 0 \
6.3175e+05 3.9701e-07 24 6 0 0 "" green red "" 43 33 0
simundump text /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305/notes 0 \
"Rates from autocamtide phosphorylation, from Hanson and Schulman JBC 267:24 17216-17224 1992. Jan 1 1998: Speed up 12x to match fig 5."
call /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305/notes LOAD \
"Rates from autocamtide phosphorylation, from Hanson and Schulman JBC 267:24 17216-17224 1992. Jan 1 1998: Speed up 12x to match fig 5."
simundump kenz /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286 0 0 0 0 0 \
6.3175e+05 3.3084e-08 2 0.5 0 0 "" green red "" 43 34 0
simundump text /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286/notes 0 ""
call /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286/notes LOAD \
""
simundump kpool /kinetics/CaMKII/tot_autonomous_CaMKII 0 0.0 0 0 0 0 0 0 \
6.3175e+05 0 /kinetics/geometry green 9 51 31 0
simundump text /kinetics/CaMKII/tot_autonomous_CaMKII/notes 0 ""
call /kinetics/CaMKII/tot_autonomous_CaMKII/notes LOAD \
""
simundump kenz /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305 0 0 0 0 0 \
6.3175e+05 2.5771e-07 24 6 0 0 "" green red "" 51 33 0
simundump text /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305/notes 0 \
"See Hanson and Schulman again, for afterburst rates of phosph."
call /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305/notes LOAD \
"See Hanson and Schulman again, for afterburst rates of phosph."
simundump kenz /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286 0 0 0 0 0 \
6.3175e+05 2.1477e-08 2 0.5 0 0 "" green red "" 51 34 0
simundump text /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286/notes 0 ""
call /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286/notes LOAD \
""
simundump kreac /kinetics/CaMKII/CaMKII_bind_CaM 0 7.9145e-05 5 "" white 9 41 \
23 0
simundump text /kinetics/CaMKII/CaMKII_bind_CaM/notes 0 \
"This is tricky. There is some cooperativity here arising from interactions between the subunits of the CAMKII holoenzyme. However, the stoichiometry is 1. Kb/Kf = 6e4 #/cell. Rate is fast (see Hanson et al Neuron 12 943-956 1994) so lets say kb = 10. This gives kf = 1.6667e-4 H&S AnnRev Biochem 92 give tau for dissoc as 0.2 sec at low Ca, 0.4 at high. Low Ca = 100 nM = physiol."
call /kinetics/CaMKII/CaMKII_bind_CaM/notes LOAD \
"This is tricky. There is some cooperativity here arising from interactions between the subunits of the CAMKII holoenzyme. However, the stoichiometry is 1. Kb/Kf = 6e4 #/cell. Rate is fast (see Hanson et al Neuron 12 943-956 1994) so lets say kb = 10. This gives kf = 1.6667e-4 H&S AnnRev Biochem 92 give tau for dissoc as 0.2 sec at low Ca, 0.4 at high. Low Ca = 100 nM = physiol."