A pharmacophore docking algorithm and its application to the cross‐docking of 18 HIV‐NNRTI's in their binding pockets

The docking of small molecules into the binding site of a target protein is an important but difficult step in structure‐based drug design. The performance of a docking algorithm is usually evaluated by re‐docking ligands into their native binding sites. We have explored the cross‐docking of 18 HIV‐NNRTIs (non‐nucleoside inhibitors of HIV reverse transcriptase) of which the ligand–protein structure has been determined: each of the 18 ligands was docked into each of the 18 binding sites. The docking algorithms studied are an energy‐based simulated annealing algorithm and a novel pharmacophore docking algorithm. It turns out that the energy‐based docking of the ligands into non‐native pockets is far less successful than the docking into their native pockets. The results can be improved by using explicit pharmacophore information, and by docking a ligand into a panel of protein structures and selecting the ligand–proteincombination with the lowest interaction energy as the final result. Proteins 2004;54:000–000. © 2003 Wiley‐Liss, Inc.

[1]  Eiji Osawa,et al.  Corner flapping: a simple and fast algorithm for exhaustive generation of ring conformations , 1989 .

[2]  D. Stuart,et al.  The structure of HIV-1 reverse transcriptase complexed with 9-chloro-TIBO: lessons for inhibitor design. , 1995, Structure.

[3]  Todd J. A. Ewing,et al.  Critical evaluation of search algorithms for automated molecular docking and database screening , 1997, J. Comput. Chem..

[4]  W. Kabsch.  A discussion of the solution for the best rotation to relate two sets of vectors , 1978 .

[5]  D I Stuart,et al.  Crystal structures of HIV-1 reverse transcriptase in complex with carboxanilide derivatives. , 1998, Biochemistry.

[6]  Thomas Lengauer,et al.  Flexible docking under pharmacophore type constraints , 2002, J. Comput. Aided Mol. Des..

[7]  D I Stuart,et al.  Structural basis for the resilience of efavirenz (DMP-266) to drug resistance mutations in HIV-1 reverse transcriptase. , 2000, Structure.

[8]  Thomas A. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[9]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[10]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[11]  Richard T. Walker,et al.  Complexes of HIV-1 reverse transcriptase with inhibitors of the HEPT series reveal conformational changes relevant to the design of potent non-nucleoside inhibitors. , 1996, Journal of medicinal chemistry.

[12]  R Abagyan,et al.  High-throughput docking for lead generation. , 2001, Current opinion in chemical biology.

[13]  Jonathan W. Essex,et al.  A review of protein-small molecule docking methods , 2002, J. Comput. Aided Mol. Des..

[14]  Henri Moereels,et al.  Structure of HIV-1 RT/TIBO R 86183 complex reveals similarity in the binding of diverse nonnucleoside inhibitors , 1995, Nature Structural Biology.

[15]  A. D. Clark,et al.  Structures of Tyr188Leu mutant and wild-type HIV-1 reverse transcriptase complexed with the non-nucleoside inhibitor HBY 097: inhibitor flexibility is a useful design feature for reducing drug resistance. , 1998, Journal of molecular biology.

[16]  Jan Balzarini,et al.  Phenylethylthiazolylthiourea (PETT) Non-nucleoside Inhibitors of HIV-1 and HIV-2 Reverse Transcriptases , 2000, The Journal of Biological Chemistry.

[17]  Steven G. Louie,et al.  A Monte carlo simulated annealing approach to optimization over continuous variables , 1984 .

[18]  Eddy Arnold,et al.  On the detection of multiple-binding modes of ligands to proteins, from biological, structural, and modeling data , 2003, J. Comput. Aided Mol. Des..

[19]  D I Stuart,et al.  Binding of the Second Generation Non-nucleoside Inhibitor S-1153 to HIV-1 Reverse Transcriptase Involves Extensive Main Chain Hydrogen Bonding* , 2000, The Journal of Biological Chemistry.

[20]  J. Gasteiger,et al.  FROM ATOMS AND BONDS TO THREE-DIMENSIONAL ATOMIC COORDINATES : AUTOMATIC MODEL BUILDERS , 1993 .

[21]  Brendan Larder,et al.  A Rapid Method for Simultaneous Detection of Phenotypic Resistance to Inhibitors of Protease and Reverse Transcriptase in Recombinant Human Immunodeficiency Virus Type 1 Isolates from Patients Treated with Antiretroviral Drugs , 1998, Antimicrobial Agents and Chemotherapy.

[22]  Shuichi Hirono,et al.  Estimation of active conformations of drugs by a new molecular superposing procedure , 1999, J. Comput. Aided Mol. Des..

[23]  Thomas Lengauer,et al.  Evaluation of the FLEXX incremental construction algorithm for protein–ligand docking , 1999, Proteins.

[24]  Pieter F. W. Stouten,et al.  A molecular mechanics/grid method for evaluation of ligand–receptor interactions , 1995, J. Comput. Chem..

[25]  Didier Rognan,et al.  ConsDock: A new program for the consensus analysis of protein–ligand interactions , 2002, Proteins.

[26]  Gennady M Verkhivker,et al.  Exploring the energy landscapes of molecular recognition by a genetic algorithm: analysis of the requirements for robust docking of HIV-1 protease and FKBP-12 complexes. , 1996, Proteins.

[27]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[28]  D I Stuart,et al.  Design of MKC-442 (emivirine) analogues with improved activity against drug-resistant HIV mutants. , 1999, Journal of medicinal chemistry.

[29]  X Fradera,et al.  Similarity‐driven flexible ligand docking , 2000, Proteins.

[30]  Thomas Lengauer,et al.  A fast flexible docking method using an incremental construction algorithm. , 1996, Journal of molecular biology.

[31]  H. Vahaboğlu,et al.  Activities of cefepime and five other antibiotics against nosocomial PER-1-type and/or OXA-10-type beta-lactamase-producing Pseudomonas aeruginosa and Acinetobacter spp. , 1998, The Journal of antimicrobial chemotherapy.

[32]  Yvonne Jones,et al.  High resolution structures of HIV-1 RT from four RT–inhibitor complexes , 1995, Nature Structural Biology.

[33]  Gennady M Verkhivker,et al.  Exploring the energy landscapes of molecular recognition by a genetic algorithm: Analysis of the requirements for robust docking of HIV‐1 protease and FKBP‐12 complexes , 1996 .