Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase.

[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in the complex with the wild-type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-random accelerated molecular dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122 and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2, and H2 from the wild-type [NiFe] hydrogenase, and we observed that while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.

New light on ancient enzymes - in vitro CO2 Fixation by Pyruvate Synthase of Desulfovibrio africanus and Sulfolobus acidocaldarius.

Two variants of the enzyme family pyruvate:ferredoxin oxidoreductase (PFOR), derived from the anaerobic sulfate-reducing bacterium Desulfovibrio africanus and the extremophilic crenarchaeon Sulfolobus acidocaldarius, respectively, were evaluated for their capacity to fixate CO2 in vitro. PFOR reversibly catalyzes the conversion of acetyl-CoA and CO2 to pyruvate using ferredoxin as redox partner. The oxidative decarboxylation of pyruvate is thermodynamically strongly favored, and most previous studies only considered the oxidative direction of the enzyme. To assay the pyruvate synthase function of PFOR during reductive carboxylation of acetyl-CoA is more challenging and requires to maintain the reaction far from equilibrium. For this purpose, a biochemical assay was established where low-potential electrons were introduced by photochemical reduction of EDTA/deazaflavin and the generated pyruvate was trapped by chemical derivatization with semicarbazide. The product of CO2 fixation could be detected as pyruvate semicarbazone by HPLC-MS. In a combinatorial approach, both PFORs were tested with ferredoxins from different sources. The pyruvate semicarbazone product could be detected with low-potential ferredoxins of the green sulfur bacterium Chlorobium tepidum and of S. acidocaldarius whereas CO2 fixation was not supported by the native ferredoxin of D. africanus. Methylviologen as an artificial electron carrier also allowed CO2 fixation. For both enzymes, the results are the first demonstration of CO2 fixation in vitro. Both enzymes exhibited high stability in the presence of oxygen during purification and storage. In conclusion, the employed PFOR enzymes in combination with non-native ferredoxin cofactors might be promising candidates for further incorporation in biocatalytic CO2 conversion. ENZYMES: EC1.2.7.1. Pyruvate:Ferredoxin Oxidoreductase.

Hydrogenases and H2 metabolism in sulfate-reducing bacteria of the Desulfovibrio genus.

Hydrogen metabolism plays a central role in sulfate-reducing bacteria of the Desulfovibrio genus and is based on hydrogenases that catalyze the reversible conversion of protons into dihydrogen. These metabolically versatile microorganisms possess a complex hydrogenase system composed of several enzymes of both [FeFe]- and [NiFe]-type that can vary considerably from one Desulfovibrio species to another. This review covers the molecular and physiological aspects of hydrogenases and H2 metabolism in Desulfovibrio but focuses particularly on our model bacterium Desulfovibrio fructosovorans. The search of hydrogenase genes in more than 30 sequenced genomes provides an overview of the distribution of these enzymes in Desulfovibrio. Our discussion will consider the significance of the involvement of electron-bifurcation in H2 metabolism.

Adaptation of Desulfovibrio alaskensis G20 to perchlorate, a specific inhibitor of sulfate reduction.

Hydrogen sulfide produced by sulfate-reducing microorganisms (SRM) poses significant health and economic risks, particularly during oil recovery. Previous studies identified perchlorate as a specific inhibitor of SRM. However, constant inhibitor addition to natural systems results in new selective pressures. Consequently, we investigated the ability of Desulfovibrio alaskensis G20 to evolve perchlorate resistance. Serial transfers in increasing concentrations of perchlorate led to robust growth in the presence of 100 mM inhibitor. Isolated adapted strains demonstrated a threefold increase in perchlorate resistance compared to the wild-type ancestor. Whole genome sequencing revealed a single base substitution in Dde_2265, the sulfate adenylyltransferase (sat). We purified and biochemically characterized the Sat from both wild-type and adapted strains, and showed that the adapted Sat was approximately threefold more resistant to perchlorate inhibition, mirroring whole cell results. The ability of this mutation to confer resistance across other inhibitors of sulfidogenesis was also assayed. The generalizability of this mutation was confirmed in multiple evolving G20 cultures and in another SRM, D. vulgaris Hildenborough. This work demonstrates that a single nucleotide polymorphism in Sat can have a significant impact on developing perchlorate resistance and emphasizes the value of adaptive laboratory evolution for understanding microbial responses to environmental perturbations.

A new mechanistic model for an O2-protected electron-bifurcating hydrogenase, Hnd from Desulfovibrio fructosovorans.

The genome of the sulfate-reducing and anaerobic bacterium Desulfovibrio fructosovorans encodes different hydrogenases. Among them is Hnd, a tetrameric cytoplasmic [FeFe] hydrogenase that has previously been described as an NADP-specific enzyme (Malki et al., 1995). In this study, we purified and characterized a recombinant Strep-tagged form of Hnd and demonstrated that it is an electron-bifurcating enzyme. Flavin-based electron-bifurcation is a mechanism that couples an exergonic redox reaction to an endergonic one allowing energy conservation in anaerobic microorganisms. One of the three ferredoxins of the bacterium, that was named FdxB, was also purified and characterized. It contains a low-potential (Em = -450 mV) [4Fe4S] cluster. We found that Hnd was not able to reduce NADP+, and that it catalyzes the simultaneous reduction of FdxB and NAD+. Moreover, Hnd is the first electron-bifurcating hydrogenase that retains activity when purified aerobically due to formation of an inactive state of its catalytic site protecting against O2 damage (Hinact). Hnd is highly active with the artificial redox partner (methyl viologen) and can perform the electron-bifurcation reaction to oxidize H2 with a specific activity of 10 µmol of NADH/min/mg of enzyme. Surprisingly, the ratio between NADH and reduced FdxB varies over the reaction with a decreasing amount of FdxB reduced per NADH produced, indicating a more complex mechanism than previously described. We proposed a new mechanistic model in which the ferredoxin is recycled at the hydrogenase catalytic subunit.

A cell-free platform for rapid synthesis and testing of active oligosaccharyltransferases.

Protein glycosylation, or the attachment of sugar moieties (glycans) to proteins, is important for protein stability, activity, and immunogenicity. However, understanding the roles and regulations of site-specific glycosylation events remains a significant challenge due to several technological limitations. These limitations include a lack of available tools for biochemical characterization of enzymes involved in glycosylation. A particular challenge is the synthesis of oligosaccharyltransferases (OSTs), which catalyze the attachment of glycans to specific amino acid residues in target proteins. The difficulty arises from the fact that canonical OSTs are large (>70 kDa) and possess multiple transmembrane helices, making them difficult to overexpress in living cells. Here, we address this challenge by establishing a bacterial cell-free protein synthesis platform that enables rapid production of a variety of OSTs in their active conformations. Specifically, by using lipid nanodiscs as cellular membrane mimics, we obtained yields of up to 420 µg/ml for the single-subunit OST enzyme, "Protein glycosylation B" (PglB) from Campylobacter jejuni, as well as for three additional PglB homologs from Campylobacter coli, Campylobacter lari, and Desulfovibrio gigas. Importantly, all of these enzymes catalyzed N-glycosylation reactions in vitro with no purification or processing needed. Furthermore, we demonstrate the ability of cell-free synthesized OSTs to glycosylate multiple target proteins with varying N-glycosylation acceptor sequons. We anticipate that this broadly applicable production method will advance glycoengineering efforts by enabling preparative expression of membrane-embedded OSTs from all kingdoms of life.

Desulfovibrio DA2_CueO is a novel multicopper oxidase with cuprous, ferrous and phenol oxidase activity.

Desulfovibrio sp. A2 is a novel Gram-negative sulfate-reducing bacterium that was isolated from sediments of the Norilsk mining/smelting area in Russia. The organism possesses a monocistronic operon encoding a 71 kDa periplasmic multicopperoxidase, which we call DA2_CueO. Histidine-tagged DA2_CueO expressed from a plasmid in Escherichia coli and purified by Ni-NTA affinity chromatography oxidizes Cu+ and Fe2+, and exhibits phenol oxidase activity with 2,2-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid), 2,3-dihydroxybenzoic acid and 2,6-dimethoxyphenol as substrates, using O2 as the oxidant. When expressed in an E. coli cueO knock-out strain, DA2_CueO exhibits phenol oxidase activity in vivo and enhances the copper tolerance of the strain. These findings indicate that the DA2_CueO gene of Desulfovibrio sp. A2 encodes a multicopperoxidase with a role in metal ion resistance. The enzyme displays some novel structural features, which are discussed.

