@@ -674,6 +674,8 @@ In this execution mode the HADDOCK3 job should be submitted to the batch system
674674### Scenario 1: 1000 rigidbody docking models, selection of top 200 and flexible refinement + EM
675675
676676Now that we have all data ready, and know about execution modes of HADDOCK3 it is time to setup the docking for the first scenario. The restraint file to use for this is ` e2a-hpr_air.tbl ` . We proceed to produce 1000 rigidbody docking models, from which 200 will be selected and refined through flexible refinement and energy minimization.
677+ This corresponds to the default docking scenario of the HADDOCK2.X version.
678+
677679For the analysis following the docking results, we are using the solved complex [ 1GGR] ( https://www.rcsb.org/structure/1GGR ) , named e2a-hpr_1GGR.pdb.
678680The configuration file for this scenario is:
679681
@@ -756,6 +758,8 @@ submitting it to the batch system requesting in this local run mode a full node
756758### Scenario 2: 1000 rigidbody docking models, FCC clustering and selection of max 20 models per cluster followed by flexible refinement and EM
757759
758760In scenario 2, we proceed to produce 1000 rigidbody docking models, from which we proceed to do a first clustering analysis. From the top clusters a flexible refinement then energy minization is done.
761+ This scenario illustrates the new flexibility of HADDOCK3, adding a clustering step after rigid-body docking, which is not possible in the HADDOCK2.X version.
762+
759763For the analysis following the docking results, we are using the solved complex [ 1GGR] ( https://www.rcsb.org/structure/1GGR ) , named e2a-hpr_1GGR.pdb.
760764The configuration file for this scenario is:
761765
@@ -968,7 +972,7 @@ We are providing in the `scripts` directory a simple script that extract some cl
968972To use is simply call the script with as argument the run directory you want to analyze, e.g.:
969973
970974<a class =" prompt prompt-cmd " >
971- ./scripts/extract-capri-stats-clt.sh ./scenario1
975+ ./scripts/extract-capri-stats-clt.sh ./run1-full
972976</a >
973977
974978<details style =" background-color :#DAE4E7 " >
@@ -977,7 +981,7 @@ To use is simply call the script with as argument the run directory you want to
977981 </summary >
978982<pre >
979983==============================================
980- == scenario1 -full/02_caprieval/capri_clt.tsv
984+ == run1 -full/02_caprieval/capri_clt.tsv
981985==============================================
982986Total number of acceptable or better clusters: 0 out of 1
983987Total number of medium or better clusters: 0 out of 1
@@ -987,7 +991,7 @@ First acceptable cluster - rank: i-RMSD: Fnat: DockQ:
987991First medium cluster - rank: i-RMSD: Fnat: DockQ:
988992Best cluster - rank: - i-RMSD: 6.407 Fnat: 0.202 DockQ: 0.264
989993==============================================
990- == scenario1 -full/04_caprieval/capri_clt.tsv
994+ == run1 -full/04_caprieval/capri_clt.tsv
991995==============================================
992996Total number of acceptable or better clusters: 0 out of 1
993997Total number of medium or better clusters: 0 out of 1
@@ -997,7 +1001,7 @@ First acceptable cluster - rank: i-RMSD: Fnat: DockQ:
9971001First medium cluster - rank: i-RMSD: Fnat: DockQ:
9981002Best cluster - rank: - i-RMSD: 6.407 Fnat: 0.202 DockQ: 0.264
9991003==============================================
1000- == scenario1 -full/06_caprieval/capri_clt.tsv
1004+ == run1 -full/06_caprieval/capri_clt.tsv
10011005==============================================
10021006Total number of acceptable or better clusters: 1 out of 1
10031007Total number of medium or better clusters: 0 out of 1
@@ -1007,7 +1011,7 @@ First acceptable cluster - rank: - i-RMSD: 2.976 Fnat: 0.611 DockQ: 0.601
10071011First medium cluster - rank: i-RMSD: Fnat: DockQ:
10081012Best cluster - rank: - i-RMSD: 2.976 Fnat: 0.611 DockQ: 0.601
10091013==============================================
1010- == scenario1 -full/08_caprieval/capri_clt.tsv
1014+ == run1 -full/08_caprieval/capri_clt.tsv
10111015==============================================
10121016Total number of acceptable or better clusters: 1 out of 1
10131017Total number of medium or better clusters: 1 out of 1
@@ -1017,7 +1021,7 @@ First acceptable cluster - rank: - i-RMSD: 1.673 Fnat: 0.736 DockQ: 0.727
10171021First medium cluster - rank: - i-RMSD: 1.673 Fnat: 0.736 DockQ: 0.727
10181022Best cluster - rank: - i-RMSD: 1.673 Fnat: 0.736 DockQ: 0.727
10191023==============================================
1020- == scenario1 -full/11_caprieval/capri_clt.tsv
1024+ == run1 -full/11_caprieval/capri_clt.tsv
10211025==============================================
10221026Total number of acceptable or better clusters: 3 out of 4
10231027Total number of medium or better clusters: 1 out of 4
@@ -1035,7 +1039,7 @@ Similarly some simple statistics can be extracted from the single model `capriev
10351039
10361040
10371041<a class =" prompt prompt-cmd " >
1038- ./scripts/extract-capri-stats.sh ./runs/scenario1-surface
1042+ ./scripts/extract-capri-stats.sh ./runs/run1-full
10391043</a >
10401044
10411045<details style =" background-color :#DAE4E7 " >
@@ -1044,7 +1048,7 @@ Similarly some simple statistics can be extracted from the single model `capriev
10441048 </summary >
10451049<pre >
10461050==============================================
1047- == scenario1 -full/02_caprieval/capri_ss.tsv
1051+ == run1 -full/02_caprieval/capri_ss.tsv
10481052==============================================
10491053Total number of acceptable or better models: 365 out of 1000
10501054Total number of medium or better models: 199 out of 1000
@@ -1054,7 +1058,7 @@ First acceptable model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
10541058First medium model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
10551059Best model - rank: 46 i-RMSD: 1.145 Fnat: 0.556 DockQ: 0.713
10561060==============================================
1057- == scenario1 -full/04_caprieval/capri_ss.tsv
1061+ == run1 -full/04_caprieval/capri_ss.tsv
10581062==============================================
10591063Total number of acceptable or better models: 144 out of 200
10601064Total number of medium or better models: 137 out of 200
@@ -1064,7 +1068,7 @@ First acceptable model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
10641068First medium model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
10651069Best model - rank: 46 i-RMSD: 1.145 Fnat: 0.556 DockQ: 0.713
10661070==============================================
1067- == scenario1 -full/06_caprieval/capri_ss.tsv
1071+ == run1 -full/06_caprieval/capri_ss.tsv
10681072==============================================
10691073Total number of acceptable or better models: 147 out of 200
10701074Total number of medium or better models: 118 out of 200
@@ -1074,7 +1078,7 @@ First acceptable model - rank: 2 i-RMSD: 1.221 Fnat: 0.694 DockQ: 0.727
10741078First medium model - rank: 2 i-RMSD: 1.221 Fnat: 0.694 DockQ: 0.727
10751079Best model - rank: 30 i-RMSD: 0.883 Fnat: 0.750 DockQ: 0.823
10761080==============================================
1077- == scenario1 -full/08_caprieval/capri_ss.tsv
1081+ == run1 -full/08_caprieval/capri_ss.tsv
10781082==============================================
10791083Total number of acceptable or better models: 147 out of 200
10801084Total number of medium or better models: 118 out of 200
@@ -1084,7 +1088,7 @@ First acceptable model - rank: 1 i-RMSD: 1.219 Fnat: 0.833 DockQ: 0.787
10841088First medium model - rank: 1 i-RMSD: 1.219 Fnat: 0.833 DockQ: 0.787
10851089Best model - rank: 39 i-RMSD: 0.807 Fnat: 0.833 DockQ: 0.862
10861090==============================================
1087- == scenario1 -full/11_caprieval/capri_ss.tsv
1091+ == run1 -full/11_caprieval/capri_ss.tsv
10881092==============================================
10891093Total number of acceptable or better models: 141 out of 185
10901094Total number of medium or better models: 116 out of 185
@@ -1232,7 +1236,7 @@ cluster_rank cluster_id n under_eval score score_std
12321236Use the ` extract-capri-stats-clt.sh ` script to extract some simple cluster statistics for this run.
12331237
12341238<a class =" prompt prompt-cmd " >
1235- ./scripts/extract-capri-stats-clt.sh runs/scenario2 /
1239+ ./scripts/extract-capri-stats-clt.sh runs/run2-full /
12361240</a >
12371241
12381242
@@ -1242,7 +1246,7 @@ Use the `extract-capri-stats-clt.sh` script to extract some simple cluster stati
12421246 </summary >
12431247<pre >
12441248==============================================
1245- == scenario2-cltsel -full/02_caprieval/capri_clt.tsv
1249+ == run2 -full/02_caprieval/capri_clt.tsv
12461250==============================================
12471251Total number of acceptable or better clusters: 0 out of 1
12481252Total number of medium or better clusters: 0 out of 1
@@ -1252,7 +1256,7 @@ First acceptable cluster - rank: i-RMSD: Fnat: DockQ:
12521256First medium cluster - rank: i-RMSD: Fnat: DockQ:
12531257Best cluster - rank: - i-RMSD: 6.407 Fnat: 0.202 DockQ: 0.264
12541258==============================================
1255- == scenario2-cltsel -full/05_caprieval/capri_clt.tsv
1259+ == run2 -full/05_caprieval/capri_clt.tsv
12561260==============================================
12571261Total number of acceptable or better clusters: 6 out of 33
12581262Total number of medium or better clusters: 1 out of 33
@@ -1262,7 +1266,7 @@ First acceptable cluster - rank: 2 i-RMSD: 1.193 Fnat: 0.563 DockQ: 0.701
12621266First medium cluster - rank: 2 i-RMSD: 1.193 Fnat: 0.563 DockQ: 0.701
12631267Best cluster - rank: 2 i-RMSD: 1.193 Fnat: 0.563 DockQ: 0.701
12641268==============================================
1265- == scenario2-cltsel -full/07_caprieval/capri_clt.tsv
1269+ == run2 -full/07_caprieval/capri_clt.tsv
12661270==============================================
12671271Total number of acceptable or better clusters: 0 out of 1
12681272Total number of medium or better clusters: 0 out of 1
@@ -1272,7 +1276,7 @@ First acceptable cluster - rank: i-RMSD: Fnat: DockQ:
12721276First medium cluster - rank: i-RMSD: Fnat: DockQ:
12731277Best cluster - rank: - i-RMSD: 8.237 Fnat: 0.104 DockQ: 0.121
12741278==============================================
1275- == scenario2-cltsel -full/09_caprieval/capri_clt.tsv
1279+ == run2 -full/09_caprieval/capri_clt.tsv
12761280==============================================
12771281Total number of acceptable or better clusters: 0 out of 1
12781282Total number of medium or better clusters: 0 out of 1
@@ -1282,7 +1286,7 @@ First acceptable cluster - rank: i-RMSD: Fnat: DockQ:
12821286First medium cluster - rank: i-RMSD: Fnat: DockQ:
12831287Best cluster - rank: - i-RMSD: 4.840 Fnat: 0.361 DockQ: 0.400
12841288==============================================
1285- == scenario2-cltsel -full/12_caprieval/capri_clt.tsv
1289+ == run2 -full/12_caprieval/capri_clt.tsv
12861290==============================================
12871291Total number of acceptable or better clusters: 4 out of 25
12881292Total number of medium or better clusters: 2 out of 25
@@ -1299,7 +1303,7 @@ Best cluster - rank: 4 i-RMSD: 0.980 Fnat: 0.812 DockQ: 0.819
12991303Similarly some simple statistics can be extracted from the single model ` caprieval ` ` capri_ss.tsv ` files with the ` extract-capri-stats.sh ` script:
13001304
13011305<a class =" prompt prompt-cmd " >
1302- ./scripts/extract-capri-stats.sh ./runs/scenario2
1306+ ./scripts/extract-capri-stats.sh ./runs/run2-full
13031307</a >
13041308
13051309<details style =" background-color :#DAE4E7 " >
@@ -1308,7 +1312,7 @@ Similarly some simple statistics can be extracted from the single model `capriev
13081312 </summary >
13091313<pre >
13101314==============================================
1311- == scenario2-cltsel -full/02_caprieval/capri_ss.tsv
1315+ == run2 -full/02_caprieval/capri_ss.tsv
13121316==============================================
13131317Total number of acceptable or better models: 365 out of 1000
13141318Total number of medium or better models: 199 out of 1000
@@ -1318,7 +1322,7 @@ First acceptable model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
13181322First medium model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
13191323Best model - rank: 46 i-RMSD: 1.145 Fnat: 0.556 DockQ: 0.713
13201324==============================================
1321- == scenario2-cltsel -full/05_caprieval/capri_ss.tsv
1325+ == run2 -full/05_caprieval/capri_ss.tsv
13221326==============================================
13231327Total number of acceptable or better models: 62 out of 375
13241328Total number of medium or better models: 22 out of 375
@@ -1328,7 +1332,7 @@ First acceptable model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
13281332First medium model - rank: 3 i-RMSD: 1.153 Fnat: 0.556 DockQ: 0.711
13291333Best model - rank: 46 i-RMSD: 1.145 Fnat: 0.556 DockQ: 0.713
13301334==============================================
1331- == scenario2-cltsel -full/07_caprieval/capri_ss.tsv
1335+ == run2 -full/07_caprieval/capri_ss.tsv
13321336==============================================
13331337Total number of acceptable or better models: 74 out of 375
13341338Total number of medium or better models: 27 out of 375
@@ -1338,7 +1342,7 @@ First acceptable model - rank: 6 i-RMSD: 1.081 Fnat: 0.750 DockQ: 0.771
13381342First medium model - rank: 6 i-RMSD: 1.081 Fnat: 0.750 DockQ: 0.771
13391343Best model - rank: 36 i-RMSD: 0.930 Fnat: 0.778 DockQ: 0.822
13401344==============================================
1341- == scenario2-cltsel -full/09_caprieval/capri_ss.tsv
1345+ == run2 -full/09_caprieval/capri_ss.tsv
13421346==============================================
13431347Total number of acceptable or better models: 74 out of 375
13441348Total number of medium or better models: 27 out of 375
@@ -1348,7 +1352,7 @@ First acceptable model - rank: 1 i-RMSD: 3.718 Fnat: 0.333 DockQ: 0.382
13481352First medium model - rank: 3 i-RMSD: 0.991 Fnat: 0.806 DockQ: 0.821
13491353Best model - rank: 60 i-RMSD: 0.896 Fnat: 0.778 DockQ: 0.828
13501354==============================================
1351- == scenario2-cltsel -full/12_caprieval/capri_ss.tsv
1355+ == run2 -full/12_caprieval/capri_ss.tsv
13521356==============================================
13531357Total number of acceptable or better models: 65 out of 317
13541358Total number of medium or better models: 27 out of 317
0 commit comments