Scoring functions for docking

In the fields of computational chemistry and molecular modelling, scoring functions are approximate mathematical methods used to predict the strength of the non-covalent interaction (also referred to as binding affinity) between two molecules. Most commonly one of the molecules is a small organic compound such as a drug and the second is the drug's biological target such as a protein receptor.^[1] Scoring functions have also been developed to predict the strength of other types of intermolecular interactions, for example between two proteins^[2] or between protein and DNA.^[3]

Additional recommended knowledge

Guide to balance cleaning: 8 simple steps

Better weighing performance in 6 easy steps

Correct Test Weight Handling Guide: 12 Practical Tips

Utility

Scoring functions are useful for:^[4]

• virtual screening of small molecule databases to identify novel small molecules that bind to a protein target of interest and therefore are useful starting points for drug discovery^[5]
• in lead optimization of screening hits to optimize their affinity and selectivity^[6]

Prerequisites

Scoring functions are normally parameterized (or trained) against a data set consisting of experimentally determined binding affinities between molecular species similar to the species that one wishes to predict.

For predictions of affinities of ligands for proteins, the following must be known or predicted:

• tertiary structure (arrangement of atoms in three dimensional space) of the protein. Protein structures may be determined by experimental techniques such as X-ray crystallography or solution phase NMR methods or predicted by homology modelling.
• the three dimensional shape (referred to as the "active" conformation) of the other binding partner
• the relative orientation of the two binding partners (often referred to as the "binding-mode") in the complex

The information above yields the three dimensional structure of the complex. Based on this three dimensional structure, the scoring function can then estimate the strength of the association between the two partners in the complex using one of the methods outlined below. Finally the scoring function itself may be used to help predict both the binding mode and the active conformation of the small molecule binding partner in the complex.

Classes

There are three general classes of scoring functions:

• Force-field - affinities are estimated by summing the strength of intermolecular van der Waals and electrostatic interactions between all atoms of the two molecules in the complex. The intramolecular energies (also referred to as strain energy) of the two binding partners are also frequently included.
• Empirical - based on counting the number off various types of interactions between the two binding partners.^[7] These interactions may include for example
  • number of hydrogen bonds (favorable electrostatic contribution to affinity),
  • hydrophobic — hydrophobic contacts (favorable),
  • hydrophilic — hydrophobic contacts (unfavorable),
  • number of rotatable bonds immobilized in complex formation (unfavorable entropic contribution).
• Knowledge - based on statistical observations of intermolecular close contacts in large 3D databases (such as the Cambridge Structural Database or Protein Data Bank) which are used to derive "potentials of mean force". This method is founded on the assumption that close intermolecular interactions between certain types of atoms or functional groups that occur more frequently than one would expect by a random distribution are likely to be energetically favorable and therefore contribute favorably to binding affinity.^[8]

Finally hybrid scoring functions have also been developed in which the terms from two or more of the above types of scoring functions are combined into one function.

References