Transcription of [FeFe]-Hydrogenase Genes during H2 Production in Clostridium and Desulfovibrio spp. Isolated from a Paddy Field Soil.

Changes in the relative abundances of the transcripts of hydA gene paralogs for [FeFe]-hydrogenase in Clostridium sp. strain H2 and Desulfovibrio sp. strain A1 isolated from paddy field soil were analyzed during H2 production. Strains H2 and A1 had at least five and two phylogenetically different hydA genes, respectively. The relative abundances of their hydA transcripts differed among the paralogs and H2 production activity changed in a manner that depended on the growth phase and conditions. Increases or decreases in the relative abundances of the transcripts of two out of five hydA genes in strain H2 correlated with changes in H2 production rates, whereas those of the others remained unchanged or decreased. In strain A1, the relative abundances of the transcripts of two hydA genes differed between monoculture, sulfate-reducing, and syntrophic, methanogenic conditions. The relative abundance of the transcripts of one hydA gene, predicted to encode a cytosolic [FeFe]-hydrogenase, was higher under syntrophic, methanogenic conditions than sulfate-reducing conditions, while that of the transcripts of the other hydA gene decreased with time under both conditions. This study showed that the transcription of the hydA gene during growth with active H2 production was differently regulated among the paralogs in H2 producers isolated from paddy field soil.

