First-Principles Calculations on Ni,Fe-Containing Carbon Monoxide Dehydrogenases Reveal Key Stereoelectronic Features for Binding and Release of CO2 to/from the C-Cluster.

In view of the depletion of fossil fuel reserves and climatic effects of greenhouse gas emissions, Ni,Fe-containing carbon monoxide dehydrogenase (Ni-CODH) enzymes have attracted increasing interest in recent years for their capability to selectively catalyze the reversible reduction of CO2 to CO (CO2 + 2H+ + 2e- ⇌ CO + H2O). The possibility of converting the greenhouse gas CO2 into useful materials that can be used as synthetic building blocks or, remarkably, as carbon fuels makes Ni-CODH a very promising target for reverse-engineering studies. In this context, in order to provide insights into the chemical principles underlying the biological catalysis of CO2 activation and reduction, quantum mechanics calculations have been carried out in the framework of density functional theory (DFT) on different-sized models of the Ni-CODH active site. With the aim of uncovering which stereoelectronic properties of the active site (known as the C-cluster) are crucial for the efficient binding and release of CO2, different coordination modes of CO2 to different forms and redox states of the C-cluster have been investigated. The results obtained from this study highlight the key role of the protein environment in tuning the reactivity and the geometry of the C-cluster. In particular, the protonation state of His93 is found to be crucial for promoting the binding or the dissociation of CO2. The oxidation state of the C-cluster is also shown to be critical. CO2 binds to Cred2 according to a dissociative mechanism (i.e., CO2 binds to the C-cluster after the release of possible ligands from Feu) when His93 is doubly protonated. CO2 can also bind noncatalytically to Cred1 according to an associative mechanism (i.e., CO2 binding is preceded by the binding of H2O to Feu). Conversely, CO2 dissociates when His93 is singly protonated and the C-cluster is oxidized at least to the Cint redox state.

Rational Design of Fe2 (µ-PR2 )2 (L)6 Coordination Compounds Featuring Tailored Potential Inversion.

It was recently discovered that some redox proteins can thermodynamically and spatially split two incoming electrons towards different pathways, resulting in the one-electron reduction of two different substrates, featuring reduction potential respectively higher and lower than the parent reductant. This energy conversion process, referred to as electron bifurcation, is relevant not only from a biochemical perspective, but also for the ground-breaking applications that electron-bifurcating molecular devices could have in the field of energy conversion. Natural electron-bifurcating systems contain a two-electron redox centre featuring potential inversion (PI), i. e. with second reduction easier than the first. With the aim of revealing key factors to tailor the span between first and second redox potentials, we performed a systematic density functional study of a 26-molecule set of models with the general formula Fe2 (µ-PR2 )2 (L)6 . It turned out that specific features such as i) a Fe-Fe antibonding character of the LUMO, ii) presence of electron-donor groups and iii) low steric congestion in the Fe's coordination sphere, are key ingredients for PI. In particular, the synergic effects of i)-iii) can lead to a span between first and second redox potentials larger than 700 mV. More generally, the "molecular recipes" herein described are expected to inspire the synthesis of Fe2 P2 systems with tailored PI, of primary relevance to the design of electron-bifurcating molecular devices.

Reactivation of the Ready and Unready Oxidized States of [NiFe]-Hydrogenases: Mechanistic Insights from DFT Calculations.

The apparently simple dihydrogen formation from protons and electrons (2H+ + 2e- ⇄ H2) is one of the most challenging reactions in nature. It is catalyzed by metalloenzymes of amazing complexity, called hydrogenases. A better understanding of the chemistry of these enzymes, especially that of the [NiFe]-hydrogenases subgroup, has important implications for production of H2 as alternative sustainable fuel. In this work, reactivation mechanism of the oxidized and inactive Ni-B and Ni-A states of the [NiFe]-hydrogenases active site has been investigated using density functional theory. Results obtained from this study show that one-electron reduction and protonation of the active site promote the removal of the bridging hydroxide ligand contained in Ni-B and Ni-A. However, this process is sufficient to activate only the Ni-B state. H2 binding to the active site is required to convert Ni-A to the active Ni-SIa state. Here, we also propose a reasonable structure for the spectroscopically well-characterized Ni-SIr and Ni-SU species, formed respectively from the one-electron reduction of Ni-B and Ni-A. Ni-SIr, depending on the pH at which the reaction occurs, features a bridging hydroxide ligand or a water molecule terminally coordinated to the Ni atom, whereas in Ni-SU a water molecule is terminally coordinated to the Fe atom, and the Cys64 residue is oxidized to sulfenate. The sulfenate oxygen atom in the Ni-A state affects the stereoelectronic properties of the binuclear cluster by modifying the coordination geometry of Ni, and consequently, by switching the regiochemistry of H2O and H2 binding from the Ni to the Fe atom. This effect is predicted to be at the origin of the different reactivation kinetics of the oxidized and inactive Ni-B and Ni-A states.

