Tungsten-dependent formaldehyde ferredoxin oxidoreductase: Reaction mechanism from quantum chemical calculations

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.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

A bacterial tungsten-containing aldehyde oxidoreductase forms an enzymatic decorated protein nanowire

Aldehyde oxidoreductases (AORs) are tungsten enzymes catalyzing the oxidation of many different aldehydes to the corresponding carboxylic acids. In contrast to other known AORs, the enzyme from the denitrifying betaproteobacterium Aromatoleum aromaticum (AORAa) consists of three different subunits (AorABC) and uses nicotinamide adenine dinucleotide (NAD) as an electron acceptor. Here, we reveal that the enzyme forms filaments of repeating AorAB protomers that are capped by a single NAD-binding AorC subunit, based on solving its structure via cryo-electron microscopy. The polyferredoxin-like subunit AorA oligomerizes to an electron-conducting nanowire that is decorated with enzymatically active and W-cofactor (W-co) containing AorB subunits. Our structure further reveals the binding mode of the native substrate benzoate in the AorB active site. This, together with quantum mechanics:molecular mechanics (QM:MM)-based modeling for the coordination of the W-co, enables formulation of a hypothetical catalytic mechanism that paves the way to further engineering for applications in synthetic biology and biotechnology.

The Development of Tungsten Biochemistry-A Personal Recollection

The development of tungsten biochemistry is sketched from the viewpoint of personal participation. Following its identification as a bio-element, a catalogue of genes, enzymes, and reactions was built up. EPR spectroscopic monitoring of redox states was, and remains, a prominent tool in attempts to understand tungstopterin-based catalysis. A paucity of pre-steady-state data remains a hindrance to overcome to this day. Tungstate transport systems have been characterized and found to be very specific for W over Mo. Additional selectivity is presented by the biosynthetic machinery for tungstopterin enzymes. Metallomics analysis of hyperthermophilic archaeon Pyrococcus furiosus indicates a comprehensive inventory of tungsten proteins.

An unprecedented function for a tungsten-containing oxidoreductase

Five tungstopterin-containing oxidoreductases were characterized from the hyperthermophile Pyrococcus furiosus. Each enzyme catalyzes the reversible conversion of one or more aldehydes to the corresponding carboxylic acid, but they have different specificities. The physiological functions of only two of these enzymes are known: one, termed GAPOR, is a glycolytic enzyme that oxidizes glyceraldehyde-3-phosphate, while the other, termed AOR, oxidizes multiple aldehydes generated during peptide fermentation. Two of the enzymes have known structures (AOR and FOR). Herein, we focus on WOR5, the fifth tungstopterin enzyme to be discovered in P. furiosus. Expression of WOR5 was previously shown to be increased during cold shock (growth at 72 ℃), although the physiological substrate is not known. To gain insight into WOR5 function, we sought to determine both its structure and identify its intracellular substrate. Crystallization experiments were performed with a concentrated cytoplasmic extract of P. furiosus grown at 72 ℃ and the structure of WOR5 was deduced from the crystals that were obtained. In contrast to a previous report, WOR5 is heterodimeric containing an additional polyferredoxin-like subunit with four [4Fe-4S] clusters. The active site structure of WOR5 is substantially different from that of AOR and FOR and the significant electron density observed adjacent to the tungsten cofactor of WOR5 was modeled as an aliphatic sulfonate. Biochemical assays and product analysis confirmed that WOR5 is an aliphatic sulfonate ferredoxin oxidoreductase (ASOR). A catalytic mechanism for ASOR is proposed based on the structural information and the potential role of ASOR in the cold-shock response is discussed.

A QM/MM investigation of the catalytic mechanism of acetylene hydratase: insights into engineering a more effective enzyme

In the present investigation, a QM/MM approach was used to better understand the effect of the second environmental shell of the active site on the catalytic conversion of acetylene to acetaldehyde by acetylene hydratase (AH). In addition, the effect of substituting W-coordinating sulfur atoms with selenium atoms was done to provide insight into the influence of the W-coordinating atoms on the catalytic reaction. From the results, it found that the presence of the second shell environment had a significant effect on the reaction. Specifically, in the absence of the MM second shell environment (i.e., QM-cluster model), the rate-determining step is defined by the first proton transfer step. In contrast, for the QM/MM model, the rate-determining step is defined by the water attacking step. Moreover, with the presence of the MM second shell environment, a key intermediate found in the DFT-cluster investigation does not exist in the QM/MM investigation. Rather, what was a two-step process in the DFT-cluster study was calculated to occur in a single step for the QM/MM study. Regarding the sulfur to selenium substitutions, it was found that Gibbs energy for the acetylene binding phase was significantly affected. Notably, the trans-position selenium made the binding of acetylene 65.6 kJ mol−1 less endergonic. Moreover, the overall reaction became 38.2 kJ mol−1 less endergonic compared with the wild type (WT) AH model. Thus, the substitution of key W-coordinating sulfur atoms with selenium atoms may offer a means to enhance the catalytic mechanism of AH considerably.

Steady-state kinetics of the tungsten containing aldehyde: ferredoxin oxidoreductases from the hyperthermophilic archaeon Pyrococcus furiosus

