Abstract

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da that weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF. We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability by RosettaLigand, and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigands ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 {micro}M. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve (AUC) of the receiver operating characteristics (ROC) of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.