Mechanism of Nitrate Reduction byDesulfovibrio desulfuricans Nitrate Reductase—A Theoretical Investigation

Microbial Denitrification: Active Site and Reaction Path Models Predict New Isotopic Fingerprints

The study of isotopic fingerprints in nitrate (δ15N, δ18O, Δ17O) has enabled pivotal insights into the global nitrogen cycle and revealed new knowledge gaps. Measuring populations of isotopic homologs of intact NO3 - ions (isotopologues) shows promise to advance the understanding of nitrogen cycling processes; however, we need new theory and predictions to guide laboratory experiments and field studies. We investigated the hypothesis that the isotopic composition of the residual nitrate pool is controlled by the N-O bond-breaking step in Nar dissimilatory nitrate reductase using molecular models of the enzyme active sites and associated kinetic isotope effects (KIEs). We integrated the molecular model results into reaction path models representing the reduction of nitrate under either closed-system or steady-state conditions. The predicted intrinsic KIE (15ε and 18ε) of the Nar active site matches observed fractionations in both culture and environmental studies. This is what would be expected if the isotopic composition of marine nitrate were controlled by dissimilatory nitrate reduction by Nar. For a closed system, the molecular models predict a pronounced negative 15N-18O clumping anomaly in residual nitrate. This signal could encode information about the amount of nitrate consumed in a closed system and thus constrain initial nitrate concentration and its isotopic composition. Similar clumped isotope anomalies can potentially be used to distinguish whether a system is open or closed to new nitrate addition. These mechanistic predictions can be tested and refined in combination with emerging ESI-Orbitrap measurements.

Mechanism and Inhibitor Exploration with Binuclear Mg Ketol-Acid Reductoisomerase: Targeting the Biosynthetic Pathway of Branched-Chain Amino Acids

Binuclear Mg ketol-acid reductoisomerase (KARI), which converts (S)-2-acetolactate into (R)-2,3-dihydroxyisovalerate, is responsible for the second step of the biosynthesis of branched-chain amino acids in plants and microorganisms and thus serves as a key inhibition target potentially without effects on mammals. Here, through the use of density functional calculations and a chemical model, the KARI-catalyzed reaction has been demonstrated to include the initial deprotonation of the substrate C2 hydroxy group, bridged by the two Mg ions, alkyl migration from the C2-alkoxide carbon atom to the C3-carbonyl carbon atom, and hydride transfer from a nicotinamide adenine dinucleotide phosphate [NAD(P)H] cofactor to C2. A dead-end mechanism with a hydride transferred to the C3 carbonyl group has been ruled out. The nucleophilicity (migratory aptitude) of the migrating carbon atom and the provision of additional negative charge to the di-Mg coordination sphere have significant effects on the steps of alkyl migration and hydride transfer, respectively. Other important mechanistic characteristics are also revealed. Inspired by the mechanism, an inhibitor (2-carboxylate-lactic acid) was designed and predicted by barrier analysis to be effective in inactivating KARI, hence probably enriching the antifungal and antibacterial library. Two types of slow substrate analogues (2-trihalomethyl acetolactic acids and 2-glutaryl lactic acid) were also found.

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.

Mechanistic investigation of trimethylamine-N-oxide reduction catalysed by biomimetic molybdenum enzyme models

In this paper, we report a theoretical investigation of the reduction reaction mechanism of Me3NO using molybdenum containing systems that are functional and structural analogues of trimethylamine N-oxide reductase mononuclear molybdenum enzyme. The reactivity of the monooxomolybdenum(IV) benzenedithiolato complex and its derivatives to carbamoyl (t-BuNHCO) and acylamino (t-BuCONH) substituents on the benzene rings in both cis and trans arrangements was explored. The calculated energy profiles describing the steps of two mechanisms of attack considered viable (named cis- and trans-attack) by the Me3NO substrate at cis and trans positions with respect to the oxo ligand show that the attack on cis is energetically more favourable than the attack on trans. Along the pathway for the cis-attack the first step of the reaction, that is rate-determining for all the studied compounds, is the approach of the substrate to the Mo centre in cis to the oxo ligand that causes a distortion of the initial square-pyramidal geometry of the complex. The reaction steps involved in the trans position attack were also explored. Calculations confirm that, as previously suggested, the introduction of ligands able to form intramolecular NH···S hydrogen bonds accelerates the reduction of the Me3NO substrate and contributes to the tuning of the reactivity of molybdoenzyme models.

