A fast flexible docking method using an incremental construction algorithm.

We present an automatic method for docking organic ligands into protein binding sites. The method can be used in the design process of specific protein ligands. It combines an appropriate model of the physico-chemical properties of the docked molecules with efficient methods for sampling the conformational space of the ligand. If the ligand is flexible, it can adopt a large variety of different conformations. Each such minimum in conformational space presents a potential candidate for the conformation of the ligand in the complexed state. Our docking method samples the conformation space of the ligand on the basis of a discrete model and uses a tree-search technique for placing the ligand incrementally into the active site. For placing the first fragment of the ligand into the protein, we use hashing techniques adapted from computer vision. The incremental construction algorithm is based on a greedy strategy combined with efficient methods for overlap detection and for the search of new interactions. We present results on 19 complexes of which the binding geometry has been crystallographically determined. All considered ligands are docked in at most three minutes on a current workstation. The experimentally observed binding mode of the ligand is reproduced with 0.5 to 1.2 A rms deviation. It is almost always found among the highest-ranking conformations computed.

[1]  E. Shaw,et al.  STUDIES ON THE ACTIVE CENTER OF TRYPSIN. THE BINDING OF AMIDINES AND GUANIDINES AS MODELS OF THE SUBSTRATE SIDE CHAIN. , 1965, The Journal of biological chemistry.

[2]  H A Scheraga,et al.  CONFORMATIONAL ENERGY CALCULATIONS OF ENZYME‐SUBSTRATE INTERACTIONS. I. Computation of Preferred Conformations of Some Substrates of α‐Chymotrypsin , 2009 .

[3]  Richard O. Duda,et al.  Pattern classification and scene analysis , 1974, A Wiley-Interscience publication.

[4]  POSSIBLE TRANSITION‐STATE ANALOGS FOR RIBONUCLEASE, THE COMPLEXES OF URIDINE WITH OXOVANADIUM(IV) ION AND VANADIUM(V) ION , 1974 .

[5]  John W. Bunting,et al.  Reversible Inhibition of Carboxypeptidase A. IV. Inhibition of Specific Esterase Activity by Hippuric Acid and Related Species and other Amino Acid Derivatives and a Comparison with Substrate Inhibition , 1975 .

[6]  M G Rossmann,et al.  A comparison of the structures of apo dogfish M4 lactate dehydrogenase and its ternary complexes. , 1977, Journal of molecular biology.

[7]  W. Kabsch.  A solution for the best rotation to relate two sets of vectors , 1976 .

[8]  V Massey,et al.  Flavin-oxygen derivatives involved in hydroxylation by p-hydroxybenzoate hydroxylase. , 1976, The Journal of biological chemistry.

[9]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[10]  R. Huber,et al.  The transition of bovine trypsinogen to a trypsin-like state upon strong ligand binding. The refined crystal structures of the bovine trypsinogen-pancreatic trypsin inhibitor complex and of its ternary complex with Ile-Val at 1.9 A resolution. , 1978, Journal of molecular biology.

[11]  F. Allen,et al.  The Cambridge Crystallographic Data Centre: computer-based search, retrieval, analysis and display of information , 1979 .

[12]  W. Bode,et al.  The transition of bovine trypsinogen to a trypsin-like state upon strong ligand binding. II. The binding of the pancreatic trypsin inhibitor and of isoleucine-valine and of sequentially related peptides to trypsinogen and to p-guanidinobenzoate-trypsinogen. , 1979, Journal of molecular biology.

[13]  B. Matthews,et al.  Binding of hydroxamic acid inhibitors to crystalline thermolysin suggests a pentacoordinate zinc intermediate in catalysis. , 1982, Biochemistry.

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

[15]  J. Bolin,et al.  Crystal structures of Escherichia coli and Lactobacillus casei dihydrofolate reductase refined at 1.7 A resolution. I. General features and binding of methotrexate. , 1982, The Journal of biological chemistry.

[16]  D C Rees,et al.  Crystallographic studies on apocarboxypeptidase A and the complex with glycyl-L-tyrosine. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[17]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[18]  R. Huber,et al.  The Geometry of the Reactive Site and of the Peptide Groups in Trypsin, Trypsinogen and its Complexes with Inhibitors , 1983 .

[19]  B. Matthews,et al.  Binding of N-carboxymethyl dipeptide inhibitors to thermolysin determined by X-ray crystallography: a novel class of transition-state analogues for zinc peptidases. , 1984, Biochemistry.