Molecular Basis of C-N Bond Cleavage by the Glycyl Radical Enzyme Choline Trimethylamine-Lyase.

Deamination of choline catalyzed by the glycyl radical enzyme choline trimethylamine-lyase (CutC) has emerged as an important route for the production of trimethylamine, a microbial metabolite associated with both human disease and biological methane production. Here, we have determined five high-resolution X-ray structures of wild-type CutC and mechanistically informative mutants in the presence of choline. Within an unexpectedly polar active site, CutC orients choline through hydrogen bonding with a putative general base, and through close interactions between phenolic and carboxylate oxygen atoms of the protein scaffold and the polarized methyl groups of the trimethylammonium moiety. These structural data, along with biochemical analysis of active site mutants, support a mechanism that involves direct elimination of trimethylamine. This work broadens our understanding of radical-based enzyme catalysis and will aid in the rational design of inhibitors of bacterial trimethylamine production.

Spatial Abundance and Distribution of Potential Microbes and Functional Genes Associated with Anaerobic Mineralization of Pentachlorophenol in a Cylindrical Reactor.

Functional interplays of microbial activity, genetic diversity and contaminant transformation are poorly understood in reactors for mineralizing halogenated aromatics anaerobically. Here, we investigated abundance and distribution of potential microbes and functional genes associated with pentachlorophenol (PCP) anaerobic mineralization in a continuous-flow cylindrical reactor (15 cm in length). PCP dechlorination and the metabolite (phenol) were observed at segments 0-8 cm from inlet, where key microbes, including potential reductive dechlorinators (Dehalobacter, Sulfurospirillum, Desulfitobacterium and Desulfovibrio spp.) and phenol degraders (Cryptanaerobacter and Syntrophus spp.), as well as putative functional genes, including putative chlorophenol reductive dehalogenase (cprA) and benzoyl-CoA reductase (bamB), were highly enriched simultaneously. Five types of putative cprAs, three types of putative bamBs and seven types of putative nitrogenase reductase (nifHs) were determined, with their copy numbers decreased gradually from inlet to outlet. Distribution of chemicals, bacteria and putative genes confirmed PCP dechlorination and phenol degradation accomplished in segments 0-5 cm and 0-8 cm, respectively, contributing to a high PCP mineralization rate of 3.86 µM d(-1). Through long-term incubation, dechlorination, phenol degradation and nitrogen fixation bacteria coexisted and functioned simultaneously near inlet (0-8 cm), verified the feasibility of anaerobic mineralization of halogenated aromatics in the compact reactor containing multiple functional microbes.

