Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which 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 (UniProtKB Q9X0C6). 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 using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's 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 μ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 of the receiver operating characteristics 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.