Dataset Information

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets.

ABSTRACT: Virtual screening by molecular docking has become a widely used approach to lead discovery in the pharmaceutical industry when a high-resolution structure of the biological target of interest is available. The performance of three widely used docking programs (Glide, GOLD, and DOCK) for virtual database screening is studied when they are applied to the same protein target and ligand set. Comparisons of the docking programs and scoring functions using a large and diverse data set of pharmaceutically interesting targets and active compounds are carried out. We focus on the problem of docking and scoring flexible compounds which are sterically capable of docking into a rigid conformation of the receptor. The Glide XP methodology is shown to consistently yield enrichments superior to the two alternative methods, while GOLD outperforms DOCK on average. The study also shows that docking into multiple receptor structures can decrease the docking error in screening a diverse set of active compounds.

SUBMITTER: Zhou Z

PROVIDER: S-EPMC2547888 | biostudies-literature | 2007 Jul-Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets.

Zhou Zhiyong Z Felts Anthony K AK Friesner Richard A RA Levy Ronald M RM

Journal of chemical information and modeling 20070623 4

Virtual screening by molecular docking has become a widely used approach to lead discovery in the pharmaceutical industry when a high-resolution structure of the biological target of interest is available. The performance of three widely used docking programs (Glide, GOLD, and DOCK) for virtual database screening is studied when they are applied to the same protein target and ligand set. Comparisons of the docking programs and scoring functions using a large and diverse data set of pharmaceutica ...[more]

PMID: 17585856

Similar Datasets

Project description:Atropisomerism is a conformational chirality that occurs when there is hindered rotation about a σ-bond. While atropisomerism is exemplified by biaryls, it is observed in many other pharmaceutically relevant scaffolds including heterobiaryls, benzamides, diarylamines, and anilides. As bond rotation leads to racemization, atropisomers span the gamut of stereochemical stability. LaPlante has classified atropisomers based on their half-life of racemization at 37 °C: class 1 (t1/2 < 60 s), class 2 (60 s < t1/2 < 4.5 years), and class 3 (t1/2 > 4.5 years). In general, class-3 atropisomers are considered to be suitable for drug development. There are currently four FDA-approved drugs that exist as stable atropisomers, and many others are in clinical trials or have recently appeared in the drug discovery literature. Class-1 atropisomers are more prevalent, with ∼30% of recent FDA-approved small molecules possessing at least one class-1 axis. While class-1 atropisomers do not possess the requisite stereochemical stability to meet the classical definition of atropisomerism, they often bind a given target in a specific set of chiral conformations.Over the past decade, our laboratory has embarked on a research program aimed at leveraging atropisomerism as a design feature to improve the target selectivity of promiscuous lead compounds. Our studies initially focused on introducing class-3 atropisomerism into promiscuous kinase inhibitors, resulting in a proof of principle in which the different atropisomers of a compound can have different selectivity profiles with potentially improved target selectivity. This inspired a careful analysis of the binding conformations of diverse ligands bound to different target proteins, resulting in the realization that the sampled dihedral conformations about a prospective atropisomeric axis played a key role in target binding and that preorganizing the prospective atropisomeric axis into a desired target's preferred conformational range can lead to large gains in target selectivity.As atropisomerism is becoming more prevalent in modern drug discovery, there is an increasing need for strategies for atropisomerically pure samples of pharmaceutical compounds. This has led us and other groups to develop catalytic atroposelective methodologies toward pharmaceutically privileged scaffolds. Our laboratory has contributed examples of atroposelective methodologies toward heterobiaryl systems while also exploring the chirality of less-studied atropisomers such as diarylamines and related scaffolds.This Account will detail recent encounters with atropisomerism in medicinal chemistry and how atropisomerism has transitioned from a "lurking menace" into a leverageable design strategy in order to modulate various properties of biologically active small molecules. This Account will also discuss recent advances in atroposelective synthesis, with a focus on methodologies toward pharmaceutically privileged scaffolds. We predict that a better understanding of the effects of conformational restriction about a prospective atropisomeric axis on target binding will empower chemists to rapidly "program" the selectivity of a lead molecule toward a desired target.

Project description:BACKGROUND:In drug design, an efficient structure-based optimization of a ligand needs the precise knowledge of the protein-ligand interactions. In the absence of experimental information, docking programs are necessary for ligand positioning, and the choice of a reliable program is essential for the success of such an optimization. The performances of four popular docking programs, Gold, Glide, Surflex and FlexX, were investigated using 100 crystal structures of complexes taken from the Directory of Useful Decoys-Enhanced database. RESULTS:The ligand conformational sampling was rather efficient, with a correct pose found for a maximum of 84 complexes, obtained by Surflex. However, the ranking of the correct poses was not as efficient, with a maximum of 68 top-rank or 75 top-4 rank correct poses given by Glidescore. No relationship was found between either the sampling or the scoring performance of the four programs and the properties of either the targets or the small molecules, except for the number of ligand rotatable bonds. As well, no exploitable relationship was found between each program performance in docking and in virtual screening; a wrong top-rank pose may obtain a good score that allows it to be ranked among the most active compounds and vice versa. Also, to improve the results of docking, the strengths of the programs were combined either by using a rescoring procedure or the United Subset Consensus (USC). Oddly, positioning with Surflex and rescoring with Glidescore did not improve the results. However, USC based on docking allowed us to obtain a correct pose in the top-4 rank for 87 complexes. Finally, nine complexes were scrutinized, because a correct pose was found by at least one program but poorly ranked by all four programs. Contrarily to what was expected, except for one case, this was not due to weaknesses of the scoring functions. CONCLUSIONS:We conclude that the scoring functions should be improved to detect the correct poses, but sometimes their failure may be due to other varied considerations. To increase the chances of success, we recommend to use several programs and combine their results. Graphical abstract Summary of the results obtained by semi-rigid docking of crystallographic ligands. The docking was done on 100 protein-ligand X-ray structures, taken from the DUD-E database, and using four programs, Glide, Gold, Surflex and FlexX. Based on the docking results, we applied our United Subset Consensus method (USC), for which only the top4-rank poses are relevant. The number of complexes for which the best pose is correct, is represented by the gray boxes, the blue and red boxes correspond to the number of complexes with a correct pose ranked as the top 1 or within the top 4. A pose is considered correct when its root-mean-square deviation from the crystal structure is less than 2 Å.

Dataset Information

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets.

Publications

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets