Structural Dynamics of the Methyl-Coenzyme M Reductase Active Site Are Influenced by Coenzyme F430 Modifications

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.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Mononuclear nickel(ii)-flavonolate complexes of tetradentate tripodal 4N ligands as structural and functional models for quercetin 2,4-dioxygenase: structures, spectra, redox and dioxygenase activity

Three new nickel(ii)-flavonolate complexes of the type [Ni(L)(fla)](ClO4) 1-3, where L is the tripodal 4N ligand tris(pyrid-2-ylmethyl)amine (tpa, L1) or (pyrid-2-ylmethyl)bis(6-methylpyrid-2-ylmethyl)amine (6-Me2-tpa, L2) or tris(N-Et-benzimidazol-2-ylmethyl)amine (Et-ntb, L3), have been isolated as functional models for Ni(ii)-containing quercetin 2,4-dioxygenase. Single crystal X-ray structures of 1 and 3 reveal that Ni(ii) is involved in π-back bonding with flavonolate (fla-), as evident from enhancement in C[double bond, length as m-dash]O bond length upon coordination [H(fla), 1.232(3); 1, 1.245(7); 3, 1.262(8) Å]. More asymmetric chelation of fla- in 3 than in 1 [Δd = (Ni-Ocarbonyl - Ni-Oenolate): 1, 0.126; 3, 0.182 Å] corresponds to lower π-delocalization in 3 with electron-releasing N-Et substituent. The optimized structures of 1-3 and their geometrical isomers have been computed by DFT methods. The HOMO and LUMO, both localized on Ni(ii)-bound fla-, are highly conjugated bonding π- and antibonding π*-orbitals respectively. They are located higher in energy than the Ni(ii)-based MOs (HOMO-1, dx2-y2; HOMO-2/6, dz2), revealing that the Ni(ii)-bound fla- rather than Ni(ii) would undergo oxidation upon exposure to dioxygen. The results of computational studies, in combination with spectral and electrochemical studies, support the involvement of redox-inactive Ni(ii) in π-back bonding with fla-, tuning the π-delocalization in fla- and hence its activation. Upon exposure to dioxygen, all the flavonolate adducts in DMF solution decompose to produce CO and depside, which then is hydrolyzed to give the corresponding acids at 70 °C. The highest rate of dioxygenase reactivity of 3 (kO2: 3 (29.10 ± 0.16) > 1 (16.67 ± 0.70) > 2 (1.81 ± 0.04 × 10-1 M-1 s-1)), determined by monitoring the disappearance of the LMCT band in the range 440-450 nm, is ascribed to the electron-releasing N-Et substituent on bzim ring, which decreases the π-delocalization in fla- and enhances its activation.

Reaction mechanism of the PuDddK dimethylsulfoniopropionate lyase and cofactor effects of various transition metal ions

The microbial cleavage of dimethylsulfoniopropionate (DMSP) produces volatile dimethyl sulfide (DMS) via the lyase pathway, playing a crucial role in the global sulfur cycle. Herein, the DMSP decomposition catalyzed by PuDddK (a DMSP lyase) devised with various transition metal ion cofactors are investigated using density functional calculations. The PuDddK reaction has been demonstrated to employ a concerted β-elimination mechanism, where the substrate α-proton abstraction by the deprotonated Tyr64 occurs simultaneously with the C_β-S bond cleavage and C_α = C_β double bond formation. The PuDddK enzymes with diverse divalent metal ions (Ni²⁺, Mn²⁺, Fe²⁺, Co²⁺, Zn²⁺, and Cu²⁺) incorporated prefer DMSP as a monodentate ligand. The cases of Ni²⁺, Mn²⁺, Fe²⁺, Co²⁺, and Zn²⁺ with the same 3His-1Glu ligands have close reaction energy barriers, indicating that the lyase activity may be hardly affected by the divalent transition metal type with the same ligand type and number. The coordination loss of one histidine in Cu²⁺, forming a 2His-1Glu architecture, leads to a lower activity, revealing that the 3His-1Glu ligand set used by DddK appears to be a scaffold capable of more efficiently catalyzing the DMSP decomposition. Further analysis reveals that the inactivation of Fe³⁺-dependent PuDddK is derived from an electron transfer from the Tyr64 phenolate to Fe³⁺, with the implication that the PuDddK activity may be primarily affected by the redox effects induced by a strongly oxidizing transition metal ion (like Fe³⁺).

Overview of Diverse Methyl/Alkyl-Coenzyme M Reductases and Considerations for Their Potential Heterologous Expression

Methyl-coenzyme M reductase (MCR) is an archaeal enzyme that catalyzes the final step of methanogenesis and the first step in the anaerobic oxidation of methane, the energy metabolisms of methanogens and anaerobic methanotrophs (ANME), respectively. Variants of MCR, known as alkyl-coenzyme M reductases, are involved in the anaerobic oxidation of short-chain alkanes including ethane, propane, and butane as well as the catabolism of long-chain alkanes from oil reservoirs. MCR is a dimer of heterotrimers (encoded by mcrABG) and requires the nickel-containing tetrapyrrole prosthetic group known as coenzyme F430. MCR houses a series of unusual post-translational modifications within its active site whose identities vary depending on the organism and whose functions remain unclear. Methanogenic MCRs are encoded in a highly conserved mcrBDCGA gene cluster, which encodes two accessory proteins, McrD and McrC, that are believed to be involved in the assembly and activation of MCR, respectively. The requirement of a unique and complex coenzyme, various unusual post-translational modifications, and many remaining questions surrounding assembly and activation of MCR largely limit in vitro experiments to native enzymes with recombinant methods only recently appearing. Production of MCRs in a heterologous host is an important step toward developing optimized biocatalytic systems for methane production as well as for bioconversion of methane and other alkanes into value-added compounds. This review will first summarize MCR catalysis and structure, followed by a discussion of advances and challenges related to the production of diverse MCRs in a heterologous host.

Focusing on a nickel hydrocorphinoid in a protein matrix: methane generation by methyl-coenzyme M reductase with F430 cofactor and its models

Methyl-coenzyme M reductase (MCR) containing a nickel hydrocorphinoid cofactor, F430, is an essential enzyme that catalyzes anaerobic methane generation and oxidation. The active Ni(I) species in MCR converts methyl-coenzyme M (CH₃S-CoM) and coenzyme B (HS-CoB) to methane and heterodisulfide (CoM-S-S-CoB). Extensive experimental and theoretical studies focusing on the substrate-binding cavity including the F430 cofactor in MCR have suggested two principally different reaction mechanisms involving an organonickel CH₃-Ni(III) species or a transient methyl radical species. In parallel with research on native MCR itself, the functionality of MCR has been investigated in the context of model complexes of F430 and recent protein-based functional models, which include a nickel complex. In the latter case, hemoproteins reconstituted with tetradehydro- and didehydrocorrinoid nickel complexes have been found to represent useful model systems that are responsible for methane generation. These efforts support the proposed mechanism of the enzymatic reaction and provide important insight into replicating the MCR-like methane-generation process. Furthermore, the modeling of MCR described here is expected to lead to understanding of protein-supported nickel porphyrinoid chemistry as well as the creation of MCR-inspired catalysis.

Handling methane: a Ni(i) F430-like cofactor derived from VB12 is active in methyl-coenzyme M reductase

An Ni(i) F₄₃₀-like cofactor derived from vitamin B₁₂ can catalyze methane formation in the active site of methyl-coenzyme M reductase.

Methane Generation and Reductive Debromination of Benzylic Position by Reconstituted Myoglobin Containing Nickel Tetradehydrocorrin as a Model of Methyl-coenzyme M Reductase

Methyl-coenzyme M reductase (MCR), which contains the nickel hydrocorphinoid cofactor F430, is responsible for biological methane generation under anaerobic conditions via a reaction mechanism which has not been completely elucidated. In this work, myoglobin reconstituted with an artificial cofactor, nickel(I) tetradehydrocorrin (Ni^I(TDHC)), is used as a protein-based functional model for MCR. The reconstituted protein, rMb(Ni^I(TDHC)), is found to react with methyl donors such as methyl p-toluenesulfonate and trimethylsulfonium iodide with methane evolution observed in aqueous media containing dithionite. Moreover, rMb(Ni^I(TDHC)) is found to convert benzyl bromide derivatives to reductively debrominated products without homocoupling products. The reactivity increases in the order of primary > secondary > tertiary benzylic carbons, indicating steric effects on the reaction of the nickel center with the benzylic carbon in the initial step. In addition, Hammett plots using a series of para-substituted benzyl bromides exhibit enhancement of the reactivity with introduction of electron-withdrawing substituents, as shown by the positive slope against polar substituent constants. These results suggest a nucleophilic S_N2-type reaction of the Ni(I) species with the benzylic carbon to provide an organonickel species as an intermediate. The reaction in D₂O buffer at pD 7.0 causes a complete isotope shift of the product by +1 mass unit, supporting our proposal that protonation of the organonickel intermediate occurs during product formation. Although the turnover numbers are limited due to inactivation of the cofactor by side reactions, the present findings will contribute to elucidating the reaction mechanism of MCR-catalyzed methane generation from activated methyl sources and dehalogenation.

Myoglobin Reconstituted with Ni Tetradehydrocorrin as a Methane-Generating Model of Methyl-coenzyme M Reductase

