Unspecific peroxygenases (UPOs) can
selectively oxyfunctionalize
unactivated hydrocarbons by using peroxides under mild conditions.
They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit
greater stability than the latter. As such, they hold great potential
for industrial applications. A thorough understanding of their catalysis
is needed to improve their catalytic performance. However, it remains
elusive how UPOs effectively convert peroxide to Compound I (CpdI),
the principal oxidizing intermediate in the catalytic cycle. Previous
computational studies of this process primarily focused on heme peroxidases
and P450s, which have significant differences in the active site from
UPOs. Additionally, the roles of peroxide unbinding in the kinetics
of CpdI formation, which is essential for interpreting existing experiments,
have been understudied. Moreover, there has been a lack of free energy
characterizations with explicit sampling of protein and hydration
dynamics, which is critical for understanding the thermodynamics of
the proton transport (PT) events involved in CpdI formation. To bridge
these gaps, we employed multiscale simulations to comprehensively
characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive
free energy and potential energy calculations were performed in a
quantum mechanics/molecular mechanics setting. Our results indicate
that substrate-binding dehydrates the active site, impeding the PT
from H2O2 to a nearby catalytic base (Glu196).
Furthermore, the PT is coupled with considerable hydrogen bond network
rearrangements near the active site, facilitating subsequent O–O
bond cleavage. Finally, large unbinding free energy barriers kinetically
stabilize H2O2 at the active site. These findings
reveal a delicate balance among PT, hydration dynamics, hydrogen bond
rearrangement, and cosubstrate unbinding, which collectively enable
efficient CpdI formation. Our simulation results are consistent with
kinetic measurements and offer new insights into the CpdI formation
mechanism at atomic-level details, which can potentially aid the design
of next-generation biocatalysts for sustainable chemical transformations
of feedstocks.