Abstract
DockIT is a tool for interactive molecular docking that can model both the local and global conformational response of the receptor to the docking of a ligand based on information from a molecular dynamics simulation. Using DockIT we have investigated the binding process of two histone deacetylase (HDAC) inhibitors to HDAC6: the nonselective approved drug belinostat and the preclinical HPOB. To model HDAC6's conformational response to the binding of the inhibitors we performed a 200-nanosecond explicit-solvent molecular dynamics simulation on HDAC6. Unexpectedly the simulation revealed a domain movement that affects the size and shape of the binding pocket. Using automated docking and a rigid model for the inhibitors, the domain movement continuously adapts the pocket to the presence of the inhibitor. For both inhibitors, an intermediate binding site was found where it was partially inserted, with a hydrogen bond formed between the inhibitor's hydroxamic acid and the Tyr745 side chain. Pushing the inhibitor deeper into the pocket over an energy barrier and re-engaging automated docking, a final binding pose resulted with a root-mean square deviation with its respective crystallographic pose of 1.0 Å for belinostat and 1.4 Å for HPOB. We believe our results mimic substrate recognition by the enzyme, with an initial partial binding of the acetyllysine residue with Tyr745. During binding a relay of hydrogen bonds occurs coordinating the orientation of the cap and the hydroxamic acid inside the pocket. The interaction between the cap and the surface of HDAC6 explains the reason for the hydroxamic acid warhead in HPOB binding in a flipped orientation compared to belinostat.