Effective handling of induced‐fit motion in flexible docking

For structure‐based drug design, where various ligand structures need to be docked to a target protein structure, a docking method that can handle conformational flexibility of not only the ligand, but also the protein, is indispensable. We have developed a simple and effective approach for dealing with the local induced‐fit motion of the target protein, and implemented it in our docking tool, ADAM. Our approach efficiently combines the following two strategies: a vdW‐offset grid in which the protein cavity is enlarged uniformly, and structure optimization allowing the motion of ligand and protein atoms. To examine the effectiveness of our approach, we performed docking validation studies, including redocking in 18 test cases and foreign‐docking, in which various ligands from foreign crystal structures of complexes are docked into a target protein structure, in 22 cases (on five target proteins). With the original ADAM, the correct docking modes (RMSD < 2.0 Å) were not present among the top 20 models in one case of redocking and four cases of foreign‐docking. When the handling of induced‐fit motion was implemented, the correct solutions were acquired in all 40 test cases. In foreign‐docking on thymidine kinase, the correct docking modes were obtained as the top‐ranked solutions for all 10 test ligands by our combinatorial approach, and this appears to be the best result ever reported with any docking tool. The results of docking validation have thus confirmed the effectiveness of our approach, which can provide reliable docking models even in the case of foreign‐docking, where conformational change of the target protein cannot be ignored. We expect that this approach will contribute substantially to actual drug design, including virtual screening. Proteins 2006. © 2006 Wiley‐Liss, Inc.

[1]  Richard A. Lewis,et al.  Lessons in molecular recognition: the effects of ligand and protein flexibility on molecular docking accuracy. , 2004, Journal of medicinal chemistry.

[2]  Thomas Lengauer,et al.  FlexE: efficient molecular docking considering protein structure variations. , 2001, Journal of molecular biology.

[3]  I. Kuntz,et al.  Molecular docking to ensembles of protein structures. , 1997, Journal of molecular biology.

[4]  Nobuo Tomioka,et al.  GREEN: A program package for docking studies in rational drug design , 1994, J. Comput. Aided Mol. Des..

[5]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[6]  I. Kuntz,et al.  Docking flexible ligands to macromolecular receptors by molecular shape. , 1986, Journal of medicinal chemistry.

[7]  Philip M Dean,et al.  Efficient conformational sampling of local side-chain flexibility. , 2003, Journal of molecular biology.

[8]  Akiko Itai,et al.  Application and evaluation of the automated docking method. , 1993 .

[9]  L. Kuhn,et al.  Virtual screening with solvation and ligand-induced complementarity , 2000 .

[10]  M. Mizutani,et al.  Efficient method for high-throughput virtual screening based on flexible docking: discovery of novel acetylcholinesterase inhibitors. , 2004, Journal of medicinal chemistry.

[11]  M. Mizutani,et al.  Rational automatic search method for stable docking models of protein and ligand. , 1994, Journal of molecular biology.

[12]  H. Wolfson,et al.  Access the most recent version at doi: 10.1110/ps.21302 References , 2001 .

[13]  H. Broughton,et al.  A method for including protein flexibility in protein-ligand docking: improving tools for database mining and virtual screening. , 2000, Journal of molecular graphics & modelling.

[14]  I. Kuntz,et al.  Conformational analysis of flexible ligands in macromolecular receptor sites , 1992 .

[15]  A. Itai,et al.  DEVELOPMENT OF AN EFFICIENT AUTOMATED DOCKING METHOD , 1993 .

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

[17]  M. Mizutani,et al.  Discovery of novel aldose reductase inhibitors using a protein structure-based approach: 3D-database search followed by design and synthesis. , 2001, Journal of medicinal chemistry.

[18]  B. Rupp,et al.  Structure of bovine pancreatic trypsin inhibitor at 125 K definition of carboxyl-terminal residues Gly57 and Ala58. , 1996, Acta crystallographica. Section D, Biological crystallography.

[19]  Thomas M Frimurer,et al.  Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities. , 2003, Biophysical journal.

[20]  Rafael Najmanovich,et al.  Side‐chain flexibility in proteins upon ligand binding , 2000, Proteins.

[21]  Didier Rognan,et al.  Comparative evaluation of eight docking tools for docking and virtual screening accuracy , 2004, Proteins.

