Calculate DOCK6 RMSD

Jump to navigation Jump to search

Calculate ligand RMSDs using DOCK6

Calculating ligand RMSDs is a tricky thing, but DOCK6 does this quite well. All you need is two ligands in mol2 format (with the same number of atoms and same Sybil atom types) and a DOCK6 binary.

Then you run the following command:

~eidamo/work/scripts/calc_dock6_rmsd.csh xtal.mol2 dock.mol2

calc_dock6_rmsd.csh (see code below) outputs three lines:

HA_RMSDs:		4.2878 
HA_RMSDh: 2.7406
HA_RMSDm: 0.9625

The most meaningful line is the one in the middle (HA_RMSDh), which outputs the RMSD using the Hungarian Algorithm (for more information:


 #!/bin/csh -f
 if ($#argv != 2) then
         echo "Usage: $0 mol.mol2 ref.mol2"
         echo "calculates rmsd between 2 molecules using hungarian algorithm implemented in DOCK 6"
 setenv dock66_dir "/raid1/people/tbalius/zzz.programs/dock6_2012-10-09.stonybrook"
 echo "ligand_atom_file                                             $1" >
 echo "limit_max_ligands                                            no" >>
 echo "skip_molecule                                                no" >>
 echo "read_mol_solvation                                           no" >>
 echo "calculate_rmsd                                               yes" >>
 echo "use_rmsd_reference_mol                                       yes" >>
 echo "rmsd_reference_filename                                      $2" >>
 echo "use_database_filter                                          no" >>
 echo "orient_ligand                                                no" >>
 echo "use_internal_energy                                          no" >>
 echo "flexible_ligand                                              no" >>
 echo "bump_filter                                                  no" >>
 echo "score_molecules                                              no" >>
 echo "atom_model                                                   all" >>
 echo "vdw_defn_file                                                $dock66_dir/parameters/vdw_AMBER_parm99.defn" >>
 echo "flex_defn_file                                               $dock66_dir/parameters/flex.defn" >>
 echo "flex_drive_file                                              $dock66_dir/parameters/flex_drive.tbl" >>
 echo "ligand_outfile_prefix                                        output" >>
 echo "write_orientations                                           no" >>
 echo "num_scored_conformers                                        1" >>
 echo "rank_ligands                                                 no" >>
 $dock66_dir/bin/dock6 -i > dock6.log
 grep "HA_RMSD" output_scored.mol2


Thanks to Trent for pointing out this great new feature in DOCK 6.6.