PDBsum entry 2rb2

Hydrolase

PDB id

2rb2

Contents

			Protein chain

					162 a.a.

			Ligands

					PO4 ×3


					263

			Waters ×320

References listed in PDB file

Key reference

Title

Rescoring docking hit lists for model cavity sites: predictions and experimental testing.

Authors

A.P.Graves, D.M.Shivakumar, S.E.Boyce, M.P.Jacobson, D.A.Case, B.K.Shoichet.

Ref.

J Mol Biol, 2008, 377, 914-934. [DOI no: 10.1016/j.jmb.2008.01.049]

PubMed id

18280498

Abstract

Molecular docking computationally screens thousands to millions of organic molecules against protein structures, looking for those with complementary fits. Many approximations are made, often resulting in low "hit rates." A strategy to overcome these approximations is to rescore top-ranked docked molecules using a better but slower method. One such is afforded by molecular mechanics-generalized Born surface area (MM-GBSA) techniques. These more physically realistic methods have improved models for solvation and electrostatic interactions and conformational change compared to most docking programs. To investigate MM-GBSA rescoring, we re-ranked docking hit lists in three small buried sites: a hydrophobic cavity that binds apolar ligands, a slightly polar cavity that binds aryl and hydrogen-bonding ligands, and an anionic cavity that binds cationic ligands. These sites are simple; consequently, incorrect predictions can be attributed to particular errors in the method, and many likely ligands may actually be tested. In retrospective calculations, MM-GBSA techniques with binding-site minimization better distinguished the known ligands for each cavity from the known decoys compared to the docking calculation alone. This encouraged us to test rescoring prospectively on molecules that ranked poorly by docking but that ranked well when rescored by MM-GBSA. A total of 33 molecules highly ranked by MM-GBSA for the three cavities were tested experimentally. Of these, 23 were observed to bind--these are docking false negatives rescued by rescoring. The 10 remaining molecules are true negatives by docking and false positives by MM-GBSA. X-ray crystal structures were determined for 21 of these 23 molecules. In many cases, the geometry prediction by MM-GBSA improved the initial docking pose and more closely resembled the crystallographic result; yet in several cases, the rescored geometry failed to capture large conformational changes in the protein. Intriguingly, rescoring not only rescued docking false positives, but also introduced several new false positives into the top-ranking molecules. We consider the origins of the successes and failures in MM-GBSA rescoring in these model cavity sites and the prospects for rescoring in biologically relevant targets.



	Figure 3. Fig. 3. Predicted and experimental ligand orientations for the hydrophobic L99A cavity. The carbon atoms of the crystallographic pose, the DOCK predicted pose, the AMBERDOCK predicted pose, and the PLOP predicted pose are colored gray, yellow, cyan, and magenta, respectively. The F[o] − F[c] omit electron density maps (green mesh) are contoured at 2.5–3.0σ (a) β-chlorophenetole (1), (b) 4-(methylthio)nitrobenzene (2), (c) 2,6-difluorobenzylbromide (4), (d) 2-ethoxyphenol (5), and (e) 3-methylbenzylazide (6) bound to L99A.		Figure 4. Fig. 4. Predicted and experimental ligand orientations for the polar L99A/M102Q cavity site. The carbon atoms of the crystallographic, DOCK, AMBERDOCK, and PLOP predicted poses are colored gray, yellow, cyan, and magenta, respectively. Hydrogen bonds are depicted with dashed lines. The F[o] − F[c] electron density omit maps (green mesh) are contoured at 2.5–3.0σ. (a) n-Phenylglycinonitrile (10), (b) 2-nitrothiophene (11), (c) 2-(n-propylthio)ethanol (12), (d) 3-methylbenzylazide (6), (e) 2-phenoxyethanol (9), and (f) (R)-(+)-3-chloro-1-phenyl-1-propanol (13) bound to L99A/M102Q.

PROCHECK

	Headers