Theoretical studies of the proton transfer behaviors in molecular complexes analogous to catalytic triad of serine protease: Toward understanding the existence and significance of the low-barrier hydrogen-bond in enzymatic catalysis

A representative acetate-(5-methylimidazole)-methanol system has been employed as a model of catalytic triad in serine protease to validate the formation processes of low-barrier H-bonds (LBHB) at the B3LYP/6-311++G** level of theory, and variable H-bonding characters from conventional ones to LBHBs have been represented along with the proceedings of proton transfer. Solvent effect is an important factor in modulation of the existence of an LBHB, where an LBHB (or a conventional H-bond) in the gas phase can be changed into a non-LBHB (an LBHB) upon solvation. The origin of the additional stabilization energy arising from the LBHB may be attributed to the H-bonding energy difference before and after proton transfer because the shared proton can freely move between the proton donor and proton acceptor. Most importantly, the order of magnitude of the stabilization energy depends on the studied systems. Furthermore, the nonexistence of LBHBs in the catalytic triad of serine proteases has been verified in a more sophisticated model treated using the ONIOM method. As a result, only the single proton transfer mechanism in the catalytic triad has been confirmed and the origin of the powerful catalytic efficiency of serine proteases should be attributed to other factors rather than the LBHB.