Biological and Catalytic Properties of Selenoproteins

Selenocysteine is a catalytic residue at the active site of all selenoenzymes in bacteria and mammals, and it is incorporated into the polypeptide backbone by a co-translational process that relies on the recoding of a UGA termination codon into a serine/selenocysteine codon. The best-characterized selenoproteins from mammalian species and bacteria are discussed with emphasis on their biological function and catalytic mechanisms. A total of 25 genes coding for selenoproteins have been identified in the genome of mammals. Unlike the selenoenzymes of anaerobic bacteria, most mammalian selenoenzymes work as antioxidants and as redox regulators of cell metabolism and functions. Selenoprotein P contains several selenocysteine residues and serves as a selenocysteine reservoir for other selenoproteins in mammals. Although extensively studied, glutathione peroxidases are incompletely understood in terms of local and time-dependent distribution, and regulatory functions. Selenoenzymes take advantage of the nucleophilic reactivity of the selenolate form of selenocysteine. It is used with peroxides and their by-products such as disulfides and sulfoxides, but also with iodine in iodinated phenolic substrates. This results in the formation of Se-X bonds (X = O, S, N, or I) from which a selenenylsulfide intermediate is invariably produced. The initial selenolate group is then recycled by thiol addition. In bacterial glycine reductase and D-proline reductase, an unusual catalytic rupture of selenium-carbon bonds is observed. The exchange of selenium for sulfur in selenoproteins, and information obtained from model reactions, suggest that a generic advantage of selenium compared with sulfur relies on faster kinetics and better reversibility of its oxidation reactions.

Anion Binding and Oxidative Modification at the Molybdenum Cofactor of Formate Dehydrogenase from Rhodobacter capsulatus Studied by X-ray Absorption Spectroscopy

Formate dehydrogenase (FDH) enzymes are versatile catalysts for CO₂ conversion. The FDH from Rhodobacter capsulatus contains a molybdenum cofactor with the dithiolene functions of two pyranopterin guanine dinucleotide molecules, a conserved cysteine, and a sulfido group bound at Mo(VI). In this study, we focused on metal oxidation state and coordination changes in response to exposure to O₂, inhibitory anions, and redox agents using X-ray absorption spectroscopy (XAS) at the Mo K-edge. Differences in the oxidative modification of the bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor relative to samples prepared aerobically without inhibitor, such as variations in the relative numbers of sulfido (Mo═S) and oxo (Mo═O) bonds, were observed in the presence of azide (N₃^-) or cyanate (OCN^-). Azide provided best protection against O₂, resulting in a quantitatively sulfurated cofactor with a displaced cysteine ligand and optimized formate oxidation activity. Replacement of the cysteine ligand by a formate (HCO₂^-) ligand at the molybdenum in active enzyme is compatible with our XAS data. Cyanide (CN^-) inactivated the enzyme by replacing the sulfido ligand at Mo(VI) with an oxo ligand. Evidence that the sulfido group may become protonated upon molybdenum reduction was obtained. Our results emphasize the role of coordination flexibility at the molybdenum center during inhibitory and catalytic processes of FDH enzymes.

Reaction mechanism of formate dehydrogenase studied by computational methods

Formate dehydrogenases (FDHs) are metalloenzymes that catalyse the reversible conversion of formate to carbon dioxide. Since such a process may be used to combat the greenhouse effect, FDHs have been extensively studied by experimental and theoretical methods. However, the reaction mechanism is still not clear; instead five putative mechanisms have been suggested. In this work, the reaction mechanism of FDH was studied by computational methods. Combined quantum mechanical and molecular mechanic (QM/MM) optimisations were performed to obtain the geometries. To get more accurate energies and obtain a detailed account of the surroundings, big-QM calculations with a very large (1121 atoms) QM region were performed. Our results indicate that the formate substrate does not coordinate directly to Mo when it enters the oxidised active site of the FDH, but instead resides in the second coordination sphere. The sulfido ligand abstracts a hydride ion from the substrate, giving a Mo(IV)-SH state and a thiocarbonate ion attached to Cys196. The latter releases CO₂ when the active site is oxidised back to the resting (Mo^VI) state. This mechanism is supported by recent experimental studies.

Discovery of a new metal and NAD+-dependent formate dehydrogenase from Clostridium ljungdahlii

