-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathai_ld.txt
9078 lines (9078 loc) · 318 KB
/
ai_ld.txt
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
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.8116883]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 3, 'min_samples_split': 3, 'n_estimators': 10}
[[0.6348709 0.3651291 ]
[0.65538214 0.34461786]
[0.61376532 0.38623468]
...
[0.81211507 0.18788493]
[0.43033497 0.56966503]
[0.79319501 0.20680499]]
########################
hmmScoreWT 0.1326346081955458
hmmScoreMUT 0.07978838548768871
Phosphomimic 0.03796691887414842
ncontacts 0.07460360241153044
nresidues 0.022261858322754687
burr 0.056422702323868354
acc 0.012710931124527338
ChargesMUT 0.06922836188841949
p_0_pfam 0.01193794071367812
ub_0_pfam 0.01125930610121326
p_1_pfam 0.010202162504343984
ub_1_pfam 0.02170071321209211
K_WT 0.0597531094384513
S_WT 0.01866547376975793
V_WT 0.012094037592429628
A_MUT 0.027516662706600043
H_MUT 0.011799336943382203
R_MUT 0.04141281836911343
allHomologs 0.011602643538301646
exclParalogs 0.018910739907638818
specParalogs 0.05294633414941931
orthologs 0.02923952046967753
bpso 0.028487869707095215
bpsh 0.09749156744601353
Feature Importance
0 hmmScoreWT 0.132635
23 bpsh 0.097492
1 hmmScoreMUT 0.079788
3 ncontacts 0.074604
7 ChargesMUT 0.069228
12 K_WT 0.059753
5 burr 0.056423
20 specParalogs 0.052946
17 R_MUT 0.041413
2 Phosphomimic 0.037967
21 orthologs 0.029240
22 bpso 0.028488
15 A_MUT 0.027517
4 nresidues 0.022262
11 ub_1_pfam 0.021701
19 exclParalogs 0.018911
13 S_WT 0.018665
6 acc 0.012711
14 V_WT 0.012094
8 p_0_pfam 0.011938
16 H_MUT 0.011799
18 allHomologs 0.011603
9 ub_0_pfam 0.011259
10 p_1_pfam 0.010202
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.665 resistance
FGFR2/L617M/train 0.568 resistance
FGFR2/V564F/train 0.7 resistance
FGFR2/K659M/train 0.448 resistance
FGFR2/L617V/train 0.379 resistance
FGFR2/M537I/train 0.68 resistance
ALK/L1196M/train 0.628 resistance
ALK/G1202R/train 0.309 resistance
ALK/I1171N/train 0.432 resistance
ALK/C1156Y/train 0.606 resistance
ALK/S1206F/train 0.605 resistance
MET/Y1230H/train 0.578 resistance
MET/D1228N/train 0.551 resistance
MET/Y1230C/train 0.395 resistance
MET/D1228V/train 0.575 resistance
MET/Y1230S/train 0.464 resistance
MET/D1228H/train 0.592 resistance
FLT3/F691L/train 0.628 resistance
FLT3/D835I/train 0.475 resistance
EGFR/L792H/train 0.47 resistance
EGFR/G796R/train 0.429 resistance
EGFR/L792R/train 0.391 resistance
EGFR/C797G/train 0.45 resistance
EGFR/L792F/train 0.597 resistance
EGFR/L718Q/train 0.392 resistance
EGFR/C797S/train 0.665 resistance
EGFR/G796D/train 0.437 resistance
EGFR/L718V/train 0.479 resistance
EGFR/G719A/train 0.427 resistance
EGFR/L747P/train 0.302 resistance
EGFR/L792V/train 0.552 resistance
EGFR/G796S/train 0.471 resistance
EGFR/L792P/train 0.406 resistance
EGFR/D761Y/train 0.563 resistance
FGFR4/V550L/train 0.7 resistance
FGFR4/V550M/train 0.7 resistance
ABL1/T315I/train 0.456 resistance
ABL1/G250E/train 0.711 resistance
ABL1/Y253H/train 0.514 resistance
ABL1/E255K/train 0.57 resistance
ABL1/F317L/train 0.57 resistance
ABL1/D276G/train 0.577 resistance
ABL1/F311I/train 0.57 resistance
BRAF/L505H/train 0.517 resistance
KIT/D816E/train 0.749 resistance
KIT/A829P/train 0.642 resistance
KIT/T670I/train 0.62 resistance
KIT/V654A/train 0.394 resistance
KIT/Y823D/train 0.443 resistance
KIT/D820Y/train 0.461 resistance
KIT/N822Y/train 0.723 resistance
KIT/D820G/train 0.647 resistance
KIT/D820E/train 0.717 resistance
KIT/K642E/train 0.756 resistance
KIT/C809G/train 0.52 resistance
KIT/N655T/train 0.606 resistance
KIT/D816G/train 0.551 resistance
ERBB2/L869R/train 0.489 resistance
BTK/C481R/train 0.514 resistance
BTK/C481S/train 0.68 resistance
BTK/C481Y/train 0.56 resistance
MAP2K1/P124L/train 0.422 resistance
MAP2K1/G128V/train 0.488 resistance
MAP2K1/C121S/train 0.542 resistance
MAP2K1/G128D/train 0.646 resistance
MAP2K1/P124S/train 0.448 resistance
MAP2K1/E203K/train 0.57 resistance
MAP2K2/N126D/train 0.734 resistance
MAP2K2/C125S/train 0.542 resistance
MAP2K2/E207K/train 0.57 resistance
KDR/L840F/train 0.463 resistance
########################
MAP2K3/A84T/test 0.665 A
PIM1/S97N/test 0.702 A
########################
CDK5/F82I/test 0.445 TBD
MAP2K1/V211D/test 0.328 TBD
########################
CHEK2/K373E/test 0.531 Inconclusive
MAP2K3/R94L/test 0.608 Inconclusive
MAP2K3/R96W/test 0.451 Inconclusive
MAP2K3/L215W/test 0.365 Inconclusive
MAP2K3/T222M/test 0.408 Inconclusive
PIM1/Q127E/test 0.662 Inconclusive
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.82050484]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 3, 'min_samples_split': 3, 'n_estimators': 50}
[[0.65490436 0.34509564]
[0.64932213 0.35067787]
[0.62678829 0.37321171]
...
[0.77990976 0.22009024]
[0.34691614 0.65308386]
[0.78479466 0.21520534]]
########################
hmmScoreWT 0.12810998510144475
hmmScoreMUT 0.08736397062905636
Phosphomimic 0.024751942734314712
Acetylmimic 0.021508216770495924
ncontacts 0.0309386586567444
nresidues 0.01442464155509491
phi_psi 0.011992762256382839
burr 0.03672838635895645
acc 0.023613003637461265
ChargesMUT 0.03467883798583897
m1_0_pfam 0.012293825082172593
ub_1_pfam 0.012936028108680402
K_WT 0.03632808592441564
A_MUT 0.0433903033660689
E_MUT 0.029925541452835165
R_MUT 0.017043429800486815
allHomologs 0.06716360693035313
exclParalogs 0.04050339613294779
specParalogs 0.042517191020995576
orthologs 0.028661378938028143
bpso 0.05358737948856675
bpsh 0.08841787609717368
Feature Importance
0 hmmScoreWT 0.128110
21 bpsh 0.088418
1 hmmScoreMUT 0.087364
16 allHomologs 0.067164
20 bpso 0.053587
13 A_MUT 0.043390
18 specParalogs 0.042517
17 exclParalogs 0.040503
7 burr 0.036728
12 K_WT 0.036328
9 ChargesMUT 0.034679
4 ncontacts 0.030939
14 E_MUT 0.029926
19 orthologs 0.028661
2 Phosphomimic 0.024752
8 acc 0.023613
3 Acetylmimic 0.021508
15 R_MUT 0.017043
5 nresidues 0.014425
11 ub_1_pfam 0.012936
10 m1_0_pfam 0.012294
6 phi_psi 0.011993
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.665 resistance
FGFR2/L617M/train 0.6 resistance
FGFR2/V564F/train 0.678 resistance
FGFR2/K659M/train 0.434 resistance
FGFR2/L617V/train 0.495 resistance
FGFR2/M537I/train 0.662 resistance
ALK/L1196M/train 0.654 resistance
ALK/G1202R/train 0.38 resistance
ALK/I1171N/train 0.505 resistance
ALK/C1156Y/train 0.592 resistance
ALK/S1206F/train 0.576 resistance
MET/Y1230H/train 0.567 resistance
MET/D1228N/train 0.573 resistance
MET/Y1230C/train 0.412 resistance
MET/D1228V/train 0.557 resistance
MET/Y1230S/train 0.545 resistance
MET/D1228H/train 0.531 resistance
FLT3/F691L/train 0.663 resistance
FLT3/D835I/train 0.469 resistance
EGFR/L792H/train 0.517 resistance
EGFR/G796R/train 0.492 resistance
EGFR/L792R/train 0.457 resistance
EGFR/C797G/train 0.461 resistance
EGFR/L792F/train 0.629 resistance
EGFR/L718Q/train 0.339 resistance
EGFR/C797S/train 0.67 resistance
EGFR/G796D/train 0.519 resistance
EGFR/L718V/train 0.506 resistance
EGFR/G719A/train 0.494 resistance
EGFR/L747P/train 0.378 resistance
EGFR/L792V/train 0.591 resistance
EGFR/G796S/train 0.496 resistance
EGFR/L792P/train 0.455 resistance
EGFR/D761Y/train 0.531 resistance
FGFR4/V550L/train 0.683 resistance
FGFR4/V550M/train 0.697 resistance
ABL1/T315I/train 0.51 resistance
ABL1/G250E/train 0.681 resistance
ABL1/Y253H/train 0.51 resistance
ABL1/E255K/train 0.641 resistance
ABL1/F317L/train 0.601 resistance
ABL1/D276G/train 0.581 resistance
ABL1/F311I/train 0.611 resistance
BRAF/L505H/train 0.551 resistance
KIT/D816E/train 0.727 resistance
KIT/A829P/train 0.634 resistance
KIT/T670I/train 0.622 resistance
KIT/V654A/train 0.456 resistance
KIT/Y823D/train 0.479 resistance
KIT/D820Y/train 0.523 resistance
KIT/N822Y/train 0.688 resistance
KIT/D820G/train 0.679 resistance
KIT/D820E/train 0.723 resistance
KIT/K642E/train 0.704 resistance
KIT/C809G/train 0.613 resistance
KIT/N655T/train 0.613 resistance
KIT/D816G/train 0.594 resistance
ERBB2/L869R/train 0.526 resistance
BTK/C481R/train 0.47 resistance
BTK/C481S/train 0.671 resistance
BTK/C481Y/train 0.519 resistance
MAP2K1/P124L/train 0.457 resistance
MAP2K1/G128V/train 0.569 resistance
MAP2K1/C121S/train 0.52 resistance
MAP2K1/G128D/train 0.633 resistance
MAP2K1/P124S/train 0.524 resistance
MAP2K1/E203K/train 0.634 resistance
MAP2K2/N126D/train 0.698 resistance
MAP2K2/C125S/train 0.52 resistance
MAP2K2/E207K/train 0.634 resistance
KDR/L840F/train 0.421 resistance
########################
MAP2K3/A84T/test 0.662 A
PIM1/S97N/test 0.665 A
########################
CDK5/F82I/test 0.439 TBD
MAP2K1/V211D/test 0.434 TBD
########################
CHEK2/K373E/test 0.61 Inconclusive
MAP2K3/R94L/test 0.59 Inconclusive
MAP2K3/R96W/test 0.389 Inconclusive
MAP2K3/L215W/test 0.393 Inconclusive
MAP2K3/T222M/test 0.41 Inconclusive
PIM1/Q127E/test 0.662 Inconclusive
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.8232734]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 3, 'min_samples_split': 3, 'n_estimators': 100}
[[0.64587872 0.35412128]
[0.65362716 0.34637284]
[0.62198891 0.37801109]
...
[0.78263989 0.21736011]
[0.3750708 0.6249292 ]
[0.78815628 0.21184372]]
########################
hmmScoreWT 0.12866257361914754
hmmScoreMUT 0.08351307729903659
Phosphomimic 0.0410014645861479
Acetylmimic 0.028305503351252037
IUPRED 0.021149195445525917
ncontacts 0.020741201035696842
nresidues 0.019727589140536118
phi_psi 0.023815150772727535
sec 0.012649716197382201
burr 0.03018140192865479
acc 0.012982288773619583
ChargesMUT 0.03272284541572592
K_WT 0.03920544327372399
A_MUT 0.0356616364661129
E_MUT 0.022314812338358588
R_MUT 0.014762703486828819
allHomologs 0.05824141813507941
exclParalogs 0.04233296702496067
specParalogs 0.046126075158182636
orthologs 0.03602267829973058
bpso 0.056362822693421855
bpsh 0.06840855565616207
Feature Importance
0 hmmScoreWT 0.128663
1 hmmScoreMUT 0.083513
21 bpsh 0.068409
16 allHomologs 0.058241
20 bpso 0.056363
18 specParalogs 0.046126
17 exclParalogs 0.042333
2 Phosphomimic 0.041001
12 K_WT 0.039205
19 orthologs 0.036023
13 A_MUT 0.035662
11 ChargesMUT 0.032723
9 burr 0.030181
3 Acetylmimic 0.028306
7 phi_psi 0.023815
14 E_MUT 0.022315
4 IUPRED 0.021149
5 ncontacts 0.020741
6 nresidues 0.019728
15 R_MUT 0.014763
10 acc 0.012982
8 sec 0.012650
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.656 resistance
FGFR2/L617M/train 0.598 resistance
FGFR2/V564F/train 0.673 resistance
FGFR2/K659M/train 0.42 resistance
FGFR2/L617V/train 0.499 resistance
FGFR2/M537I/train 0.674 resistance
ALK/L1196M/train 0.657 resistance
ALK/G1202R/train 0.409 resistance
ALK/I1171N/train 0.49 resistance
ALK/C1156Y/train 0.579 resistance
ALK/S1206F/train 0.57 resistance
MET/Y1230H/train 0.566 resistance
MET/D1228N/train 0.567 resistance
MET/Y1230C/train 0.435 resistance
MET/D1228V/train 0.541 resistance
MET/Y1230S/train 0.542 resistance
MET/D1228H/train 0.524 resistance
FLT3/F691L/train 0.664 resistance
FLT3/D835I/train 0.494 resistance
EGFR/L792H/train 0.531 resistance
EGFR/G796R/train 0.503 resistance
EGFR/L792R/train 0.48 resistance
EGFR/C797G/train 0.464 resistance
EGFR/L792F/train 0.625 resistance
EGFR/L718Q/train 0.34 resistance
EGFR/C797S/train 0.668 resistance
EGFR/G796D/train 0.519 resistance
EGFR/L718V/train 0.538 resistance
EGFR/G719A/train 0.496 resistance
EGFR/L747P/train 0.375 resistance
EGFR/L792V/train 0.603 resistance
EGFR/G796S/train 0.525 resistance
EGFR/L792P/train 0.484 resistance
EGFR/D761Y/train 0.535 resistance
FGFR4/V550L/train 0.67 resistance
FGFR4/V550M/train 0.673 resistance
ABL1/T315I/train 0.539 resistance
ABL1/G250E/train 0.662 resistance
ABL1/Y253H/train 0.508 resistance
ABL1/E255K/train 0.619 resistance
ABL1/F317L/train 0.6 resistance
ABL1/D276G/train 0.587 resistance
ABL1/F311I/train 0.619 resistance
BRAF/L505H/train 0.54 resistance
KIT/D816E/train 0.713 resistance
KIT/A829P/train 0.609 resistance
KIT/T670I/train 0.634 resistance
KIT/V654A/train 0.461 resistance
KIT/Y823D/train 0.467 resistance
KIT/D820Y/train 0.535 resistance
KIT/N822Y/train 0.678 resistance
KIT/D820G/train 0.646 resistance
KIT/D820E/train 0.704 resistance
KIT/K642E/train 0.653 resistance
KIT/C809G/train 0.604 resistance
KIT/N655T/train 0.584 resistance
KIT/D816G/train 0.565 resistance
ERBB2/L869R/train 0.546 resistance
BTK/C481R/train 0.495 resistance
BTK/C481S/train 0.666 resistance
BTK/C481Y/train 0.507 resistance
MAP2K1/P124L/train 0.474 resistance
MAP2K1/G128V/train 0.589 resistance
MAP2K1/C121S/train 0.499 resistance
MAP2K1/G128D/train 0.626 resistance
MAP2K1/P124S/train 0.541 resistance
MAP2K1/E203K/train 0.624 resistance
MAP2K2/N126D/train 0.692 resistance
MAP2K2/C125S/train 0.499 resistance
MAP2K2/E207K/train 0.624 resistance
KDR/L840F/train 0.439 resistance
########################
MAP2K3/A84T/test 0.65 A
PIM1/S97N/test 0.655 A
########################
CDK5/F82I/test 0.47 TBD
MAP2K1/V211D/test 0.459 TBD
########################
CHEK2/K373E/test 0.595 Inconclusive
MAP2K3/R94L/test 0.592 Inconclusive
MAP2K3/R96W/test 0.393 Inconclusive
MAP2K3/L215W/test 0.417 Inconclusive
MAP2K3/T222M/test 0.441 Inconclusive
PIM1/Q127E/test 0.65 Inconclusive
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.80561944]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 4, 'min_samples_split': 3, 'n_estimators': 10}
[[0.63377024 0.36622976]
[0.66049186 0.33950814]
[0.6199757 0.3800243 ]
...
[0.80121612 0.19878388]
[0.43654535 0.56345465]
[0.76149061 0.23850939]]
########################
hmmScoreWT 0.1361752544826065
hmmScoreMUT 0.08044656598735037
ncontacts 0.08464678292907128
nresidues 0.02062274265836489
burr 0.07429643696725505
acc 0.013332866576955921
ChargesMUT 0.06967812962556595
p_0_pfam 0.011947360571989299
ub_0_pfam 0.011321945396237602
p_1_pfam 0.014745912138524383
ub_1_pfam 0.020787527487184848
K_WT 0.059753109438451316
S_WT 0.018699196801251407
V_WT 0.012094037592429632
A_MUT 0.027516662706600043
H_MUT 0.011818329164687398
R_MUT 0.04159408890489486
allHomologs 0.011602643538301646
exclParalogs 0.019236930788302267
specParalogs 0.05372919125774371
orthologs 0.03577788887820925
bpso 0.028487869707095215
bpsh 0.09761789994393241
Feature Importance
0 hmmScoreWT 0.136175
22 bpsh 0.097618
2 ncontacts 0.084647
1 hmmScoreMUT 0.080447
4 burr 0.074296
6 ChargesMUT 0.069678
11 K_WT 0.059753
19 specParalogs 0.053729
16 R_MUT 0.041594
20 orthologs 0.035778
21 bpso 0.028488
14 A_MUT 0.027517
10 ub_1_pfam 0.020788
3 nresidues 0.020623
18 exclParalogs 0.019237
12 S_WT 0.018699
9 p_1_pfam 0.014746
5 acc 0.013333
13 V_WT 0.012094
7 p_0_pfam 0.011947
15 H_MUT 0.011818
17 allHomologs 0.011603
8 ub_0_pfam 0.011322
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.665 resistance
FGFR2/L617M/train 0.562 resistance
FGFR2/V564F/train 0.7 resistance
FGFR2/K659M/train 0.458 resistance
FGFR2/L617V/train 0.373 resistance
FGFR2/M537I/train 0.68 resistance
ALK/L1196M/train 0.621 resistance
ALK/G1202R/train 0.309 resistance
ALK/I1171N/train 0.433 resistance
ALK/C1156Y/train 0.606 resistance
ALK/S1206F/train 0.605 resistance
MET/Y1230H/train 0.588 resistance
MET/D1228N/train 0.551 resistance
MET/Y1230C/train 0.397 resistance
MET/D1228V/train 0.575 resistance
MET/Y1230S/train 0.464 resistance
MET/D1228H/train 0.602 resistance
FLT3/F691L/train 0.621 resistance
FLT3/D835I/train 0.475 resistance
EGFR/L792H/train 0.461 resistance
EGFR/G796R/train 0.429 resistance
EGFR/L792R/train 0.386 resistance
EGFR/C797G/train 0.451 resistance
EGFR/L792F/train 0.591 resistance
EGFR/L718Q/train 0.387 resistance
EGFR/C797S/train 0.665 resistance
EGFR/G796D/train 0.437 resistance
EGFR/L718V/train 0.452 resistance
EGFR/G719A/train 0.419 resistance
EGFR/L747P/train 0.297 resistance
EGFR/L792V/train 0.546 resistance
EGFR/G796S/train 0.471 resistance
EGFR/L792P/train 0.4 resistance
EGFR/D761Y/train 0.564 resistance
FGFR4/V550L/train 0.7 resistance
FGFR4/V550M/train 0.7 resistance
ABL1/T315I/train 0.45 resistance
ABL1/G250E/train 0.712 resistance
ABL1/Y253H/train 0.517 resistance
ABL1/E255K/train 0.563 resistance
ABL1/F317L/train 0.564 resistance
ABL1/D276G/train 0.577 resistance
ABL1/F311I/train 0.564 resistance
BRAF/L505H/train 0.527 resistance
KIT/D816E/train 0.749 resistance
KIT/A829P/train 0.642 resistance
KIT/T670I/train 0.62 resistance
KIT/V654A/train 0.394 resistance
KIT/Y823D/train 0.437 resistance
KIT/D820Y/train 0.456 resistance
KIT/N822Y/train 0.72 resistance
KIT/D820G/train 0.641 resistance
KIT/D820E/train 0.711 resistance
KIT/K642E/train 0.756 resistance
KIT/C809G/train 0.514 resistance
KIT/N655T/train 0.606 resistance
KIT/D816G/train 0.551 resistance
ERBB2/L869R/train 0.489 resistance
BTK/C481R/train 0.515 resistance
BTK/C481S/train 0.68 resistance
BTK/C481Y/train 0.56 resistance
MAP2K1/P124L/train 0.422 resistance
MAP2K1/G128V/train 0.483 resistance
MAP2K1/C121S/train 0.542 resistance
MAP2K1/G128D/train 0.585 resistance
MAP2K1/P124S/train 0.448 resistance
MAP2K1/E203K/train 0.563 resistance
MAP2K2/N126D/train 0.735 resistance
MAP2K2/C125S/train 0.542 resistance
MAP2K2/E207K/train 0.563 resistance
KDR/L840F/train 0.457 resistance
########################
MAP2K3/A84T/test 0.665 A
PIM1/S97N/test 0.687 A
########################
CDK5/F82I/test 0.438 TBD
MAP2K1/V211D/test 0.322 TBD
########################
CHEK2/K373E/test 0.506 Inconclusive
MAP2K3/R94L/test 0.602 Inconclusive
MAP2K3/R96W/test 0.451 Inconclusive
MAP2K3/L215W/test 0.36 Inconclusive
MAP2K3/T222M/test 0.408 Inconclusive
PIM1/Q127E/test 0.6 Inconclusive
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.8201002]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 4, 'min_samples_split': 3, 'n_estimators': 50}
[[0.65435532 0.34564468]
[0.649739 0.350261 ]
[0.62774155 0.37225845]
...
[0.77679927 0.22320073]
[0.34805296 0.65194704]
[0.77790917 0.22209083]]
########################
hmmScoreWT 0.12631977936478522
hmmScoreMUT 0.08744226367137527
Phosphomimic 0.017159607781018563
Acetylmimic 0.021852356713609893
ncontacts 0.03302285772477565
nresidues 0.01736074537600189
phi_psi 0.012064162568260142
burr 0.04330268211184773
acc 0.02363031632705435
ChargesMUT 0.03523533787240301
m1_0_pfam 0.012773579167634797
ub_1_pfam 0.013463370385646878
K_WT 0.0363941342080859
A_MUT 0.04358335155491667
E_MUT 0.022102065538456766
R_MUT 0.016223281578957954
allHomologs 0.06710830926374443
exclParalogs 0.04125701147383701
specParalogs 0.043176932199550024
orthologs 0.02916353973449406
bpso 0.053017437141738434
bpsh 0.09088938842652393
Feature Importance
0 hmmScoreWT 0.126320
21 bpsh 0.090889
1 hmmScoreMUT 0.087442
16 allHomologs 0.067108
20 bpso 0.053017
13 A_MUT 0.043583
7 burr 0.043303
18 specParalogs 0.043177
17 exclParalogs 0.041257
12 K_WT 0.036394
9 ChargesMUT 0.035235
4 ncontacts 0.033023
19 orthologs 0.029164
8 acc 0.023630
14 E_MUT 0.022102
3 Acetylmimic 0.021852
5 nresidues 0.017361
2 Phosphomimic 0.017160
15 R_MUT 0.016223
11 ub_1_pfam 0.013463
10 m1_0_pfam 0.012774
6 phi_psi 0.012064
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.667 resistance
FGFR2/L617M/train 0.596 resistance
FGFR2/V564F/train 0.68 resistance
FGFR2/K659M/train 0.436 resistance
FGFR2/L617V/train 0.494 resistance
FGFR2/M537I/train 0.664 resistance
ALK/L1196M/train 0.655 resistance
ALK/G1202R/train 0.38 resistance
ALK/I1171N/train 0.506 resistance
ALK/C1156Y/train 0.592 resistance
ALK/S1206F/train 0.577 resistance
MET/Y1230H/train 0.574 resistance
MET/D1228N/train 0.574 resistance
MET/Y1230C/train 0.413 resistance
MET/D1228V/train 0.557 resistance
MET/Y1230S/train 0.549 resistance
MET/D1228H/train 0.533 resistance
FLT3/F691L/train 0.663 resistance
FLT3/D835I/train 0.469 resistance
EGFR/L792H/train 0.516 resistance
EGFR/G796R/train 0.492 resistance
EGFR/L792R/train 0.455 resistance
EGFR/C797G/train 0.462 resistance
EGFR/L792F/train 0.63 resistance
EGFR/L718Q/train 0.344 resistance
EGFR/C797S/train 0.672 resistance
EGFR/G796D/train 0.52 resistance
EGFR/L718V/train 0.501 resistance
EGFR/G719A/train 0.49 resistance
EGFR/L747P/train 0.378 resistance
EGFR/L792V/train 0.588 resistance
EGFR/G796S/train 0.496 resistance
EGFR/L792P/train 0.455 resistance
EGFR/D761Y/train 0.532 resistance
FGFR4/V550L/train 0.679 resistance
FGFR4/V550M/train 0.699 resistance
ABL1/T315I/train 0.509 resistance
ABL1/G250E/train 0.674 resistance
ABL1/Y253H/train 0.514 resistance
ABL1/E255K/train 0.64 resistance
ABL1/F317L/train 0.595 resistance
ABL1/D276G/train 0.581 resistance
ABL1/F311I/train 0.608 resistance
BRAF/L505H/train 0.553 resistance
KIT/D816E/train 0.721 resistance
KIT/A829P/train 0.634 resistance
KIT/T670I/train 0.62 resistance
KIT/V654A/train 0.457 resistance
KIT/Y823D/train 0.477 resistance
KIT/D820Y/train 0.525 resistance
KIT/N822Y/train 0.69 resistance
KIT/D820G/train 0.671 resistance
KIT/D820E/train 0.71 resistance
KIT/K642E/train 0.693 resistance
KIT/C809G/train 0.612 resistance
KIT/N655T/train 0.613 resistance
KIT/D816G/train 0.594 resistance
ERBB2/L869R/train 0.526 resistance
BTK/C481R/train 0.471 resistance
BTK/C481S/train 0.673 resistance
BTK/C481Y/train 0.52 resistance
MAP2K1/P124L/train 0.464 resistance
MAP2K1/G128V/train 0.566 resistance
MAP2K1/C121S/train 0.52 resistance
MAP2K1/G128D/train 0.628 resistance
MAP2K1/P124S/train 0.524 resistance
MAP2K1/E203K/train 0.633 resistance
MAP2K2/N126D/train 0.7 resistance
MAP2K2/C125S/train 0.52 resistance
MAP2K2/E207K/train 0.633 resistance
KDR/L840F/train 0.42 resistance
########################
MAP2K3/A84T/test 0.661 A
PIM1/S97N/test 0.665 A
########################
CDK5/F82I/test 0.438 TBD
MAP2K1/V211D/test 0.431 TBD
########################
CHEK2/K373E/test 0.6 Inconclusive
MAP2K3/R94L/test 0.587 Inconclusive
MAP2K3/R96W/test 0.39 Inconclusive
MAP2K3/L215W/test 0.393 Inconclusive
MAP2K3/T222M/test 0.43 Inconclusive
PIM1/Q127E/test 0.648 Inconclusive
columns to consider ['Acc' 'Gene' 'Mutation' 'hmmScoreWT' 'hmmScoreMUT' 'Phosphomimic'
'Acetylmimic' 'IUPRED' 'ncontacts' 'nresidues' 'phi_psi' 'sec' 'burr'
'acc' 'ChargesWT' 'ChargesMUT' 'ac_-1' 'gl_-1' 'm1_-1' 'm2_-1' 'm3_-1'
'me_-1' 'p_-1' 'sm_-1' 'ub_-1' 'ac_-1_pfam' 'gl_-1_pfam' 'm1_-1_pfam'
'm2_-1_pfam' 'm3_-1_pfam' 'me_-1_pfam' 'p_-1_pfam' 'sm_-1_pfam'
'ub_-1_pfam' 'ac_0' 'gl_0' 'm1_0' 'm2_0' 'm3_0' 'me_0' 'p_0' 'sm_0'
'ub_0' 'ac_0_pfam' 'gl_0_pfam' 'm1_0_pfam' 'm2_0_pfam' 'm3_0_pfam'
'me_0_pfam' 'p_0_pfam' 'sm_0_pfam' 'ub_0_pfam' 'ac_1' 'gl_1' 'm1_1'
'm2_1' 'm3_1' 'me_1' 'p_1' 'sm_1' 'ub_1' 'ac_1_pfam' 'gl_1_pfam'
'm1_1_pfam' 'm2_1_pfam' 'm3_1_pfam' 'me_1_pfam' 'p_1_pfam' 'sm_1_pfam'
'ub_1_pfam' 'A_WT' 'C_WT' 'D_WT' 'E_WT' 'F_WT' 'G_WT' 'H_WT' 'I_WT'
'K_WT' 'L_WT' 'M_WT' 'N_WT' 'P_WT' 'Q_WT' 'R_WT' 'S_WT' 'T_WT' 'V_WT'
'W_WT' 'Y_WT' 'A_MUT' 'C_MUT' 'D_MUT' 'E_MUT' 'F_MUT' 'G_MUT' 'H_MUT'
'I_MUT' 'K_MUT' 'L_MUT' 'M_MUT' 'N_MUT' 'P_MUT' 'Q_MUT' 'R_MUT' 'S_MUT'
'T_MUT' 'V_MUT' 'W_MUT' 'Y_MUT' 'allHomologs' 'exclParalogs'
'specParalogs' 'orthologs' 'bpso' 'bpsh' 'A_-1' 'D_-1' 'R_-1' 'A_-1_pfam'
'D_-1_pfam' 'R_-1_pfam' 'A_0_pfam' 'D_0_pfam' 'R_0_pfam' 'A_1' 'D_1'
'R_1' 'A_1_pfam' 'D_1_pfam' 'R_1_pfam' 'MUT_TYPE']
mean_test_score [0.82186219]
-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-#-
Best model found during the CV
{'max_depth': 3, 'max_features': 'log2', 'min_samples_leaf': 4, 'min_samples_split': 3, 'n_estimators': 100}
[[0.63984863 0.36015137]
[0.65625457 0.34374543]
[0.62154099 0.37845901]
...
[0.78804362 0.21195638]
[0.37332265 0.62667735]
[0.78919547 0.21080453]]
########################
hmmScoreWT 0.12945362364724156
hmmScoreMUT 0.0837826265248172
Phosphomimic 0.041453586826247044
Acetylmimic 0.02779672423055593
IUPRED 0.02246622141489417
ncontacts 0.023019487807961467
nresidues 0.02099401733450277
phi_psi 0.02378103147284423
sec 0.012402400828476844
burr 0.03346442784335046
acc 0.012990945118416121
ChargesMUT 0.03322850733668184
K_WT 0.039238467415559114
A_MUT 0.035811116573955644
E_MUT 0.014734751820939258
R_MUT 0.015569173449450923
allHomologs 0.0618339604845394
exclParalogs 0.04297373200539816
specParalogs 0.04630727736784269
orthologs 0.032695996527537255
bpso 0.054640080217575646
bpsh 0.06947829969729515
Feature Importance
0 hmmScoreWT 0.129454
1 hmmScoreMUT 0.083783
21 bpsh 0.069478
16 allHomologs 0.061834
20 bpso 0.054640
18 specParalogs 0.046307
17 exclParalogs 0.042974
2 Phosphomimic 0.041454
12 K_WT 0.039238
13 A_MUT 0.035811
9 burr 0.033464
11 ChargesMUT 0.033229
19 orthologs 0.032696
3 Acetylmimic 0.027797
7 phi_psi 0.023781
5 ncontacts 0.023019
4 IUPRED 0.022466
6 nresidues 0.020994
15 R_MUT 0.015569
14 E_MUT 0.014735
10 acc 0.012991
8 sec 0.012402
########################
activatingresistance results ( 14 )
REC: 0.857
########################
increaseresistance results ( 1 )
REC: 1.0
########################
FGFR2/N549K/train 0.66 resistance
FGFR2/L617M/train 0.598 resistance
FGFR2/V564F/train 0.678 resistance
FGFR2/K659M/train 0.416 resistance
FGFR2/L617V/train 0.499 resistance
FGFR2/M537I/train 0.675 resistance
ALK/L1196M/train 0.66 resistance
ALK/G1202R/train 0.406 resistance
ALK/I1171N/train 0.496 resistance
ALK/C1156Y/train 0.585 resistance
ALK/S1206F/train 0.577 resistance
MET/Y1230H/train 0.568 resistance
MET/D1228N/train 0.564 resistance
MET/Y1230C/train 0.435 resistance
MET/D1228V/train 0.542 resistance
MET/Y1230S/train 0.542 resistance
MET/D1228H/train 0.526 resistance
FLT3/F691L/train 0.666 resistance
FLT3/D835I/train 0.491 resistance
EGFR/L792H/train 0.529 resistance
EGFR/G796R/train 0.503 resistance
EGFR/L792R/train 0.479 resistance
EGFR/C797G/train 0.466 resistance
EGFR/L792F/train 0.628 resistance
EGFR/L718Q/train 0.34 resistance
EGFR/C797S/train 0.673 resistance
EGFR/G796D/train 0.52 resistance
EGFR/L718V/train 0.533 resistance
EGFR/G719A/train 0.491 resistance
EGFR/L747P/train 0.373 resistance
EGFR/L792V/train 0.603 resistance
EGFR/G796S/train 0.526 resistance
EGFR/L792P/train 0.481 resistance
EGFR/D761Y/train 0.536 resistance
FGFR4/V550L/train 0.669 resistance
FGFR4/V550M/train 0.675 resistance
ABL1/T315I/train 0.539 resistance
ABL1/G250E/train 0.658 resistance
ABL1/Y253H/train 0.516 resistance
ABL1/E255K/train 0.62 resistance
ABL1/F317L/train 0.599 resistance
ABL1/D276G/train 0.591 resistance
ABL1/F311I/train 0.616 resistance
BRAF/L505H/train 0.543 resistance
KIT/D816E/train 0.709 resistance
KIT/A829P/train 0.612 resistance
KIT/T670I/train 0.631 resistance
KIT/V654A/train 0.462 resistance
KIT/Y823D/train 0.472 resistance
KIT/D820Y/train 0.54 resistance
KIT/N822Y/train 0.685 resistance
KIT/D820G/train 0.645 resistance
KIT/D820E/train 0.695 resistance
KIT/K642E/train 0.648 resistance
KIT/C809G/train 0.606 resistance
KIT/N655T/train 0.585 resistance
KIT/D816G/train 0.563 resistance
ERBB2/L869R/train 0.549 resistance
BTK/C481R/train 0.497 resistance
BTK/C481S/train 0.671 resistance
BTK/C481Y/train 0.517 resistance
MAP2K1/P124L/train 0.479 resistance
MAP2K1/G128V/train 0.592 resistance
MAP2K1/C121S/train 0.499 resistance
MAP2K1/G128D/train 0.627 resistance
MAP2K1/P124S/train 0.541 resistance
MAP2K1/E203K/train 0.625 resistance
MAP2K2/N126D/train 0.696 resistance
MAP2K2/C125S/train 0.499 resistance
MAP2K2/E207K/train 0.625 resistance
KDR/L840F/train 0.436 resistance
########################
MAP2K3/A84T/test 0.65 A
PIM1/S97N/test 0.659 A
########################
CDK5/F82I/test 0.466 TBD
MAP2K1/V211D/test 0.459 TBD
########################
CHEK2/K373E/test 0.581 Inconclusive
MAP2K3/R94L/test 0.595 Inconclusive
MAP2K3/R96W/test 0.396 Inconclusive
MAP2K3/L215W/test 0.42 Inconclusive