Kinetic Properties of Pyruvate Ferredoxin Oxidoreductase of Intestinal Sulfate-Reducing Bacteria Desulfovibrio piger Vib-7 and Desulfomicrobium sp. Rod-9.

Intestinal sulfate-reducing bacteria reduce sulfate ions to hydrogen sulfide causing inflammatory bowel diseases of humans and animals. The bacteria consume lactate as electron donor which is oxidized to acetate via pyruvate in process of the dissimilatory sulfate reduction. Pyruvate-ferredoxin oxidoreductase activity and the kinetic properties of the enzyme from intestinal sulfate-reducing bacteria Desulfovibrio piger and Desulfomicrobium sp. have never been well-characterized and have not been yet studied. In this paper we present for the first time the specific activity of pyruvate-ferredoxin oxidoreductase and the kinetic properties of the enzyme in cell-free extracts of both D. piger Vib-7 and Desulfomicrobium sp. Rod-9 intestinal bacterial strains. Microbiological, biochemical, biophysical and statistical methods were used in this work. The optimal temperature (+35°C) and pH 8.5 for enzyme reaction were determined. The spectral analysis of the puri- fied pyruvate-ferredoxin oxidoreductase from the cell-free extracts was demonstrated. Analysis of the kinetic properties of the studied enzyme was carried out. Initial (instantaneous) reaction velocity (V0), maximum amount of the product of reaction (Pmax), the reaction time (half saturation period) and maximum velocity of the pyruvate-ferredoxin oxidoreductase reaction (V ) were defined. Michaelis constants (Km) of the enzyme reaction were calculated for both intestinal bacterial strains. The studies of the kinetic enzyme properties in the intestinal sulfate-reducing bacteria strains in detail can be prospects for clarifying the etiological role of these bacteria in the development of inflammatory bowel diseases.

Rates and Routes of Electron Transfer of [NiFe]-Hydrogenase in an Enzymatic Fuel Cell.

Hydrogenase enzymes are being used in enzymatic fuel cells immobilized on a graphite or carbon electrode surface, for example. The enzyme is used for the anodic oxidation of molecular hydrogen (H2) to produce protons and electrons. The association and orientation of the enzyme at the anode electrode for a direct electron transfer is not completely resolved. The distal FeS-cluster in [NiFe]-hydrogenases contains a histidine residue which is known to play a critical role in the intermolecular electron transfer between the enzyme and the electrode surface. The [NiFe]-hydrogenase graphite electrode association was investigated using Brownian Dynamics simulations. Residues that were shown to be in proximity to the electrode surface were identified (His184, Ser196, Glu461, Glu464), and electron transfer routes connecting the distal FeS-cluster with the surface residues were investigated. Several possible pathways for electron transfer between the distal FeS-cluster and the terminal amino acid residues were probed in terms of their rates of electron transfer using DFT methods. The reorganization energies λ of the distal iron-sulfur cluster and coronene as a molecular model for graphite were calculated. The reorganization energy of the distal (His)(Cys)3 cluster was found to be not very different from that of a standard cubane clusters with a (Cys)4 coordination. Electronic coupling matrix elements and rates of electron transfer for the different pathways were calculated according to the Marcus equation. The rates for glutamate-mediated electrode binding were found to be incompatible with experimental data. A direct electron transfer from the histidine ligand of the distal FeS-cluster to the electrode yielded rates of electron transfer in excellent agreement with experiment. A second pathway, however, from the distal FeS-cluster to the Ser196 residue was found to be equally efficient and feasible.

[NiFe]-hydrogenases revisited: nickel-carboxamido bond formation in a variant with accrued O2-tolerance and a tentative re-interpretation of Ni-SI states.

