Skip to main content


We're creating a new version of this page. See preview

  • Poster presentation
  • Open Access

Rescoring of docking poses using force field-based methods

  • 1,
  • 1 and
  • 1
Journal of Cheminformatics20102 (Suppl 1) :P52

  • Published:


  • Fast Filter
  • Born Surface Area
  • Complementary Shape
  • Generalize Born Surface
  • Generalize Born Surface Area

Existing protein-ligand docking methods computationally screen thousands to millions of organic molecules against protein structures, trying to find those with complementary shapes and highest binding free energies. To allow large molecular databases to be screened rapidly, simple and approximative scoring functions are used as a fast filter, resulting in low hit rates. Therefore, docking hit lists are commonly rescored in a final step by more rigorous and time-consuming methods to gain a more accurate final list of ranked compounds.

Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) or generalized Born surface area (MM-GBSA) methods are currently considered to be suitable techniques for rescoring. These physically realistic approaches incorporate more sophisticated models for solvation and electrostatic interactions than most scoring functions. Hence, they can discriminate more reliably between correct and incorrect docking poses.

A high-throughput rescoring protocol using force field-based methods has been proposed by Brown and Muchmore [1]. They used 18 in-house urokinase-ligand crystal structures and their corresponding experimentally determined binding free energies as a test set. On the basis of this rescoring protocol we tested several molecular dynamics simulation protocols in combination with different MM-PBSA, and MM-GBSA calculation procedures using Amber 10 [2] and Gromacs 4 [3]. Considering performance and accuracy, our best rescoring protocol performs similarly to the one described by Brown and Muchmore [1]. It has a comparable run-time and achieves a correlation between experimental and rescored values of 0.88 compared to 0.87.

However, we used urokinase-ligand complexes generated using the docking program Glide [4] instead of crystal structures. Additionally, this shows that our rescoring protocol improves the correlation of 0.57 between experimental values and Glide scores significantly to 0.88, thereby achieving a more accurate list of ranked compounds.

The protocol will be incorporated into BALLView [5], an open-source molecular viewer and modeling tool. Thus, it will be available free of charge and can be conveniently used to rescore (docked) protein-ligand complexes.

Authors’ Affiliations

Center for Bioinformatics Tübingen, Eberhard Karls University, Sand 14, 72076 Tübingen, Germany


  1. Brown SP: J Chem Inf Model. 2007, 47: 1493-10.1021/ci700041j.View ArticleGoogle Scholar
  2. Case DA: 2008, University of California, San FranciscoGoogle Scholar
  3. Hess B: J Chem Theory Comput. 2008, 4: 435-10.1021/ct700301q.View ArticleGoogle Scholar
  4. Friesner RA: J Med Chem. 2004, 47: 1739-10.1021/jm0306430.View ArticleGoogle Scholar
  5. Moll A: J Comput Aided Mol Des. 2005, 19: 791-10.1007/s10822-005-9027-x.View ArticleGoogle Scholar


© Nina M et al; licensee BioMed Central Ltd. 2010

This article is published under license to BioMed Central Ltd.