[22]  D. Rognan,et al.  Protein-based virtual screening of chemical databases. 1. Evaluation of different docking/scoring combinations. , 2000, Journal of medicinal chemistry.

[23]  Robin Taylor,et al.  A new test set for validating predictions of protein–ligand interaction , 2002, Proteins.

[24]  A. Leach,et al.  Ligand docking to proteins with discrete side-chain flexibility. , 1994, Journal of molecular biology.

[25]  W Patrick Walters,et al.  A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance , 2004, Proteins.

[26]  J M Blaney,et al.  A geometric approach to macromolecule-ligand interactions. , 1982, Journal of molecular biology.

[27]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[28]  Frank Wien,et al.  Exploring the active site of herpes simplex virus type‐1 thymidine kinase by X‐ray crystallography of complexes with aciclovir and other ligands , 1998, Proteins.

[29]  R Nussinov,et al.  Flexible docking allowing induced fit in proteins: Insights from an open to closed conformational isomers , 1998, Proteins.

[30]  Heinz Sklenar,et al.  Harmonic modes as variables to approximately account for receptor flexibility in ligand-receptor docking simulations: Application to DNA minor groove ligand complex , 1999, J. Comput. Chem..

[31]  Charles L. Brooks,et al.  Detailed analysis of grid‐based molecular docking: A case study of CDOCKER—A CHARMm‐based MD docking algorithm , 2003, J. Comput. Chem..

[32]  R Abagyan,et al.  Flexible protein–ligand docking by global energy optimization in internal coordinates , 1997, Proteins.

[33]  J. Thornton,et al.  Conformational changes observed in enzyme crystal structures upon substrate binding. , 2005, Journal of molecular biology.

[34]  Nobuo Tomioka,et al.  A method for fast energy estimation and visualization of protein-ligand interaction , 1987, J. Comput. Aided Mol. Des..

[35]  Charles L. Brooks,et al.  Assessing energy functions for flexible docking , 1998, J. Comput. Chem..

[36]  David S. Goodsell,et al.  Distributed automated docking of flexible ligands to proteins: Parallel applications of AutoDock 2.4 , 1996, J. Comput. Aided Mol. Des..

[37]  Irwin D. Kuntz,et al.  Automated flexible ligand docking method and its application for database search , 1997 .

[38]  Ajay N. Jain Surflex: fully automatic flexible molecular docking using a molecular similarity-based search engine. , 2003, Journal of medicinal chemistry.

[39]  Maria Kontoyianni,et al.  Evaluation of docking performance: comparative data on docking algorithms. , 2004, Journal of medicinal chemistry.

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

[41]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[42]  Noriaki Hirayama,et al.  Ph4Dock: pharmacophore-based protein-ligand docking. , 2004, Journal of medicinal chemistry.

[43]  J. Irwin,et al.  Lead discovery using molecular docking. , 2002, Current opinion in chemical biology.

[44]  D. Goodsell,et al.  Automated docking to multiple target structures: Incorporation of protein mobility and structural water heterogeneity in AutoDock , 2002, Proteins.

[45]  D. Moras,et al.  A 'specificity' pocket inferred from the crystal structures of the complexes of aldose reductase with the pharmaceutically important inhibitors tolrestat and sorbinil. , 1997, Structure.

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

[47]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[48]  J. Tainer,et al.  Screening a peptidyl database for potential ligands to proteins with side‐chain flexibility , 1998, Proteins.

[49]  Amedeo Caflisch,et al.  Docking small ligands in flexible binding sites , 1998, J. Comput. Chem..

[50]  Brian K Shoichet,et al.  Testing a flexible-receptor docking algorithm in a model binding site. , 2004, Journal of molecular biology.

[51]  B. Shoichet,et al.  Soft docking and multiple receptor conformations in virtual screening. , 2004, Journal of medicinal chemistry.

[52]  I. Kuntz Structure-Based Strategies for Drug Design and Discovery , 1992, Science.

[53]  David A. Agard,et al.  The Structural Basis of Estrogen Receptor/Coactivator Recognition and the Antagonism of This Interaction by Tamoxifen , 1998, Cell.

[54]  Gisbert Schneider,et al.  Virtual screening and fast automated docking methods. , 2002, Drug discovery today.