[ Receptor Considerations | Ligand Considerations | Chirality Considerations | Protonation Considerations | Internal Hydrogen Bonds | Pyramidal Sampling | Membrane Exposed Docking | Waters | Scoring | Ligand Strain ]

For accurate ligand docking, the goal is to have an adequate three-dimensional model of the receptor pocket you are planning to dock ligands to. If this is the case then ICM docking has been shown to be very accurate in a number of independent assessments.

However, there are a number of pitfalls which need to be overcome to achieve accurate ligand docking. The pitfalls are that your model is not accurate overall, does not reflect the induced fit, or alternative conformations of the receptor binding pocket are missed. MolSoft has a number of in built tools to overcome some of these challenges

Some key points about ICM Ligand Docking:

• An average docking time is 0.1 seconds to 30 seconds per ligand per processor depending on the size of the ligand. The time per ligand was chosen to be the smallest possible to allow screening of very large data sets. To increase the time spent per ligand, change the thoroughness/effort parameter. For even faster docking, ideal for libraries more than a billion chemicals see our GPU enhanced methods such as RIDGE, GigaScreen and CombiRIDGE.
• ICM docking is the most accurate predictive tool of the binding geometry today. ICM docking has consistently ranked first place compared to other leading docking software in terms of accuracy. ICM has been successful in many drug design applications by scientists in academia and industry.
• ICM ligand docking procedure performs docking of the fully flexible small-molecule ligand to a known receptor 3D structure.
• The goal of the flexible docking calculation is prediction of correct binding geometry for each binder.
• ICM stochastic global optimization algorithm attempts to find the global minimum of the energy function that includes five grid potentials describing interaction of the flexible ligand with the receptor and internal conformational energy of the ligand. During the docking process a stack of alternative low energy conformations is saved.

Please click here to read more about ICM Docking.

• X-ray structure is preferable to an NMR structure.
• High resolution X-ray structure ( less than 2.1A ) is much better than, say 2.5A .
• Watch out for high-B-factor regions and avoid them; sometimes crystallographers deposit fantasy coordinates with high-B-factors.
• Check the pocket for correct placement of polar hydrogens and choose correct form of histidine.
• A bound conformation of the receptor is preferable, however if you use an apo-model, an NMR structure or a model by homology, the side-chains in a pocket may be incorrect. Frequently they stick out and prevent a ligand from binding. Those stubborn side-chains can be 'tamed', (i) manually; (ii) by a side chain simulation with elevated surfaceTension; or (iii) using a method such as Dual Alanine Scanning and Refinement ( SCARE).
• A model by homology can be built with the build model command (see molecular modeling section of this manual) and used for docking.

Usually a good place to start is to try to dock the known ligand(s) to the receptor model. You may also want to dock a library of compounds in order to identify lead candidates. In this case the main pitfalls are that the library is too restricted, molecules are not chemically feasible or not drug-like. For peptide docking please use the peptide docking protocol.

There is no need to convert the ligands to 3D, this is done "on-the-fly" during the docking process.

To sample linear chiral centers during docking you can check the option "Sample Racemic Centers" in Docking/Preferences/General. To sample chiral centers in rings it is best to generate all 2D enantiomers using Chemistry/Generate Stereoisomers. For most linear cases 'sample racemic centers' is acceptable. For most rigorous treatment one should pre-generate all stereoisomers (Chemistry/Generate Stereoisomers). If you do not select anything you get one (essentially arbitrary) stereoisomer docked.

The protonation state can be set in the Docking/Preferences/General dialog box. Choose the drop down option in the field labeled Charge Groups. To use MolSoft's pKa model choose the option auto, which uses built-in prediction of Ka/Kb to charge and protonate/deprotonate appropriate groups. If you do not want to use the pKa model then charges of some ionizable basic groups can be set, it currently accepts the following values: NH2, NH, NT for primary secondary and tertiary aliphatic amines. There are also options for imidazole and amidine.The none option still charges acids unless neutralAcids flag is set to no in the docking table file .dtb or .tab.

To determine protonation at a range of pH values you can use the options in Chemistry/Set Formal Charges/ and choose the option Auto using Pka Model.

Internal hydrogen bonds are taken into account during docking by default. You can check in the docking_project_name.dtb table file to make sure that l_internalHB field is 'yes'.

If the parameter l_samplePyramid is set to yes in the Docking Project .dtb file:

