Investigation of the hydroxylation mechanism of noncoupled copper oxygenases by ab initio molecular dynamics simulations.

In Nature, the family of copper monooxygenases comprised of peptidylglycine α-hydroxylating monooxygenase (PHM), dopamine β-monooxygenase (DβM), and tyramine β-monooxygenase (TβM) is known to perform dioxygen-dependent hydroxylation of aliphatic C-H bonds by using two uncoupled metal sites. In spite of many investigations, including biochemical, chemical, and computational, details of the C-H bond oxygenation mechanism remain elusive. Herein we report an investigation of the mechanism of hydroxylation by PHM by using hybrid quantum/classical potentials (i.e., QM/MM). Although previous investigations using hybrid QM/MM techniques were restricted to geometry optimizations, we have carried out ab initio molecular dynamics simulations in order to include the intrinsic flexibility of the active sites in the modeling protocol. The major finding of this study is an extremely fast rebound step after the initial hydrogen-abstraction step promoted by the cupric-superoxide adduct. The hydrogen-abstraction/rebound sequence leads to the formation of an alkyl hydroperoxide intermediate. Long-range electron transfer from the remote copper site subsequently triggers its reduction to the hydroxylated substrate. We finally show two reactivity consequences inherent in the new mechanistic proposal, the investigation of which would provide a means to check its validity by experimental means.