-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathLoadParameters.m
828 lines (725 loc) · 9.15 KB
/
LoadParameters.m
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
function modelparameters=LoadParameters(cond)
%% Signaling parameters
p_1_1 = 1;
p_1_2 = 2;
p_2_1 = 1;
p_2_2 = 1.75;
p_3_1 = 1;
p_3_2 = 0.5;
p_4_1 = 1;
p_4_2 = 0.06;
p_5_1 = 1;
p_5_2 = 0.003;
p_6_1 = 0.014999625;
p_6_2 = 0.006;
p_7_1 = 0.41482;
p_7_2 = 45;
p_8_1 = 80.357;
p_8_2 = 180;
p_9_1 = 0.096786;
p_9_2 = 10.5;
p_10_1 = 18.75;
p_10_2 = 42.001;
p_11_1 = 6;
p_12_1 = 0.01;
p_13_1 = 0.0001;
p_14_1 = 0.0013;
p_14_2 = 1;
p_15_1 = 90.6537658442;
p_15_2 = 2.466;
p_16_1 = 60;
p_16_2 = 0.01;
p_17_1 = 3;
p_17_2 = 1;
p_18_1 = 3600; % 36 for high Ca baseline condition
p_18_2 = 144;
p_19_1 = 1;
p_19_2 = 6;
p_20_1 = 0.465;
p_20_2 = 10;
p_21_1 = 72;
p_21_2 = 72;
p_22_1 = 3.6;
p_22_2 = 10;
p_23_1 = 0.6;
p_23_2 = 0.5;
p_24_1 = 30;
p_24_2 = 1;
p_25_1 = 288;
p_26_1 = 8;
p_27_1 = 500;
p_28_1 = 4.89E-04;
p_28_2 = 1;
p_29_1 = 4.89E-04;
p_29_2 = 1;
p_30_1 = 4.89E-04;
p_30_2 = 1;
p_31_1 = 4.89E-04;
p_31_2 = 1;
p_32_1 = 4.89E-04;
p_32_2 = 1;
p_33_1 = 4.89E-04;
p_33_2 = 1;
p_34_1 = 4.89E-04;
p_34_2 = 1;
p_35_1 = 4.89E-04;
p_35_2 = 1;
p_36_1 = 0.11059;
p_36_2 = 1;
p_37_1 = 25;
p_38_1 = 1800;
p_38_2 = 360;
p_39_1 = 50;
p_39_2 = 0.1;
p_40_1 = 19.231;
p_40_2 = 1;
p_41_1 = 1000.2;
p_41_2 = 0.1;
p_42_1 = 50;
p_42_2 = 5;
p_43_1 = 0.3431978192;
p_43_2 = 1.4;
p_44_1 = 0.35;
p_45_1 = 0.003;
p_46_1 = 0.35;
p_47_1 = 0.35;
p_48_1 = 0.3431978192;
p_48_2 = 1.4;
p_49_1 = 0.3431978192;
p_49_2 = 1.4;
p_50_1 = 6;
p_51_1 = 6;
p_52_1 = 599995.2000384;
p_52_2 = 2;
p_53_1 = 924009.46185689;
p_53_2 = 2;
p_54_1 = 0.35;
p_55_1 = 7199942.4004608;
p_55_2 = 24;
p_56_1 = 11088113.542283;
p_56_2 = 24;
p_57_1 = 0.3431978192;
p_57_2 = 1.4;
p_58_1 = 0.35;
p_59_1 = 0.3431978192;
p_59_2 = 1.4;
p_60_1 = 0.5;
p_61_1 = 0.5;
p_62_1 = 25;
p_63_1 = 144;
p_64_1 = 0.15;
p_65_1 = 0.0079998;
p_65_2 = 8.6348;
p_66_1 = 0.0006;
p_66_2 = 0.1;
p_67_1 = 0.003;
p_68_1 = 2.5;
p_69_1 = 10;
p_70_1 = 75;
p_71_1 = 160;
p_72_1 = 1;
p_72_2 = 1.6;
p_73_1 = 1;
p_73_2 = 0.2;
p_74_1 = 75;
p_74_2 = 25.5;
p_75_1 = 90.6537658442;
p_75_2 = 3.0825;
p_76_1 = 0.0133;
p_77_1 = 1.6667;
p_78_1 = 25.2;
p_78_2 = 1;
p_79_1 = 2.52;
p_79_2 = 1;
p_80_1 = 0.6;
p_80_2 = 1;
p_81_1 = 0.006;
p_81_2 = 0.0001;
p_82_1 = 0.35719;
p_82_2 = 45.72;
p_83_1 = 0.11;
p_84_1 = 117.33;
p_85_1 = 39.2;
p_86_1 = 35;
p_86_2 = 0.336;
p_87_1 = 55;
p_88_1 = 1.035;
p_88_2 = 157.32;
p_89_1 = 2;
p_90_1 = 0.00078582;
p_90_2 = 11.43;
p_91_1 = 0.00207;
p_91_2 = 39.33;
p_92_1 = 733.3216665208;
p_92_2 = 469.32;
p_93_1 = 52.972972973;
p_93_2 = 156.8;
p_94_1 = 20;
p_95_1 = 1.9999600008;
p_95_2 = 8;
p_96_1 = 1.1;
p_97_1 = 10.667;
p_98_1 = 11.9570416105;
p_98_2 = 220;
p_99_1 = 4.125;
p_100_1 = 4.2666;
p_100_2 = 1.024;
p_101_1 = 1.0311468853;
p_101_2 = 16.5;
p_102_1 = 4.9997500125;
p_102_2 = 80;
p_103_1 = 0.09;
p_104_1 = 0.003;
p_105_1 = 1;
p_105_2 = 0.001;
p_106_1 = 0.013;
p_107_1 = 10.8839285714;
p_107_2 = 12.19;%15.2375
p_108_1 = 10.8839285714;
p_108_2 = 12.19;
p_109_1 = 2.7209642857;
p_109_2 = 3.04748;
p_110_1 = 1.8333333333;
p_110_2 = 4.4;
p_111_1 = 2.9629962944;
p_111_2 = 42.6669333;
p_112_1 = 0.0012;
p_113_1 = 1;
p_113_2 = 0.001;
p_114_1 = 4.675;
p_115_1 = 0.0060335;
p_115_2 = 0.070136;
p_116_1 = 0.256;
p_117_1 = 31.2;
p_117_2 = 2.496;
p_118_1 = 0.056;
p_118_2 = 0.624;
p_119_1 = 0.05;
p_119_2 = 1;
p_120_1 = 1.5625;
p_120_2 = 2;%2.5
p_121_1 = 1;
p_121_2 = 2;%2.5
p_122_1 = 0.6666666667;
p_122_2 = 16;%20
p_123_1 = 0.5;
p_124_1 = 0.5;
p_125_1 = 4;
p_126_1 = 1.5625;
p_126_2 = 2;
p_127_1 = 1;
p_127_2 = 2;
p_128_1 = 0.5;
p_129_1 = 0.5;
p_130_1 = 1;
p_130_2 = 0.001;
p_131_1 = 0.014999625;
p_131_2 = 0.0048;
p_132_1 = 0.5624859379;
p_132_2 = 0.36;
p_133_1 = 3.0475;
p_134_1 = 3.0475;
p_135_1 = 0.76187;
p_136_1 = 0.001;
p_136_2 = 1;
p_137_1 = 179.8076923077;
p_137_2 = 18.7;
p_138_1 = 0.15;
p_139_1 = 0.003;
p_140_1 = 53.75;
p_140_2 = 17.2;
p_141_1 = 0.013;
p_142_1 = 1;
p_142_2 = 0.001;
p_143_1 = 0.000335;
p_144_1 = 0.001;
p_144_2 = 1;
p_145_1 = 0.2321;
p_146_1 = 0.01355;
p_147_1 = 0.26875;
p_147_2 = 4.3;
p_148_1 = 0.001;
p_148_2 = 1;
p_149_1 = 0.0132;
p_150_1 = 1;
p_150_2 = 45;
p_151_1 = 1;
p_152_1 = 1.92;
p_153_1 = 1.75E-005;
p_154_1 = 1;
p_154_2 = 0.001;
p_155_1 = 0.1746268657;
p_155_2 = 0.936;
p_156_1 = 0.034047619;
p_156_2 = 0.2288;
p_157_1 = 0.001;
p_157_2 = 1;
p_158_1 = 29.0125;
p_158_2 = 0.9284;
p_159_1 = 1;
p_159_2 = 16;
p_160_1 = 0.0038;
p_160_2 = 2.376;
p_161_1 = 1;
p_161_2 = 0.001;
p_162_1 = 1;
p_162_2 = 19;
p_163_1 = 2.7216666667;
p_163_2 = 6.532;
p_164_1 = 5.1;
p_165_1 = 0.0012;
p_165_2 = 1.26;
p_166_1 = 0.0108695652;
p_166_2 = 0.25;
p_167_1 = 0.001;
p_167_2 = 1;
p_168_1 = 220;
p_168_2 = 0.0528;
p_169_1 = 0.0016666667;
p_169_2 = 0.05;
p_170_1 = 1923;
p_170_2 = 10;
p_171_1 = 0.00012;
p_171_2 = 0.1;
p_172_1 = 1;
p_173_1 = 50;
p_174_1 = 0.0012;
p_174_2 = 0.1;
p_175_1 = 1.2705;
p_176_1 = 3.5026;
p_177_1 = 0.6;
p_177_2 = 0.5;
p_178_1 = 1;
p_179_1 = 0.1;
p_180_1 = 0.018;
p_180_2 = 2;
p_181_1 = 2;
p_182_1 = 0.2;
p_183_1 = 0.6250312516;
p_183_2 = 12.5;
p_184_1 = 1.6667;
p_185_1 = 2.5200342725;
p_185_2 = 50;
p_186_1 = 160.003200064;
p_186_2 = 800;
p_187_1 = 75;
p_187_2 = 375;
p_188_1 = 0.234;
p_189_1 = 0.0572;
p_190_1 = 0.0027;
p_190_2 = 2.5;
p_191_1 = 0.6165;
p_192_1 = 1.633;
p_193_1 = 3.00E-004;
p_194_1 = 1923;
p_194_2 = 10;
p_195_1 = 25;
p_196_1 = 586.65;
p_197_1 = 196;
p_198_1 = 15.2375;
p_199_1 = 15.2375;
p_200_1 = 3.80935;
p_201_1 = 53.3339333;
p_202_1 = 1.17;
p_203_1 = 0.286;
p_204_1 = 8.165;
p_205_1 = 10;
p_206_1 = 20.625;
p_207_1 = 275;
p_208_1 = 100;
p_209_1 = 5.5;
p_210_1 = 0.45;
p_211_1 = 23.375;
p_212_1 = 1.1605;
p_213_1 = 0.066;
p_214_1 = 0.2;
p_215_1 = 0.04;
p_216_1 = 20.4;
p_217_1 = 1.75;
p_218_1 = 1.75;
p_219_1 = 1.75;
p_220_1 = 1.75;
p_221_1 = 1.75;
p_222_1 = 2.5;
p_223_1 = 30;
p_224_1 = 2.5;
p_225_1 = 2.5;
p_226_1 = 2.5;
p_227_1 = 2.5;
p_228_1 = 2.5;
p_229_1 = 30;
p_230_1 = 0.006;
p_231_1 = 20;
g.na = 3*cond(2);% Na, g,max
g.kdr = 0.9*cond(3);% Kdr, g,max
g.ka = 0.15*cond(4);% Ka, g,max
g.kahp = 0.15*cond(5); % Kahp, g,max
g.cal = 0.0015*cond(6); % CaL, g,max
g.leak = 0.01; % gleak, g,max
%phosMK1.max = 1.35;
%phosPKC1.max = 0.5200;
phosMK1.k_p = 1; % 0 for CaMKII blocking condition
phosPKC1.k_p = 10; % 0 for PKC blocking condition
vmax = 1.2e-1;
%% initializing all params
modelparameters = [
p_1_1
p_1_2
p_2_1
p_2_2
p_3_1
p_3_2
p_4_1
p_4_2
p_5_1
p_5_2
p_6_1
p_6_2
p_7_1
p_7_2
p_8_1
p_8_2
p_9_1
p_9_2
p_10_1
p_10_2
p_11_1
p_12_1
p_13_1
p_14_1
p_14_2
p_15_1
p_15_2
p_16_1
p_16_2
p_17_1
p_17_2
p_18_1
p_18_2
p_19_1
p_19_2
p_20_1
p_20_2
p_21_1
p_21_2
p_22_1
p_22_2
p_23_1
p_23_2
p_24_1
p_24_2
p_25_1
p_26_1
p_27_1
p_28_1
p_28_2
p_29_1
p_29_2
p_30_1
p_30_2
p_31_1
p_31_2
p_32_1
p_32_2
p_33_1
p_33_2
p_34_1
p_34_2
p_35_1
p_35_2
p_36_1
p_36_2
p_37_1
p_38_1
p_38_2
p_39_1
p_39_2
p_40_1
p_40_2
p_41_1
p_41_2
p_42_1
p_42_2
p_43_1
p_43_2
p_44_1
p_45_1
p_46_1
p_47_1
p_48_1
p_48_2
p_49_1
p_49_2
p_50_1
p_51_1
p_52_1
p_52_2
p_53_1
p_53_2
p_54_1
p_55_1
p_55_2
p_56_1
p_56_2
p_57_1
p_57_2
p_58_1
p_59_1
p_59_2
p_60_1
p_61_1
p_62_1
p_63_1
p_64_1
p_65_1
p_65_2
p_66_1
p_66_2
p_67_1
p_68_1
p_69_1
p_70_1
p_71_1
p_72_1
p_72_2
p_73_1
p_73_2
p_74_1
p_74_2
p_75_1
p_75_2
p_76_1
p_77_1
p_78_1
p_78_2
p_79_1
p_79_2
p_80_1
p_80_2
p_81_1
p_81_2
p_82_1
p_82_2
p_83_1
p_84_1
p_85_1
p_86_1
p_86_2
p_87_1
p_88_1
p_88_2
p_89_1
p_90_1
p_90_2
p_91_1
p_91_2
p_92_1
p_92_2
p_93_1
p_93_2
p_94_1
p_95_1
p_95_2
p_96_1
p_97_1
p_98_1
p_98_2
p_99_1
p_100_1
p_100_2
p_101_1
p_101_2
p_102_1
p_102_2
p_103_1
p_104_1
p_105_1
p_105_2
p_106_1
p_107_1
p_107_2
p_108_1
p_108_2
p_109_1
p_109_2
p_110_1
p_110_2
p_111_1
p_111_2
p_112_1
p_113_1
p_113_2
p_114_1
p_115_1
p_115_2
p_116_1
p_117_1
p_117_2
p_118_1
p_118_2
p_119_1
p_119_2
p_120_1
p_120_2
p_121_1
p_121_2
p_122_1
p_122_2
p_123_1
p_124_1
p_125_1
p_126_1
p_126_2
p_127_1
p_127_2
p_128_1
p_129_1
p_130_1
p_130_2
p_131_1
p_131_2
p_132_1
p_132_2
p_133_1
p_134_1
p_135_1
p_136_1
p_136_2
p_137_1
p_137_2
p_138_1
p_139_1
p_140_1
p_140_2
p_141_1
p_142_1
p_142_2
p_143_1
p_144_1
p_144_2
p_145_1
p_146_1
p_147_1
p_147_2
p_148_1
p_148_2
p_149_1
p_150_1
p_150_2
p_151_1
p_152_1
p_153_1
p_154_1
p_154_2
p_155_1
p_155_2
p_156_1
p_156_2
p_157_1
p_157_2
p_158_1
p_158_2
p_159_1
p_159_2
p_160_1
p_160_2
p_161_1
p_161_2
p_162_1
p_162_2
p_163_1
p_163_2
p_164_1
p_165_1
p_165_2
p_166_1
p_166_2
p_167_1
p_167_2
p_168_1
p_168_2
p_169_1
p_169_2
p_170_1
p_170_2
p_171_1
p_171_2
p_172_1
p_173_1
p_174_1
p_174_2
p_175_1
p_176_1
p_177_1
p_177_2
p_178_1
p_179_1
p_180_1
p_180_2
p_181_1
p_182_1
p_183_1
p_183_2
p_184_1
p_185_1
p_185_2
p_186_1
p_186_2
p_187_1
p_187_2
p_188_1
p_189_1
p_190_1
p_190_2
p_191_1
p_192_1
p_193_1
p_194_1
p_194_2
p_195_1
p_196_1
p_197_1
p_198_1
p_199_1
p_200_1
p_201_1
p_202_1
p_203_1
p_204_1
p_205_1
p_206_1
p_207_1
p_208_1
p_209_1
p_210_1
p_211_1
p_212_1
p_213_1
p_214_1
p_215_1
p_216_1
p_217_1
p_218_1
p_219_1
p_220_1
p_221_1
p_222_1
p_223_1
p_224_1
p_225_1
p_226_1
p_227_1
p_228_1
p_229_1
p_230_1
p_231_1
g.na
g.kdr
g.ka
g.kahp
g.cal
g.leak
phosMK1.k_p
phosPKC1.k_p
vmax
];