Catalysis by methyl-coenzyme M reductase: a theoretical study for heterodisulfide product formation

Hybrid density functional theory has been used to investigate the catalytic mechanism of methyl-coenzyme M reductase (MCR), an essential enzyme in methanogenesis. In a previous study of methane formation, a scheme was suggested involving oxidation of Ni(I) in the starting square-planar coordination to the high-spin Ni(II) form in the CoM-S-Ni(II)F(430) octahedral intermediate. The methyl radical, concomitantly released by methyl-coenzyme M (CoM), is rapidly quenched by hydrogen atom transfer from the coenzyme B (CoB) thiol group, yielding methane as the first product of the reaction. The present investigation primarily concerns the second and final step of the reaction: oxidation of CoB and CoM to the CoB-S-S-CoM heterodisulfide product and reduction of nickel back to the Ni(I) square-planar form. The activation energy for the second step is found to be around 10 kcal/mol, implying that the first step of methane formation with an activation energy of 20 kcal/mol should be rate-limiting. An oxygen of the Gln147 residue, occupying the rear axial position in the oxidized Ni(II) state, is shown to stabilize the intermediate by 6 kcal/mol, thereby slightly decreasing the barrier for the preceding rate-limiting transition state. The mechanism suggested is discussed in the context of available experimental data. An analysis of the flexibility of the F(430) cofactor during the reaction cycle is also given.