Probing the conformational space available to inhibitors in the thermolysin active site using Monte Carlo/energy minimization techniques

The conformational space available to four inhibitors of the bacterial enzyme thermolysin has been searched in the enzyme binding site using a method that combines Monte Carlo type techniques with energy minimization for exploration of the conformational potential energy hypersurface. Molecular mechanics methodology using the AMBER force field was employed for computation of the molecular energetics. Solvation energies were also included in the calculations by employing a technique that estimates hydration energies based on the exposed solvent accessible surface area for each atom of the inhibitor and active site. It was found that in each case, the crystallographically observed conformation was among the low energy conformers discovered. In fact, in three of the calculations it was the lowest energy conformation. The methodology described in this article is expected to be quite useful for studies involving computer aided design and evaluation of enzyme inhibitors.