Probing the importance of hydrogen bonds in the active site of the subtilisin nattokinase by site-directed mutagenesis and molecular dynamics simulation

Hydrogen bonds occurring in the catalytic triad (Asp 32 , His 64 and Ser 221 ) and the oxyanion hole (Asn 155 ) are very important to the catalysis of peptide bond hydrolysis by serine proteases. For the subtilisin NK (nattokinase), a bacterial serine protease, construction and analysis of a three-dimensional structural model suggested that several hydrogen bonds formed by four residues function to stabilize the transition state of the hydrolysis reaction. These four residues are Ser 33 , Asp 60 , Ser 62 and Thr 220 . In order to remove the effect of these hydrogen bonds, four mutants (Ser 33 →Ala 33 , Asp 60 →Ala 60 , Ser 62 →Ala 62 , and Thr 220 →Ala 220 ) were constructed by site-directed mutagenesis. The results of enzyme kinetics indicated that removal of these hydrogen bonds increases the free-energy of the transition state (ΔΔ G T ). We concluded that these hydrogen bonds are more important for catalysis than for binding the substrate, because removal of these bonds mainly affects the k cat but not the K m values. A substrate, SUB1 (succinyl-Ala-Ala-Pro-Phe- p -nitroanilide), was used during enzyme kinetics experiments. In the present study we have also shown the results of FEP (free-energy perturbation) calculations with regard to the binding and catalysis reactions for these mutant subtilisins. The calculated difference in FEP also suggested that these four residues are more important for catalysis than binding of the substrate, and the simulated values compared well with the experimental values from enzyme kinetics. The results of MD (molecular dynamics) simulations further demonstrated that removal of these hydrogen bonds partially releases Asp 32 , His 64 and Asn 155 so that the stability of the transition state decreases. Another substrate, SUB2 (H-D-Val-Leu-Lys- p -nitroanilide), was used for FEP calculations and MD simulations.