Abstract
Activation and reaction energies for four model systems capturing the essential physicochemical features of the hydrolysis of the peptide bond have been calculated at various level of theory, including the presumably accurate CCSD(T) calculations. The models studied covered a part of the spectrum encountered in biological systems: the hydrolysis in the absence of metal ions (represented by formamide and Ala-Ala) and the hydrolysis in the presence of one and two zinc(II) ions, mimicking the active sites of mono- and dizinc metallopeptidases, respectively (by using thermolysin and glutamate carboxypeptidaseII as the model catalytic systems and formamide as the model substrate).
The results obtained using CCSD(T)/def2-TZVP and CCSD(T)/aug-cc-pVTZ calculations were used as the benchmark values to which the set of cheaper methods, such as (RI-)DFT, (RI-)MP2, and SCS-MP2, were referenced. It was shown that deviations of 3-5kcalmol(-1) (translating to 2-3 orders in reaction constants) with respect to the reference CCSD(T) barriers are frequently encountered for many correlated methods and most of studied DFT functionals.
It has been concluded that from the set of wave-function methods, both MP2 and SCS-MP2 methods can be recommended for smaller models (measured by the mean absolute deviation of the activation barriers over the four systems studied), whereas among the popular DFT functionals, B3LYP and especially M06-2X are likely to be reasonable choices for calculating the activation barriers of zinc metallopeptidases. Finally, with the model of glutamate carboxypeptidaseII, issues related to the convergence of the calculated barriers with the size of the model system used as the representative of the enzyme active site were addressed.