The ligands were docked against the selected targets by using MVD 6.0. The scoring function used was MolDock Score (GRID) with a grid resolution of 0.30Å that gives a 4-5times enhancement by speeding up precalculating potential-energy values on an evenly spaced cubic grid. Tri-linear interpolation between the relevant grid points were considered for evaluating the energy potential. Ligand evaluation for each pose was performed by considering internal electrostatic interactions, internal hydrogen bonds and sp2-sp2 torsions into count for dihedral angles. The search algorithm
used was MolDock SE (Simplex Evolution) with 50 runs for a maximum of 1500 iterations and a population size of 50. The
poses were built incrementally by the pose generator from their rigid root point for the initial population size (=50) which tests a number of different torsions angles (either sp2-sp2, sp2-sp3 or sp3-sp3), rotations and translations, evaluation of the affected part of the molecule and chooses the value which results in the lowest energy contribution. The energy threshold was kept 100 for pose generation with a value of 10 as minimum, 30 for maximum and 10 for quick tries. A combined global/local search by the algorithm is done for a maximum of 300 steps and a Neighbor distance factor
of 0.1as the point of the initial simplex for random selection of individuals in the population [32]. Clustering of the poses were done by the Tabu clustering technique of MVD with an RMSD threshold of 2.00Å and an extra energy penalty of 100.00 which was performed by considering the intrinsic ligand symmetries.