Home LiteratureArticle Details
PMID: 9385547 Published · ppublish English Journal Article

Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes.

Journal of computer-aided molecular design ·Vol. 11 ·No. 5 ·1997-09-00 ·Pages 425-45

Eldridge MD, Murray CW, Auton TR, Paolini GV, Mee RP

Abstract

This paper describes the development of a simple empirical scoring function designed to estimate the free energy of binding for a protein-ligand complex when the 3D structure of the complex is known or can be approximated. The function uses simple contact terms to estimate lipophilic and metal-ligand binding contributions, a simple explicit form for hydrogen bonds and a term which penalises flexibility. The coefficients of each term are obtained using a regression based on 82 ligand-receptor complexes for which the binding affinity is known. The function reproduces the binding affinity of the complexes with a cross-validated error of 8.68 kJ/mol. Tests on internal consistency indicate that the coefficients obtained are stable to changes in the composition of the training set. The function is also tested on two test sets containing a further 20 and 10 complexes, respectively. The deficiencies of this type of function are discussed and it is compared to approaches by other workers.

MeSH Terms
Binding Sites Drug Design Endopeptidases/chemistry Ligands Protein Binding Proteins/chemistry Regression Analysis Software Thermodynamics
Chemicals
Ligands Proteins Endopeptidases
Authors & Affiliations
5 authors, click to expand affiliations / ORCID
Eldridge M D
Proteus Molecular Design Ltd., Macclesfield, Cheshire, U.K.
Murray C W
Auton T R
Paolini G V
Mee R P
References (75)
75 references, click to expand
  1. Hydroxyethylene isostere inhibitors of human immunodeficiency virus-1 protease: structure-activity analysis using enzyme kinetics, X-ray crystallography, and infected T-cell assays.
    Biochemistry. 1992 Jul 28;31(29):6646-59 PMID: 1637805
  2. Structure of human neutrophil collagenase reveals large S1' specificity pocket.
    Nat Struct Biol. 1994 Feb;1(2):119-23 PMID: 7656015
  3. Crystallographic analysis of a complex between human immunodeficiency virus type 1 protease and acetyl-pepstatin at 2.0-A resolution.
    J Biol Chem. 1990 Aug 25;265(24):14209-19 PMID: 2201682
  4. Molecular basis of crossreactivity and the limits of antibody-antigen complementarity.
    Nature. 1993 Oct 28;365(6449):859-63 PMID: 8413674
  5. Crystal structures of native and inhibited forms of human cathepsin D: implications for lysosomal targeting and drug design.
    Proc Natl Acad Sci U S A. 1993 Jul 15;90(14):6796-800 PMID: 8393577
  6. Hammerhead: fast, fully automated docking of flexible ligands to protein binding sites.
    Chem Biol. 1996 Jun;3(6):449-62 PMID: 8807875
  7. Local and transmitted conformational changes on complexation of an anti-sweetener Fab.
    J Mol Biol. 1994 Feb 11;236(1):247-74 PMID: 7893280
  8. A new method for predicting binding affinity in computer-aided drug design.
    Protein Eng. 1994 Mar;7(3):385-91 PMID: 8177887
  9. Analyses of ligand binding in five endothiapepsin crystal complexes and their use in the design and evaluation of novel renin inhibitors.
    J Med Chem. 1993 Nov 26;36(24):3809-20 PMID: 8254610
  10. Effect of conformational flexibility and solvation on receptor-ligand binding free energies.
    Biochemistry. 1994 Nov 29;33(47):13977-88 PMID: 7947806
  11. Recognition of a cell-surface oligosaccharide of pathogenic Salmonella by an antibody Fab fragment.
    Science. 1991 Jul 26;253(5018):442-5 PMID: 1713710
  12. Refined structure for the complex of 1-deoxynojirimycin with glucoamylase from Aspergillus awamori var. X100 to 2.4-A resolution.
    Biochemistry. 1993 Feb 16;32(6):1618-26 PMID: 8431441
  13. X-ray crystallographic studies of a series of penicillin-derived asymmetric inhibitors of HIV-1 protease.
    Biochemistry. 1994 Jul 19;33(28):8417-27 PMID: 8031777
  14. Correlation of binding affinities with non-bonded interaction energies of thrombin-inhibitor complexes.
    Acta Crystallogr D Biol Crystallogr. 1995 Jul 1;51(Pt 4):560-6 PMID: 15299844
  15. The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure.
    J Comput Aided Mol Des. 1994 Jun;8(3):243-56 PMID: 7964925
  16. Design, activity, and 2.8 A crystal structure of a C2 symmetric inhibitor complexed to HIV-1 protease.
    Science. 1990 Aug 3;249(4968):527-33 PMID: 2200122
  17. Structure of complex of synthetic HIV-1 protease with a substrate-based inhibitor at 2.3 A resolution.
    Science. 1989 Dec 1;246(4934):1149-52 PMID: 2686029
  18. X-ray crystallographic analysis of inhibition of endothiapepsin by cyclohexyl renin inhibitors.
    Biochemistry. 1992 Sep 8;31(35):8142-50 PMID: 1525155
  19. Atomic structure of adenosine deaminase complexed with a transition-state analog: understanding catalysis and immunodeficiency mutations.
    Science. 1991 May 31;252(5010):1278-84 PMID: 1925539
  20. Binding of N-carboxymethyl dipeptide inhibitors to thermolysin determined by X-ray crystallography: a novel class of transition-state analogues for zinc peptidases.
    Biochemistry. 1984 Nov 20;23(24):5724-9 PMID: 6395881
  21. Principles and pitfalls in designing site-directed peptide ligands.
    Proteins. 1993 Jul;16(3):246-67 PMID: 8346191
  22. Calculation of relative binding affinities of purine nucleoside phosphorylase inhibitors.
    Acta Crystallogr D Biol Crystallogr. 1995 Jul 1;51(Pt 4):536-40 PMID: 15299841
  23. 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.
    J Mol Biol. 1992 Aug 20;226(4):1085-99 PMID: 1518046
  24. X-ray analyses of aspartic proteinases. III Three-dimensional structure of endothiapepsin complexed with a transition-state isostere inhibitor of renin at 1.6 A resolution.
    J Mol Biol. 1990 Dec 20;216(4):1017-29 PMID: 2266553
  25. High resolution X-ray analyses of renin inhibitor-aspartic proteinase complexes.
    Nature. 1987 May 28-Jun 3;327(6120):349-52 PMID: 3295561
  26. Electron paramagnetic resonance detectable states of cytochrome P-450cam.
    Biochemistry. 1980 Jul 22;19(15):3590-9 PMID: 6250573
  27. Crystal structure of the complex between carboxypeptidase A and the biproduct analog inhibitor L-benzylsuccinate at 2.0 A resolution.
    J Mol Biol. 1992 Jan 20;223(2):573-8 PMID: 1738164
  28. X-ray crystallographic structure of a complex between a synthetic protease of human immunodeficiency virus 1 and a substrate-based hydroxyethylamine inhibitor.
    Proc Natl Acad Sci U S A. 1990 Nov;87(22):8805-9 PMID: 2247451
  29. Influenza B virus neuraminidase can synthesize its own inhibitor.
    Structure. 1993 Sep 15;1(1):19-26 PMID: 8069621
  30. Crystallographic refinement of human serum retinol binding protein at 2A resolution.
    Proteins. 1990;8(1):44-61 PMID: 2217163
  31. Rational design of potent, bioavailable, nonpeptide cyclic ureas as HIV protease inhibitors.
    Science. 1994 Jan 21;263(5145):380-4 PMID: 8278812
  32. Sugar-binding and crystallographic studies of an arabinose-binding protein mutant (Met108Leu) that exhibits enhanced affinity and altered specificity.
    Biochemistry. 1991 Jul 16;30(28):6861-6 PMID: 2069949
  33. A Pro to Gly mutation in the hinge of the arabinose-binding protein enhances binding and alters specificity. Sugar-binding and crystallographic studies.
    J Biol Chem. 1990 Sep 25;265(27):16592-603 PMID: 2204627
  34. The Protein Data Bank: a computer-based archival file for macromolecular structures.
    J Mol Biol. 1977 May 25;112(3):535-42 PMID: 875032
  35. 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-thrombin. X-ray crystallographic determination of the NAPAP-trypsin complex and modeling of NAPAP-thrombin and MQPA-thrombin.
    Eur J Biochem. 1990 Oct 5;193(1):175-82 PMID: 2226434
  36. Solvent accessibility as a predictive tool for the free energy of inhibitor binding to the HIV-1 protease.
    Protein Sci. 1995 Jul;4(7):1356-64 PMID: 7670378
  37. Hydrogen bonding and biological specificity analysed by protein engineering.
    Nature. 1985 Mar 21-27;314(6008):235-8 PMID: 3845322
  38. Three-dimensional structure of the tetragonal crystal form of egg-white avidin in its functional complex with biotin at 2.7 A resolution.
    J Mol Biol. 1993 Jun 5;231(3):698-710 PMID: 8515446
  39. Crystal structure of the complex of carboxypeptidase A with a strongly bound phosphonate in a new crystalline form: comparison with structures of other complexes.
    Biochemistry. 1990 Jun 12;29(23):5546-55 PMID: 2386784
  40. Scoring noncovalent protein-ligand interactions: a continuous differentiable function tuned to compute binding affinities.
    J Comput Aided Mol Des. 1996 Oct;10(5):427-40 PMID: 8951652
  41. Relocating a negative charge in the binding pocket of trypsin.
    J Mol Biol. 1993 Apr 5;230(3):934-49 PMID: 8478942
  42. Crystal structure and ligand-binding studies of a screened peptide complexed with streptavidin.
    Biochemistry. 1992 Oct 6;31(39):9350-4 PMID: 1390720
  43. A preference-based free-energy parameterization of enzyme-inhibitor binding. Applications to HIV-1-protease inhibitor design.
    Protein Sci. 1995 Sep;4(9):1881-903 PMID: 8528086
  44. Comparison of the structures of three carboxypeptidase A-phosphonate complexes determined by X-ray crystallography.
    Biochemistry. 1991 Aug 20;30(33):8171-80 PMID: 1868092
  45. CNAD: a potent and specific inhibitor of alcohol dehydrogenase.
    J Med Chem. 1994 Feb 4;37(3):392-9 PMID: 8308865
  46. Allosteric modifiers of hemoglobin. 2. Crystallographically determined binding sites and hydrophobic binding/interaction analysis of novel hemoglobin oxygen effectors.
    J Med Chem. 1991 Feb;34(2):758-67 PMID: 1995898
  47. High-resolution X-ray diffraction study of the complex between endothiapepsin and an oligopeptide inhibitor: the analysis of the inhibitor binding and description of the rigid body shift in the enzyme.
    EMBO J. 1989 Aug;8(8):2179-88 PMID: 2676515
  48. Crystallographic study of coenzyme, coenzyme analogue and substrate binding in 6-phosphogluconate dehydrogenase: implications for NADP specificity and the enzyme mechanism.
    Structure. 1994 Jul 15;2(7):651-68 PMID: 7922042
  49. Prediction of drug binding affinities by comparative binding energy analysis.
    J Med Chem. 1995 Jul 7;38(14 ):2681-91 PMID: 7629807
  50. A check on rational drug design: crystal structure of a complex of human immunodeficiency virus type 1 protease with a novel gamma-turn mimetic inhibitor.
    J Med Chem. 1995 Aug 18;38(17):3246-52 PMID: 7650677
  51. A strategy for theoretical binding constant, Ki, calculations for neuraminidase aromatic inhibitors designed on the basis of the active site structure of influenza virus neuraminidase.
    Proteins. 1995 Oct;23(2):264-77 PMID: 8592707
  52. Interactions of nucleotides with fully unadenylylated glutamine synthetase from Salmonella typhimurium.
    Biochemistry. 1994 Sep 20;33(37):11184-8 PMID: 7727369
  53. Definition and display of steric, hydrophobic, and hydrogen-bonding properties of ligand binding sites in proteins using Lee and Richards accessible surface: validation of a high-resolution graphical tool for drug design.
    J Med Chem. 1992 May 15;35(10):1671-84 PMID: 1588550
  54. A fast flexible docking method using an incremental construction algorithm.
    J Mol Biol. 1996 Aug 23;261(3):470-89 PMID: 8780787
  55. X-ray crystallographic investigation of substrate binding to carboxypeptidase A at subzero temperature.
    Proc Natl Acad Sci U S A. 1986 Oct;83(20):7568-72 PMID: 3463986
  56. Binding of hydroxamic acid inhibitors to crystalline thermolysin suggests a pentacoordinate zinc intermediate in catalysis.
    Biochemistry. 1981 Nov 24;20(24):6912-20 PMID: 7317361
  57. Crystallographic structural analysis of phosphoramidates as inhibitors and transition-state analogs of thermolysin.
    Eur J Biochem. 1986 Jun 2;157(2):261-8 PMID: 3709536
  58. Structure-based understanding of ligand affinity using human thrombin as a model system.
    Biochemistry. 1996 Jul 30;35(30):9690-9 PMID: 8703940
  59. Affinity and specificity of serine endopeptidase-protein inhibitor interactions. Empirical free energy calculations based on X-ray crystallographic structures.
    J Mol Biol. 1993 Dec 5;234(3):661-79 PMID: 8254666
  60. A priori prediction of activity for HIV-1 protease inhibitors employing energy minimization in the active site.
    J Med Chem. 1995 Jan 20;38(2):305-17 PMID: 7830273
  61. Three-dimensional structure of Escherichia coli dihydrodipicolinate reductase.
    Biochemistry. 1995 Mar 21;34(11):3502-12 PMID: 7893645
  62. 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.
    Protein Eng. 1995 Jul;8(7):677-91 PMID: 8577696
  63. Direct observation by X-ray analysis of the tetrahedral "intermediate" of aspartic proteinases.
    Protein Sci. 1992 Mar;1(3):322-8 PMID: 1304340
  64. Chelation of serine 39 to Mg2+ latches a gate at the active site of enolase: structure of the bis(Mg2+) complex of yeast enolase and the intermediate analog phosphonoacetohydroxamate at 2.1-A resolution.
    Biochemistry. 1994 Aug 9;33(31):9333-42 PMID: 8049235
  65. Geometry of binding of the N alpha-tosylated piperidides of m-amidino-, p-amidino- and p-guanidino phenylalanine to thrombin and trypsin. X-ray crystal structures of their trypsin complexes and modeling of their thrombin complexes.
    FEBS Lett. 1991 Aug 5;287(1-2):133-8 PMID: 1879520
  66. PRO_SELECT: combining structure-based drug design and combinatorial chemistry for rapid lead discovery. 1. Technology.
    J Comput Aided Mol Des. 1997 Mar;11(2):193-207 PMID: 9089436
  67. Computational methods to predict binding free energy in ligand-receptor complexes.
    J Med Chem. 1995 Dec 22;38(26):4953-67 PMID: 8544170
  68. Dehydrogenase binding by tiazofurin anabolites.
    J Med Chem. 1990 Apr;33(4):1123-7 PMID: 1969483
  69. Conformational changes of small molecules binding to proteins.
    Bioorg Med Chem. 1995 Apr;3(4):411-28 PMID: 8581425
  70. Crystallographic analysis of transition-state mimics bound to penicillopepsin: phosphorus-containing peptide analogues.
    Biochemistry. 1992 Jun 9;31(22):5201-14 PMID: 1606144
  71. Calculation of the free energy of association for protein complexes.
    Protein Sci. 1992 Jan;1(1):169-81 PMID: 1339024
  72. Prediction of new serine proteinase inhibitors.
    Nat Struct Biol. 1994 Oct;1(10):735-43 PMID: 7634078
  73. Major enhancement of the affinity of an enzyme for a transition-state analog by a single hydroxyl group.
    Science. 1989 Mar 24;243(4898):1591-3 PMID: 2928795
  74. Refined 1.89-A structure of the histidine-binding protein complexed with histidine and its relationship with many other active transport/chemosensory proteins.
    Biochemistry. 1994 Apr 26;33(16):4769-79 PMID: 8161536
  75. A relationship between protein stability and protein function.
    Proc Natl Acad Sci U S A. 1995 Jan 17;92(2):452-6 PMID: 7831309
Article Info
Journal
Journal of computer-aided molecular design
Abbr.
J Comput Aided Mol Des
ISSN
0920-654X
Published
1997-09-00
Pages
425-45
Language
English
Region
Netherlands
NLM ID
8710425
Subset
IM
Analysis Services
Analysis Services

Contact

No. 2 Wenbo Road, Zhangqiu District, Jinan, Shandong

Qilu Normal University · Genelibs Bioinformatics Lab

750 Shunhua Rd, Jinan

2F, Bldg F, University Science Park

Tel: 0531-88819269

WeChat Official Account

Follow our WeChat subscription account for real-time updates and the latest in medical and biological research.


Business Email

E-mail: product@genelibs.com