Abstract
Protein kinases have high structural plasticity: their structure can change significantly, depending on what ligands are bound to them. Rigid‐protein docking methods are not capable of describing such effects. Here, we present a new flexible‐ligand flexible‐protein docking model in which the protein can adopt conformations between two extremes observed experimentally. The model utilized a molecular dynamics‐based simulated annealing cycling protocol and a distance‐dependent dielectric model to perform docking. By testing this model on docking four diverse ligands to protein kinase A, we found that the ligands were able to dock successfully to the protein with the proper conformations of the protein induced. By imposing relatively soft conformational restraints to the protein during docking, this model reduced computational costs yet permitted essential conformational changes that were essential for these inhibitors to dock properly to the protein. For example, without adequate movement of the glycine‐rich loop, it was difficult for the ligands to move from the surface of the protein to the binding site. In addition, these simulations called for better ways to compare simulation results with experiment other than using the popular root‐mean‐square deviation between the structure of a ligand in a docking pose and that in experiment because the structure of the protein also changed. In this work, we also calculated the correlation coefficient between protein–ligand/protein–protein distances in the docking structure and those in the crystal structure to check how well a ligand docked into the binding site of the protein and whether the proper conformation of the protein was induced. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009