The tungsten containing Aldehyde:ferredoxin oxidoreductases (AOR) offer interesting opportunities for biocatalytic approaches towards aldehyde oxidation and carboxylic acid reduction. The hyperthermophilic archaeon Pyrococcus furiosus encodes five different AOR family members: glyceraldehyde-3-phosphate oxidoreductase (GAPOR), aldehyde oxidoreductase (AOR), and formaldehyde oxidoreductase (FOR), WOR4 and WOR5. GAPOR functions as a glycolytic enzyme and is highly specific for the substrate glyceraldehyde-3-phosphate (GAP). AOR, FOR and WOR5 have a broad substrate spectrum, and for WOR4 no substrate has been identified to date. As ambiguous kinetic parameters have been reported for different AOR family enzymes the steady state kinetics under different physiologically relevant conditions was explored. The GAPOR substrate GAP was found to degrade at 60 °C by non-enzymatic elimination of the phosphate group to methylglyoxal with a half-life t_1/2 = 6.5 min. Methylglyoxal is not a substrate or inhibitor of GAPOR. D-GAP was identified as the only substrate oxidized by GAPOR, and the kinetics of the enzyme was unaffected by the presence of L-GAP, which makes GAPOR the first enantioselective enzyme of the AOR family. The steady-state kinetics of GAPOR showed partial substrate inhibition, which assumes the GAP inhibited form of the enzyme retains some activity. This inhibition was found to be alleviated completely by a 1 M NaCl resulting in increased enzyme activity at high substrate concentrations. GAPOR activity was strongly pH dependent, with the optimum at pH 9. At pH 9, the substrate is a divalent anion and, therefore, positively charged amino acid residues are likely to be involved in the binding of the substrate. FOR exhibited a significant primary kinetic isotope effect of the apparent Vmax for the deuterated substrate, formaldehyde-d₂, which shows that the rate-determining step involves a CH bond break from the aldehyde. The implications of these results for the reaction mechanism of tungsten-containing AORs, are discussed.

How does Mo-dependent perchlorate reductase work in the decomposition of oxyanions?

The mechanisms of Mo-dependent perchlorate reductase (PcrAB)-catalyzed decomposition of perchlorate, bromate, iodate, and nitrate were revealed by density functional calculations.

Computational Understanding of the Selectivities in Metalloenzymes

Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.

DFT Study on the Interaction of Tris(benzene-1,2-dithiolato)molybdenum Complex with Water. A Hydrolysis Mechanism Involving a Feasible Seven-Coordinate Aquomolybdenum Intermediate

In the present work, the reactivity of the tris(benzene-1,2-dithiolato)molybdenum complex ([Mo(bdt)₃]) toward water is studied by means of the density functional theory (DFT). DFT calculations were performed using the M06, B3P86, and B3PW91 hybrid functionals for comparison purposes. The M06 method was employed to elucidate the reaction pathway, relative stability of the intermediate products, nature of the Mo-S bond cleavage, and electronic structure of the involved molybdenum species. This functional was also used to study the transference of electrons from the molybdenum center toward the ligands. The reaction pathway confirms that [Mo(bdt)₃] undergoes hydrolysis, yielding dihydroxo-bis(benzene-1,2-dithiolato)molybdenum complex ([Mo(OH)₂(bdt)₂]) and benzenedithiol. The reaction takes place through seven transition structures, one of them involving an aquo seven-coordinate molybdenum intermediate stabilized by a lone pair (LP) LP_O→LP_Mo hyperconjugative interaction. This heptacoordinate species allows understanding of the observed oxygen atom exchange between water and tertiary phosphines mediated by these complexes. Calculations also show that [Mo(C₂H₄S₂)₃] and [Mo(OH)₂(C₂H₄S₂)₂] have d² and d⁰ electronic configuration, and hence an electron pair must be transferred during the course of the hydrolysis. The frontier molecular orbital (FMO) analysis concludes that the electron pair is transferred in the rupture of the second Mo-S bond, from the occupied donating Mo d_x²-y² orbital to the unoccupied C₂H₄(SH)₂ S-C σ* ligand orbital. This result is supported by the bond dissociation energy calculations, which demonstrate that the neutral dissociation of the second Mo-S bond is energetically the more favorable.

Theoretically exploring the key role of the Lys412 residue in the conversion of N2O to N2by nitrous oxide reductase from Achromobacter cycloclastes

The active Cu_Zcluster of NOR provides strong back-donation to coordinated N₂O and activates the O atom of the N₂O group facilitating H-bonding and protonationviathe Lys412 residue.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Why is the molybdenum-substituted tungsten-dependent formaldehyde ferredoxin oxidoreductase not active? A quantum chemical study

Formaldehyde ferredoxin oxidoreductase is a tungsten-dependent enzyme that catalyzes the oxidative degradation of formaldehyde to formic acid. The molybdenum ion can be incorporated into the active site to displace the tungsten ion, but is without activity. Density functional calculations have been employed to understand the incapacitation of the enzyme caused by molybdenum substitution. The calculations show that the enzyme with molybdenum (Mo-FOR) has higher redox potential than that with tungsten, which makes the formation of the Mo(VI)=O complex endothermic by 14 kcal/mol. Following our previously suggested mechanism for this enzyme, the formaldehyde substrate oxidation was also investigated for Mo-FOR using the same quantum-mechanics-only model, except for the displacement of tungsten by molybdenum. The calculations demonstrate that formaldehyde oxidation occurs via a sequential two-step mechanism. Similarly to the tungsten-catalyzed reaction, the Mo(VI)=O species performs the nucleophilic attack on the formaldehyde carbon, followed by proton transfer in concert with two-electron reduction of the metal center. The first step is rate-limiting, with a total barrier of 28.2 kcal/mol. The higher barrier is mainly due to the large energy penalty for the formation of the Mo(VI)=O species.