The Decarboxylation of α,β-Unsaturated Acid Catalyzed by Prenylated FMN-Dependent Ferulic Acid Decarboxylase and the Enzyme Inhibition

Ferulic acid decarboxylase (Fdc1) is able to catalyze the decarboxylation of α,β-unsaturated acids using a novel cofactor, prenylated flavin mononucleotide (PrFMN). Using density functional theory calculations, we here have investigated the Fdc1 reaction mechanism with the substrate of α-methylcinnamic acid. It is demonstrated that Fdc1 employs a 1,3-dipolar cycloaddition mechanism involving four concerted steps, where the Glu282 acts as a crucial proton donor to protonate the α carbon (C_α). The last step, the decomposition of a pyrrolidine species, is rate-limiting with an overall barrier of 18.9 kcal mol^-1. Furthermore, when α-hydroxycinnamic acid is used, the Glu282 is found to have another face to transport the hydroxyl proton to the C_β atom to promote the tautomerization from enol intermediate to ketone species leading to the inhibition of the Fdc1 enzyme. The PrFMN roles are also discussed in detail.

Equilibrium Geometries of Noncovalently Bound Intermolecular Complexes Derived from Subsystem Formulation of Density Functional Theory

The subsystem formulation of density functional theory is used to obtain equilibrium geometries and interaction energies for a representative set of noncovalently bound intermolecular complexes. The results are compared with literature benchmark data. The range of applicability of two considered approximations to the exchange-correlation- and nonadditive kinetic energy components of the total energy is determined. Local density approximation, which does not involve any empirical parameters, leads to excellent intermolecular equilibrium distances for hydrogen-bonded complexes (maximal error 0.13 Å for NH3-NH3). It is a method of choice for a wide class of weak intermolecular complexes including also dipole-bound and the ones formed by rare gas atoms or saturated hydrocarbons. The range of applicability of the chosen generalized gradient approximation, which was shown in our previous works to lead to good interaction energies in such complexes, where π-electrons are involved in the interaction, remains limited to this group because it improves neither binding energies nor equilibrium geometries in the wide class of complexes for which local density approximation is adequate. An efficient energy minimization procedure, in which optimization of the geometry and the electron density of each subsystem is made simultaneously, is proposed and tested.

Periplasmic nitrate reductase and formate dehydrogenase: similar molecular architectures with very different enzymatic activities

It is remarkable how nature has been able to construct enzymes that, despite sharing many similarities, have simple but key differences that tune them for completely different functions in living cells. Periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) from the DMSOr family are representative examples of this. Both enzymes share almost identical three-dimensional protein foldings and active sites, in terms of coordination number, geometry and nature of the ligands. The substrates of both enzymes (nitrate and formate) are polyatomic anions that also share similar charge and stereochemistry. In terms of the catalytic mechanism, both enzymes have a common activation mechanism (the sulfur-shift mechanism) that ensures a constant coordination number around the metal ion during the catalytic cycle. In spite of these similarities, they catalyze very different reactions: Nap abstracts an oxygen atom from nitrate releasing nitrite, whereas FdH catalyzes a hydrogen atom transfer from formate and releases carbon dioxide. In this Account, a critical analysis of structure, function, and catalytic mechanism of the molybdenum enzymes periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) is presented. We conclude that the main structural driving force that dictates the type of reaction, catalyzed by each enzyme, is a key difference on one active site residue that is located in the top region of the active sites of both enzymes. In both enzymes, the active site is centered on the metal ion of the cofactor (Mo in Nap and Mo or W in Fdh) that is coordinated by four sulfur atoms from two pyranopterin guanosine dinucleotide (PGD) molecules and by a sulfido. However, while in Nap there is a Cys directly coordinated to the Mo ion, in FdH there is a SeCys instead. In Fdh there is also an important His that interacts very closely with the SeCys, whereas in Nap the same position is occupied by a Met. The role of Cys in Nap and SeCys in FdH is similar in both enzymes; however, Met and His have different roles. His participates directly on catalysis, and it is therefore detrimental for the catalytic cycle of FdH. Met only participates in substrate binding. We concluded that this small but key difference dictates the type of reaction that is catalyzed by each enzyme. In addition, it allows explaining why formate can bind in the Nap active site in the same way as the natural substrate (nitrate), but the reaction becomes stalled afterward.