Myoglobin reconstituted with Ni tetradehydrocorrin was investigated as a model of F430-containing methyl-coenzyme M reductase, which catalyzes anaerobic methane generation. The Ni^II tetradehydrocorrin complex has a Ni^II /Ni^I redox potential of -0.34 V vs. SHE and EPR spectroscopy indicates the formation of a Ni^I species upon reduction by dithionite. This redox potential is approximately 0.31 V more positive than that of F430. The Ni^I tetradehydrocorrin moiety is bound to the apo-form of myoglobin to yield the reconstituted protein. Methane gas is generated in the reaction of the model with methyl iodide in the presence of the reconstituted protein under reductive conditions, whereas the Ni^I complex itself does not produce methane gas. This is the first example of a protein-based functional model of F430-containing methyl-coenzyme M reductase.

Theoretical Studies of Nickel-Dependent Enzymes

The advancements of quantum chemical methods and computer power allow detailed mechanistic investigations of metalloenzymes. In particular, both quantum chemical cluster and combined QM/MM approaches have been used, which have been proven to successfully complement experimental studies. This review starts with a brief introduction of nickel-dependent enzymes and then summarizes theoretical studies on the reaction mechanisms of these enzymes, including NiFe hydrogenase, methyl-coenzyme M reductase, nickel CO dehydrogenase, acetyl CoA synthase, acireductone dioxygenase, quercetin 2,4-dioxygenase, urease, lactate racemase, and superoxide dismutase.

Methyl (Alkyl)-Coenzyme M Reductases: Nickel F-430-Containing Enzymes Involved in Anaerobic Methane Formation and in Anaerobic Oxidation of Methane or of Short Chain Alkanes

Methyl-coenzyme M reductase (MCR) catalyzes the methane-forming step in methanogenic archaea. The active enzyme harbors the nickel(I) hydrocorphin coenzyme F-430 as a prosthetic group and catalyzes the reversible reduction of methyl-coenzyme M (CH₃–S-CoM) with coenzyme B (HS-CoM) to methane and CoM-S–S-CoB. MCR is also involved in anaerobic methane oxidation in reverse of methanogenesis and most probably in the anaerobic oxidation of ethane, propane, and butane. The challenging question is how the unreactive CH₃–S thioether bond in methyl-coenzyme M and the even more unreactive C–H bond in methane and the other hydrocarbons are anaerobically cleaved. A key to the answer is the negative redox potential (E_o′) of the Ni(II)F-430/Ni(I)F-430 couple below −600 mV and the radical nature of Ni(I)F-430. However, the negative one-electron redox potential is also the Achilles heel of MCR; it makes the nickel enzyme one of the most O₂-sensitive enzymes known to date. Even under physiological conditions, the Ni(I) in MCR is oxidized to the Ni(II) or Ni(III) states, e.g., when in the cells the redox potential (E′) of the CoM-S–S-CoB/HS-CoM and HS-CoB couple (E_o′ = −140 mV) gets too high. Methanogens therefore harbor an enzyme system for the reactivation of inactivated MCR in an ATP-dependent reduction reaction. Purification of active MCR in the Ni(I) oxidation state is very challenging and has been achieved in only a few laboratories. This perspective reviews the function, structure, and properties of MCR, what is known and not known about the catalytic mechanism, how the inactive enzyme is reactivated, and what remains to be discovered.

Insights into the dioxygen activation and catalytic mechanism of the nickel-containing quercetinase

The catalytic mechanism of Ni-QueD^FLA was elucidated by QM/MM calculations, and the different reactivities of nickel and iron were illuminated.

From NAD+ to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase

Lactate racemase (LarA), a new nickel enzyme discovered recently, catalyzes the racemization between d- and l-lactates with a novel nickel pincer cofactor (Ni-PTTMN) derived from nicotinic acid. In this study, by using DFT and a 200-atom active-site model, LarA is revealed to employ a modified proton-coupled hydride-transfer mechanism in which a hydride is transferred to a cofactor pyridine carbon from the substrate α-carbon along with proton transfer from the substrate hydroxy group to a histidine, and then moved back from the opposite side. Tyr294 and Lys298 provide significant acceleration effects by orientating substrates and stabilizing the negative charge developing at the substrate hydroxy oxygen. The barrier was determined to be 12.0 kcal mol^-1 , which reveals enhanced racemase activity relative to the LarA reaction using NAD⁺ -like cofactors. Compared with NAD⁺ , Ni-PTTMN has a stronger hydride-addition reactivity in moderate and high environmental polarity and may fit perfectly the moderately polar active site of LarA.

Biochemistry of Methyl-Coenzyme M Reductase

