Abstract
Tobacco etch virus protease (TEVp) is wildly exploited for various biotechnological applications. These applications take advantage of TEVp's ability to cleave specific substrate sequences to study protein function and interactions. A major limitation of this enzyme is its relatively slow catalytic rate. In this study, MD simulations were conducted on TEV enzymes and known highly active mutants (eTEV and uTEV3) to explore the relationship between mutation, conformation, and catalytic function. The results suggest that mutations distant from the active site can influence the substrate-binding pocket through interaction networks. MD analysis of eTEV demonstrates that, by stabilizing the orientation of the substrate at the catalytic site, mutations that appropriately enlarge the substrate-binding pocket will be beneficial for Kcat, enhancing the catalytic efficiency of the enzyme. On the contrary, mutations in uTEV3 reduced the flexibility of the active pocket and increased the hydrogen bonding between the substrate and enzyme, resulting in higher affinity. At the same time, the MD simulation demonstrates that mutations outside of the active site residues could affect the dynamic movement of the binding pocket by altering residue networks and communication pathways, thereby having a profound impact on reactivity. These findings not only provide a molecular mechanistic explanation for the excellent mutants, but also serve as a guiding framework for rational computational design.