Basic instructions to screen huge chemical collections with rDock

rDock is a fast & free docking program designed to be used in High-Throughput Virtual Screening. This is a summary of the steps necessary to do so.

1) You mst have previously defined yor system & cavity (consult the manual)

2) Execute 50 docking runs for a set of ligands representative og your entire collection:

rbdock -i $f -o ${f}_out -r REC.prm -p dock.prm -n 50 > ${f}_out}.log

So, your output files are named *_out.sd. Create a report file:

sdreport -t *.sd > results_table.txt

Now execute 'ht_protocol_finder.pl' to identify the optimal threshold for HT-VS. If you execute the program without arguments, help is printed. The basic idea is that there is no need to perform exhaustive docking for every single molecule. We can start with a small number of runs (Genetic Algorithm optimizations). If the molecule reaches a reasonable docking score it will be worth continuing until completeness (50 docking runs). Otherwise, the molecule is unlikely to attain a good score and can be rejected earlier. 
This is how you execute it. To obtain a first guess of the values, you may take a look at the results table.

ht_protocol_finder.pl results_table.txt protocol.out -12 -15 5 15

After several trials, on of the predictions is:

 7.103, 21.880,  3.624,    5,  -15,   15,  -20 ***

This means that this 2-step process should be fairly efficient:
  • Will need 7.1% of the time needed to perfomr exhaustive docking
  • Will do 5 runs, and will continue if score reaches -15 (this will occur in 21.88% of cases)
  • Then complete 15 runs and will continue if a score fo -20 is reached (will happen in 3.6% of cases)
If you are satisfied with these conditions, you have to create a 'filter file' with this information. 

3)