Redox Potentials of Small Inorganic Radicals and Hexa-Aquo Complexes of First-Row Transition Metals in Water: A DFT Study Based on the Grand Canonical Ensemble.

The potentials of redox systems involving nitrogen, oxygen, and metal ions of the first-row transition series have been computed according to the general approach of the grand canonical ensemble, which leads to the equilibrium value of the reduction potential via a (complete) sampling of configuration space at a given temperature. The approach is a single configuration approach in the sense that identical molecular structures are sampled for both the oxidized and reduced species considered in water solution. In this study, the solute and a cluster of 11-12 water molecules are treated explicitly at the same level of theory and embedded in a continuum solvent. The molecular energies are computed in the framework of the density functional theory. Our approach is basically different from the approach based on the ThermoDynamic Cycle involving gas-phase calculations of the electron affinity of the oxidized species, corrected by the differential hydration energy (obtained from continuum solvent models only) between oxidized and reduced forms. The calculated redox potentials are in agreement with the available experimental data much closer than other results so far presented in the literature. Our results are very satisfactory also in the case of the 3+/2+ redox states of the first-row transition metals, i.e., systems with a high positive charge for which enhanced effects of the solvent are expected.

Catalytic Mechanism of Fungal Lytic Polysaccharide Monooxygenases Investigated by First-Principles Calculations.

Lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes that facilitate the degradation of recalcitrant polysaccharides by the oxidative cleavage of glycosidic bonds. They are gaining rapidly increasing attention as key players in biomass conversion, especially for the production of second-generation biofuels. Elucidation of the detailed mechanism of the LPMO reaction is a major step toward the assessment and optimization of LPMO efficacy in industrial biotechnology, paving the way to utilization of sustainable fuel sources. Here, we used density functional theory calculations to study the reaction pathways suggested to date, exploiting a very large active-site model for a fungal AA9 LPMO and using a celloheptaose unit as a substrate mimic. We identify a copper oxyl intermediate as being responsible for H-atom abstraction from the substrate, followed by a rapid, water-assisted hydroxyl rebound, leading to substrate hydroxylation.

Theoretical investigation of aerobic and anaerobic oxidative inactivation of the [NiFe]-hydrogenase active site.

The extraordinary capability of [NiFe]-hydrogenases to catalyse the reversible interconversion of protons and electrons into dihydrogen (H2) has stimulated numerous experimental and theoretical studies addressing the direct utilization of these enzymes in H2 production processes. Unfortunately, the introduction of these natural H2-catalysts in biotechnological applications is limited by their inhibition under oxidising (aerobic and anaerobic) conditions. With the aim of contributing to overcome this limitation, we studied the oxidative inactivation mechanism of [NiFe]-hydrogenases by performing Density Functional Theory (DFT) calculations on a very large model of their active site in which all the amino acids forming the first and second coordination spheres of the NiFe cluster have been explicitly included. We identified an O2 molecule and two H2O molecules as sources of the two oxygen atoms that are inserted at the active site of the inactive forms of the enzyme (Ni-A and Ni-B) under aerobic and anaerobic conditions, respectively. Furthermore, our results support the experimental evidence that the Ni-A-to-Ni-B ratio strongly depends on the number of reducing equivalents available for the process and on the oxidizing conditions under which the reaction takes place.

Theoretical insights into [NiFe]-hydrogenases oxidation resulting in a slowly reactivating inactive state.

[NiFe]-hydrogenases catalyse the relevant H2 â†’ 2H+ + 2e- reaction. Aerobic oxidation or anaerobic oxidation of this enzyme yields two inactive states called Ni-A and Ni-B. These states differ for the reactivation kinetics which are slower for Ni-A than Ni-B. While there is a general consensus on the structure of Ni-B, the nature of Ni-A is still controversial. Indeed, several crystallographic structures assigned to the Ni-A state have been proposed, which, however, differ for the nature of the bridging ligand and for the presence of modified cysteine residues. The spectroscopic characterization of Ni-A has been of little help due to small differences of calculated spectroscopic parameters, which does not allow to discriminate among the various forms proposed for Ni-A. Here, we report a DFT investigation on the nature of the Ni-A state, based on systematic explorations of conformational and configurational space relying on accurate energy calculations, and on comparisons of theoretical geometries with the X-ray structures currently available. The results presented in this work show that, among all plausible isomers featuring various protonation patterns and oxygenic ligands, the one corresponding to the crystallographic structure recently reported by Volbeda et al. (J Biol Inorg Chem 20:11-22, 19)-featuring a bridging hydroxide ligand and the sulphur atom of Cys64 oxidized to bridging sulfenate-is the most stable. However, isomers with cysteine residues oxidized to terminal sulfenate are very close in energy, and modifications in the network of H-bond with neighbouring residues may alter the stability order of such species.

