Computational analysis for selectivity of histone deacetylase inhibitor by replica-exchange umbrella sampling molecular dynamics simulations

We performed protein-ligand docking simulations with a ligand T247, which has been reported as a selective inhibitor of a histone deacetylase HDAC3, by the replica-exchange umbrella sampling method in order to estimate the free energy profiles along ligand docking pathways of HDAC3-T247 and HDAC2-T247 systems. The simulation results showed that the docked state of the HDAC3-T247 system is more stable than that of the HDAC2-T247 system although the amino-acid sequences and structures of HDAC3 and HDAC2 are very similar. By comparing structures obtained from the simulations of both systems, we found the difference between structures of hydrophobic residues at the entrance of the catalytic site. Moreover, we performed conventional molecular dynamics simulations of HDAC3 and HDAC2 systems without T247, and the results also showed the same difference of the hydrophobic structures. Therefore, we consider that this hydrophobic structure contributes to the stabilization of the docked state of the HDAC3-T247 system. Furthermore, we show that Tyr209, which is one of the hydrophobic residues in HDAC2, plays a key role in the instability from the simulation results of a mutated-HDAC2 system.