How to analyse docking results from HADDOCK or refine models?
Best practice guide
There are various way how can your docking run be compared to a reference structure:
local HADDOCK scripts
If you run HADDOCK locally, you can use the scripts that were installed in the tools directory during HADDOCK installation. Such analysis is described in the local HADDOCK tutorial and in the HADDOCK2.4 manual These scripts can be also found here Github.
After using the HADDOCK webserver, one can analyse their results locally with our analysis scripts in the ana_scripts directory. This analysis is described in the HADDOCK ab-intio tutorial. However they need to be modified for each individual use. These scripts can be also found here Github.
DockQ estimates protein-protein model quality in a quantitative way using all parameters as the CAPRI classification e.g. fraction of native contacts (fcc), interface root-mean-square-deviation (i-RMSD), ligand RMSD (l-RMSD) DockQ can be found on Github in both python or fortran.
Clustering and scoring
A short section about key analysis parameters is written here.
The HADDOCK2.4 server provides a dedicated web interface to run a refinement on a molecular complex (still experimental). As input, only a PDB file for each partner of the complex is needed. In case one wants to tune the default parameters, it is possible to run the refinement also locally or using the regular submission interface. Then following settings need to be adjusted:
Settings to run water refinement locally
|Parameter||run.cns name||default value||optimal value|
|Center of mass restraints||
|Surface contact restraints||
|Number of structures for rigid body docking (it0)||
||1000||same as itw structures|
|Number of structures for semi-flexible refinement (it1)||
||200||same as itw structures|
|Sample 180 degrees rotated solutions during rigid body EM||
|Refine with short molecular dynamics in explicit solvent?||
|Multiply the number of calculated structures by all combination||
|Randomize starting orientations||
|Perform initial rigid body minimisation||
|Allow translation in rigid body minimisation||
|Number of MD steps for rigid body high temperature TAD||
|Number of MD steps during first rigid body cooling stage||
|Number of MD steps during second cooling stage with flexible side-chains at interface||
|Number of MD steps during third cooling stage with fully flexible interface||
* - only in json files, needs to be modified by hand in run.cns
Any more questions about analysis of the HADDOCK run? Have a look at our HADDOCK bioexcel forum hosted by . There is a very high chance that your problem has already been addressed.