• During Monte Carlo sampling, the algorithm includes steps to invert the pyramid geometry at nitrogen atoms.
• This specifically targets uncharged aliphatic amines and sulfonamides.

(l_membraneScore)

If you are docking to a membrane-exposed site, set the l_membraneScore flag in the docking project .dtb file to yes.

Before docking, you will need to convert your PDB structure into an ICM object. During this step, you can choose how water molecules are handled. The options are to keep all waters, delete all waters, or remove only loose waters.

Loose waters are defined as water molecules that do not make any hydrogen bonds with the protein or a co-crystallized ligand. Tight waters are those that make three or more hydrogen bonds and are often more structurally stable.

Before docking, crystallographic water molecules in or near the binding pocket should be reviewed carefully. In most cases, bulk or weakly ordered waters should be deleted, as docking methods implicitly account for solvent effects and these waters may block ligand placement.

Waters that may be safely removed include:

• Waters located at the pocket entrance or solvent-exposed regions
• Waters with high B-factors or partial occupancy
• Waters not involved in stable hydrogen-bond networks
• Waters displaced by known ligands in co-crystal structures

Waters that may be worth retaining include:

• Conserved waters observed across multiple crystal structures
• Waters that bridge key protein-ligand interactions
• Waters forming stable hydrogen-bond networks with low B-factors
• Waters deeply buried in the pocket with limited exchange

When waters are retained, they should be treated consistently during docking and scoring. If the role of a water molecule is uncertain, it can be useful to run docking with and without selected waters and compare docking scores and poses.

• Physics-based Score. The Score is the ICM VLS Score which is a GBSA/MM-type scoring function augmented with a directional hydrogen bonding term (see Neves et al 2012). A ICM score -of 32 and lower are generally considered good scores - but it depends on the receptor (e.g. exposed pockets or pockets with metal ions may have higher scores than -32).

  An important parameter in a VLS run is the score threshold. A docked conformation for a given ligand is stored by the ICM VLS procedure only if its binding score is below this threshold.

  The threshold can be chosen in several ways:

  1. If a co-crystal structure is available, remove the ligand, re-dock it, and use the resulting score as a reference.
  2. If known ligands exist, dock them and examine their scores. A threshold set slightly above the typical scores of these ligands is usually a good starting point.
  3. If no ligands are known, run a short pre-simulation using approximately 1000 compounds from the screening database. Based on the score distribution, set the threshold to retain roughly 1 percent of the ligands.

  You can edit the Score threshold in the Docking/Preferences dialog box or edit the .dtb file.

• RTCNN Score is a Neural Network Score. RTCNN is a Radial Convolutional Neural Net including layers that do Topological (chemical graph) convolutions and 3D Radial convolutions. The RTCNN score is different from the regular ICM Score (see above) - it does not use any molecular mechanics or physical energy terms. Instead, the score is trained to recognize native-like complexes versus decoys directly, based only on geometries of putative complexes. So, for example, if internally the neural net recognizes a hydrogen bod, it is because it is observed during training that the configurations of atoms that we call hydrogen bonds are present to a greater extend in native-like structures. This is the same for other favorable or unfavourable interactions - it is all learnt from examples, not from any input terms.

  The RTCNN score performs extremely well on our benchmarks and the CASF benchmark (Su et al J.Chem.Inf.Model 2019). Please find details about this benchmark here: RTCNN_performance.pdf and read more about RTCNN Score here.

  When analyzing a hitlist from virtual screening we would recommend to plot Score vs RTCNN and prioritize those in the quadrant that have low Score and low RTCNN score. One key difference you may observe between RTCNN score and the regular Score is that the RTCNN score has a lower penalty for small induced fit ligand-receptor clashes. The lower the RTCNN score the better the predicted ligand-receptor interaction.

• mfScore is the statistical potential of mean force which was derived according to the published methodology (see Muegge J.Med.Chem 2006)

To calculate Strain, an ensemble of low-energy ligand conformers is generated in solution using the MMFF94s force field. The lowest-energy conformer is used as the unstrained reference. Strain is then defined as the energy difference between this reference conformation and the docked pose, with both energies computed using the same force field.

Analysis of specific torsion can be performed via torsion profile calculation at force-field level with MolMechanics/Torsion Scan, or via Crystallography Open Database search with appropriate substructure (Search Tab, select COD search, draw substructure) and statistics calculation with Chemistry/Torsion analysis.