A Bioinspired Molybdenum Complex as a Catalyst for the Photo- and Electroreduction of Protons

A molybdenum-dithiolene-oxo complex was prepared as a model of some active sites of Mo/W-dependent enzymes. The ligand, a quinoxaline-pyran-fused dithiolene, mimics molybdopterin present in these active sites. For the first time, this type of complex was shown to be active as a catalyst for the photoreduction of protons with excellent turnover numbers (500) and good stability in aqueous/organic media and for the electroreduction of protons in acetonitrile with remarkable rate constants (1030 s(-1) at -1.3 V versus Ag/AgCl). DFT calculations provided insight into the catalytic cycle of the reaction, suggesting that the oxo ligand plays a key role in proton exchange. These results provide a basis to optimize this new class of H2 -evolving catalysts.

DFT study on the oxygen transfer mechanism in nitroethenediamine based H2-receptor antagonists using the bis-dithiolene complex as the model catalyst for N-oxide reductase enzyme

Nitroethenediamine is an important functional unit, which is present in H2-receptor antagonists. These drugs show low bioavailability due to the bacterial degradation caused by the N-oxide reductase type of enzymes present in the human colon. Quantum chemical studies have been carried out to elucidate the mechanism of metabolic degradation of nitroethenediamine in the active site of N-oxide reductase. Three different pathways have been explored for the N-oxide bond cleavage by the model system, Mo(IV) bis-dithiolene complex [Mo(OMe)(mdt)2](-), (where mdt=1,2-dimethyl-ethene-1,2-dithiolate) using B3LYP/6-311+G(d,p) and M06/6-311+G(d,p) Density Functional Theory methods. The oxygen atom transfer from the nitrogen atom of nitroethenediamine to the Mo(IV) complex, involves simultaneous weakening of the N-oxide bond and the formation of Mo-O bond through a least motion path. During this transfer, Mo center is converted from a square pyramidal geometry to a distorted octahedral geometry, to facilitate the process of oxygen atom transfer. The energy barrier for the oxygen atom transfer from the imine tautomer has been estimated to be 25.9kcal/mol however, the overall reaction has been found to be endothermic. On the other hand, oxygen transfer reaction from the nitronic acid tautomer requires 30.5kcal/mol energy leading to a highly exothermic metabolite (M-1) directly hence, this path can be considered thermodynamically favorable for this metabolite. The alternative path involving the oxygen atom transfer from the enamine tautomer requires comparatively a higher energy barrier (32.6kcal/mol) and leads to a slightly endothermic metabolite. This study established the structural and energetic details associated with the Mo(IV) bis-dithiolene complex that catalyzes the degradation of nitroethenediamine based drug molecules.

Proposed toxic and hypoxic impairment of a brainstem locus in autism

Electrophysiological findings implicate site-specific impairment of the nucleus tractus solitarius (NTS) in autism. This invites hypothetical consideration of a large role for this small brainstem structure as the basis for seemingly disjointed behavioral and somatic features of autism. The NTS is the brain's point of entry for visceral afference, its relay for vagal reflexes, and its integration center for autonomic control of circulatory, immunological, gastrointestinal, and laryngeal function. The NTS facilitates normal cerebrovascular perfusion, and is the seminal point for an ascending noradrenergic system that modulates many complex behaviors. Microvascular configuration predisposes the NTS to focal hypoxia. A subregion--the "pNTS"--permits exposure to all blood-borne neurotoxins, including those that do not readily transit the blood-brain barrier. Impairment of acetylcholinesterase (mercury and cadmium cations, nitrates/nitrites, organophosphates, monosodium glutamate), competition for hemoglobin (carbon monoxide, nitrates/nitrites), and higher blood viscosity (net systemic oxidative stress) are suggested to potentiate microcirculatory insufficiency of the NTS, and thus autism.

The sulfur shift: an activation mechanism for periplasmic nitrate reductase and formate dehydrogenase

