Sampling Methods

From Rizzo_Lab
Revision as of 14:36, 23 February 2009 by Sudipto (talk | contribs)
(diff) ←Older revision | view current revision (diff) | Newer revision→ (diff)
Jump to: navigation, search

Sampling the correct bound pose is the first requirement for a successful docking simulation. No matter how good our scoring function is, we cannot report the correct bound pose unless we find it first.

Anchor and Grow

DOCK uses the Anchor and Growth method to sample the torsions in flexible molecules.

Sampling01.png

Clash Checking

During growth, some of the molecule conformations produced may be invalid. Clash checks removes unphysical poses. The clash_overlap parameter controls this. By default, the clash_penalty parameter is set to 0.5

  if nonbonded_atoms(a1, a2) 
    ref = clash_penalty * (atom_radius[a1] + atom_radius[a2])
    report clash if distance(a1,a2) < ref

Torsion Drive Libraries

Flexdrv01.png

max_orientations

The max number of orientations generated in the orienting code, before geometrical clustering.

min_anchor_size

Determines multiple anchor settings

pruning_max_orients

The number of anchor orients passed on to the anchor and grow code. A second clustering is performed in grow_periphery after minimization of the anchors which uses both energy and RMSD diversity. Presumably setting this equal to max_orientations turns off the geometric clustering in the orienting code.

pruning_clustering_cutoff

A parameter that roughly determines the number of conformers retained at each stage of growth after clustering. Default is 100. Greater numbers mean a greater # of conformers retained. Not a hard cutoff though, the number of conformers retained can be much greater at times.

pruning_conformer_score_cutoff

The max score (usually +ve) that will be allowed while pruning minimized conformers during growth. The same value is used for both grid score and internal energy. Typically 50 works well.