[20]  P. D. Clercq,et al.  Systematic conformational analysis. A microcomputer method for the semiquantitative evaluation of polycyclic systems containing five-, six- and seven-membered rings. 1. Program characteristics. , 1984 .

[21]  Systematic conformational analysis. A microcomputer method for the semiquantitative evaluation of polycyclic systems containing five-, six- and seven-membered rings. 2. Scope and limitations. , 1984 .

[22]  A. Ghose,et al.  Geometrically feasible binding modes of a flexible ligand molecule at the receptor site , 1985 .

[23]  A Wlodawer,et al.  Nuclear magnetic resonance and neutron diffraction studies of the complex of ribonuclease A with uridine vanadate, a transition-state analogue. , 1985, Biochemistry.

[24]  A. Fersht,et al.  Use of binding energy in catalysis analyzed by mutagenesis of the tyrosyl-tRNA synthetase. , 1986, Biochemistry.

[25]  Udo Heinemann,et al.  Restrained least‐squares refinement of the crystal structure of the ribonuclease T1*2'‐guanylic acid complex at 1·9 Å resolution , 1987 .

[26]  J L Sussman,et al.  Refined crystal structure of dogfish M4 apo-lactate dehydrogenase. , 1989, Journal of molecular biology.

[27]  Robert P. Sheridan,et al.  Docking Flexible Ligands to Macromolecular Receptors by Molecular Shape. , 1987 .

[28]  Timothy F. Havel,et al.  A new approach to the problem of docking two molecules: The ellipsoid algorithm , 1987, Biopolymers.

[29]  P. Brick,et al.  Crystal structure of a deletion mutant of a tyrosyl-tRNA synthetase complexed with tyrosine. , 1987, Journal of molecular biology.

[30]  Brian W. Matthews,et al.  Structural basis of the action of thermolysin and related zinc peptidases , 1988 .

[31]  Gordon M. Crippen,et al.  Distance Geometry and Molecular Conformation , 1988 .

[32]  Larry S. Davis,et al.  Pose Determination of a Three-Dimensional Object Using Triangle Pairs , 1988, IEEE Trans. Pattern Anal. Mach. Intell..

[33]  I. Kuntz,et al.  Using shape complementarity as an initial screen in designing ligands for a receptor binding site of known three-dimensional structure. , 1988, Journal of medicinal chemistry.

[34]  J. Wendoloski,et al.  Structural origins of high-affinity biotin binding to streptavidin. , 1989, Science.

[35]  J. Drenth,et al.  The coenzyme analogue adenosine 5-diphosphoribose displaces FAD in the active site of p-hydroxybenzoate hydroxylase. An x-ray crystallographic investigation. , 1989, Biochemistry.

[36]  H. Zollner,et al.  Handbook of Enzyme Inhibitors , 1989 .

[37]  The coenzyme analogue adenosine 5-diphosphoribose displaces FAD in the active site of p-hydroxybenzoate hydroxylase. An x-ray crystallographic investigation. , 1989 .

[38]  D. Goodsell,et al.  Automated docking of substrates to proteins by simulated annealing , 1990, Proteins.

[39]  Valerie J. Gillet,et al.  Automated structure design in 3D , 1990 .

[40]  D. Turk,et al.  Geometry of binding of the benzamidine- and arginine-based inhibitors Nα-(2-naphthyl-sulphonyl-glycyl)-dl-p-amidinophenylalanyl-piperidine (NAPAP) and (2R,4R)-4-methyl-1-[Nα-(3-methyl-1,2,3,4-tetrahydro-8-quinolinesulphonyl)-l-arginyl]-2-piperidine carboxylic acid (MQPA) to human α-thrombin , 1990 .

[41]  W Bode,et al.  Geometry of binding of the benzamidine- and arginine-based inhibitors N alpha-(2-naphthyl-sulphonyl-glycyl)-DL-p-amidinophenylalanyl-pipe ridine (NAPAP) and (2R,4R)-4-methyl-1-[N alpha-(3-methyl-1,2,3,4-tetrahydro-8- quinolinesulphonyl)-L-arginyl]-2-piperidine carboxylic acid (MQPA) to human alpha-t , 1990, European journal of biochemistry.

[42]  M. Katharine Holloway,et al.  X-Ray Crystal Structure of the HIV Protease Complex with L-700,417, an Inhibitor with Pseudo C2 Symmetry , 1991 .

[43]  D. Banner,et al.  Crystallographic analysis at 3.0-A resolution of the binding to human thrombin of four active site-directed inhibitors. , 1994, The Journal of biological chemistry.