A structural rearrangement known as sulfur shift occurs in some Mo-containing enzymes of the DMSO reductase family. This mechanism is characterized by the displacement of a coordinating cysteine thiol (or SeCys in Fdh) from the first to the second shell of the Mo-coordination sphere metal. The hexa-coordinated Mo ion found in the as-isolated state cannot bind directly any exogenous ligand (substrate or inhibitors), while the penta-coordinated ion, attained upon sulfur shift, has a free binding site for direct coordination of the substrate. This rearrangement provides an efficient mechanism to keep a constant coordination number throughout an entire catalytic pathway. This mechanism is very similar to the carboxylate shift observed in Zn-dependent enzymes, and it has been recently detected by experimental means. In the present paper, we calculated the geometries and energies involved in the sulfur-shift mechanism using QM-methods (M06/(6-311++G(3df,2pd),SDD)//B3LYP/(6-31G(d),SDD)). The results indicated that the sulfur-shift mechanism provides an efficient way to enable the metal ion for substrate coordination.

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.

The prokaryotic Mo/W-bisPGD enzymes family: a catalytic workhorse in bioenergetic

Over the past two decades, prominent importance of molybdenum-containing enzymes in prokaryotes has been put forward by studies originating from different fields. Proteomic or bioinformatic studies underpinned that the list of molybdenum-containing enzymes is far from being complete with to date, more than fifty different enzymes involved in the biogeochemical nitrogen, carbon and sulfur cycles. In particular, the vast majority of prokaryotic molybdenum-containing enzymes belong to the so-called dimethylsulfoxide reductase family. Despite its extraordinary diversity, this family is characterized by the presence of a Mo/W-bis(pyranopterin guanosine dinucleotide) cofactor at the active site. This review highlights what has been learned about the properties of the catalytic site, the modular variation of the structural organization of these enzymes, and their interplay with the isoprenoid quinones. In the last part, this review provides an integrated view of how these enzymes contribute to the bioenergetics of prokaryotes. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Mechanism of tungsten-dependent acetylene hydratase from quantum chemical calculations

Acetylene hydratase is a tungsten-dependent enzyme that catalyzes the nonredox hydration of acetylene to acetaldehyde. Density functional theory calculations are used to elucidate the reaction mechanism of this enzyme with a large model of the active site devised on the basis of the native X-ray crystal structure. Based on the calculations, we propose a new mechanism in which the acetylene substrate first displaces the W-coordinated water molecule, and then undergoes a nucleophilic attack by the water molecule assisted by an ionized Asp13 residue at the active site. This is followed by proton transfer from Asp13 to the newly formed vinyl anion intermediate. In the subsequent isomerization, Asp13 shuttles a proton from the hydroxyl group of the vinyl alcohol to the α-carbon. Asp13 is thus a key player in the mechanism, but also W is directly involved in the reaction by binding and activating acetylene and providing electrostatic stabilization to the transition states and intermediates. Several other mechanisms are also considered but the energetic barriers are found to be very high, ruling out these possibilities.

On the inhibitor effects of bergamot juice flavonoids binding to the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) enzyme

Density functional theory was applied to study the binding mode of new flavonoids as possible inhibitors of the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), an enzyme that catalyzes the four-electron reduction of HMGCoA to mevalonate, the committed step in the biosynthesis of sterols. The investigated flavonoid conjugates brutieridin and melitidin were recently quantified in the bergamot fruit extracts and identified to be structural analogues of statins, lipids concentration lowering drugs that inhibit HMGR. Computations allowed us to perform a detailed analysis of the geometrical and electronic features affecting the binding of these compounds, as well as that of the excellent simvastatin drug, to the active site of the enzyme and to give better insight into the inhibition process.

Human insulin-degrading enzyme working mechanism

The possible mechanism by which the insulin-degrading enzyme (IDE) zinc-binding protease carries out its catalytic function toward two peptides of different length, simulating a portion of B chain of insulin, was investigated on an enzymatic model consisting of 130 /159 atoms, using the density functional theory method and the hybrid exchange-correlation functional B3LYP in gas phase and in the protein environment. Based on the geometry and relative stabilities of minima and transition states on the potential energy profiles, we determined that proteolysis reaction is exothermic and proceeds quickly as the barrier in the rate-limiting step falls widely within the range of values expected for an enzymatic catalysis, both in vacuum and in protein medium.