[NiFe]-hydrogenases are well-studied enzymes capable of oxidizing molecular hydrogen and reducing protons. EPR and FTIR spectroscopic studies have shown that these enzymes can be isolated in several redox states that include paramagnetic oxidized inactive Ni-A and Ni-B species and a reduced Ni-C form. The latter and the diamagnetic respectively more oxidized Ni-SI and more reduced Ni-R forms are generally thought to be involved in the catalytic cycle of [NiFe]-hydrogenases. With the exception of Ni-SI, these different stable states have been well characterized. Here, based on the crystal structure of a partially reduced Desulfovibrio fructosovorans (Df) enzyme and data from the literature we propose that at least one of the Ni-SI sub-states contains an unexpected combination of hydride and sulfenic acid moieties. We have also determined the structure of the less oxygen-sensitive Df [NiFe]-hydrogenase V74C mutant and found that more than half of the active site nickel occupies a novel position, called Ni'. In this new position, the metal ion is coordinated by two cysteine thiolates, a bridging species modeled as SH(-) and a main chain carboxamido N atom. The Ni' coordination is similar to the one found in Ni superoxide dismutase, an enzyme that operates at significantly more positive potentials than [NiFe]-hydrogenases. We propose that the oxygen-tolerance of the V74C variant results from a high potential stabilization of a Ni'(iii) species induced by the change in the metal ion coordination sphere. We also propose that transient Ni'(iii) species can rapidly attract successive electrons from the Fe4S4 proximal cluster accelerating the reduction of oxygen to water and hydroxide. The naturally occurring oxygen-tolerant [NiFe]-hydrogenases have an unusual proximal cluster that has been shown to be exceptionally plastic and capable of undergoing two successive one-electron oxidations. This double oxidation is modulated by the migration of one of the iron atoms in the cluster to the main chain where, as Fe(iii), it forms a bond with a carboxamido N ligand. Like in the Df V74C variant the electrons from the proximal cluster help reducing O2 to H2O and OH(-). In conclusion, in both cases a metal-carboxamido bond may explain, at least partially, the observed oxygen tolerance.

A threonine stabilizes the NiC and NiR catalytic intermediates of [NiFe]-hydrogenase.

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

Multiscale simulations give insight into the hydrogen in and out pathways of [NiFe]-hydrogenases from Aquifex aeolicus and Desulfovibrio fructosovorans.

[NiFe]-hydrogenases catalyze the cleavage of molecular hydrogen into protons and electrons and represent promising tools for H2-based technologies such as biofuel cells. However, many aspects of these enzymes remain to be understood, in particular how the catalytic center can be protected from irreversible inactivation by O2. In this work, we combined homology modeling, all-atom molecular dynamics, and coarse-grain Brownian dynamics simulations to investigate and compare the dynamic and mechanical properties of two [NiFe]-hydrogenases: the soluble O2-sensitive enzyme from Desulfovibrio fructosovorans, and the O2-tolerant membrane-bound hydrogenase from Aquifex aeolicus. We investigated the diffusion pathways of H2 from the enzyme surface to the central [NiFe] active site, and the possible proton pathways that are used to evacuate hydrogen after the oxidation reaction. Our results highlight common features of the two enzymes, such as a Val/Leu/Arg triad of key residues that controls ligand migration and substrate access in the vicinity of the active site, or the key role played by a Glu residue for proton transfer after hydrogen oxidation. We show specificities of each hydrogenase regarding the enzymes internal tunnel network or the proton transport pathways.

ArsC3 from Desulfovibrio alaskensis G20, a cation and sulfate-independent highly efficient arsenate reductase.

Desulfovibrio alaskensis G20, a sulfate-reducing bacterium, contains an arsRBC2C3 operon that encodes two putative arsenate reductases, DaG20_ArsC2 and DaG20_ArsC3. In this study, resistance assays in E. coli transformed with plasmids containing either of the two recombinant arsenate reductases, showed that only DaG20_ArsC3 is functional and able to confer arsenate resistance. Kinetic studies revealed that this enzyme uses thioredoxin as electron donor and therefore belongs to Staphylococcus aureus plasmid pI258 and Bacillus subtilis thioredoxin-coupled arsenate reductases family. Both enzymes from this family contain a potassium-binding site, but only in Sa_ArsC does potassium actually binds resulting in a lower K m. Important differences between the S. aureus and B. subtilis enzymes and DaG20_ArsC3 are observed. DaG20_ArsC3 contains only two (Asn10, Ser33) of the four (Asn10, Ser33, Thr63, Asp65) conserved amino acid residues that form the potassium-binding site and the kinetics is not significantly affected by the presence of either potassium or sulfate ions. Isothermal titration calorimetry measurements confirmed nonspecific binding of K(+) and Na(+), corroborating the non-relevance of these cations for catalysis. Furthermore, the low K m and high k cat values determined for DaG20_ArsC3 revealed that this enzyme is the most catalytically efficient potassium-independent arsenate reductase described so far and, for the first time indicates that potassium binding is not essential to have low K m, for Trx-arsenate reductases.

