Running DOCK

Starting a DOCK Run

You are now ready to run DOCK . Since dock can use a substantial amount of CPU time, it is a good idea to check whether there are other jobs running on the same machine. Consider any other users sharing your computers when deciding whether to start more than one run at a time. Be aware of any policies your site has regarding submitting background jobs.

The easiest way to select dock parameters is to run dock interactively. Do this by typing dock -i dock.in. You will be prompted for a value for each parameter. Any value you enter will be stored in the file dock.in. This file does not need to exist beforehand. If it does exist, then dock will extract all the relevant parameters it can find from the file. For each parameter, dock will supply a default value. If you want to use the default value, just hit return. The following tables list recommended values for running dock in two different ways: first to dock a single ligand, and second to dock a database of ligands. If you are viewing this manual on-line, then click on any of the keywords to view the reference entry for it.

Table 2. Recommended dock parameters for a new user.

2A: General

Keyword

Suggestions

flexible_ligand

no; try later

orient_ligand

yes; searches ligand orientations

score_ligand

yes; scores each ligand orientation

minimize_ligand

no; try later

multiple_ligands

no

random_seed

0 is fine

 

2B: Orientation Search

Keyword

Suggestions

match_receptor_sites

yes; matches ligand to site points

random_search

no

ligand_centers

no

automated_matching

yes; otherwise need to enter geometric match parameters

maximum_orientations

500; number of orientations to try

write_orientations

yes; to write out multiple orientations for single molecule

rank_orientations

yes; to save the top scoring orientations

rank_orientation_total

10; to save the top 10 orientations

 

2C: Scoring  

Keyword

Suggestions

intermolecular_score

yes

gridded_score

yes

grid_version

4

bump_filter

yes

bump_maximum

3

contact_score

no

chemical_score

no

energy_score

yes

atom_model

u; for united atom model

vdw_scale

1

electrostatic_scale

1

 

2D: Input  

Keyword

Suggestions

ligand_atom_file

Enter the ligand MOL2 file name here (including the directory path if this file is not in the current directory).

receptor_site_file

Enter the sphgen site point file name here.

score_grid_prefix

Enter the grid file name prefix here.

vdw_definition_file

~dock/parameter/vdw.defn

 

2E: Output

Keyword

Suggestions

ligand_energy_file

dock_nrg.mol2

Table 3. Recommended beginner's Dock parameters for a database search run

3A: Parameters to Modify from Database Run

Keyword

Previous

Suggestions

multiple_ligands

no

yes; to consider multiple molecules

maximum_orientations

500

50; just so that the run doesn't take too long

write_orientations

yes

no; to write only the best orientation for each molecule

 

3B: Multiple Ligands  

Keyword

Suggestions

parallel_jobs

no

ligands_maximum

first try 10 to make sure everything is working, then set it to a number larger than the number of molecules in input file

initial_skip

0; n>0 will skip the first n molecules in input file

interval_skip

0; n>0 will skip n molecules for each molecule processed

heavy_atoms_minimum

4; to discard small molecules

heavy_atoms_maximum

50; to discard large molecules

rank_ligands

yes; to save a top score list

rank_ligand_total

50; to save the top 50 molecules

restart_interval

100; to save restart info every 100th molecule processed

 

3C: Additional Input and Output

Keyword

Suggestions

quit_file

dock.quit

dump_file

dock.dump

info_file

dock.info

restart_file

dock.rst

If you happen to enter the wrong value for any parameter and wish to change it, then you may edit the dock.in file directly and modify the parameter value. Once all parameters have been entered, dock should begin the calculation and the dock.in file is complete. You may kill the process with a CTRL-C and relaunch the process in background by typing "dock -i dock.in -o dock.out&". If you would like to run dock multiple times from the same directory, then you may use a different name for dock.in and dock.out. Just be sure to change the names of the output files inside the new dock.in file so that two processes don't end up overwriting each other's output files.

Check a few minutes after you start the run to be sure that it is still going; if it has stopped, look for mistakes in the input. Beginners should check disk usage occasionally while the job is running, just in case the program is creating incredibly large files which might overflow the available space.

During a database search run (which can take anywhere from hours to days to weeks to finish), you can follow dock's progress through the database by inspecting the *.info file.

Restarting a Search Run

In database search mode of dock (when multiple_ligands, orient_ligand, and rank_ligands are selected), dock periodically saves information necessary to restart the search from its current location in the database in a *.rst file. If there is a power failure or the system crashes, you can set up a new run to start where the last one was stopped. First, make a copy of dock.out so that status of the previous run are saved. Then relaunch the job using the -r flag at the command line. (Do not change the remaining files, since dock needs them to restart successfully.) When the restarted run finishes, the sorted list of ligands in the output file will include the top scorers from the entire database.

Looking at the Results

Dock puts the resulting molecule orientations in a file for each type of scoring function used. The scores are given in the comment records at the beginning of each molecule record. If you have selected MOL2 format for your output files and your graphical viewer does not read this format, then convert the file to PDB by typing dock -i dock_nrg.mol2 -o dock_nrg.pdb.

Other Post-Docking Tasks

Depending on your particular project, you might be interesting in any one or several of the following post-docking techniques:

 

© UC Regents 1998
Top / Up / Previous / Next