Prev | 23.6 FAQ-Docking | Next |
[ Energy Units | Calculate Binding Energy | Faq ringflex dock | Faq guided dock | Faq reload dock | Faq dock repeat | SCORE | Faq hitlist | Faq score | Ligand and box | Pocket Finder | Dock Time | Thoroughness | Probe | Faq receptor selection | Flexible Receptor Dock | Dock Command Line | Faq scanScoreExternal | Faq scanScoreExternal2 | Background Job | flexible rings | Dock Racemic | Dock charge groups | Working directory ]
Frequently asked questions regarding small molecule and protein-protein docking.
23.6.1 What are the units of the energy values displayed after docking? |
The energy units are kcal/mol
23.6.2 I do not have ICM-VLS but I would like to calculate the binding energy of my docked complex - how can I do this? |
For more in depth information on this topic please see the command line manual www.molsoft.com/man however the basic approach is this: Calculuate the energy of the receptor (e1, a_1), energy of the ligand (e2, a_2) and the energy of the complex (e12) then the binding energy = e12 - e1 - e2
Here is a script to do this whereby the ligand is a_2 and the receptor is a_1:
electroMethod="boundary element" surfaceMethod="constant tension" surfaceTension=0.020 dielConst = 12.7 set terms "sf,el,en" read object s_icmhome+"2ptc" show energy a_1 a_1 mute e1 =Energy("el,sf,en") show energy a_2 a_2 mute e2 =Energy("el,sf,en") show energy mute e12 =Energy("el,sf,en") print "Binding energy = ", e12 - e1 - e2
There are many different approaches to the evaluation of binding energy. One of the reasonable approximations has the following features: van der Waals/hydrogen bonding interaction is excluded since it has close magnitudes for protein-protein and for protein-solvent interactions; electrostatic free energy change is calculated by the REBEL method (see also the section "How to calculate the electrostatic free energy ... ") above); side-chain entropy change is calculated by standard ICM entropic term based on exposed surface area of flexible side-chains; hydrophobic energy change is calculated using surface term with constant surface tension of 20. cal/Angstrom.
23.6.3 How do I sample conformations of flexible rings in docked ligands, for example, a "chair-boat" transition? |
Find your docking project file ( yourDockingProjectName.tab ) and set the ringFlexibilityLevel to 1 or 2.
In GUI, this parameter can be reset in the "Docking Preferences" dialog. It is called "Flexible ring sampling level".
23.6.4 How can I guide my docking to a known conformation of a smilar ligand? |
Use Docking/Template as described here.
23.6.5 How do I reload a docking project? |
To reload a docking project.
/Docking/Set Project - Type in the Docking Project Name (Case Sensitive)
Now you can browse scan solutions etc.... and use the maps to dock another ligand.
23.6.6 In a VLS run how many times should I run the docking? |
Generally we suggest the docking should be repeated 2-3 times and the lowest ICM score pose should be taken.
23.6.7 Which score value should I use for analysis. |
The value under the heading SCORE relates to the ICM docking score and is the best one to use for docking result analysis. The other score we provide - potential of mean force score (mfscore - http://www.molsoft.com/man/terms.html#term-mf) provides an independent score of the strength of ligand-receptor interaction.
23.6.8 Some compounds are missing from my HITLIST. |
The hitlist is filtered according to a score cutoff defined in the PROJECTNAME.TAB file. Therefore poor scoring compounds are not reported in the HITLIST - this can be changed by opening the .tab file in a suitable text editor e.g. notepad in windows or vi and changing the DOCK1.r_ScoreThreshold value (by default it is -32). The scores for all compounds in a VLS screen are available in the PROJECTNAME.OU file.
You can also change this value in Graphical User Interface:
23.6.9 What constitutes a good docking score? |
Generally a score below -32 is regarded as a good docking score. A good score depends on the system into which you are docking. For example is the pocket open or closed and are there metal ions interacting with the ligand. If the pocket is open scores higher than -32 may indicate potential binders. If a crystal structure is available remove the ligand and re-dock it to get an indication of approximately what is a good score for the receptor you are interested in.
23.6.10 When I view my docking run my ligand never jumps into the box - what did I do wrong? |
Here are some reasons and some solutions for why your ligand is sampling outside of the binding pocket:
23.6.11 How do I identify the binding pockets in my receptor? |
To do this go to:
More information can be found in the section entitled Identifying Binding Pockets
23.6.12 How long does it take to dock one ligand using ICM-VLS? |
It takes approximately 30-60 seconds per ligand depending on the size of the ligand and the nature of the pocket. ICM ranked first place compared to other leading docking software in terms of accuracy in a recent analysis undertaken by Astra Zeneca scientists.
See: On Evaluating Molecular-Docking Methods for Pose Prediction and Enrichment Factors Hongming Chen, Paul D. Lyne, Fabrizio Giordanetto, Timothy Lovell, and Jin Li J. Chem. Inf. Model.; 2006; 46(1) pp 401 - 415
23.6.13 What does thoroughness mean? |
When you send a docking job either using Docking/Interactive or Docking/Run Docking Batch you are asked to enter a thoroughness value. This value represents the length of the docking simulation. The default value is 1. and this works well with nearly every kind of docking scenario, however in certain circumstances such as if you have a very large pocket this value should be increased slightly to a range of between 5. and 10.
23.6.14 When I setup the receptor I am asked to move a probe - what is this? |
The probe which you see after Receptor Setup (Docking/ Receptor Setup) represents the initial starting position for the ligand. Usually this does not have to be changed as ICM by default places it into the center of the pocket. However if you do wish to move it to position closer to a critical region of the pocket you can do this using the middle mouse button when prompted in receptor setup.
23.6.15 I want to dock to the receptor and include other molecules in the receptor such as a tightly bound water molecule - how can I do this? |
All molecules that you wish to dock to need to be stated in the Docking/Receptor Setup/ Receptor molecule (s) data entry box using the ICM selection language. For example if you wish to dock to the protein with PDB code 1m17 and water molecule number 20 you need to enter the information as shown below:
23.6.16 How can I run docking with a flexible receptor? |
23.6.17 How can I run the docking simulation from the command line? |
Use dockScan as described here: http://www.molsoft.com/man/ligand-fit.html#dockScan
23.6.18 I have a complex I wish to generate an ICM VLS Score for, however I did not dock it using VLS. How can I do this? |
type in the command line:
scanScoreExternal
OR
23.6.19 Why is there always a small difference between the score calculated interactively by scanScoreExternal and that obtained by docking (VLS)? |
The reason for this is that ICM score has terms that require calculations on the reference 'free' state of the ligand, in particular solvation electrostatics and internal force-field strain energy change are calculated as a difference between free and bound state. VLS uses as a free state the lowest energy conformation found by MC search for the unbound ligand. Interactive score just minimizes the ligand locally. To ensure consistency we recomend you use one method or the other for scoring or you could recalculate the interactive score for your ligand from VLS before modifying/minimizing.
23.6.20 How do I monitor and terminate a background docking job? |
If a background job is running you will see a message saying "1 bgrnd job" at the top of the gui interface (gui blue title panel).
To monitor the progress of the job:
To terminate a background job:
To view the current output of a background job:
When a background job has finished a message will appear in the graphical user interface
23.6.21 How do I sample flexible ring conformations (boat, chair etc..) during docking? |
MolSoft's ICM docking algorithm has flexible ring sampling included on the fly. Just set ring sampling level to 1 (flex ring only in pre-sampling step) or 2 (throughout the simulation).
To do this:
OR,
If you want to generate the conformations before docking and you have ICM-Pro + ICM-Chemistry then you can use the conformation generator algorithm described here: http://www.molsoft.com/gui/conf-gen.html
23.6.22 I am docking a racemic compound how can I sample both R and S states during docking? |
To sample both R and S states of a compound during docking. Edit the project_name.tab file and edit l_sampleRacemin to yes
If l_sampleRacemic is 'yes', R and S states are sampled for racemic centers and best-fitting one is chosen. If it is 'no', they are kept fixed (in an aribtrary R or S state). Note that stereo centers that centers with pre-assigned R or S state are never sampled, if sampling is desired they need to be reassigned as racemic.
23.6.23 What do the options mean in Docking/Preferences/Charge Groups? |
This relates only to the ligands:
In the very latest version of ICM (3.7-2a) we have a prediction of pKA which will automatically charge each ligand according to the pKA prediction. If you want to change the charge of the residues on the receptor please right click on the ligand select mutate and their are uncharged versions of each residue. If you do this then you would have to re-make the maps.
23.6.24 How do I change the working directory? |
This is a common question from people who use Windows especially for docking. By default all the docking files are stored in a directory in “My Documents”. This may not be convenient and so the way to change this is to use the command line and type: set directory "c:/path_to_directory/"' In the very latest version of the software there is now a GUI button to do this in the Tools menu.
Prev Closest Neighbor | Home Up | Next FAQ-Cheminformatics |