[44]  S E Ealick,et al.  Application of crystallographic and modeling methods in the design of purine nucleoside phosphorylase inhibitors. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[45]  W. Howe,et al.  Computer design of bioactive molecules: A method for receptor‐based de novo ligand design , 1991, Proteins.

[46]  G Vriend,et al.  Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex. , 1991, Journal of molecular biology.

[47]  C L Verlinde,et al.  Anion binding at the active site of trypanosomal triosephosphate isomerase. Monohydrogen phosphate does not mimic sulphate. , 1991, European journal of biochemistry.

[48]  Akiko Itai,et al.  Automatic creation of drug candidate structures based on receptor structure. Starting point for artificial lead generation , 1991 .

[49]  Andrew Smellie,et al.  Fast drug-receptor mapping by site-directed distances: a novel method of predicting new pharmacological leads , 1991, J. Chem. Inf. Comput. Sci..

[50]  Hans-Joachim Böhm,et al.  The computer program LUDI: A new method for the de novo design of enzyme inhibitors , 1992, J. Comput. Aided Mol. Des..

[51]  Hans-Joachim Böhm,et al.  LUDI: rule-based automatic design of new substituents for enzyme inhibitor leads , 1992, J. Comput. Aided Mol. Des..

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

[53]  J. J. Wendoloski,et al.  Crystallographic and thermodynamic comparison of natural and synthetic ligands bound to streptavidin , 1992 .

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

[55]  Mark A. Murcko,et al.  GenStar: A method for de novo drug design , 1993, J. Comput. Aided Mol. Des..

[56]  J. Scott Dixon,et al.  A good ligand is hard to find: Automated docking methods , 1993 .

[57]  S Mangani,et al.  High-resolution structure of the complex between carboxypeptidase A and L-phenyl lactate. , 1993, Acta crystallographica. Section D, Biological crystallography.

[58]  THE HIGH RESOLUTION CRYSTAL STRUCTURE OF THE COMPLEX BETWEEN CARBOXYPEPTIDASE A AND L-PHENYL LACTATE , 1994 .

[59]  K Gubernator,et al.  Design and synthesis of potent and highly selective thrombin inhibitors. , 1994, Journal of medicinal chemistry.

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

[61]  H J Berendsen,et al.  Molecular dynamics simulation of the docking of substrates to proteins , 1994, Proteins.

[62]  G. Klebe.  The use of composite crystal-field environments in molecular recognition and the de novo design of protein ligands. , 1994, Journal of molecular biology.

[63]  Valerie J. Gillet,et al.  SPROUT: Recent developments in the de novo design of molecules , 1994, J. Chem. Inf. Comput. Sci..

[64]  Clark F. Olson.  Time and space efficient pose clustering , 1994, 1994 Proceedings of IEEE Conference on Computer Vision and Pattern Recognition.

[65]  P. Colman,et al.  Structure-based drug design. , 1994, Current opinion in structural biology.

[66]  I. Kuntz,et al.  Structure-Based Molecular Design , 1994 .

[67]  Gerhard Klebe,et al.  A fast and efficient method to generate biologically relevant conformations , 1994, J. Comput. Aided Mol. Des..

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

[69]  Hans-Joachim Böhm,et al.  The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[70]  Struktur Und Guanosintriphosphat-Hydrolysemechanismus Des C-Terminal Verkuerzten Menschlichen Krebsproteins P21-H-Ras , 1994 .

[71]  Andrew R. Leach,et al.  Current methods for site-directed structure generation , 1994, J. Comput. Aided Mol. Des..

[72]  R. Glen,et al.  Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. , 1995, Journal of molecular biology.

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

[74]  Thomas Lengauer,et al.  Time-Efficient Docking of Flexible Ligands into Active Sites of Proteins , 1995, ISMB.

[75]  Richard S. Judson,et al.  Docking flexible molecules: A case study of three proteins , 1995, J. Comput. Chem..

[76]  R. Nussinov,et al.  A geometry-based suite of molecular docking processes. , 1995, Journal of molecular biology.

[77]  J. Scott Dixon,et al.  Flexible ligand docking using a genetic algorithm , 1995, J. Comput. Aided Mol. Des..

[78]  Ruth Nussinov,et al.  An automated computer vision and robotics-based technique for 3-D flexible biomolecular docking and matching , 1995, Comput. Appl. Biosci..

[79]  C. Frömmel,et al.  The automatic search for ligand binding sites in proteins of known three-dimensional structure using only geometric criteria. , 1996, Journal of molecular biology.

[80]  Thomas Lengauer,et al.  Placement of medium-sized molecular fragments into active sites of proteins , 1996, J. Comput. Aided Mol. Des..