The large subunit of the regulatory [NiFe]-hydrogenase from Ralstonia eutropha - a minimal hydrogenase?

Chemically synthesized compounds that are capable of facilitating the reversible splitting of dihydrogen into protons and electrons are rare in chemists' portfolio. The corresponding biocatalysts - hydrogenases - are, however, abundant in the microbial world. [NiFe]-hydrogenases represent a major subclass and display a bipartite architecture, composed of a large subunit, hosting the catalytic NiFe(CO)(CN)2 cofactor, and a small subunit whose iron-sulfur clusters are responsible for electron transfer. To analyze in detail the catalytic competence of the large subunit without its smaller counterpart, we purified the large subunit HoxC of the regulatory [NiFe]-hydrogenase of the model H2 oxidizer Ralstonia eutropha to homogeneity. Metal determination and infrared spectroscopy revealed a stoichiometric loading of the metal cofactor. This enabled for the first time the determination of the UV-visible extinction coefficient of the NiFe(CO)(CN)2 cofactor. Moreover, the absence of disturbing iron-sulfur clusters allowed an unbiased look into the low-spin Fe2+ of the active site by Mössbauer spectroscopy. Isolated HoxC was active in catalytic hydrogen-deuterium exchange, demonstrating its capacity to activate H2. Its catalytic activity was drastically lower than that of the bipartite holoenzyme. This was consistent with infrared and electron paramagnetic resonance spectroscopic observations, suggesting that the bridging position between the active site nickel and iron ions is predominantly occupied by water-derived ligands, even under reducing conditions. In fact, the presence of water-derived ligands bound to low-spin Ni2+ was reflected by the absorption bands occurring in the corresponding UV-vis spectra, as revealed by time-dependent density functional theory calculations conducted on appropriate in silico models. Thus, the isolated large subunits indeed represent simple [NiFe]-hydrogenase models, which could serve as blueprints for chemically synthesized mimics. Furthermore, our data point to a fundamental role of the small subunit in preventing water access to the catalytic center, which significantly increases the H2 splitting capacity of the enzyme.

A theoretical study on the reactivity of the Mo/Cu-containing carbon monoxide dehydrogenase with dihydrogen.

The Mo/Cu-dependent CO dehydrogenase from Oligotropha carboxidovorans is an enzyme that is able to catalyze CO oxidation to CO2; moreover, it can also oxidize H2, thus eliciting a characteristic EPR signal. Interestingly, the Ag-substituted enzyme form proved unable to catalyze H2 oxidation. In the present contribution, we characterized the reactivity of the enzyme with H2 by quantum-chemical calculations. It was found that dihydrogen binding to the wild-type enzyme requires significant structural rearrangements of the active site Theoretical EPR spectra for plausible H2-bound models of the partially reduced, paramagnetic active site are also presented and compared with the experimental counterpart. Finally, density functional theory modeling shows that Ag substitution impairs H2 binding at the active site.

Copper ion interaction with the RNase catalytic site fragment of the angiogenin protein: an experimental and theoretical investigation.

The angiogenin protein (Ang) is a member of the vertebrate-specific secreted ribonucleases and one of the most potent angiogenic factors known. Ang is a normal constituent of human plasma and its concentration increases under some physiological and pathological conditions to promote neovascularization. Ang was originally identified as an angiogenic tumour factor, but its biological activity has been found to extend from inducing angiogenesis to promoting cell survival in different neurodegenerative diseases. Ang exhibits weak ribonucleolytic activity, which is critical for its biological functions. The RNase catalytic sites are two histidine residues, His-13 and His-114, and the lysine Lys-40. Copper is also an essential cofactor in angiogenesis and influences angiogenin's biological properties. The main Cu(ii) anchoring site of Ang is His-114, where metal binding inhibits RNase activity of the protein. To reveal the Cu(ii) coordination environment in the C-terminal domain of the Ang protein, we report on the characterization, by means of potentiometric, voltammetric, and spectroscopic (CD, UV-Vis and EPR) methods and DFT calculations, of Cu(ii) complexes formed with a peptide fragment including the Ang sequence 112-117 (PVHLDQ). Potentiometric titrations indicated that [CuLH-2] is the predominant species at physiological pH. EPR, voltammetric data and DFT calculations are consistent with a CuN3O2 coordination mode in which a distorted square pyramidal arrangement of the peptide was observed with the equatorial positions occupied by the nitrogen atoms of the deprotonated amides of the Asp and Leu residues, the δ-N atom of histidine and the oxygen atom of the aspartic carboxylic group. Moreover, two analogous peptides encompassing the PVHLNQ and LVHLDQ sequences were also characterized by using thermodynamic, spectroscopic and DFT studies to reveal the role they play in Cu(ii) complex formation by the carboxylate side chain of the Asp and Pro residues, a known breaking-point in metal coordination.