Dock
Metaphorics Logo

Dock


 
Introduction
Overview
Installation
Quickstart
Hints and tips
Utilities
PVM version

usage: cat cmd.cex protein.cex spheres.cex molecules.cex | dock > docked.cex

Dock docks small molecules from "molecules.cex" into the protein from "protein.cex" in two stages:

  1. The small molecule's coordinates are generated by distance geometry into the spheres in "spheres.cex".

  2. The small molecule is optimized with a simple molecular mechanics function into the binding site in "protein.cex" which contains the spheres.

During the optimization, optional user specified distance constraints can be applied.

A reference conformation (e.g. the crystallographically determined bound ligand) is specified by adding coordinates for the molecule and a REMARK field of "1 LIGAND REFERENCE" (this is done automatically by pdb2site). The RMS difference between each docked conformation and the reference conformation will be reported, along with the reference conformation's score.

Including hydrogens (setting ADDH to TRUE) will typically at least double the amount of CPU time required, but hydrogens produce higher-quality ligand geometries, especially for molecules with flexible rings and aliphatic chains.

If the output CEX file is too large during large docking runs (for example, if it exceeds the 2GB file limit on Linux), the simplest solution is to pipe the output of dock into either gzip or split. The CEX output compresses well, so gzip is usually a good solution.

Dock options are specified in "cmd.cex", which is a simple text file. For now, the only way to change these options is to edit this file. The first section of "cmd.cex" defines CEX datatypes for each command option and should not be changed. Editable command options begin after $MSG<"dock options">. Frequently changed commands are:

VERBOSE

write summary output to stdout (TRUE/FALSE)

ADDH

generate the ligand(s) with all hydrogens included (TRUE/FALSE)

ZAPH

remove ligand hydrogens before output (TRUE/FALSE)

NUM_CONFORMATIONS

# dockings per ligand

BEST_CONFORMATIONS

# best dockings to output (0 = all)

SEED

random number seed. Set to RANDOM to get a different seed for each run

BUMP14

applies VDW bumping to 1-4 distances (TRUE/FALSE)

RANDCOORD

uses random coordinates method instead of metric matrix (TRUE/FALSE)

USE_REFERENCE_CONF

refine starting from the reference conformation (TRUE/FALSE)

SCORE_FUNCTION

scoring function which is used in docking refinement. Default is dockscore (the DockIt molecular mechanics score)

RESCORE_FUNCTIONS

scoring functions used in single point rescore. Default pmfscore, plpscore.

The following commands are rarely modified (an additional third stage of refinement is added if hydrogens are included using ADDH):

DO_SPHERES_FIRST

spheres are visible in the first stage (TRUE/FALSE)

DO_SPHERES_SECOND

spheres are visible in the second stage (TRUE/FALSE)

DO_SPHERES_THIRD

spheres are visible in the third stage (TRUE/FALSE)

DO_ENERGY_FIRST

site/ligand energy is included in the first stage (TRUE/FALSE)

DO_ENERGY_SECOND

site/ligand energy is included in the second stage (TRUE/FALSE)

DO_ENERGY_THIRD

site/ligand energy is included in the third stage (TRUE/FALSE)

TIGHT_SPHERE_SCORE

selects sphere subset by LIGSITE score for initial positioning of ligand by LIGSITE score

SPHERE_NEIGHBORS

# nearest spheres seen by each ligand atom

SPHERE_NEIGHBOR_UPDATE_RATE

frequency (# function calls) for updating sphere neighbors list

ENERGY_CUTOFF

distance threshold for selecting nearest site atoms seen by each ligand atom

SITE_NEIGHBOR_UPDATE_RATE

frequency (# function calls) for updating site atom nearest neighbors list

TRIALS

# attempts per docking before giving up

DIELECTRIC

dielectric constant

SITE_CUTOFF

maximum distance between protein atom and spheres' surface for including protein atom in docking site

MXDV

maximum allowed intra-ligand distance violation

MXVV

maximum allowed intra-ligand chiral violation

CONVERGENCE_1ST

gradient norm required for first stage

CONVERGENCE_2ND

gradient norm required for second stage

CONVERGENCE_3RD

gradient norm required for third stage

EVAL_LIMIT_1ST

max # function evaluations for first stage

EVAL_LIMIT_2ND

max # function evaluations for second stage

EVAL_LIMIT_3RD

max # function evaluations for third stage

SCALE_4D_1ST

4th-dimensional penalty function weight 1st stage

SCALE_4D_2ND

4th-dimensional penalty function weight 2nd stage

SCALE_4D_3RD

4th-dimensional penalty function weight 3rd stage

SCALE_SPHERE_1ST

sphere function weight for 1st stage

SCALE_SPHERE_2ND

sphere function weight for 2nd stage

SCALE_SPHERE_3RD

sphere function weight for 3rd stage

SCALE_ENERGY_1ST

energy function weight for 1st stage

SCALE_ENERGY_2ND

energy function weight for 2nd stage

SCALE_ENERGY_3RD

energy function weight for 3rd stage

© Metaphorics, LLC
info@metaphorics.com