Over the next decades, with the growing concern of rising atmospheric carbon dioxide (CO₂) levels, the importance of investigating new approaches for its reduction becomes crucial. Reclamation of CO₂ for conversion into biofuels represents an alternative and attractive production method that has been studied in recent years, now with enzymatic methods gaining more attention. Formate dehydrogenases (FDHs) are NAD(P)H-dependent oxidoreductases that catalyze the conversion of formate into CO₂ and have been extensively used for cofactor recycling in chemoenzymatic processes. A new FDH from Clostridium ljungdahlii (ClFDH) has been recently shown to possess activity in the reverse reaction: the mineralization of CO₂ into formate. In this study, we show the successful homologous expression of ClFDH in Escherichia coli. Biochemical and kinetic characterization of the enzyme revealed that this homologue also demonstrates activity toward CO₂ reduction. Structural analysis of the enzyme through homology modeling is also presented.

Anaerobic Formate and Hydrogen Metabolism

Numerous recent developments in the biochemistry, molecular biology, and physiology of formate and H2 metabolism and of the [NiFe]-hydrogenase (Hyd) cofactor biosynthetic machinery are highlighted. Formate export and import by the aquaporin-like pentameric formate channel FocA is governed by interaction with pyruvate formate-lyase, the enzyme that generates formate. Formate is disproportionated by the reversible formate hydrogenlyase (FHL) complex, which has been isolated, allowing biochemical dissection of evolutionary parallels with complex I of the respiratory chain. A recently identified sulfido-ligand attached to Mo in the active site of formate dehydrogenases led to the proposal of a modified catalytic mechanism. Structural analysis of the homologous, H2-oxidizing Hyd-1 and Hyd-5 identified a novel proximal [4Fe-3S] cluster in the small subunit involved in conferring oxygen tolerance to the enzymes. Synthesis of Salmonella Typhimurium Hyd-5 occurs aerobically, which is novel for an enterobacterial Hyd. The O2-sensitive Hyd-2 enzyme has been shown to be reversible: it presumably acts as a conformational proton pump in the H2-oxidizing mode and is capable of coupling reverse electron transport to drive H2 release. The structural characterization of all the Hyp maturation proteins has given new impulse to studies on the biosynthesis of the Fe(CN)2CO moiety of the [NiFe] cofactor. It is synthesized on a Hyp-scaffold complex, mainly comprising HypC and HypD, before insertion into the apo-large subunit. Finally, clear evidence now exists indicating that Escherichia coli can mature Hyd enzymes differentially, depending on metal ion availability and the prevailing metabolic state. Notably, Hyd-3 of the FHL complex takes precedence over the H2-oxidizing enzymes.

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.

Reduction of Carbon Dioxide by a Molybdenum-Containing Formate Dehydrogenase: A Kinetic and Mechanistic Study

Carbon dioxide accumulation is a major concern for the ecosystems, but its abundance and low cost make it an interesting source for the production of chemical feedstocks and fuels. However, the thermodynamic and kinetic stability of the carbon dioxide molecule makes its activation a challenging task. Studying the chemistry used by nature to functionalize carbon dioxide should be helpful for the development of new efficient (bio)catalysts for atmospheric carbon dioxide utilization. In this work, the ability of Desulfovibrio desulfuricans formate dehydrogenase (Dd FDH) to reduce carbon dioxide was kinetically and mechanistically characterized. The Dd FDH is suggested to be purified in an inactive form that has to be activated through a reduction-dependent mechanism. A kinetic model of a hysteretic enzyme is proposed to interpret and predict the progress curves of the Dd FDH-catalyzed reactions (initial lag phase and subsequent faster phase). Once activated, Dd FDH is able to efficiently catalyze, not only the formate oxidation (kcat of 543 s(-1), Km of 57.1 μM), but also the carbon dioxide reduction (kcat of 46.6 s(-1), Km of 15.7 μM), in an overall reaction that is thermodynamically and kinetically reversible. Noteworthy, both Dd FDH-catalyzed formate oxidation and carbon dioxide reduction are completely inactivated by cyanide. Current FDH reaction mechanistic proposals are discussed and a different mechanism is here suggested: formate oxidation and carbon dioxide reduction are proposed to proceed through hydride transfer and the sulfo group of the oxidized and reduced molybdenum center, Mo(6+)═S and Mo(4+)-SH, are suggested to be the direct hydride acceptor and donor, respectively.

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.

Sulfido and cysteine ligation changes at the molybdenum cofactor during substrate conversion by formate dehydrogenase (FDH) from Rhodobacter capsulatus