TupA: a tungstate binding protein in the periplasm of Desulfovibrio alaskensis G20.

The TupABC system is involved in the cellular uptake of tungsten and belongs to the ABC (ATP binding cassette)-type transporter systems. The TupA component is a periplasmic protein that binds tungstate anions, which are then transported through the membrane by the TupB component using ATP hydrolysis as the energy source (the reaction catalyzed by the ModC component). We report the heterologous expression, purification, determination of affinity binding constants and crystallization of the Desulfovibrio alaskensis G20 TupA. The tupA gene (locus tag Dde_0234) was cloned in the pET46 Enterokinase/Ligation-Independent Cloning (LIC) expression vector, and the construct was used to transform BL21 (DE3) cells. TupA expression and purification were optimized to a final yield of 10 mg of soluble pure protein per liter of culture medium. Native polyacrylamide gel electrophoresis was carried out showing that TupA binds both tungstate and molybdate ions and has no significant interaction with sulfate, phosphate or perchlorate. Quantitative analysis of metal binding by isothermal titration calorimetry was in agreement with these results, but in addition, shows that TupA has higher affinity to tungstate than molybdate. The protein crystallizes in the presence of 30% (w/v) polyethylene glycol 3350 using the hanging-drop vapor diffusion method. The crystals diffract X-rays beyond 1.4 Å resolution and belong to the P21 space group, with cell parameters a = 52.25 Å, b = 42.50 Å, c = 54.71 Å, ß = 95.43°. A molecular replacement solution was found, and the structure is currently under refinement.

FTIR spectroscopy of metalloproteins.

Absorption of infrared radiation by proteins gives important information about their structure and function. The most intense infrared bands correspond to the overlap of all the peptide bond absorption. Additionally, in many metalloproteins their prosthetic groups have intrinsic ligands or bind substrates/inhibitors that absorb intensively in the infrared. Here, we describe thoroughly several Fourier transform infrared methods for studying structure-function relationships in metalloproteins, using hydrogenases as an example.

Accurate calculations of geometries and singlet-triplet energy differences for active-site models of [NiFe] hydrogenase.

We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni(2+) models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22 electrons in 22 orbitals. This is so large that we need to use restricted active-space theory (RASPT2). The calculations predict that the singlet state is 48-57 kJ mol(-1) more stable than the triplet state for this model of the Ni-SIa state. However, in the [NiFe] hydrogenase protein, the structure around the Ni ion is far from the square-planar structure preferred by the singlet state. This destabilises the singlet state so that it is only ∼24 kJ mol(-1) more stable than the triplet state. Finally, we have studied how various density functional theory methods compare to the experimental, CCSD(T), CASPT2, and RASPT2 results. Semi-local functionals predict the best singlet-triplet energy differences, with BP86, TPSS, and PBE giving mean unsigned errors of 12-13 kJ mol(-1) (maximum errors of 25-31 kJ mol(-1)) compared to CCSD(T). For bond lengths, several methods give good results, e.g. TPSS, BP86, and M06, with mean unsigned errors of 2 pm for the bond lengths if relativistic effects are considered.

O2migration rates in [NiFe] hydrogenases. A joint approach combining free-energy calculations and kinetic modeling.

Hydrogenases are promising candidates for the catalytic production of green energy by means of biological ways. The major impediment to such a production is rooted in their inhibition under aerobic conditions. In this work, we model dioxygen migration rates in mutants of a hydrogenase of Desulfovibrio fructusovorans. The approach relies on the calculation of the whole potential of mean force for O2 migration within the wild-type as well as in V74M, V74F, and V74Q mutant channels. The three free-energy barriers along the entire migration pathway are converted into chemical rates through modeling based on Transition State Theory. The use of such a model recovers the trend of O2 migration rates among the series.