Specific interactions and binding free energies between thermolysin and dipeptides: molecular simulations combined with ab initio molecular orbital and classical vibrational analysis

Thermolysin (TLN) is a metalloprotease widely used as a nonspecific protease for sequencing peptide and synthesizing many useful chemical compounds by the chemical industry. It was experimentally shown that the activity and functions of TLN are inhibited by the binding of many types of amino acid dipeptides. However, the binding mechanisms between TLN and dipeptides have not been clarified at the atomic and electronic levels. In this study, we investigated the binding mechanisms between TLN and four dipeptides. Specific interactions and binding free energies (BFEs) between TLN and the dipeptides were calculated using molecular simulations based on classical molecular dynamics and ab initio fragment molecular orbital (FMO) methods. The molecular systems were embedded in solvating water molecules during calculations. The calculated BFEs were qualitatively consistent with the trend of the experimentally observed inhibition of TLN activity by binding of the dipeptides. In addition, the specific interactions between the dipeptides and each amino acid residue of TLN or solvating water molecules were elucidated by the FMO calculations.