Theoretical modeling of the hydroxylation of methane as mediated by the particulate methane monooxygenase.

We present here the results of density functional theory (DFT) calculations directed toward elucidation of the CH bond activation mechanism that might be adopted by the particulate methane monooxygenase (pMMO) in the hydroxylation of methane and related small alkanes. In these calculations, we considered three of the most probable models for the transition metal active site mediating the "oxo-transfer": (i) the trinuclear copper cluster bis(mu(3)-oxo)trinuclear copper(II, II, III) complex 1, recently proposed by Chan et al. [S.I. Chan, K.H.-C. Chen, S.S.-F. Yu, C.-L. Chen, S.S.-J. Kuo, Biochemistry 43 (2004) 4421-4430.]; (ii) the most frequently used model complex, bis(mu-oxo)Cu(III)(2) complex 2; and (iii) the mixed-valence bis(mu-oxo)Cu(II)Cu(III) complex 3. The results obtained indicate that the methane hydroxylation chemistry mediated by the trinuclear copper cluster bis(mu(3)-oxo)trinuclear copper(II, II, III) complex 1 offers the most facile pathway for methane hydroxylation, and this model yields KIE values that are in good agreement with experiment. In this mechanism, the reaction proceeds along a "singlet" potential surface and a "singlet oxene" is directly inserted across a CH bond in a concerted manner. Kinetic isotope effects (k(H)/k(D) or KIE) associated with the concerted oxene insertion process mediated by complex 1 are calculated to be 5.2 at 300K when tunneling effects are included. Overall rate constants for the methane hydroxylation by the three models have been calculated as a function of temperature, and the rates are at least 5-6 orders of magnitude more facile when the chemistry is mediated by complex 1 compared to complex 2 or complex 3.