Methanogens are masters of CO2 reduction. They conserve energy by coupling H2 oxidation to the reduction of CO2 to CH4, the primary constituent of natural gas. They also generate methane by the reduction of acetic acid, methanol, methane thiol, and methylamines. Methanogens produce 109 tons of methane per year and are the major source of the earth’s atmospheric methane. Reverse methanogenesis or anaerobic methane oxidation, which is catalyzed by methanotrophic archaea living in consortia among bacteria that can act as an electron acceptor, is responsible for annual oxidation of 108 tons of methane to CO2. This chapter briefly describes the overall process of methanogenesis and then describes the enzymatic mechanism of the nickel enzyme, methyl-CoM reductase (MCR), the key enzyme in methane synthesis and oxidation. MCR catalyzes the formation of methane and the heterodisulfide (CoBSSCoM) from methyl-coenzyme M (methyl-CoM) and coenzyme B (HSCoB). Uncovering the mechanistic and molecular details of MCR catalysis is critical since methane is an abundant and important fuel and is the second (to CO2) most prevalent greenhouse gas.

The radical mechanism of biological methane synthesis by methyl-coenzyme M reductase

Methyl-coenzyme M reductase, the rate-limiting enzyme in methanogenesis and anaerobic methane oxidation, is responsible for the biological production of more than 1 billion tons of methane per year. The mechanism of methane synthesis is thought to involve either methyl-nickel(III) or methyl radical/Ni(II)-thiolate intermediates. We employed transient kinetic, spectroscopic, and computational approaches to study the reaction between the active Ni(I) enzyme and substrates. Consistent with the methyl radical-based mechanism, there was no evidence for a methyl-Ni(III) species; furthermore, magnetic circular dichroism spectroscopy identified the Ni(II)-thiolate intermediate. Temperature-dependent transient kinetics also closely matched density functional theory predictions of the methyl radical mechanism. Identifying the key intermediate in methanogenesis provides fundamental insights to develop better catalysts for producing and activating an important fuel and potent greenhouse gas.

My Lifelong Passion for Biochemistry and Anaerobic Microorganisms

Early parental influence led me first to medical school, but after developing a passion for biochemistry and sensing the need for a deeper foundation, I changed to chemistry. During breaks between semesters, I worked in various biochemistry labs to acquire a feeling for the different areas of investigation. The scientific puzzle that fascinated me most was the metabolism of the anaerobic bacterium Clostridium kluyveri, which I took on in 1965 in Karl Decker's lab in Freiburg, Germany. I quickly realized that little was known about the biochemistry of strict anaerobes such as clostridia, methanogens, acetogens, and sulfate-reducing bacteria and that these were ideal model organisms to study fundamental questions of energy conservation, CO2 fixation, and the evolution of metabolic pathways. My passion for anaerobes was born then and is unabated even after 50 years of study.

An inhibitor persistently decreased enteric methane emission from dairy cows with no negative effect on milk production

A quarter of all anthropogenic methane emissions in the United States are from enteric fermentation, primarily from ruminant livestock. This study was undertaken to test the effect of a methane inhibitor, 3-nitrooxypropanol (3NOP), on enteric methane emission in lactating Holstein cows. An experiment was conducted using 48 cows in a randomized block design with a 2-wk covariate period and a 12-wk data collection period. Feed intake, milk production, and fiber digestibility were not affected by the inhibitor. Milk protein and lactose yields were increased by 3NOP. Rumen methane emission was linearly decreased by 3NOP, averaging about 30% lower than the control. Methane emission per unit of feed dry matter intake or per unit of energy-corrected milk were also about 30% less for the 3NOP-treated cows. On average, the body weight gain of 3NOP-treated cows was 80% greater than control cows during the 12-wk experiment. The experiment demonstrated that the methane inhibitor 3NOP, applied at 40 to 80 mg/kg feed dry matter, decreased methane emissions from high-producing dairy cows by 30% and increased body weight gain without negatively affecting feed intake or milk production and composition. The inhibitory effect persisted over 12 wk of treatment, thus offering an effective methane mitigation practice for the livestock industries.

Methane oxidation by anaerobic archaea for conversion to liquid fuels

Given the recent increases in natural gas reserves and associated drawbacks of current gas-to-liquids technologies, the development of a bioconversion process to directly convert methane to liquid fuels would generate considerable industrial interest. Several clades of anaerobic methanotrophic archaea (ANME) are capable of performing anaerobic oxidation of methane (AOM). AOM carried out by ANME offers carbon efficiency advantages over aerobic oxidation by conserving the entire carbon flux without losing one out of three carbon atoms to carbon dioxide. This review highlights the recent advances in understanding the key enzymes involved in AOM (i.e., methyl-coenzyme M reductase), the ecological niches of a number of ANME, the putative metabolic pathways for AOM, and the syntrophic consortia that they typically form.