Formate dehydrogenase (FDH) enzymes are attractive catalysts for potential carbon dioxide conversion applications. The FDH from Rhodobacter capsulatus (RcFDH) binds a bis-molybdopterin-guanine-dinucleotide (bis-MGD) cofactor, facilitating reversible formate (HCOO(-)) to CO2 oxidation. We characterized the molecular structure of the active site of wildtype RcFDH and protein variants using X-ray absorption spectroscopy (XAS) at the Mo K-edge. This approach has revealed concomitant binding of a sulfido ligand (Mo=S) and a conserved cysteine residue (S(Cys386)) to Mo(VI) in the active oxidized molybdenum cofactor (Moco), retention of such a coordination motif at Mo(V) in a chemically reduced enzyme, and replacement of only the S(Cys386) ligand by an oxygen of formate upon Mo(IV) formation. The lack of a Mo=S bond in RcFDH expressed in the absence of FdsC implies specific metal sulfuration by this bis-MGD binding chaperone. This process still functioned in the Cys386Ser variant, showing no Mo-S(Cys386) ligand, but retaining a Mo=S bond. The C386S variant and the protein expressed without FdsC were inactive in formate oxidation, supporting that both Mo-ligands are essential for catalysis. Low-pH inhibition of RcFDH was attributed to protonation at the conserved His387, supported by the enhanced activity of the His387Met variant at low pH, whereas inactive cofactor species showed sulfido-to-oxo group exchange at the Mo ion. Our results support that the sulfido and S(Cys386) ligands at Mo and a hydrogen-bonded network including His387 are crucial for positioning, deprotonation, and oxidation of formate during the reaction cycle of RcFDH.

Assembly and catalysis of molybdenum or tungsten-containing formate dehydrogenases from bacteria

The global carbon cycle depends on the biological transformations of C1 compounds, which include the reductive incorporation of CO₂into organic molecules (e.g. in photosynthesis and other autotrophic pathways), in addition to the production of CO₂from formate, a reaction that is catalyzed by formate dehydrogenases (FDHs). FDHs catalyze, in general, the oxidation of formate to CO₂and H⁺. However, selected enzymes were identified to act as CO₂reductases, which are able to reduce CO₂to formate under physiological conditions. This reaction is of interest for the generation of formate as a convenient storage form of H₂for future applications. Cofactor-containing FDHs are found in anaerobic bacteria and archaea, in addition to facultative anaerobic or aerobic bacteria. These enzymes are highly diverse and employ different cofactors such as the molybdenum cofactor (Moco), FeS clusters and flavins, or cytochromes. Some enzymes include tungsten (W) in place of molybdenum (Mo) at the active site. For catalytic activity, a selenocysteine (SeCys) or cysteine (Cys) ligand at the Mo atom in the active site is essential for the reaction. This review will focus on the characterization of Mo- and W-containing FDHs from bacteria, their active site structure, subunit compositions and its proposed catalytic mechanism. We will give an overview on the different mechanisms of substrate conversion available so far, in addition to providing an outlook on bio-applications of FDHs. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications.

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.

Quantum Chemical Modeling of Enzymatic Reactions: The Case of Decarboxylation

We present a systematic study of the decarboxylation step of the enzyme aspartate decarboxylase with the purpose of assessing the quantum chemical cluster approach for modeling this important class of decarboxylase enzymes. Active site models ranging in size from 27 to 220 atoms are designed, and the barrier and reaction energy of this step are evaluated. To model the enzyme surrounding, homogeneous polarizable medium techniques are used with several dielectric constants. The main conclusion is that when the active site model reaches a certain size, the solvation effects from the surroundings saturate. Similar results have previously been obtained from systematic studies of other classes of enzymes, suggesting that they are of a quite general nature.

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.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p).

Genome sequence of Azotobacter vinelandii, an obligate aerobe specialized to support diverse anaerobic metabolic processes

Azotobacter vinelandii is a soil bacterium related to the Pseudomonas genus that fixes nitrogen under aerobic conditions while simultaneously protecting nitrogenase from oxygen damage. In response to carbon availability, this organism undergoes a simple differentiation process to form cysts that are resistant to drought and other physical and chemical agents. Here we report the complete genome sequence of A. vinelandii DJ, which has a single circular genome of 5,365,318 bp. In order to reconcile an obligate aerobic lifestyle with exquisitely oxygen-sensitive processes, A. vinelandii is specialized in terms of its complement of respiratory proteins. It is able to produce alginate, a polymer that further protects the organism from excess exogenous oxygen, and it has multiple duplications of alginate modification genes, which may alter alginate composition in response to oxygen availability. The genome analysis identified the chromosomal locations of the genes coding for the three known oxygen-sensitive nitrogenases, as well as genes coding for other oxygen-sensitive enzymes, such as carbon monoxide dehydrogenase and formate dehydrogenase. These findings offer new prospects for the wider application of A. vinelandii as a host for the production and characterization of oxygen-sensitive proteins.