An example of a protein ligand found by database mining: description of the docking method and its verification by a 2.3 A X-ray structure of a thrombin-ligand complex.

A computer program (SANDOCK) has been developed for the automated docking of small ligands to a target protein. It uses a guided matching algorithm to fit ligand atoms into the protein binding pocket. The protein is described by a modified Lee-Richard's dotted surface with each dot coded by chemical property and accessibility. Orientations of the ligand in the active site are generated such that a chemical and a shape complementary between the ligand and the active site cavity have to be fulfilled. The generated fits are evaluated with scoring functions which account for van der Waals, hydrophobic and hydrogen bonding interactions. This newly developed docking program can efficiently screen very large databases in a reasonable time and has been used to successfully identify novel ligands. The X-ray structure of a thrombin-ligand complex predicted by SANDOCK is described. The ligand binds to thrombin with a Kd of 65 microM and has an rmsd of 0.7 A for all ligand atoms from the predicted binding mode by SANDOCK.

[1]  Y. Martin,et al.  3D database searching in drug design. , 1992, Journal of medicinal chemistry.

[2]  S. Evans,et al.  p-Aminobenzamidine as a fluorescent probe for the active site of serine proteases. , 1982, The Journal of biological chemistry.

[3]  D S Goodsell,et al.  Automated docking of flexible ligands: Applications of autodock , 1996, Journal of molecular recognition : JMR.

[4]  C. Ehrhardt,et al.  Synthetic low-molecular weight thrombin inhibitors: molecular design and pharmacological profile. , 1993, Trends in pharmacological sciences.

[5]  M. L. Connolly.  Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

[6]  T. Stouch,et al.  Affinity and specificity of serine endopeptidase-protein inhibitor interactions. Empirical free energy calculations based on X-ray crystallographic structures. , 1993, Journal of molecular biology.

[7]  J. Dunitz.  The entropic cost of bound water in crystals and biomolecules. , 1994, Science.

[8]  R C Wade,et al.  Further development of hydrogen bond functions for use in determining energetically favorable binding sites on molecules of known structure. 2. Ligand probe groups with the ability to form more than two hydrogen bonds. , 1993, Journal of medicinal chemistry.

[9]  A. Berger,et al.  On the size of the active site in proteases. I. Papain. , 1967, Biochemical and biophysical research communications.

[10]  Hans-Joachim Böhm,et al.  Towards the automatic design of synthetically accessible protein ligands: Peptides, amides and peptidomimetics , 1996, J. Comput. Aided Mol. Des..

[11]  Gerhard Klebe,et al.  What Can We Learn from Molecular Recognition in Protein–Ligand Complexes for the Design of New Drugs? , 1996 .

[12]  W. Bode,et al.  Crystal structure of the thrombin‐hirudin complex: a novel mode of serine protease inhibition. , 1990, The EMBO journal.

[13]  PatrickY.-S. Lam,et al.  Rational design of potent, bioavailable, nonpeptide cyclic ureas as HIV protease inhibitors. , 1994, Science.

[14]  A. McPherson,et al.  Current approaches to macromolecular crystallization. , 1990, European journal of biochemistry.

[15]  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.

[16]  Ajay,et al.  Computational methods to predict binding free energy in ligand-receptor complexes. , 1995, Journal of medicinal chemistry.

[17]  Hans-Joachim Böhm,et al.  Site-directed structure generation by fragment-joining , 1995 .

[18]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[19]  K. Mann,et al.  Human d-Phe-Pro-Arg-CH2-α-thrombin crystallization and diffraction data , 1989 .

[20]  Mary L. Westbrook,et al.  Structure of the hirugen and hirulog 1 complexes of alpha-thrombin. , 1991 .

[21]  Daniel A. Gschwend,et al.  Molecular docking towards drug discovery , 1996, Journal of molecular recognition : JMR.

[22]  H Brandstetter,et al.  Enzyme flexibility, solvent and 'weak' interactions characterize thrombin-ligand interactions: implications for drug design. , 1996, Structure.

[23]  Irwin D. Kuntz,et al.  A fast and efficient method for 2D and 3D molecular shape description , 1992, J. Comput. Aided Mol. Des..

[24]  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..

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

[26]  Randy J. Read,et al.  Monte Carlo algorithms for docking to proteins , 1995 .

[27]  H Brandstetter,et al.  Refined 2.3 A X-ray crystal structure of bovine thrombin complexes formed with the benzamidine and arginine-based thrombin inhibitors NAPAP, 4-TAPAP and MQPA. A starting point for improving antithrombotics. , 1992, Journal of molecular biology.

