Methyl-coenzyme M reductase (MCR) is a central player in methane biogeochemistry, governing methanogenesis and the anaerobic oxidation of methane (AOM) in methanogens and anaerobic methanotrophs (ANME), respectively. The prosthetic group of MCR is coenzyme F430, a nickel-containing tetrahydrocorphin. Several modified versions of F430 have been discovered, including the 172-methylthio-F430 (mtF430) used by ANME-1 MCR. Here, we employ molecular dynamics (MD) simulations to investigate the active site dynamics of MCR from Methanosarcina acetivorans and ANME-1 when bound to the canonical F430 compared to 172-thioether coenzyme F430 variants and substrates (methyl-coenzyme M and coenzyme B) for methane formation. Our simulations highlight the importance of the Gln to Val substitution in accommodating the 172 methylthio modification in ANME-1 MCR. Modifications at the 172 position disrupt the canonical substrate positioning in M. acetivorans MCR. However, in some replicates, active site reorganization to maintain substrate positioning suggests that the modified F430 variants could be accommodated in a methanogenic MCR. We additionally report the first quantitative estimate of MCR intrinsic electric fields that are pivotal in driving methane formation. Our results suggest that the electric field aligned along the CH3-S-CoM thioether bond facilitates homolytic bond cleavage, coinciding with the proposed catalytic mechanism. Structural perturbations, however, weaken and misalign these electric fields, emphasizing the importance of the active site structure in maintaining their integrity. In conclusion, our results deepen the understanding of MCR active site dynamics, the enzyme's organizational role in intrinsic electric fields for catalysis, and the interplay between active site structure and electrostatics.
Read full abstract