# Difference between revisions of "DOCK3.7 INDOCK Minimization Parameter"

Enkhjargal (talk | contribs) |
|||

(6 intermediate revisions by 2 users not shown) | |||

Line 22: | Line 22: | ||

sim_itmax 500 | sim_itmax 500 | ||

− | This is the initial distance in angstroms the molecule is translated. | + | This is the initial distance in angstroms the molecule is translated (In Angstrom). (This is used to initiliaze the simplex) |

sim_trnstep 0.2 | sim_trnstep 0.2 | ||

− | How many degrees of initial rotation are done. | + | How many degrees of initial rotation are done (In Degree). (This is used to initiliaze the simplex) |

sim_rotstep 5.0 | sim_rotstep 5.0 | ||

− | If the energy changes by this much, restart the minimizer from this newest position. | + | If the energy changes by this much, restart the minimizer from this newest position. Between 2 subseq runs: if energy changes by more than 1kcal restart minimizer from latest position for another round. |

sim_need_to_restart 1.0 | sim_need_to_restart 1.0 | ||

− | How much the total energy can changed to be considered converged. Setting this higher will stop faster, setting it lower will cause it to do more iterations before converging (or potentially hitting the iteration max above). | + | How much the total energy can changed to be considered converged. within one run: considered converged when energies of lowest and highest energy state differ by only 0.1kcal/mol. Setting this higher will stop faster, setting it lower will cause it to do more iterations before converging (or potentially hitting the iteration max above). |

sim_cnvrge 0.1 | sim_cnvrge 0.1 | ||

− | Don't minimize molecules that score above | + | Don't minimize molecules that score above this large number. |

min_cut 1.0e15 | min_cut 1.0e15 | ||

− | To initiliaze the simplex, this | + | To initiliaze the simplex, this iseed number used to generate a random number. e.g. used for trn and rotstep calc |

iseed 777 | iseed 777 | ||

+ | |||

+ | |||

+ | == Information on number_save and number_write parameters == | ||

+ | |||

+ | This is true for DOCK v3.7.1rc1 (and later versions). | ||

+ | |||

+ | "number_save" is the number (for each hierarchy) that are minimized and written to OUTDOCK. | ||

+ | |||

+ | "number_write" are the number (for each hierarchy) that are written to the mol2 file. | ||

+ | |||

+ | For example, if you set number_save to 10 and number_write to 2 then 10 poses will be minimized but only the top two scoring minimized poses will be written out to the mol2 file. | ||

+ | |||

+ | |||

+ | The parameters (number_save/number_write) are applied to each db2 hierarchy. Each protomer/tautomer has one or more db2 hierarchies. So if there are multiple protonation state each will have the "number_write" poses written out to the mol2 file (caveat, dock might not generate that many in which case it will write all that are possible). | ||

+ | |||

+ | We recommended setting these numbers to 1 for production run discovery campaigns. | ||

+ | we topically only modify number_save/number_write to 10 if you want to perform more exhaustive sampling for a small number (e.g. less than 1000) of molecules. | ||

+ | |||

+ | [[File:Miminizer.2017.11.14.png|thumb|center|500px|how resulting poses will very with different values of number_save and number_write]] |

## Latest revision as of 00:19, 15 November 2017

The best scored ligand conformation can be minimized using the following parameters.

Add following section to existing INDOCK after SCORING section.

##################################################### # MINIMIZATION minimize yes sim_itmax 500 sim_trnstep 0.2 sim_rotstep 5.0 sim_need_to_restart 1.0 sim_cnvrge 0.1 min_cut 1.0e15 iseed 777

The default for minimize is no. Yes turns on the minimization.

minimize yes

How many iterations of minimization to do. More means longer run times, but potentially better poses.

sim_itmax 500

This is the initial distance in angstroms the molecule is translated (In Angstrom). (This is used to initiliaze the simplex)

sim_trnstep 0.2

How many degrees of initial rotation are done (In Degree). (This is used to initiliaze the simplex)

sim_rotstep 5.0

If the energy changes by this much, restart the minimizer from this newest position. Between 2 subseq runs: if energy changes by more than 1kcal restart minimizer from latest position for another round.

sim_need_to_restart 1.0

How much the total energy can changed to be considered converged. within one run: considered converged when energies of lowest and highest energy state differ by only 0.1kcal/mol. Setting this higher will stop faster, setting it lower will cause it to do more iterations before converging (or potentially hitting the iteration max above).

sim_cnvrge 0.1

Don't minimize molecules that score above this large number.

min_cut 1.0e15

To initiliaze the simplex, this iseed number used to generate a random number. e.g. used for trn and rotstep calc

iseed 777

## Information on number_save and number_write parameters

This is true for DOCK v3.7.1rc1 (and later versions).

"number_save" is the number (for each hierarchy) that are minimized and written to OUTDOCK.

"number_write" are the number (for each hierarchy) that are written to the mol2 file.

For example, if you set number_save to 10 and number_write to 2 then 10 poses will be minimized but only the top two scoring minimized poses will be written out to the mol2 file.

The parameters (number_save/number_write) are applied to each db2 hierarchy. Each protomer/tautomer has one or more db2 hierarchies. So if there are multiple protonation state each will have the "number_write" poses written out to the mol2 file (caveat, dock might not generate that many in which case it will write all that are possible).

We recommended setting these numbers to 1 for production run discovery campaigns. we topically only modify number_save/number_write to 10 if you want to perform more exhaustive sampling for a small number (e.g. less than 1000) of molecules.