QM methods in structure based design: utility in probing protein-ligand interactions

Small changes in ligand structure can lead to large unexpected changes in activity yet it is often not possible to rationalize these effects using empirical modeling techniques, suggesting more effective methods are required. In this study we investigate the use of high level QM methods to study the interactions found within protein-ligand complexes as improved understanding of these could help in the design of new, more active molecules. We study aspects of ligand binding in a set of protein ligand complexes containing ligand efficient, fragment-like inhibitors as these structures are often challenging to determine experimentally. To assess the reliability of our theoretical models we compare the MP2/6-31+G** QM results to the original X-ray coordinates and to QM/MM B3LYP/6-31G*//UFF results which we have previously reported. We also contrast these results with data obtained from an analysis of the distribution of comparable interactions found in (a) high resolution kinase complexes (≤ 1.8Å) from the PDB and (b) more generic, small molecule crystal structures from the CSD.