The protein-protein interactions between SMPI and thermolysin studied by molecular dynamics and MM/PBSA calculations

Thermolysin is a zinc-metalloendopeptidase secreted by the gram-positive thermophilic bacterium Bacillus thermoproteolyticus. Thermolysin belongs to the gluzinicin family of enzymes, which is selectively inhibited by Steptomyces metalloproteinase inhibitor (SMPI). Very little is known about the interaction between SMPI and thermolysin. Knowledge about the protein-protein interactions is very important for designing new thermolysin inhibitors with possible industrial or pharmaceutical applications. In the present study, two binding modes between SMPI and thermolysin were studied by 2300 picoseconds (ps) of comparative molecular dynamics (MD) simulations and calculation of the free energy of binding using the molecular mechanics-Poisson-Boltmann surface area (MM/PBSA) method. One of the positions, the 'horizontal arrow head docking' (HAHD) was similar to the previously proposed binding mode by Tate et al. (Tate, S., Ohno, A., Seeram, S. S., Hiraga, K., Oda, K., and Kainosho, M. J. Mol. Biol. 282, 435-446 (1998)). The other position, the 'vertical arrow head docking' (VAHD) was obtained by a manual docking guided by the shape and charge distribution of SMPI and the binding pocket of thermolysin. The calculations showed that SMPI had stronger interactions with thermolysin in the VAHD than in the HAHD complex, and the VAHD complex was considered more realistic than the HAHD complex. SMPI interacted with thermolysin not only at the active site but had auxiliary binding sites contributing to proper interactions. The VAHD complex can be used for designing small molecule inhibitors mimicking the SMPI-thermolysin binding interfaces.