The hydrolysis by thermolysin of a Gly-Phe-Leu peptide, considered as a mod
el substrate of the enkephalin family, has been studied with a mixed QM/MM
method with the AM1/AMBER parameterization. This study is based on the mech
anism proposed by Matthews in which the Glu-143 residue plays the role of a
proton shuttle in the course of the reaction. The study focused on the des
cription of every step of the process, reaction intermediates and transitio
n states, and on the influence, both energetical and structural, of the who
le protein on these stationary points. The overall mechanism appears to be
quite realistic, but the study shows that some reaction steps that were ass
umed to be concerted should occur in two phases. Analysis of the role of th
e amino-acids surrounding the active site has shown their important partici
pation in the fluctuations of the energy. In particular, the major role of
His-231 on the overall mechanism has been confirmed. This study shows that
modeling reaction mechanisms for enzymes is quite feasible and opens the wa
y for computer experiments that may be helpful in devising and interpreting
detailed experimental investigations.