How Modification of Accessible Lysines to Phenylalanine Modulates the Structural and Functional Properties of Horseradish Peroxidase: A Simulation Study

Horseradish Peroxidase (HRP) is one of the most studied peroxidases and a great number of chemical modifications and genetic manipulations have been carried out on its surface accessible residues to improve its stability and catalytic efficiency necessary for biotechnological applications. Most of the stabilized derivatives of HRP reported to date have involved chemical or genetic modifications of three surface-exposed lysines (K174, K232 and K241). In this computational study, we altered these lysines to phenylalanine residues to model those chemical modifications or genetic manipulations in which these positively charged lysines are converted to aromatic hydrophobic residues. Simulation results implied that upon these substitutions, the protein structure becomes less flexible. Stability gains are likely to be achieved due to the increased number of stable hydrogen bonds, improved heme-protein interactions and more integrated proximal Ca²⁺ binding pocket. We also found a new persistent hydrogen bond between the protein moiety (F174) and the heme prosthetic group as well as two stitching hydrogen bonds between the connecting loops GH and F′F″ in mutated HRP. However, detailed analysis of functionally related structural properties and dynamical features suggests reduced reactivity of the enzyme toward its substrates. Molecular dynamics simulations showed that substitutions narrow the bottle neck entry of peroxide substrate access channel and reduce the surface accessibility of the distal histidine (H42) and heme prosthetic group to the peroxide and aromatic substrates, respectively. Results also demonstrated that the area and volume of the aromatic-substrate binding pocket are significantly decreased upon modifications. Moreover, the hydrophobic patch functioning as a binding site or trap for reducing aromatic substrates is shrunk in mutated enzyme. Together, the results of this simulation study could provide possible structural clues to explain those experimental observations in which the protein stability achieved concurrent with a decrease in enzyme activity, upon manipulation of charge/hydrophobicity balance at the protein surface.