The role of hydrogen bond in catalytic triad of serine proteases ^†

In order to investigate the origin of catalytic power for serine proteases, the role of the hydrogen bond in the catalytic triad was studied in the proteolysis process of the peptides chymotrypsin inhibitor 2 (CI2), MCTI-A, and a hexapeptide (SUB), respectively. We first calculated the free energy profile of the proton transfer between His and Asp residues of the catalytic triad in the enzyme-substrate state and transition state by employing QM/MM molecular dynamics simulations. The results show that a low-barrier hydrogen bond (LBHB) only forms in the transition state of the acylation of CI2, while it is a normal hydrogen bond in the acylation of MCTI-A or SUB. In addition, the change of the hydrogen bond strength is much larger in CI2 and SUB systems than in MCTI-A system, which decreases the acylation energy barrier significantly for CI2 and SUB. Clearly, a LBHB formed in the transition state region helps accelerate the acylation reaction. But to our surprise, a normal hydrogen bond can also help to decrease the energy barrier. The key to reducing the reaction barrier is the increment of hydrogen bond strength in the transition state state, whether it is a LBHB or not. Our studies cast new light on the role of the hydrogen bond in the catalytic triad, and help to understand the catalytic triad of serine proteases.