This paper is devoted to results obtained by the docking program SOL and the post-processing program DISCORE at the CSAR benchmark. SOL and DISCORE programs are described. SOL is the original docking program developed on the basis of the genetic algorithm, MMFF94 force field, rigid protein, precalculated energy grid including desolvation in the frame of simplified GB model, vdW, and electrostatic interactions and taking into account the ligand internal strain energy. An important SOL feature is the single- or multi-processor performance for up to hundreds of CPUs. DISCORE improves the binding energy scoring by the local energy optimization of the ligand docked pose and a simple linear regression on the base of available experimental data. The docking program SOL has demonstrated a good ability for correct ligand positioning in the active sites of the tested proteins in most cases of CSAR exercises. SOL and DISCORE have not demonstrated very exciting results on the protein-ligand binding free energy estimation. Nevertheless, for some target proteins, SOL and DISCORE were among the first in prediction of inhibition activity. Ways to improve SOL and DISCORE are discussed.

The findings of this review confirm the great demand in docking programs for discovery of new medicine substances with the help of molecular modeling. New trends in docking programs design are revealed. These trends are focused on the increase of the docking accuracy at the expense of more accurate molecular energy calculations without any fitting parameters, including quantum-chemical methods and implicit solvent models, and by using new global optimization algorithms which make it possible to treat flexibility of ligands and mobility of protein atoms simultaneously. Finally, it is shown that all the necessary prerequisites for increasing the docking accuracy can be accomplished in practice.

The adequate choice of the docking target function impacts the accuracy of the ligand positioning as well as the accuracy of the protein-ligand binding energy calculation. To evaluate a docking target function we compared positions of its minima with the experimentally known pose of the ligand in the protein active site. We evaluated five docking target functions based on either the MMFF94 force field or the PM7 quantum-chemical method with or without implicit solvent models: PCM, COSMO, and SGB. Each function was tested on the same set of 16 protein-ligand complexes. For exhaustive low-energy minima search the novel MPI parallelized docking program FLM and large supercomputer resources were used. Protein-ligand binding energies calculated using low-energy minima were compared with experimental values. It was demonstrated that the docking target function on the base of the MMFF94 force field in vacuo can be used for discovery of native or near native ligand positions by finding the low-energy local minima spectrum of the target function. The importance of solute-solvent interaction for the correct ligand positioning is demonstrated. It is shown that docking accuracy can be improved by replacement of the MMFF94 force field by the new semiempirical quantum-chemical PM7 method.

No abstract

Results of the combined use of the classical force field and the recent quantum chemical PM7 method for docking are presented. Initially the gridless docking of a flexible low molecular weight ligand into the rigid target protein is performed with the energy function calculated in the MMFF94 force field with implicit water solvent in the PCM model. Among several hundred thousand local minima, which are found in the docking procedure, about eight thousand lowest energy minima are chosen and then energies of these minima are recalculated with the recent quantum chemical semiempirical PM7 method. This procedure is applied to 16 test complexes with different proteins and ligands. For almost all test complexes such energy recalculation results in the global energy minimum configuration corresponding to the ligand pose near the native ligand position in the crystalized protein-ligand complex. A significant improvement of the ligand positioning accuracy comparing with MMFF94 energy calculations is demonstrated.