[28]  P. Grootenhuis,et al.  Correlation of binding affinities with non-bonded interaction energies of thrombin-inhibitor complexes. , 1995, Acta crystallographica. Section D, Biological crystallography.

[29]  G Klebe,et al.  Energetic and entropic factors determining binding affinity in protein-ligand complexes. , 1997, Journal of receptor and signal transduction research.

[30]  Brian K. Shoichet,et al.  Molecular docking using shape descriptors , 1992 .

[31]  R F Standaert,et al.  Atomic structures of the human immunophilin FKBP-12 complexes with FK506 and rapamycin. , 1993, Journal of molecular biology.

[32]  J. Zou,et al.  Improved methods for building protein models in electron density maps and the location of errors in these models. , 1991, Acta crystallographica. Section A, Foundations of crystallography.

[33]  Randy J. Read,et al.  A multiple‐start Monte Carlo docking method , 1992 .

[34]  Robert Huber,et al.  The refined 1.9 A crystal structure of human alpha‐thrombin: interaction with D‐Phe‐Pro‐Arg chloromethylketone and significance of the Tyr‐Pro‐Pro‐Trp insertion segment. , 1989 .

[35]  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..

[36]  M. Walkinshaw,et al.  Hydrophilicity of proteins and DNA , 1990 .

[37]  M J Sternberg,et al.  Empirical scale of side-chain conformational entropy in protein folding. , 1993, Journal of molecular biology.

[38]  S. Anderson,et al.  Water molecules participate in proteinase‐inhibitor interactions: Crystal structures of Leu18, Ala18, and Gly18 variants of turkey ovomucoid inhibitor third domain complexed with Streptomyces griseus proteinase B , 1995, Protein science : a publication of the Protein Society.

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

[40]  V. Mikol,et al.  The role of water molecules in the structure-based design of (5-hydroxynorvaline)-2-cyclosporin: synthesis, biological activity, and crystallographic analysis with cyclophilin A. , 1995, Journal of medicinal chemistry.

[41]  I. Kuntz,et al.  Automated docking with grid‐based energy evaluation , 1992 .

[42]  Axel T. Brunger,et al.  X-PLOR Version 3.1: A System for X-ray Crystallography and NMR , 1992 .

[43]  B. Matthews,et al.  Specificity of ligand binding in a buried nonpolar cavity of T4 lysozyme: linkage of dynamics and structural plasticity. , 1995, Biochemistry.

[44]  Gennady M Verkhivker,et al.  Empirical free energy calculations of ligand-protein crystallographic complexes. I. Knowledge-based ligand-protein interaction potentials applied to the prediction of human immunodeficiency virus 1 protease binding affinity. , 1995, Protein engineering.

[45]  M. Karplus,et al.  Molecular recognition in proteins. Simulation analysis of substrate binding by a tyrosyl-tRNA synthetase mutant. , 1994, Journal of molecular biology.

[46]  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.

[47]  M Karplus,et al.  HOOK: A program for finding novel molecular architectures that satisfy the chemical and steric requirements of a macromolecule binding site , 1994, Proteins.

[48]  B. Matthews,et al.  Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4 lysozyme. , 1995, Biochemistry.

[49]  K D Wilkinson,et al.  The neuron-specific protein PGP 9.5 is a ubiquitin carboxyl-terminal hydrolase. , 1989, Science.

[50]  Andrew Smellie,et al.  Poling: Promoting conformational variation , 1995, J. Comput. Chem..

[51]  J. Pflugrath,et al.  Crystal orientation and X-ray pattern prediction routines for area-detector diffractometer systems in macromolecular crystallography , 1987 .

[52]  W. Bode,et al.  The clot thickens: clues provided by thrombin structure. , 1995, Trends in biochemical sciences.

[53]  F. Quiocho,et al.  Crystallographic evidence of a large ligand-induced hinge-twist motion between the two domains of the maltodextrin binding protein involved in active transport and chemotaxis. , 1992, Biochemistry.

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

[55]  R. Wade,et al.  New hydrogen-bond potentials for use in determining energetically favorable binding sites on molecules of known structure. , 1989, Journal of medicinal chemistry.

[56]  M J Sternberg,et al.  Side‐chain conformational entropy in protein folding , 1995, Protein science : a publication of the Protein Society.

[57]  I. Kuntz,et al.  Structure-based discovery of inhibitors of thymidylate synthase. , 1993, Science.

[58]  G. Claeson,et al.  In vitro and in vivo characterization of a neutral boron-containing thrombin inhibitor. , 1993, The Journal of biological chemistry.

[59]  Diana C. Roe,et al.  BUILDER v.2: Improving the chemistry of a de novo design strategy , 1995, J. Comput. Aided Mol. Des..